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,801 records · Page 100
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively.

Allgood, Daniel C.↗
DOE OSTI2022

Finite orbit width effects in large aspect ratio stellarators

New orbit-averaged equations for low collisionality neoclassical fluxes in large aspect ratio stellarators with mirror ratios close to unity are derived. The equations retain finite orbit width effects by employing the second adiabatic invariant $J$ as a velocity-space coordinate and they have been implemented in the orbit-averaged neoclassical code KNOSOS (Velasco et al. , J. Comput. Phys. , vol. 418, 2020, 109512; Velasco et al. , Nucl. Fusion , vol. 61, 2021, 116013). The equations are used to study the $1/\nu$ regime and the lower collisionality regimes. For generic large aspect ratio stellarators with mirror ratios close to unity, as the collision frequency decreases, the $1/\nu$ regime transitions directly into the $\nu$ regime, without passing through a $\sqrt {\nu }$ regime. An explicit formula for the neoclassical fluxes in the $\nu$ regime is obtained. The formula includes the effect of particles that transition between different types of wells. While these transitions produce stochastic scattering independent of the value of the collision frequency in velocity space, the diffusion in real space remains proportional to the collision frequency. The $\sqrt {\nu }$ regime is only recovered in large aspect ratio stellarators close to omnigeneity: large aspect ratio stellarators with large mirror ratios and optimized large aspect ratio stellarators with mirror ratios close to unity. Neoclassical transport in large aspect ratio stellarators with large mirror ratios can be calculated with the orbit-averaged equations derived by Calvo et al. ( Plasma Phys. Control. Fusion , vol. 59, 2017, 055014). In these stellarators, the $\sqrt {\nu }$ regime exists in the collisionality interval $(a/R) \ln (R/a) \ll \nu _{ii} R a/\rho _i v_{ti} \ll R/a$ . In optimized large aspect ratio stellarators with mirror ratios close to unity, the $\sqrt {\nu }$ regime occurs in an interval of collisionality that depends on the deviation from omnigeneity $\delta$ : $\delta ^{2} |\ln \delta | \ll \nu _{ii} R a/\rho _i v_{ti} \ll 1$ . Here, $\nu _{ii}$ is the ion–ion collision frequency, $\rho _i$ and $v_{ti}$ are the ion gyroradius and thermal speed, and $a$ and $R$ are the minor and major radii.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Long-time integration of parametric evolution equations with physics-informed DeepONets

Ordinary and partial differential equations (ODEs/PDEs) play a paramount role in analyzing and simulating complex dynamic processes across all corners of science and engineering. In recent years machine learning tools are aspiring to introduce new effective ways of simulating such equations, however existing approaches are not able to reliably return stable and accurate predictions across long temporal horizons. We aim to address this challenge by introducing an effective framework for learning evolution operators that map random initial conditions to associated ODE/PDE solutions within a short time interval. Such operators can be parametrized by deep neural networks that are trained in an entirely self-supervised manner without requiring one to generate any paired input-output observations. Global long-time predictions across a range of initial conditions can be then obtained by iteratively evaluating the trained model using each prediction as the initial condition for the next evaluation step. Here, this introduces a new approach to temporal domain decomposition that is shown to be effective in performing accurate long-time simulations for a wide range of parametric ODE and PDE systems, from wave propagation, to reaction-diffusion dynamics and stiff chemical kinetics, introducing a new way of rapidly emulating non-equilibrium processes in science and engineering.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Multi-physics Peridynamic Modeling of Damage Processes in Protective Coatings

Protective coatings can revolutionize the aerospace industry by providing a way to build a fuel-efficient and low-emission gas turbine engine operating at elevated temperature. Currently, significant challenges persist in the multi-physics modeling of damage processes in protective coatings exposed to the hot corrosive environment of gas turbine engine. Examples of such detrimental processes include an interface damage between different material layers forming protective coating, a formation of the through-the-thickness vertical cracks, a diffusion of oxygen and moisture through those cracks, an oxidation, and a delamination. In this paper, a feasibility of a multi-physics peridynamic model to provide insight to these detrimental processes has been demonstrated. 3-Dimensional peridynamic equations which are valid everywhere, including dynamically evolving discontinuities, have been presented. In conclusion, validation and demonstration have been shown for conduction heat transfer, reactive oxidation, and delamination of protective coating with a composition borrowed from the open literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Electrochemical properties of poly(ethylene oxide) electrolytes above the entanglement threshold

