SEARCH · Engineering Papers

Results for “mass diffusivity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Thermoacoustic Soret separation

In a fluid mixture in a channel with an axial time-averaged temperature gradient, high-amplitude oscillating flow can greatly increase the axial flux of thermal diffusion (Soret) separation of the components of the mixture. The enhancement occurs when the oscillating lateral temperature gradient greatly exceeds the axial gradient, causing a large oscillating concentration that can be favorably time-phased with the oscillating flow. This process can occur even with a negligible pressure oscillation or with a negligible temperature response to pressure, as is the case in most liquid solutions. The thermal boundary condition imposed by realistic solids on thermoacoustic liquids is imperfect, adding mathematical complications that are absent for typical gases, for which the solid surface is temporally isothermal. Compared with gas mixtures, the high Lewis number in typical liquid solutions reduces the separation flux associated with the time-averaged temperature gradient, but it also reduces the remixing associated with the time-averaged mole-fraction gradient. For large enough channels, the second-law separation efficiency is only slightly reduced from that of steady liquid Soret separation.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Pit rim decomposition into multiple quantum dots on surfaces of epitaxial thin films grown on pit-patterned substrates

Here, we report results of dynamical simulations according to an experimentally validated surface morphological evolution model on the formation of regular arrays of quantum dot molecules (QDMs) consisting of 1D arrays of smaller interacting quantum dots (QDs). These QD arrays form along the sides of each pit rim on the surface of a coherently strained thin film epitaxially deposited on a semiconductor substrate, the surface of which consists of a periodic pattern of inverted prismatic pits with rectangular pit openings. We find that this complex QDM pattern results from the decomposition of the pit rim from a “quantum fortress” with four elongated QDs into four 1D arrays of multiple smaller QDs arranged along each side of the pit rim. Systematic parametric analysis indicates that varying the pit opening dimensions and the pit wall inclination directly impacts the number of QDs in the resulting QDM pattern, while varying the pit depth only affects the dimensions of the QDs in the nanostructure pattern. Therefore, the number, arrangement, and sizes of QDs in the resulting pattern of QDMs on the film surface can be engineered precisely by proper tuning of the pit design parameters. Our simulation results are supported by predictions of morphological stability analysis, which explains the pit rim decomposition into multiple QDs as the outcome of a tip-splitting instability and provides a fundamental characterization of the post-instability nanostructure pattern. Our theoretical findings can play a vital role in designing optimal semiconductor surface patterns toward enabling future nanofabrication technologies.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Fluid modeling of low-temperature plasmas

Fluid models are essential for understanding and predicting low-temperature plasma (LTP) behavior in various scientific and industrial settings. This paper provides an introductory tutorial on fluid modeling of LTPs, covering model formulation, implementation, and computational simulations. The tutorial focuses on five main components of the formulation of LTP fluid models: fluid flow, energy, chemistry, electromagnetism, and material properties, as well as in essential aspects of model implementations, including multiscale phenomena, multiphysics coupling, and numerical convergence. Designed for students and early-career researchers, this work offers a practical foundation for developing and using fluid models, from in-house computational codes to commercial software, bridging fundamental theory with real-world applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Decoherence and Brownian motion of a polarizable particle near a medium

Optically levitated nanoparticles are ideal experimental testbeds for investigating macroscopic superpositions and microscopic thermodynamics. Integrating such levitated nanoparticles with photonic structures can enable strong coupling between their center-of-mass motion and guided photonic modes, facilitating enhanced control and probing of their motion. When coupling a particle to a photonic structure, such as a waveguide, the effects of fluctuations become prominent at nanoscales. In this work, we analyze the classical and quantized center-of-mass motion of a polarizable particle interacting with the fluctuations of the electromagnetic field in the presence of a medium. We derive a position localization master equation for the particle's quantized center of mass, and examine its classical center-of-mass momentum diffusion, elucidating correspondences between classical and quantum Brownian motion of polarizable particles near media. We study the decoherence rate of the particle in the presence of a planar surface as a function of temperature and distance from the surface, comparing it to common sources of decoherence. Our results are pertinent to experiments aimed at preparing levitated nanospheres in macroscopic quantum states and investigating their Brownian dynamics.

Brownian motion↗

DOE OSTI2022

The Influence of Shallow Cloud Populations on Transitions to Deep Convection in the Amazon

