Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational multiscale”

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

Multiscale plasticity of geomaterials predicted via constrained optimization‐based granular micromechanics

Abstract A general framework to derive nonlinear elastic and elastoplastic material models from granular micromechanics is proposed, where a constraint‐based variational structure is introduced to classical grain contact‐based homogenization methods of hyperelasticity. Like the classical hyperelastic methods, reference solutions for closed‐form hyperelastic material models are analytically derived from the grain‐scale contact mechanics. However, unlike prior methods, the proposed homogenization framework defines closed‐form hyperelastoplastic material models that extend multiscale variational methods to granular plasticity. The proposed framework is used to develop novel granular micromechanics‐based macroscopic models for a Mises type solid, Drucker–Prager type plasticity, and grain‐contact cohesive‐debonding with a deviatorically and volumetrically coupled nonlinearly elastic response. Macroscopic plastic parameters and yield criteria are explicitly related to their microscale counterparts, for example, the friction coefficient governing intergranular slip. Numerical examples and comparison to measurements from the literature, including triaxial compaction of concrete, are provided to investigate model predictions and demonstrate calibration to experimental data.

Bryant, E. C.↗
DOE OSTI2023

Variational quantum simulation of the critical Ising model with symmetry averaging

Here we investigate the use of deep multiscale entanglement renormalization ansatz (DMERA) circuits as a variational ansatz. We use the exactly solvable one-dimensional critical transverse-field Ising model as a test bed. Numerically exact simulation of the quantum circuit ansatz can in this case be carried out to hundreds of qubits by exploiting efficient classical algorithms for simulating matchgate circuits. We find that, for this system, the DMERA strongly outperforms a standard quantum approximate optimization algorithm (QAOA)–style ansatz, and that a major source of systematic error in correlation functions approximated using the DMERA is the breaking of the translational and Kramers-Wannier symmetries of the transverse-field Ising model. We are able to reduce this error by up to four orders of magnitude by symmetry averaging, without incurring additional cost in qubits or circuit depth. Here, we propose that this technique for mitigating systematic error could be applied to noisy intermediate-scale quantum (NISQ) simulations of physical systems with other symmetries.

1-dimensional spin chains↗
DOE OSTI2020

The Strong Impact of Weak Horizontal Convergence on Continental Shallow Convection

Idealized large-eddy simulation (LES) is a basic tool for studying three-dimensional turbulence in the planetary boundary layer. LES is capable of providing benchmark solutions for parameterization development efforts. However, real small-scale atmospheric flows develop in heterogeneous and transient environments with locally varying vertical motions inherent to open multiscale interactive dynamical systems. These variations are often too subtle to detect them by state-of-the-art remote and in situ measurements, and are typically excluded from idealized simulations. The present study addresses the impact of weak [i.e., O(10-6) s-1] short-lived low-level large-scale convergence/divergence perturbations on continental shallow convection. The results show a strong response of shallow nonprecipitating convection to the applied weak large-scale dynamical forcing. Evolutions of CAPE, mean liquid water path, and cloud-top heights are significantly affected by the imposed convergence/divergence. In contrast, evolving cloud-base properties, such as the area coverage and mass flux, are only weakly affected. To contrast those impacts with microphysical sensitivity, the baseline simulations are perturbed assuming different observationally based cloud droplet number concentrations and thus different rainfall. For the tested range of microphysical perturbations, the imposed convergence/divergence provides significantly larger impact than changes in the cloud microphysics. Simulation results presented here provide a stringent test for convection parameterizations, especially important for large-scale models progressing toward resolving some nonhydrostatic effects.

Meteorology & Atmospheric Sciences↗
DOE OSTI2026

Polarons are probes of the dynamic nanoscale environments found in electrochemically doped π-conjugated polymers

