Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistent scheme”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Development of a Hybrid RANS/LES Method for Turbulent Mixing Layers

Significant research has been underway for several years in NASA Glenn Research Center's nozzle branch to develop advanced computational methods for simulating turbulent flows in exhaust nozzles. The primary efforts of this research have concentrated on improving our ability to calculate the turbulent mixing layers that dominate flows both in the exhaust systems of modern-day aircraft and in those of hypersonic vehicles under development. As part of these efforts, a hybrid numerical method was recently developed to simulate such turbulent mixing layers. The method developed here is intended for configurations in which a dominant structural feature provides an unsteady mechanism to drive the turbulent development in the mixing layer. Interest in Large Eddy Simulation (LES) methods have increased in recent years, but applying an LES method to calculate the wide range of turbulent scales from small eddies in the wall-bounded regions to large eddies in the mixing region is not yet possible with current computers. As a result, the hybrid method developed here uses a Reynolds-averaged Navier-Stokes (RANS) procedure to calculate wall-bounded regions entering a mixing section and uses a LES procedure to calculate the mixing-dominated regions. A numerical technique was developed to enable the use of the hybrid RANS-LES method on stretched, non-Cartesian grids. With this technique, closure for the RANS equations is obtained by using the Cebeci-Smith algebraic turbulence model in conjunction with the wall-function approach of Ota and Goldberg. The LES equations are closed using the Smagorinsky subgrid scale model. Although the function of the Cebeci-Smith model to replace all of the turbulent stresses is quite different from that of the Smagorinsky subgrid model, which only replaces the small subgrid turbulent stresses, both are eddy viscosity models and both are derived at least in part from mixing-length theory. The similar formulation of these two models enables the RANS and LES equations to be solved with a single solution scheme and computational grid. The hybrid RANS-LES method has been applied to a benchmark compressible mixing layer experiment in which two isolated supersonic streams, separated by a splitter plate, provide the flows to a constant-area mixing section. Although the configuration is largely two dimensional in nature, three-dimensional calculations were found to be necessary to enable disturbances to develop in three spatial directions and to transition to turbulence. The flow in the initial part of the mixing section consists of a periodic vortex shedding downstream of the splitter plate trailing edge. This organized vortex shedding then rapidly transitions to a turbulent structure, which is very similar to the flow development observed in the experiments. Although the qualitative nature of the large-scale turbulent development in the entire mixing section is captured well by the LES part of the current hybrid method, further efforts are planned to directly calculate a greater portion of the turbulence spectrum and to limit the subgrid scale modeling to only the very small scales. This will be accomplished by the use of higher accuracy solution schemes and more powerful computers, measured both in speed and memory capabilities.

Georgiadis, Nicholas J.↗
NASA NTRS1973

Fuel cell technology program

A fuel cell technology program was established to advance the state-of-the-art of hydrogen-oxygen fuel cells using low temperature, potassium hydroxide electrolyte technology as the base. Program tasks are described consisting of baseline cell design and stack testing, hydrogen pump design and testing, and DM-2 powerplant testing and technology extension efforts. A baseline cell configuration capable of a minimum of 2000 hours of life was defined. A 6-cell prototype stack, incorporating most of the scheme cell features, was tested for a total of 10,497 hours. A 6-cell stack incorporating all of the design features was tested. The DM-2 powerplant with a 34 cell stack, an accessory section packaged in the basic configuration anticipated for the space shuttle powerplant and a powerplant control unit, was defined, assembled, and tested. Cells were used in the stack and a drag-type hydrogen pump was installed in the accessory section. A test program was established, in conjunction with NASA/JSC, based on space shuttle orbiter mission. A 2000-hour minimum endurance test and a 5000-hour goal were set and the test started on August 8, 1972. The 2000-hour milestone was completed on November 3, 1972. On 13 March 1973, at the end of the thirty-first simulated seven-day mission and 5072 load hours, the test was concluded, all goals having been met. At this time, the DM-2 was in excellent condition and capable of additional endurance.

Source record↗
NASA NTRS1973

Determination of normal points for the relativity test of the celestial mechanics experiment

