Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

FDNS CFD Code Benchmark for RBCC Ejector Mode Operation: Continuing Toward Dual Rocket Effects

Computational Fluid Dynamics (CFD) analysis results are compared with benchmark quality test data from the Propulsion Engineering Research Center's (PERC) Rocket Based Combined Cycle (RBCC) experiments to verify fluid dynamic code and application procedures. RBCC engine flowpath development will rely on CFD applications to capture the multi -dimensional fluid dynamic interactions and to quantify their effect on the RBCC system performance. Therefore, the accuracy of these CFD codes must be determined through detailed comparisons with test data. The PERC experiments build upon the well-known 1968 rocket-ejector experiments of Odegaard and Stroup by employing advanced optical and laser based diagnostics to evaluate mixing and secondary combustion. The Finite Difference Navier Stokes (FDNS) code [2] was used to model the fluid dynamics of the PERC RBCC ejector mode configuration. Analyses were performed for the Diffusion and Afterburning (DAB) test conditions at the 200-psia thruster operation point, Results with and without downstream fuel injection are presented.

West, Jeff↗
NASA NTRS2014

Study of Aerothermodynamic Modeling Issues Relevant to High-Speed Sample Return Vehicles

This paper examines the application of state-of-the-art coupled ablation and radiation simulations to highspeed sample return vehicles, such as those returning from Mars or an asteroid. A defining characteristic of these entries is that the surface recession rates and temperatures are driven by nonequilibrium convective and radiative heating through a boundary layer with significant surface blowing and ablation products. Measurements relevant to validating the simulation of these phenomena are reviewed and the Stardust entry is identified as providing the best relevant measurements. A coupled ablation and radiation flowfield analysis is presented that implements a finite-rate surface chemistry model. Comparisons between this finite-rate model and a equilibrium ablation model show that, while good agreement is seen for diffusion-limited oxidation cases, the finite-rate model predicts up to 50% lower char rates than the equilibrium model at sublimation conditions. Both the equilibrium and finite rate models predict significant negative mass flux at the surface due to sublimation of atomic carbon. A sensitivity analysis to flowfield and surface chemistry rates show that, for a sample return capsule at 10, 12, and 14 km/s, the sublimation rates for C and C3 provide the largest changes to the convective flux, radiative flux, and char rate. A parametric uncertainty analysis of the radiative heating due to radiation modeling parameters indicates uncertainties ranging from 27% at 10 km/s to 36% at 14 km/s. Applying the developed coupled analysis to the Stardust entry results in temperatures within 10% of those inferred from observations, and final recession values within 20% of measurements, which improves upon the 60% over-prediction at the stagnation point obtained through an uncoupled analysis. Emission from CN Violet is shown to be over-predicted by nearly and order-of-magnitude, which is consistent with the results of previous independent analyses. Finally, the coupled analysis is applied to a 14 km/s Earth entry representative of a Mars sample return. Although the radiative heating provides a larger fraction of the total heating, the influence of ablation and radiation on the flowfield are shown to be similar to Stardust.

Johnston, Christopher O.↗
DOE OSTI2013

Investigation of the possibility to use a fine-mesh solver for resolving coupled neutronics and thermal-hydraulics

The development of a fine-mesh coupled neutronic/thermal-hydraulic solver is touched upon in this paper. The reported work investigates the feasibility of using finite volume techniques to discretize a set of conservation equations modeling neutron transport, fluid dynamics, and heat transfer within a single numerical tool. With the long-term objective of developing fine-mesh computing capabilities for a few selected fuel assemblies in a nuclear core, this preliminary study considers an infinite array of a single fuel assembly having a finite height. Thermal-hydraulic conditions close to the ones existing in PWRs are taken as a first test case. The neutronic modeling relies on the diffusion approximation in a multi-energy group formalism, with cross-sections pre-calculated and tabulated at the sub-pin level using a Monte Carlo technique. The thermal-hydraulics is based on the Navier-Stokes equations, complemented by an energy conservation equation. The non-linear coupling terms between the different conservation equations are fully resolved using classical iteration techniques. Early tests demonstrate that the numerical tool provides an unprecedented level of details of the coupled solution estimated within the same numerical tool and thus avoiding any external data transfer, using fully consistent models between the neutronics and the thermal-hydraulics. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

