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 217 records · Page 12

Study of turbulence-induced refraction of lower hybrid waves using synthetic scrape-off layer filaments

Here, turbulence-induced refraction effects to lower hybrid (LH) wave propagation and current drive are studied using synthetic scrape-off layer (SOL) blob/filament fields. A synthetic 3D, field-following, blob turbulence model is implemented in the ray-tracing/Fokker-Planck (RTFP) codes GENRAY/CQL3D. In Alcator C-Mod, the blob field is shown to significantly affect LH ray-trajectories, leading to increased on-axis damping and smoother current profiles. This effect depends on the average blob size and amplitude. In addition, the diffusion of ray-trajectories in phase-space caused by turbulence increases the robustness of the RTFP model. A modified $N_{||}$ launch spectrum, acting as a proxy for parametric decay instability (PDI) effects, is included in simulations with the blob model. A synergy between the modified launch spectrum and turbulence-induced refraction results in synthetic hard x-ray profiles that agree with experiment. Lastly, the blob model is used to predict the effect of SOL turbulence on DIII-D high-field side (HFS) LH launch. Assuming low turbulence amplitude in the HFS SOL (~5%), turbulence-induced refraction is predicted to have little effect on current drive efficiency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Particle-in-cell modeling of electron beam generated plasma

Plasmas generated using energetic electron beams are well known for their low electron temperature (T e ) and plasma potential, which makes them attractive for atomic-precision plasma processing applications such as atomic layer etch and deposition. A 2-dimensional particle-in-cell model for an electron beam-generated plasma in argon confined by a constant applied magnetic field is described here in this article. Plasma production primarily occurs in the path of the beam electrons in the center of the chamber. The resulting plasma spreads out in the chamber through non-ambipolar diffusion with a short-circuiting effect allowing unequal electron and ion fluxes to different regions of the bounding conductive chamber walls. The cross-field transport of the electrons (and thus the steady-state characteristics of the plasma) are strongly impacted by the magnetic field. T e is anisotropic in the electron beam region, but low and isotropic away from the plasma production zone. The plasma density increases and the plasma becomes more confined near the region of production when the magnetic field strengthens. The magnetic field reduces both electron physical and energy transport perpendicular to the magnetic field. T e is uniform along the magnetic field lines and slowly decreases perpendicular to it. Electrons are less energetic in the sheath regions where the sheath electric field repels and confines the low-energy electrons from the bulk plasma. Even though electron and ion densities are similar in the bulk plasma due to quasi-neutrality, electron and ion fluxes on the grounded chamber walls are unequal at most locations. Electron confinement by the magnetic field weakens with increasing pressure, and the plasma spread out farther from the electron beam region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fabrication of high-performance gas-diffusion-electrode based membrane-electrode assemblies

This paper demonstrates the fabrication and processing steps required to produce high performance fuel cell membrane electrode assemblies (MEAs) based on spray-coated gas-diffusion electrodes (GDEs). It is demonstrated that coating the catalyst layer with a thin layer of ionomer and then hot pressing the GDEs to the membrane is required to achieve comparable catalyst activity and air performance to catalyst-coated-membrane MEAs. We show that there is a critical amount of ionomer required to achieve maximized performance. Using electron microscopy, we show that the combination of the ionomer overlayer and hot-pressing bonds the catalyst layer to the membrane, increasing the interfacial contact area and quality of this interface. We also find that the ionomer overlayer smooths the surface of the GDE and provides increased contact area between the GDE and the membrane. Additionally, we demonstrate that much less ionomer is required for high-performance than has been previously reported. Through model fitting of electrochemical impedance spectroscopy, we show that this improvement in the catalyst layer - membrane interface reduces the effective catalyst layer resistance, which reduces Ohmic losses and increases catalyst utilization.

33 ADVANCED PROPULSION SYSTEMS↗

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↗

Evaluation of Model Summertime Boundary Layer Cloud Development over Complex Terrain in New York State

