Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion equation”

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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

An explicit Runge-Kutta method for turbulent reacting flow calculations

The paper presents a numerical method for the solution of the conservation equations governing steady, reacting, turbulent viscous flow in two-dimensional geometries, in both Cartesian and axisymmetric coordinates. These equations are written in Favre-averaged form and closed with a first order model. A two-equation K-epsilon model, where low Reynolds number and compressibility effects are included, and a modified eddy-break up model are used to simulate fluid mechanics turbulence, chemistry and turbulence-combustion interaction. The solution is obtained by using a pseudo-unsteady method with improved perturbation propagation properties. The equations are discretized in space by using a finite volume formulation. An explicit multi-stage dissipative Runge-Kutta algorithm is then used to advance the flow equations in the pseudo-time. The method is applied to the computation of both diffusion and premixed turbulent reacting flows. The computed temperature distributions compare favorably with experimental data.

Boretti, A. A.↗
NASA NTRS1991

Computations of the three-dimensional flow and heat transfer within a coolant passage of a radial turbine blade

A numerical code is developed for computing three-dimensional, turbulent, compressible flow within coolant passages of turbine blades. The code is based on a formulation of the compressible Navier-Stokes equations in a rotating frame of reference in which the velocity dependent variable is specified with respect to the rotating frame instead of the inertial frame. The algorithm employed to obtain solutions to the governing equation is a finite-volume LU algorithm that allows convection, source, as well as diffusion terms to be treated implicitly. In this study, all convection terms are upwind differenced by using flux-vector splitting, and all diffusion terms are centrally differenced. This paper describes the formulation and algorithm employed in the code. Some computed solutions for the flow within a coolant passage of a radial turbine are also presented.

Shih, T. I.-P.↗
NASA NTRS1991

Cosmic-ray heating of cooling flows - A critical analysis

It is shown that a combination of MHD wave-mediated cosmic ray heating and thermal conduction could balance cooling in intracluster media and substantially reduce the rate of inflow. The appropriate system of steady state equations is solved, including a new self-consistent formulation for the cosmic-ray diffusivity. Models which can produce substantial positive temperature gradients in static configurations are found when conduction is reduced by a factor of 10 or more. These models have too-flat thermal pressure profiles compared with observations. It is found that cosmic-ray heating is unlikely either to stabilize positive density perturbations against condensation or to contribute appreciably to the powering of the optical filaments.

Loewenstein, Michael↗
NASA NTRS1993

A numerical analysis of high-temperature heat pipe startup from the frozen state

Continuum and rarefied vapor flows co-exist along the heat pipe length for most of the startup period. A two-region model is proposed in which the vapor flow in the continuum region is modeled by the compressible Navier-Stokes equations, and the vapor flow in the rarefied region is simulated by a self-diffusion model. The two vapor regions are linked with appropriate boundary conditions, and heat pipe wail, wick, and vapor flow are solved as a conjugate problem. The numerical solutions for the entire heat pipe startup process from the frozen state are compared with the corresponding experimental data with good agreement.

Cao, Y.↗
NASA NTRS1995

Boundary layer blockage in expansion tube nozzles

The results of a first order perfect gas correction for the effects of the boundary layer formation within expansion tubes with nozzles are presented. The analytical model developed to describe the boundary layer formation within the expansion tube and an expansion nozzle located at the end of the acceleration tube is based on the Karman integral equations. The results of this analytical model are compared with experimental data from an expansion diffuser. The model provides a useful tool for the preliminary design of nozzles for such facilities.

Sudnitsin, Olga↗
NASA NTRS1996

Self-Consistent Simulation of the Brownian Stage of Dust Growth