Towards High Performance Oxide-Dispersion-Strengthened Alloys

The current first-generation candidate swelling-resistant alloys envisioned to replace austenitics are ferritic and ferritic-martensitic steels, with second-generation alloys being yttria dispersion-hardened variants of these steels. In ferritic-martensitic alloys, ferrite grains always swell earlier than martensite. Currently, most oxide dispersion-strengthened (ODS) alloys have only a ferrite matrix, thus giving up the first line of defense against swelling provided by tempered martensite. Towards development of better swelling-resistant ODS alloys, dual phase (ferrite + martensite) ODS alloys with oxide particles uniformly distributed in both ferrite and martensite phases is needed. Furthermore, fundamental understanding on mechanisms governing dispersoid stability and its correlation with void swelling is greatly needed. We reported here further alloy improvement by exploring new 9Cr ODS alloy with majority of oxide dispersoids uniformly distributed in tempered martensite grains. The alloy was irradiated by 3.5 MeV Fe self ions up to 600 dpa (displacements per atom) and at temperatures ranging from 325 to 625 deg. C. According to SRIM code, the dpa profile peaks at about 1000 nm. First of all, ion irradiation reduces the dispersoid sizes. The dispersoid sizes begin to increase beyond ∼1000 nm. Second, higher irradiation temperatures result in larger dispersoids on average. Third, there is no significant dispersoid size difference between 100 and 200 dpa irradiation. The equilibrium sizes are quickly saturated. Size and coherency of a dispersoid is closely related. In martensite grains, fine dispersoids tend to be coherent while coarse dispersoids tend to incoherent. We found that at a low irradiation temperature, coherent small dispersoids become smaller and their interface coherency remains. In a comparison, incoherent large dispersoids disappear after irradiation. At an intermediate irradiation temperature, the final equilibrium dispersoid sizes increase. At a high irradiation temperature, the final dispersoid sizes further increase and some large semi-coherent dispersoids survive. In ferrite grains, as-received dispersoids are relatively smaller in comparison with martensite grains. Their equilibrium size changes after irradiation follow a similar trend but large dispersoids completely disappear even under high temperature ion irradiation. We propose a model to explain the above dispersoid size changes under various irradiation conditions. The model considers the competing effects between ballistic-collision- assisted dispersoid dissolution into matrix and thermally driven back-diffusion for recovery. When the competing reaches a balance, an equilibrium size distribution curve is determined. For dispersoids smaller than the equilibrium size, they grow upon ion irradiation. Large dispersoids may shrink and reach the equilibrium size but ultra-large ones will completely disappear if they shrinkage paths are far away from the curve. The model further considers the complexity from the dpa rates, and damage cascade morphology.

36 MATERIALS SCIENCE↗
NASA NTRS2002

Simultaneous Raman-Rayleigh-LIF Measurements and Numerical Modeling Results of a Lifted H2/N2 Turbulent Jet Flame in a Vitiated Coflow

An experimental and numerical investigation is presented of a H2/N2 turbulent jet flame burner that has a novel vitiated coflow. The vitiated coflow emulates the recirculation region of most combustors, such as gas turbines or furnaces. Additionally, since the vitiated gases are coflowing, the burner allows for exploration of recirculation chemistry without the corresponding fluid mechanics of recirculation. Thus the vitiated coflow burner design facilitates the development of chemical kinetic combustion models without the added complexity of recirculation fluid mechanics. Scalar measurements are reported for a turbulent jet flame of H2/N2 in a coflow of combustion products from a lean ((empty set) = 0.25) H2/Air flame. The combination of laser-induced fluorescence, Rayleigh scattering, and Raman scattering is used to obtain simultaneous measurements of the temperature, major species, as well as OH and NO. Laminar flame calculation with equal diffusivity do agree when the premixing and preheating that occurs prior to flame stabilization is accounted for in the boundary conditions. Also presented is an exploratory pdf model that predicts the flame's axial profiles fairly well, but does not accurately predict the lift-off height.

