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 667 records · Page 37

Gravitational Influences on Flame Propagation through Non-Uniform, Premixed Gas Systems

Flame propagation through non-uniformly premixed (or layered) gases has importance both in useful combustion systems and in unintentional fires. As summarized previously, non-uniform premixed gas combustion receives scant attention compared to the more usual limiting cases of diffusion or uniformly premixed flames, especially regarding the role gravity plays. This paper summarizes our progress on furthering the knowledge of layered combustion, in which a fuel concentration gradient exists normal to the direction of flame spread. We present experimental and numerical results for flame spread through propanol-air layers formed near the flash point temperature (25 C) or near the stoichiometric temperature (33 C). Both the model and experimental results show that the removal of gravity results in a faster spreading flame, by as much as 80% depending on conditions. This is exactly the opposite effect as that predicted by an earlier model reported. We also found that having a gallery lid results in faster flame spread, an effect more pronounced at normal gravity, demonstrating the importance of enclosure geometry. Also reported here is the beginning of our spectroscopic measurements of fuel vapor.

Miller, Fletcher J.↗

MATBOX, an Open-Source Microstructure Analysis Toolbox for Meshing, Generation, Segmentation, and Characterization of 3D Heterogenous Volumes

Battery performance is strongly correlated with electrode microstructural properties. To account for its impact, lithium-ion battery (LIB) models either abstract the microstructural heterogeneity of composite electrodes using effective macroscopic properties (macro- or meso- scale models) or directly solve the system of equations on the microstructure geometry or mesh (microstructure-scale models). Therefore, to be adequate, both families of models require information from the microstructure geometry, which can be provided by the numerical tool presented in this work. MATBOX is a MATLAB open-source application [1] developed by NREL for performing various microstructure-related tasks including microstructure numerical generation, image filtering and microstructure segmentation, microstructure characterization and correlation, visualization, and microstructure meshing. MATBOX was originally developed for the analysis of LIB electrode microstructures; however, the algorithms provided by the toolbox are widely applicable to other heterogeneous materials. The toolbox provides a user-friendly experience thanks to a Graphical-User Interface, requires no coding by the user, and is well documented. This presentation will illustrate various MATBOX features for the characterization of a LIB electrode, including a fully automated Representative Volume Element (RVE) analysis, the numerical generation of complex 'virtual' microstructure, including dual-layer electrodes and carbon-binder additive phase, and the meshing of a complex NMC/graphite full cell microstructure suitable for 3D finite-element modeling. Other modules (segmentation, visualization, and correlation) will be briefly presented. Thanks to its modular, open-source approach, MATBOX can easily incorporate third-party algorithms to eventually build a standard in the field that will benefit the whole scientific community. Effective diffusion coefficient [2], additive phase numerical generation [3], and meshing [4] third-party algorithms have been already integrated in the toolbox with more to come.

DIRECT ENERGY CONVERSION,MATHEMATICS AND COMPUTING↗

Effect of non-solvent from the phase inversion method on the morphology and performance of the anode supported microtubular solid oxide fuel cells

The microstructure of the anode in anode-supported solid oxide fuel cells has significant influence on the cell performance. Here, microtubular Ni-yttria stabilized zircona (Zr 0.8 Y 0.2 O 2 , YSZ) anode support has been prepared by the phase inversion method. Different compositions of non-solvent have been used for the fabrication of the Ni-YSZ anode support, and the correlation between non-solvent composition and characteristics of the microstructure of the anode support has been investigated. The presence of ethanol or isopropanol in the non-solvent can inhibit the growth of the finger-like pores in the anode support. With the increase of the concentration of ethanol or isopropanol in the non-solvent, a thin dense layer can be observed on the top of the prepared tubular anode support. In addition, the mechanism of pore formation is explained based on the inter-diffusivity between the solvent and the non-solvent. Lastly, the prepared microtubular solid oxide fuel cells (MT-SOFCs) have been tested, and the influence of the anode microstructure on the cell electrochemical performance is analyzed based on a polarization model.

