Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical diffusion”

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 721 records · Page 40

Constraints on the thermal evolution of Venus inferred from Magellan data

The impact craters with diameters from 1.5 to 280 km compiled from Magellan observations indicate that the crater population on Venus has a completely spatially random distribution and the size/density distribution of craters with diameters greater than or equal to 35 km is consistent with a 'production' population with an age of 500 plus or minus 250 m.y. The similarity in size distribution from area to area indicates that the crater distribution is independent of crater size. Also, the forms of the modified craters are virtually identical to those of the pristine craters. These observations imply that Venus reset its cratering record by global resurfacing 500 m.y. ago, and resurfacing declined relatively fast. The fact that less than 40 percent of all craters have been modified and that the few volcanically embayed craters are located on localized tectonic regions indicate that only minor and localized volcanism and tectonism have occurred since the latest vigorous resurfacing event approximately 500 m.y. ago and the interior of Venus has been solid and possibly colder than Earth's. This is because the high-temperature lithosphere of Venus would facilitate upward ascending of mantle plumes and result in extensive volcanism if the venusian upper mantle were as hot as or hotter than Earth's. Therefore, the present surface morphology of Venus may provide useful constraints on the pattern of that vigorous convection, and possibly on the thermal state of the venusian mantle. We examine this possibility through numerical calculations of three-dimensional thermal convection models in a spherical shell with temperature- and pressure-dependent Newtonian viscosity, temperature-dependent thermal diffusivity, pressure-dependent thermal expansion coefficient, and time-dependent internal heat production rate solar magnitude.

Arkani-Hamed, Jafar↗

A numerical treatment of geodynamic viscous flow problems involving the advection of material interfaces

Effective numerical treatment of multicomponent viscous flow problems involving the advection of sharp interfaces between materials of differing physical properties requires correction techniques to prevent spurious diffusion and dispersion. We develop a particular algorithm, based on modern shock-capture techniques, employing a two-step nonlinear method. The first step involves the global application of a high-order upwind scheme to a hyperbolic advection equation used to model the distribution of distinct material components in a flow field. The second step is corrective and involves the application of a global filter designed to remove dispersion errors that result from the advection of discontinuities (e.g., material interfaces) by high-order, minimally dissipative schemes. The filter introduces no additional diffusion error. Nonuniform viscosity across a material interface is allowed for by the implementation of a compositionally weighted-inverse interface viscosity scheme. The combined method approaches the optimal accuracy of modern shock-capture techniques with a minimal increase in computational time and memory. A key advantage of this method is its simplicity to incorporate into preexisting codes be they finite difference, element, or volume of two or three dimensions.

Lenardic, A.↗

Uncertain Representations of Sub-Grid Pollutant Transport in Chemistry-Transport Models and Impacts on Long-Range Transport and Global Composition

Sub-grid transport, by convection and turbulence, is known to play an important role in lofting pollutants from their source regions. Consequently, the long-range transport and climatology of simulated atmospheric composition are impacted. This study uses the Goddard Earth Observing System, Version 5 (GEOS-5) atmospheric model to study pollutant transport. The baseline model uses a Relaxed Arakawa-Schubert (RAS) scheme that represents convection through a sequence of linearly entraining cloud plumes characterized by unique detrainment levels. Thermodynamics, moisture and trace gases are transported in the same manner. Various approximate forms of trace-gas transport are implemented, in which the box-averaged cloud mass fluxes from RAS are used with different numerical approaches. Substantial impacts on forward-model simulations of CO (using a linearized chemistry) are evident. In particular, some aspects of simulations using a diffusive form of sub-grid transport bear more resemblance to space-biased CO observations than do the baseline simulations with RAS transport. Implications for transport in the real atmosphere will be discussed. Another issue of importance is that many adjoint/inversion computations use simplified representations of sub-grid transport that may be inconsistent with the forward models: implications will be discussed. Finally, simulations using a complex chemistry model in GEOS-5 (in place of the linearized CO model) are underway: noteworthy results from this simulation will be mentioned.

Pawson, Steven↗

Investigation of Unsteady Flow Interaction Between an Ultra-Compact Inlet and a Transonic Fan

