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,477 records · Page 82

Solar modulation of galactic cosmic rays

The current theory for the solar modulation of galactic cosmic rays is reviewed briefly. The four basic mechanisms which are important for modulation-diffusion, effects associated with the large-scale magnetic field, convection, and energy change - are discussed in simple terms. The equation which describes cosmic-ray behavior in the heliosphere is presented and some sample solutions are considered. It is pointed out that, although the basic mechanisms which describe the cosmic-ray behavior are understood, the actual cause of the variation in the cosmic-ray flux over the solar cycle, or over longer intervals, is beyond current comprehension.

Fisk, L. A.↗

The coupling between flame surface dynamics and species mass conservation in premixed turbulent combustion

Current flamelot models based on a description of the flame surface dynamics require the closure of two inter-related equations: a transport equation for the mean reaction progress variable, (tilde)c, and a transport equation for the flame surface density, Sigma. The coupling between these two equations is investigated using direct numerical simulations (DNS) with emphasis on the correlation between the turbulent fluxes of (tilde)c, bar(pu''c''), and Sigma, (u'')(sub S)Sigma. Two different DNS databases are used in the present work: a database developed at CTR by A. Trouve and a database developed by C. J. Rutland using a different code. Both databases correspond to statistically one-dimensional premixed flames in isotropic turbulent flow. The run parameters, however, are significantly different, and the two databases correspond to different combustion regimes. It is found that in all simulated flames, the correlation between bar(pu''c'') and (u'')(sub S)Sigma is always strong. The sign, however, of the turbulent flux of (tilde)c or Sigma with respect to the mean gradients, delta(tilde)c/delta(x) or delta(Sigma)/delta(x), is case-dependent. The CTR database is found to exhibit gradient turbulent transport of (tilde)c and Sigma, whereas the Rutland DNS features counter-gradient diffusion. The two databases are analyzed and compared using various tools (a local analysis of the flow field near the flame, a classical analysis of the conservation equation for (tilde)(u''c''), and a thin flame theoretical analysis). A mechanism is then proposed to explain the discrepancies between the two databases and a preliminary simple criterion is derived to predict the occurrence of gradient/counter-gradient turbulent diffusion.

Trouve, A.↗

A space-time discretization procedure for wave propagation problems

Higher order compact algorithms are developed for the numerical simulation of wave propagation by using the concept of a discrete dispersion relation. The dispersion relation is the imprint of any linear operator in space-time. The discrete dispersion relation is derived from the continuous dispersion relation by examining the process by which locally plane waves propagate through a chosen grid. The exponential structure of the discrete dispersion relation suggests an efficient splitting of convective and diffusive terms for dissipative waves. Fourth- and eighth-order convection schemes are examined that involve only three or five spatial grid points. These algorithms are subject to the same restrictions that govern the use of dispersion relations in the constructions of asymptotic expansions to nonlinear evolution equations. A new eighth-order scheme is developed that is exact for Courant numbers of 1, 2, 3, and 4. Examples are given of a pulse and step wave with a small amount of physical diffusion.

Davis, Sanford↗

On the modelling of non-reactive and reactive turbulent combustor flows

A study of non-reactive and reactive axisymmetric combustor flows with and without swirl is presented. Closure of the Reynolds equations is achieved by three models: kappa-epsilon, algebraic stress and Reynolds stress closure. Performance of two locally nonequilibrium and one equilibrium algebraic stress models is analyzed assuming four pressure strain models. A comparison is also made of the performance of a high and a low Reynolds number model for combustor flow calculations using Reynolds stress closures. Effects of diffusion and pressure-strain models on these closures are also investigated. Two models for the scalar transport are presented. One employs the second-moment closure which solves the transport equations for the scalar fluxes, while the other solves the algebraic equations for the scalar fluxes. In addition, two cases of non-premixed and one case of premixed combustion are considered. Fast- and finite-rate chemistry models are applied to non-premixed combustion. Both show promise for application in gas turbine combustors. However, finite rate chemistry models need to be examined to establish a suitable coupling of the heat release effects on turbulence field and rate constants.