36 MATERIALS SCIENCE↗

Bulky Cation-Modified Interfaces for Thermally Stable Lead Halide Perovskite Solar Cells

Charged conjugated organic molecules offer promising prospects for reducing nonradiative recombination at interfaces in perovskite solar cells, while protecting the active layer from moisture. However, several studies have shown that the heat induced diffusion of these cations leads to irreversible solar cell degradation. Passivation molecules for perovskite can reconstruct the film surface into lower-dimensional phases when exposed to thermal stress, impeding charge extraction and affecting the photoconversion efficiency (PCE) of devices. In this work, we study how molecular interactions between passivation molecules and 3D CsFAPbI 3 perovskite impact stability and charge extraction at the perovskite/hole transport layer interfaces. Two model π- conjugated molecules are studied: 2-([2,2′-bithiophen]-5-yl)ethan-1-aminium iodide (2TI) and 2-(3‴,4′-dimethyl-[2,2′:5′,2″:5″,2‴- quaterthiophen]-5-yl)ethan-1-ammonium iodide (4TmI). We demonstrate that the speed of surface layer reconstruction under thermal stress can be controlled by the cation size and correlate these structural changes with the solar cell performance and stability. Devices treated with 2TI and 4TmI achieve PCEs over 21% and maintain their performance under thermal stress. Our findings demonstrate that thermal stability in PSCs can be achieved via the design engineering of passivation agents, offering a blueprint for developing next-generation passivation molecules.

36 MATERIALS SCIENCE↗

Specific features of fluorescence transfer in multiply scattering randomly inhomogeneous layers under intense laser pumping

Based on the analysis of experimental data on the effect of the pulsed laser pump intensity on the spectral properties and the size of the fluorescent response zone in randomly inhomogeneous fluorescent layers, we found that the amplification of spontaneous and stimulated emission significantly affects the statistical properties of the propagation lengths of the fluorescent field partial components in the layers. The experiments are performed with layers of SiO{sub 2} and TiO{sub 2} nanoparticles saturated with rhodamine 6G, pumped by 532-nm laser radiation in the intensity range corresponding to the transient regime from excitation of spontaneous fluorescence to random lasing in the layer. The experimental data are compared with the results of statistical modelling of fluorescence transfer. It is shown that, even at a pump intensity below the random lasing threshold, the spontaneous emission amplification in a layer leads to a significant increase in the contributions to the fluorescence response from partial components with propagation lengths much larger than the layer thickness. This can be interpreted as a manifestation of the quasi-waveguide effect, in which the probability of propagation of diffuse fluorescence components along the layer over distances many times greater than its thickness and the size of the pumped region increases significantly with a decrease in the characteristic radiation amplification length in the layer. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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

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

Johnston, Christopher O.↗

Enhancing Unknown Waveform Detection by Learning Intra and Inter-domain Dependencies with Advanced Attention Fusion Mechanisms