Cabra, R.↗
DOE OSTI2020

The Numerical Computation of Clear-Sky Transmittance and Reflectance of Solar Radiation for the Modeling of PV System Performance

Modeling of photovoltaics (PV) system performance requires numerical computation of plane-of-array (POA) irradiances that are influenced by atmospheric absorption and scattering. Clear-sky transmittance and reflectance of diffuse irradiance are essential for the understanding of solar radiation related to multiple reflection between land surface and clouds. This study analytically solves the radiative transfer equation for clear-sky conditions and numerically computes the transmittance and reflectance of diffuse irradiance as functions of the optical thickness of the atmosphere. The computation will be implemented in all-sky radiative transfer models to improve their computational efficiency.

14 SOLAR ENERGY↗
DOE OSTI2023

Investigation of desorption of hydrogen gas from polymer matrix using thermal desorption analysis and finite element modeling

Polymer materials are used widely as sealing materials and liners in hydrogen storage and transportation infrastructure. During operating conditions, polymers are exposed to high pressure hydrogen gas periodically. When in contact, the hydrogen gas diffuses through the polymer material with ease due to its smaller sized molecules compared to the polymer molecules. In the event of rapid decompression of outside hydrogen gas, the hydrogen molecules within the polymer diffuse out. This process is known as the desorption of polymer and can cause blisters or cracking causing permanent damage to the polymer. Here, we used a combined experimental and modeling approach to investigate the desorption of the hydrogen gas in the polymer. Thermal desorption analysis (TDA) experiments were performed on the polymer sample exposed to the high-pressure hydrogen gas to measure the diffusion coefficient and equilibrium hydrogen concentration. Experiments also provided the quantity of hydrogen within the polymer during the desorption process which followed Fick’s law. However, for higher pressurization values, some anomaly was observed in the quantity of hydrogen remaining in the polymer sample. A continuum mechanics based coupled diffusion–deformation-damage model was deployed to simulate the gas desorption process in the material using the finite element method to understand its effect on the material. It was found that the anomaly in TDA results were due to the blister formation on the polymer surface which resulted in longer hydrogen retention. Simulations also showed that blister formation can be affected by multiple parameters such as the addition of filler particles, and pre-existing cavity size.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermal power plant upgrade via a rotating detonation combustor and retrofitted turbine with optimized endwalls

Over the past decade, pressure gain combustion research has promised over 10 percentage-points of increase in power plant thermal efficiency. Alas, to realize such potential gain, one must effectively couple the turbine with the detonation combustor, whose exhaust conditions differ substantially from current state of the art gas turbines. This paper presents a modeling approach that enables a superior thermodynamic cycle with a rotating detonation combustor and a retrofitted gas turbine by including a diffuser downstream of the combustor and by contouring the turbine endwall, while preserving the airfoil geometry. We propose a multi-step optimization strategy, parametrizing the endwall geometry with a few control points, without altering the airfoil geometry, and with the stage turbine efficiency as objective function. In a first step the turbine performance is assessed with steady inlet conditions by solving the steady three-dimensional Reynolds-Averaged Navier-Stokes equations. In a second step, the inlet conditions are unsteady, as predicted from a detonation combustor and a diffuser, and three-dimensional full unsteady simulations are performed with an unsteady Reynolds-Averaged Navier-Stokes solver. By altering the vane endwall, the steady optimization yielded an efficiency increase of 12% relative to the baseline, while the unsteady optimization resulted in 21% increase compared to the datum turbine. Finally, a full engine analysis demonstrated the superiority of pressure gain combustion which included a realistic thermodynamic cycle of the combustor, the diffuser and the optimized turbine components.

42 ENGINEERING↗
DOE OSTI2024

Subsurface Redox Interactions Regulate Ebullitive Methane Flux in Heterogeneous Mississippi River Deltaic Wetland

