Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale simulation”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Impact of Tropical Cyclone Wind Forcing on the Global Climate in a Fully Coupled Climate Model

Tropical cyclones (TCs) alter upper-ocean temperature and influence ocean heat content via enhanced turbulent mixing. A better understanding of the role of TCs within the climate system requires a fully coupled modeling framework, where TC-induced ocean responses feed back to the atmosphere and subsequently to the climate mean state and variability. Here, we investigate the impacts of TC wind forcing on the global ocean and the associated feedbacks within the climate system using the fully coupled Community Earth System Model version 1.3 (CESM1.3). Using the low-resolution version of CESM1.3 (1° atmosphere and ocean grid spacing) with no intrinsic TCs, we conduct a suite of sensitivity experiments by inserting TC winds extracted from a high-resolution (0.25° atmosphere grid spacing) TC-permitting simulation into the low-resolution model. Results from the low-resolution TC experiment are compared to a low-resolution control simulation to diagnose TCs’ impact. We found that the added TC winds can increase ocean heat content by affecting ocean vertical mixing, air–sea enthalpy fluxes, and cloud amount. Here, the added TCs can influence mean SST, precipitation, ocean subsurface temperature, and ocean mixed layer depth. We found a strengthening of the wind-driven subtropical cells and a weakening of the Atlantic meridional overturning circulation due to the changes of surface buoyancy fluxes. TCs in the model cause anomalous equatorward ocean heat convergence in the deep tropics and an increase of poleward ocean heat transport out of the subtropics. Our modeling results provide new insights into the multiscale interactions between TCs and the coupled climate system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Micro-Mechanically Guided High-Throughput Alloy Design Exploration Towards Metastability-Induced H Embrittlement Resistance

We develop a high-throughput approach for studying H embrittlement (HE)-resistance in alloys, which is based on combinatorial compositional screening of metastability effects by in situ scanning electron microscopy H-analyses. The project objective included: (i) technique development of high-throughput screening (HTS) for HE-resistance, (ii) discovery of new metallic materials with superior HE-resistance, (iii) Multiscale verification of HE-resistance of the new alloys and H-barrier layers, from atomic scale to an engineering scale. The investigation focused on a model system that enables exploration of different metastable states in Fe-based complex-concentrated alloys. Composition spread islands with hundreds of varying compositions are fabricated on a single substrate using a combinatorial co-sputtering technique. To screen these alloys, we developed and employed a home-built SEM-based integrated analysis system capable of characterizing H-permeability, H-trapping, H-influence on mechanical properties, and other material properties, as well as atomistic to continuum simulations to study the underlying physics.

08 HYDROGEN↗
NASA NTRS2008

New Developments in the Embedded Statistical Coupling Method: Atomistic/Continuum Crack Propagation

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain has been enhanced. The concurrent MD-FEM coupling methodology uses statistical averaging of the deformation of the atomistic MD domain to provide interface displacement boundary conditions to the surrounding continuum FEM region, which, in turn, generates interface reaction forces that are applied as piecewise constant traction boundary conditions to the MD domain. The enhancement is based on the addition of molecular dynamics-based cohesive zone model (CZM) elements near the MD-FEM interface. The CZM elements are a continuum interpretation of the traction-displacement relationships taken from MD simulations using Cohesive Zone Volume Elements (CZVE). The addition of CZM elements to the concurrent MD-FEM analysis provides a consistent set of atomistically-based cohesive properties within the finite element region near the growing crack. Another set of CZVEs are then used to extract revised CZM relationships from the enhanced embedded statistical coupling method (ESCM) simulation of an edge crack under uniaxial loading.

Saether, E.↗
DOE OSTI2023

Multiscale Electricity Modeling for Evaluating Carbon Capture and Sequestration Technologies (Final Report)