The relativity test is based on Mariner 9 tracking data acquired by the Deep Space Network. These data consist of two types of measurements: (1) two-way frequency difference (Doppler) measurements that are proportional to the tracking station-spacecraft range rate, and (2) signal round trip time-delay measurements that are proportional to the station-spacecraft range. A data compression scheme has been employed to alleviate accuracy and cost difficulties by (1) using the Doppler data to solve for the spacecraft orbit and to relate spacecraft position to the center of mass of Mars; and (2) combining this result with the station-spacecraft range measurements to obtain measurements of the Earth-Mars distance, called normal points, that contain all of the information pertinent to a detailed knowledge of the Earth-Mars motion.

Jordan, J. F.↗
NASA NTRS1982

Seasat synthetic-aperture radar data reduction using parallel programmable array processors

This paper presents a digital processing system that produces the Seasat synthetic-aperture radar (SAR) imagery. The system consists of a SEL 32/77 host minicomputer and three AP-120B array processors. The partitioning of the SAR processing functions and the design of software modules is described. The rationale for selecting the parallel array processor architecture and the methodology for developing the parallel processing scheme on this system is described. This system attains a Seasat SAR data reduction speed of 2.5 h per 25-m resolution 4-look and 100 km x 100 km image frame. A preliminary performance evaluation of this parallel processing system and potential future applications for remote sensing data reduction are described.

Wu, C.↗
NASA NTRS1983

Spectral consequences of photoreceptor sampling in the rhesus retina

Optical transforms were used to compute the power spectra of rhesus cones treated as arrays of image sampling points. Spectra were obtained for the central fovea, parafovea, periphery, and far periphery. All were consistent with a novel spatial sampling principle that introduces minimal noise for spatial frequencies below the Nyquist limits implied by local receptor densities, while frequencies above the nominal Nyquist limits are not converted into conspicuous moiree patterns, but instead are scattered into broadband noise. This sampling scheme allows the visual system to escape aliasing distortion despite a large mismatch between retinal image bandwidth and the Nyquist limits implied by extrafoveal cone densities.

Yellott, J. I., Jr.↗
NASA NTRS1995

AC motor controller with 180 degree conductive switches

An ac motor controller is operated by a modified time-switching scheme where the switches of the inverter are on for electrical-phase-and-rotation intervals of 180.degree. as opposed to the conventional 120.degree.. The motor is provided with three-phase drive windings, a power inverter for power supplied from a dc power source consisting of six switches, and a motor controller which controls the current controlled switches in voltage-fed mode. During full power, each switch is gated continuously for three successive intervals of 60.degree. and modulated for only one of said intervals. Thus, during each 60.degree. interval, the two switches with like signs are on continuously and the switch with the opposite sign is modulated.

Oximberg, Carol A.↗
NASA NTRS2013

A Hybrid Fiber/Solid-State Regenerative Amplifier with Tunable Pulse Widths for Satellite Laser Ranging

A fiber/solid-state hybrid seeded regenerative amplifier, capable of achieving high output energy with tunable pulse widths, has been developed for satellite laser ranging applications. The regenerative amplifier cavity uses a pair of Nd:YAG zigzag slabs oriented orthogonally to one another in order to make thermal lensing effects symmetrical and simplify optical correction schemes. The seed laser used is a fiber-coupled 1,064-nm narrowband (<0.02 nm) diode laser that is discretely driven in a new short-pulsed mode, enabling continuously tunable seed pulse widths in the 0.2-to-0.4-ns range. The amplifier gain unit consists of a pair of Brewster-cut 6-bounce zigzag Nd:YAG laser slabs, oriented 90deg relative to each other in the amplifier head. This arrangement creates a net-symmetrical thermal lens effect (an opposing singleaxis effect in each slab), and makes thermo-optical corrections simple by optimizing the curvature of the nearest cavity mirror. Each slab is pumped by a single 120-W, pulsed 808-nm laser diode array. In this configuration, the average pump beam distribution in the slabs had a 1-D Gaussian shape, which matches the estimated cavity mode size. A half-wave plate between the slabs reduces losses from Fresnel reflections due to the orthogonal slabs Brewster-cut end faces. Successful "temporal" seeding of the regenerative amplifier cavity results in a cavity Q-switch pulse envelope segmenting into shorter pulses, each having the width of the input seed, and having a uniform temporal separation corresponding to the cavity round-trip time of approx. =10 ns. The pulse energy is allowed to build on successive passes in the regenerative amplifier cavity until a maximum is reached, (when cavity gains and losses are equal), after which the pulse is electro- optically switched out on the next round trip The overall gain of the amplifier is approx. =82 dB (or a factor of 1.26 million). After directing the amplified output through a LBO frequency doubling crystal, approx. = 2.1 W of 532-nm output (>1 mJ) was measured. This corresponds to a nonlinear conversion efficiency of >60%. Furthermore, by pulse pumping this system, a single pulse per laser shot can be created for the SLR (satellite laser ranging) measurement, and this can be ejected into the instrument. This is operated at the precise frequency needed by the measurement, as opposed to commercial short-pulsed, mode-locked systems that need to operate in a continuous fashion, or CW (continuous wave), and create pulses at many MHz. Therefore, this design does not need to throw away or dump 99% of the laser energy to produce what is required; this system can be far smaller, more efficient, cheaper, and readily deployed in the field when packaged efficiently. Finally, by producing custom diode seed pulses electronically, two major advantages over commercial systems are realized: First, this pulse shape is customizable and not affected by the cavity length or gain of the amplifier cavity, and second, it can produce adjustable (selectable) pulse widths by simply adding multiple seed diodes and coupling each into commercial, low-cost fiber-optic combiners.

