Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale simulations”

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 541 records · Page 30

Electro-reduction of organics on metal cathodes: A multiscale-modeling study of benzaldehyde on Au (111)

Here, we present a multiscale modeling study on the electro-reduction of benzaldehyde (BZY) on Au(111) under realistic electrochemical conditions. To model the electrochemical cell, we adopt a capacitor model in which complex solvents are confined between a Au cathode and a carbon anode. Classical molecular dynamics simulations reveal that electrode charge density and the presence of alcohol show strong effects on the density, adsorption geometry and dynamics of benzaldehyde on the Au electrode. Under charging conditions, the surface concentration of benzaldehyde on the Au electrode decreases, while the content of other species increases. Finally, we proposed a scheme that correlates the electric current running through the Au/solvent interface with the applied bias. This study provides a molecular level understanding of how solvent composition, in this case water/alcohol content, controls the activity of electrocatalytic hydrogenation

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A multiscale anisotropic polymer network model coupled with phase field fracture

The study of polymers has continued to gain substantial attention due to their expanding range of applications, spanning essential engineering fields to emerging domains like stretchable electronics, soft robotics, and implantable sensors. These materials exhibit remarkable properties, primarily stemming from their intricate polymer chain network, which, in turn, increases the complexity of precisely modeling their behavior. Especially for modeling elastomers and their fracture behavior, accurately accounting for the deformations of the polymer chains is vital for predicting the rupture in highly stretched chains. Despite the importance, many robust multiscale continuum frameworks for modeling elastomer fracture tend to simplify network deformations by assuming uniform behavior among chains in all directions. Recognizing this limitation, our study proposes a multiscale fracture model that accounts for the anisotropic nature of elastomer network responses. At the microscale, damage in the chains is assumed to be driven by both the chain's entropy and the internal energy due to molecular bond distortions. In order to bridge the stretching in the chains to the macroscale deformation, we employ the maximal advance path constraint network model, inherently accommodating anisotropic network responses. As a result, chains oriented differently can be predicted to exhibit varying stretch and, consequently, different damage levels. To drive macroscale fracture based on damages in these chains, we utilize the micromorphic regularization theory, which involves the introduction of dual local-global damage variables at the macroscale. The macroscale local damage variable is obtained through the homogenization of the chain damage values, resulting in the prediction of an isotropic material response. The macroscale global damage variable is subjected to nonlocal effects and boundary conditions in a thermodynamically consistent phase field continuum formulation. Moreover, the total dissipation in the system is considered to be mainly due to the breaking of the molecular bonds at the microscale. To validate our model, we employ the double-edge notched tensile test as a benchmark, comparing simulation predictions with existing experimental data. Additionally, to enhance our understanding of the fracturing process, we conduct uniaxial tensile experiments on a square film made up of polydimethylsiloxane (PDMS) rubber embedded with a hole and notches and then compare our simulation predictions with the experimental observations. Furthermore, we visualize the evolution of stretch and damage values in chains oriented along different directions to assess the predictive capacity of the model. In conclusion, the results are also compared with another existing model to evaluate the utility of our model in accurately simulating the fracture behavior of rubber-like materials.

42 ENGINEERING↗

Sensitivity of pore collapse heating to the melting temperature and shear viscosity of HMX

A multiscale modeling strategy is used to quantify factors governing the temperature rise in hot spots formed by pore collapse from supported and unsupported shock waves in the high explosive HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Two physical aspects are examined in detail, namely the melting temperature and liquid shear viscosity. All-atom molecular dynamics simulations of phase coexistence are used to predict the pressure-dependent melting temperature up to 5 GPa. Equilibrium simulations and the Green–Kubo formalism are used to obtain the temperature- and pressure-dependent liquid shear viscosity. Starting from a simplified continuum-based grain-scale model of HMX, in this study we systematically increase the complexity of treatments for the solid–liquid phase transition and liquid shear viscosity in simulations of pore collapse. Using a realistic pressure-dependent melting temperature completely suppresses melting for supported shocks, which is otherwise predicted when treating it as a constant determined at atmospheric pressure. Alternatively, melt pools form around collapsed pores when the pressure (and melting temperature) are reduced during the release stage of unsupported shocks. Capturing the pressure dependence of the shear viscosity increases the peak temperature of melt pools by hundreds of Kelvin through viscous work. The complicated interplay of the solid-phase plastic work, solid–liquid phase transition, and liquid-phase viscous work identified here motivate taking a systematic approach to building increasingly complex grain-scale models.

36 MATERIALS SCIENCE↗

