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

Design model for the baffle aperture region of a hollow cathode thruster

A current balance on the main discharge chamber of an electron bombardment thruster is carried out. It is shown that the discharge current is the sum of the net cathode emission current, the beam current, and the ion current to cathode potential surfaces. Measurements of ion currents outside the virtual anode surface in a SERT II divergent magnetic field thruster show that relatively few ions cross this surface. This implies that ions are directed preferentially toward the anode and cathode pole pieces. A simple theoretical model useful as an aid in the design of the baffle aperture region of a hollow cathode equipped thruster is developed. The electron diffusion process through the aperture is modeled in accordance with the Bohm diffusion theory. The model is shown to yield consistent results for a given thruster geometry over substantial changes in operating conditions. The model's design usefulness is limited by a factor of two uncertainty observed over different thruster/cathode region geometries and by the accuracy to which the plasma parameters required as inputs to the model can be specified.

Brophy, J. R.↗
NASA NTRS1991

Turbulence and stall in plane diffusers - Computational study

The effect of free-stream turbulence level on stall conditions in plane diffusers is investigated with a finite-volume numerical method which utilizes a k-epsilon turbulence model or a differential Reynolds stress model for closure. Computed results show qualitative agreement with measured data. The ability to eliminate stall and thereby increase pressure recovery in diffusers at large included angles by utilizing high free-stream turbulence is correctly predicted.

Demuren, A. O.↗
NASA NTRS1995

Calculations of Diffuser Flows with an Anisotropic K-Epsilon Model

A newly developed anisotropic K-epsilon model is applied to calculate three axisymmetric diffuser flows with or without separation. The new model uses a quadratic stress-strain relation and satisfies the realizability conditions, i.e., it ensures both the positivity of the turbulent normal stresses and the Schwarz' inequality between any fluctuating velocities. Calculations are carried out with a finite-element method. A second-order accurate, bounded convection scheme and sufficiently fine grids are used to ensure numerical credibility of the solutions. The standard K-epsilon model is also used in order to highlight the performance of the new model. Comparison with the experimental data shows that the anisotropic K-epsilon model performs consistently better than does the standard K-epsilon model in all of the three test cases.

Zhu, J.↗
NASA NTRS2008

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Aero-Acoustic Analyses and Experimental Validation - Part 3

A unique assessment of acoustic similarity scaling laws and acoustic analogy methodologies in predicting the far-field acoustic signature from a sub-scale altitude rocket test facility at the NASA Stennis Space Center was performed. A directional, point-source similarity analysis was implemented for predicting the acoustic far-field. In this approach, experimental acoustic data obtained from "similar" rocket engine tests were appropriately scaled using key geometric and dynamic parameters. The accuracy of this engineering-level method is discussed by comparing the predictions with acoustic far-field measurements obtained. In addition, a CFD solver was coupled with a Lilley's acoustic analogy formulation to determine the improvement of using a physics-based methodology over an experimental correlation approach. In the current work, steady-state Reynolds-averaged Navier-Stokes calculations were used to model the internal flow of the rocket engine and altitude diffuser. These internal flow simulations provided the necessary realistic input conditions for external plume simulations. The CFD plume simulations were then used to provide the spatial turbulent noise source distributions in the acoustic analogy calculations. Preliminary findings of these studies will be discussed.

Allgood, Daniel C.↗
DOE OSTI2023

Simulations of self- and Xe diffusivity in uranium mononitride including chemistry and irradiation effects

A combination of density functional theory and empirical potential atomic scale simulations have been used to determine a model for defect stability and mobility in uranium mononitride (UN), as a function of temperature (T) and N 2 partial pressure (p N 2 ). Using the model, predictions of hypo-stoichiometry under U-rich conditions compare favorably to CALPHAD calculations using the TAF-ID database. Furthermore, our predictions of U and N self-diffusivity are in good agreement with experiments carried out as a function of T at specific partial pressures under thermal equilibrium. The validated atomic scale data have then been implemented within a cluster dynamics method to simulate irradiation-enhanced defect concentrations. All defects and clusters studied have significantly enhanced concentrations, with respect to thermal equilibrium, as T is lowered. The irradiation-enhanced Xe diffusivity is compared to post-irradiation annealing and in-pile experiments. In conclusion, the contributions of various defects and clusters to non-stoichiometry, self-diffusivity, and Xe diffusivity are discussed.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2016

