SEARCH · Search NASA

Results for “Gradient Limiters”

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.

SourceResearch connectionFrom yearThrough year

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

At least 19 records

All sources NASA NTRS DOE OSTI

NASA NTRS2002

Stereo transparency and the disparity gradient limit

Several studies (Vision Research 15 (1975) 583; Perception 9 (1980) 671) have shown that binocular fusion is limited by the disparity gradient (disparity/distance) separating image points, rather than by their absolute disparity values. Points separated by a gradient >1 appear diplopic. These results are sometimes interpreted as a constraint on human stereo matching, rather than a constraint on fusion. Here we have used psychophysical measurements on stereo transparency to show that human stereo matching is not constrained by a gradient of 1. We created transparent surfaces composed of many pairs of dots, in which each member of a pair was assigned a disparity equal and opposite to the disparity of the other member. For example, each pair could be composed of one dot with a crossed disparity of 6' and the other with uncrossed disparity of 6', vertically separated by a parametrically varied distance. When the vertical separation between the paired dots was small, the disparity gradient for each pair was very steep. Nevertheless, these opponent-disparity dot pairs produced a striking appearance of two transparent surfaces for disparity gradients ranging between 0.5 and 3. The apparent depth separating the two transparent planes was correctly matched to an equivalent disparity defined by two opaque surfaces. A test target presented between the two transparent planes was easily detected, indicating robust segregation of the disparities associated with the paired dots into two transparent surfaces with few mismatches in the target plane. Our simulations using the Tsai-Victor model show that the response profiles produced by scaled disparity-energy mechanisms can account for many of our results on the transparency generated by steep gradients.

NASA Discipline Space Human Factors↗

NASA NTRS2019

Efficient and Robust Weighted Least-Squares Cell-Average Gradient Construction Methods for the Simulation of Scramjet Flows

White, Jeffery A.↗

NASA NTRSDate not supplied

A 3-D Nodal-Averaged Gradient Approach for Unstructured-Grid Cell-Centered Finite-Volume Methods for Application to Turbulent Hypersonic Flow

A 2-D nodal weighted least-squares gradient method and a related face-averaged nodal gradient approach that were developed for use with triangular grids are extended to 3-D for use with tetrahedral grids. In addition, a method, developed in 2-D, to stabilize the iterative convergence of these methods on quadrilateral cells is described and extended to 3-D and remedies are investigated to determine the nodal gradient averaging approach most suitable for use with grids made up of hexahedral, prismatic, pyramidal and tetrahedral cells. Moreover, due to an interest in hypersonic flow, a robust multidimensional gradient limiter procedure that is consistent with the stencil used to construct the nodal gradients is described. Finally, we demonstrate that the resulting 3-D methods are sufficiently robust for use in scramjet computations through the solution of three canonical turbulent hypersonic flow problems as well as a physically realistic 3-D scramjet inlet geometry.

Jeffery A White↗

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 NTRS2019

Weighted Least-Squares Cell-Average Gradient Construction Methods for the VULCAN-CFD Second-Order Accurate Unstructured-Grid Cell-Centered Finite-Volume Solver

White, Jeffery A.↗

NASA NTRS2013

Atomic Force Microscope Mediated Chromatography

The atomic force microscope (AFM) is used to inject a sample, provide shear-driven liquid flow over a functionalized substrate, and detect separated components. This is demonstrated using lipophilic dyes and normal phase chromatography. A significant reduction in both size and separation time scales is achieved with a 25-micron-length column scale, and one-second separation times. The approach has general applications to trace chemical and microfluidic analysis. The AFM is now a common tool for ultra-microscopy and nanotechnology. It has also been demonstrated to provide a number of microfluidic functions necessary for miniaturized chromatography. These include injection of sub-femtoliter samples, fluidic switching, and sheardriven pumping. The AFM probe tip can be used to selectively remove surface layers for subsequent microchemical analysis using infrared and tip-enhanced Raman spectroscopy. With its ability to image individual atoms, the AFM is a remarkably sensitive detector that can be used to detect separated components. These diverse functional components of microfluidic manipulation have been combined in this work to demonstrate AFM mediated chromatography. AFM mediated chromatography uses channel-less, shear-driven pumping. This is demonstrated with a thin, aluminum oxide substrate and a non-polar solvent system to separate a mixture of lipophilic dyes. In conventional chromatographic terms, this is analogous to thin-layer chromatography using normal phase alumina substrate with sheardriven pumping provided by the AFM tip-cantilever mechanism. The AFM detection of separated components is accomplished by exploiting the variation in the localized friction of the separated components. The AFM tip-cantilever provides the mechanism for producing shear-induced flows and rapid pumping. Shear-driven chromatography (SDC) is a relatively new concept that overcomes the speed and miniaturization limitations of conventional liquid chromatography. SDC is based on a sliding plate system, consisting of two flat surfaces, one of which has a recessed channel. A fluid flow is produced by axially sliding one plate past another, where the fluid has mechanical shear forces imposed at each point along the channel length. The shear-induced flow rates are very reproducible, and do not have pressure or voltage gradient limitations. SDC opens up a new range of enhanced separation kinetics by permitting the sample confinement with submicron dimensions. Small, highly confined liquid is advantageous for chromatographic separation because the separation rate is known to scale according to the square of the confined sample diameter. In addition, because shear-driven flows are not limited by fluid velocity, shear-driven liquid chromatography may provide up to 100,000 plate efficiency.

