Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stress 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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Fluidization of Group A Glass Particles: Experiments and Preliminary Validation

This work is aimed at providing reliable high-quality data from fluidization experiments. Geldart Group A glass beads are used as bed material. The fluidizing medium is air at atmospheric pressure which enters the system at 40% relative humidity. The air inlet velocity is varied over a broad range from 2 U mf to 26 U mf , where the minimum fluidization velocity, U mf , is 0.608 cm/s. There is not a significant change in the mean values of differential pressure; however, there are noticeable trends in their standard deviation values. Preliminary validation studies are performed using the open-source software Multiphase Flow with Interphase eXchanges (MFiX) developed by the National Energy Technology Laboratory (NETL). The results from Eulerian-Eulerian (MFiX-TFM, two-fluid model) and Eulerian-Lagrangian (MFiX-PIC, particle-in-cell) formulations are compared. The work presented in this report highlights the applicability of MFiX-PIC for large-scale applications, and reveals that PIC simulation offers a suitable alternative to TFM simulation when effects due to multiple components or polydispersity are significant. PIC modeling provides a readily extendible framework for such systems. However, PIC models do rely on an empirical closure for inter-particle stress. Some solution discrepancies were evident in MFiX-PIC predictions which require further analysis. As such, the study underlines the need for a systematic approach to quantify the sources of numerical uncertainty within the PIC model. However, even with known empiricism, MFiX-PIC has demonstrated considerable potential in analyzing gas-solid systems. The results obtained are encouraging and warrant further investigation to improve the predictive capability of MFiX-PIC.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unraveling fundamental mechanisms of silicon nitride crystallization in microelectronics manufacturing

This research project investigates the fundamental mechanisms of silicon nitride (SiN) crystallization, aiming to enhance the understanding of this critical material in microelectronics manufacturing. Through a collaborative effort between Sandia National Laboratories, the University of Tennessee, and the University of Florida, we developed a comprehensive framework that integrates experimental techniques, atomistic modeling, meso-scale simulations, and an integrated multi-scale model to capture this physical phenomenon on multiple time and length scales . The project developed a new machine learning based atomistic potential and utilized advanced phase field modeling to capture the complexities of polycrystalline growth and the influence of mechanical stresses on crystallization dynamics. By employing a grain tracker algorithm, the meso-scale model effectively identified and tracked individual crystal grains, enabling the simulation of anisotropic growth behaviors reflective of SiN’s physical properties. The integration of atomistic simulations with meso-scale modeling created a powerful multi-scale framework that validated atomistic inputs and enhanced predictive accuracy for crystallization dynamics at larger scales, validated experimentally. This adaptable modeling capability not only accelerates development times by informing manufacturing processes but also serves as a valuable starting point for understanding crystallization physics in similar materials. The insights gained from this research unlock new opportunities for the development of advanced materials tailored for future microelectronics and photonics applications. Overall, this project represents a significant advancement in understanding of fundamental physics of SiN and establishes a foundation for future research in material science, bridging the gap between atomic-level phenomena and macroscopic material behaviors for practical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Acute Ozone-Induced Transcriptional Changes in Markers of Oxidative Stress and Glucocorticoid Signaling in the Rat Hippocampus and Hypothalamus Are Sex-Specific

Exposure to a prototypic air pollutant ozone (O 3 ) has been associated with the activation of neuroendocrine stress response along with neural changes in oxidative stress (OS), inflammation, and Alzheimer’s disease-like pathologies in susceptible animal models. We hypothesized that neural oxidative and transcriptional changes induced by O 3 in stress responsive regions are sex-dependent. Male and female adult Long–Evans rats were exposed to filtered air or O 3 for two consecutive days (0.8 ppm, 4 h/day) and brain regions were flash-frozen. Activities of cerebellar OS parameters and mitochondrial complex I, II, and IV enzymes were assessed to confirm prior findings. We assessed transcriptional changes in hypothalamus (HYP) and hippocampus (HIP) for markers of OS, microglial activity and glucocorticoid signaling using qPCR. Although there were no O 3 or sex-related differences in the cerebellar activities of OS and mitochondrial enzymes, the levels of protein carbonyls and complex II activities were higher in females regardless of O 3 . There were no statistical differences in baseline expression of genes related to OS (Cat, Dhcr24, Foxm1, Gpx1, Gss, Nfe2l2, Sod1) except for lower HYP Sod1 expression in air-exposed females than males, and higher HIP Gss expression in O 3 -exposed females relative to matched males. Microglial marker Aif1 expression was higher in O 3 -exposed females relative to males; O 3 inhibited Itgam only in males. The expression of Bdnf in HIP and HYP was inhibited by O 3 in both sexes. Genes related to glucocorticoid signaling (Fkbp4, Fkbp5, Hsp90aa1, Hspa4, nr3c1, nr3c2) showed sex-specific effects due to O 3 exposure. Baseline expression of HIP Fkbp4 was higher in females relative to males. O 3 inhibited Nr3c1 in female HIP and male HYP, but Nr3c2 was inhibited in male HYP. Fkbp4 expression was higher in O 3 -exposed females when compared to matched males, whereas Fkbp5 was expressed at higher levels in both brain regions of males and females. These results indicate that sex-specific brain region responses to O 3 might, in part, be caused by OS and regulation of glucocorticoid signaling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Rewiring the unfolded protein response for plant growth recovery after stress