Nikjooy, Mohammad↗

Charging of radioactive and environmental airborne particles

The charging of various airborne particles was investigated using single-particle levitation and charge-balance equations. Though radioactive decay and triboelectrification can induce charging, it is typically assumed that the aerosols in a radioactive plume will not carry significant charge at steady state since atmospheric particles can have their charge neutralized through the capture of adjacent counter-ions (i.e., diffusion charging). Here to assess this assumption, we directly measured the surface charge and charge density of various triboelectrically charged aerosols including radioactive uranium oxide (<1 μm), urban dust, Arizona desert dust, hydrophilic and hydrophobic silica nanoparticles, and graphene oxide powders using an electric field-assisted particle levitator in air. Of these particles, uranium oxide aerosols exhibited the highest surface charge density. Charge balance equations were employed to predict the average charge gained from radioactive decay as a function of time and to evaluate the effects of diffusion charging on triboelectrically charged radioactive and non-radioactive particles in the atmosphere. Simulation results show that particles, initially charged through triboelectrification, can be quickly discharged by diffusion charging in the absence of radioactive decay. Nevertheless, simulation results also indicate that particles can be strongly charged when they carry radionuclides. These experimental and simulation results suggest that radioactive decay can induce strong particle charging that may potentially affect atmospheric transport of airborne radionuclides.

54 ENVIRONMENTAL SCIENCES↗

Transcritical diffuse-interface hydrodynamics of propellants in high-pressure combustors of chemical propulsion systems

Rocket engines and high-power new generations of gas-turbine jet engines and diesel engines oftentimes involve the injection of one or more reactants at subcritical temperatures into combustor environments at high pressures, and more particularly at pressures higher than those corresponding to the critical points of the individual components of the mixture, which typically range from 13 to 50 bars for most propellants. This class of trajectories in the thermodynamic space has been traditionally referred to as transcritical. However, the fundamental understanding of fuel atomization, vaporization, mixing, and combustion processes at such high pressures remains elusive. In particular, whereas fuel sprays are relatively well characterized at normal pressures, analyses of dispersion of fuel in high-pressure combustors are hindered by the limited experimental diagnostics and theoretical formulations available. The description of the thermodynamics of hydrocarbon-fueled mixtures employed in chemical propulsion systems is complex and involves mixing-induced phenomena, including an elevation of the critical point whereby the coexistence region of the mixture extends up to pressures much larger than the critical pressures of the individual components. As a result, interfaces subject to surface-tension forces may persist in multicomponent systems despite the high pressures, and may give rise to unexpected spray-like atomization dynamics that are otherwise absent in monocomponent systems above their critical point. In this article, the current understanding of this phenomenon is reviewed within the context of propulsion systems fueled by heavy hydrocarbons. Emphasis is made on analytical descriptions at mesoscopic scales of interest for computational fluid dynamics. In particular, a set of modifications of the constitutive laws in the Navier–Stokes equations for multicomponent flows, supplemented with a high-pressure equation of state and appropriate redefinitions of the thermodynamic potentials, are introduced in this work based on an extended version of the diffuse-interface theory of van der Waals. The resulting formulation involves revisited forms of the stress tensor and transport fluxes of heat and species, and enables a description of the mesoscopic volumetric effects induced by transcritical interfaces consistently with the thermodynamic phase diagram of the mixture at high pressures. Additionally, applications of the theory are illustrated in canonical problems, including dodecane/nitrogen transcritical interfaces in non-isothermal systems. The results indicate that a transcritical interface is formed between the propellant streams that persists downstream of the injection orifice over distances of the same order as the characteristic thermal-entrance length of the fuel stream. The transcritical interface vanishes at an edge that gives rise to a fully supercritical mixing layer.

02 PETROLEUM↗

On the eigenvector bias of Fourier feature networks: From regression to solving multi-scale PDEs with physics-informed neural networks