Numerical investigation of thermomechanical response of multiscale porous Ultra-High Temperature Ceramics

Recent advances in Ultra-High Temperature Ceramics (UHTC) manufacturing have permitted the development of multiscale porous UHTC microstructures. Within the target application of hypersonic vehicles, dense UHTCs are suitable for thermal protection on leading edges, whereas porous UHTCs may find a role in providing thermally insulated interfaces for temperature-sensitive interior components. Designing vehicles incorporating porous UHTCs requires a characterization of their thermomechanical properties across the full range of expected operating temperatures spanning –20 °C–2500 °C. This research represents a preliminary study in performing this characterization. Several numerical experiments are performed using a coupled thermomechanical implementation of the Material Point Method to determine the temperature dependence of effective material properties both with and without damage. Furthermore, complex time-dependent boundary conditions derived from known hypersonic flight profiles are simulated in order to probe the various couplings between deformation, damage and heat transfer. It is shown that the model reveals the importance of micro-buckling in determining effective material stiffness and thermal conductivity.

36 MATERIALS SCIENCE↗

Identifying the Growth Phase of Magnetic Reconnection Using Pressure‐Strain Interaction

Abstract Magnetic reconnection often initiates abruptly and then rapidly progresses to a nonlinear quasi‐steady state. While satellites frequently detect reconnection events, ascertaining whether the system has achieved steady‐state or is still evolving in time remains challenging. Here, we propose that the relatively rapid opening of the reconnection separatrices within the electron diffusion region serves as an indicator of the growth phase of reconnection. The opening of the separatrices is produced by electron flows diverging away from the neutral line downstream of the X‐line and flowing around a dipolarization front. This flow pattern leads to characteristic spatial structures in the electron pressure‐strain interaction that could be a useful indicator for the growth phase of a reconnection event. We employ two‐dimensional particle‐in‐cell numerical simulations of anti‐parallel magnetic reconnection to validate this prediction. We find that the signature discussed here, alongside traditional reconnection indicators, can serve as a marker of the growth phase. This signature is potentially accessible using multi‐spacecraft single‐point measurements, such as with NASA's Magnetospheric Multiscale satellites in Earth's magnetotail. Applications to other settings where reconnection occurs are also discussed.