Abstract In this study, a pair of convection-permitting (2-km grid spacing), month-long, wet-season Weather Research and Forecasting (WRF) Model simulations with and without the eddy-diffusivity mass-flux (EDMF) scheme are performed for a portion of the Green Ocean Amazon (GoAmazon) 2014/15 field campaign period. EDMF produces an ensemble of subgrid-scale convective plumes that evolve in response to the boundary layer meteorological conditions and can develop into shallow clouds. The objective of this study is to determine how different treatments of shallow cumulus clouds (i.e., with and without EDMF) impact the total cloud population and precipitation across the Amazonian rain forest, with emphasis on impacts on the likelihood of shallow-to-deep convection transitions. Results indicate that the large-scale synoptic conditions in the EDMF and control simulations are nearly identical; however, on the local scale their rainfall patterns diverge drastically and the biases decrease in EDMF. The EDMF scheme significantly increases the frequency of shallow clouds, but the frequencies of deep clouds are similar between the simulations. Deep convective clouds are tracked using a cloud-tracking algorithm to examine the impact of shallow cumulus on the surrounding ambient environment where deep convective clouds initiate. Results suggest that a rapid increase of low-level cloudiness acts to cool and moisten the low to midtroposphere during the day, favoring the transition to deep convection.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Diagnosing Near-Surface Model Errors with Candidate Physics Parameterization Schemes for the Multiphysics Rapid Refresh Forecast System (RRFS) Ensemble during Winter over the Northeastern United States and Southern Great Plains

Abstract During the winter of 2020/21 an ensemble of FV3-LAM forecasts was produced over the contiguous United States for the Winter Weather Experiment using five physics suites. These forecasts are evaluated with the goal of optimizing physics parameterizations within the future operational Rapid Refresh Forecast System (RRFS) in the Unified Forecast System (UFS) realm and for selecting suitable physics suites for a multiphysics RRFS ensemble. The five physics suites have different combinations of land surface models (LSMs), planetary boundary layer (PBL) parameterizations, and surface layer schemes, chosen from those used in current and possible future operational systems and likely to be supported in the operational UFS. Full-season evaluation reveals a persistent near-surface cold bias in the U.S. Northeast from one suite and a nighttime warm bias in the southern Great Plains in another suite, while other suites have smaller biases. A representative case is chosen to diagnose the cause for each of these biases using sensitivity simulations with different physics combinations or modified parameters and verified with additional mesonet observations. The cold bias in the Northeast is attributed to aspects of the Noah-MP LSM over snow cover, where Noah-MP simulates lower soil water content, and thus lower thermal conductivity than other LSMs, leading to less upward ground heat flux during nighttime and consequently lower surface temperature. The nighttime warm bias found in the southern Great Plains is attributed to overestimation of vertical mixing in the K -profile-based eddy-diffusivity mass-flux (K-EDMF) PBL scheme and insufficient land–atmospheric coupling from the GFS surface layer scheme over short vegetation. A few key parameters driving these systematic biases are identified.

Meteorology & Atmospheric Sciences↗

DOE OSTI2024

Large eddy simulation of a supersonic lifted hydrogen flame: Impacts of Lewis, turbulent Schmidt and Prandtl numbers

Parametric large eddy simulations (LES) of a supersonic lifted hydrogen flame are reported. The emphases are on two aspects: impacts of (1) Lewis number (Lei of the ith species) and (2) turbulent Schmidt and Prandtl numbers (⁠Sc t and Pr t ⁠) on supersonic turbulent flame and flow structures. Five cases are considered: species-specific Le i , Sc t = Sr t = 1.0 (C 0 ); unity Le i , Sc t = Pr t = 1.0 (C 1 ); species-specific Sc t = 0.5, Pr t = 1.0 (C 2 ); species-specific Le i , Sc t = 1.0, Pr t = 0.5 (C 3 ); and species-specific Le i , Sc t = Pr t = 0.5 (C 4 ). Numerical results of instantaneous and/or time-averaged species mole fractions, mixture fraction, heat release rate, flame base location, and mixed modes of premixed and diffusion combustion are compared between cases C 0 and C 1 . Differences in auto-ignition locations and strengths and flame structures and stabilization specify the impacts of Lewis number. They are triggered by different predictions of species mass and thermal diffusions at fuel-coflow and/or coflow-ambient air mixing layers. These differences are rationalized by a scale analysis of mass/thermal diffusion and convection for case C 0 , which suggests the relatively low but non-negligible former against the latter. Cases C 0 and C 2 –C 4 barely see differences in terms of instantaneous and/or time-averaged temperature, velocity, and mixed combustion modes except for further downstream areas where combustion occurs. Both Sc t and Pr t impose less significant influences than Lewis number, as sub-grid scale mass/thermal diffusion is subordinate to its resolved counterpart according to their scale analysis for case C 4 .