Abstract As interfaces connecting terrestrial and ocean ecosystems, coastal wetlands develop temporally and spatially complex redox conditions, which drive uncertainties in greenhouse gas emission as well as the total carbon budget of the coastal ecosystem. To evaluate the role of complex redox reactions in methane emission from coastal wetlands, a coupled reactive‐transport model was configured to represent subsurface biogeochemical cycles of carbon, nitrogen, and sulfur, along with production and transport of multiple gas species through diffusion and ebullition. This model study was conducted at multiple sites along a salinity gradient in the Barataria Basin at the Mississippi River Deltaic Plain. Over a freshwater to saline gradient, simulated total flux of methane was primarily controlled by its subsurface production and consumption, which were determined by redox reactions directly (e.g., methanogenesis, methanotrophy) and indirectly (e.g., competition with sulfate reduction) under aerobic and/or anaerobic conditions. At fine spatiotemporal scales, surface methane fluxes were also strongly dependent on transport processes, with episodic ebullitive fluxes leading to higher spatial and temporal variability compared to the gradient‐driven diffusion flux. Ebullitive methane fluxes were determined by methane fraction in total ebullitive gas and the frequency of ebullitive events, both of which varied with subsurface methane concentrations and other gas species. Although ebullition thresholds are constrained by local physical factors, this study indicates that redox interactions not only determine gas composition in ebullitive fluxes but can also regulate ebullition frequency through gas production.

03 NATURAL GAS↗
DOE OSTI2022

Modeling the Dynamics of Energetic Protons in Earth's Inner Magnetosphere

Recently, Van Allen Probes data have revealed the interesting energy-dependent dynamics of energetic protons in Earth's inner magnetosphere. To quantify the role of radial diffusion and charge exchange to the observed energy-dependent dynamics of ring current protons, a 1D radial diffusion model with charge exchange loss is implemented. The observed proton flux over the long-term period of November 2012 to September 2013 is first converted to phase space densities, which are then simulated by our model driven by an outer boundary condition at L* = 5.5 derived from the data. The simulation results show that our model generally captures the transport and acceleration of energetic protons at μ = 30, 50, and 80 MeV/G and K = 0.11 G 1/2 R E , suggesting that radial diffusion is the dominant source mechanism for >75 keV protons at L* = 3.5–5.5. Additionally, the observed fast decay of protons at lower μ and slow decay at higher μ are well captured by the model, demonstrating the dominant role of charge exchange in explaining the observed energy-dependent proton decay. For higher μ protons, prompt losses of protons on the time scale of hours are observed over a wide range of L*, which is too fast to be explained by charge exchange. Some of the prompt losses at high L* regions are reproduced by the model with outward radial diffusion to the outer boundary. However, many prompt losses observed at lower L* regions are not captured by the model, which could be due to other loss mechanisms including electromagnetic ion cyclotron wave scattering and field line curvature scattering.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2000

Measurements and Modeling of Nitric Oxide Formation in Counterflow, Premixed CH4/O2/N2 Flames