Sorption Modeling and Verification for Off-Gas Treatment

This project was successfully executed to provide valuable adsorption data and improve a comprehensive model developed in previous work by the authors. Data obtained were used in an integrated computer program to predict the behavior of adsorption columns. The model is supported by experimental data and has been shown to predict capture of off gas similar to that evolving during the reprocessing of nuclear waste. The computer program structure contains (a) equilibrium models of off-gases with the adsorbate; (b) mass-transfer models to describe off-gas mass transfer to a particle, diffusion through the pores of the particle, and adsorption on the active sites of the particle; and (c) incorporation of these models into fixed bed adsorption modeling, which includes advection through the bed. These models are being connected with the MOOSE (Multiphysics Object-Oriented Simulation Environment) software developed at the Idaho National Laboratory through DGOSPREY (Discontinuous Galerkin Off-gas SeParation and REcoverY) computer codes developed in this project. Experiments for iodine and water adsorption have been conducted on reduced silver mordenite (Ag0Z) for single layered particles. Adsorption apparatuses have been constructed to execute these experiments over a useful range of conditions for temperatures ranging from ambient to 250°C and water dew points ranging from -69 to 19°C. Experimental results were analyzed to determine mass transfer and diffusion of these gases into the particles and to determine which models best describe the single and binary component mass transfer and diffusion processes. The experimental results were also used to demonstrate the capabilities of the comprehensive models developed to predict single-particle adsorption and transients of the adsorption-desorption processes in fixed beds. Models for adsorption and mass transfer have been developed to mathematically describe adsorption kinetics and transport via diffusion and advection processes. These models were built on a numerical framework for solving conservation law problems in one-dimensional geometries such as spheres, cylinders, and lines. Coupled with the framework are specific models for adsorption in commercial adsorbents, such as zeolites and mordenites. Utilizing this modeling approach, the authors were able to accurately describe and predict adsorption kinetic data obtained from experiments at a variety of different temperatures and gas phase concentrations. A demonstration of how these models, and framework, can be used to simulate adsorption in fixed- bed columns is provided. The CO 2 absorption work involved modeling with supportive experimental information. A dynamic model was developed to simulate CO 2 absorption using high alkaline content water solutions. The model is based upon transient mass and energy balances for chemical species commonly present in CO 2 absorption. A computer code was developed to implement CO 2 absorption with a chemical reaction model. Experiments were conducted in a laboratory scale column to determine the model parameters. The influence of geometric parameters and operating variables on CO 2 absorption was studied over a wide range of conditions. Continuing work could employ the model to control column operation and predict the absorption behavior under various input conditions and other prescribed experimental perturbations. The value of the validated models and numerical frameworks developed in this project is that they can be used to predict the sorption behavior of off-gas evolved during the reprocessing of nuclear waste and thus reduce the cost of the experiments. They can also be used to design sorption processes based on concentration limits and flow-rates determined at the plant level.

42 ENGINEERING↗
NASA NTRS1976

Post sunset behavior of the 6300 A atomic oxygen airglow emission