It is a widely accepted view that in proto-planetary accretion disks the collision and following sticking of dust particles embedded in the gas eventually leads to the formation of planetesimals (coagulation). For the smallest dust grains, Brownian motion is assumed to be the dominant source of their relative velocities leading to collisions between these dust grains. As the dust grains grow they eventually couple to the turbulent motion of the gas which then drives the coagulation much more efficiently. Many numerical coagulation simulations have been carried out to calculate the fractal dimension of the aggregates, which determines the duration of the ineffective Brownian stage of growth. Predominantly on-lattice and off-lattice methods were used. However, both methods require simplification of the astrophysical conditions. The aggregates found by those methods had a fractal dimension of approximately 2 which is equivalent to a constant, mass-independent friction time. If this value were valid for the conditions in an accretion disk, this would mean that the coagulation process would finally 'freeze out' and the growth of a planetesimal would be impossible within the lifetime of an accretion disk. In order to investigate whether this fractal dimension is model independent, we simulate self-consistently the Brownian stage of the coagulation by an N-particle code. This method has the advantage that no further assumptions about homogeneity of the dust have to be made. In our model, the dust grains are considered as aggregates built up of spheres. The equation of motion of the dust grains is based on the probability density for the diffusive transport within the gas atmosphere. Because of the very low number density of the dust grains, only 2-body-collisions have to be considered. As the Brownian stage of growth is very inefficient, the system is to be simulated over long periods of time. In order to find close particle pairs of the system which are most likely to undergo a collision, we use a particle-in-cell (PIC) method for the early stages of the simulation where the system is still very homogeneous and a tree method later when the particles are more clustered.

Kempf, S.↗
NASA NTRS1998

Tape-Drop Transient Model for In-Situ Automated Tape Placement of Thermoplastic Composites

Composite parts of nonuniform thickness can be fabricated by in-situ automated tape placement (ATP) if the tape can be started and stopped at interior points of the part instead of always at its edges. This technique is termed start/stop-on-the-part, or, alternatively, tape-add/tape-drop. The resulting thermal transients need to be managed in order to achieve net shape and maintain uniform interlaminar weld strength and crystallinity. Starting-on-the-part has been treated previously. This paper continues the study with a thermal analysis of stopping-on-the-part. The thermal source is switched off when the trailing end of the tape enters the nip region of the laydown/consolidation head. The thermal transient is determined by a Fourier-Laplace transform solution of the two-dimensional, time-dependent thermal transport equation. This solution requires that the Peclet number Pe (the dimensionless ratio of inertial to diffusive heat transport) be independent of time and much greater than 1. Plotted isotherms show that the trailing tape-end cools more rapidly than the downstream portions of tape. This cooling can weaken the bond near the tape end; however the length of the affected region is found to be less than 2 mm. To achieve net shape, the consolidation head must continue to move after cut-off until the temperature on the weld interface decreases to the glass transition temperature. The time and elapsed distance for this condition to occur are computed for the Langley ATP robot applying PEEK/carbon fiber composite tape and for two upgrades in robot performance. The elapsed distance after cut-off ranges from about 1 mm for the present robot to about 1 cm for the second upgrade.

Costen, Robert C.↗
NASA NTRS2008

Innermost Van Allen Radiation Belt for High Energy Protons at Saturn

The high energy proton radiation belts of Saturn are energetically dominated by the source from cosmic ray albedo neutron decay (CRAND), trapping of protons from beta decay of neutrons emitted from galactic cosmic ray nuclear interactions with the main rings. These belts were originally discovered in wide gaps between the A-ring, Janus/Epimetheus, Mimas, and Enceladus. The narrow F and G rings significant affected the CRAND protons but did not produce total depletion. Voyager 2 measurements subsequently revealed an outermost CRAND proton belt beyond Enceladus. Although the source rate is small, the trapping times limited by radial magnetospheric diffusion are very long, about ten years at peak measured flux inwards of the G ring, so large fluxes can accumulate unless otherwise limited in the trapping region by neutral gas, dust, and ring body interactions. One proposed final extension of the Cassini Orbiter mission would place perikrone in a 3000-km gap between the inner D ring and the upper atmosphere of Saturn. Experience with CRAND in the Earth's inner Van Allen proton belt suggests that a similar innermost belt might be found in this comparably wide region at Saturn. Radial dependence of magnetospheric diffusion, proximity to the ring neutron source, and northward magnetic offset of Saturn's magnetic equator from the ring plane could potentially produce peak fluxes several orders of magnitude higher than previously measured outside the main rings. Even brief passes through such an intense environment of highly penetrating protons would be a significant concern for spacecraft operations and science observations. Actual fluxes are limited by losses in Saturn's exospheric gas and in a dust environment likely comparable to that of the known CRAND proton belts. The first numerical model of this unexplored radiation belt is presented to determine limits on peak magnitude and radial profile of the proton flux distribution.