The unfolded protein response (UPR) is a highly coordinated signaling network that alleviates endoplasmic reticulum (ER) stress, a condition induced by diverse environmental challenges in plants. Over the past two decades, substantial progress has been made in elucidating the genetic and molecular mechanisms of ER stress sensing and signal transduction in plants, largely through studies in the model plant Arabidopsis thaliana . These advances have established the UPR as a central regulator of proteostasis and underscored its broader relevance to plant growth and development and crop productivity under stress conditions. Despite this progress, critical knowledge gaps remain, particularly concerning the downstream biological processes required for growth recovery once ER stress has subsided and how these processes are coordinated by UPR regulators. Recent systems-level and integrative studies have begun to reveal critical roles of UPR signaling in pathways governing growth re-establishment and homeostasis of nutrient allocation and energy metabolism. In this review, we highlight recent findings on the functional roles of the plant UPR in recovery from ER stress, with a focus on mechanisms mediated by UPR regulators and downstream biological pathways that enable the transition from stress mitigation to growth restoration. Although this research area is still emerging, accumulating evidence supports a model in which the UPR functions as a dynamic regulatory network that actively coordinates post-stress physiological recovery to support plant fitness.

ER stress↗
DOE OSTI2025

Phase-field modeling and experiments of dynamic fracture in single crystal quartz

Predicting the onset and characteristics of brittle fracture is important for a wide range of engineering and geological material applications. In this paper, we study important aspects of brittle fracture in α-quartz by phase-field modeling and experiments using a top-down approach. In the modeling framework, the work term in the Griffith energy balance is replaced with internal energy contributions that represent surface energy, thermal energy, and elastic strain energy stored in defects. This allows parametrization of individual energy contributions in terms of internal state variables and keeps track of energy partitioning after the onset of fracture. The path and history dependence of fracture is included in evolution laws for internal state variables, e.g., entropy evolution, while the energy remains a true potential. In the experimental part, dynamic compression experiments coupled with X-ray phase contrast imaging are performed on cube-like samples with a hole. In the top-down analysis, dynamic compression and three point bending experiments from the literature are simulated with the developed phase-field damage model. In conclusion, the fitted model highlights the strain rate, size, and stress state dependence of damage nucleation and evolution in single crystal α-quartz.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Predicting fatal heat and humidity using the heat index model

A unique wet-bulb temperature of 35°C is often used as the threshold for human survivability, but recent experiments have shown that a person’s core temperature starts to rise at a wide range of critical wet-bulb temperatures. Here, it is shown that the model underlying the heat index correctly predicts those critical wet-bulb temperatures, explaining 95% of the variance in the values observed in laboratory heat-stress experiments. This is the first time the heat-index model has been validated against physiological data from laboratory experiments. For light and moderate exertion in an indoor setting, the heat index model predicts that the critical wet-bulb temperature ranges from 20 to 32°C, depending on the relative humidity, consistent with experimental results. Furthermore, for the same setting and exertion, the heat index model predicts fatal wet-bulb temperatures ranging from 24 to 37°C. Recent experiments have identified the critical combinations of heat and humidity, in an indoor setting, above which an individual is unable to maintain a standard core temperature, indicating severe heat stress. It is shown here why this state of severe heat stress cannot be predicted using the wet-bulb temperature. Instead, it is shown that the recently extended heat index model can explain nearly all of the variance in the observed critical combinations of temperature and humidity, and can be used to calculate fatal combinations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

High Throughput Computational Framework of Materials Properties for Extreme Environments