Charge carriers (i.e., polarons) in electrochemically doped organic (semi)conductors are proposed to be regulated by several physicochemical features, including the chemical compositions and structures of the π-conjugated frameworks of the semiconductor building blocks, the chemistry of the electrolyte (considering both the salt and solvent), multiscale and time-dependent morphology variations across the material as a function of electrochemical processes, and assorted permutations of these and other factors. To address these hypotheses, we investigate the energetic, optoelectronic, chemical, and local structural properties of negative polarons (radical anions) as a function of electrochemical doping in the donor–acceptor, π-conjugated redox copolymer P(NDI2OD-T2), also referred to as N2200. A critical finding is that there is not just “one type” of polaron in electrochemically doped P(NDI2OD-T2). Rather, an ensemble of polarons exists, with the polarons having optoelectronic signatures that are defined by their nanoscale environments. Importantly, the polaron optical signatures serve as local probes for how the operando electrochemical environments are dynamically working in concert to facilitate charge transport. Collectively, the distinctive and extensive integration of theory and measurement science presented here establishes a baseline for the roles that semiconductor and electrolyte chemistries and dynamic structural features have on polarons in electrochemically doped organic semiconductors and how these factors influence the energetics and rates of polaron transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2022

Thermospheric Neutral Density Variation During the “SpaceX” Storm: Implications From Physics-Based Whole Geospace Modeling

The Starlink satellites launched on 3 February 2022 were lost before they fully arrived in their designated orbits. The loss was attributed to two moderate geomagnetic storms that occurred consecutively on February 3-4. We investigate the thermospheric neutral mass density variation during these storms with the Multiscale Atmosphere-Geospace Environment (MAGE) model, a first-principles, fully coupled geospace model. Simulated neutral density enhancements are validated by Swarm satellite measurements at the altitude of 400-500 km. Comparison with standalone TIEGCM and empirical NRLMSIS 2.0 and DTM-2012 models suggests better performance by MAGE in predicting the maximum density enhancement and resolving the gradual recovery process. Along the Starlink satellite orbit in the middle thermosphere (∼ 200 km altitude), MAGE predicts up to 150% density enhancement near the second storm peak while standalone TIEGCM, NRLMSIS 2.0 and DTM-2012 suggest only ∼ 50% increase. MAGE also suggests altitudinal, longitudinal, and latitudinal variability of storm-time percentage density enhancement due to height dependent Joule heating deposition per unit mass, thermospheric circulation changes, and travelling atmospheric disturbances. This study demonstrates that a moderate storm can cause substantial density enhancement in the middle thermosphere. Thermospheric mass density strongly depends on the strength, timing, and location of high-latitude energy input, which cannot be fully reproduced with empirical models. A physics-based, fully coupled geospace model that can accurately resolve the high-latitude energy input and its variability is critical to modeling the dynamic response of thermospheric neutral density during storm time.

Starlink↗
DOE OSTI2025

Convergence and Quantum Advantage of Trotterized MERA for Strongly-Correlated Systems

Strongly-correlated quantum many-body systems are difficult to study and simulate classically. We recently proposed a variational quantum eigensolver (VQE) based on the multiscale entanglement renormalization ansatz (MERA) with tensors constrained to certain Trotter circuits. Here, we determine the scaling of computation costs for various critical spin chains which substantiates a polynomial quantum advantage in comparison to classical MERA simulations based on exact energy gradients or variational Monte Carlo. Algorithmic phase diagrams suggest an even greater separation for higher-dimensional systems. Hence, the Trotterized MERA VQE is a promising route for the efficient investigation of strongly-correlated quantum many-body systems on quantum computers. Furthermore, we show how the convergence can be substantially improved by building up the MERA layer by layer in the initialization stage and by scanning through the phase diagram during optimization. For the Trotter circuits being composed of single-qubit and two-qubit rotations, it is experimentally advantageous to have small rotation angles. We find that the average angle amplitude can be reduced considerably with negligible effect on the energy accuracy. Benchmark simulations suggest that the structure of the Trotter circuits for the TMERA tensors is not decisive; in particular, brick-wall circuits and parallel random-pair circuits yield very similar energy accuracies.