Weather forecasting over complex terrain with diverse land cover is challenging. Utilizing the high-resolution observations from New York State Mesonet (NYSM), we are able to evaluate the surface processes of the Weather Research Forecast (WRF) Model in a detailed, scale-dependent manner. Here, in the study, possible impacts of land–atmosphere interaction on surface meteorology and boundary layer cloud development are investigated with different model resolutions, land surface models (LSMs), and planetary boundary layer (PBL) physical parameterizations. The High-Resolution Rapid Refresh, version 3 (HRRR), forecasting model is used as a reference for the sensitivity evaluation. Results show that over complex terrain, the high-resolution simulations (1 km × 60 vertical levels) generally perform better compared to low-resolution (3 km × 50 levels) in both surface meteorology and cloud fields. LSMs play a more important role in surface meteorology compared to PBL schemes. The NoahMP land surface model exhibits daytime warmer and drier biases compared to the Rapid Update Cycle (RUC) due to better prediction of the Bowen ratio in RUC. The PBL schemes would affect the convective strength in the boundary layer. The Shin–Hong (SH) scale-aware scheme tends to produce the strongest convective strength in the PBL, while the ACM2 PBL scheme rarely resolved convection even at 1-km resolution. By considering the radiation effect of subgrid-scale (SGS) clouds, the Mellor–Yamada–Nakanishi–Niino eddy diffusivity mass flux (MYNN-EDMF) predicted the highest cloud coverage and lowest surface solar radiation bias. The configuration of SGS clouds in MYNN-EDMF would not only significantly reduce shortwave radiation bias, but also affect the convection behaviors through land surface–cloud–radiation interaction.

54 ENVIRONMENTAL SCIENCES↗

Local wave-number model for inhomogeneous two-fluid mixing

We analyze the local wave-number (LWN) model, a two-point spectral closure model for turbulence, as applied to the Rayleigh-Taylor (RT) instability, the flow induced by the relaxation of a statically-unstable density stratification. Model outcomes are validated against data from 3D simulations of the RT instability. In the first part of the study we consider the minimal model terms required to capture inhomogeneous mixing and show that this version, with suitable model coefficients, is sufficient to capture the evolution of important mean global quantities including mix-width, turbulent mass flux velocity, and Reynolds stress, if the start time is chosen such that the earliest transitions are avoided. However, this simple model does not permit the expected finite asymptote of the density-specific-volume covariance b. In the second part of the study, we investigate two forms for a source term for the evolution of the spectrum of density-specific-volume covariance for the LWN model. The first includes an empirically motivated calibration of the source to achieve the final asymptotic state of constant b. Here, the second form does not require calibration but, in conjunction with enhanced diffusion and drag captures the full evolution of all the dynamical quantities, namely, the mix-layer growth, turbulent mass-flux velocity, Reynolds stress, as well as the desired behavior of b.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The structure, composition, and performance impact of a YSZ-GDC interdiffusion layer in solid oxide electrolysis cells

This study provides a combined experimental and computational investigation into the structure and impact of the cation interdiffusion layer that appears at the gadolinium doped ceria (GDC)/yttria stabilized zirconia (YSZ) interface in solid oxide electrolysis cells (SOECs). Scanning transmission electron microscopy (STEM) illustrates that a ∼0.4 μm interdiffusion layer (IDL) with an intermixed cation distribution and fine grain size forms upon sintering. STEM identifies that the interdiffusion layer exists in the cubic fluorite structure despite changes in cation composition. The interdiffusion layer microstructure formed during sintering does not change during SOEC testing at either 1.3V or heightened voltage pulse testing. Modeling predicts that ionic conductivity may decrease in the interdiffusion layer due to Coulombic trapping between mobile oxygen vacancies and excess Gd 3+ acceptor dopants. Yet, the density and continuous nature of the layer should benefit cell stability by substantially reducing the formation of SrZrO 3 , which is corroborated by STEM and Synchrotron X-ray diffraction (XRD). We conclude that the interdiffusion layer acts as a beneficial barrier to Sr diffusion, when operating in a regime where electrolyte void formation is not observed.

