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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Properties of the +70 kilometers per second cloud toward HD 203664

I present high-resolution International Ultraviolet Explorer (IUE) spectra of the ultraviolet absorption in an intermediate-velocity interstellar cloud (nu(sub LSR) approximately equal to +70 km/s) toward HD 203664. The combined, multiple IUE images result in spectra with S/N = 15-40 and resolutions of approximately 20-25 km/s. The intermediate-velocity cloud absorption is present in ultraviolet lines of C II, C II(sup *), C IV, N I, O I, Mg I, Mg II, Al II, Al III, Si II, Si III, Si IV, S II, Cr II, Mn II, Fe II, and Zn II. The relative abundances of low-ionization species suggest an electron density of 0.15-0.34/cu cm and a temperature of 5300-6100 K in the neutral and weakly ionized gas. Given the presence of high-ionization gas tracers such as Si IV and C IV, ionized portions of the cloud probably contribute to the relatively large values of n(sub e) derived from measurements of the lower ionization species. The high-ionization species in the cloud have an abundance ratio, N(C IV)/N(Si IV) approximately equal to 4.5, similar to that inferred for collisionally ionized cloud interfaces at temperatures near 10(exp 5) K along other sight lines. When referenced to sulfur, the abundances of most elements in the cloud are within a factor of 5 of their solar values, which suggests that the +70 km/s gas has a previous origin in the Galactic disk despite a recent determination by Little et al. that the cloud lies at a distance of 200-1500 pc below the Galactic plane. I have checked this result against a model of the ionization for the diffuse ionized gas layer of the Galaxy and find that this conclusion is essentially unchanged as long as the ionization parameter is low as implied by the abundances of adjoining ionization states of aluminum and silicon. The processes responsible for the production of highly ionized gas in the +70 km/s cloud appear to be able to account for the inferred dust grain destruction as well.

Sembach, Kenneth R.↗
DOE OSTI2021

Demystifying the Stern layer at a metal–electrolyte interface: Local dielectric constant, specific ion adsorption, and partial charge transfer

Electric double layer (EDL) represents one of the most basic concepts in electrochemistry and is pertinent to diverse engineering applications ranging from electrocatalysis to energy storage. Whereas phenomenological and coarse-grained models have been long established to describe ionic distributions in the diffuse layer, a faithful prediction of the physicochemical properties of the electrode–electrolyte interface from a molecular perspective remains a daunting challenge. In this work, we investigate the charging behavior of an Ag (111) electrode in NaF aqueous solutions leveraging experimental results and theoretical calculations based on the classical density functional theory for ion distributions in the diffuse layer and on the joint density functional theory (JDFT) for the electronic structure. When the Ag electrode is applied with a negative voltage, the surface charge density can be reasonably described by assuming a neutral Stern layer with the dielectric constant dependent on the local electric field as predicted by the Kirkwood equation. However, the specific adsorption of F− ions must be considered when the electrode is positively charged and the fluoride adsorption can be attributed to both physical and chemical interactions. Qualitatively, F− binding and partial charge transfer are supported by JDFT calculations, which predict an increased binding energy as the voltage increases. Our findings shed insight on the molecular characteristics of the Stern layer and the charge behavior of adsorbed species not specified by conventional EDL models.

Wang, Xuepeng (ORCID:0000000285704546)↗
NASA NTRS2010

Multiple Scattering of Laser Pulses in Snow Over Ice: Modeling the Potential Bias in ICESat Altimetry

The primary goal of NASA's current ICESat and future ICESat2 missions is to map the altitude of the Earth's land ice with high accuracy using laser altimetry technology, and to measure sea ice freeboard. Ice however is a highly transparent optical medium with variable scattering and absorption properties. Moreover, it is often covered by a layer of snow with varying depth and optical properties largely dependent on its age. We describe a modeling framework for estimating the potential altimetry bias caused by multiple scattering in the layered medium. We use both a Monte Carlo technique and an analytical diffusion model valid for optically thick media. Our preliminary numerical results are consistent with estimates of the multiple scattering delay from laboratory measurements using snow harvested in Greenland, namely, a few cm. Planned refinements of the models are described.