Miao, Qiang [Duke Quantum Center, Duke University,↗
DOE OSTI2025

Meta Biome: a multiscale model integrating agent-based and metabolic networks to reveal spatial regulation in gut mucosal microbial communities

ABSTRACT Mucosal microbial communities (MMCs) are complex ecosystems near the mucosal layers of the gut essential for maintaining health and modulating disease states. Despite advances in high-throughput omics technologies, current methodologies struggle to capture the dynamic metabolic interactions and spatiotemporal variations within MMCs. In this work, we presentMetaBiome, a multiscale model integrating agent-based modeling (ABM), finite volume methods, and constraint-based models to explore the metabolic interactions within these communities. Integrating ABM allows for the detailed representation of individual microbial agents each governed by rules that dictate cell growth, division, and interactions with their surroundings. Through a layered approach—encompassing microenvironmental conditions, agent information, and metabolic pathways—we simulated different communities to showcase the potential of the model. Using ourin-silicoplatform, we explored the dynamics and spatiotemporal patterns of MMCs in the proximal small intestine and the cecum, simulating the physiological conditions of the two gut regions. Our findings revealed how specific microbes adapt their metabolic processes based on substrate availability and local environmental conditions, shedding light on spatial metabolite regulation and informing targeted therapies for localized gut diseases.MetaBiome provides a detailed representation of microbial agents and their interactions, surpassing the limitations of traditional grid-based systems. This work marks a significant advancement in microbial ecology, as it offers new insights into predicting and analyzing microbial communities. IMPORTANCE Our study presents a novel multiscale model that combines agent-based modeling, finite volume methods, and genome-scale metabolic models to simulate the complex dynamics of mucosal microbial communities in the gut. This integrated approach allows us to capture spatial and temporal variations in microbial interactions and metabolism that are difficult to study experimentally. Key findings from our model include the following: (i) prediction of metabolic cross-feeding and spatial organization in multi-species communities, (ii) insights into how oxygen gradients and nutrient availability shape community composition in different gut regions, and (iii) identification of spatiallyregulated metabolic pathways and enzymes inE. coli. We believe this work represents a significant advance in computational modeling of microbial communities and provides new insights into the spatial regulation of gut microbiome metabolism. The multiscale modeling approach we have developed could be broadly applicable for studying other complex microbial ecosystems.

Microbiology↗
DOE OSTI2024

Multiscale simulation of spatially correlated microstructure via a latent space representation

When deformation gradients act on the scale of the microstructure of a part due to geometry and loading, spatial correlations and finite-size effects in simulation cells cannot be neglected. We propose a multiscale method that accounts for these effects using a variational autoencoder to encode the structure–property map of the stochastic volume elements making up the statistical description of the part. In this paradigm the autoencoder can be used to directly encode the microstructure or, alternatively, its latent space can be sampled to provide likely realizations. Furthermore, we demonstrate the method on three examples using the common additively manufactured material AlSi10Mg in: (a) a comparison with direct numerical simulation of the part microstructure, (b) a push forward of microstructural uncertainty to performance quantities of interest, and (c) a simulation of functional gradation of a part with stochastic microstructure.

Elastoplasticity↗
DOE OSTI2023

Microstructural characterization of nuclear graphite: from the microscale to the nanoscale

Multiple reactor designs use graphite as a moderator of the nuclear reactions and as structural support. During the lifetime of the reactor, multiple aging factors such as neutron irradiation, oxidation, and temperature along with others induce changes in the microstructure and crystal lattice of graphite components. The pore morphology and crystal structure of some phases in graphite can be used to trace the evolution of irradiation defects and mechanical properties of graphite. We present a combination of results from several microscopy techniques to investigate the differences between nuclear graphite grades and the effects of neutron irradiation and oxidation at multiple length scales. This multiscale approach is needed to understand the microstructural variations caused by the raw materials and manufacturing processes as well as how the different phases of graphite are affected by the reactor environment. The results provide insight into the oxidation- and radiation-induced changes of graphite and create a robust baseline of microstructure information that can be used for the selection of materials for the next generation of nuclear power stations. Moreover, the experiments conducted in this work provide an overview of the advantages and limitations of the most common techniques used to characterize nuclear graphite and how these techniques might be applied to study other carbon-based materials used in the nuclear industry.

Arregui Mena, Jose'↗
NASA NTRS2013

The Effect of Fiber Strength Stochastics and Local Fiber Volume Fraction on Multiscale Progressive Failure of Composites

Continuous fiber unidirectional polymer matrix composites (PMCs) can exhibit significant local variations in fiber volume fraction as a result of processing conditions that can lead to further local differences in material properties and failure behavior. In this work, the coupled effects of both local variations in fiber volume fraction and the empirically-based statistical distribution of fiber strengths on the predicted longitudinal modulus and local tensile strength of a unidirectional AS4 carbon fiber/ Hercules 3502 epoxy composite were investigated using the special purpose NASA Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC); local effective composite properties were obtained by homogenizing the material behavior over repeating units cells (RUCs). The predicted effective longitudinal modulus was relatively insensitive to small (~8%) variations in local fiber volume fraction. The composite tensile strength, however, was highly dependent on the local distribution in fiber strengths. The RUC-averaged constitutive response can be used to characterize lower length scale material behavior within a multiscale analysis framework that couples the NASA code FEAMAC and the ABAQUS finite element solver. Such an approach can be effectively used to analyze the progressive failure of PMC structures whose failure initiates at the RUC level. Consideration of the effect of local variations in constituent properties and morphologies on progressive failure of PMCs is a central aspect of the application of Integrated Computational Materials Engineering (ICME) principles for composite materials.