Anderson, Mark S.↗

NASA NTRS1956

Investigation of Heat Transfer From

The convective heat transfer from the surface of an ellipsoidal forebody of fineness ratio 3 and 20-inch maximum diameter was investigated in clear air for both stationary and rotating operation over a range of conditions including air speeds up to 240 knots, rotational speeds up to 1200 rpm, and angles of attack of 0 deg, 3 deg, and 6 deg. The results are presented in the form of heat-transfer coefficients and the correlation of Nusselt and Reynolds numbers. Both a uniform surface temperature and a uniform input heater density distribution were used. The experimental results agree well with theoretical predictions for uniform surface temperature distribution. Complete agreement was not obtained with uniform input heat density in the laminar-flow region because of conduction effects. No significant effects of rotation were obtained over the range of airstream and rotational speeds investigated. Operation at angle of attack had only minor effects on the local heat transfer. Transition from laminar to turbulent heat transfer occurred over a wide range of Reynolds numbers. The location of transition depended primarily on surface roughness and pressure and temperature gradients. Limited transient heating data indicate that the variation of surface temperature with time followed closely an exponential relation.

Lewis, James P.↗

NASA NTRS1988

Explicit upwind algorithm for the parabolized Navier-Stokes equations

A new explicit upwind algorithm based on Roe's flux-difference splitting (FDS) method has been developed for the three-dimensional Parabolized Navier-Stokes (PNS) equations. For three-dimensional flows, FDS's are determined separately for the two nonmarching directions and modified to account for the calculated shock angle in the crossflow plane. Second-order FDS is applied to the pressure and convection terms with the streamwise pressure gradient limited in the subsonic region to maintain a hyperbolic inviscid equation set. Second-order central differencing is obtained in the two-step algorithm for the shear and heat flux terms. The new algorithm is demonstrated for three laminar flow test cases: supersonic flow over a flat plate, hypersonic flow over a 15 deg ramp, and hypersonic flow past a 10 deg cone at a 24 deg angle of attack. The computed results agree well with experimental measurements.

Korte, John J.↗

NASA NTRS1994

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells are created using polygon-clipping algorithms. The grid is stored in a binary-tree structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded: a gradient-limited, linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The more robust of a series of viscous flux functions is used to provide the viscous fluxes at the cell interfaces. Adaptively-refined solutions of the Navier-Stokes equations using the Cartesian, cell-based approach are obtained and compared to theory, experiment, and other accepted computational results for a series of low and moderate Reynolds number flows.

Coirier, William J.↗

NASA NTRS1995

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells are created using polygon-clipping algorithms. The grid is stored in a binary-tree data structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded: A gradient-limited, linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The more robust of a series of viscous flux functions is used to provide the viscous fluxes at the cell interfaces. Adaptively-refined solutions of the Navier-Stokes equations using the Cartesian, cell-based approach are obtained and compared to theory, experiment and other accepted computational results for a series of low and moderate Reynolds number flows.

Coirier, William J.↗

NASA NTRS1988

On the origin of steep edges and filaments in vorticity and potential vorticity fields

High-resolution numerical calculations are shown which capture the fundamental process responsible for the intensification of vorticity gradients in an isolated vortex subject to externally imposed distrubances. Imposition of almost any weak strain or shear field in stripping away the relatively weak vorticity at the edge of the vortex and leaves it with gradients four to six orders of magnitude greater than in the initial state. Calculations displaying such enormous gradients have never been reported previously, because of the artificial eddy diffusivities that always limit such gradients in standard numerical models. The present calculations, which have no such limitations, have been made possible by the development of a novel and robust new numerical technique for vortex dynamics called contour surgery.