Carbon capture and sequestration (CCS) technologies that can operate with a high degree of operating flexibility could provide necessary electric grid flexibility in a system with high shares of variable renewables. This project examines the deployment and dispatch potential of twelve unique flexible CCS (FLECCS) technologies that encompass post-combustion carbon dioxide (CO 2 ) capture designs, concepts using a storage media to enable energy arbitrage, and hybrid processes that integrate CCS with direct air capture (DAC) for flexibility with net zero or negative CO 2 emissions. FLECCS technology potential is explored with a multi-model, multi-scale framework including the Regional Energy Deployment System (ReEDS) electric sector capacity expansion model (CEM) and the PLEXOS production cost model (PCM). Innovative methods were developed to represent FLECCS technology operating modes, performance, and cost in the two models. ReEDS was then used to simulate nine scenarios for each FLECCS technology, three CO 2 emissions price futures reaching $\$$150, $\$$225, and $\$$300/tCO 2 in 2050; and three scenarios for FLECCS technology deployment favorability relative to competing technologies. For each CO 2 price and reference FLECCS favorability, the 2050 infrastructures from ReEDS model are downscaled and implemented in PLEXOS to examine hourly dispatch under detailed operational constraints that are not included in ReEDS. FLECCS technologies exhibited a wide range of deployment potential ranging from none to several hundred gigawatts of capacity, with outcomes highly sensitive to input cost and performance parameters that are inherently highly uncertain. When deployed, FLECCS tended to displace a combination of wind, solar, and natural gas-based technologies rather than supporting increased renewable deployment. As a result, CO 2 emissions reductions facilitated by FLECCS deployment tended to come with higher overall system costs and electricity prices. When economically competitive, FLECCS technologies can contribute significant flexible generation and firm capacity to the grid, but continued technology development and an expanded analytical scope are necessary to fully understand FLECCS deployment potential its impact on the electric power sector. Follow-on analysis incorporating captured CO 2 tax credit value from the Inflation Reduction Act (IRA) and other potential policy scenarios could be particularly valuable, as this policy can substantially change the relative competitiveness of FLECCS technologies.

03 NATURAL GAS↗
DOE OSTI2025

Designing Physicochemically‐Ordered Interphases for High‐Performance Composites

To enhance the mechanical properties of carbon fiber‐reinforced polymer composites, a physicochemical scaffold is designed incorporating microscopically architected chemically reactive nanofibers that act as a multiscale bridge between the carbon fibers and the matrix. Thermally activated nanofibers leverage their morphologically driven mechanochemical properties to form covalent bonds with adjacent polymer molecules, creating a co‐continuous network that dramatically enhances fiber‐matrix load transfer. By meticulously controlling the nanofiber architecture through variable surface area, functional group availability, and polymer chain alignment effects, the extent of covalent bonding between nanofibers and the matrix is manipulated ultimately resulting in improved carbon fiber‐matrix adhesion. Further, the concept was validated using polyacrylonitrile nanofibers within an acrylonitrile butadiene styrene matrix in a discontinuous carbon fiber‐reinforced composite system. Nanomechanical studies using atomic force microscopy and low‐field nuclear magnetic resonance spectroscopy confirmed immobilized, chemically transferred, and ordered nanostructures at the interphase. The resulting composites demonstrated ≈56% and ≈175% improvements in tensile strength and toughness, respectively, compared to composites without nanofiber. Comprehensive thermal, rheological, and X‐ray scattering analyses, along side all‐atomic molecular dynamics simulations, revealed the fundamental mechanisms behind these improvements in mechanical behavior. The versatility and efficacy of the approach have the potential to address longstanding interphase challenges in the composite industry.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Denoising diffusion algorithm for inverse design of microstructures with fine-tuned nonlinear material properties

