Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum scale modeling”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Bridging molecular-scale interfacial science with continuum-scale models

Solid–water interfaces are crucial for clean water, conventional and renewable energy, and effective nuclear waste management. However, reflecting the complexity of reactive interfaces in continuum-scale models is a challenge, leading to oversimplified representations that often fail to predict real-world behavior. This is because these models use fixed parameters derived by averaging across a wide physicochemical range observed at the molecular scale. Recent studies have revealed the stochastic nature of molecular-level surface sites that define a variety of reaction mechanisms, rates, and products even across a single surface. To bridge the molecular knowledge and predictive continuum-scale models, we propose to represent surface properties with probability distributions rather than with discrete constant values derived by averaging across a heterogeneous surface. This conceptual shift in continuum-scale modeling requires exponentially rising computational power. By incorporating our molecular-scale understanding of solid–water interfaces into continuum-scale models we can pave the way for next generation critical technologies and novel environmental solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

A data-driven multiscale model for reactive wetting simulations

Here, we describe a data-driven, multiscale technique to model reactive wetting of a silver–aluminum alloy on a Kovar™ (Fe-Ni-Co alloy) surface. We employ molecular dynamics simulations to elucidate the dependence of surface tension and wetting angle on the drop’s composition and temperature. A design of computational experiments is used to efficiently generate training data of surface tension and wetting angle from a limited number of molecular dynamics simulations. The simulation results are used to parameterize models of the material’s wetting properties and compute the uncertainty in the models due to limited data. The data-driven models are incorporated into an engineering-scale (continuum) model of a silver–aluminum sessile drop on a Kovar™ substrate. Model predictions of the wetting angle are compared with experiments of pure silver spreading on Kovar™ to quantify the model-form errors introduced by the limited training data versus the simplifications inherent in the molecular dynamics simulations. The paper presents innovations in the determination of “convergence” of noisy MD simulations before they are used to extract the wetting angle and surface tension, and the construction of their models which approximate physio-chemical processes that are left unresolved by the engineering-scale model. Together, these constitute a multiscale approach that integrates molecular-scale information into continuum scale models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Temperature-dependent model of helium bubble expansion and bursting in tungsten

Continuum-scale models that can reliably predict the behavior of helium in tungsten are of interest to the fusion community due to the projected impact of these materials on fusion reactor operation. In this work, we perform molecular dynamics simulations of spherical helium bubbles of various sizes in tungsten at different temperatures and depths with the goal of determining a mathematical model of the pressure and density at which the bubbles initially expand or burst as a function of depth, size, temperature, and surface orientation. The bubbles are small enough that their loop-punching pressures cannot be accurately predicted with continuum mechanics, and their expansion behavior is important, as it appears to cause many of the features observed on helium-irradiated tungsten surfaces. We vary the temperature, bubble size, bubble depth, and surface orientation in each case, recording the bubble pressure and density that result in bubble expansion. An exponential function with three adjustable parameters is found to fit the results well; the parameters that best fit our results are provided.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Density functional theory-based surrogate kinetic models for heterogeneous reactions of hydrocarbon intermediates on silicon carbide

The increasing demand for high-performance materials in advanced technologies highlights the importance of achieving a fundamental understanding and potential control of silicon carbide (SiC) deposition processes. However, existing models often lack sufficient theoretical detail, relying heavily on empirical data and offering limited predictive capability. In particular, the complex surface chemistry governing SiC growth remains poorly understood. This study addresses these challenges by employing density functional theory (DFT) to investigate key heterogeneous reactions involving hydrocarbon intermediates on SiC surfaces, including dehydrogenation, hydrogenation, and carbon deposition. Transition state searches were conducted to identify reaction pathways and energy barriers. While first-principles calculations offer high accuracy, they are computationally intensive. To extend the utility of these first-principles results, vibrational analyses were performed using phonon-based statistical thermochemistry to compute temperature-dependent reaction rates which were used to develop Arrhenius-type surrogate kinetic models. Furthermore, the resulting framework provides a more rigorous, physically grounded basis for integrating atomistic insights into continuum-scale modeling, ultimately enabling improved prediction and optimization of SiC film growth in high-performance material systems.

Density Functional Theory↗
DOE OSTI2026

Microstructure, Transport, and Mechanics of Compacted Clay Simulated at the 0.1 μm Scale (1400 Smectite Clay Particles) Using a Coarse-Grained Model with Explicit Counterions