The ion transport in electrolytes depends on three transport coefficients, conductivity (κ), salt diffusion coefficient (D), and the cation transference number with respect to the solvent velocity ($t_+^0$), and the thermodynamic factor ($T_f$). Current methods for determining these parameters involve four separate experiments, and the coupled nature of the equations used to determine them generally results in large experimental uncertainty. We present data obtained from 64 independent polymer electrolytes comprising poly(ethylene oxide) (PEO) and lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) salt. The molecular weights of PEO ranged from 5 to 275 kg mol -1 ; these samples are all above the entanglement threshold. We minimize the experimental uncertainty in transport and thermodynamic measurements by exploiting the fact that ion transport in entangled polymer electrolytes should be independent of molecular weight. The dependence of $κ, D, t_+^0$, and $T_f$ as a function of salt concentration in the range 0.035 ≤ r ≤ 0.30 are presented with a 95% confidence interval, where r is the molar ratio of lithium ions to ethylene oxide monomer units. While κ, D, and $T_f$ are all positive as required by thermodynamic constraints, there is no constraint on the sign of $t_+^0$. We find that $t_+^0$ is negative in the salt concentration range of 0.093 ≤ r ≤ 0.189.

25 ENERGY STORAGE↗
DOE OSTI2021

Nuclear energy density functionals grounded in ab initio calculations

Here, we discuss the construction of a nuclear energy density functional (EDF) from ab initio computations and advocate the need for a methodical approach that is free from ad hoc assumptions. The equations of state (EoSs) of symmetric nuclear and pure neutron matter are computed using the chiral NNLO sat and the phenomenological AV4' + UIX c Hamiltonians as inputs to self-consistent Green's function (SCGF) and auxiliary field diffusion Monte Carlo (AFDMC) methods. We propose a convenient parametrization of the EoS as a function of the Fermi momentum and fit it on the SCGF and AFDMC calculations. We apply the ab initio based EDF to carry out an analysis of the binding energies and charge radii of different nuclei in the local density approximation. The NNLO sat -based EDF produces encouraging results, whereas the AV4' + UIX c -based one is farther from experiment. Possible explanations of these different behaviors are suggested, and the importance of gradient and spin-orbit terms is analyzed. Our paper paves the way for a practical and systematic way to merge ab initio nuclear theory and density functional theory, while shedding light on some critical aspects of this procedure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Assimilating partial observation to enhance feedback control of stochastic dynamical systems

Here, in this paper, we present a novel methodology to tackle feedback optimal control problems in scenarios where the exact state of the controlled process is unknown. It integrates data assimilation techniques and optimal control solvers to manage partial observation of the state process, a common occurrence in practical scenarios. Traditional stochastic optimal control methods assume full state observation, which is often not feasible in real-world fluid dynamics control problems. Our approach underscores the significance of utilizing observational data to inform control policy design. Specifically, we introduce a kernel learning backward stochastic differential equation (SDE) filter to enhance data assimilation efficiency and propose a sample-wise stochastic optimization method within the stochastic maximum principle framework. We demonstrate the efficacy and accuracy of our method in the control of advection-diffusion-reaction flow problem and the Dubins airplane maneuvering problem with model uncertainty.

data driven↗
NASA NTRS1973

A study of the force-field equation for the propagation of galactic cosmic rays

A new development is given of the solution of the equation of the force-field approximation for the propagation of galactic cosmic rays in the interplanetary region. It leads to simpler methods for determining the force-field parameters. A method is given for determining the separable diffusion coefficient from observations of galactic electron spectrum and near-earth electron spectra; it is shown that this diffusion coefficient is not unique but may have a periodic-like dependence upon rigidity; and the method is used to obtain diffusion coefficients for 1965 and 1968. Approximate formulae relating small changes in intensity and diffusion coefficient are developed and some applications of these noted; in one it is shown that the form of, and changes in, diffusion coefficient deduced previously for a neutron monitor event during June-September 1969 are unnecessarily constrained and therefore probably not correct.

Gleeson, L. J.↗
NASA NTRS1974

Analytic approximations in the study of the solar modulation of electrons

Numerical solutions to the transport equation of galactic cosmic rays in the interplanetary medium have been used to investigate the applicability of commonly used approximate analytic solutions for electrons (10 MeV to 10 GeV). We find that for a given cosmic-ray diffusion coefficient the force-field approximation is in reasonable agreement with the numerical solution at energies above 200 MeV, but deviates significantly at lower energies, depending on the shape of the interstellar electron spectrum. The diffusion-convection approximation agrees generally with the numerical solution within a factor of 2 over the entire energy range.