Barbhuiya, M. Hasan [Department of Physics and Ast↗

Topology induced crossover between Langevin, subdiffusion, and Brownian diffusion regimes in supercooled water

Despite extensive studies of supercooled water, it remains challenging to understand its peculiar dynamic anomalous properties. In this work, we integrated full atomistic simulations of supercooled water over the temperature range of room temperature to 200 K using a quantum-mechanics-based polarizable force field with the dressed dynamics method that couples fast collision events and slow reorganization dynamics of hydrogen-bond networks. Our analysis unveils the salient multiscale features in the transient relaxation dynamics of supercooled water. A classical Langevin behavior dominates at fast timescales, while long-time relaxations unveil two different activation barriers in two temperature regions: below and above 230 K. The modulation of the entropy spectrum by temperature is elucidated in terms of a three-state model underlined by the complexity of the water dynamics associated with a topological transition of a strong hydrogen-bond network. This state-dependent network topology is quantitatively characterized by power-law exponents of inverse network connectivity from 200 to 298 K. Finally, this work provides valuable guidance for further studies on the transient relaxation dynamics of supercooled water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impacts of Biomass Feedstock Pre-Processing on Heat and Mass Transfer During Pyrolysis Using X-Ray Computed Tomography and Multiscale Modeling

Knowledge of the transport properties of biomass particles such as porosity, tortuosity, and permeability is paramount for high-fidelity modeling of biomass pyrolysis due to the heat and mass transfer limitations imposed by particle microstructure. X-ray computed tomography (XCT) is a non-destructive imaging method that enables full 3D reconstructions of the biomass particle microstructure with high resolution, permitting direct calculation of porosity, tortuosity, and permeability from real particle geometries. In this study, XCT imaging revealed the 3D microstructures of particles and chars from pyrolytic conversion of cylindrically cut or milled/pelletized loblolly pine samples. The porosity, tortuosity, and permeability were calculated directly from the XCT geometries via open-source microstructural analysis tool MATBOX+TauFactor (https://github.com/NREL/MATBOX_Microstructure_analysis_toolbox) and computational fluid dynamics (CFD) simulations using our solver, Mesoflow (https://github.com/NREL/mesoflow). These properties were used in a reactor scale model developed in COMSOL of the single particle reactor at NREL to investigate the impact of feedstock pre-processing on biomass conversion during pyrolysis with rigorous experimental validation.

biomass↗

Hydrological Controls on Riverbed Methane Emissions: A Numerical Investigation of Hydrodynamic and Ebullitive Mechanisms from Site to Basin Scales

Methane production and emission from riverbed sediments constitute a significant yet underexplored component of the river methane budget. Despite their importance, a mechanistic and quantitative understanding of these processes under dynamic flow conditions remains elusive, particularly at the basin scale. In this study, we investigate the multiscale (site-to-basin) mechanisms governing methane emissions from riverbeds, focusing on the interplay between biogeochemical processes and hydrological dynamics. We developed a reactive transport model that integrates methane production with hydrodynamic and ebullitive pathways and coupled it with a basin-scale groundwater flow model to simulate methane emissions across scales. By incorporating key factors such as vertical hydrological exchange flow (VHEF), particulate organic carbon availability, sediment properties, and temperature sensitivity, our model captures the spatiotemporal variability of methane ebullition at both the site and basin scales. Our results reveal that methane ebullition is strongly modulated by VHEF, with lower flow rates promoting methane accumulation and subsequent ebullition. Sensitivity analyses underscore nonlinear responses of methane fluxes to environmental drivers. Furthermore, these findings emphasize the critical roles of sediment composition, hydrological exchange, and environmental conditions in shaping methane dynamics. Overall, this work advances our understanding of riverine methane emissions and offers insights for guiding climate models and mitigation strategies.

Fluxes↗

The Method of Finite Averages

The Method of Finite Averages (MoFA) is a rigorous multiscale modeling methodology for efficiently modeling multi-physical phenomena in heterogeneous porous media. The code developed in this project aims to perform the numerical calculations required to formulate, implement, and verify MoFA models for Earth and Energy systems (i.e., model verification refers to performing fully-resolved simulations of the systems and comparing their results to those of the models). In general, MoFA transforms partial differential equations (PDEs) describing the fine-scale physics of a system into coupled ordinary differential equations (ODEs)---in time---that describe the coarse-scale---or "average"---physical behaviors of the system. This transformation significantly expedites system simulation, as the coarse-scale ODEs involve vastly fewer degrees of freedom than the fine-scale PDEs. The code developed under this project will allow users to 1.) generate system geometries and numerical meshes, 2.) solve the PDE and ODE systems required for MoFA model formulation and implementation, 3.) solve the PDE systems required to obtain fully-resolved simulation results for model verification, and 4.) compare and plot results (e.g., the model and fully-resolved simulation solutions, the error between the solutions, etc.).

Pietrzyk, KyleM [Lawrence Livermore National Labor↗

Nanosecond Laser Drying in Wafer Foundries for Significant Energy Efficiency Improvement

Surface drying is a critical process in microelectronics wafer fabrication. In order to reduce and eliminate watermarks, many surface drying processes have been studied and used. Most of the wafer foundries currently use drying processes based on Isopropyl Alcohol (IPA), which causes serious safety (fire), health, environmental, and energy efficiency issues. The long-term goal of this project is to develop a fast, effective, chemical-free, and extremely energy efficient wafer drying technology that employs nanosecond laser-induced sub-surface evaporation and explosion. The objectives of this project are to 1) experimentally develop and verify the laser surface drying method, explore and identify how and to what extent various physical parameters influence the drying efficiency and effectiveness; and 2) understand and optimize the underlying thermodynamics in the proposed laser surface drying, including micro/nanoscale sub-surface superheating, evaporation, explosion, and photo/thermochemical bond breaking via atomistic modeling and experimental characterization. Toward these objectives, five tasks have been accomplished successfully, including 1) develop in-lab laser surface drying setup and study the effect of varied laser parameters, 2) use multiscale hybrid-modeling to study and understand the physics of water behavior during drying toward process optimization, 3) investigate the effect of drying conditions on the drying effectiveness, 4) develop a large-scale laser drying scanning system with automatic scanning, and 5) conduct large-scale simulation to study the effect of dopant level and develop strategies for system development toward industry deployment.

42 ENGINEERING↗

Tailoring ductile-phase toughened tungsten hierarchical microstructures for plasma-facing materials