Clay-rich geomaterials play a critical role in many subsurface systems. The macroscale properties of these materials (low permeability, high ionic conductivity, high swelling pressure, etc.) are sensitive to molecular-level adsorption and hydration interactions at clay−water interfaces. Efforts to develop multiscale simulation approaches to predict these properties reveal a scale gap between atomistic simulations (typically limited to systems smaller than 10 nm) and continuum-scale models (which use computational grid elements with dimensions ≳ 10 μm). In this study, we present a coarse-grained (CG) framework that partly bridges this gap by simulating compacted smectite clay assemblages with dimensions of 0.1 μm containing 1,400 clay particles across a range of dry densities (1,050 to 1,850 kg·m −3 ) and Na/Ca counterion compositions (Na fraction ranging from 0.2 to 1). The simulated systems, along with their reconstructed binary three-dimensional pore networks, are used to evaluate the microstructure, pore size distribution, tortuosity, ion diffusivity, and swelling pressure of compacted smectite clay. Results show that our approach captures important features of the mesoscale heterogeneity of compacted clays, including tactoid formation, hierarchical porosity, and anisotropic pore networks. Results also reveal how compaction and counterion composition govern emergent behaviors, including dominant pore sizes, directional transport, and electrochemical response. This work highlights the potential of CG simulations to bridge molecular and continuum scales and to advance geotechnical and environmental applications involving clay-rich materials as well as related nanoporous media such as geopolymers and calcium-silicate-hydrate. However, the results also suggest that accurate prediction of certain microstructural and mechanical properties (e.g., swelling pressure) may require even larger-scale systems on the order of 1 μm.

02 PETROLEUM↗
DOE OSTI2021

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Large-scale molecular dynamics simulations of bubble collapse in water: Effects of system size, water model, and nitrogen

Molecular dynamics simulations in the microcanonical ensemble are performed to study the collapse of a bubble in liquid water using the single-site mW and the four-site TIP4P/2005 water models. To study system size effects, simulations for pure water systems are performed using periodically replicated simulation boxes with linear dimensions, L, ranging from 32 to 512 nm with the largest systems containing 8.7 × 10 6 and 4.5 × 10 9 molecules for the TIP4P/2005 and mW water models, respectively. The computationally more efficient mW water model allows us to reach converging behavior when the bubble dynamics results are plotted in reduced units, and the limiting behavior can be obtained through linear extrapolation in L –1 . Qualitative differences are observed between simulations with the mW and TIP4P/2005 water models, but they can be explained by the models’ differences in predicted viscosity and surface tension. Although bubble collapse occurs on time scales of only hundreds of picoseconds, the system sizes used here are sufficiently large to obtain bubble dynamics consistent with the Rayleigh–Plesset equation when using the models’ thermophysical properties as input. For the conditions explored here, extreme heating of the interfacial water molecules near the time of collapse is observed for the larger mW water systems (but the model underpredicts the viscosity), whereas heating is less pronounced for the TIP4P/2005 water systems because its larger viscosity contribution slows the collapse dynamics. The presence of nitrogen within the bubble only starts to affect bubble dynamics near the very end of the initial collapse, leading to an incomplete collapse and strong rebound for the mW water model. Although nitrogen is non-condensable at 300 K, it becomes highly compressed and reaches a liquid-like density near the collapse point. We find that the dissolution of nitrogen is much slower than the movement of the collapsing water front, and the re-expansion of the dense nitrogen droplet gives rise to bubble rebound. The incompatibility of the collapse and dissolution time scales should be considered for continuum-scale modeling of bubble dynamics. Finally, we also confirm that the diffusion coefficient for dissolved nitrogen is insensitive to pressure as the liquid transitions from a compressed to a stretched state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

High-Fidelity Arc-Discharge Model for Hydrogen-Plasma-Smelting-Reduction of Iron Ore