Coyle, Barry↗
NASA NTRS1991

Hybrid solution for high-speed target acquisition and identification systems

A typical hierarchy for a general object recognition problem consists of object detection, classification and identification. This paper establishes necessary building blocks required for high-speed object recognition applications. An architecture that combines digital and optical processing, exploiting current image processing techniques for detection and classification, and optical processing hardware is described. An optical processing scheme is suggested for the identification aspect. Numerical results of each proposed concept are presented.

Udomkesmalee, Suraphol↗
NASA NTRS1995

Sensitivity analysis of frictional contact response of axisymmetric composite structures

A computational procedure is presented for evaluating the sensitivity coefficients of the static frictional contact response of axisymmetric composite structures. The structures are assumed to consist of an arbitrary number of perfectly bonded homogeneous anisotropic layers. The material of each layer is assumed to be hyperelastic, and the effect of geometric nonlinearity is included. The sensitivity coefficients measure the sensitivity of the response variations in different material, lamination and geometric parameters of the structure. A displacement finite element model is used for the discretization. The normal contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of nodal displacements, and Lagrange multipliers associated with the contact conditions. The Lagrange multipliers are allowed to be discontinuous at interelement boundaries. Tangential contact conditions are incorporated by using a penalty method in conjunction with the classical Coulomb's friction model. The Newton-Raphson iterative scheme is used for the solution of the resulting nonlinear algebraic equations, and for the determination of the contact region, contact conditions (sliding or sticking), and the contact pressures. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented for the frictional contact of a composite spherical cap pressed against a rigid plate.

Karaoglan, I.↗
NASA NTRS1992

A numerical study of three-dimensional surface tension driven convection with fre surface deformation

The steady three-dimensional thermocapillary motion with a deformable free surface is studied numerically in both normal and zero gravity environments. Flow configurations consist of a square cavity heated from the side. In the analysis, the free surface is allowed to deform and the grid distribution is adapted to the surface deformation. The divergence-free condition is satisfied by using a dual time-stepping approach in the numerical scheme. Convective flux derivatives are evaluated using a third-order accurate upwind-biased flux-split differencing technique. The numerical solutions at the midplane of the square cavity are compared with the results from two-dimensional calculations. In addition, numerial results for cases under zero and normal gravity conditions are compared. Significantly different flow structures and surface deformation have been observed. The comparison of calculated results will be compared with experimental data in the updated version of this paper.

Hsieh, Kwang-Chung↗
NASA NTRS1995

Calculations of Diffuser Flows with an Anisotropic K-Epsilon Model