Here we introduce a denoising diffusion algorithm to discover microstructures with nonlinear fine-tuned properties. Denoising diffusion probabilistic models are generative models that use diffusion-based dynamics to gradually denoise images and generate realistic synthetic samples. By learning the reverse of a Markov diffusion process, we design an artificial intelligence to efficiently manipulate the topology of microstructures to generate a massive number of prototypes that exhibit constitutive responses sufficiently close to designated nonlinear constitutive behaviors. To identify the subset of microcstructures with sufficiently precise fine-tuned properties, a convolutional neural network surrogate is trained to replace high-fidelity finite element simulations to filter out prototypes outside the admissible range. Results of this study indicate that the denoising diffusion process is capable of creating microstructures of fine-tuned nonlinear material properties within the latent space of the training data. More importantly, this denoising diffusion algorithm can be easily extended to incorporate additional topological and geometric modifications by introducing high-dimensional structures embedded in the latent space. Numerical experiments are conducted on the open-source mechanical MNIST data set (Lejeune, 2020). Consequently, this algorithm is not only capable of performing inverse design of nonlinear effective media, but also learns the nonlinear structure–property map to quantitatively understand the multiscale interplay among the geometry, topology, and their effective macroscopic properties.

42 ENGINEERING↗
DOE OSTI2023

Unveiling the catalytic mechanism of GTP hydrolysis in microtubules

Microtubules (MTs) are large cytoskeletal polymers, composed of αβ-tubulin heterodimers, capable of stochastically converting from polymerizing to depolymerizing states and vice versa. Depolymerization is coupled with hydrolysis of guanosine triphosphate (GTP) within β-tubulin. Hydrolysis is favored in the MT lattice compared to a free heterodimer with an experimentally observed rate increase of 500- to 700-fold, corresponding to an energetic barrier lowering of 3.8 to 4.0 kcal/mol. Mutagenesis studies have implicated α-tubulin residues, α:E254 and α:D251, as catalytic residues completing the β-tubulin active site of the lower heterodimer in the MT lattice. The mechanism for GTP hydrolysis in the free heterodimer, however, is not understood. Additionally, there has been debate concerning whether the GTP-state lattice is expanded or compacted relative to the GDP state and whether a “compacted” GDP-state lattice is required for hydrolysis. In this work, extensive quantum mechanics/molecular mechanics simulations with transition-tempered metadynamics free-energy sampling of compacted and expanded interdimer complexes, as well as a free heterodimer, have been carried out to provide clear insight into the GTP hydrolysis mechanism. α:E254 was found to be the catalytic residue in a compacted lattice, while in the expanded lattice, disruption of a key salt bridge interaction renders α:E254 less effective. The simulations reveal a barrier decrease of 3.8 ± 0.5 kcal/mol for the compacted lattice compared to a free heterodimer, in good agreement with experimental kinetic measurements. Additionally, the expanded lattice barrier was found to be 6.3 ± 0.5 kcal/mol higher than compacted, demonstrating that GTP hydrolysis is variable with lattice state and slower at the MT tip.

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

Machine-learning modeling of magnetization dynamics in quasi-equilibrium and driven metallic spin systems

Here, we present a perspective on recent progress in machine-learning (ML) force-field approaches for large-scale Landau–Lifshitz–Gilbert (LLG) simulations of metallic spin systems. Building on a generalization of the Behler–Parrinello (BP) architecture originally developed for quantum molecular dynamics, we develop scalable and transferable ML models that faithfully capture the complex, environment-dependent electron-mediated exchange fields characteristic of itinerant magnets. A central ingredient of this framework is the implementation of symmetry-aware magnetic descriptors based on group-theoretical bispectrum formalisms. Leveraging these ML force fields, LLG simulations faithfully reproduce hallmark non-collinear magnetic orders—such as the 120° and tetrahedral states—on the triangular lattice, and successfully capture the complex spin textures emerging in the mixed-phase states of a square-lattice double-exchange model under thermal quench. We further discuss a generalized potential theory that extends the BP formalism to incorporate both conservative and nonconservative electronic torques, thereby enabling ML models to learn nonequilibrium exchange fields from computationally demanding microscopic approaches such as nonequilibrium Green’s-function techniques. This extension yields quantitatively accurate predictions of voltage-driven domain-wall motion and establishes a foundation for quantum-accurate, multiscale modeling of nonequilibrium spin dynamics and spintronic functionalities.