multiscale analysis↗
DOE OSTI2025

Autonomous Nanoparticle Synthesis Guided by In Situ Multiscale Structural Characterization

Autonomous synthesis platforms promise rapid exploration of vast parameter spaces; yet, integrating in situ structural characterization in closed-loop synthesis optimization remains challenging. We demonstrate a realization of such a closed-loop platform coupled with a droplet-flow microreactor, in situ X-ray scattering methods (SAXS/WAXS), and Gaussian process optimization to synthesize citrate-reduced Au nanoparticles with targeted characteristics. The system efficiently explored ∼19,000 synthesis recipes through 365 experiments, achieving precise control over size (4–60 nm) and polydispersity (σ < 0.11) across large citrate/gold ratios, exceeding traditional synthesis boundaries (1–10). Beyond confirming classical Turkevich–Frens trends, partial-dependence analysis revealed strong nonlinear coupling among precursor, citrate, and pH effects. Combining quantitative SAXS/WAXS analysis with electron microscopy characterization, we uncovered that crystallite size (d c ) and particle size (d) follow d c = 0.18d + β, where synthesis chemistry controls the intercept β while maintaining a universal slope. This parallel-band structure enables independent tuning of crystallite domain size at fixed particle diameter through a combination of chloride, gold precursor, citrate, and pH contributions (cross-validated Spearman ρ = 0.7 ± 0.1). High-resolution electron microscopy shows multiple lattice-fringe orientations within single particles, directly confirming polycrystalline domains and the ability to tune d c at the fixed d. The platform’s validation includes indistinguishable static versus flowing measurements, stable droplet transport at 100 °C, and <5% run-to-run variation, establishing a robust framework for mapping and controlling multiscale nanoparticle structure across expansive chemical spaces. In conclusion, the developed closed-loop platform can be applied to a borad range of nanosyntheis processes.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Variations in Creep Performance of P92 Steel Within Its Composition Range: Influence of N Solubility