A newly developed anisotropic K-epsilon model is applied to calculate three axisymmetric diffuser flows with or without separation. The new model uses a quadratic stress-strain relation and satisfies the realizability conditions, i.e., it ensures both the positivity of the turbulent normal stresses and the Schwarz' inequality between any fluctuating velocities. Calculations are carried out with a finite-element method. A second-order accurate, bounded convection scheme and sufficiently fine grids are used to ensure numerical credibility of the solutions. The standard K-epsilon model is also used in order to highlight the performance of the new model. Comparison with the experimental data shows that the anisotropic K-epsilon model performs consistently better than does the standard K-epsilon model in all of the three test cases.

Zhu, J.↗
NASA NTRS2001

A Coupled Plasma and Sheath Model for High Density Reactors

We present a coupled plasma and collisionless; sheath model for the simulation of high density plasma processing reactors. Due to inefficiencies in numerical schemes and the resulting computational burden, a coupled multidimensional plasma and sheath simulation has not been possible model for gas mixtures and high density reactors of practical interest. In this work we demonstrate that with a fully implicit algorithm and a refined computational mesh, a self-consistent plasma and sheath simulation is feasible. We discuss the details of the model equations, the importance of ion inertia, and the resulting sheath profiles for argon and chlorine plasmas. We find that at low operating pressures (10-30 mTorr), the charge separation occurs only within a 0.5 mm layer near the surface in a 300 mm inductively coupled plasma etch reactor. A unified model eliminates the use of off-line or loosely coupled sheath models with simplifying assumptions which generally lead to uncertainties in ion flux and sheath electrical properties.

Deepak, Bose↗
NASA NTRS2001

Formation of Electrostatic Potential Drops in the Auroral Zone

In order to examine the self-consistent formation of large-scale quasi-static parallel electric fields in the auroral zone on a micro/meso scale, a particle in cell simulation has been developed. The code resolves electron Debye length scales so that electron micro-processes are included and a variable grid scheme is used such that the overall length scale of the simulation is of the order of an Earth radii along the magnetic field. The simulation is electrostatic and includes the magnetic mirror force, as well as two types of plasmas, a cold dense ionospheric plasma and a warm tenuous magnetospheric plasma. In order to study the formation of parallel electric fields in the auroral zone, different magnetospheric ion and electron inflow boundary conditions are used to drive the system. It has been found that for conditions in the primary (upward) current region an upward directed quasi-static electric field can form across the system due to magnetic mirroring of the magnetospheric ions and electrons at different altitudes. For conditions in the return (downward) current region it is shown that a quasi-static parallel electric field in the opposite sense of that in the primary current region is formed, i.e., the parallel electric field is directed earthward. The conditions for how these different electric fields can be formed are discussed using satellite observations and numerical simulations.

Schriver, D.↗
NASA NTRS2017

Geometrically Flexible and Efficient Flow Analysis of High Speed Vehicles Via Domain Decomposition, Part 1: Unstructured-Grid Solver for High Speed Flows

The ability to solve the equations governing the hypersonic turbulent flow of a real gas on unstructured grids using a spatially-elliptic, 2nd-order accurate, cell-centered, finite-volume method has been recently implemented in the VULCAN-CFD code. This paper describes the key numerical methods and techniques that were found to be required to robustly obtain accurate solutions to hypersonic flows on non-hex-dominant unstructured grids. The methods and techniques described include: an augmented stencil, weighted linear least squares, cell-average gradient method, a robust multidimensional cell-average gradient-limiter process that is consistent with the augmented stencil of the cell-average gradient method and a cell-face gradient method that contains a cell skewness sensitive damping term derived using hyperbolic diffusion based concepts. A data-parallel matrix-based symmetric Gauss-Seidel point-implicit scheme, used to solve the governing equations, is described and shown to be more robust and efficient than a matrix-free alternative. In addition, a y+ adaptive turbulent wall boundary condition methodology is presented. This boundary condition methodology is deigned to automatically switch between a solve-to-the-wall and a wall-matching-function boundary condition based on the local y+ of the 1st cell center off the wall. The aforementioned methods and techniques are then applied to a series of hypersonic and supersonic turbulent flat plate unit tests to examine the efficiency, robustness and convergence behavior of the implicit scheme and to determine the ability of the solve-to-the-wall and y+ adaptive turbulent wall boundary conditions to reproduce the turbulent law-of-the-wall. Finally, the thermally perfect, chemically frozen, Mach 7.8 turbulent flow of air through a scramjet flow-path is computed and compared with experimental data to demonstrate the robustness, accuracy and convergence behavior of the unstructured-grid solver for a realistic 3-D geometry on a non-hex-dominant grid.

