Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffuse layer model”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Verification of Unstructured Grid Adaptation Components

Adaptive unstructured grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic grid adaptation mechanics. Known differences and previously unknown differences in grid adaptation components and their integrated processes are identified here for study. Unstructured grid adaptation tools are verified using analytic functions and the Code Comparison Principle. Three analytic functions with different smoothness properties are adapted to show the impact of smoothness on implementation differences. A scalar advection-diffusion problem with an analytic solution that models a boundary layer is adapted to test individual grid adaptation components. The scalar problems illustrate known differences in a grid adaptation component implementation and a previously unknown interaction between components. Laminar flow over a delta wing is verified with multiple, independent grid adaptation procedures to show consistent convergence to fine-grid forces and pitching moment.

Park, Michael A.↗
DOE OSTI2020

Surrogate Model Development of Spent Fuel Degradation for Repository Performance Assessment

In model simulations of deep geologic repositories, UO 2 fuel matrix degradation typically begins as soon as the waste package breaches and groundwater contacts the fuel surface. The initial degradation rate depends on the timing of these events, burnup of the fuel, temperature, and concentrations of dissolved reactants. Estimating the initial rate of degradation is fairly straightforward, but as UO 2 corrosion products precipitate on the fuel surface and the movement of dissolved species between the fuel surface and environment is impeded by the precipitated solids, the rate is more difficult to quantify. At that point, calculating the degradation rate becomes a reactive-transport problem in which a large number of equations must be solved by iteration for a large number of grid cells at each time step. The consequence is that repository simulations, which are already expensive, become much more expensive, especially when hundreds or thousands of waste packages breach. The Fuel Matrix Degradation (FMD) model is the process model of the Spent Fuel and Waste Science and Technology (SFWST) campaign of the US Department of Energy (DOE). It calculates spent fuel degradation rates as a function of radiolysis, redox reactions, electrochemical reactions, alteration layer growth, and diffusion of reactants through the alteration layer. Like other similar fuel degradation process models, it is a complicated model requiring a large number of calculations and iterations at each time step.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Electric Field Effects on Water and Ion Structure and Diffusion at the Orthoclase (001)–Water Interface

Understanding the electrochemical properties of mineral–water interfaces tends to rely upon electrical double layer (EDL) models, but these models are based on the assumption that electrostatic equilibrium is constantly maintained. In reality, interfacial reactions, ion diffusion, and their electrochemical signatures are based in nonequilibrium conditions of locally or globally imbalanced electrical fields where current EDL models have limited purview. In this work, we performed molecular dynamics (MD) simulations of the orthoclase (001) surface in contact with a 1 M NaCl aqueous solution under various electric fields, to explore the interplay between EDL structure and dynamics when perturbed by electric fields of different direction and strength, by confinement, and by different distributions of structural surface charge. The simulations showed that confinement between two opposing (001) surfaces led to the development of an induced field when the applied field was perpendicular to the surfaces and, as a result, to ionic diffusion coefficients that were independent of electric field strength. In contrast, when the applied field was parallel to the surfaces, confinement resulted in ionic diffusion coefficients that were more strongly dependent on the magnitude of the electric field than in bulk water. Differences in the density and distribution of aluminol groups on the two surfaces had a significant impact on how the interfacial structure and dynamics varied in the presence of an electric field. Notably, these differences resulted in an electro-osmotic flow with opposite directions at the two surfaces under parallel applied electric field. Overall, the MD simulations highlighted the importance of considering atomic-level structure and heterogeneities when developing models of the electrochemical properties of mineral–water interfaces.

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

Solutions for Reacting and Nonreacting Viscous Shock Layers with Multicomponent Diffusion and Mass Injection

Numerical solutions are presented for the viscous shocklayer equations where the chemistry is treated as being either frozen, equilibrium, or nonequilibrium. Also the effects of the diffusion model, surface catalyticity, and mass injection on surface transport and flow parameters are considered. The equilibrium calculations for air species using multicomponent: diffusion provide solutions previously unavailable. The viscous shock-layer equations are solved by using an implicit finite-difference scheme. The flow is treated as a mixture of inert and thermally perfect species. Also the flow is assumed to be in vibrational equilibrium. All calculations are for a 45 deg hyperboloid. The flight conditions are those for various altitudes and velocities in the earth's atmosphere. Data are presented showing the effects of the chemical models; diffusion models; surface catalyticity; and mass injection of air, water, and ablation products on heat transfer; skin friction; shock stand-off distance; wall pressure distribution; and tangential velocity, temperature, and species profiles.