Detection of unknown waveforms in mission-critical communications is a crucial area of interest for the Department of Energy (DoE). Traditional methods and recent deep learning-based approaches often assume that the training set includes all possible classes, which is impractical for detecting new waveforms. This limitation gives rise to the problem of open-set recognition (OSR), which involves correctly identifying known classes while detecting and rejecting unknown or unseen classes. To address this limitation, we propose a novel dual-domain complex-valued neural architecture that jointly processes time-domain and frequency-domain signal representations using transformer mechanisms. A transformer model is a deep learning architecture that uses self-attention mechanisms to process and learn relationships in sequential data. Our model employs a cosine similarity loss to extract domain-specific features and incorporates a transformer architecture in the latent space to weigh the importance of different features from the time and frequency domains. The transformer layer includes stacked self-attention and cross-attention modules to learn intra-domain and inter-domain dependencies, creating a more holistic signal representation. An attention-based fusion module intelligently combines the time and frequency-domain features using multi-head attention, enabling the network to learn the optimal feature for each domain in each input signal. Quantitative results demonstrate the impact of these architectural choices on overall performance, showing significant improvement after incorporating self and cross-attention modules and using complex attention fusion over simple weighted fusion. Our ongoing work will focus on addressing the limitations of threshold-based OSR methods by developing a novel generative framework that integrates a conditional diffusion probabilistic model (DPM). DPM is a generative framework that learns to synthesize complex data by reversing a gradual noising process using a neural network trained to denoise step-by-step. Our goal is to leverage the inherent strengths of DPMs for identifying unknown signals more robustly. One primary advantage of using a DPM is its ability to provide a more reliable anomaly score based on the model's reconstruction error, rather than relying solely on classifier confidence. Additionally, the iterative denoising process of DPMs makes this approach naturally resilient to low Signal-to-Noise Ratio (SNR) conditions, where traditional methods often fail. By implementing this generative framework, we aim to enhance the model's capability to accurately detect unknown waveforms and maintain performance in challenging environments.

99 - GENERAL AND MISCELLANEOUS↗

Estimating viscoelastic compliance of desiccating cementitious materials using drying prism tests

Highlights: • Determining relative humidity and degree of reaction of desiccating cementitious materials via an integrated method • Extracting moisture diffusion coefficient of cement mortar using the mass loss data of the prism specimens • Estimating uniaxial viscoelastic compliance of cement mortar using the shrinkage strain data of the prism specimens An integrated modeling approach to determining internal relative humidity and degree of reaction of desiccating cementitious materials has been developed. By simultaneously accounting for both self-desiccation and external drying, the model is capable of predicting the evolution of material properties of desiccating cementitious materials. Using the mass loss data of the prism specimens, the moisture diffusion coefficient of a nonlinear model for cement mortar was determined. Additionally, a previous poroviscoelastic shrinkage model was modified and extended. Through comparison of predicted and measured axial strain of the prism specimens, the viscoelastic compliance of cement mortar was extracted. The present work revealed that it is possible to use a simple drying prism test to fit the moisture diffusion coefficient and creep compliance of cementitious materials, though the comparison of the extracted viscoelastic compliance with measured compliance was only reasonable if the outer layer of the drying mortar was presumed to have microcracking.

36 MATERIALS SCIENCE↗

Systematic Evaluation of Counterpoise Correction in Density Functional Theory

A widespread belief persists that the Boys–Bernardi function counterpoise (CP) procedure “overcorrects” supramolecular interaction energies for the effects of basis-set superposition error. To the extent that this is true for correlated wave function methods, it is usually an artifact of low-quality basis sets. The question has not been considered systematically in the context of density functional theory, however, where basis-set convergence is generally less problematic. We present a systematic assessment of the CP procedure for a representative set of functionals and basis sets, considering both benchmark data sets of small dimers and larger supramolecular complexes. The latter include layered composite polymers with ~150 atoms and ligand–protein models with ~300 atoms. Provided that CP correction is used, we find that intermolecular interaction energies of nearly complete-basis quality can be obtained using only double-ζ basis sets. Furthermore, this is less expensive as compared to triple-ζ basis sets without CP correction. CP-corrected interaction energies are less sensitive to the presence of diffuse basis functions as compared to uncorrected energies, which is important because diffuse functions are expensive and often numerically problematic for large systems. Our results upend the conventional wisdom that CP “overcorrects” for basis-set incompleteness. In small basis sets, CP correction is mandatory in order to demonstrate that the results do not rest on error cancellation.

74 ATOMIC AND MOLECULAR PHYSICS↗

Cooperative mechanisms behind nanoscale smoothing on metal surfaces: From adatom diffusion to step nucleation