Biological materials, such as bones, nacre, etc., have been found to show attractive and unique combinations of stiffness, strength and fracture toughness. Research has been conducted in various engineering areas to mimic the naturally hierarchical microstructures or nanostructures of these materials to produce materials with optimum stiffness, high strength and fracture toughness for their intended structural applications. In the same objective, this work applies a multiscale microstructural approach recently developed to investigate the deformation and fracture behavior of ductile phase toughened tungsten (W) materials such as tungsten-nickel/iron (W-Ni/Fe) composites that possess lamellar-like and hierarchical “brick-and-mortar” (BAM) microstructures. First, the approach is used to simulate tensile loading of W-Ni/Fe specimens cut out from hot-rolled W-Ni/Fe plates which possess lamellar-like microstructures. The finite element model of the gage section of the specimen consists of a W-Ni/Fe dual-phase microstructural domain in which the constitutive behaviors of W and Ni/Fe phases are described by an elastic-plastic damage model. The predicted material stress-strain response and crack pattern development as a function of loading are then compared to the corresponding experimental results to determine the constitutive model parameters for W and Ni/Fe. Subsequently, the model parameters are used to analyze W-Ni/Fe BAM microstructures that are artificially created to investigate the effects of microstructural features and morphologies on the composite stress-strain response, damage, and fracture patterns. The modeling shows that regular BAM microstructures that experience bridging mechanism combined with crack penetration across W-phase regions exhibit significantly higher strengths than the rather random lamellar-like microstructures. More significantly, however adjusting the brick’s length-to-height ratio, the BAM microstructure can be designed to allow a more distributed damage zone which leads to increased strength, ductility and fracture energy.

36 MATERIALS SCIENCE↗

Learning nonlinear operators in latent spaces for real-time predictions of complex dynamics in physical systems

Abstract Predicting complex dynamics in physical applications governed by partial differential equations in real-time is nearly impossible with traditional numerical simulations due to high computational cost. Neural operators offer a solution by approximating mappings between infinite-dimensional Banach spaces, yet their performance degrades with system size and complexity. We propose an approach for learning neural operators in latent spaces, facilitating real-time predictions for highly nonlinear and multiscale systems on high-dimensional domains. Our method utilizes the deep operator network architecture on a low-dimensional latent space to efficiently approximate underlying operators. Demonstrations on material fracture, fluid flow prediction, and climate modeling highlight superior prediction accuracy and computational efficiency compared to existing methods. Notably, our approach enables approximating large-scale atmospheric flows with millions of degrees, enhancing weather and climate forecasts. Here we show that the proposed approach enables real-time predictions that can facilitate decision-making for a wide range of applications in science and engineering.

97 MATHEMATICS AND COMPUTING↗

Development of Composite Photocatalyst Materials that are Highly Selective for Solar Hydrogen Production and their Evaluation in Z-Scheme Reactor Designs

The key technology gap preventing a vertically stacked dual-bed particle suspension reactor from achieving the DOE MYRD&D ultimate cost target for H 2 production remains the lack of materials in particle form factor that exhibit ≥10% solar-to-H 2 energy conversion (STH) efficiency as a suspension. Therefore, our project goals centered around strategies to increase the STH efficiency by enhancing photophysical properties of perovskite oxide particles including increased visible-light absorption, increased selectivity for electrocatalysis of the H 2 evolution reaction (HER) and the O 2 evolution reaction (OER) through development of ultrathin oxide coatings, correlating composition and structure to function, and improving understanding of multiscale transport and kinetic processes.

08 HYDROGEN↗

Energy transfer and scale dynamics in 2D and 3D laser-driven jets

We demonstrate a methodology for diagnosing the multiscale dynamics and energy transfer in complex HED flows with realistic driving and boundary conditions. The approach separates incompressible, compressible, and baropycnal contributions to energy scale-transfer and quantifies the direction of these transfers in (generalized) wavenumber space. We use this to compare the kinetic energy (KE) transfer across scales in simulations of 2D axisymmetric vs fully 3D laser-driven plasma jets. Using the FLASH code, we model a turbulent jet ablated from an aluminum cone target in the configuration outlined. We show that, in addition to its well known bias for underestimating hydrodynamic instability growth, 2D modeling suffers from significant spurious energization of the bulk flow by a turbulent upscale cascade. In 2D, this arises as vorticity and strain from instabilities near the jet's leading edge transfer KE upscale, sustaining a coherent circulation that helps propel the axisymmetric jet farther (⁠ ≈25% by 3.5 ns) and helps keep it collimated. In 3D, the coherent circulation and upscale KE transfer are absent. Here, the methodology presented here may also help with inter-model comparison and validation, including future modeling efforts to alleviate some of the 2D hydrodynamic artifacts highlighted in this study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Lower length scale model for palladium attack of silicon carbide in TRISO fuel