This project aims to establish a framework capable of efficiently predicting the properties of structural materials for service in harsh environments over a wide range of temperatures and over long periods of time. The approach is to develop and integrate high throughput first-principles calculations in combination with machine learning (ML) methods, perform high throughput CALPHAD (calculations of phase diagrams) modeling, and carry out finite element method (FEM) simulations. Relevant to high temperature service in fossil power system, nickel-based superalloys such as Inconel 740 and Haynes 282 as well as the associated (Ni-Cr-Co)-Al-C-Fe-Mn-Mo-Nb-Si-Ti system, were investigated. The present framework was built on the concept of phase-based property data, in which properties of individual phases are modeled as a function of internal and external independent variables. This project established an open-source infrastructure with the following capabilities: (1) High throughput implementation of first-principles calculations at finite temperatures and variable compositions using both accurate phonon calculations and the efficient Debye model for thermodynamic properties, elastic constants, diffusion coefficients, vacancy formation, stacking and twin faults, and dislocation mobility; i.e., using the developed code DFTTK; (2) Machine learning capabilities to predict the above properties so that the number of first-principles calculations can be significantly reduced; e.g., using the developed code SIPFENN; (3) High throughput CALPHAD modeling of the above properties as a function of temperature and composition using our unique capability based on ESPEI and PyCalphad; (4) New capabilities to predict the stress-strain behavior of individual phases; and (5) New models for tensile strength prediction in common FEM software with the crystal plasticity finite element simulations (CPFEM).

, Ni-based superalloys↗
DOE OSTI2026

X-ray tomography of damage dynamics in advanced materials using a laser wakefield accelerator

Additively manufactured (AM) metals offer the potential for customizable, cost-effective components, but qualification and certification are crucial. Key to this process is understanding pore dynamics under stress, typically analyzed using micro-computed tomography. This study introduces laboratory-scale “betatron” x-rays from laser wakefield acceleration as a high-throughput alternative for x-ray tomography of advanced materials, such as AM AlSi10Mg alloys. Coupled with 3D finite element modeling, this method provides detailed insights into stress-porosity interactions. The approach delivers high-resolution scans, revealing that pore shape and local triaxiality significantly influence fracture dynamics, supporting advanced material characterization. This work also demonstrates the potential and versatility of laser-betatron x-ray μCT for generating large datasets to accelerate our understanding of the stochastic, process-specific nature of pore formation in AM alloys.

Senthilkumaran, Vigneshvar↗
DOE OSTI2022

A directional microcrack damage theory for brittle solids based on continuous hyperplasticity

We present a novel microcrack-damage theory for brittle solids under compression. Instead of using internal variables like zeroth, second or fourth rank damage tensors, the state of material damage is represented by an internal function that encapsulates the information of direction, density and size of microcracks. Just like other internal variables, the evolution of this state function must obey the second law of thermodynamics for arbitrary loading paths. This is done by casting the model in the framework of continuous hyperplasticity and enforcing a non-negative dissipation rate functional. The proposed framework offers predictions on the continuous evolution of microcrack density and the induced material anisotropy along with the macroscopic stress-strain curves. The use of continuous damage function grants the model significantly enhanced resolution in characterizing the direction-dependent response of cracked solids compared to classical models that are based on damage tensors. Two scenarios are considered in developing the theory, one assumes frictionless cracks and the other incorporates friction between crack surfaces. The results highlight that inelasticity, pressure dependence, and loading-unloading hysteresis exhibited by brittle solids are natural consequences of frictional microcracks. The proposed theory offers a generic and versatile framework to upscale micromechanical processes operating at individual crack scale to explain the macroscopic behavior of cracked solids.

Materials Science↗
DOE OSTI2022

Scaling of Shear Rheology of Concentrated Charged Colloidal Suspensions across Glass Transition

Electrostatic interparticle interactions are a key component in controlling and designing rheological characteristics of concentrated charged colloidal suspensions. Herein, we investigate electroviscous effects on shear rheology using highly charged silica particles. By fixing the volume fraction but varying the salinity, the system undergoes a glass transition as evidenced by the evolution of the yield stress and zero-shear viscosity. We show that the steady shear viscosities obey a critical scaling relation that scales the flow curves into a super- and a sub-critical branch with glass transition salinity serving as the bifurcation point; we also demonstrate an isoviscosity scaling that collapses all isoviscosity lines into a single master curve that exhibits no singularity. Based on each scaling relation, in conjunction with common modeling equations, the quantitative relationships between the shear viscosity, stress, and salinity are established. This study demonstrates a new framework to model the steady shear rheology of concentrated charged colloids.

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