Many technologies require stable or metastable surface morphology. Here, in this paper we study the factors that control the metastability of a common feature of rough surfaces: “hillocks.” We use low energy electron microscopy to follow the evolution of the individual atomic steps in hillocks on Pd(111). We show that the uppermost island in the stack often adopts a static, metastable configuration. Modeling this result shows that the degree of the metastability depends on the configuration of steps dozens of atomic layers lower. Our model allows us to link surface metastability to the atomic processes of surface evolution.

Fernandes Cauduro, Andre Luis↗

Transit time and charge storage measurements in heavily doped emitters

A first direct measurement of the minority-carrier transit time in a transparent heavily doped emitter layer is reported. The value was obtained by a high-frequency conductance method recently developed and used for low-doped Si. The transit time coupled with the steady-state current enables the determination of the quasi-static charge stored in the emitter and the quasi-static emitter capacitance. Using a transport model, from the measured transit time, the value for the minority-carrier diffusion coefficient and mobility is estimated. The measurements were done using a heavily doped emitter of the Si p(+)-n-p bipolar transistor. The new result indicates that the position-averaged minority-carrier diffusion coefficients may be much smaller than the corresponding majority-carrier values for emitters having a concentration ranging from about 3 x 10 to the 19th per cu cm to 10 to the 20th per cu cm.

Neugroschel, A.↗

Polar and Non-Polar Layers on Mars: A Single Mechanism for Formation?

The recent discovery of vast quantities of near-subsurface ice in both polar regions of Mars by the Mars Odyssey Gamma Ray Spectrometer (GRS) has presented us with an interesting quandary. On one hand, these deposits, found poleward of 60 deg in both hemispheres, are consistent with thermal models suggesting ice will be best protected in these regions during periods of high obliquity. On the other hand, the current paradigm regarding the placement of these deposits, i.e., diffusive deposition of water vapor, appears to be inconsistent with the large volume mixing ratios (approx. 90%) inferred from the GRS data. This incongruity argues that diffusion alone cannot be the primary mechanism for the creation of these reservoirs, and that an alternate, large-scale process should be considered.

Mischna, M. A.↗

Uncovering grain and subgrain microstructure at the scale of additive manufacturing melt tracks with a scalable cellular automaton solidification model

Metal additive manufacturing, characterized by rapid solidification, yields refined grains with a distinctive cellular subgrain microstructure that plays a pivotal role in determining material properties. Due to the significant computational expense demanded to simulate the required physics with submicron spatial resolution, their numerical simulations have been limited to proof-of-concept studies to either 2D or small subregions of a melt pool. In this study, an open-source, scalable, solidification code, muMatScale, based on the cellular automaton method, has been developed to predict the grain and the underlying subgrain microstructure over an entire melt pool. The model incorporates flexible parallelization schemes, utilizing MPI and OpenMP GPU Offloading, in addition to appropriate multi-physics specific to non-equilibrium rapid solidification in AM. The impact of nucleation parameters on grain microstructures was investigated with a focus on grain size variations and morphology transitions. With selected nucleation parameters, the simulation predicted the grain size, subgrain morphology, crystallographic orientation, and microsegregation aligned with experimental measurements. The model demonstrates that epitaxial grain growth is a dominant factor at the melt pool boundary, influencing grain size variation under different grain sizes in the build plate while maintaining consistent primary dendrite arm spacing under identical thermal conditions. Here, the highly efficient numerical model enables large-scale simulations with a spatial resolution of 100 nm or less, unveiling unprecedented insights into thermal and solutal diffusion driven grain growth, and the subgrains with microsegregation within grains in 3D across scales. muMatScale will enable the linking of submicron length-scale microstructure to part-level material behavior by investigating fundamental solidification problems at the intercellular scale in many-track and many-layer builds.

36 MATERIALS SCIENCE↗

Diffusion of Sn donors in β-Ga 2 O 3