Electrification and use of renewable hydrogen is currently a necessity for decarbonizing the iron-and-steel industry. In this regard, hydrogen plasma smelting reduction (HPSR) is a novel pathway that is being explored for reduction of iron ore. HPSR provides several decarbonization merits compared to conventional blast furnaces. Firstly, the use of renewable hydrogen drastically reduces the CO2 emissions compared to the use of coke. Secondly, renewable electricity in the form of a thermal plasma for making reactive hydrogen species (radicals, ions) are more efficient at reducing iron ore compared to neutral H2. Thirdly, a molten product compatible with downstream processes is obtained from the intense heat transfer from the plasma. However, the scale-up of this technology requires fundamental exploration of hydrogen plasma dynamics and its interaction with complex solid material that include phase changing iron-ore and slag. In this work, we present a first principles continuum scale model for thermal plasmas in Ar/H2 gas mixtures typically used for HPSR. The thermal plasma governing equations for mass, momentum and energy with Lorentz force and Joule heating source terms are solved along with electromagnetic equations for electrostatic and magnetic vector potential. Our solver will be based on Pele, a suite of reacting flow solvers designed for advanced scientific computing architectures (Henry De Frahan et al., Proceedings of SIAM Parallel Processing, 13-25, 2024), and will utilize adaptive mesh generation for enhanced resolutions at locations of intense physicochemical interactions. This study will present the impact of Ar to H2 ratios on excited/dissociated hydrogen species concentrations, plasma temperature and conductivity along with the impact of outgassed species (water, metal vapor, O, OH radicals) from ore surface on gas phase chemistry. Furthermore, the heat and species flux to the surface will be quantified as a function of applied voltages in a transferred arc configuration.

hydrogen plasma↗
DOE OSTI2020

Enabling 6C Fast Charging of Li–Ion Batteries with Graphite/Hard Carbon Hybrid Anodes

Li-ion batteries that can simultaneously achieve high-energy density and fast charging are essential for electric vehicles. Graphite anodes enable a high-energy density, but suffer from an inhomogeneous reaction current and irreversible Li plating during fast charging. In contrast, hard carbon exhibits superior rate performance but lower energy density owing to its lower initial coulombic efficiency and higher average voltage. In this work, these tradeoffs are overcome by fabricating hybrid anodes with uniform mixtures of graphite and hard carbon, using industrially-relevant multi-layer pouch cells (>1 Ah) and electrode loadings (3 mAh cm –2 ). By controlling the graphite/hard carbon ratio, this study shows that battery performance can be systematically tuned to achieve both high-energy density and efficient fast charging. Pouch cells with optimized hybrid anodes retain 87% and 82% of their initial specific energy after 500 cycles of 4C and 6C fast-charge cycling, respectively. This is significantly higher than the 61% and 48% specific energy retention with graphite anodes under the same conditions. The enhanced performance is attributed to improved homogeneity of the reaction current throughout the hybrid anode, which is supported by continuum-scale modeling. Furthermore, this process is directly compatible with existing roll-to-roll battery manufacturing, representing a scalable pathway to fast charging.

25 ENERGY STORAGE↗
DOE OSTI2022

Computational Design of a Simple Flyer Plate Launcher

The response of materials to shock loading is important to understand for a variety of applications. When shock physics emerged during and after WWII, direct explosive loading or explosively driven plate impact was the primary tool for these studies. Subsequent decades have seen the widespread use of large caliber guns for plate impact studies, laser-shock facilities and pulsed power facilities. INL currently lacks a gun suitable for plate impact or explosives casting and machining facilities; however, it does possess explosives use and handling capabilities. An option for performing plate impact experiments was needed, therefore continuum scale models were utilized to explore a few simple donor-acceptor explosive plane wave lens designs, one of which could be hand packed with plastic explosives to launch flyer plates. 2D simulations were performed to study different geometries in an effort to minimize the difference in shock arrival across the central portion of a small copper flyer plate. A shock wave arrival time difference under 50 ns across 50 percent of the center of the flyer was achieved with a few designs. This work summarizes the computational models and results.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Establishing the Role of Operating Potential and Mass Transfer in Multicarbon Product Generation for Photoelectrochemical CO 2 Reduction Cells Using a Cu Catalyst

There is increasing interest in the possibility of photoelectrochemical (PEC) reduction of CO 2 to C 2+ products; however, the criteria for maximizing PEC solar-to-C 2+ (STC 2+ ) rates are not well understood. We report here a continuum-scale model of PEC CO 2 reduction (CO 2 R) on Cu in 0.1 M CsHCO 3 and use it to optimize the design and operating conditions for generating C 2+ products. Furthermore, we demonstrate that the potential-dependent product distribution of CO 2 R on Cu requires operating near the potential that maximizes C 2+ generation rates ($V$ id ), unlike PEC water splitting, which desires operation at the maximum photocurrent density. Because of this requirement, the criterion for a high STC 2+ rate includes high-photocurrent semiconductors with photovoltages near $V$ id and low series resistance. The STC 2+ rate in these systems is enhanced by optimal CO 2 transport and exhibits low sensitivity to dirunal solar irradiance variations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Hole formation effect on surface morphological response of plasma-facing tungsten