42 ENGINEERING↗

DOE OSTI2022

Influence of the Soret effect on binary-species boundary layers at high pressure

Direct numerical simulations of binary-species temporal boundary layers at high pressure are performed. The main objective is to investigate the influence of the Soret effect on flow physics of binary-species boundary layers where the fluid has a uniform composition. The working fluid is a mixture of 25 % methane and 75 % nitrogen in mass fraction. Although the fluid composition is uniform at the initial condition, the mass fraction of methane increases near the wall when the wall temperature is hotter than the free stream temperature, whereas it decreases when the wall temperature is colder. The non-uniform mass fraction indicates that the uphill diffusion occurs near the wall. Investigation of fluctuations of the mass fraction reveals that the mass fraction fluctuates in the whole boundary layer, indicating that the uphill diffusion occurs even far from the wall. Examination of the species-mass diffusion balance for mean flow fields clarified that the Soret effect flux becomes large near the wall, and the large flux causes the non-uniform profile of the mass fraction near the wall.

Binary Diffusion Coefficient↗

DOE OSTI2023

Mass transfer in catalytic depolymerization: External effectiveness factors and serendipitous processivity in stagnant and stirred melts

Several heterogeneous catalysts are being developed to recycle plastics. Most operate in viscous polymer melts, where external mass transfer effects could limit the supply of co-reactants to active sites. External mass transfer can also impede the diffusion of long chain products away from the catalyst after each cut. Product egress limitations could potentially confer unintentional processivity to catalyst operation, i.e. a tendency for the catalyst to repeatedly cut the same chain after an initial encounter. We formulate reaction–diffusion equations to quantify mass transfer effects on the co-reactant transport to the catalyst and the degree of serendipitous processivity. Results are developed for catalysts in stagnant or stirred melts, with simple expressions involving Damkohler, Peclet, and Sherwood numbers, i.e. dimensionless combinations of rate constants, catalyst particle size, polymer diffusivities, and shear rates (where applicable). In conclusion, we estimate the impact of these effects for a spherical core–shell catalyst.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Heavy Quark Diffusion from 2+1 Flavor Lattice QCD with 320 MeV Pion Mass

We present the first calculations of the heavy flavor diffusion coefficient using lattice QCD with light dynamical quarks corresponding to a pion mass of around 320 MeV. For temperatures 195 MeV< 𝑇 <352 MeV, the heavy quark spatial diffusion coefficient is found to be significantly smaller than previous quenched lattice QCD and recent phenomenological estimates. The result implies very fast hydrodynamization of heavy quarks in the quark-gluon plasma created during ultrarelativistic heavy-ion collision experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Pore connectivity influences mass transport in natural rocks: Pore structure, gas diffusion and batch sorption studies

For this work, six rocks (one granodiorite, one limestone, two chalks, one mudstone, and one dolostone) with different extents of heterogeneity at six different particle sizes (from 75 to 8000 μm) were studied to describe the effects of pore connectivity on mass transport. The methods applied were (i) porosity measurement of granular rocks, (ii) analyses of gas-phase diffusive transport in a bed of packed particles, along with a solid quartz method at these six particle sizes being developed to identify the contribution of intraparticle diffusion, and (iii) batch sorption tests of multiple ions (anions and cations) with subsequent analyses of inductively coupled plasma-mass spectrometry. Granular porosity measurement results reveal that with decreasing particle sizes, the effective porosities for the “heterogenous” group of rocks (Grimsel granodiorite and Edwards limestone) increase, whereas the porosities of another “homogeneous” group (two Israel chalk samples, Japan mudstone, and Wyoming dolostone) remain constant. Gas diffusion results show that the intraparticle gas diffusion coefficient among these two sample groups, varying in the magnitude of 10 -8 to 10 -6 m 2 /s, are not directly correlated to the porosity differences. Moreover, the batch sorption work displays a different affinity of rocks for various tracers. For Grimsel granodiorite, Japan mudstone, and Wyoming dolostone, the adsorption capacity of Sm 3+ and Eu 3+ increases as the particle size decreases. In general, this integrated research of grain size distribution, granular rock porosity, intraparticle diffusivity, and ionic sorption capacity gives insights into the pore connectivity effect on both physical and chemical transport behaviors for different lithologies and/or different particle sizes.

58 GEOSCIENCES↗

DOE OSTI2023