In the present study, unsteady flow interaction between an ultra-compact inlet and a transonic fan stage is investigated. Future combat aircraft require ultra-compact inlet ducts as part of an integrated, advanced propulsion system to improve air vehicle capability and effectiveness to meet future mission needs. The main purpose of the study is to advance the current understanding of the flow interaction between two different ultra-compact inlets and a transonic fan for future design applications. Both URANS and LES approaches are used to calculate the unsteady flow field and are compared with the available measured data. The present study indicates that stall inception is mildly affected by the distortion pattern generated by the inlet with the current test set-up. The numerical study indicates that the inlet distortion pattern decays significantly before it reaches the fan face for the current configuration. Numerical results with a shorter distance between the inlet and fan show that counter-rotating vortices near the rotor tip due to the serpentine diffuser affects fan characteristics significantly.

Inlet↗

Refinement of Vortex Generators in a Streamline-Traced, External-Compression Supersonic Inlet

Computational simulations of the flow within a streamline-traced, external-compression supersonic inlet for Mach 1.664 without and with vortex generators were performed to refine the characterization of the inlet performance as measured by the total pressure recovery and the radial and circumferential total pressure distortion indices at the engine face. The refinement of the simulations concerned two aspects: 1) refinement of the grid for the simulations to evaluate and reduce uncertainty, and 2) refinement of the modeling and design of the vortex generators. The vortex generators studied were rectangular vane-type vortex generators arranged in co-rotating arrays within the subsonic diffuser. The vortex generator geometric factors of interest included the height, circumferential spacing, and angle-of-incidence. The flow through the inlet was simulated numerically through the solution of the steady-state, Reynolds-averaged Navier-Stokes equations on multi-block, structured grids using the Wind-US flow solver. The vortex generators were simulated using either a vortex generator model or with grids generated about each vortex generator. Statistical methods were used to compute confidence intervals and grid convergence indices to establish the uncertainties of the analyses with respect to grid refinement. Design-of-experiments methods were applied to quantify the effects of the geometric factors of the vortex generators. The analyses of the computed results illustrate the complexities of quantifying the uncertainties of the inlet performance and the implications of the uncertainties for the design of a vortex generator array for the STEX inlet.

Slater, John W.↗

Numerical Simulation of the High-Speed Leg of the National Transonic Facility

Numerical simulations for the flow inside the high-speed leg of the National Transonic Facility was conducted. The NASA Tetrahedral Unstructured Software System (TetrUSS) with its USM3D_ME solver was used to perform the numerical simulations. USM3D_ME is developed and maintained by NASA Langley Research Center. Simulations were conducted for three configurations: empty tunnel, body of revolution installed, and NASA Common Research Model installed in the test section. Simulations were performed for a test section Mach number of 0.7 and 0.85 and a corresponding Reynolds number of 8 million per ft. A controller was developed that automated dynamic outflow boundary and streamlined the process of running multiple simulations. The use of a dynamic outflow boundary was the key parameter to drive simulation to the desired tunnel conditions. The numerical simulations captured expected flow features in NTF test section and in the tunnel plenum. The simulations depicted that separation of the flow inside the diffuser appears to be asymmetric and more extensive toward the top and bottom walls. Analysis of the flow field were conducted, and the computed drag coefficient compared to wind tunnel data.

NTF↗

Numerical analysis of confined recirculating gas-solid turbulent flows

Particle dispersion in confined recirculating turbulent flows has been investigated numerically. The present computational scheme utilizes Eulerian forms of the governing equations and allows two-way coupling between gas and solid phases. A recently developed two-phase closure model based on the multiple-scale turbulence model was used for the estimation of turbulent viscosities and diffusivities. For the particle size and loading considered in this study, the effect of particulate phase on the conveying gas is small, the nonequilibrium between the two phases is significant. Effects of recirculation, expanded chamber size and secondary annular jet momentum on the particle mixing rate are also investigated. In general, the present numerical results are in reasonably good agreement with the available experimental data.

Chen, C. P.↗

STAN5: A program for numerical computation of two-dimensional internal and external boundary layer flows

A large variety of two dimensional flows can be accommodated by the program, including boundary layers on a flat plate, flow inside nozzles and diffusers (for a prescribed potential flow distribution), flow over axisymmetric bodies, and developing and fully developed flow inside circular pipes and flat ducts. The flows may be laminar or turbulent, and provision is made to handle transition.