Deformation twinning and detwinning in extruded Mg-4Al: In-situ experiment and crystal plasticity simulation

Deformation twinning and detwinning in extruded Mg-4Al were investigated using in-situ SEM-DIC experiments and crystal plasticity finite element (CPFE) simulation. In this study, the in-situ SEM-DIC method was used to provide a unique set of data including twin/detwin characteristics and twin area fraction in addition to strain maps. A statistical analysis of the activation of twin variants and twin area fraction during both twinning and detwinning was conducted. A strong correlation was found between twin growth/shrinkage and the Schmid Factor (SF) for individual twin variants Higher twin SF during loading and unloading led to higher twin growth and shrinkage, respectively. However, after the applied compressive strain was removed, the pattern of the twin area fractions of the residual twin variants versus their nominal SFs did not follow the trend observed at the maximum compressive strain. Using a systematic methodology and an advanced twin/detwin model, the PRISMS-Plasticity CPFE simulation was calibrated using experimentally determined stress versus strain and twin area fraction versus strain information. A comprehensive evaluation of the CPFE model was conducted to determine its ability to capture the statistics of twin variants activation and twin area fraction. CPFE accurately captured the statistical aspects of both twinning and detwinning. Finally, it also predicted the first dominant twin variant for 47.5% of the grains and at least one of the two dominant twin variants for 80% of the grains at maximum compressive strain.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Numerical Modeling Framework for Flocculation and Cohesive Sediment Transport in the Wave Bottom Boundary Layer

Flocculation, a critical process in coastal and estuarine systems, plays a significant role in sediment transport, nutrient cycling, and ecological health. This study develops a cohesive sediment transport modeling framework tailored to the wave bottom boundary layer under dilute and equilibrium conditions, explicitly incorporating flocculation effects via a Population Balance Equation (PBE). Using Direct Numerical Simulation, six baseline cases, each with a distinct sediment concentration profile resulting from a constant settling velocity and critical erosion shear stress, are generated to drive the PBE flocculation model for given floc yield strength and stickiness. Results reveal that flocculation significantly influences sediment concentration profiles promoting three distinct stages, well‐mixed, transition to lutocline, and well‐developed lutocline. At low concentrations with well‐mixed profiles, cohesive floc properties are less significant, and turbulence is a main flocculation driver. In contrast, as concentration increases, cohesive floc properties become crucial, facilitating lutocline formation. Here, the analysis also highlights limitations of depth‐averaged settling velocity as a parameterization. It is suitable for well‐mixed and transitional profiles but fails in well‐developed lutoclines, where empirical formulations that explicitly incorporate turbulent shear rate and sediment concentration better capture variability. This study underscores the necessity of incorporating flocculation effects into sediment transport models to enhance predictions of sediment dynamics in wave bottom boundary layers.

