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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Adjoint-Based Mesh Adaptation and Shape Optimization for Simulations with Propulsion

We demonstrate a well-posed formulation of permeable boundary conditions and mass- flow-rate functionals for adjoint-based mesh refinement and shape optimization governed by the steady Euler equations. The boundary conditions are used to model propulsion- system effects of inlets and nozzles. A two-shock diffuser with an analytic solution is used to verify the implementation. Numerical examples show that the adjoint solution is smooth at the boundary, indicating that the discretization is adjoint consistent when exit pressure is specified at subsonic outflow, and stagnation temperature and pressure at subsonic inflow. The results focus on improving simulation techniques for low-boom aircraft analysis and design. By including mass-flow-rate outputs, we obtain reliable estimates of engine flow rates concurrently with nearfield pressure signatures without increasing simulation cost. We also demonstrate the importance of mass-flow-rate constraints in shape optimization by examining trade-offs between maximizing performance of a shrouded supersonic nozzle and minimizing shocks in its nearfield.

Nemec, Marian↗
NASA NTRS1999

Melt Stabilization of PbSnTe in a Magnetic Field

Both the experimental observation and numerical simulation indicate that the Bridgman growth of PbSnTe under the microgravity environment in space is still greatly influenced by buoyancy-induced convection. The application of a magnetic field during the semiconductor growth can dampen the convective flow in the metal-like melt. However, for Bridgman growth of PbSnTe on earth (with either vertical or horizontal configuration), both experimental observation and numerical modeling suggest that even with a strong magnetic furnace (5-Tesla constant axial magnetic field), the convective flow in the melt still cannot be sufficiently suppressed to reach the diffusion-controlled level. In order to completely dampen the buoyancy-induced convection on earth, estimates based on scaling analysis indicate that for common experimental conditions, an extremely high magnetic field is required, far beyond the capacity of the experimental apparatus currently available. Therefore, it is proposed that only the combination of microgravity environment and magnetic damping will produce the desired diffusion-controlled growth state for this particular material. The primary objectives of this study are to provide a quantitative understanding of the complex transport phenomena during solidification of non-dilute binarys, to furnish a numerical tool for furnace design and growth condition optimization, to provide estimates of the required magnetic field strength for low gravity growth, and to assess the role of magnetic damping for space and earth control of the double-diffusive convection. As an integral part of a NASA research program, our numerical simulation supports both the flight and ground-based experiments in an effort to bring together a complete picture of the complex physical phenomena involved in the crystal growth process. For Bridgman growth of PbSnTe under microgravity (with both vertical and horizontal configurations), the simulations suggest that a moderate axial magnetic field of only a few kilo-Gauss in strength could effectively eliminate buoyancy-induced convection in the melt and control solute segregation. Therefore, this work confirms the idea that the combination of microgravity environment and the magnetic damping will indeed be sufficient to produce the desired diffusion-controlled growth state for PbSnTe.

Fripp, Archibald L.↗
NASA NTRS1994

New Developments in the Method of Space-Time Conservation Element and Solution Element-Applications to Two-Dimensional Time-Marching Problems

A new numerical discretization method for solving conservation laws is being developed. This new approach differs substantially in both concept and methodology from the well-established methods, i.e., finite difference, finite volume, finite element, and spectral methods. It is motivated by several important physical/numerical considerations and designed to avoid several key limitations of the above traditional methods. As a result of the above considerations, a set of key principles for the design of numerical schemes was put forth in a previous report. These principles were used to construct several numerical schemes that model a 1-D time-dependent convection-diffusion equation. These schemes were then extended to solve the time-dependent Euler and Navier-Stokes equations of a perfect gas. It was shown that the above schemes compared favorably with the traditional schemes in simplicity, generality, and accuracy. In this report, the 2-D versions of the above schemes, except the Navier-Stokes solver, are constructed using the same set of design principles. Their constructions are simplified greatly by the use of a nontraditional space-time mesh. Its use results in the simplest stencil possible, i.e., a tetrahedron in a 3-D space-time with a vertex at the upper time level and other three at the lower time level. Because of the similarity in their design, each of the present 2-D solvers virtually shares with its 1-D counterpart the same fundamental characteristics. Moreover, it is shown that the present Euler solver is capable of generating highly accurate solutions for a famous 2-D shock reflection problem. Specifically, both the incident and the reflected shocks can be resolved by a single data point without the presence of numerical oscillations near the discontinuity.