A theoretical model of the 6300 A OI airglow emission was developed based on the assumptions that both the charged and neutral portions of the Earth's upper atmosphere are in steady state conditions of diffusive equilibrium. Intensities of 6300 A OI emission line were calculated using electron density true height profiles from a standard C-4 ionosonde and exospheric temperatures derived from Fabry-Perot interferometer measurements of the Doppler broadened 6300 A emission line shape as inputs to the model. Reaction rate coefficient values, production mechanism efficiencies, solar radiation fluxes, absorption cross sections, and models of the neutral atmosphere were varied parametrically to establish a set of acceptable inputs which will consistently predict 6300 A emission intensities that closely agree with intensities observed during the post-sunset twilight period by an airglow observatory consisting of a Fabry-Perot interferometer and a turret photometer. Emission intensities that can only result from the dissociative recombination of molecular oxygen ions were observed during the latter portion of the observational period. Theoretical calculations indicate that contamination of the 6300 A OI emission should be on the order of or less than 3 percent; however, these results are very sensitive to the wavelengths of the individual lines and their intensities relative to the 6300 A OI intensity. This combination of a model atmosphere, production mechanism efficiencies, and quenching coefficient values was used when the dissociative photoexcitation and direct impact excitation processes were contributing to the intensity to establish best estimates of solar radiation fluxes in the Schumann--Runge continuum and associated absorption cross sections. Results show that the Jacchia 1971 model of the upper atmosphere combined with the Ackerman recommended solar radiation fluxes and associated absorption cross sections produces theoretically calculated intensities that more closely agree with the observed intensities than all the other combinations.

Smith, R. E.↗
DOE OSTI2021

Particle Size Distribution Dynamics Can Help Constrain the Phase State of Secondary Organic Aerosol

Particle phase state is a property of atmospheric aerosols that has important implications for the formation, evolution, and gas/particle partitioning of secondary organic aerosol (SOA). In this work, we use a size-resolved chemistry and microphysics model (SOM-TOMAS), updated to include an explicit treatment of particle phase state, to constrain the bulk diffusion coefficient (Db) of SOA produced from ??- pinene ozonolysis. By leveraging data from laboratory experiments performed in the absence of a seed and under dry conditions, we find that the Db for SOA can be constrained (1-7 ×10-15 cm2 s-1 in these experiments) by simultaneously reproducing the time-varying SOA mass concentrations and the evolution of the particle size distribution. Another version of our model that used the predicted SOA composition to calculate the glass transition temperature, viscosity, and, ultimately, Db (~10-15 cm2 s-1) of the SOA was able to reproduce the mass and size distribution measurements when we included oligomer formation (oligomers accounted for about a fifth of the SOA mass). Our work highlights the potential of a size resolved SOA model to constrain the particle phase state of SOA by utilizing historical measurements of the evolution of the particle size distribution.

He, Yicong↗
DOE OSTI2022

Modeling Ni redistribution in the hydrogen electrode of solid oxide cells through Ni(OH) 2 diffusion and Ni-YSZ wettability change

Ni redistribution in the hydrogen electrodes of solid oxide cells is an important degradation mechanism. Its driving force is still under debate. This work focuses on the interplay between Ni(OH) 2 diffusion and Ni-YSZ wettability change. With a Ni(OH) 2 diffusion model developed in our previous work, we further employ three models for the Ni-YSZ contact angle to quantify the Ni-YSZ wettability change. The microstructure evolutions in a reconstructed Ni-YSZ electrode are simulated under selected experimental operating conditions. It is shown that the phenomenological model can capture the driving force of Ni spreading/detachment on YSZ surface and the Ni migration reported in experiments can be qualitatively reproduced through the competition between Ni(OH) 2 diffusion and Ni-YSZ wettability change. It is also found that the initial microstructure of the active layer and the microstructure characteristics in the support layer can strongly affect the distribution of steam partial pressure and overpotential, hence the Ni redistribution, in the active layer, which may explain the inconsistencies in experiments. Further, our results show that there are two missing pieces in current theory of Ni redistribution: a mechanism of fast Ni diffusion under humid atmosphere and the physical reason behind the Ni spreading/detachment.

25 ENERGY STORAGE↗
NASA NTRS1996

A Computational Fluid Dynamic and Heat Transfer Model for Gaseous Core and Gas Cooled Space Power and Propulsion Reactors