Crawford, M. E.↗

Time-dependent Computational Studies of Premixed Flames in Microgravity

This report describes the research performed at the Center for Reactive Flow and Dynamical Systems in the Laboratory for Computational Physics and Fluid Dynamics, at the Naval Research Laboratory, in support of NASA Microgravity Science and Applications Program. The primary focus of this research is on investigating fundamental questions concerning the propagation and extinction of premixed flames in earth gravity and in microgravity environments. Our approach is to use detailed time-dependent, multispecies, numerical models as tools to simulate flames in different gravity environments. The models include a detailed chemical kinetics mechanism consisting of elementary reactions among the eight reactive species involved in hydrogen combustion, coupled to algorithms for convection, thermal conduction, viscosity, molecular and thermal diffusion, and external forces. The external force, gravity, can be put in any direction relative to flame propagation and can have a range of values. Recently more advanced wall boundary conditions such as isothermal and no-slip have been added to the model. This enables the simulation of flames propagating in more practical systems than before. We have used the numerical simulations to investigate the effects of heat losses and buoyancy forces on the structure and stability of flames, to help resolve fundamental questions on the existence of flammability limits when there are no external losses or buoyancy forces in the system, to understand the interaction between the various processes leading to flame instabilities and extinguishment, and to study the dynamics of cell formation and splitting. Our studies have been able to bring out the differences between upward- and downward-propagating flames and predict the zero-gravity behavior of these flames. The simulations have also highlighted the dominant role of wall heat losses in the case of downward-propagating flames. The simulations have been able to qualitatively predict the formation of multiple cells and the cessation of cell-splitting. Our studies have also shown that some flames in a microgravity environment can be extinguished due to a chemical instability and without any external losses. However, further simulations are needed to more completely understand upward-propagating and zero-gravity flames as well as to understand the potential effect of radiative heat losses.

Kailasanath, K.↗

Dynamic response of shock waves in transonic diffuser and supersonic inlet - An analysis with the Navier-Stokes equations and adaptive grid

An existing method which solves the multi-dimensional ensemble-averaged compressible time-dependent Navier-Stokes equations in conjunction with mixing length turbulence model and shock capturing technique has been extended to include the shock-tracking adaptive grid systems. The numerical scheme for solving the governing equations is based on a linearized block implicit approach. The effects of grid-motion and grid-distribution on the calculated flow solutions have been studied in relative detail and this is carried out in the context of physically steady, shocked flows computed with non-stationary grids. Subsequently, the unsteady dynamics of the flows occurring in a supercritically operated transonic diffuser and a mixed compression supersonic inlet have been investigated with the adaptive grid systems by solving the Navier-Stokes equations.

Liu, N.-S.↗

A pseudo-spectral collocation method applied to the problem of convective diffusive transport in fluids subject to unsteady residual accelerations

The problem of determining the sensitivity of Bridgman-Stockbarger directional solidification experiments to residual accelerations of the type associated with spacecraft in low earth orbit is analyzed numerically using a pseudo-spectral collocation method. The approach employs a novel iterative scheme combining the method of artificial compressibility and a generalized ADI method. The results emphasize the importance of the consideration of residual accelerations and careful selection of the operating conditions in order to take full advantages of the low gravity conditions.

Alexander, J. Iwan↗

Application of the Space-Time Conservation Element and Solution Element Method to One-Dimensional Convection-Diffusion Problems

In the space-time conservation element and solution element (CE/SE) method, the independent marching variables used comprise not only the mesh value of the physical dependent variables but also, in contrast to it typical numerical method, the Mesh values of the spatial derivatives of the physical variables The use of the extra marching variables results from the need to construct the two-level explicit and nondissipative schemes which are at the core of the CE/SE development. It also results from the need to minimize the stencil while maintaining accuracy. In this paper using the 1D(sub (alpha)-mu) scheme as an example, the effect of this added complication on consistency, accuracy and operation count is assessed. As part of this effort, an equivalent yet more efficient form of the alpha-mu scheme in which the independent marching variables are the local fluxes tied to each mesh point is introduced. Also, the intriguing relations that exist among the alpha-mu. Leapfrog, and DuFort-Frankel schemes are further explored. In addition, the redundance of the Leapfrog, DUFort-Frankel, and Lax scheme and the remedy for this redundance are discussed. This paper is concluded with the construction and evaluation of a CE/SE solver for the inviscid Burger equation.