organic↗

Study of interface reaction in a B 4 C/Cr mirror at elevated temperature using soft X-ray reflectivity

Boron carbide is a prominent material for high-brilliance synchrotron optics as it remains stable up to very high temperatures. The present study shows a significant change taking place at 550°C in the buried interface region formed between the Cr adhesive layer and the native oxide layer present on the silicon substrate. An in situ annealing study is carried out at the Indus-1 Reflectivity beamline from room temperature to 550°C (100°C steps). The studied sample is a mirror-like boron carbide thin film of 400 Å thickness deposited with an adhesive layer of 20 Å Cr on a silicon substrate. The corresponding changes in the film structure are recorded using angle-dependent soft X-ray reflectivity measurements carried out in the region of the boron K-edge after each annealing temperature. Analyses performed using the Parratt recursive formalism reveal that the top boron carbide layer remains intact but interface reactions take place in the buried Cr–SiO2 region. After 300°C the Cr layer diffuses towards the substrate. At higher temperatures of 500°C and 550°C the Cr reacts with the native oxide layer and tends to form a low-density compound of chromium oxysilicide (CrSiOx). Depth profiling of Si and Cr distributions obtained from secondary ion mass spectroscopy measurements corroborate the layer model obtained from the soft X-ray reflectivity analyses. Details of the interface reaction taking place near the substrate region of boron carbide/Cr sample are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

15.3% AM1.5G Efficiency GaAs Solar Cells Fabricated via an Epitaxy-Free Process

Here, we report simple and potentially low-cost techniques for creating high-quality n-type gallium arsenide (GaAs) and GaAs p/n junctions and fabricate GaAs p/n junction solar cells. Detailed-balance modeling suggests that 20% AM1.5G efficiency p/n homojunction devices may be possible if the surface doping concentration can be limited to values less than ∼ 5 × 10 19 cm −3 . Our process exploits an open-tube, vapor-phase, deposition-free, zinc diffusion technique for forming p-type layers in melt-grown n-GaAs substrates that results in sheet resistances less than 1 kΩ/$\square$. In addition, we have improved the minority carrier diffusion lengths of melt-grown GaAs from less than one micron to over five microns using an open-tube, vacuum-free, annealing process which reduces the density of EL2 midgap defects. Finally, we have combined these advances to fabricate epitaxy-free, GaAs solar cells with a validated AM1.5G efficiency of 15.3%.

14 SOLAR ENERGY↗

Influence of carrier trapping on radiation detection properties in $\mathrm{CVD}$ grown 4H-SiC epitaxial layers with varying thickness up to 250 µm

We report low doped n-type 4H-SiC epitaxial layers of thickness 50, 150, and 250 µm grown by hot wall chemical vapor deposition were used to fabricate Ni/4H-SiC Schottky radiation detectors. The epitaxial layers were grown on the (0001) face of highly conductive bulk 4H-SiC substrates 8° offcut towards the $\langle 11\bar{2}1\rangle$ direction. The 50, 150, and 250 µm thick epilayer detectors, under optimized settings, showed energy resolutions of 2.0%, 0.78%, and 0.63%, respectively for 5.48 MeV alpha particles. Deep level transient spectroscopy studies showed that the observed variation in the detector resolution is linked to the defect parameters in the devices. Least-squares fitting of the bias dependence of the charge collection efficiency according to a drift-diffusion model, revealed minority carrier diffusion lengths of 16, 10, and 9.2 µm, respectively, implying that the detectors are not limited by minority carrier trapping. The detector performance was observed to be primarily dependent on the concentration and capture cross-sections of the lifetime killing electron traps Z 1/2 and EH 6/7 .

36 MATERIALS SCIENCE↗