Thermokinetics of point defects in α-Fe 2 O 3

Abstract Point defect formation and migration in oxides governs a wide range of phenomena from corrosion kinetics and radiation damage evolution to electronic properties. In this study, we examine the thermodynamics and kinetics of anion and cation point defects using density functional theory in hematite ( α -Fe 2 O 3 ), an important iron oxide highly relevant in both corrosion of steels and water-splitting applications. These calculations indicate that the migration barriers for point defects can vary significantly with charge state, particularly for cation interstitials. Additionally, we find multiple possible migration pathways for many of the point defects in this material, related to the low symmetry of the corundum crystal structure. The possible percolation paths are examined, using the barriers to determine the magnitude and anisotropy of long-range diffusion. Our findings suggest highly anisotropic mass transport in hematite, favoring diffusion along the c -axis of the crystal. In addition, we have considered the point defect formation energetics using the largest Fe 2 O 3 supercell reported to date.

36 MATERIALS SCIENCE↗

DOE OSTI2022

On the Spectral Evolution of Hot White Dwarf Stars. II. Time-dependent Simulations of Element Transport in Evolving White Dwarfs with STELUM

White dwarf stars are subject to various element transport mechanisms that can cause their surface composition to change radically as they cool, a phenomenon known as spectral evolution. In this paper, we undertake a comprehensive theoretical investigation of the spectral evolution of white dwarfs. First, we introduce STELUM, a new implementation of the stellar evolutionary code developed at the Université de Montréal. We provide a thorough description of the physical content and numerical techniques of the code, covering the treatment of both stellar evolution and chemical transport. Then, we present two state-of-the-art numerical simulations of element transport in evolving white dwarfs. Atomic diffusion, convective mixing, and mass loss are considered simultaneously as time-dependent diffusive processes and are fully coupled to the cooling. We first model the PG 1159-DO-DB-DQ evolutionary channel: a helium-, carbon-, and oxygen-rich PG 1159 star transforms into a pure-helium DB white dwarf due to gravitational settling and then into a helium-dominated, carbon-polluted DQ white dwarf through convective dredge-up. We also compute for the first time the full DO-DA-DC evolutionary channel: a helium-rich DO white dwarf harboring residual hydrogen becomes a pure-hydrogen DA star through the float-up process and then a helium-dominated, hydrogen-bearing DC star due to convective mixing. We demonstrate that our results are in excellent agreement with available empirical constraints. In particular, our DO-DA-DC simulation perfectly reproduces the lower branch of the bifurcation observed in the Gaia color–magnitude diagram, which can therefore be interpreted as a signature of spectral evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Sublimation and oxidation measurements of graphite and carbon black at high temperatures in a shock tube using absorption imaging and thermal emission

Surface mass loss rates due to sublimation and oxidation at temperatures of 3000–7000 K have been measured in a shock tube for graphite and carbon black (CB) particles. Diagnostics are presented for measuring surface mass loss rates by diffuse backlit illumination extinction imaging and thermal emission. Here, the surface mass loss rate is found by regression fitting extinction and emission signals with an independent spherical primary particle assumption. Measured graphite sublimation and oxidation rates are reported to be an order of magnitude greater than CB sublimation and oxidation rates. It is speculated that the difference between CB and graphite surface mass loss rates is largely due to the primary particle assumption of the presented technique which misrepresents the effective surface area of an aggregate particle where primary particles overlap and shield inner particles. Measured sublimation rates are compared to sublimation models in the literature, and it is seen graphite shows fair agreement with the models while CB underestimates, likely a result of the particle shielding affect not being considered in the sublimation model.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Simulations of flaming combustion and flaming-to-smoldering transition in wildland fire spread at flame scale

Our objective in the present study is to provide basic insights into the coupling between external-gas and solid biomass vegetation processes that control the dynamics of flame spread in wildland fire problems. Here we focus on a modeling approach that resolves processes occurring at vegetation and flame scales, i.e., the formation of flammable vapors due to the thermal degradation of the solid biomass, the subsequent combustion in ambient air, the thermal feedback to the biomass through radiative and convective heat transfer, and the possible transition from flaming combustion (taking place outside of the solid biomass) to smoldering combustion (taking place inside the solid biomass). The capability uses a multiphase combustion framework and treats external-gas processes through a Large Eddy Simulation solver and solid biomass processes through a discrete particle model. The discrete particle model adopts a one-dimensional porous medium formulation, includes descriptions of drying, thermal pyrolysis, oxidative pyrolysis, and char oxidation, as well as a description of the external-gas-to-solid-biomass diffusion of oxygen mass; the discrete particle model thereby provides a treatment of in-depth oxidative processes and allows the simulation of smoldering combustion. The modeling capability is applied to the simulation of fire spread across a surrogate biomass vegetation bed corresponding to a discrete array of cylindrical-shaped, vertically-oriented, pine wood sticks, characterized by a monomodal size distribution, in horizontal flat terrain and under wind-aided conditions. The numerical results demonstrate that the model can simulate successful flaming-to-smoldering transition followed by complete biomass consumption.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2026