Dritschel, D. G.↗

NASA NTRS1999

Ames Optimized TCA Configuration

Configuration design at Ames was carried out with the SYN87-SB (single block) Euler code using a 193 x 49 x 65 C-H grid. The Euler solver is coupled to the constrained (NPSOL) and the unconstrained (QNMDIF) optimization packages. Since the single block grid is able to model only wing-body configurations, the nacelle/diverter effects were included in the optimization process by SYN87's option to superimpose the nacelle/diverter interference pressures on the wing. These interference pressures were calculated using the AIRPLANE code. AIRPLANE is an Euler solver that uses a unstructured tetrahedral mesh and is capable of computations about arbitrary complete configurations. In addition, the buoyancy effects of the nacelle/diverters were also included in the design process by imposing the pressure field obtained during the design process onto the triangulated surfaces of the nacelle/diverter mesh generated by AIRPLANE. The interference pressures and nacelle buoyancy effects are added to the final forces after each flow field calculation. Full details of the (recently enhanced) ghost nacelle capability are given in a related talk. The pseudo nacelle corrections were greatly improved during this design cycle. During the Ref H and Cycle 1 design activities, the nacelles were only translated and pitched. In the cycle 2 design effort the nacelles can translate vertically, and pitch to accommodate the changes in the lower surface geometry. The diverter heights (between their leading and trailing edges) were modified during design as the shape of the lower wing changed, with the drag of the diverter changing accordingly. Both adjoint and finite difference gradients were used during optimization. The adjoint-based gradients were found to give good direction in the design space for configurations near the starting point, but as the design approached a minimum, the finite difference gradients were found to be more accurate. Use of finite difference gradients was limited by the CPU time limit available on the Cray machines. A typical optimization run using finite difference gradients can use only 30 to 40 design variables and one optimization iteration within the 8 hour queue limit for the chosen grid size and convergence level. The efficiency afforded by the adjoint method allowed for 50-120 design variables and 5-10 optimization iterations in the 8 hour queue. Geometric perturbations to the wing and fuselage were made using the Hicks/Henne (HH) shape functions. The HH functions were distributed uniformly along the chords of the wing defining sections and lofted linearly. During single-surface design, constraints on thickness and volume at selected wing stations were imposed. Both fuselage camber and cross-sectional area distributions were permitted to change during design. The major disadvantage to the use of these functions is the inherent surface waviness produced by repeated use of such functions. Many smoothing operations were required following optimization runs to produce a configuration with reasonable smoothness. Wagner functions were also used on the wing sections but were never used on the fuselage. The Wagner functions are a family of increasingly oscillatory functions that have also been used extensively in airfoil design. The leading and trailing edge regions of the wing were designed by use of polynomial and monomial functions respectively. Twist was attempted but was abandoned because of little performance improvement available from changing the baseline twist.

Cliff, Susan E.↗

NASA NTRS1994

Multidisciplinary design optimization using genetic algorithms

Multidisciplinary design optimization (MDO) is an important step in the conceptual design and evaluation of launch vehicles since it can have a significant impact on performance and life cycle cost. The objective is to search the system design space to determine values of design variables that optimize the performance characteristic subject to system constraints. Gradient-based optimization routines have been used extensively for aerospace design optimization. However, one limitation of gradient based optimizers is their need for gradient information. Therefore, design problems which include discrete variables can not be studied. Such problems are common in launch vehicle design. For example, the number of engines and material choices must be integer values or assume only a few discrete values. In this study, genetic algorithms are investigated as an approach to MDO problems involving discrete variables and discontinuous domains. Optimization by genetic algorithms (GA) uses a search procedure which is fundamentally different from those gradient based methods. Genetic algorithms seek to find good solutions in an efficient and timely manner rather than finding the best solution. GA are designed to mimic evolutionary selection. A population of candidate designs is evaluated at each iteration, and each individual's probability of reproduction (existence in the next generation) depends on its fitness value (related to the value of the objective function). Progress toward the optimum is achieved by the crossover and mutation operations. GA is attractive since it uses only objective function values in the search process, so gradient calculations are avoided. Hence, GA are able to deal with discrete variables. Studies report success in the use of GA for aircraft design optimization studies, trajectory analysis, space structure design and control systems design. In these studies reliable convergence was achieved, but the number of function evaluations was large compared with efficient gradient methods. Applicaiton of GA is underway for a cost optimization study for a launch-vehicle fuel-tank and structural design of a wing. The strengths and limitations of GA for launch vehicle design optimization is studied.