Chang, Sin-Chung↗
NASA NTRS1993

Numerical simulations of three-dimensional laminar flow over a backward facing step; flow near side walls

This paper reports the results of numerical simulations of steady, laminar flow over a backward-facing step. The governing equations used in the simulations are the full 'compressible' Navier-Stokes equations, solutions to which were computed by using a cell-centered, finite volume discretization. The convection terms of the governing equations were discretized by using the Advection Upwind Splitting Method (AUSM), whereas the diffusion terms were discretized using central differencing formulas. The validity and accuracy of the numerical solutions were verified by comparing the results to existing experimental data for flow at identical Reynolds numbers in the same back step geometry. The paper focuses attention on the details of the flow field near the side wall of the geometry.

Steinthorsson, Erlendur↗
NASA NTRS1984

Differential rotation in a solar-driven quasi-axisymmetric circulation

The concept of a quasi-axisymmetric circulation is used to explore the global scale dynamics of planetary atmospheres. A numerical circulation model applicable to Jupiter is presented, and an analytical study is performed elucidating the conditions leading to differential rotation in an atmosphere which is convectively unstable. A linear system forced by solar differential heating is considered, with nonlinear effects arising from advection being represented in the form of eddy diffusion. An empirical, latitudinal spectrum of the observed zonal wind field on Jupiter is discussed. Numerical solutions are presented which reveal banded wind fields with alternating and equatorial zonal jets and a multicellular Ferrel-Thomson meridional circulation consistent with the observed cloud striations on Jupiter. The vertical derivatives are parameterized to construct a simplified one-layer model.

Mayr, H. G.↗
NASA NTRS1986

Transitions to chaos in two-dimensional double-diffusive convection

The partial differential equations governing two-dimensional thermosolutal convection in a Boussinesq fluid with free boundary conditions have been solved numerically in a regime where oscillatory solutions can be found. A systematic study of the transition from nonlinear periodic oscillations to temporal chaos has revealed sequences of period-doubling bifurcations. Overstability occurs if the ratio of the solutal to the thermal diffusivity tau is less than 1 and the solutal Rayleigh number Rs is sufficiently large. Solutions have been obtained for two representative values of tau. For tau = 0.316, R(s) = 10,000, symmetrical oscillations undergo a bifurcation to asymmetry, followed by a cascade of period-doubling bifurcations leading to aperiodicity, as the thermal Rayleigh number R(T) is increased. At higher values of R(T), the bifurcation sequence is repeated in reverse, restoring simple periodic solutions. As R(T) is further increased more period-doubling cascades, followed by chaos, can be identified. Within the chaotic regions there are narrow periodic windows, and multiple branches of oscillatory solutions coexist. Eventually the oscillatory branch ends and only steady solutions can be found. The development of chaos has been investigated for tau = 0.1 by varying R(T) for several different values of R(s). When R(s) is sufficiently small there are periodic solutions whose period becomes infinite at the end of the oscillatory branch. As R(s) is increased, chaos appears in the neighborhood of these heteroclinic orbits. At higher values of R(s), chaos is found for a broader range in R(T). A truncated fifth-order model suggest that the appearance of chaos is associated with heteroclinic bifurcations.

Knobloch, E.↗
NASA NTRS1991

Ignition and structure of a laminar diffusion flame in a compressible mixing layer with finite rate chemistry