Simulating Mixed‐Phase Open Cellular Clouds Observed During COMBLE: Evaluation of Parameterized Turbulence Closure

Marine cold-air outbreaks, or CAOs, are airmass transformations whereby relatively cold boundary layer (BL) air is transported over relatively warm water. To more deeply understand BL and mixed-phase cloud properties during CAO conditions, the Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) took place from late 2019 into early 2020. During COMBLE, the U.S. Department of Energy's first Atmospheric Radiation Measurement Mobile Facility (AMF1) was deployed to Andenes, Norway, far downstream (~1,000 km) from the Arctic pack ice. This study examines the two most intense CAOs sampled at the AMF1 site. The observed BL structures are open cellular with high (~3–5 km) and cold (–30 to –50 °C) cloud tops, and they often have pockets of high liquid water paths (LWPs; up to ~1,000 g m –2 ) associated with strong updrafts and enhanced turbulence. We use a high-resolution mesoscale model to explore how well four turbulence closure methods represent open cellular clouds. After applying a radar simulator to model outputs for direct evaluation, cloud top properties agree well with AMF1 observations (within ~10%), but radar reflectivity and LWP agreement is more variable. Results suggest that the turbulent Prandtl number may play an important role for the simulated BL and cloud properties. All simulations produce enhanced precipitation rates that are well-correlated with a cloud transition. Finally, the eddy-diffusivity/mass-flux approach produces the deepest cloud layer and therefore the largest and most coherent cellular structures. Furthermore, we recommend the use of a non-local turbulence closure approach to better capture turbulent processes in intense CAOs.

54 ENVIRONMENTAL SCIENCES↗

Lithium-ion battery physics and statistics-based state of health model

A pseudo-2d model using COMSOL Multiphysics® software is developed to simulate performance and performance degradation of Li-ion batteries consisting of layered and olivine cathodes with graphite anode when subjected to peak shaving grid service. Multiple degradation pathways are considered, including solid electrolyte interphase (SEI) formation and breakdown at the anode, cathode dissolution and its synergistic effect on SEI formation at the anode. The model is validated by simulating commercial cylindrical cell performance. A global model is developed to simulate performance across all chemistries, along with individual chemistry models using global model parameters as initial values. There is good agreement between these models for various optimization parameters such as SEI equilibrium potential, cathode dissolution exchange current density, solvent diffusivity in the SEI and SEI ionic conductivity. To circumvent time constraints related to the COMSOL model, a 0d global model is developed which fits data well and provides more clarity on differences in cathode dissolution exchange current density. Again, good agreement for various optimization parameters is obtained among the COMSOL global & individual chemistry models and the 0-d model. The lessons learned from the physics-based model is used to develop a top down statistics-based model using current, voltage and anode volumetric change per mole lithium intercalated, along with their interactions as degradation predictors. This model predicts out of sample degradation for multiple grid services and electric vehicle drive cycle with high accuracy and provides the pathway to develop an efficient battery management system combining machine learning and findings from physics-based computationally intensive algorithms.

Crawford, Aladsair J.↗

Mesoscopic analyses of the impact of morphology and operating conditions on the transport resistances in a proton-exchange-membrane fuel-cell catalyst layer

Exploring the origins of local transport resistance and characterizing the oxygen transport resistances in the catalyst layer (RCL) are critical for cost reduction. In this paper, a comprehensive mesoscopic model for simulating coupled transport processes of oxygen and water vapor for different structural parameters under different operating conditions in reconstructed microstructures is proposed. The local transport resistance is calculated after achieving the limiting current density and the effective diffusivity of oxygen. The results demonstrate that RCL increases greatly with decreasing platinum loading (LPt) and the transport resistances in other components of the cell dominate for high-loadings. Both the reduced oxygen permeation coefficient in the ionomer thin-film and the adsorption resistance account for the origins of local transport resistance. The local transport resistance increases with the bare carbon ratio for a constant LPt and Pt/C ratio due to the decreased effective ionomer surface, and increases with the I/C ratio due to the increased ionomer thickness and decreased Knudsen diffusivity. Due to the presence of liquid water, a slight decrease followed by an increase of the local transport resistance versus relative humidity is obtained. The contribution of ionomer thin-films to RCL is more sensitive to liquid saturation compared with that of pores. This journal is