Chang, Sin-Chung↗

Stability aspects of diverging subsonic flow

The paper describes the computation of two-dimensional, subsonic, diverging internal flows and how they differ from the corresponding converging flows. Such diverging or decelerating flows occur in such obvious places as subsonic diffusers and inlets; however, such flows also occur in supersonic nozzles in the presence of a normal shock. The flow instability and its relation to the numerical method used, boundary conditions, and viscous effects are assessed both analytically and numerically. The inviscid flow is shown to be physically unstable and a poor representation of the true viscous flow.

Cline, M. C.↗

Numerical simulation of a supersonic reacting mixing layer

In order to arrive at physical models that can adequately describe supersonic combustion, and develop accurate and efficient numerical techniques for the solution of such models' governing equations, a computer program has been developed for the study of reacting flows which considers the multicomponent diffusion and convection of important chemical species, as well as their finite state reaction and the interaction of the fluid mechanics and the chemistry that occurred. The code employs a hybrid Chebyshev pseudospectral technique for integration of the models' resulting governing equations; the program is here used to study a spatially developing and reacting mixing layer.

Drummond, J. Philip↗

Bipolar supernova remnants and the obliquity dependence of shock acceleration

The diffusive shock acceleration mechanism proposed to explain the bipolarity observed in the synchrotron radio emission of young adiabatically expanding shell SNRs is investigated by means of numerical simulations. The theoretical basis of the SNR models and the numerical computation methods are explained, and the results are presented in graphs and synthetic radio maps and discussed in detail. It is found that the efficiency of the acceleration process depends on the obliquity angle theta(Bn) between the shock normal and the uniform magnetic field: models with theta(Bn) of about 90 deg can reproduce the observed azimuthal intensity ratios in most cases, but models with theta(Bn) near 0 deg cannot.

Fulbright, Michael S.↗

Phase-field model for isothermal phase transitions in binary alloys

A new phase field model is described which models isothermal phase transitions between ideal binary alloy solution phases. Equations are developed for the temporal and spatial variation of the phase field, which describes the identity of the phase, and of the composition. An asymptotic analysis, as the gradient energy coefficient of the phase field becomes small, was conducted. From the analysis, it is shown that the model recovers classical sharp interface models of this situation when the interfacial layers are thin, and they show how to relate the parameters appearing in the phase field model to material and growth parameters in real systems. Further, three stages of temporal evolution are identified: the first corresponding to interfacial genesis which occurs very rapidly; the second to interfacial motion controlled by the local energy difference across the interface and diffusion; the last taking place on a long time scale in which curvature effects are important and which correspond to Ostwald ripening. The results of the numerical calculations are presented.

Wheeler, A. A.↗

Flow effects in a vertical CVD reactor

A model is presented to simulate the non-Boussinesq flow in a vertical, two-dimensional, chemical vapor deposition reactor under atmospheric pressure. Temperature-dependent conductivity, mass diffusivity, viscosity models, and reactive species mass transfer to the substrate are incorporated. In the limits of small Mach number and small aspect ratio, asymptotic expressions for the flow, temperature, and species fields are developed. Soret diffusion effects are also investigated. Analytical solutions predict an inverse relationship between temperature field and concentration field due to Soret effects. This finding is consistent with numerical simulations, assisting in the understanding of the complex interactions amongst the flow, thermal, and species fields in a chemically reacting system.

Young, G. W.↗

Design of two-dimensional channels with prescribed velocity distributions along the channel walls

A general method of design is developed for two-dimensional unbranched channels with prescribed velocities as a function of arc length along the channel walls. The method is developed for both compressible and incompressible, irrotational, nonviscous flow and applies to the design of elbows, diffusers, nozzles, and so forth. In part I solutions are obtained by relaxation methods; in part II solutions are obtained by a Green's function. Five numerical examples are given in part I including three elbow designs with the same prescribed velocity as a function of arc length along the channel walls but with incompressible, linearized compressible, and compressible flow. One numerical example is presented in part II for an accelerating elbow with linearized compressible flow, and the time required for the solution by a Green's function in part II was considerably less than the time required for the same solution by relaxation methods in part I.

Stanitz, John D↗