A computational model based on the axisymmetric, thin-layer Navier-Stokes equations is developed to predict the convective, radiation and conductive heat transfer in high temperature space nuclear reactors. An implicit-explicit, finite volume, MacCormack method in conjunction with the Gauss-Seidel line iteration procedure is utilized to solve the thermal and fluid governing equations. Simulation of coolant and propellant flows in these reactors involves the subsonic and supersonic flows of hydrogen, helium and uranium tetrafluoride under variable boundary conditions. An enthalpy-rebalancing scheme is developed and implemented to enhance and accelerate the rate of convergence when a wall heat flux boundary condition is used. The model also incorporated the Baldwin and Lomax two-layer algebraic turbulence scheme for the calculation of the turbulent kinetic energy and eddy diffusivity of energy. The Rosseland diffusion approximation is used to simulate the radiative energy transfer in the optically thick environment of gas core reactors. The computational model is benchmarked with experimental data on flow separation angle and drag force acting on a suspended sphere in a cylindrical tube. The heat transfer is validated by comparing the computed results with the standard heat transfer correlations predictions. The model is used to simulate flow and heat transfer under a variety of design conditions. The effect of internal heat generation on the heat transfer in the gas core reactors is examined for a variety of power densities, 100 W/cc, 500 W/cc and 1000 W/cc. The maximum temperature, corresponding with the heat generation rates, are 2150 K, 2750 K and 3550 K, respectively. This analysis shows that the maximum temperature is strongly dependent on the value of heat generation rate. It also indicates that a heat generation rate higher than 1000 W/cc is necessary to maintain the gas temperature at about 3500 K, which is typical design temperature required to achieve high efficiency in the gas core reactors. The model is also used to predict the convective and radiation heat fluxes for the gas core reactors. The maximum value of heat flux occurs at the exit of the reactor core. Radiation heat flux increases with higher wall temperature. This behavior is due to the fact that the radiative heat flux is strongly dependent on wall temperature. This study also found that at temperature close to 3500 K the radiative heat flux is comparable with the convective heat flux in a uranium fluoride failed gas core reactor.

Anghaie, S.↗
DOE OSTI2016

Activation volumes of oxygen self-diffusion in fluorite structured oxides

In this study, fluorite structured oxides are used in numerous applications and as such it is necessary to determine their materials properties over a range of conditions. In the present study we employ molecular dynamics calculations to calculate the elastic and expansivity data, which are then used in a thermodynamic model (the cBΩ model) to calculate the activation volumes of oxygen self-diffusion coefficient in ThO 2 , UO 2 and PuO 2 fluorite structured oxides over a wide temperature range. We present relations to calculate the activation volumes of oxygen self-diffusion coefficient in ThO 2 , UO 2 and PuO 2 for a wide range of temperature (300–1700 K) and pressure (–7.5 to 7.5 GPa).

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Development of an Advanced Clothing Model that Considers Heat and Moisture Transport

Despite recent advancements in the development of human thermal physiology models, clothing is often considered in only a rudimentary manner, that is, as a boundary condition at the skin characterized by its bulk thermal and evaporative impedances. While such an approach may be sufficient for thermoregulation modeling of humans operating in relatively low thermal burden scenarios, e.g., within an office environment, a clothing thermal model that does not explicitly consider the layer-to-layer transport of vapor can yield incorrect results, especially under conditions corresponding to high thermal stress or extreme environments. For example, vapor diffusing through layers of clothing can reach its dew point, in which case a phase change will occur that can fundamentally affect the transport of moisture and heat within the clothing ensemble. Specifically, the evaporation and condensation of moisture within the fabric layers can affect liquid content and temperature and, ultimately, affect thermo-physiological response and thermal comfort perception. Development of an advanced clothing model that considers heat and moisture transport within the fabric layers can improve the accuracy of human thermal predictions in scenarios where clothing-based condensation and evaporation occur. This paper describes how a general heat transfer software modeling tool was modified to consider the transport of vapor, the storage of liquid, and the thermal effects of phase change within a clothing ensemble. The advanced clothing moisture model was tested by comparing its predictions to measurements derived from sweating hot plate experiments and human physiological experiments. The model was able to reproduce the thermal and moisture results from tests performed under conditions in which condensation within the clothing layers was known to occur. This study demonstrates the importance of considering layer-to-layer vapor transport in clothing models and its impact on human thermal modeling predictions.

Human thermal modeling↗
DOE OSTI2014

Exponential time differencing methods with Chebyshev collocation for polymers confined by interacting surfaces