Davis, A. B.↗
NASA NTRS2019

ACC Subduction by Mesoscales

The mesoscale contribution to Subduction in the Southern Ocean was recently studied by Salle'e and Rintoul (2011, SR) using the following mesoscale model. The adiabatic A-regime was modeled with the GM stream function, the diabatic D-regime was modeled with tapering functions,the D-A interface was taken to be at the mixed layer depth and the mesoscale diffusivity was either a constant or given by a 2D model. Since the resulting subductions were an order of magnitude smaller than the data of ±200 m/yr (Mazloff et al., 2010), SR showed that if instead of the above model-dependent mesoscale diffusivities, they employed the ones by Salle'e et al. (2008) from surface drifter observations, the subductions compared significantly better with the data. On those grounds, SR suggested a tenfold increase of the diffusivity. In this work, we suggest that since the mesoscale diffusivity is but one component of a muc large mesoscale parameterization, one should first assess the latter's overall performance followed by the assessment of the predicted ACC subduction. We employ the mesoscale model formulated31 in Canuto et al. (2018; 2019, that includes recent theoretical and observational advances and that was assessed against a variety of data including the output of 17 other OGCMs. The ACC diffusivities compare well with drifter data by Salle'e et al. (2008) and the ACC subduction rates are in agreement with the data.

V M Canuto↗
NASA NTRS1982

Three dimensional flow field inside compressor rotor, including blade boundary layers

The Reynolds stress equation, pressure strain correlation, and dissipative terms and diffusion are discussed in relation to turbulence modelling using the Reynolds stress model. Algebraic modeling of Reynolds stresses and calculation of the boundary layer over an axial cylinder are examined with regards to the kinetic energy model for turbulence modelling. The numerical analysis of blade and hub wall boundary layers, and an experimental study of rotor blade boundary layer in an axial flow compressor rotor are discussed. The Patankar-Spalding numerical method for two dimensional boundary layers is included.

Galmes, J. M.↗
NASA NTRS1984

Analytic modeling of a spray diffusion flame

A detailed model for a spray diffusion flame is described. The model is based on the boundary layer form of the equations of motion, with droplet transport accounted for using a discretized droplet size distribution function. Interphase transport of mass and energy are accounted for, with a flame-sheet model used to describe the combustion process on a droplet scale. Near dynamic equilibrium is assumed for the description of droplet transport; droplets can diffuse relative to the gas phase. Gas-phase mixing is accounted for using a two-equation turbulence model; buoyancy effects are included, with a temperature fluctuation equation used to account for buoyancy effects on turbulence structure. Thermal radiation from gas-phase CO2 and H2O is included. Gas-phase chemical kinetics are modeled using a 20-reaction, 10-species version of the advanced quasi-global chemical kinetics formulation. Results are compared with data for a vaporizing Freon spray and a pentane spray flame. It is shown that the computational approach provides a reasonably valid picture of the overall development of a spray diffusion flame, and, furthermore, provides a useful tool for the parametric examination of the spray combustion process.

Harsha, P. T.↗
NASA NTRS2017

Evaluation of Full Reynolds Stress Turbulence Models in FUN3D

Full seven-equation Reynolds stress turbulence models are a relatively new and promising tool for todays aerospace technology challenges. This paper uses two stress-omega full Reynolds stress models to evaluate challenging flows including shock-wave boundary layer interactions, separation and mixing layers. The Wilcox and the SSGLRR full second-moment Reynolds stress models are evaluated for four problems: a transonic two-dimensional diffuser, a supersonic axisymmetric compression corner, a compressible planar shear layer, and a subsonic axisymmetric jet. Simulation results are compared with experimental data and results using the more commonly used Spalart-Allmaras (SA) one-equation and the Menter Shear Stress Transport (SST) two-equation models.

computational fluid dynamics↗
NASA NTRS1976

Wind enhanced planetary escape: Collisional modifications