Cooper, John F.↗
NASA NTRS1988

The role of gravity wave generated advection and diffusion in transport of tracers in the mesosphere

The roles of advection and turbulent diffusion in controlling tracer transport in the meridional plane of the mesosphere are compared using a version of the global primitive equation model for the middle atmosphere developed by Holton and Wehrbein (1980). The results of model sensitivity tests demonstrated that transport in the mesosphere for long-lived species with large-scale heights, such as H2O, CO, and NO(x), is primarily an advective process; on a global scale, small-scale turbulent diffusion produced by gravity wave breaking plays a secondary role. It is concluded that, for such tracers, the eddy diffusion coefficient used in one-dimensional chemical models of the mesosphere should be regarded as a parameterization of the transport by the meridional circulation, not as a parameterization of mixing by small-scale diffusion.

Holton, James R.↗
NASA NTRS1986

Thermofluid analysis of the SSME preburner using a gas-gas diffusion model for oxygen and hydrogen combustion at supercritical pressures

The paper discusses the thermofluid analysis of the Space Shuttle Main Engine (SSME) fuelside preburner. The governing equations have been solved numerically to predict flow, heat transfer, mixing, and combustion. A two-fluid approach is adopted in which oxygen is regarded as one fluid and hydrogen is regarded as the other fluid. The chemical kinetics is assumed to be very fast so that combustion is primarily controlled by the rate of mixing between oxygen and hydrogen. The preburner pressure is much greater than the critical pressures of oxygen and hydrogen; hence, a gas-gas diffusion model (rather than an evaporation model) has been developed to compute the rate of interphase mixing. Empirical correlations have been incorporated to account for the effect of slip on the interphase exchange. A sensitivity study has been performed with various model parameters. It is observed that the model can predict possibility of incomplete combustion and local regions of high temperatures under steady operating conditions. Some of these anomalies have been observed in actual tests, and the numerical model is useful for understanding possible causes and remedies. At least some measurements are needed for quantitative verification of the model.

Prakash, C.↗
NASA NTRS1988

A two-dimensional numerical simulation of a supersonic, chemically reacting mixing layer

Research has been undertaken to achieve an improved understanding of physical phenomena present when a supersonic flow undergoes chemical reaction. A detailed understanding of supersonic reacting flows is necessary to successfully develop advanced propulsion systems now planned for use late in this century and beyond. In order to explore such flows, a study was begun to create appropriate physical models for describing supersonic combustion, and to develop accurate and efficient numerical techniques for solving the governing equations that result from these models. From this work, two computer programs were written to study reacting flows. Both programs were constructed to consider the multicomponent diffusion and convection of important chemical species, the finite rate reaction of these species, and the resulting interaction of the fluid mechanics and the chemistry. The first program employed a finite difference scheme for integrating the governing equations, whereas the second used a hybrid Chebyshev pseudospectral technique for improved accuracy.

Drummond, J. Philip↗
DOE OSTI2025

From disorganized data to emergent dynamic models: Questionnaires to partial differential equations