Mu, YT↗

Machine Learning Self-Diffusion Prediction for Lennard-Jones Fluids in Pores

Predicting the diffusion coefficient of fluids under nanoconfinement is important for many applications including the extraction of shale gas from kerogen and product turnover in porous catalysts. Due to the large number of important variables, including pore shape and size, fluid temperature and density, and the fluid–wall interaction strength, simulating diffusion coefficients using molecular dynamics (MD) in a systematic study could prove to be prohibitively expensive. Here, we use machine learning models trained on a subset of MD data to predict the self-diffusion coefficients of Lennard-Jones fluids in pores. Our MD data set contains 2280 simulations of ideal slit pore, cylindrical pore, and hexagonal pore geometries. We use the forward feature selection method to determine the most useful features (i.e., descriptors) for developing an artificial neutral network (ANN) model with an emphasis on easily acquired features. Our model shows good predictive ability with a coefficient of determination (i.e., R 2 ) of ~0.99 and a mean squared error of ~2.9 × 10 –5 . Finally, we propose an alteration to our feature set that will allow the ANN model to be applied to nonideal pore geometries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Prediction of divertor heat flux width for ITER pre-fusion power operation using BOUT++ transport code

Prediction of divertor heat flux width is performed for the first and the second pre-fusion power operation (PFPO) phases specified in the new ITER research plan using BOUT++ transport code (Li et al 2018 Comput. Phys. Commun.228 69–82). Here the initial plasma profiles inside the separatrix are taken from CORSICA scenario studies. Transport coefficients in transport code are calculated by inverting the plasma profiles inside the separatrix and are assumed to be constants in the scrape-off-layer. An anomalous thermal diffusivity scan is performed with E × B and magnetic drifts. The results in two scenarios identifying two distinct regimes: a drift-dominant regime when diffusivity is smaller than the respective critical thermal diffusivity χ c and a turbulence-dominant regime when diffusivity is larger than it. The Goldston heuristic drift model and the ITPA multi-machine experimental scaling yield a lower limit of the width λ q . From transport simulations, we obtain the critical thermal diffusivity χ c = 0.5 m 2 s –1 for the PFPO-1 scenario with toroidal magnetic field B = 1.77 T and plasma current I p = 5 MA, and χ c = 0.3 m 2 s –1 for the PFPO-2 scenario with toroidal magnetic field B = 2.65 T and plasma current I p = 7.5 MA. Separatrix temperature and collisionality also have a significant impact on the heat flux width in the drift-dominant regime. The investigation clearly yields a scaling for critical thermal diffusivity ${\chi }_{\text{c}}\propto {A}^{1/2}/(Z{\left(1+Z\right)}^{\frac{1}{2}}{B}_{\text{p}}^{2})$ using ITER scenarios with fixed safety factor q 95 , major radius R, aspect ratio R/a, and the separatrix temperature T sep , and establishes the connection with CFETR and C-Mod discharges. This scaling implies that for a given tokamak device with q 95 , R, R/a, and T sep fixed, a reduction of poloidal magnetic field by a factor of 3 leads to a 9 times higher critical value of thermal diffusivity χ c , possibly yielding a transition from turbulence- to drift-dominant regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Extraction Kinetics of Rare Earth Elements from Ion-Adsorbed Underclays