We present a fast and accurate numerical method for the self-consistent field theory calculations of confined polymer systems. It introduces an exponential time differencing method (ETDRK4) based on Chebyshev collocation, which exhibits fourth-order accuracy in temporal domain and spectral accuracy in spatial domain, to solve the modified diffusion equations. Similar to the approach proposed by Hur et al. [Macromolecules 45, 2905 (2012)], non-periodic boundary conditions are adopted to model the confining walls with or without preferential interactions with polymer species, avoiding the use of surface field terms and the mask technique in a conventional approach. The performance of ETDRK4 is examined in comparison with the operator splitting methods with either Fourier collocation or Chebyshev collocation. Numerical experiments show that our exponential time differencing method is more efficient than the operator splitting methods in high accuracy calculations. This method has been applied to diblock copolymers confined by two parallel flat surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Controlling kinetic and diffusive length-scales during absorptive hydrogen removal in methane dehydroaromatization on MoC x /H-ZSM-5 catalysts

Addition of Zr metal absorbent to MoC x /H-ZSM-5 in the form of staged-bed, stratified-bed, and interpellet physical mixtures effectively scavenges H 2 from catalyst proximity, enhancing maximum single-pass benzene + naphthalene yield during methane dehydroaromatization (DHA) reactions to 14-16% compared to 8% in formulations without zirconium. The coupling of spatially-distinct catalytic and absorptive functions is achieved by dispersive/diffusive transport which conveys H 2 to staged Zr both co- and counter-current to bulk advection, thereby suppressing axial H 2 partial pressure profiles along the catalyst bed and enhancing net aromatization rates. We evince hitherto unreported significance of dispersive hydrogen transport during methane DHA by measurement of Péclet number, Pe = 1.32, in H 2 tracer studies with step-change or impulse input to inert catalyst proxies. Kinetic limits to methane pyrolysis are quantified by Damköhler number, Da, for synthesis of benzene, Da B = 0.15, and naphthalene, Da N = 0.03, determined from kinetic studies which rigorously account for reversibility of DHA reactions. Here, detailed reaction-transport models synthesize interplay of kinetic, diffusive, and convective length-scales captured by Péclet and Damköhler number to predict influence of catalyst-absorbent proximity and process flow-conditions on aromatization rates. Systematic control of catalyst bed-length, L, or linear flow velocity, u, predictably alters Pe and Da to effect improvements in methane conversion with and without Zr metal, corroborating results from simulation of the reaction-transport model.

03 NATURAL GAS↗
DOE OSTI2014

FY 2014 Sodium-Supercritical CO 2 Interactions in the SNAKE Experiment Facility

One appealing feature of the supercritical carbon dioxide Brayton cycle energy conversion system is the small footprint that the hardware requires, which is in part due to the use of compact diffusion-bonded heat exchangers such as Printed Circuit Heat Exchanger (PCHEs) or Hybrid Heat Exchangers (H 2 X) as the heat source heat exchanger (sodium-to-CO 2 ) as well as the recuperator and cooler modules. Although PCHEs and H 2 Xs have a high degree of structural integrity, the potential for leaks to develop between the sodium and CO 2 coolant channels in the heat exchanger cannot be ruled out, and this would lead to discharge of high pressure CO 2 into the secondary coolant circuit. Due to the robustness of the PCHE design, catastrophic failure leading to CO 2 jet blowdown into the secondary sodium loop is not deemed likely. Rather, small cracks (or micro-leaks) may develop in which CO 2 will bleed into the secondary system at a relatively low rate and chemically react with the sodium. The goal of the sodium-CO 2 interaction tests is to gain a fundamental understanding of sodium-CO 2 interactions under prototypical conditions of compact diffusion-bonded heat exchanger failure, a fundamental understanding of self-plugging if it occurs, and the development of one-dimensional phenomenological models for the interactions between high-pressure CO 2 issuing into liquid sodium from a micro-leak across a stainless steel pressure boundary. These models will be validated using experiment data. Therefore, an experiment program at Argonne was initiated in Fiscal Year 2010 to investigate the reaction behavior between sodium and CO 2 under micro-leak conditions. Several reports have described the facility scaling rationale and design. Assembly of the SNAKE (S-CO 2 , Na Kinetics Experiment) began in Fiscal Year 2011 and was completed in July 2012. The sodium inventory in this facility consists of approximately 44 lbs (20 kg; ~21 L/5.5 gal).

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2013