The problem of thermal escape is considered in which both the effects of thermospheric winds at the exobase and collisions below the exobase are included in a Monte Carlo calculation. The collisions are included by means of a collisional relaxation layer of a background gas which models the transition region between the exosphere and the thermosphere. The wind effects are considered in the limiting cases of vertical and horizontal flows. Two species are considered: terrestrial hydrogen and terrestrial helium. In the cases of terrestrial hydrogen the escape fluxes were found to be strongly filtered or throttled by collisions at high exospheric temperatures. The model is applied to molecular hydrogen diffusing through a methane relaxation layer under conditions possible on Titan. The results are similar to the case of terrestrial hydrogen with wind enhanced escape being strongly suppressed by collisions. It is concluded that wind enhanced escape is not an important process on Titan.

Curtis, S. A.↗
DOE OSTI2025

Comparative Evaluation of the Ability of the MYNN‐EDMF PBL Scheme in WRF Model to Reproduce Near Surface Wind Speed Over Different Topographical Types

Abstract This study systematically evaluates the performance of the Mellor‐Yamada‐Nakanishi‐Niino‐Eddy‐Diffusion‐Mass‐Flux planetary boundary layer (PBL) scheme within the Weather Research and Forecasting (WRF) model in simulating near‐surface wind speeds across various topographies in New York State (NYS). Simulated wind speeds are compared with in‐situ measurements from 22 surface sites, grouped into six topographic categories: continental plain (CT), lakeside (LS), river valley (RV), Long Island (LI), Block Island (BI), and offshore ocean (OO). A quantitative evaluation based on Relative Euclidean Distance shows that wind speeds at the OO site are the most accurately reproduced, followed by those at LI sites, while the model performs less accurately for the remaining topographic groups. Wind speeds over CT sites tend to be overestimated by approximately 1 m/s, although their diurnal variability (DV) is well captured. In contrast, the model underestimates wind DV at LS, RV, LI, and BI sites, with the largest biases occurring at LI and BI, resulting in underestimated daytime wind speed and/or overestimated nighttime wind speed. The OO winds exhibit minimal diurnal variation, accurately captured by our WRF model. The surface wind diurnal variation is closely linked to PBL development. Among the indicators of PBL development, surface potential temperature biases most strongly correlate with wind speed biases. Our WRF model faces challenges in capturing the distinctions between winds influenced by local circulations and those over continental plains, and the significantly stronger winds at OO compared to BI. Potential causes for these biases are discussed, offering pathways for improving surface wind simulations in future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The nanoscale distribution of copper and its influence on charge collection in CdTe solar cells

For decades, Cu has been the primary dopant in CdTe photovoltaic absorbers. Typically, Cu acceptor concentrations in these devices are on the order of 1 1014 cm-3, which has made it notoriously difficult to directly correlate nanoscale Cu distributions to the local charge transport properties of these devices. To measure and correlate these properties, measurement techniques require high sensitivity to elemental concentration, large penetration depth, and operando compatibility. Techniques such as secondary-ion mass spectroscopy and X-ray energy dispersive spectroscopy are widely adopted to measure Cu concentrations, but they are limited by penetration depth, sensitivity, or spatial resolution. Additionally, they lack the operando capabilities required to correlate one-to-one Cu concentrations to electrical performance. In this work, correlative X-ray microscopy is used to investigate the spatial distribution of Cu and its impact on charge collection through the depth and breadth of CdTe photovoltaic devices. Plan-view, nanoscale X-ray fluorescence maps clearly demonstrate the spatial segregation of copper around regions thought to be CdTe grain boundaries. Complementary cross-section imaging unveils the transition of the maximum charge-collection efficiency from the ZnTe-CdTe interface to the CdS-CdTe interface as a function of Cu incorporation. The copper concentration through the depth of the CdTe layer is characterized by slow and fast diffusion components, and cross-section charge-transport modeling shows that the experimentally obtained charge collection can be explained by the modeled acceptor distribution through the depth of the CdTe layer.

14 SOLAR ENERGY↗
DOE OSTI2025