Physics-informed neural networks (PINNs) are demonstrating remarkable promise in integrating physical models with gappy and noisy observational data, but they still struggle in cases where the target functions to be approximated exhibit high-frequency or multi-scale features. Here in this work we investigate this limitation through the lens of Neural Tangent Kernel (NTK) theory and elucidate how PINNs are biased towards learning functions along the dominant eigen-directions of their limiting NTK. Using this observation, we construct novel architectures that employ spatio-temporal and multi-scale random Fourier features, and justify how such coordinate embedding layers can lead to robust and accurate PINN models. Numerical examples are presented for several challenging cases where conventional PINN models fail, including wave propagation and reaction–diffusion dynamics, illustrating how the proposed methods can be used to effectively tackle both forward and inverse problems involving partial differential equations with multi-scale behavior.

42 ENGINEERING↗

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗

Lithospheric thinning associated with rifting in East Africa

The rate of ascent of the thermal perturbation of eruptive rock types in the Kenya Dome region can be calculated for depths of origin ranging from greater than 170 km to the lower crust (25-30 km) as combined with dated occurrences. An equation is derived and solved iteratively for ascent velocity at 10 km depth intervals from 150-30 km, assuming a thermal diffusivity of 1 sq mm/sec. For radii of 50 and 100 km and constant gradient increases within the range 55-60 C/km over the lithosphere-asthenosphere boundary (LAB) ascent depth range, the calculated rate of upwelling corresponds well with the results of a simple thermal model of lithospheric thinning. A 1 C/km change in gradient results approximately in a 0.55 Myr change in ascent time from 150 to 30 km, due to the increasing differential between the geotherm and the solidus with decreasing depths and the increasing rate of heat loss to the surface. It is suggested by the temporal-spatial-compositional relations of eruptives in the Kenya Dome region that there is a shallowing of magma source regions with time, implying encroachment of a thermal anomaly.

Wendlandt, R. F.↗

Experimental study of the thermal-acoustic efficiency in a long turbulent diffusion-flame burner

An acoustic source/propagation model is used to interpret measured noise spectra from a long turbulent burner. The acoustic model is based on the perturbation solution of the equations describing the unsteady one-dimensional flow of an inviscid ideal gas with a distributed heat source. The model assumes that the measured noise spectra are due uniquely to the unsteady component of combustion heat release. The model was applied to a long cylindrical hydrogen burner operating over a range of power levels between 4.5 kW and 22.3 kW. Acoustic impedances at the inlet to the burner and at the exit of the tube downstream of the burner were measured and are used as boundary conditions for the model. These measured impedances are also presented.

Mahan, J. R.↗

Application of program LAURA to three-dimensional AOTV flowfields

Program LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) is a robust, finite volume, single-level storage, implicit upwind differencing algorithm which has been documented and tested on several three-dimensional blunt-body flows. The algorithm can run at unlimited Courant numbers (relaxing the steady-state equations) but requires the inversion of only a 5 x 5 matrix per computational cell. An alternating directional sweep Gauss-Seidel substitution strategy is used to relax the governing equations. At present, the Euler and thin-layer Navier-Stokes equations using Sutherland's law for viscosity have been modeled for a perfect gas, equilibrium air, and nonequilibrium air chemistry neglecting diffusion. The equilibrium and nonequilibrium air chemistry options have been described in a companion paper. Good comparisons with experimental data and another calculation method for pressure distributions, aerodynamic coefficients, and heat-transfer distributions have been demonstrated for three-dimensional blunt-body flows.

Gnoffo, P. A.↗

Evaluation of transport conditions during physical vapor transport growth of opto-electronic crystals

Transport conditions were evaluated during the vapor phase growth of mercurous chloride crystals in a closed tube. Experimentally observed growth rates were much smaller than those calculated by the Hertz-Knudsen (H-K) equation. The Arrhenius behavior of growth rate with the temperature was used to derive the sticking coefficient. A one-dimensional diffusion model was used to calculate the total mass flux and was compared with the condensing flux. It was predicted that growth occurred in the convecto-diffusive range.

Singh, N. B.↗

Numerical solution of a three-dimensional cubic cavity flow by using the Boltzmann equation