The ignition and structure of a reacting compressible mixing layer is considered using finite rate chemistry lying between two streams of reactants with different freestream speeds and temperatures. Numerical integration of the governing equations show that the structure of the reacting flow can be quite complicated depending on the magnitude of the Zeldovich number. An analysis of both the ignition a diffusion flame regimes is presented using a combination of large Zeldovich number asymptotics and numerics. This allows to analyze the behavior of these regimes as a function of the parameters of the problem.

Grosch, C. E.↗
NASA NTRS1991

Ignition and structure of a laminar diffusion flame in a compressible mixing layer with finite rate chemistry

The ignition and structure of a reacting compressible mixing layer is considered using finite rate chemistry lying between two streams of reactants with different freestream speeds and temperatures. Numerical integration of the governing equations show that the structure of the reacting flow can be quite complicated depending on the magnitude of the Zeldovich number. An analysis of both the ignition and diffusion flame regimes is presented using a combination of large Zeldovich number asymptotics and numerics. This allows to analyze the behavior of these regimes as a function of the parameters of the problem.

Grosch, C. E.↗
NASA NTRS1972

A generalized theory for the design of contraction cones and other low speed ducts

A generalization of the Tsien method of contraction cone design is described. The design velocity distribution is expressed in such a form that the required high order derivatives can be obtained by recursion rather than by numerical or analytic differentiation. The method is applicable to the design of diffusers and converging-diverging ducts as well as contraction cones. The computer program is described and a FORTRAN listing of the program is provided.

Barger, R. L.↗
NASA NTRS1975

Numerical solution of the Navier-Stokes equations with topography

A finite difference scheme for solving the equations of fluid motion in a generalized coordinate system has been constructed. The scheme conserves mass and all the first integral moments of the motion. The scheme also advectively 'almost conserves' second moments, in that the magnitude of implicit numerical smoothing is typically about an order smaller than explicit viscosity and diffusion. Calculations with the model support the theoretical conjecture that the difference scheme is stable whenever the analogous Cartesian scheme is stable. The scheme has been used to calculate dry atmospheric convection due to differential heating between top and bottom of mountainous terrain. The general small-scale characteristics of mountain up-slope winds have been simulated. In addition, the results have demonstrated the crucial role played by the eddy diffusivities and the environmental stability, in determining both the quantitative and the qualitative features of the circulation.

Gal-Chen, T.↗
NASA NTRS1982

Dissipation of ionospheric irregularities by wave-particle and collisional interactions

The nonlinear dissipation of plasma irregularities aligned parallel to an ambient magnetic field is studied numerically using a model which employs both wave-particle and collisional diffusion. A wave-particle diffusion coefficient derived from a local theory of the universal drift instability is used. This coefficient is effective in regions of nonzero plasma gradients and produces triangular-shaped irregularities with spectra which vary as f to the -4th, where f is the spatial frequency. Collisional diffusion acts rapidly on the vertices of the irregularities to reduce their amplitude. The simultaneous action of the two dissipative processes is more efficient than collisions acting alone. In this model, wave-particle diffusion mimics the forward cascade process of wave-wave coupling.

Bernhardt, P. A.↗
NASA NTRS1982

Magnetograph response to canopy-type fields

The response of longitudinal-field magnetographs to magnetic fields which are semi-infinite or confined to a horizontal layer is discussed with respect to the interpretation of solar diffuse fields, observed towards the limb, in terms of magnetic canopy models. Numerical results are presented for several reference solar models and typical 'calibration' curves are shown for the CI 9111 A, Fe I 8688 A, and Ca II 8542 A lines in magnetostatic atmospheres derived from a mean model. A procedure is developed for determining the base heights of magnetic canopies from observations with an uncertainty not exceeding the order of a pressure scale height. Until definitive information regarding atmospheric structure inside flux tubes can be developed from theory or observation, reliable field strengths cannot be derived from the data.

Jones, H. P.↗
NASA NTRS1986

Applications of a new wall function to turbulent flow computations