Descriptors↗
NASA NTRSDate not supplied

Hierarchical Multiscale Process Modeling of A Textile Composite Y-Joint for the Aurora D8 Aircraft

Woven polymer matrix composites (PMCs) are widely used in aerospace applications for their favorable specific properties. However, the behavior of the inter- and intra-tow matrix during manufacturing of woven PMCs is difficult to isolate experimentally due to geometric and processing complexities. Classical closed-form homogenization approaches allow for composite thermos-chemo-mechanical properties to be characterized at any fiber volume fraction. In this work, virtual curing is applied to the inter- and intra-tow matrix materials of a plain weave repeating unit cell (RUC). Thermo-chemo-mechanical property evolution within the tows is calculated using the Concentric Cylinder Model (CCM) and the Rule of Mixtures (ROM). Curing is implemented through user-written subroutines within the commercial finite element method software Abaqus. Virtual curing results are compared with those from a linear-elastic thermal cool down simulation in Abaqus. Results indicate that tow property evolution during curing has a non-negligible effect on the final stress state.

curing↗
DOE OSTI2025

Influences of shock imprinting on mix in a 3D-printed porous media

Mixing of materials in porous media can cause a significant impact on fusion yield as previously demonstrated by the National Ignition Facility (NIF) MARBLE Campaign. Initially, the reactants are separated, with deuterium in the lattice struts and a tritium gas fill in the voids. Lattice parameters such as the strut thickness and relative pitch, provide a control for the mix parameters in the experiment. Los Alamos National Laboratory’s (LANL) BOSQUE project looks to better understand how the mix of the reactants and shell materials impact the fusion burn and resultant yield on various laser platforms. xRAGE’s Eulerian hydrodynamics and adaptive mesh refinement (AMR) provide the unique ability to study the impacts of multiscale features of complex lattice structures. This modeling provides the ability to measure shock front variations as the wave progress’ through a given media. Initial conditions of the lattice are essential to accurately model mix and burn measured by experiment. By varying the initial orientation and densities of these lattice regions the early time dynamics of how the shock is launched into the system is changed and advocates for the study of resulting effects. In this work, we will study the sensitivities of shock effects in varying 3D printed geometric systems and how these shocks alter the structure and mix in the lattice. We will discuss both preliminary experimental results and simulations to help plan and constrain future experiments where we will study the impact of different lattice geometries and lattice bulk densities. This work concludes with the relative impacts of lattice geometries on shock speeds at different bulk densities and the resultant mix widths due to those shock interactions. We see agreement with theory at the higher end of our bulk density study, however, as we approach lower bulk density systems the dynamics of these interactions begin to change.

2 Photon Polymerization↗
DOE OSTI2024

Driving macro-scale transformations in three-dimensional-printed biopolymers through controlled induction of molecular anisotropy at the nanoscale

Motivated by the need to harness the properties of renewable and biodegradable polymers for the design and manufacturing of multi-scale structures with complex geometries, we have employed our additive manufacturing platform that leverages molecular self-assembly for the production of metre-scale structures characterized by complex geometries and heterogeneous material composition. As a precursor material, we used chitosan, a chemically modified form of chitin, an abundant and sustainable structural polysaccharide. We demonstrate the ability to control concentration-dependent crystallization as well as the induction of the preferred orientation of the polymer chains through the combination of extrusion-based robotic fabrication and directional toolpathing. Anisotropy is demonstrated and assessed through high-resolution micro-X-ray diffraction in conjunction with finite element simulations. Using this approach, we can leverage controlled and user-defined small-scale propagation of residual stresses to induce large-scale folding of the resulting structures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Determining the Length Scale of Transport Impedances in Li-Ion Electrodes: Li(Ni 0.33 Mn 0.33 Co 0.33 )O 2