Penaloza‐Giraldo, Jorge A. [Oak Ridge National Lab↗
DOE OSTI2022

Hydromechanical impact of basement rock on injection-induced seismicity in Illinois Basin

Abstract The common explanation of observed injection-induced microseismicity is based on the measured stress state at the injection interval and the assumption that it remains the same in the vicinity. We argue here that representing the stress state in different geologic formations over the injection site with the single Mohr’s circle is insufficient due to local stratigraphic features and contrast in compressibilities of the involved formations. The role of hydromechanical coupling in the microseismic response is also crucial for the proper assessment of the problem. Thoroughly monitored Illinois Basin Decatur Project revealed the majority of CO 2 injection-associated microseismic events being originated in the crystalline basement. Even though basement faults can serve as the conduits for fluid flow—the predicted pressure increase seems to be insufficient to trigger seismicity. To address this issue, accurate laboratory measurements of rock properties from the involved formations are conducted. The pre-injection stress state and its evolution are evaluated with the hydromechanically coupled numerical model. It appears that the presence of an offset in a stiff competent layer affects the stress state in its vicinity. Therefore, both the pre-injection stress state and its evolution during the fluid injection should be addressed during the induced seismicity assessment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Data from: "Warming and the dependence of limber pine (Pinus flexilis) establishment on summer soil moisture within and above its current elevation range"

This data package contains data that were used for analysis in “Warming and the dependence of limber pine (Pinus flexilis) establishment on summer soil moisture within and above its current elevation range”, by Moyes et al. 2013. All data collection and field research were completed on Niwot Ridge, Colorado, USA.This data package contains nine comma-separated-values (.csv) files, one text (.txt) file, and two zipped seedling folders that were used for leaf area analysis. One zipped folder contains 468 .jpg photographs of seedlings, and the second contains 468 corresponding Image J-processed .jpg images that include silhouette leaf area values. .csv and .txt files can be opened using any compatible simple text-editor software such as TextEdit (Mac) and Notepad (Windows); .csv’s can also be opened using R and Microsoft Excel. Image files can be opened using Preview (Mac) and Photos (Windows). In addition, there are a total of 31 Microsoft Excel files: three .xlsx files, and 28 raw Li-Cor output .xls files. This data user’s guide is available in .pdf format, and can be opened using Adobe Acrobat Reader, or any other compatible file viewing software. Geospatial data showing field site locations are also included in the archive for use and reference. There are two geospatial formats in this archive: ESRI shapefiles (.shp) and keyhole markup-language (.kml) files. Both file types contain bounding box information, with the former being polygons, and the latter containing corner coordinates for each site. ESRI shapefiles can be opened using any geospatial software compatible with the file type (such as ESRI’s ArcGIS suite and QGIS), and .kml files are compatible with Google Earth and Google Maps.--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------Continued changes in climate are projected to alter the geographic distributions of plant species, in part by affecting where individuals can establish from seed. We tested the hypothesis that warming promotes uphill redistribution of subalpine tree populations by reducing cold limitation at high elevation and enhancing drought stress at low elevation. We seeded limber pine (Pinus flexilis) into plots with combinations of infrared heating and water addition treatments, at sites positioned in lower subalpine forest, the treeline ecotone, and alpine tundra. In 2010, first-year seedlings were assessed for physiological performance and survival over the snow-free growing season. Seedlings emerged in midsummer, about 5–8 weeks after snowmelt. Low temperature was not observed to limit seedling photosynthesis or respiration between emergence and October, and thus experimental warming did not appear to reduce cold limitation at high elevation. Instead, gas exchange and water potential from all sites indicated a prevailing effect of summer moisture stress on photosynthesis and carbon balance. Infrared heaters raised soil growing degree days (base 5 °C, p < 0.001) and August–September mean soil temperature (p < 0.001). Despite marked differences in vegetation cover and meteorological conditions across sites, volumetric soil moisture content (θ) at 5–10 cm below 0.16 and 0.08 m^3 m^(−3) consistently corresponded with moderate and severe indications of drought stress in midday stem water potential, stomatal conductance, photosynthesis, and respiration. Seedling survival was greater in watered plots than in heated plots (p = 0.01), and negatively related to soil growing degree days and duration of exposure to θ < 0.08 m^3 m^(−3) in a stepwise linear regression model (p < 0.0001). We concluded that seasonal moisture stress and high soil surface temperature imposed a strong limitation to limber pine seedling establishment across a broad elevation gradient, including at treeline, and that these limitations are likely to be enhanced by further climate warming.--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------Maintenance log:Amended file name error in Data User's Guide May 19 2022

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Comparative transcriptomics and metabolomics reveal specialized metabolite drought stress responses in switchgrass ( Panicum virgatum )

Summary Switchgrass ( Panicum virgatum ) is a bioenergy model crop valued for its energy efficiency and drought tolerance. The related monocot species rice ( Oryza sativa ) and maize ( Zea mays ) deploy species‐specific, specialized metabolites as core stress defenses. By contrast, specialized chemical defenses in switchgrass are largely unknown. To investigate specialized metabolic drought responses in switchgrass, we integrated tissue‐specific transcriptome and metabolite analyses of the genotypes Alamo and Cave‐in‐Rock that feature different drought tolerance. The more drought‐susceptible Cave‐in‐Rock featured an earlier onset of transcriptomic changes and significantly more differentially expressed genes in response to drought compared to Alamo. Specialized pathways showed moderate differential expression compared to pronounced transcriptomic alterations in carbohydrate and amino acid metabolism. However, diterpenoid‐biosynthetic genes showed drought‐inducible expression in Alamo roots, contrasting largely unaltered triterpenoid and phenylpropanoid pathways. Metabolomic analyses identified common and genotype‐specific flavonoids and terpenoids. Consistent with transcriptomic alterations, several root diterpenoids showed significant drought‐induced accumulation, whereas triterpenoid abundance remained predominantly unchanged. Structural analysis verified select drought‐responsive diterpenoids as oxygenated furanoditerpenoids. Drought‐dependent transcriptome and metabolite profiles provide the foundation to understand the molecular mechanisms underlying switchgrass drought responses. Accumulation of specialized root diterpenoids and corresponding pathway transcripts supports a role in drought stress tolerance.

transcriptomics↗
DOE OSTI2025

Development and Validation of MALAMUTE model for Electric Field Assisted Sintering of Structural Materials

Fusion power plant designs feature extreme material performance requirements for structural material candidates. In addition to conventional alloys, more advanced composites and oxide dispersion strengthened (ODS) alloys are being explored, however, achieving the desired microstructures to maximize performance using traditional manufacturing methods can be challenging. The advanced manufacturing (AM) electric field-assisted sintering (EFAS) technique offers improved control over the final microstructure through higher heating and cooling rates and moderate pressures. Modeling and simulation tools show promise in elucidating the process-structure-property-performance (PSPP) correlation for AM-produced parts, including the EFAS process. An inherently multiscale process, the EFAS technique aligns well with the multiscale modeling capability of the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE)[cite]. We present here an electro-thermo-mechanical approach to modeling the EFAS process using the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE) code. Prediction of the field and gradient distributions across the EFAS tooling is required to accurately describe the conditions for the lower-scale microstructural evolution models. In this work we present the MALAMUTE model developed to predict the electrical potential, temperature, and mechanical stress distribution across the EFAS graphite tooling and part at the larger engineering-scale. Validation of the MALAMUTE engineering-scale model is completed using data from experimental densification and pre-densified runs of iron powder via EFAS at 1000oC. These runs were conducted using a Thermal Technology DCS-5 EFAS system. Data collected during the experiment runs include the direct current (DC) supplied to the graphite tooling, the temperature of the graphite tooling as measured with a pyrometer, and the force applied to the top of the graphite tooling stack, and the data were recorded every 10 seconds. Our validation approach used the current and force data from the EFAS run as boundary condition inputs to the MAMALUTE simulation; the temperature data were used to evaluate the MALAMUTE EFAS model prediction. Results of the MALAMUTE simulations are employed to connect the external pyrometer temperature measurement to the temperature profile across the part undergoing consolidation. We investigate the impact of material property variation and mesh deformation on the temperature profile as predicted by MALAMUTE. We conclude by highlighting projects where the MALAMUTE EFAS modeling and simulation capabilities will be used to assist experimental design.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Multiscale Modeling Framework for Lithium Nucleation in 3D Porous Carbon Anodes