Starting with sets of disorganized observations of spatially varying and temporally evolving systems, obtained at different (also disorganized) sets of parameters, we demonstrate the data-driven derivation of parameter dependent, evolutionary partial differential equation (PDE) models capable of generating the data. This tensor type of data is reminiscent of shuffled (multidimensional) puzzle tiles. The independent variables for the evolution equations (their “space” and “time”) as well as their effective parameters are all emergent , i.e. determined in a data-driven way from our disorganized observations of behavior in them. We use a diffusion map based questionnaire approach to build a smooth parametrization of our emergent space/time/parameter space for the data. This approach iteratively processes the data by successively observing them on the “space,” the “time” and the “parameter” axes of a tensor. Once the data become organized, we use machine learning (here, neural networks) to approximate the operators governing the evolution equations in this emergent space. Our illustrative examples are based (i) on a simple advection–diffusion model; (ii) on a previously developed vertex-plus-signaling model of Drosophila embryonic development; and (iii) on two complex dynamic network models (one neuronal and one coupled oscillator model) for which no obvious smooth embedding geometry is known a priori. This allows us to discuss features of the process like symmetry breaking, translational invariance, and autonomousness of the emergent PDE model, as well as its interpretability.

generative models↗
NASA NTRS1991

Application of a Reynolds Stress turbulence model to a supersonic hydrogen-air diffusion flame

A second-order differential Reynolds Stress turbulence model has been applied to the Favre-averaged Navier-Stokes equations for the study of supersonic flows undergoing hydrogen-air chemical reactions. An assumed Beta Probability Density Function is applied to account for the chemical source terms in the conservation equations. An algebraic Reynolds Flux model is used for the fluctuating density-velocity as well as the species mass fraction-velocity correlations. The variances of temperature and species fluctuations are also modelled using an algebraic flux technique. A seven-species, seven-reaction finite rate chemistry mechanism is used to simulate the combustion processes. The resulting formulation is validated by comparison with experimental data on reacting supersonic axisymmetric jets. Results obtained for specific conditions indicate that the effect of chemical reaction on the turbulence is significant.

Chandrasekhar, R.↗
NASA NTRS1993

Application of a Reynolds stress turbulence model to a supersonic radiating hydrogen-air diffusion flame

A second-order differential Reynolds Stress turbulence model has been applied to the Favre-averaged Navier-Stokes equations for the study of supersonic flows with finite-rate chemistry and radiation. An assumed Beta Probability Density Function is applied to account for the chemical source terms and the radiative flux terms in the conservation equations. A seven-species, seven-reaction finite rate chemistry mechanism is used to simulate the combustion process. The tangent slab approximation is used in radiative flux formulation. A pseudo-gray gas model is used to represent the absorption-emission characteristics of the participating species. The turbulence/radiation interaction is achieved via a new formulation. The resulting formulation is validated by comparison with experimental data on reacting supersonic axisymmetric jets. Results obtained for specific conditions indicate that the effect of chemical reaction on the turbulence is significant. Also, the radiative heat transfer is enhanced by the turbulence.

Chandrasekhar, R.↗
NASA NTRS2024

Detection of Extensive Equatorial Plasma Depletions after the 2022 Tongan Volcanic Eruption from Multiple Geodetic Satellite Ranging Systems

We present a number of unique observations of ionospheric anomalies following the Hunga-Tonga Hunga-Ha'apai (HTHH) volcanic eruption on 15 January 2022. All are based on non-dedicated geodetic satellite systems: Global Positioning System tracking of Low Earth Orbit (LEO) CubeSats, intersatellite tracking between two GRACE Follow-On satellites, satellite radar altimeters to the ocean surface, and Doppler radio beacons from ground stations to LEO geodetic satellites. Their observations revealed the development of anomalously large trough-like plasma depletions, along with plasma bubbles, in the equatorial regions of the Pacific and East Asian sectors. Trough-like plasma depletions appeared to be confined within approximately ±20° magnetic latitude, accompanied by density enhancements just outside this latitude range. These plasma depletions and enhancements were aligned with the magnetic equator and occurred across broad longitudes. They were detected in regions where atmospheric waves from the HTHH eruption passed through around the time of the sunset terminator. We interpret these phenomena in terms of the E dynamo electric fields driven by atmospheric waves from the eruption. The uplift of the ionosphere beyond satellite altitudes, followed by subsequent plasma diffusion to higher latitudes along magnetic field lines, results in the formation of trough-like plasma depletions around the magnetic equator and density enhancement at higher latitudes. The detection of plasma bubbles in the Asian sector during the non-bubble season (January) is likely associated with the uplift of the ionosphere at the sunset terminator.