Moss, J. N.↗
DOE OSTI2022

Assessment and Improvement of Fission Product Transport Predictions of Particle Fuel in BISON

The U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel cycle systems. This program has been providing engineering-scale support for the development of BISON, a high-fidelity and high-resolution fuel performance tool. This study was motivated by the need to incorporate more physics-based models in BISON in order to foster tri-structural isotropic (TRISO) applications. This document details the integration of new modeling capabilities in BISON, including (1) development of pyrolytic carbon (PyC) and silicon carbide (SiC) layer anisotropic thermal and mass transport capabilities, (2) verification of the mass diffusion solution in TRISO modeling, (3) calibration of fission product diffusivity using Advanced Gas Reactor (AGR) experiments, (4) improved fission product release modeling by developing compact diffusion modeling capabilities, and (5) documentation of accelerated failure analysis on the BISON website. Improvements made to the diffusion models and parameters were documented and validated against AGR-1 and -2 experiment data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1979

Effects of turbulence model selection on the prediction of complex aerodynamic flows

Numerical simulations of viscous transonic flow over a circular-arc airfoil and in a diffuser are described. The simulations are made with a new computer program designed to serve as a tool in the development of improved turbulence models for complex flows. The program incorporates zero-, one-, and two-equation eddy viscosity models and includes a variety of subsonic and supersonic boundary conditions. The airfoil flow contains a shock-separated boundary-layer interaction that has resisted previous attempts at simulation. The diffuser flow also contains a shock-boundary-layer interaction, which has not been simulated previously. Calculations using standard turbulence models, developed originally for incompressible unseparated flows, are described. Results indicate that although there are interesting differences in predictions between the various models, none of them predict the flows accurately. Suggestions for improved turbulence models are discussed.

Coakley, T. J.↗
DOE OSTI2021

Understanding Carbon Dioxide Transfer in Direct Methanol Fuel Cells Using a Pore-Scale Model

Abstract The gas flow of carbon dioxide from the catalyst layer (CL) through the microporous layer (MPL) and gas diffusion layer (GDL) has great impacts on the water and fuel management in direct methanol fuel cells (DMFCs). This work has developed a liquid–vapor two-phase model considering the counter flow of carbon dioxide gas, methanol, and water liquid solution in porous electrodes of DMFC. The model simulation includes the capillary pressure as well as the pressure drop due to flow resistance through the fuel cell components. The pressure drop of carbon dioxide flow is found to be about two to three orders of magnitude higher than the pressure drop of the liquid flow. The big difference between liquid and gas pressure drops can be explained by two reasons: volume flowrate of gas is three orders of magnitude higher than that of liquid; only a small fraction of pores (<5%) in hydrophilic fuel cell components are available for gas flow. Model results indicate that the gas pressure and the mass transfer resistance of liquid and gas are more sensitive to the pore size distribution than the thickness of porous components. To buildup high gas pressure and high mass transfer resistance of liquid, the MPL and CL should avoid micro-cracks during manufacture. Distributions of pore size and wettability of the GDL and MPL have been designed to reduce the methanol crossover and improve fuel efficiency. The model results provide design guidance to obtain superior DMFC performance using highly concentrated methanol solutions or even pure methanol.

Electrochemistry↗
NASA NTRS1983

The stability of ground ice in the equatorial region of Mars

The lifetime of an unreplenished layer of ground ice lying within 30 deg of the Martian equator was examined within the context of the existing data base on Martian regolith and climate. Data on the partial pressure of H2O in the Martian atmosphere and the range of mean annual temperatures indicated the ground ice would be restricted to latitudes poleward of 40 deg. However, the ground ice near the poles may be a relic from early Martian geologic times held in place by a thin layer of regolith. Consideration of twelve model pore size distributions, similar to silt- and clay-type earth soils, was combined with a parallel pore model of gaseous diffusion to calculate the flux of H2O molecules escaping from the subsurface ground ice layer. Martian equatorial ground ice was found to be influenced by the soil structure, the magnitude of the geothermal gradient, the climatic desorption of CO2 from the regolith. It is concluded that equatorial ground ice is present on Mars only if a process of replenishment is active.