Battery cathodes are complex multiscale, multifunctional materials. The length scale at which the dominant impedance arises may be difficult to determine even with the most advanced experimental characterization efforts, and thus modeling can play an important role in analysis. Discharge and voltage relaxation curves, interrogated with theory, are used to distinguish between transport impedance that arise on the scale of the active crystal and on the scale of agglomerates (secondary particles) comprised of nanoscale crystals. Model-selection algorithms are applied to determine that the agglomerate scale is dominant in the Li Ni 0.33 Mn 0.33 Co 0.33 O 2 electrode studied here. Furthermore, conditions where the agglomerate and crystal-scale models yield distinct simulation results are demonstrated, providing approaches that can be applied to other systems.

25 ENERGY STORAGE↗
NASA NTRSDate not supplied

Porosity Model Parameter Effects on Metal Additive Manufacturing Simulations

This work investigates the effects of model parameters on thermomechanical simulations of a metal additive manufacturing process. The first half of this work explores the sensitivity and effects of the two porosity parameters used in the evolving constitutive model of the simulation. Varying the two parameters of the model show that for the simulations conducted here, the linear coefficient drives the predicted residual stress results while the constant term has little impact. The second half of the work considers the size, order, and the number of element layers in the finite element mesh. The stress distribution plots show that second-order elements along with a dense mesh marginally improve mesh convergence at a significant cost.

Additive Manufacturing↗
NASA NTRSDate not supplied

NavCube3-mini Lunar GNSS Receiver

This paper describes development and testing of NASA Goddard Space Flight Center’s new NavCube3-mini (NC3m) spaceborne, weak-signal GNSS receiver, which targets all Earth orbit regimes with special focus on lunar applications. NC3m derives from the ground-breaking Magnetospheric Multiscale (MMS) mission Navigator GPS receiver. The MMS-Navigator (launched 2015) holds a Guinness World Record for highest altitude GPS fix and is currently in a highly elliptic orbit with a 29 Earth radii apogee, nearly half lunar distance. NC3m has reduced size, weight, and power compared to the MMS-Navigator, making it suitable for smallsat applications, and adds multi-frequency and multi-GNSS capabilities, among other improvements. A NC3m engineering test unit was subjected to and successfully completed a comprehensive Technology Readiness Level 6 (system/subsystem model or prototype demonstration in a relevant environment) testing campaign in the second half of 2022. This involved high-fidelity simulations in LEO, GEO, and lunar regimes and full environmental testing, including vibration, thermal vacuum, and electromagnetic compatibility. An overview of the NC3m receiver, the test setup, and a sample of results is presented. The results include predicted performance at high-altitude and in the lunar regime.

NavCube3↗
DOE OSTI2023

The Role of Topography in Controlling Evapotranspiration Age

Abstract Evapotranspiration (ET) age is a key metric of water sustainability but a major unknown partly due to the extreme difficulty in modeling it. Groundwater is found to be important in ET age variations in small‐scale studies, yet our understanding is insufficient because groundwater systems are nested across scales. Here, we conducted GPU‐accelerated particle tracking with integrated hydrologic modeling to quantify the variations in ET age at a regional scale of ∼0.4 M km 2 . Simulation results reveal topography‐driven flow paths shaping the spatial and temporal patterns of ET age variations. On ridges, where root zone decoupling with deep subsurface storage, ET age is generally young, with seasonal variations dominated by meteorological conditions. In the valley bottom, ET age is generally old, with significant subseasonal variations caused by the convergence of subsurface flow paths. On hillslopes with water table depths ranging from 1 to 10 m, ET age shows strong seasonal variations caused by the connections with lateral groundwater regulated by ET demand. Our modeling approach provides insights into the basic linkages between ET age and topography at large scale. Our work highlights the perspective of multiscale studies of ET age, suggesting new field experiments to test these process connections and to determine if such linkages warrant inclusion in Earth System Models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A novel conditional formulation of the Vlasov–Ampère equations: a conservative, positivity, asymptotic and Gauss law preserving scheme