Porous carbon scaffolds offer a promising route for mitigating non-uniform lithium (Li) plating to enhance the safety and longevity of Li metal batteries. However, the influence of microstructural morphology on Li nucleation is not well understood. Here, we present a multiscale modeling framework to investigate how the porous microstructure of carbon materials affects Li nucleation behavior. Ab initio molecular dynamics simulations quantify the nucleation energy barriers of Li on graphene as a function of Li content, surface curvature, and applied potential, providing key parameters for a classical nucleation theory (CNT) model. From macroscale half-cell simulations, we obtained Li concentration and electrical potential profiles to define boundary conditions for mesoscopic simulations. At the mesoscale, three distinct synthetic 3D microstructures with different porosities and characteristic feature sizes are generated to resolve local distributions of Li flux, current density, and mechanical stress. These outputs are integrated into the CNT model to map spatial variation in nucleation rates. Our findings reveal trade-offs between suppressing nucleation rates and achieving spatial uniformity, offering design guidelines for optimizing porous carbon anodes to balance nucleation control and mechanical integrity.

Materials science↗
DOE OSTI2022

Deep Learning for Subgrid‐Scale Turbulence Modeling in Large‐Eddy Simulations of the Convective Atmospheric Boundary Layer

Abstract In large‐eddy simulations, subgrid‐scale (SGS) processes are parameterized as a function of filtered grid‐scale variables. First‐order, algebraic SGS models are based on the eddy‐viscosity assumption, which does not always hold for turbulence. Here we apply supervised deep neural networks (DNNs) to learn SGS stresses from a set of neighboring coarse‐grained velocity from direct numerical simulations of the convective boundary layer at friction Reynolds numbers Re τ up to 1243 without invoking the eddy‐viscosity assumption. The DNN model was found to produce higher correlation between SGS stresses compared to the Smagorinsky model and the Smagorinsky‐Bardina mixed model in the surface and mixed layers and can be applied to different grid resolutions and various stability conditions ranging from near neutral to very unstable. The DNN model can capture key statistics of turbulence in a posteriori (online) tests when applied to large‐eddy simulations of the atmospheric boundary layer.

54 ENVIRONMENTAL SCIENCES↗