Underspecified composition ranges often lead to alloys with unpredictable mechanical performance. To better understand the changes in microstructure and mechanical performance associated with variations of key elements, three versions of P92 are formulated within, or close to, the specified allowable for N, B, and C ranges. Chromium and Si are also varied to influence N solubility. Different service conditions (i.e., temperature and stress) are explored. It is observed that >80% decrease in creep life occurs at 625 °C and 155 MPa for the highest B and N containing alloy. Multiscale characterization reveals key changes due to the trace element variation. The high B and N containing alloy forms deleterious BN precipitates with morphology that promotes crack nucleation and damage accumulation, but this alloy additionally forms higher fractions of beneficial MX precipitates. The alloy with the lowest B and N concentrations but greater C content shows the best creep performance—a consequence of the refined M 23 C 6 carbide precipitate population and the absence of large-scale inclusions or BN precipitates. Calculations of creep activation energy reveal that the high B and N containing alloy is more prone to damage accumulation which causes an early onset of accelerating creep and greater minimum creep rate.

36 MATERIALS SCIENCE↗
NASA NTRS2017

Solar Activity Across the Scales: From Small-Scale Quiet-Sun Dynamics to Magnetic Activity Cycles

Observations as well as numerical and theoretical models show that solar dynamics is characterized by complicated interactions and energy exchanges among different temporal and spatial scales. It reveals magnetic self-organization processes from the smallest scale magnetized vortex tubes to the global activity variation known as the solar cycle. To understand these multiscale processes and their relationships, we use a two-fold approach: 1) realistic 3D radiative MHD simulations of local dynamics together with high resolution observations by IRIS, Hinode, and SDO; and 2) modeling of solar activity cycles by using simplified MHD dynamo models and mathematical data assimilation techniques. We present recent results of this approach, including the interpretation of observational results from NASA heliophysics missions and predictive capabilities. In particular, we discuss the links between small-scale dynamo processes in the convection zone and atmospheric dynamics, as well as an early prediction of Solar Cycle 25.

Solar↗
NASA NTRS2017

Solar Activity Across the Scales: from Small-Scale Quiet-Sun Dynamics to Magnetic Activity Cycles

Observations as well as numerical and theoretical models show that solar dynamics is characterized by complicated interactions and energy exchanges among different temporal and spatial scales. It reveals magnetic self-­‐organization processes from the smallest scale magnetized vortex tubes to the global activity variation known as the solar cycle. To understand these multiscale processes and their relationships, we use a two-­fold approach: 1) рealistic 3D radiative MHD simulations of local dynamics together with high-­‐resolution observations by IRIS, Hinode, and SDO; and 2) modeling of solar activity cycles by using simplified MHD dynamo models and mathematical data assimilation techniques. We present recent results of this approach, including the interpretation of observational results from NASA heliophysics missions and predictive capabilities. In particular, we discuss the links between small-­‐scale dynamo processes in the convection zone and atmospheric dynamics, as well as an early prediction of Solar Cycle 25.

Magneti↗
DOE OSTI2022

Exploring the Impact of N Solubility and Trace Elements on the Creep Properties of P92 Steel

Trace elements can have major consequences on the properties of creep resistant martensitic steels. To better understand the changes in microstructure and mechanical performance associated with the variation of trace elements and enable design of steels operating at more demanding conditions, we formulated three versions of P92 within the specified allowable N, B, C and Si ranges. Different service conditions were explored, and >80% decrease in creep life was observed at 625°C 155MPa for the highest B and N containing alloy. Multiscale characterization (OM, SEM, TEM) revealed key changes due to the trace element variation. The high trace content alloy formed deleterious BN precipitates with morphology that would promote crack nucleation, but also formed higher fractions of beneficial MX precipitates. The alloy with the lowest trace content showed the best creep performance – a result of the refined precipitate populations and the absence of large-scale inclusions and BN precipitates.

Antonov, Stoichko↗
DOE OSTI2023

Modeling and Simulation of Advanced Manufacturing Techniques using MOOSE and MALAMUTE