We propose a novel reformulation of the Vlasov–Ampère equations for plasmas that reveals discrete symmetries that enables simultaneous conservation of mass, momentum and energy; preservation of Gauss’s law; positivity of the distribution function; and consistency with quasi-neutral asymptotics. The approach employs variable and coordinate transformations to yield a coupled system comprising a modified Vlasov equation and associated moment–field equations. The modified Vlasov equation advances a conditional distribution function that excludes mass, momentum and energy densities, which are instead evolved through moment equations enforcing the relevant symmetries, conservation laws and involution constraints. This reformulation aligns naturally with a recent slow-manifold reduction technique, which separates fast electron time scales and simplifies the treatment of the quasi-neutral limit within the reduced moment–field subsystem. Using this framework, we develop a numerical method for the reduced 1D1V subsystem that, for the first time in the literature, satisfies all key physical constraints while maintaining a quasi-neutral asymptotic behaviour. The advantages of the method are demonstrated on canonical electrostatic test problems, including the multiscale ion acoustic shock wave.

1D1V↗
NASA NTRS2010

Multiscale Modeling of Solar Coronal Magnetic Reconnection

Magnetic reconnection is widely believed to be the primary process by which the magnetic field releases energy to plasma in the Sun's corona. For example, in the breakout model for the initiation of coronal mass ejections/eruptive flares, reconnection is responsible for the catastrophic destabilizing of magnetic force balance in the corona, leading to explosive energy release. A critical requirement for the reconnection is that it have a "switch-on' nature in that the reconnection stays off until a large store of magnetic free energy has built up, and then it turn on abruptly and stay on until most of this free energy has been released. We discuss the implications of this requirement for reconnection in the context of the breakout model for CMEs/flares. We argue that it imposes stringent constraints on the properties of the flux breaking mechanism, which is expected to operate in the corona on kinetic scales. We present numerical simulations demonstrating how the reconnection and the eruption depend on the effective resistivity, i.e., the effective Lundquist number, and propose a model for incorporating kinetic flux-breaking mechanisms into MHO calculation of CMEs/flares.

Antiochos, Spiro K.↗
NASA NTRS2018

Effects of Multiscale Phase-Mixing and Interior Conductance in the Lunar-like Pickup Ion Plasma Wake. First Results from 3-D Hybrid Kinetic Modeling

The study of multiscale pickup ion phase-mixing in the lunar plasma wake with a hybrid model is the main subject of our investigation in this paper. Photoionization and charge exchange of protons with the lunar exosphere are the ionization processes included in our model. The computational model includes the self-consistent dynamics of the light (hydrogen ions or hydrogen molecule ion or helium ion), and heavy (sodium ion ) pickup ions. The electrons are considered as a fluid. The lunar interior is considered as a weakly conducting body. In this paper we considered for the first time the cumulative effect of heavy neutrals in the lunar exosphere (e.g., Aluminum, Argon), an effect which was simulated with one species of sodium ion but with a tenfold increase in total production rates. We find that various species produce various types of plasma tail in the lunar plasma wake. Specifically, sodium ion and helium ion pickup ions form a cycloid-like tail, whereas the hydrogen ion or hydrogen molecule ion pickup ions form a tail with a high density core and saw-like periodic structures in the flank region. The length of these structures varies from 1:5 R(sub M) to 3:3 R(sub M) depending on the value of gyro radius for hydrogen ion or hydrogen molecule ion pickup ions. The light pickup ions produce more symmetrical jump in the density and magnetic field at the Mach cone which is mainly controlled by the conductivity of the interior, an effect previously unappreciated. Although other pickup ion species had little effect on the nature of the interaction of the Moon with the solar wind, the global structure of the lunar tail in these simulations appeared quite different when the hydrogen molecule ion production rate was high.

Lipatov, A.S.↗