Laser-induced fluorescence (LIF) measurements of NO concentration in a variety of CH4/O2/N2 flames are used to evaluate the chemical kinetics of NO formation. The analysis begins with previous measurements in flat, laminar, premixed CH4/O2/N2 flames stabilized on a water-cooled McKenna burner at pressures ranging from 1 to 14.6 atm, equivalence ratios from 0.5 to 1.6, and volumetric nitrogen/oxygen dilution ratios of 2.2, 3.1 and 3.76. These measured results are compared to predictions to determine the capabilities and limitations of the comprehensive kinetic mechanism developed by the Gas Research Institute (GRI), version 2.11. The model is shown to predict well the qualitative trends of NO formation in lean-premixed flames, while quantitatively underpredicting NO concentration by 30-50%. For rich flames, the model is unable to even qualitatively match the experimental results. These flames were found to be limited by low temperatures and an inability to separate the flame from the burner surface. In response to these limitations, a counterflow burner was designed for use in opposed premixed flame studies. A new LIF calibration technique was developed and applied to obtain quantitative measurements of NO concentration in laminar, counterflow premixed, CH4/O2/N2 flames at pressures ranging from 1 to 5.1 atm, equivalence ratios of 0.6 to 1.5, and an N2/O2 dilution ratio of 3.76. The counterflow premixed flame measurements are combined with measurements in burner-stabilized premixed flames and counterflow diffusion flames to build a comprehensive database for analysis of the GRI kinetic mechanism. Pathways, quantitative reaction path and sensitivity analyses are applied to the GRI mechanism for these flame conditions. The prompt NO mechanism is found to severely underpredict the amount of NO formed in rich premixed and nitrogen-diluted diffusion flames. This underprediction is traced to uncertainties in the CH kinetics as well as in the nitrogen oxidation chemistry. Suggestions are made which significantly improve the predictive capability of the GRI mechanism in near-stoichiometric, rich, premixed flames and in atmospheric-pressure, diffusion flames. However, the modified reaction mechanism is unable to model the formation of NO in ultra-rich, premixed or in high-pressure, nonpremixed flames, thus indicating the need for additional study under these conditions.

Thomsen, D. Douglas↗
NASA NTRS2001

Numerical Study of Buoyancy and Different Diffusion Effects on the Structure and Dynamics of Triple Flames

Numerical simulations of 2-D triple flames under gravity force have been implemented to identify the effects of gravity on triple flame structure and propagation properties and to understand the mechanisms of instabilities resulting from both heat release and buoyancy effects. A wide range of gravity conditions, heat release, and mixing widths for a scalar mixing layer are computed for downward-propagating (in the same direction with the gravity vector) and upward-propagating (in the opposite direction of the gravity vector) triple flames. Results of numerical simulations show that gravity strongly affects the triple flame speed through its contribution to the overall flow field. A simple analytical model for the triple flame speed, which accounts for both buoyancy and heat release, is developed. Comparisons of the proposed model with the numerical results for a wide range of gravity, heat release and mixing width conditions, yield very good agreement. The analysis shows that under neutral diffusion, downward propagation reduces the triple flame speed, while upward propagation enhances it. For the former condition, a critical Froude number may be evaluated, which corresponds to a vanishing triple flame speed. Downward-propagating triple flames at relatively strong gravity effects have exhibited instabilities. These instabilities are generated without any artificial forcing of the flow. Instead disturbances are initiated by minute round-off errors in the numerical simulations, and subsequently amplified by instabilities. A linear stability analysis on mean profiles of stable triple flame configurations have been performed to identify the most amplified frequency in spatially developed flows. The eigenfunction equations obtained from the linearized disturbance equations are solved using the shooting method. The linear stability analysis yields reasonably good agreements with the observed frequencies of the unstable triple flames. The frequencies and amplitudes of disturbances increase with the magnitude of the gravity vector. Moreover, disturbances appear to be most amplified just downstream of the premixed branches. The effects of mixing width and differential diffusion are investigated and their roles on the flame stability are studied.

Chen, Jyh-Yuan↗
NASA NTRS1997

Temperature and Radiative Heat Flux Measurements in Microgravity Jet Diffusion Flames

The objective of this project is to provide detailed measurements and modeling analyses of local soot concentration, temperature and radiation heat flux distributions in laminar and turbulent jet diffusion flames under normal (1-g) and reduced gravity (0-g) conditions. Results published to date by these co-PI's and their co-workers include: 1. thermophoretic sampling and size and morphological analyses of soot aggregates in laminar flames under normal and reduced gravity conditions; 2. full-field absorption imaging for soot volume fraction maps in laminar and turbulent flames under normal and reduced gravity conditions; 3. an accurate solver module for detailed radiation heat transfer in nongray nonhomogeneous media; 4. a complete model to include flame structure, soot formation and an energy equation to couple with radiation solver.

Ku, Jerry C.↗
DOE OSTI2020

Graphite Lithiation under Fast Charging Conditions: Atomistic Modeling Insights