We report simulation results on the effect of helium (He) bubble bursting-mediated surface hole formation on the surface morphological response of tungsten plasma-facing components (PFCs) in nuclear fusion devices. Our analysis is based on an atomistically informed, continuum-scale model, which is capable of accessing the spatiotemporal scales relevant to the fuzz nanostructure formation process on the surface of PFC tungsten. Our simulations account, in an empirical fashion, for two types of subsurface bubble dynamical phenomena in the nanobubble region of PFC tungsten during He plasma irradiation, involving bubble bursting and surface crater formation. We demonstrate that the hole formation effect on the PFC tungsten surface accelerates the growth rate of nanotendrils and the onset of fuzz formation. As a result, the predicted incubation time for surface nanotendril growth is shortened, in agreement with experimental data of incubation fluence at comparable plasma exposure conditions. We also explore systematically the dependence of the PFC surface morphological response on the areal density of holes introduced at regular time intervals onto the He-implanted tungsten surface, a parameter in our analysis that serves as a proxy for the rate of He bubble bursting. More importantly, our simulations capture fine surface features in the PFC tungsten surface morphology and predict that the average spacing between nanotendrils is on the order of 100 nm, consistent with the experimental findings.

Chen, Chao-Shou (ORCID:0000000307910042)↗
DOE OSTI2023

Energy localization efficiency in 1,3,5-trinitro-2,4,6-triaminobenzene pore collapse mechanisms

Atomistic and continuum scale modeling efforts have shown that the shock-induced collapse of porosity can occur via a wide range of mechanisms dependent on pore morphology, the shockwave pressure, and material properties. The mechanisms that occur under weaker shocks tend to be more efficient at localizing thermal energy but do not result in high, absolute temperatures or spatially large localizations compared to mechanisms found under strong shock conditions. However, the energetic material 1,3,5-trinitro-2,4,6-triaminobenzene (TATB) undergoes a wide range of collapse mechanisms that are not typical of similar materials, leaving the collapse mechanisms and the resultant energy localization from the collapse, i.e., hotspots, relatively uncharacterized. Therefore, we present the pore collapse simulations of cylindrical pores in TATB for a wide range of pore sizes and shock strengths that trigger viscoplastic collapses that occur almost entirely perpendicular to the shock direction for weak shocks and hydrodynamic-like collapses for strong shocks that do not break the strong hydrogen bonds of the TATB basal planes. The resulting hotspot temperature fields from these mechanisms follow trends that differ considerably from other energetic materials; hence, we compare them under normalized temperature values to assess the relative efficiency of each mechanism to localize energy. The local intra-molecular strain energy of the hotspots is also assessed to better understand the physical mechanisms behind the phenomena that lead to a latent potential energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Helium bubble size effects on the surface morphological response of plasma-facing tungsten

Abstract We report a simulation study on the effects of helium (He) bubble size on the morphological evolution and pattern formation on the surface of tungsten used as a plasma-facing component (PFC) in nuclear fusion devices. We have carried out a systematic investigation based on self-consistent dynamical simulations of surface morphological evolution according to an atomistically-informed, 3D continuum-scale model that captures well the relevant length and time scales of surface nanostructure formation in PFC tungsten. The model accounts for PFC surface diffusion, driven by the biaxial compressive stress originating from the over-pressurized He bubbles in the near-surface region of PFC tungsten as a result of He plasma exposure, combined with the formation of self-interstitial atoms in tungsten that diffuse toward the PFC surface and the flux of surface adatoms generated as a result of surface vacancy-adatom pair formation upon He implantation; this transport of surface adatoms contributes to the anisotropic growth of surface nanostructural features due to the different rates of adatom diffusion along and across step edges of islands on the tungsten surface. Our detailed analysis reveals that varying the average He bubble size plays an important role in the PFC surface growth kinetics as well as the resulting surface topography. Specifically, we find that an increase in the He bubble size leads to a deceleration in the growth rate of the tungsten nanotendrils that emanate from the PFC surface. We also find that the separation distance between the resulting surface features increases with increasing He bubble size, as well as over time. This coarsening effect is a thermally activated process resulting in an accurate description of the temperature dependence of the average surface feature separation by an Arrhenius relation.

36 MATERIALS SCIENCE↗