Multiscale, mechanistic modeling of cesium transport in silicon carbide for TRISO fuel performance prediction

Understanding cesium (Cs) transport in TRistructural ISOtropic (TRISO) particle fuel is crucial for predicting fission product release in high-temperature reactors. However, current challenges include significant scatter in diffusivity data and unexplained temperature-dependent diffusion regimes in the silicon carbide layer. This study addresses these challenges by developing a multiscale, mechanistic Cs transport model integrating atomistic simulations and phase field modeling. Our model quantifies temperature and grain size effects on Cs diffusivity, attributing experimentally observed regimes to a transition from bulk-dominated diffusivity at high temperatures to grain boundary-dominated diffusivity at lower temperatures. The model, validated against diffusion measurements and advanced gas reactor (AGR)-1 and AGR-2 post-irradiation fission product release data, enhances the predictive capability of the BISON fuel performance code. This study advances our understanding of Cs release from TRISO particles and its dependence on temperature and silicon carbide grain size, with implications for the safety and efficiency of high-temperature nuclear reactors.

BISON↗
DOE OSTI2023

Time dependent chlorination of CeO 2 , La 2 O 3 and Nd 2 O 3 by ZrCl 4 dissolved in eutectic LiCl–KCl

Partial chlorination of CeO 2 , La 2 O 3 , and Nd 2 O 3 was achieved using ZrCl 4 dissolved in molten LiCl-KCl at 450°C over a period of 4-5 hr. ZrCl 4 was pre-volatilized into an argon stream to maintain a constant vapor pressure for equilibration with the salt. This approach was effective at maintaining the ZrCl 4 concentration at the solubility limit of 1.3 mol% in the molten salt throughout the duration of two of the chlorination experiments and within about 60% of the solubility limit in the third experiment. Based on the analysis of salt samples, fractional conversion after four hours of reaction for CeO 2 , and five hours for La 2 O 3 , and Nd 2 O 3 were calculated to be 0.97, 0.28, and 0.65 respectively. Assuming a shrinking core diffusion-limited model, the diffusivity of ZrCl 4 through the hypothesized ZrO 2 layer was calculated for the chlorination of CeO 2 , La 2 O 3 , and Nd 2 O 3 were 1.76 x 10 -10 cm 2 /sec, 7.85 x 10 -12 cm 2 /sec, and 7.42 x 10 -12 cm 2 /sec, respectively. Given the extremely low values calculated for diffusivity, it was concluded that the process is chemical reaction rate limited and that a complete reaction of each rare earth oxide could be achieved. Furthermore, a chemically controlled shrinking core model was shown to also fit the experimental data closely.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1991

Diurnal forcing of planetary atmospheres

A free convection parameterization has been introduced into the Mars Planetary Boundary Layer Model (MPBL). Previously, the model would fail to generate turbulence under conditions of zero wind shear, even when statically unstable. This in turn resulted in erroneous results at the equator, for example, when the lack of Coriolis forcing allowed zero wind conditions. The underlying cause of these failures was the level 2 second-order turbulence closure scheme which derived diffusivities as algebraic functions of the Richardson number (the ratio of static stability to wind shear). In the previous formulation, the diffusivities were scaled by the wind shear--a convenient parameter since it is non-negative. This was the drawback that all diffusivities are zero under conditions of zero shear (viz., the free convection case). The new scheme tests for the condition of zero shear in conjunction with static instability and recalculates the diffusivities using a static stability scaling. The results for a simulation of the equatorial boundary layer at autumnal equinox are presented. (Note that after some wind shear is generated, the model reverts to the traditional diffusivity calculation.)

Houben, Howard C.↗
DOE OSTI2025

Control of Two Solid Electrolyte Interphases at the Negative Electrode of an Anode‐Free All Solid‐State Battery based on Argyrodite Electrolyte