The charging of lithium ion batteries in a fast and safe manner is critical for promoting the mass adoption of electric vehicles. Li intercalation in graphite electrodes is known to be one of the bottlenecks during the fast charging process. The mechanism of Li diffusion in highly polarized graphite anode at high current rates remains, however, not well understood. In this work, Density Functional Theory (DFT) calculations are used to gain insights into the Li diffusion process in graphite when it is far from equilibrium under fast charging conditions. The effect of uncompensated charges on Li mobility is determined in the highly polarized regions of the anode close to the interfaces. The extra charge was found to increase the interlayer spacing in the diffusion layer and adjacent channels, increasing the diffusivity and promoting the formation of Li clusters. A concerted diffusion mechanism at the edge of high-concentration Li domains is proposed to enhance the diffusion of Li.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

Meso-Beta scale numerical simulation studies of terrain-induced jet streak mass/momentum perturbations

The Control-B simulation experiment was designed to improve on the Control-A experiment performed with the GMASS model. This experiment addressed several inadequacies with the first smooth terrain numerical simulation by including: (1) increased nested-grid resolution to better define the simulated gravity waves, (2) increased horizontal diffusion to remove outflow boundary condition noise, and (3) the use of reanalyzed rawinsonde data and surface observations in the initial state to increase the definition of the observed jet streak as well as other low-level features. A smoothed-terrain dry simulation with the nested-grid GMASS model has revealed many important aspects of the processes which resulted in the generation of gravity waves in the region and time when and where they were observed. However, the vertical structure, number, and characteristics of the waves are still quite different from observed waves as diagnosed thus necessitating future improved simulations. However, this control simulation has produced substantial insight into processes which occur on many spatial scales over a 30 hour time period thus allowing one to draw promising inferences as to the mechanisms for the complex process which occurred in nature during the CCOPE case study. The theoretical aspects of the project have focus on understanding the nature of the ageostrophic circulations which are produced in idealized models of the atmosphere in which the troposphere is modeled in one of two ways. The first model assumes that the lower atmosphere can be represented as a single layer of homogeneous fluid whose upper surface is free to exhibit vertical displacement. Two-dimensional internal convergence (divergence) occurring during the adjustment to an asymptotic equilibrium state from an ageostrophic initial state whose momentum structure is representative of a midlatitude localized zonal wind anomaly will cause the free upper surface of the homogeneous atmosphere to rise (fall), and therefore the response can be viewed as being physically three-dimensional. The second model assumes that the troposphere can be represented by an unbounded continuously stratified Boussinesq fluid of constant Brunt-Vaisala frequency N, where the vertical gradient of the basic state potential temperature profile allows for the existence of vertically propagating internal inertia-gravity waves.

Lin, Yuh-Lang↗
NASA NTRS1985

The perturbed neutral circulation in the vicinity of a symmetric stable auroral arc

In the present investigation, a sophisticated high-resolution time-dependent numerical model is used to study the neutral response to changes in the E region electrodynamical forcing and particle heating associated with the appearance of a stable auroral arc. The initial condition is a nighttime ionosphere in the presence of a background precipitation of electrons associated with diffuse auroras. A simulation study is conducted of the perturbation in the neutral circulation. The simulations were performed on the basis of a model which describes nonlinear, nonhydrostatic viscous flow on a rotating atmosphere.

Walterscheid, R. L.↗
NASA NTRS1972

The effects of boundary condition asymmetries on the interplanetary magnetic field-moon interaction.

Boundary condition asymmetries inherent in the solar wind flow past the moon are included in a cylindrical model of the interplanetary magnetic field-moon interaction. Numerical examinations of the sunward side response of this model are compared in the frequency domain with those of symmetrically excited spherical and cylindrical models and two characteristic differences are observed: the response of the asymmetric model is depressed at low frequencies due to magnetic diffusion around a conducting core, and is flattened at high frequencies because of the finite application time of the incident interplanetary magnetic field. The diffusion of field lines around the core is also evident in the time response of the model in the antisolar cavity. The above features of the lunar response resulting from boundary condition asymmetries are shown to be evident in observational measurements.

Reisz, A. C.↗