Diffusion of the n-type dopant Sn in β-Ga 2 O 3 is studied using secondary-ion mass spectrometry combined with hybrid functional calculations. The diffusion of Sn from a Sn-doped bulk substrate with surface orientation (001) into an epitaxial layer is observed after heat treatments in the temperature range of 1050–1250 °C. Calculated formation energies of Sn-related and intrinsic defects show that the migration of Sn is mediated by Ga vacancies ( V Ga ) through the formation and dissociation of intermittent mobile V Ga Sn Ga complexes. The evolution of the Sn concentration vs depth profiles after heat treatments can be well described by a reaction–diffusion model. Using model parameters guided by the hybrid functional calculations, we extract a V Ga Sn Ga complex migration barrier of 3.0 ± 0.4 eV with a diffusion coefficient of 2 × 10 –1 cm 2 /s. The extracted migration barrier is consistent with our theoretical predictions using the nudged elastic band method, which shows migration barriers of 3.42, 3.15, and 3.37 eV for the [100], [010], and [001] directions, respectively.

36 MATERIALS SCIENCE↗

Two-dimensional simulations of internal gravity waves in a $\mathrm{5M_⊙}$ zero-age-main-sequence model

Main-sequence intermediate-mass stars present a radiative envelope that supports internal gravity waves (IGWs). Excited at the boundary with the convective core, IGWs propagate towards the stellar surface and are suspected to impact physical processes such as rotation and chemical mixing. Using the fully compressible time-implicit code MUSIC, we study IGWs in two-dimensional simulations of a zero-age-main-sequence 5 solar mass star model up to 91 percent of the stellar radius with different luminosity and radiative diffusivity enhancements. Our results show that low-frequency waves excited by core convection are strongly impacted by radiative effects as they propagate. This impact depends on the radial profile of radiative diffusivity which increases by almost 5 orders of magnitude between the centre of the star and the top of the simulation domain. In the upper layers of the simulation domain, we observe an increase of the temperature. Our study suggests that this is due to heat added in these layers by IGWs damped by radiative diffusion. We show that non-linear effects linked to large amplitude IGWs may be relevant just above the convective core. Both these effects are intensified by the artificial enhancement of the luminosity and radiative diffusivity, with enhancement factors up to 10 4 times the realistic values. Our results also highlight that direct comparison between numerical simulations with enhanced luminosity and observations must be made with caution. Finally, our work suggests that thermal effects linked to the damping of IGWs could have a non-negligible impact on stellar structure.

79 ASTRONOMY AND ASTROPHYSICS↗

Formation of polar stratospheric clouds simulated in a two dimensional model of the atmosphere

A microphysics code has been implemented in a two dimensional model of the atmosphere to study formation of polar stratospheric clouds containing HCl or HNO3. The model range is from pole to pole in latitude and from the ground to about 20 km in altitude. Resolution in latitude is 10 deg and about 0.8 km in altitude. This is an Eulerian model with prescribed eddy diffusion coefficients and the circulation obtained from observations. The chemistry of the model follows the family approach for NO(x), Cl(x) and HO(x) while the ozone is fixed and changed seasonally. The aerosol code is based on an assigned population of condensation, coagulation and sedimentation. Aerosol growth is simulated in nine different size bins ranging between 0.01u and 2.56u. The model has been built to study aerosol layers formation in the upper troposphere and lower stratosphere and has been validated for sulfate aerosol resulting from a rather complex sulfur chemistry.

Visconti, Guido↗

Design, analysis, and fabrication of oxide-coated iridium/rhenium combustion chambers