Sentinel-3B↗
DOE OSTI2020

Rechargeable Alkali-Ion Battery Materials: Theory and Computation

Since its development in the 1970s, the rechargeable alkali-ion battery has proven to be a truly transformative technology, providing portable energy storage for devices ranging from small portable electronics to sizable electric vehicles. Here, we present a review of modern theoretical and computational approaches to the study and design of rechargeable alkali-ion battery materials. Additionally, starting from fundamental thermodynamics and kinetics phenomenological equations, we rigorously derive the theoretical relationships for key battery properties, such as voltage, capacity, alkali diffusivity, and other electrochemically relevant computable quantities. We then present an overview of computational techniques for the study of rechargeable alkali-ion battery materials, followed by a critical review of the literature applying these techniques to yield crucial insights into battery operation and performance. Finally, we provide perspectives on outstanding challenges and opportunities in the theory and computation of rechargeable alkali-ion battery materials.

25 ENERGY STORAGE↗
DOE OSTI2023

Staggered scheme for the compressible fluctuating hydrodynamics of multispecies fluid mixtures

Here, we present a numerical formulation for the solution of nonisothermal, compressible Navier-Stokes equations with thermal fluctuations to describe mesoscale transport phenomena in multispecies fluid mixtures. The novelty of our numerical method is the use of staggered grid momenta along with a finite volume discretization of the thermodynamic variables to solve the resulting stochastic partial differential equations. The key advantages of the numerical scheme are that it significantly simplifies the discretization of diffusive and stochastic momentum fluxes into a more compact form, and it provides an unambiguous prescription of boundary conditions involving pressure. The staggered grid scheme more accurately reproduces the equilibrium static structure factor of hydrodynamic fluctuations in gas mixtures compared to a collocated scheme described previously by Balakrishnan et al. [Phys. Rev. E 89, 013017 (2014)1539-375510.1103/PhysRevE.89.013017]. The numerical method is tested for ideal noble gases mixtures under various nonequilibrium conditions, such as applied thermal and concentration gradients, to assess the role of cross-diffusion effects, such as Soret and Dufour, on the long-ranged correlations of hydrodynamic fluctuations, which are also more accurately reproduced compared to the collocated scheme. We numerically study giant nonequilibrium fluctuations driven by concentration gradients and fluctuation-driven Rayleigh-Taylor instability in gas mixtures. Wherever applicable, excellent agreement is observed with theory and measurements from the direct simulation Monte Carlo method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Application of Scalar Filtered Density Function to Turbulent Flows Under Supercritical Condition

Abstract The scalar filtered density function (FDF) methodology is extended and employed for large eddy simulation (LES) of turbulent flows under supercritical condition. To describe real fluid behavior, the extended methodology incorporates the generalized heat and mass diffusion models along with real fluid thermodynamic relations which are derived using the cubic Peng–Robinson equation of state. These models are implemented within the stochastic differential equations comprising the scalar FDF transport. Simulations are conducted of a temporally developing mixing layer under supercritical condition and the results are assessed by comparing with data generated by direct numerical simulation (DNS) of the same layer. The consistency of the proposed FDF methodology is assessed. The LES-FDF predictions are shown to agree favorably with the DNS data and exhibit several key features pertaining to supercritical turbulent flows.

Energy & Fuels↗