Cummings, A. C.↗
NASA NTRS1976

Theoretical investigations of plasma processes

System analyses are presented for electrically sustained, collision dominated plasma centrifuges, in which the plasma rotates under the influence of the Lorentz forces resulting from the interaction of the current density fields with an external magnetic field. It is shown that gas discharge centrifuges are technically feasible in which the plasma rotates at speeds up to 1 million cm/sec. The associated centrifugal forces produce a significant spatial isotope separation, which is somewhat perturbed in the viscous boundary layers at the centrifuge walls. The isotope separation effect is the more pronounced. The induced magnetic fields have negligible influence on the plasma rotation if the Hall coefficient is small. In the technical realization of collision dominated plasma centrifuges, a trade-off has to be made between power density and speeds of rotation. The diffusion of sputtered atoms to system surfaces of ion propulsion systems and the deposition of the atoms are treated theoretically by means of a simple model which permits an analytical solution. The problem leads to an inhomogeneous integral equation.

Wilhelm, H. E.↗
NASA NTRS1982

Convective thinning of the lithosphere - A mechanism for the initiation of continental rifting

A model of lithospheric thinning, in which heat is convected to the base and conducted within the lithosphere, is presented. An analytical equation for determinining the amount of thinning attainable on increasing the heat flux from the asthenosphere is derived, and a formula for lithosphere thickness approximations as a function of time is given. Initial and final equilibrium thicknesses, thermal diffusivity, transition temperature profile, and plume temperature profile are all factors considered for performing rate of thinning determinations. In addition, between initial and final equilibrium states, lithospheric thinning occurs at a rate which is inversely proportional to the square root of the time. Finally, uplift resulting from thermal expansion upon lithospheric thinning is on the order of 10 to the 2nd to 10 to the 3rd m.

Spohn, T.↗
NASA NTRS1985

Investigation of parabolic computational techniques for internal high-speed viscous flows

A feasibility study was conducted to assess the applicability of an existing parabolic analysis (ADD-Axisymmetric Diffuser Duct), developed previously for subsonic viscous internal flows, to mixed supersonic/subsonic flows with heat addition simulating a SCRAMJET combustor. A study was conducted with the ADD code modified to include additional convection effects in the normal momentum equation when supersonic expansion and compression waves were present. It is concluded from the present study that for the class of problems where strong viscous/inviscid interactions are present a global iteration procedure is required.

Anderson, O. L.↗
NASA NTRS1985

Numerical investigation of internal high-speed viscous flows using a parabolic technique

A feasibility study has been conducted to assess the applicability of an existing parabolic analysis (ADD-Axisymmetric Diffuser Duct), developed previously for subsonic viscous internal flows, to mixed supersonic/subsonic flows with heat addition simulating a SCRAMJET combustor. A study was conducted with the ADD code modified to include additional convection effects in the normal momentum equation when supersonic expansion and compression waves are present. A set of test problems with weak shock and expansion waves have been analyzed with this modified ADD method and stable and accurate solutions were demonstrated provided the streamwise step size was maintained at levels larger than the boundary layer displacement thickness. Calculations made with further reductions in step size encountered departure solutions consistent with strong interaction theory. Calculations were also performed for a flow field with a flame front in which a specific heat release was imposed to simulate a SCRAMJET combustor. In this case the flame front generated relatively thick shear layers which aggravated the departure solution problem. Qualitatively correct results were obtained for these cases using a marching technique with the convective terms in the normal momentum equation suppressed. It is concluded from the present study that for the class of problems where strong viscous/inviscid interactions are present a global iteration procedure is required.

Anderson, O. L.↗
NASA NTRS1988

Evaluation of surface energy and radiation balance systems for FIFE

The energy balance and radiation balance components were determined at six sites during the First International Satellite Land Surface Climatology Project Field Experiment (FIFE) conducted south of Manhattan, Kansas during the summer of 1987. The objectives were: to determine the effect of slope and aspect, throughout a growing season, on the magnitude of the surface energy balance fluxes as determined by the Energy Balance Method (EBM); to investigate the calculation of the soil heat flux density at the surface as calculated from the heat capacity and the thermal conductivity equations; and to evaluate the performance of the Surface Energy and Radiation Balance System (SERBS). A total of 17 variables were monitored at each site. They included net, solar (up and down), total hemispherical (up and down), and diffuse radiation, soil temperature and heat flux density, air and wet bulb temperature gradients, wind speed and direction, and precipitation. A preliminary analysis of the data, for the season, indicate that variables including net radiation, air temperature, vapor pressure, and wind speed were quite similar at the sites even though the sites were as much as 16 km apart and represented four cardinal slopes and the top of a ridge.