Iridium-coated rhenium (Ir/Re) combustion chambers provide high temperature, oxidation-resistant operation for radiation-cooled liquid-fueled rocket engines. A 22-N (5-lb(sub f)) chamber has been operated for 15 hours at 2200 C (4000 F) using nitrogen tetroxide/monomethyl hydrazine (NTO/MMH) propellant, with negligible internal erosion. The oxidation resistance of these chambers could be further increased by the addition of refractory oxide coatings, providing longer life and/or operation in more oxidizing and higher temperature environments. The oxide coatings would serve as a thermal and diffusion barrier for the iridium coating, lowering the temperature of the iridium layer while also preventing the ingress of oxygen and egress of iridium oxides. This would serve to slow the failure mechanisms of Ir/Re chambers, namely the diffusion of rhenium to the inner surface and the oxidation of iridium. Such protection could extend chamber lifetimes by tens or perhaps hundreds of hours, and allow chamber operation on stoichiometric or higher mixture ratio oxygen/hydrogen (O2/H2) propellant. Extensive thermomechanical, thermochemical, and mass transport modeling was performed as a key material/structure design tool. Based on the results of these analyses, several 22-N oxide-coated Ir/Re chambers were fabricated and delivered to NASA Lewis Research Center for hot-fire testing.

Jang, Q.↗

H2O-H2SO4 system in Venus' clouds and OCS, CO, and H2SO4 profiles in Venus' troposphere

A coupled problem of diffusion and condensation is solved for the H2SO4-H2O system in Venus' cloud layer. The position of the lower cloud boundary and profiles of the H2O and H2SO4 vapor mixing ratios and of the H2O/H2SO4 ratio of sulfuric acid aerosol and its flux are calculated as functions of the column photochemical production rate of sulfuric acid, Phi(sub H2SO4). Variations of the lower cloud boundary are considered. Our basic model, which is constrained to yield f(sub H2O) (30 km) = 30 ppm (Pollack et al. 1993), predicts the position of the lower cloud boundary at 48.4 km coinciding with the mean Pioneer Venus value, the peak H2SO4 mixing ratio of 5.4 ppm, and the H2SO4 production rate Phi(sub H2SO4) = 2.2 x 10(exp 12)/sq cm/s. The sulfur to sulfuric acid mass flux ratio in the clouds is 1:27 in this model, and the mass loading ratio may be larger than this value if sulfur particles are smaller than those of sulfuric acid. The model suggests that the extinction coefficient of sulfuric acid particles with radius 3.7 micrometers (mode 3) is equal to 0.3/km in the middle cloud layer. The downward flux of CO is equal to 1.7 x 10(exp 12)/sq cm/s in this model. Our second model, which is constrained to yield f(sub H2O) = 10 ppm at the lower cloud boundary, close to the value measured by the Magellan radiooccultations, predicts the position of this boundary to be at 46.5 km, which agrees with the Magellan data; f(sub H2O) (30 km) = 90 ppm, close to the data of Moroz et al. (1983) at this altitude; Phi(sub H2SO4) = 6.4 x 10(exp 12)/sq cm/s; and Phi(sub co) = 4.2 x 10(exp 12)/sq cm/s. The S/H2SO4 flux mass ratio is 1:18, and the extinction coefficient of the mode 3 sulfuric acid particles is equal to 0.9 km in the middle cloud layer. A strong gradient of the H2SO4 vapor mixing ratio near the bottom of the cloud layer drives a large upward flux of H2SO4, which condenses and forms the excessive downward flux of liquid sulfuric acid, which is larger by a factor of 4-7 than the flux in the middle cloud layer. This is the mechanism of formation of the lower cloud layer. Variations of the lower cloud layer are discussed. Our modeling of the OCS and CO profiles in the lower atmosphere measured by Pollack et al. (1993) provides a reasonable explanation of these data and shows that the rate coefficient of the reaction SO3 + CO yields CO2 + SO2 is equal to 10(exp -11) exp(-(13,100 +/- 1000)/T) cc/s. The main channel of the reaction between SO3 and OCS is CO2 + (SO)2, and its rate coefficient is equal to 10(exp -11) exp(-(8900 +/- 500)/T) cc/s. In the conditions of Venus' lower atmosphere, (SO)2 is removed by the reaction (SO)2 + OCS yields CO + S2 + SO2. The model predicts an OCS mixing ratio of 28 ppm near the surface.

Krasnopolsky, V. A.↗