Citric acid has been identified as an environmentally sustainable organic acid capable of leaching up to ~30% of easily accessible REEs from underclay material. An analysis of the leaching profiles was performed to discern the reaction rates, extraction efficiencies, and potential leaching mechanisms of REEs and cations of interest from ion-adsorbed underclays. The initial leaching stage follows a slow intraparticle diffusion mechanism followed by a second stage controlled by a mixed diffusion regime. The leaching profiles of Ca and P were similar to those of REEs, suggesting that REEs are most likely derived from mineral surfaces such as hydroxyapatite or crandallite rather than predominately from underclays. Fitting to a modified diffusion control model found diffusion-controlled leaching to be the primary mechanism whereas non-diffusive mechanisms made up about 22% of the extracted REEs. Gangue cations associated with underclays had less non-diffusive leaching than REE species, indicating that their leaching kinetics may be dominated by diffusion from within the material or potentially from product layer formation. Fitting to Boyd plots further indicated that REEs were leached following intraparticle diffusion control. These results have important implications for the development of more efficient and sustainable methods for extracting REEs or critical minerals from alternative feedstocks.

Prem, Priscilla↗

Grazing Incidence Wide-Angle X-ray Scattering of Water Adsorption in Polyamide Barrier Layers of Reverse Osmosis Membranes

To understand the relationship between the intermolecular structure of aromatic polyamide (PA) scaffold and the water molecules in the barrier layers of reverse osmosis (RO) membranes, a grazing incidence wide-angle X-ray scattering (GIWAXS) study was carried out on freestanding PA thin films at varying relative humidity (RH) conditions. The scattering results were analyzed by an interference scattering model, containing a phase factor between a PA chain and an adsorbed water molecule. This model yielded good fits to the GIWAXS profiles where the water adsorption was found to vary linearly with RH. Atomistic molecular dynamics (MD) simulations were also performed to complement the experimental study. Furthermore, the simulations revealed that a rapid condensation layer initially formed on the PA film surface, followed by the slow water molecule diffusion inside the PA membrane. Sparse adsorbed water, isolated in subnanopores of the PA film adjacent to the polar atoms, even in very low quantities, modifies the X-ray scattering. Atomistic simulations at the microscopic scale provide partial support for several X-ray scattering findings.

36 MATERIALS SCIENCE↗

Tools for Assessing Performance (TAP) Project: FY2023 Quarter 2 Report

LANL personnel presented TAP-sponsored work at the 103 rd Annual Meeting of the American Meteorological Society, which was held in Denver, CO in January 2023. There were four different presentations related to TAP-sponsored work that were given by LANL staff and studentsat this meeting: (1) Gloeckler, T., M. Nelson, H. Shah, M. Deshler, L. Wedell, R. Linn, N. Duboc, S. Pol, P. Bieringer, A. Annunzio, B. Martin, and H. Tinnesand, 2023. Testing Validity and Limitations of a Low-Order Diffusive Wake Model. 15 th Symposium on the Urban Environment, January 8-12, 2023, Denver, CO, U.S.A. (2) Nelson, M., H. Shah, T. Gloeckler, M. Deshler, L. Wedell, A. Ortega, R. Linn, N. Duboc, P. Conry, S. Pol, and H. Tinnesand, 2023. The Isolated Building Wake Experiment. Paper 1A.1. 24 th Symposium on Boundary Layers and Turbulence, January 8-12, 2023. Denver, CO, U.S.A. (3) Shah, H., M. Nelson, S. Pol, R. Linn, N. Duboc, and H. Tinnesand, 2023. Detailed Measurement of a Building Wake Under Real-World Atmospheric Conditions. Poster 102. 24 th Symposium on Boundary Layers and Turbulence, January 8-12, 2023, Denver, CO, U.S.A. (4) Wedell, L., M. Nelson, H. Shah, T. Gloeckler, M. Deshler, R. Linn, S. Pol, H. Tinnesand, 2023. The Isolated Building Wake Experiment: Using a Rotating Building to Mulitply the Number of Wake Measurement Locations without Moving the Instruments. Paper 1A.4A. 24 th Symposium on Boundary Layers and Turbulence, January 8-12, 2023, Denver, CO, U.S.A.

54 ENVIRONMENTAL SCIENCES↗