A multiscale packed-bed reactor model for sustainable ethylene production via chemical looping oxidative coupling of methane

The rising global warming concerns and shale gas discovery have prompted research in the direction of greenhouse gas (GHG), such as methane, reduction and conversion. Oxidative coupling of methane (OCM) offers a pathway to low carbon-intense valorization of methane while producing ethylene, a chemical regarded as central to the petrochemical industry. Even after decades of OCM discovery, researchers keep understanding the process and underlying chemical reactions in a pursuit to achieve industrial viability for OCM. Here, in general, OCM suffers from low C 2 selectivity, yield and reactor temperature runaways due to highly exothermic nature of its reactions. Computational Fluid Dynamics (CFD) tools help analyze spatial gradients within the reactor to deeply understand the diffusion of species, mass and heat transfer phenomena. Furthermore, challenges associated with scaling up such as hot spot formation and parametric sensitivity can be addressed without having to expend on costly experiments. The current paper presents a multiscale packed-bed reactor CFD model coupled with a chemical kinetic model for the chemical looping OCM. The CFD model includes two scales i.e., macroscale for catalyst bed and microscale for individual pellets. Moreover, a chemical kinetic model based on 10 gas-phase reactions is integrated with the CFD model. An additional surface reaction for the formation of gas-phase oxygen from catalyst surface is added to account for the absence of feed oxygen. The model is calibrated against experimental results. The calibrated model captures trends in CH 4 conversion, C 2 selectivity and C 2 yield within a ± 4.35 % range across a temperature range of 700-900 °C. Moreover, model fidelity is evaluated by varying key computational parameters such as mesh resolution and time step size. The model is also verified by varying the inlet methane concentration and the gas hourly space velocity (GHSV) and comparing the results with literature. A sensitivity analysis and scale-up of the current model is undergoing.

Chemical looping↗

DOE OSTI2024

Sequential Monte Carlo for Cut-Bayesian Posterior Computation

We propose a sequential Monte Carlo (SMC) method to efficiently and accurately compute cut-Bayesian posterior quantities of interest, variations of standard Bayesian approaches constructed primarily to account for model misspecification. We prove finite sample concentration bounds for estimators derived from the proposed method along with a linear tempering extension and apply these results to a realistic setting where a computer model is misspecified. We then illustrate the SMC method for inference in a modular chemical reactor example that includes submodels for reaction kinetics, turbulence, mass transfer, and diffusion. The samples obtained are commensurate with a direct-sampling approach that consists of running multiple Markov chains, with computational efficiency gains using the SMC method. Overall, the SMC method presented yields a novel, rigorous approach to computing with cut-Bayesian posterior distributions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Efficient Reformulation of Linear and Nonlinear Solid-Phase Diffusion in Lithium-ion Battery Models using Symmetric Polynomials: Mass Conservation and Computational Efficiency

Lithium-ion batteries are typically modeled using porous electrode theory coupled with various transport and reaction mechanisms, along with suitable discretization or approximations for the solid-phase diffusion equation. The solid-phase diffusion equation represents the main computational burden for typical pseudo-2-dimensional (p2D) models since these equations in the pseudo r -dimension must be solved at each point in the computational grid. This substantially increases the complexity of the model as well as the computational time. Traditional approaches towards simplifying solid-phase diffusion possess certain significant limitations, especially in modeling emerging electrode materials which involve phase changes and variable diffusivities. A computationally efficient representation for solid-phase diffusion is discussed in this paper based on symmetric polynomials using Orthogonal Collocation and Galerkin formulation (weak form). A systematic approach is provided to increase the accuracy of the approximation (p form in finite element methods) to enable efficient simulation with a minimal number of semi-discretized equations, ensuring mass conservation even for non-linear diffusion problems involving variable diffusivities. These methods are then demonstrated by incorporation into the full p2D model, illustrating their advantages in simulating high C-rates and short-time dynamic operation of Lithium-ion batteries.

25 ENERGY STORAGE↗