Unal, Resit↗

NASA NTRS1993

On the accuracy of limiters and convergence to steady state solutions

This paper addresses the practical problem of obtaining convergence to steady state solutions when limiters are used in conjunction with upwind schemes on unstructured grids. The base scheme forms a gradient and limits it by imposing monotonicity conditions in the reconstruction stage. It is shown by analysis in one dimension that such an approach leads to various schemes meeting TVD requirements in one dimension. It is further shown that these formally second order accurate schemes are less than second order accurate in practice because of the action of the limiter function in smooth regions of the solution. Modifications are proposed to the limiter that restore the second order accuracy. In multiple dimensions these schemes produce steady state solutions that are monotone and devoid of oscillations. However, convergence stalls after a few orders of reduction in the residual. With the modified limiter, on the other hand, it is shown that converged steady state solutions can be obtained.

Venkatakrishnan, V.↗

NASA NTRS2008

Spacecraft Design Thermal Control Subsystem

The Thermal Control Subsystem engineers task is to maintain the temperature of all spacecraft components, subsystems, and the total flight system within specified limits for all flight modes from launch to end-of-mission. In some cases, specific stability and gradient temperature limits will be imposed on flight system elements. The Thermal Control Subsystem of "normal" flight systems, the mass, power, control, and sensing systems mass and power requirements are below 10% of the total flight system resources. In general the thermal control subsystem engineer is involved in all other flight subsystem designs.

Jupiter Icy Moon Orbiter (JIMO)↗

NASA NTRS2005

Development of Detectability Limits for On-Orbit Inspection of Space Shuttle Wing Leading Edge

At the conclusion of the Columbia Accident Investigation, one of the recommendations of the Columbia Accident Investigation Board (CAIB) was that NASA develop and implement an inspection plan for the Reinforced Carbon-Carbon (RCC) system components of the Space Shuttle. To address these issues, a group of scientists and engineers at NASA Langley Research Center proposed the use of an IR camera to inspect the RCC. Any crack in an RCC panel changes the thermal resistance of the material in the direction perpendicular to the crack. The change in thermal resistance can be made visible by introducing a heat flow across the crack and using an IR camera to image the resulting surface temperature distribution. The temperature difference across the crack depends on the change in the thermal resistance, the length of the crack, the local thermal gradient, and the rate of radiation exchange with the environment. This paper describes how the authors derived the minimum thermal gradient detectability limits for a through crack in an RCC panel. This paper will also show, through the use of a transient, 3-dimensional, finite element model, that these minimum gradients naturally exist on-orbit. The results from the finite element model confirm that there are sufficient thermal gradient to detect a crack on 96% of the RCC leading edge.

Stephan, Ryan A.↗

NASA NTRS1986

Combustion theory for liquids with a free surface. 3: Special problems

Two special problems concerning the combustion of liquids with a free surface, i.e., flame quenching during the mixing of a burning liquid inside a container and liquid burnout from a porous layer, are analyzed using a quasi-one-dimensional model. The critical parameters corresponding to the quenching of a burning fluid with a free surface are determined. Determinations are also made of the limiting pressure gradients corresponding to the transition from the combustion mode where the liquid evaporates from the surface of a porous layer to the mode where the phase transition surface lies inside the porous layer.

Milkov, S. N.↗

NASA NTRS1992

Learning and optimization with cascaded VLSI neural network building-block chips

To demonstrate the versatility of the building-block approach, two neural network applications were implemented on cascaded analog VLSI chips. Weights were implemented using 7-b multiplying digital-to-analog converter (MDAC) synapse circuits, with 31 x 32 and 32 x 32 synapses per chip. A novel learning algorithm compatible with analog VLSI was applied to the two-input parity problem. The algorithm combines dynamically evolving architecture with limited gradient-descent backpropagation for efficient and versatile supervised learning. To implement the learning algorithm in hardware, synapse circuits were paralleled for additional quantization levels. The hardware-in-the-loop learning system allocated 2-5 hidden neurons for parity problems. Also, a 7 x 7 assignment problem was mapped onto a cascaded 64-neuron fully connected feedback network. In 100 randomly selected problems, the network found optimal or good solutions in most cases, with settling times in the range of 7-100 microseconds.

Duong, T.↗