Clifford, S. M.↗
DOE OSTI2022

Understanding the E balance for water management in hydroxide exchange membrane fuel cells

As a key challenge for hydroxide exchange membrane fuel cells (HEMFCs), understanding water management deeply is critical for performance and durability improvement. Here, in this work, we show that the balance between water generation and phase-change-induced (PCI) flow is the main factor controlling cell voltage over a wide range of conditions where the cells are vulnerable to severe flooding. The regulation of flooding by heat generation can be reflected by a stable voltage with increase of current density in the polarization curve. This PCI balance voltage (E balance ), which determines the shape of polarization curve, is modeled as a function of cell operating conditions and gas diffusion layer (GDL) properties. Temperature, pressure, and thermal conductivity of GDL were found important for tuning the level of E balance , both in HEMFCs and proton exchange membrane fuel cells. The understanding of E balance will help to generate strategies for resolving flooding in HEMFCs, and it offers a new perspective for polarization curve interpretation.

25 ENERGY STORAGE↗
NASA NTRS1987

Downstream boundary effects on the frequency of self-excited oscillations in transonic diffuser flows

An investigation of downstream boundary effects on the frequency of self-excited oscillations in two-dimensional, separated transonic diffuser flows has been conducted numerically by solving the compressible, Reynolds-averaged, thin-layer Navier-Stokes equation with a two-equation turbulence model. It was found that the unsteady diffuser flowfields are very sensitive to the location of the downstream boundary. Extension of the diffuser downstream boundary significantly reduces the frequency and amplitude of oscillations for pressure, velocity and shock. Computational results suggest that the mechanism causing the self-excited oscillation changes from viscous convective wave dominated oscillations to inviscid acoustic wave dominated oscillations when the location of downstream boundary varies from 8.66 to 134.7 throat height. The existence of a suction slot in the experimental setup obscures the physical downstream boundary and, therefore, presents a difficulty for quantitative comparisons between computation and experiment.

Hsieh, T.↗
NASA NTRS1977

Contribution to polar albedo from a mesospheric aerosol layer

An examination is made of the impact of a layer of particulate matter, assumed to be ice crystals, on the albedo of the polar region. The model is time dependent, includes the growth of the layer, and incorporates the diffuse nature of radiation reflected from the surface and atmosphere. Although the magnitude of the effect is about an order of magnitude less than previous results, the impact is one of heating instead of cooling. It is also shown that ignoring the diffuse nature of the radiation reflected from the underlying earth-atmosphere system, as has been done in many previous simple models, can result in overestimation of the climatological impact of aerosols in sign and magnitude by a factor of up to 4-6.

Hummel, J. R.↗
NASA NTRS2015

Vibroacoustic Tailoring of a Rod-Stiffened Composite Fuselage Panel with Multidisciplinary Considerations

An efficient multi-objective design tailoring procedure seeking to improve the vibroacoustic performance of a fuselage panel while maintaining or reducing weight is presented. The structure considered is the pultruded rod stitched efficient unitized structure, a highly integrated composite structure concept designed for a noncylindrical, next-generation flight vehicle fuselage. Modifications to a baseline design are evaluated within a six-parameter design space including spacing, flange width, and web height for both frame and stringer substructure components. The change in sound power radiation attributed to a design change is predicted using finite-element models sized and meshed for analyses in the 500 Hz, 1 kHz, and 2 kHz octave bands. Three design studies are carried out in parallel while considering a diffuse acoustic field excitation and two types of turbulent boundary-layer excitation. Kriging surrogate models are used to reduce the computational costs of resolving the vibroacoustic and weight objective Pareto fronts. The resulting Pareto optimal designs are then evaluated under a static pressurization ultimate load to assess structural strength and stability. Results suggest that choosing alternative configurations within the considered design space can reduce weight and improve vibroacoustic performance without compromising strength and stability of the structure under the static load condition considered, but the tradeoffs are significantly influenced by the spatial characteristics of the assumed excitation field.

Allen, Albert R.↗
NASA NTRSDate not supplied

Multiscale Modeling of the Uptake of Hydrogen Chloride on Anodized Aluminum in Relevance to Spacecraft Fire Safety