Advanced manufacturing techniques offer increased geometry complexity, energy and material usage efficiency improvements, and an expanded palette of materials as compared to conventional manufacturing approaches. Advanced-manufacturing-produced parts can experience wide variations in the final microstructure, and these microstructure variations significantly impact the parts’ performance. In this chapter, we present recent code developments within Multiphysics Object-Oriented Simulation Environment (MOOSE) and in the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE). Here we demonstrate applying these modeling and simulation codes to two advanced manufacturing process types: advanced sintering techniques and laser-based additive manufacturing techniques. The multiphysics and multiscale capabilities of these codes enable prediction of the microstructure evolution resulting from variations in the Advanced manufacturing process parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Irradiation damage concurrent challenges with RAFM and ODS steels for fusion reactor first-wall/blanket: a review

Reduced activation ferritic martensitic (RAFM) and oxide dispersion strengthened (ODS) steels are the most promising candidates for fusion first-wall/blanket (FW/B) structures. The performance of these steels will deteriorate during service due to neutron damage and transmutation-induced gases, such as helium/hydrogen, at elevated operating temperatures. Here, after highlighting the operating conditions of fusion reactor concepts and a brief overview, the main irradiation-induced degradation challenges associated with RAFM/ODS steels are discussed. Their long-term degradation scenarios such as (a) low-temperature hardening embrittlement (LTHE)—including dose-temperature dependent yield stress, tensile elongations, necking ductility, test temperature effect on hardening, Charpy impact ductile-to-brittle transition temperature and fracture toughness, (b) intermediate temperature cavity swelling, (c) the effect of helium on LTHE and cavity swelling, (d) irradiation creep and (e) tritium management issues are reviewed. The potential causes of LTHE are discussed, which highlights the need for advanced characterisation techniques. The mechanical properties, including the tensile/Charpy impact of RAFM and ODS steels, are compared to show that the current generation of ODS steels also suffers from LTHE, and shows irradiation hardening up to high temperatures of ~400 °C–500 °C. To minimise this, future ODS steel development for FW/B-specific application should target materials with a lower Cr concentration (to minimise α'), and minimise other elements that could form embrittling phases under irradiation. RAFM steel-designing activities targeting improvements in creep and LTHE are reviewed. The need to better understand the synergistic effects of helium on the thermo-mechanical properties in the entire temperature range of FW/B is highlighted. Because fusion operating conditions will be complex, including stresses due to the magnetic field, primary loads like coolant pressure, secondary loads from thermal gradients, and due to spatial variation in damage levels and gas production rates, an experimentally validated multiscale modelling approach is suggested as a pathway to future reactor component designing such as for the fusion neutron science facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Multimodal and multiscale strengthening mechanisms in Al-Ni-Zr-Ti-Mn alloy processed by laser powder bed fusion additive manufacturing

The unique thermokinetics of laser-powder bed fusion additive manufacturing (L-PBFAM) has been exploited for development of a novel high-strength Al-Ni-Ti-Zr-Mn alloy. The addition of 0.5 wt% Mn leads to extraordinary improvement in ultimate tensile strength (502 MPa) and work hardening due to the activation of two Mn-induced strengthening mechanisms. First, by a bimodal particle strengthening effect due to Al 31 Mn 6 Ni 12 nano-quasi-crystal particles rejected in inter-dendritic spaces and fibrous Al 3 Ni eutectic dendritic channels, which predominately contributes to the strength improvement, and second by solid solution strengthening from remaining Mn entrapped in Al. These mechanisms supplement the particle strengthening effect imparted by coherent and incoherent Al 3 (Ti,Zr) co-precipitates present at melt pool boundaries, dislocation strengthening due to solidification induced strain, and Hall-Petch and backstress strengthening effect due to heterogenous grain size distribution occurring at various length scales. In conclusion, the solidification dynamics and hierarchical heat distribution that are associated with L-PBFAM resulted in complex spatial variations in these strengthening phenomena and were investigated via a high-throughput multiscale structure–property correlation that involved thermokinetic simulation, X-ray diffraction, high-resolution nanoindentation mapping, and site-specific transmission electron microscopy of the alloy.

36 MATERIALS SCIENCE↗