A three-dimensional cubic cavity flow has been analyzed for diatomic gases by using the Boltzmann equation with the Bhatnagar-Gross-Krook (B-G-K) model. The method of discrete ordinate was applied, and the diffuse reflection boundary condition was assumed. The results, which show a consistent trend toward the Navier-Stokes solution as the Knudson number is reduced, give us confidence to apply the method to a three-dimensional geometry for practical predictions of rarefied-flow characteristics. The CPU time and the main memory required for a three-dimensional geometry using this method seem reasonable.

Hwang, Danny P.↗

The interaction of atomic oxygen with thin copper films

A source of thermal, ground-state atomic oxygen has been used to expose thin copper films at a flux of 1.4 x 10 exp 17 atoms/sq cm s for times up to 50 min for each of five temperatures between 140 and 200 C. Rutherford backscattering spectroscopy was used to characterize the oxide formed during exposure. The observations are consistent with the oxide phase Cu2O. The time dependence and the temperature dependence of the oxide layer thickness can be described using oxide film growth theory based on rate limitation by diffusion. Within the time and temperature ranges of this study, the growth of the oxide layers is well described by the equation L(T,t) = 3.6 x 10 to 8th exp(- 1.1/2k sub B T)t exp 1/2, where L,T, and t are measured in angstroms, degrees Kelvin, and minutes, respectively. The deduced activation energy is 1.10 +/- 0.15 eV, with the attendant oxidation rate being greater than that for the corresponding reaction in molecular oxygen.

Gibson, B. C.↗

The cosmic Doppler instability

The equations governing the behavior of perturbations of a mixture of nearly homogeneous and isotropic matter and radiation are derived, using a diffusion approximation where spatial perturbations in the radiation spectrum are allowed to vary with frequency. A simple model of line opacity leads to dispersion relations which display a new bulk instability. The model is used to derive an approximate dispersion relation for radiation interacting via resonance scattering opacity in atomic hydrogen at low density and low temperature. Possible applications to cosmology are briefly discussed.

Hogan, Craig J.↗

Development of an explicit multiblock/multigrid flow solver for viscous flows in complex geometries

A new computer program is being developed for doing accurate simulations of compressible viscous flows in complex geometries. The code employs the full compressible Navier-Stokes equations. The eddy viscosity model of Baldwin and Lomax is used to model the effects of turbulence on the flow. A cell centered finite volume discretization is used for all terms in the governing equations. The Advection Upwind Splitting Method (AUSM) is used to compute the inviscid fluxes, while central differencing is used for the diffusive fluxes. A four-stage Runge-Kutta time integration scheme is used to march solutions to steady state, while convergence is enhanced by a multigrid scheme, local time-stepping and implicit residual smoothing. To enable simulations of flows in complex geometries, the code uses composite structured grid systems where all grid lines are continuous at block boundaries (multiblock grids). Example results are shown a flow in a linear cascade, a flow around a circular pin extending between the main walls in a high aspect-ratio channel, and a flow of air in a radial turbine coolant passage.

Steinthorsson, E.↗

Development of an explicit multiblock/multigrid flow solver for viscous flows in complex geometries

A new computer program is being developed for doing accurate simulations of compressible viscous flows in complex geometries. The code employs the full compressible Navier-Stokes equations. The eddy viscosity model of Baldwin and Lomax is used to model the effects of turbulence on the flow. A cell centered finite volume discretization is used for all terms in the governing equations. The Advection Upwind Splitting Method (AUSM) is used to compute the inviscid fluxes, while central differencing is used for the diffusive fluxes. A four-stage Runge-Kutta time integration scheme is used to march solutions to steady state, while convergence is enhanced by a multigrid scheme, local time-stepping, and implicit residual smoothing. To enable simulations of flows in complex geometries, the code uses composite structured grid systems where all grid lines are continuous at block boundaries (multiblock grids). Example results shown are a flow in a linear cascade, a flow around a circular pin extending between the main walls in a high aspect-ratio channel, and a flow of air in a radial turbine coolant passage.

Steinthorsson, E.↗