Fire safety is a critical component of space missions. An electrical fire within the spacecraft often produces gaseous hydrogen chloride (HCl) from the smoldering of wire insulation. These gases tend to stick to the interior materials of the spacecraft, which is often comprised of treated aluminum. Determining the rate of uptake for HCl by aluminum helps the safe design of the spacecraft to reduce catastrophic mission risk. Our previous studies showed that a one-step global reaction model for adsorption of HCl on aluminum is incapable of predicting the experimentally measured HCl uptake of anodized aluminum, compared to other treatments with thinner oxide layers. X-ray photoelectron spectroscopy revealed that aluminum chloride is found deep in the oxide layer of the anodized aluminum, implying that the transport of HCl within the oxide layer covering the aluminum surface must be considered. This work uses a multiscale modeling approach to predict the uptake of HCl. A pore model is developed wherein the HCl reacts with the aluminum oxide pore walls to create an aluminum chloride layer. The HCl then diffuses through that layer for further reaction, until the diffusion resistance grows too large, and the surface is said to be saturated. The pore model is then coupled with the reactor-scale model to predict overall HCl uptake. Results show that this new multiscale approach predicts the uptake of HCl far more accurately than previously used single-scale models.

Modeling↗
NASA NTRS1985

Electrical coupling effects on the temporal evolution of F layer plasma structure

A time dependent model of F region structure decay by 'classical' cross field diffusion and electrical coupling along magnetic field lines to the E region is examined. The temporal behavior of the ion concentration fluctuations is determined by the electric field in the coupled system as well as by the initial perturbation spectra and the E region recombination rate. The formation of image structure in the E region ion concentration affects the lifetime of F layer structure in a scale size dependent way. Once an image is formed, the image amplitude and the driving F region structure amplitude decay at the same rate. At large scale sizes of at least lambda(2pi/k), this rate is proportional to k2 and the ratio of the temperatures in each region. At small scale sizes it depends on the E region recombination rate and the temperatures of the two regions but is only very weakly dependent on k. The background E region concentration determines the wave number beyond which the structure amplitude decay rate is almost independent of its scale size.

Heelis, R. A.↗
DOE OSTI2022

Improved Representation of Horizontal Variability and Turbulence in Mesoscale Simulations of an Extended Cold-Air Pool Event

Abstract Cold-air pools (CAPs), or stable atmospheric boundary layers that form within topographic basins, are associated with poor air quality, hazardous weather, and low wind energy output. Accurate prediction of CAP dynamics presents a challenge for mesoscale forecast models in part because CAPs occur in regions of complex terrain, where traditional turbulence parameterizations may not be appropriate. This study examines the effects of the planetary boundary layer (PBL) scheme and horizontal diffusion treatment on CAP prediction in the Weather Research and Forecasting (WRF) Model. Model runs with a one-dimensional (1D) PBL scheme and Smagorinsky-like horizontal diffusion are compared with runs that use a new three-dimensional (3D) PBL scheme to calculate turbulent fluxes. Simulations are completed in a nested configuration with 3-km/750-m horizontal grid spacing over a 10-day case study in the Columbia River basin, and results are compared with observations from the Second Wind Forecast Improvement Project. Using event-averaged error metrics, potential temperature and wind speed errors are shown to decrease both with increased horizontal grid resolution and with improved treatment of horizontal diffusion over steep terrain. The 3D PBL scheme further reduces errors relative to a standard 1D PBL approach. Error reduction is accentuated during CAP erosion, when turbulent mixing plays a more dominant role in the dynamics. Last, the 3D PBL scheme is shown to reduce near-surface overestimates of turbulence kinetic energy during the CAP event. The sensitivity of turbulence predictions to the master length-scale formulation in the 3D PBL parameterization is also explored. Significance Statement In this article, we demonstrate how a new framework for modeling atmospheric turbulence improves cold pool predictions, using a case study from January 2017 in the Columbia River basin (U.S. Pacific Northwest). Cold pools are regions of cold, stagnant air that form within valleys or basins, and improved forecasts could help to mitigate the risks they pose to air quality, transportation, and wind energy production. For the chosen case study, our tests show a reduction in temperature and wind speed errors by up to a factor of 2–3 relative to standard model options. These results strongly motivate continued development of the framework as well as its application to other complex weather events.

17 WIND ENERGY↗
DOE OSTI2025

Ab Initio Design of High-Entropy Thermal/ Environmental Barrier Coatings