A new wall function approach is developed based on a wall law suitable for incompressible turbulent boundary layers under strong adverse pressure gradients. This wall law was derived from a one-dimensional analysis of the turbulent kinetic energy equation with gradient diffusion concept employed in modeling the near-wall shear stress gradient. Numerical testing cases for the present wall functions include turbulent separating flows around an airfoil and turbulent recirculating flows in several confined regions. Improvements on the predictions using the present wall functions are illustrated. For cases of internal recirculating flows, one modification factor for improving the performance of the k-epsilon turbulence model in the flow recirculation regions is also included.

Chen, Y. S.↗
NASA NTRS1988

Effect of heat release on the spatial stability of a supersonic reacting mixing layer

A numerical study of the stability of compressible mixing layers in which a diffusion flame is embedded is described. The mean velocity profile has been approximated by a hyperbolic tangent profile and the limit of infinite activation energy taken, which reduces the diffusion flame to a flame sheet. The addition of combustion in the form of a flame sheet was found to have important, and complex, effects on the flow stability.

Jackson, T. L.↗
NASA NTRS1989

Ballistic transport in planetary ring systems due to particle erosion mechanisms. I - Theory, numerical methods, and illustrative examples

Ballistic transport, defined as the net radial transport of mass and angular momentum due to exchanges of meteoroid hypersonic-impact ejecta by neighboring planetary ring regions on time-scales orders-of-magnitude shorter than the age of the solar system, is presently considered as a problem in mathematical physics. The preliminary results of a numerical scheme for following the combined effects of ballistic transport and viscous diffusion demonstrate that ballistic transport generates structure near sharp edges already present in the ring-mass distribution; the entire ring system ultimately develops an undulatory structure whose length scale is typically of the order of the radial excursion of the impact ejecta.

Durisen, Richard H.↗
NASA NTRS1990

Supersonic combustor modeling

The physical phenomena involved when a supersonic flow undergoes chemical reaction are discussed. Detailed physical models of convective and diffusive mixing, and finite rate chemical reaction in supersonic flow are presented. Numerical algorithms used to solve the equations governing these processes are introduced. Computer programs using these algorithms are used to analyze the structure of the reacting mixing layer. It is concluded that, as in subsonic flow, exothermic heat release in unconfined supersonic flows retards fuel/air mixing. Non mixing is shown to be a potential problem in reducing the efficiency of supersonic as well as subsonic combustion. Techniques for enhancing fuel/air mixing and combustion are described.

Riggins, David W.↗
NASA NTRS1991

Multigrid for hypersonic viscous two- and three-dimensional flows

The use of a multigrid method with central differencing to solve the Navier-Stokes equations for hypersonic flows is considered. The time dependent form of the equations is integrated with an explicit Runge-Kutta scheme accelerated by local time stepping and implicit residual smoothing. Variable coefficients are developed for the implicit process that removes the diffusion limit on the time step, producing significant improvement in convergence. A numerical dissipation formulation that provides good shock capturing capability for hypersonic flows is presented. This formulation is shown to be a crucial aspect of the multigrid method. Solutions are given for two-dimensional viscous flow over a NACA 0012 airfoil and three-dimensional flow over a blunt biconic.

Turkel, E.↗
NASA NTRS1991

Multigrid for hypersonic viscous two- and three-dimensional flows

The use of a multigrid method with central differencing to solve the Navier-Stokes equations for hypersonic flows is considered. The time-dependent form of the equations is integrated with an explicit Runge-Kutta scheme accelerated by local time stepping and implicit residual smoothing. Variable coefficients are developed for the implicit process that remove the diffusion limit on the time step, producing significant improvement in convergence. A numerical dissipation formulation that provides good shock-capturing capability for hypersonic flows is presented. This formulation is shown to be a crucial aspect of the multigrid method. Solutions are given for two-dimensional viscous flow over a NACA 0012 airfoil and three-dimensional viscous flow over a blunt biconic.

Turkel, E.↗