Description of the First Observed Sodium-C0 2 Reactions in the Sodium CO 2 Interaction Experiment (SNAKE)

One appealing feature of the supercritical carbon dioxide Brayton cycle energy conversion system is the small footprint that the hardware requires, which is in part due to the use of Printed Circuit Heat Exchangers (PCHEs) as the heat source heat exchanger (sodium-to-CO 2 ) as well as the recuperator and cooler modules. Although PCHEs have a high degree of structural integrity, the potential for leaks to develop between the sodium and CO 2 coolant channels in the secondary heat exchanger cannot be ruled out, and this would lead to discharge of high pressure CO 2 into the secondary coolant circuit. Due to the robustness of the PCHE design, catastrophic failure leading to CO 2 jet blowdown into the secondary sodium loop is not deemed likely. Rather, small cracks (or micro-leaks) may develop in which CO 2 will bleed into the secondary system at a relatively low rate and chemically react with the sodium. The goal of the sodium-CO 2 interaction tests is to gain a fundamental understanding of sodium-CO 2 interactions under prototypical conditions of compact diffusion-bonded heat exchanger failure, a fundamental understanding of self-plugging if it occurs, and the development of one-dimensional phenomenological models for the interactions between high-pressure CO 2 issuing into liquid sodium from a micro-leak across a stainless steel pressure boundary. These models will be validated using experiment data. Therefore, an experiment program at Argonne was initiated in Fiscal Year 2010 to investigate the reaction characteristics between sodium and CO 2 under micro-leak conditions. Several reports have described the facility scaling rationale and design. Assembly of the SNAKE (S-CO 2 , Na Kinetics Experiment) began in Fiscal Year 2011 and was completed in July 2012. Approximately 44 lbs (20 kg; ~21 L/5.5 gal) of sodium was transferred into the SNAKE dump tank from a drum of clean sodium in July 2012.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1994

Calculation of Propulsive Nozzle Flowfields in Multidiffusing Chemically Reacting Environments

An advanced engineering model has been developed to aid in the analysis and design of hydrogen/oxygen chemical rocket engines. The complete multispecies, chemically reacting and multidiffusing Navier-Stokes equations are modelled, including the Soret thermal diffusion and the Dufour energy transfer terms. In addition to the spectrum of multispecies aspects developed, the model developed in this study is also conservative in axisymmetric flow for both inviscid and viscous flow environments and the boundary conditions employ a viscous, chemically reacting, reference plane characteristics method. Demonstration cases are presented for a 1030:1 area ratio nozzle, a 25 lbf film cooled nozzle, and a transpiration cooled plug and spool rocket engine. The results indicate that the thrust coefficient predictions of the 1030:1 and the 25 lbf film cooled nozzle are within 0.2 to 0.5 percent, respectively, of experimental measurements when all of the chemical reaction and diffusion terms are considered. Further, the model's predictions agree very well with the heat transfer measurements made in all of the nozzle test cases. The Soret thermal diffusion term is demonstrated to have a significant effect on the predicted mass fraction of hydrogen along the wall of the nozzle in both the laminar flow 1030:1 nozzle and the turbulent flow plug and spool nozzle analysis cases performed. Further, the Soret term was shown to represent an important fraction of the diffusion fluxes occurring in a transpiration cooled rocket engine.

Kacynski, Kenneth John↗
DOE OSTI2016

Laboratory Experiments on Bentonite Samples: FY16 Progress

The primary goal of this study is to improve the understanding of U(VI) sorption and diffusion behavior in sodium-montmorillonite in order to support the development of realistic conceptual models describing these processes in performance assessment models while (1) accounting for potential changes in system conditions over time and space, (2) avoiding overly conservative transport predictions, and (3) using a minimum number of fitting parameters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