Next generation thermal/environmental barrier coatings (TEBC) require carefully balancing various properties including phase stability, thermal conductivity, coefficient of thermal expansion (CTE), mechanical properties, and resistance against hot corrosion and water vapor recession. This work mainly focuses on rapid design of cost-effective high entropy rare-earth disilicates and aluminum garnets to protect SiC-based ceramic matrix composites and nickel-based superalloys in the hot section of gas turbine engines using density functional theory methods. Our calculations identify several low-cost high entropy TEBC exhibiting ultralow thermal conductivity at 1500 K and desirable CTE while maintaining good mechanical properties, including Er1/2Y3/4Yb3/4Si2O7, Gd1/4Er1/4Y3/4Yb3/4Si2O7, Eu1/4Er1/4Y3/4Yb3/4Si2O7, and (Y1/4Gd1/4Er1/4Yb1/4)3Al5O12. This work also aims to gain fundamental understanding of oxygen diffusion in model disilicates. Minimizing oxidizer (such as water vapor and oxygen) permeability through the EBC layer can significantly decrease the growth rate of thermally grown oxide and extend the service life of the coating system. Oxygen diffusion mechanisms including formation energy of defects under varying oxygen conditions and defect migration energy barriers will be presented.

coefficient of thermal expansion↗
DOE OSTI2022

Analysis of Compressibility Effects and Nonlinear Property Variations in a Supercritical CO2 Mixing Layer

Calculations are performed using the large eddy simulation technique to qualitatively assess the effects of fluid compressibility and thermodynamic nonlinearities on the dynamics of a supercritical flow field. Resulting data are also used to interrogate instantaneous subfilter-scale fields, as modeled by the mixed-dynamic Smagorinsky closure and gradient diffusion model. A three-dimensional wall-resolved calculation of a spatially evolving mixing layer composed solely of carbon dioxide is performed. This is followed by a coarse grid calculation to examine model performance as a function of resolution. Theoretical analysis indicates that the partial derivatives of density with respect to pressure and temperature are important modulators of the temperature and pressure field evolution, respectively. Results pertaining to subfilter velocity-velocity and velocity-temperature fields indicate that the models are most stressed in the regions where compressibility and thermo-fluid nonlinearities dominate the physics.

42 ENGINEERING↗
NASA NTRS2019

A Multi-Physics Study on High-Specific Power Li-O2 Batteries for Electric Aircraft

Commercialization of lithium-air batteries faces many challenges, such as electrolyte decomposition, short cycle life, low energy and power density, etc. However, commercialization of Li-O2 batteries for aeronautics is much more challenging due to additional safety constraints on cyclability and performance (high specific power and specific energy). For this presentation, we will discuss inter-related aspects of physics-based modeling of a pack: cell and battery model calibration. In addition, we will evaluate and present optimal battery designs for high discharge current density, high discharge time, and low battery mass using simulation-based optimization.The Finite Element Model (FEM) used to simulate a Li-O2 cell is based on the work of Bevara [1]. The different aspects of the model are based on: porous electrode theory and concentrated electrolyte theory; quantum tunneling model for the resistance of conformal layer of discharge product (Li2O2) [1]; Butler-Volmer kinetics for electrochemical reaction; Fick's diffusion for oxygen transport; and an oxygen dissolution model is applied at the air/electrolyte interface [2]. The electrolyte properties such as ion conductivity, ion diffusion, oxygen diffusion, and mass density of the electrolyte were taken from Molecular Dynamics (MD) simulations [3]; while the other model parameters, which includes mass of cell components, were calibrated to match experiments at high discharge current densities. The cell mass includes the anode, cathode, separator, electrolyte, and other components (such as current collector). This calibrated model is used to perform parametric studies on cathode thickness, porosity, tortuosity, carbon particle size, electrolyte transport and material properties, partial pressure of oxygen, discharge time, and discharge current density to study optimal designs for high specific power and energy. References:1. Bevara, V. & Andrei, P. (2014), J. Electrochem. Soc. 161 (14), A2068-A2079.2.Mehta, M. & Andrei, P. (2015), J. Power Sources. 286, 299-308.3.Liyana-Arachchi, T.; Haskins, J.; Burke, C.; Diederichsen, K.; McCloskey, B.; & Lawson, J. (2018), J. Phys. Chem. B. 122 (36), 8548 - 8559.4.Choi, W.; Kikumoto, H.; Choudhary, R. & Ooka, R. (2018), Applied Energy, 209, 306-321.

Mehta, Mohit↗