White, Jeffery A.↗
NASA NTRS2020

Reduction of structural impacts and distinction of photosynthetic pathways in a global estimation of GPP from space-borne solar-induced chlorophyll fluorescence

Quantifying global photosynthesis remains a challenge due to a lack of accurate remote sensing proxies. Solar-induced chlorophyll fluorescence (SIF) has been shown to be a good indicator of photosynthetic activity across various spatial scales. However, a global and spatially challenging estimate of terrestrial gross primary production (GPP) based on satellite SIF remains unresolved due to the confounding effects of species-specific physical and physiological traits and external factors, such as canopy structure or photosynthetic pathway (C3orC4). Here we analyze an ensemble of far-red SIF data from OCO-2 satellite and ground observations at multiple sites, using the spectral invariant theory to reduce the effects of canopy structure and to retrieve a structure-corrected total canopy SIF emission (SIF total). We find that the relationships between observed canopy-leavingSIF and ecosystem GPP vary significantly among biomes. In contrast, the relationships between SIF total and GPP converge around two unique models, one for C3and one for C4plants. We show that the two single empirical models can be used to globally scale satellite SIF observations to terrestrial GPP. We obtain an independent estimate of global terrestrial GPP of 129.56 ± 6.54 PgC/year for the 2015–2017 period, which is consistent with the state-of-the-art data- and process-oriented models. The new GPP product shows improved sensitivity to previously undetected ‘hotspots’ of productivity, being able to resolve the double-peak in GPP due to rotational cropping systems. We suggest that the direct scheme to estimate GPP presented here, which is based on satelliteSIF, may open up new possibilities to resolve the dynamics of global terrestrial GPP across space and time.

scale satellite SIF observations to terrestrial GP↗
NASA NTRS1973

An unsupervised classification technique for multispectral remote sensing data.

Description of a two-part clustering technique consisting of (a) a sequential statistical clustering, which is essentially a sequential variance analysis, and (b) a generalized K-means clustering. In this composite clustering technique, the output of (a) is a set of initial clusters which are input to (b) for further improvement by an iterative scheme. This unsupervised composite technique was employed for automatic classification of two sets of remote multispectral earth resource observations. The classification accuracy by the unsupervised technique is found to be comparable to that by traditional supervised maximum-likelihood classification techniques.

Su, M. Y.↗
NASA NTRS1982

Mature thunderstorm cloud top structure - Three-dimensional numerical simulation versus satellite observations

Preliminary results of four runs with a three-dimensional model of the effects of vertical wind shear on cloud top height/temperature structure and the internal properties of isolate midlatitude thunderstorms are reported. The model is being developed as an aid to analyses of GEO remote sensing satellite data. The grid is a 27 x 27 x 20 mesh with 2 km horizontal resolution and 0.9 vertical resolution. The total grid is 54 km on a side and 18 km deep. A second-order Crowley scheme for advecting momentum is extended with a third-order correction for spatial truncation error, and the earth-relative horizontal surface wind components are decreased to 50 percent of their values at 0.45 km. A temperature increase with height is included, together with an initial impulse consisting of a nonrotating cylindrical weak buoyant updraft 10 km in radius. The results of the runs are discussed in terms of the time variation of the vertical velocity extrema, the effects of strong and weak shear on a storm, the cloud top height, the Lagrangian dynamics of a thermal couplet, and data from a real storm.

Schlesinger, R. E.↗
NASA NTRS1984

A combined direct/inverse three-dimensional transonic wing design method for vector computers

A three-dimensional transonic-wing design algorithm for vector computers is developed, and the results of sample computations are presented graphically. The method incorporates the direct/inverse scheme of Carlson (1975), a Cartesian grid system with boundary conditions applied at a mean plane, and a potential-flow solver based on the conservative form of the full potential equation and using the ZEBRA II vectorizable solution algorithm of South et al. (1980). The accuracy and consistency of the method with regard to direct and inverse analysis and trailing-edge closure are verified in the test computations.

Weed, R. A.↗