Fritschen, Leo J.↗
NASA NTRS1987

Fermi acceleration at shocks with arbitrary velocity profiles

Diffusive particle acceleration in modified shock fronts is considered. A method by means of which to obtain the resulting particle spectrum for an arbitrary flow profile is outlined; this method is based on a reformulation of the transport equation as a nonlinear first-order equation. Besides the shape of the flow profile, the spectral index of particles depends on the Peclet number of the shock. Results are given for some sample profiles.

Schneider, Peter↗
NASA NTRS1994

Probability density function approach for compressible turbulent reacting flows

The objective of the present work is to extend the probability density function (PDF) tubulence model to compressible reacting flows. The proability density function of the species mass fractions and enthalpy are obtained by solving a PDF evolution equation using a Monte Carlo scheme. The PDF solution procedure is coupled with a compression finite-volume flow solver which provides the velocity and pressure fields. A modeled PDF equation for compressible flows, capable of treating flows with shock waves and suitable to the present coupling scheme, is proposed and tested. Convergence of the combined finite-volume Monte Carlo solution procedure is discussed. Two super sonic diffusion flames are studied using the proposed PDF model and the results are compared with experimental data; marked improvements over solutions without PDF are observed.

Hsu, A. T.↗
NASA NTRS1993

Particle-Gas Dynamics in the Midplane of a Protoplanetary Nebula

In this paper, we study the stage of planetary formation during which the particulate material is still dispersed as centimeter-to-meter sized primordial aggregates. During this stage, particles are able to settle toward the midplane into a layer of mass density comparable to or much greater than that of the gas. Nonlinear, coupled interactions between the particles and the nebula gas become significant and ultimately determine the vertical profiles of the particle density and the mean velocities of the particles and the nearby gas. This is the environment in which the earliest planetesimals probably form. Our numerical models rely on the Reynolds averaged Navier-Stokes equations for the gas and particles and are fully viscous and turbulent (and, for the particles. compressible). Our turbulence modeling uses a Prandtl local shear parametrization, validated by laboratory experiments. We have developed a new model for particle diffusion involving the profile of the nebula gas turbulence and the particle Schmidt number, which is a function of particle size and density. We model a cool, quiescent nebula at 1 AU (280 K) and 10 AU (90 K) and a higher temperature stage (1000 K) at 1 AU. Our main results include: (a) rapid accretion of planetesimals by gravitationally unstable fragmentation on an orbital time scale (the 'Goldreich-Ward instability') is unlikely to occur until objects have already accreted by some other process to the mass of the largest known meteorite samples, if at all; (b) from 'seeds' as small as 10 m, growth of 10- to 100-km planetesimals can proceed rapidly by drift-augmented accretion in the particle-rich midplane with orbital decay of about 1% for the growing planetesimals; (c) outward transport of vapor and small entrained chips can account for significant radial compositional and mineralogical mixing in primitive meteorite parent bodies. Other implications may be drawn concerning particle random collisional velocities and the time scales of midplane turbulence.

Cuzzi, Jeffrey N.↗
NASA NTRS2011

3D Modeling of Ultrasonic Wave Interaction with Disbonds and Weak Bonds

Ultrasonic techniques, such as the use of guided waves, can be ideal for finding damage in the plate and pipe-like structures used in aerospace applications. However, the interaction of waves with real flaw types and geometries can lead to experimental signals that are difficult to interpret. 3-dimensional (3D) elastic wave simulations can be a powerful tool in understanding the complicated wave scattering involved in flaw detection and for optimizing experimental techniques. We have developed and implemented parallel 3D elastodynamic finite integration technique (3D EFIT) code to investigate Lamb wave scattering from realistic flaws. This paper discusses simulation results for an aluminum-aluminum diffusion disbond and an aluminum-epoxy disbond and compares results from the disbond case to the common artificial flaw type of a flat-bottom hole. The paper also discusses the potential for extending the 3D EFIT equations to incorporate physics-based weak bond models for simulating wave scattering from weak adhesive bonds.

Leckey, C.↗

Refine your search to explore more results.