TRistructural ISOtropic (TRISO) particle fuels rely on silicon carbide (SiC) as the primary barrier for metallic fission product (FP) release. Palladium (Pd) generated by fission degrades the SiC layer, resulting in the formation of lamellar layers of palladium silicides (PdxSi) and carbon (C) perpendicular to the direction of attack. The Pd attack has been hypothesized to be responsible for failure of the SiC layer and enhance FP release. To better understand and quantify Pd attack of SiC in TRISO particles, a multiscale, mechanistic model of Pd transport is being developed by the NEAMS program. Previous work provided an initial hypothesis for modeling lamellar microstructure formation in SiC due to Pd attack using a phase-field model. The work described in this report builds on the previous model by using molecular dynamics (MD) simulations to parameterize the phase-field model kinetics, and build a reduced order model in BISON using the improved mesoscale Pd penetration model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Reaction-induced departures from continuum Navier–Stokes turbulence

Reactive hydrodynamic turbulence is an inherently multiscale phenomenon, characterized by the separation between energy-containing, viscous, and molecular length and time scales. The separation between the viscous scale (the Kolmogorov scale) and the molecular mean free path ostensibly justifies a macroscopic description of reactive turbulence via the Navier–Stokes (NS) equations. However, here we use molecular-level simulations to demonstrate that exothermic bimolecular reactions can cause the NS description of turbulence to break down in the near-continuum regime. Sufficiently energetic heat-releasing reactive collisions strongly distort the Maxwell–Boltzmann velocity distribution function, modifying not only the macroscopic chemical rate law but the kinetic-energy-transfer processes as well. This translational nonequilibrium ultimately introduces significant departures from the NS description in the kinetic energy spectra at scales orders of magnitude larger than both the molecular mean free path and the Kolmogorov length scale. These departures prove substantial enough to meaningfully alter integrated quantities, including the overall turbulence kinetic energy itself.

molecular gas dynamics↗

The persistence of large scale structures. Part I. Primordial non-Gaussianity

Abstract We develop an analysis pipeline for characterizing the topology of large scale structure and extracting cosmological constraints based on persistent homology . Persistent homology is a technique from topological data analysis that quantifies the multiscale topology of a data set, in our context unifying the contributions of clusters, filament loops, and cosmic voids to cosmological constraints. We describe how this method captures the imprint of primordial local non-Gaussianity on the late-time distribution of dark matter halos, using a set of N-body simulations as a proxy for real data analysis. For our best single statistic, running the pipeline on several cubic volumes of size 40 (Gpc/h) 3 , we detect f NL loc =10 at 97.5% confidence on ~ 85% of the volumes. Additionally wetest our ability to resolve degeneracies betweenthe topological signature of f NL loc and variation of σ 8 and argue that correctly identifying nonzero f NL loc in this case is possible via an optimal template method. Our method relies on information living at $\mathcal{O}$(10) Mpc/h, a complementary scale with respect to commonly used methods such as the scale-dependent bias in the halo/galaxy power spectrum. Therefore, while still requiring a large volume, our method does not require sampling long-wavelength modes to constrain primordial non-Gaussianity. Moreover, our statistics are interpretable: we are able to reproduce previous results in certain limits and we make new predictions for unexplored observables, such as filament loops formed by dark matter halos in a simulation box.

Astronomy & Astrophysics↗

Direct Imaging of Hydrogen‐Driven Dislocation and Strain Field Evolution in a Stainless Steel Grain

Hydrogen embrittlement (HE) poses a significant challenge to the durability of materials used in hydrogen production and utilization. Disentangling the competing nanoscale mechanisms driving HE often relies on simulations and electron-transparent sample techniques, limiting experimental insights into hydrogen-induced dislocation behavior in bulk materials. This study employs in situ Bragg coherent X-ray diffraction imaging to track three-dimensional (3D) dislocation and strain field evolution during hydrogen charging in a bulk grain of austenitic 316 stainless steel. Tracking a single dislocation reveals hydrogen-enhanced mobility and relaxation, consistent with dislocation dynamics simulations. Subsequent observations reveal dislocation unpinning and climb processes, likely driven by osmotic forces. Additionally, nanoscale strain analysis around the dislocation core directly measures hydrogen-induced elastic shielding. These findings experimentally validate theoretical predictions and offer mechanistic insights into hydrogen-driven dislocation behavior. The quantified nanoscale phenomena serve as critical inputs for multiscale modeling frameworks to predict bulk material responses and accelerate the development of HE-resistant alloys.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