Abstract Anode‐free all solid‐state batteries (AF‐ASSBs) employ “empty” current collector with three active interfaces that determine electrochemical stability; lithium metal – Solid electrolyte (SE) interphase (SEI‐1), lithium – current collector interface, and collector – SE interphase (SEI‐2). Argyrodite Li 6 PS 5 Cl (LPSCl) solid electrolyte (SE) displays SEI‐2 containing copper sulfides, formed even at open circuit. Bilayer of 140 nm magnesium/30 nm tungsten (Mg/W‐Cu) controls the three interfaces and allows for state‐of‐the‐art electrochemical performance in half‐cells and fullcells. AF‐ASSB with NMC811 cathode achieves 150 cycles with Coulombic efficiency (CE) above 99.8%. With high mass‐loading cathode (8.6 mAh cm −2 ), AF‐ASSB retains 86.5% capacity after 45 cycles at 0.2C. During electrodeposition of Li, gradient Li‐Mg solid solution is formed, which reverses upon electrodissolution. This promotes conformal wetting/dewetting by Li and stabilizes SEI‐1 by lowering thermodynamic driving force for SE reduction. Inert refractory W underlayer is required to prevent ongoing formation of SEI‐2 that also drives electrochemical degradation. Inert Mo and Nb layers likewise protect Cu from corroding, while Li‐alloying layers (Mg, Sn) are less effective due to ongoing volume changes and associated pulverization. Mechanistic explanation for observed Li segregation within alloying Li x Mg layer is provided through mesoscale modelling, considering opposing roles of diffusivity differences and interfacial stresses.

Wang, Yixian [Materials Science and Engineering Pr↗
DOE OSTI2023

Disorder and hydrogenation in graphene nanopowder revealed by complementary X-ray and neutron scattering

Functionalization of graphene's two-dimensional sheets can be used to modify the physical properties of graphene, such as changing its conductivity or inducing ferromagnetism, which is of broad interest for a myriad of applications. However, functionalized graphene can possess significant structural disorder that is often overlooked and is difficult to characterize quantitatively. Here we investigate hydrogenated graphene nanopowder (H-Gr) produced by Birch reduction of graphene oxide. By combining complementary x-ray and neutron diffraction, we show how to quantitatively characterize the H-Gr nanostructure, including the hydrogen content. Further, we show that the majority of the H-Gr consists of highly disordered molecular-scale carbon while a small portion of the sample contains few-layered graphene having an expanded interlayer spacing. Modeling the coherent diffuse scattering for both x-ray and neutron diffraction, and comparing it with the incoherent neutron scattering, the hydrogen to carbon ratio was measured and it was determined that most of the hydrogen resides within the disordered carbon rather than in the graphene. Evidence is presented for hydrogen clustering as well as for hydrogen bonded perpendicular to the plane of molecular-carbon bonds. Our quantitative approach to characterizing carbon nanopowders has broad relevance to understanding disorder and bonding of hydrogen in graphitic carbons and other nanomaterials.

36 MATERIALS SCIENCE↗
NASA NTRS1986

A three-dimensional axisymmetric calculation procedure for turbulent flows in a radial vaneless diffuser

An analytical model is proposed to calculate the three-dimensional axisymmetric turbulent flowfield in a radial vneless diffuser. The model assumes that the radial and tangential boundary layer profiles are approximated by power law profiles. Then, using the integrated radial and tangential momentum and continuity equations for the boundary layer and corresponding inviscid equations for the core flow, there results six ordinary differential equation in six unknowns which are easily solved using a Runge-Kutta technique. A model is also proposed for fully developed flow. The results using this technique were compared with the results from a three-dimensional viscous, axisymmetric duct code and with experimental data and good quantitative agreement was obtained.

Schumann, L. F.↗
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. Laminar flow over a delta wing and turbulent flow over an ONERA M6 wing are verified with multiple, independent grid adaptation procedures to show consistent convergence to fine-grid forces and a moment. The scalar problems illustrate known differences in a grid adaptation component implementation and a previously unknown interaction between components. The wing adaptation cases in the current study document a clear improvement to existing grid adaptation procedures. The stage is set for the infusion of verified grid adaptation into production fluid flow simulations.

Park, Michael A.↗
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.↗