Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiscale systems”

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 307 records · Page 17

Multiscale Crystalline Structure of Confined Polypeptoid Films: The Effect of Alkyl Side Chain Branching

Here we report the effect of alkyl side chain branching on melt-recrystallization of nanoconfined polypeptoid films using poly(N-octyl glycine) (PNOG) and poly(N-2-ethyl-1-hexyl glycine) (PNEHG) as model systems. Upon cooling from the isotropic melt, confined PNOG molecules recrystallize into a near-perfect orthorhombic crystal structure with the board-like molecules stacked face-to-face in the substrate-parallel direction, resulting in long-range ordered wormlike lamellae that occupy the entire film. By contrast, rod-like PNEHG molecules bearing branched N-2-ethyl-1-hexyl side chains stack into a columnar hexagonal mesophase with their backbones oriented parallel to the substrates, forming micron-sized sheaf-like superstructures under confinement, exposing large areas of empty spaces in the film. These findings highlight the effect of alkyl side chain branching on the packing motif and multiscale crystalline structure of polypeptoids under a nanoconfined geometry.

36 MATERIALS SCIENCE↗

Permeation of CO 2 and N 2 through glassy poly(dimethyl phenylene) oxide under steady- and presteady-state conditions

Glassy polymers are often used for gas separations because of their high selectivity. Although the dual-mode permeation model correctly fits their sorption and permeation isotherms, its physical interpretation is disputed, and it does not describe permeation far from steady state, a condition expected when separations involve intermittent renewable energy sources. To develop a more comprehensive permeation model, we combine experiment, molecular dynamics, and multiscale reaction–diffusion modeling to characterize the time-dependent permeation of N 2 and CO 2 through a glassy poly(dimethyl phenylene oxide) membrane, a model system. Simulations of experimental time-dependent permeation data for both gases in the presteady-state and steady-state regimes show that both single- and dual-mode reaction–diffusion models reproduce the experimental observations, and that sorbed gas concentrations lag the external pressure rise. The results point to environment-sensitive diffusion coefficients as a vital characteristic of transport in glassy polymers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Numerical evaluation of photosensitive tracers as a strategy for separating surface and subsurface transient storage in streams

For this work, we numerically evaluated photosensitive tracers as a potential strategy for separating the effects of surface and hyporheic storage zones (SSZs and HSZs, respectively) on stream corridor transport. Correctly separating HSZ and SSZ effects is critical to estimating the hydro-biogeochemical function of a stream because HSZs and SSZs expose solutes to significantly different biogeochemical conditions, like sunlight exposure, microbial processes, and oxygen availability. Our numerical experiments used a multiscale river-corridor transport model implemented in the ATS code, which accommodates multiple storage zones with distinct travel time distributions and biogeochemical reactions. For parameter inferences, we used Bayesian inverse modeling. We found that breakthrough curves for photo-decaying tracers from day and night injection can delineate surface and hyporheic transient storage contributions, but only when interpreted jointly through a two-storage zone model. Numerical experiments that used only daytime injection or interpreted breakthrough curves with a single storage zone model yielded good fit to breakthrough curves, but parameter estimates were biased and controlling processes misattributed, examples of good model fits for the wrong reasons. Using those biased parameter estimates in reactive transport simulations resulted in significantly different projections of denitrification, which underscores the potential for stream function to be mischaracterized if tracer tests are interpreted through an inappropriately simplified model for transient storage. More generally, this study highlights the role of modeling in evaluating the experimental design and identifying the potential of system mischaracterization and its implications.

54 ENVIRONMENTAL SCIENCES↗

Mineralized sclerites in the gorgonian coral Leptogorgia chilensis as a natural jamming system

The soft corals (Cnidaria, Octocorallia), a diverse group of colonial marine invertebrates, can reversibly tune their body stiffness in response to external stimuli. This capability is attributed to their dynamic skeletal systems, which consist of thousands of mineralized skeletal elements, called sclerites, embedded within a gel-like matrix that swells/deswells and unjams/jams the sclerites, thus modulating skeletal stiffness. While sclerite morphology is widely used for species identification, its role in the mechanical performance of a soft coral’s skeletal system is largely unknown. Here, we investigated structure-jamming relationships in sclerite-based skeletal architectures using the red gorgonian octocoral Leptogorgia chilensis as a model system. The sclerites of L. chilensis exhibit a shaft-like geometry with two axial branches and two sets of triradiate side branches, which are aligned with the crystallographic symmetry of the constituent magnesium-containing calcite. By combining multiscale three-dimensional (3D) structural characterization, parametric geometrical modeling, 3D printing, mechanical testing, and discrete element simulations, we demonstrate how sclerite geometry achieves a balanced jamming performance in terms of stiffness, weight, strength, and fracture resistance in comparison to alternative geometries parametrically modified from the native sclerites (e.g., changes in the length and number of side branches). Here, we also found that these performance metrics are achieved through the effective interlocking among side and axial branches, which is further enhanced by the fractal-like microscopic spikes on the branch tips. The findings in this natural jamming system offer insights for designing synthetic mechanotunable material architectures for a wide range of applications, from soft robotics to mechanical dampeners.

36 MATERIALS SCIENCE↗

Development of a MOOSE thermal model of the MPC-32 canister and HI-STORM overpack

Nuclear power is a significant source of electricity in the United States, but the average age of nuclear power plants is around 40 years old. Safe management of spent nuclear fuel (SNF) is a key aspect of the back end of the nuclear fuel cycle, and SNF dry storage systems are becoming a popular, effective solution in this area, given the absence of a final disposal system. The spent fuel cask system (dry cask method) provides a feasible solution for maintaining SNF (~60 years) prior to final disposal. This project aims to develop a thermal model of the MPC-32 canister and HI-STORM overpack, using the Multiphysics Object-Oriented Simulation Environment (MOOSE). MOOSE is an open-source framework developed by Idaho National Laboratory (INL) for multiscale, multiphysics simulations. This study will investigate and demonstrate the thermal-hydraulics capabilities of the MOOSE framework, including natural circulation, heat transfer, porous flows, etc. The ultimate goal of the project is to verify whether MOOSE tools (including Pronghorn) can be used to study the thermal performance of the SNF dry cask storage system. This study provides reliable and inclusive solving strategy for dry cask problems. The detailed information about the solving scheme and the governing equations related to the physics of the system is provided in the report. The results for thermal-hydraulic analysis of the HI-STORM system is produced with using open source modules of the MOOSE framework. This results highlights the flexibility and modularity of the MOOSE which makes it a unique candidate for the frameworks and code packages. Therefore, integration of the MOOSE to UNF ST&DARDS will improve the thermal-hydraulic capability of the system while providing distinctive features to users.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Advanced Materials for Energy-Water Systems: The Central Role of Water/Solid Interfaces in Adsorption, Reactivity, and Transport

The structure, chemistry, and charge of interfaces between materials and aqueous fluids play a central role in determining properties and performance of numerous water systems. Sensors, membranes, sorbents, and heterogeneous catalysts almost uniformly rely on specific interactions between their surfaces and components dissolved or suspended in the water-and often the water molecules themselves-to detect and mitigate contaminants. Deleterious processes in these systems such as fouling, scaling (inorganic deposits), and corrosion are also governed by interfacial phenomena. Despite the importance of these interfaces, much remains to be learned about their multiscale interactions. Developing a deeper understanding of the molecular- and mesoscale phenomena at water/solid interfaces will be essential to driving innovation to address grand challenges in supplying sufficient fit-for-purpose water in the future. In this Review, we examine the current state of knowledge surrounding adsorption, reactivity, and transport in several key classes of water/solid interfaces, drawing on a synergistic combination of theory, simulation, and experiments, and finally, provide an outlook for prioritizing strategic research directions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

GLUE Code: A framework handling communication and interfaces between scales

Many scientific applications are inherently multiscale in nature. Such complex physical phenomena often require simultaneous execution and coordination of simulations spanning multiple time and length scales. This is possible by combining expensive small-scale simulations (such as molecular dynamics simulations) with larger scale simulations (such continuum limit/hydro solvers) to allow for considerably larger systems using task and data parallelism. However, the granularity of the tasks can be very large and often leads to load imbalance. Traditionally, we use approximations to streamline the computation of the more costly interactions and this introduces trade-offs between simulation cost and accuracy. In recent years, the available computational power and the advances in machine learning have made computing these scale-bridging interactions and multiscale simulations more feasible. One driving application has been plasma modeling in inertial confinement fusion (ICF), which is fundamentally multiscale in nature. This requires deep understanding of how to extrapolate microscopic information into macroscopically relevant scales. For example, in ICF one needs an accurate understanding of the connection between experimental observables and the underlying microphysics. The properties of the larger scales are often affected by the microscale behavior incorporated usually into the equations of state and ionic and electronic transport coefficients (Liboff, 1959; Rinderknecht et al., 2014; Rosenberg et al., 2015; Ross et al., 2017). Instead of incorporating this information using reliable molecular dynamics (MD) simulations, one often needs to use theoretical models, due to the inability of MD to reach engineering scales (Glosli et al., 2007; Marinak et al., 1998). One approach to resolve this issue is by coupling two MD simulations of different scales via force interpolation, e.g., the AdResS method (Krekeler et al., 2018; Nagarajan et al., 2013). Another approach, which we will pursue in the scope of this work, is by enabling scale bridging between MD simulations and meso/macro-scale models through the development and support of application programming interfaces that these different applications can interact through.

54 ENVIRONMENTAL SCIENCES↗

Revisiting Discharge Mechanism of CF x as a High Energy Density Cathode Material for Lithium Primary Battery

Lithium/fluorinated graphite (Li/CFx) primary batteries show great promise for applications in a wide range of energy storage systems due to their high energy density (>2100 Wh kg –1 ) and low self-discharge rate (<0.5% per year at 25 °C). While the electrochemical performance of the CF x cathode is indeed promising, the discharge reaction mechanism is not thoroughly understood to date. In this article, a multiscale investigation of the CF x discharge mechanism is performed using a novel cathode structure to minimize the carbon and fluorine additives for precise cathode characterizations. Titration gas chromatography, X-ray diffraction, Raman spectroscopy, X-ray photoelectron spectroscopy, scanning electron microscopy, cross-sectional focused ion beam, high-resolution transmission electron microscopy, and scanning transmission electron microscopy with electron energy loss spectroscopy are utilized to investigate this system. Results show no metallic lithium deposition or intercalation during the discharge reaction. Crystalline lithium fluoride particles uniformly distributed with <10 nm sizes into the CF x layers, and carbon with lower sp2 content similar to the hard-carbon structure are the products during discharge. This article deepens the understanding of CF x as a high energy density cathode material and highlights the need for future investigations on primary battery materials to advance performance.

25 ENERGY STORAGE↗

Enabling probabilistic learning on manifolds through double diffusion maps

Here, we present a generative learning framework for probabilistic sampling that extends Probabilistic Learning on Manifolds (PLoM), which is designed to generate statistically consistent realizations of a random vector in a finite-dimensional Euclidean space, informed by a (representative) set of observations. In its original form, PLoM constructs a reduced-order probabilistic model by combining three main components: (a) kernel density estimation to approximate the underlying probability measure, (b) Diffusion Maps to characterize the manifold of the data, and (c) a reduced-order Itô Stochastic Differential Equation (ISDE) to sample from the learned distribution. However, its sampling dynamics are posed in the ambient space and the retained number of reduced coordinates is chosen by projection-reconstruction error. In practice, this often (i) requires more coordinates than the data’s intrinsic dimension to achieve stable sampling and (ii) lacks a smooth, basis-independent lifting back to the data domain; moreover, standard Diffusion Maps emphasize harmonic eigenfunctions and can miss non-harmonic latent structure. We address these limitations by decoupling geometry learning from sampling: a first Diffusion Maps pass identifies non-harmonic coordinates on which we formulate a full-order ISDE directly in the latent space, while Double Diffusion Maps captures multiscale geometric features and Geometric Harmonics (GH) learns a smooth lifting map to the ambient variables that is independent of the particular diffusion basis. This hybrid design preserves the system’s dynamical richness with a compact geometric representation and enables principled out-of-sample inference. The effectiveness and robustness of the proposed method are illustrated through two numerical studies: one based on data generated from two-dimensional Hermite polynomial functions and another based on high-fidelity simulations of a detonation wave in a reactive flow.

Double diffusion maps↗

Cross-Scale Catalyst Modeling Applied to H 2 Storage and Release via Formic Acid

Here, we propose the Systems-to-Atoms (S2A) modeling framework that integrates the kinetics of reaction chemistry and structural configurations across various length scales with the aim of establishing a versatile template for multiscale modeling of reactive flow problems and to predict the operando activity of catalyst materials. The approach encompasses a microkinetic model to analyze surface reactions on individual facets of catalyst nanoparticles coupled with the computation of average surface reaction rates for catalyst nanoparticles of specific size distributions. Macro-homogeneous surface reaction kinetics are derived as a function of catalyst loading and used as input parameters for the continuum-scale reactor model. The cross-scale framework enables the optimization of catalyst utilization through reactor design and operating strategy. To demonstrate the framework, we studied the storage and release of hydrogen from formic acid, a promising liquid organic hydrogen carrier (LOHC), over Pd, Pt, and Cu catalysts. The framework predicts observed trends in formic acid dehydrogenation activity for catalysts with comparable weight loadings and metal particle diameters, demonstrating satisfactory quantitative alignment. Finally, the seamless transmission of parameter uncertainties between scales is also discussed.

08 HYDROGEN↗

Ion Correlations and Partial Ionicities in the Lamellar Phases of Block Copolymeric Ionic Liquids

Recently, significant interest has arisen on the impact of dynamical ion correlations on the conductivity and transport properties of polymeric electrolyte materials. It has been hypothesized that confining ion motion to narrow channels may reduce such ion correlations and enhance the resulting ionic conductivity. Motivated by such considerations, in this study we used a multiscale simulation framework to study the dynamical ion correlations in the microphase-separated lamella phase of block copolymeric ionic liquids and compare with the corresponding results for homopolymeric systems. We probed the influence of ion correlations through the partial ionicity, Δ, which quantifies the ratio of true conductivity to the ideal, Nernst–Einstein conductivity for the anion-related contributions. Consistent with our original hypothesis, our results demonstrate that the partial ionicity relating to the mobile anions is much larger in the lamella phases of block copolymers compared to that in homopolymers. Analysis of the distinct conductivity contributions demonstrates that such results arise as a result of an intricate compensation among the nonideal dynamical correlations relating to anions in lamella phases. Overall, our results suggest that self-assembled phases of block copolymers may provide an avenue to tune the dynamical ion correlations in polymer electrolyte systems.

36 MATERIALS SCIENCE↗

Load‐Balancing Intense Physics Calculations to Embed Regionalized High‐Resolution Cloud Resolving Models in the E3SM and CESM Climate Models

Abstract We design a new strategy to load‐balance high‐intensity sub‐grid atmospheric physics calculations restricted to a small fraction of a global climate simulation's domain. We show why the current parallel load balancing infrastructure of Community Earth System Model (CESM) and Energy Exascale Earth Model (E3SM) cannot efficiently handle this scenario at large core counts. As an example, we study an unusual configuration of the E3SM Multiscale Modeling Framework (MMF) that embeds a binary mixture of two separate cloud‐resolving model grid structures that is attractive for low cloud feedback studies. Less than a third of the planet uses high‐resolution (MMF‐HR; sub‐km horizontal grid spacing) relative to standard low‐resolution (MMF‐LR) cloud superparameterization elsewhere. To enable MMF runs with Multi‐Domain cloud resolving models (CRMs), our load balancing theory predicts the most efficient computational scale as a function of the high‐intensity work's relative overhead and its fractional coverage. The scheme successfully maximizes model throughput and minimizes model cost relative to precursor infrastructure, effectively by devoting the vast majority of the processor pool to operate on the few high‐intensity (and rate‐limiting) high‐resolution (HR) grid columns. Two examples prove the concept, showing that minor artifacts can be introduced near the HR/low‐resolution CRM grid transition boundary on idealized aquaplanets, but are minimal in operationally relevant real‐geography settings. As intended, within the high (low) resolution area, our Multi‐Domain CRM simulations exhibit cloud fraction and shortwave reflection convergent to standard baseline tests that use globally homogenous MMF‐LR and MMF‐HR. We suggest this approach can open up a range of creative multi‐resolution climate experiments without requiring unduly large allocations of computational resources.

54 ENVIRONMENTAL SCIENCES↗

Variational nonlinear WKB in the Eulerian frame

Nonlinear WKB is a multiscale technique for studying locally plane-wave solutions of nonlinear partial differential equations (PDEs). Its application comprises two steps: (1) replacement of the original PDE with an extended system separating the large scales from the small and (2) reduction of the extended system to its slow manifold. In the context of variational fluid theories with particle relabeling symmetry, nonlinear WKB in the mean Eulerian frame is known to possess a variational structure. This much has been demonstrated using, for instance, the theoretical apparatus known as the generalized Lagrangian mean. On the other hand, the variational structure of nonlinear WKB in the conventional Eulerian frame remains mysterious. By exhibiting a variational principle for the extended equations from step (1) above, here we demonstrate that nonlinear WKB in the Eulerian frame is in fact variational. Remarkably, the variational principle for the extended system admits loops of relabeling transformations as a symmetry group. Noether’s theorem therefore implies that the extended Eulerian equations possess a family of circulation invariants parameterized by S 1 . As an illustrative example, we use our results to systematically deduce a variational model of high-frequency acoustic waves interacting with a larger-scale compressible isothermal flow.

97 MATHEMATICS AND COMPUTING↗

Coupled SAM/Griffin Model of a Reference Fluoride-Salt-Cooled High-Temperature Reactor for Multi-Physics Simulations

A multi-physics coupled simulation model of a reference pebble bed fluoride-salt-cooled high-temperature reactor (PB-FHR) has been developed with SAM and Griffin computer codes for transient safety analysis. The reference problem of a prototypical reactor design serves as the foundation for the U.S. NRC (Nuclear Regulatory Commission) to verify the adequacy of computer codes and evaluation models for specific reactor types. In this work, the previously developed SAM model for PB-FHR primary loop has been updated for the coupled simulation. The updated SAM PB-FHR model includes a 2-D axial symmetric core region and external core components in 0-D/1-D. In addition to the primary loop, a detailed model of the RCCS (reactor cavity cooling system) is added. The 2-D and 1-D domains are tightly coupled using the single-solve approach developed in SAM. In the pebble bed region, the SAM multiscale explicit pebble model is applied to calculate the pebble and TRISO fuel kernel temperatures. The Griffin model used in this work is based on a model developed at Idaho National Laboratory in collaboration with the U.S. NRC. The Griffin neutronics model and SAM thermal hydraulics model is coupled with the Comprehensive Reactor Analysis Bundle (CRAB or alternately BlueCRAB) application. Both steady-state and transient scenarios are simulated to demonstrate the model's suitability for multi-physics simulations of PB-FHR transients.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Evaluating grid stress and reliability in future electricity grids across a range of demand, generation mix, and weather trends

The reliability of power grids in the future will depend on how system planners account for the integration of new technologies, extreme weather events, and uncertainties in demand growth from increased electrification and data centers. This study introduces an open-source, multisectoral, multiscale modeling framework that projects grid stress and reliability trends between 2020 and 2055 in the Western Interconnection of the United States. The framework integrates global to national energy-water-land dynamics with power plant siting and hourly grid operations modeling. We analyze future wholesale electricity price shocks and unserved energy events across eight scenarios spanning a range of population growth and economic change, generation mixes, and weather conditions. Our results show future grids with high percentage of non-renewable generation and strong economic growth are characterized by higher reliability and lower wholesale electricity prices than lower growth scenarios because of larger reliance on dispatchable generators and lower fossil fuel extraction costs. Scenarios with high percentage of renewable resources have lower median but more volatile wholesale electricity prices as well as more frequent and severe unserved energy events compared to scenarios relying more on dispatchable generators. These events occur because higher proportion of solar and wind energy causes net demand curves to deepen during midday (duck curves get progressively severe), exacerbating the challenge of meeting demand during summer evening peaks. This study suggests that robust and co-optimized transmission and energy storage planning could help maintain low wholesale electricity prices and high reliability levels in future electricity grids across uncertainties in generation mixes.

Electric grid reliability↗

Dissecting the Genetic Architecture of Carbon Partitioning in Sorghum Using Multiscale Phenotypes

Carbon partitioning in plants may be viewed as a dynamic process composed of the many interactions between sources and sinks. The accumulation and distribution of fixed carbon is not dictated simply by the sink strength and number but is dependent upon the source, pathways, and interactions of the system. As such, the study of carbon partitioning through perturbations to the system or through focus on individual traits may fail to produce actionable developments or a comprehensive understanding of the mechanisms underlying this complex process. Using the recently published sorghum carbon-partitioning panel, we collected both macroscale phenotypic characteristics such as plant height, above-ground biomass, and dry weight along with microscale compositional traits to deconvolute the carbon-partitioning pathways in this multipurpose crop. Multivariate analyses of traits resulted in the identification of numerous loci associated with several distinct carbon-partitioning traits, which putatively regulate sugar content, manganese homeostasis, and nitrate transportation. Using a multivariate adaptive shrinkage approach, we identified several loci associated with multiple traits suggesting that pleiotropic and/or interactive effects may positively influence multiple carbon-partitioning traits, or these overlaps may represent molecular switches mediating basal carbon allocating or partitioning networks. Conversely, we also identify a carbon tradeoff where reduced lignin content is associated with increased sugar content. The results presented here support previous studies demonstrating the convoluted nature of carbon partitioning in sorghum and emphasize the importance of taking a holistic approach to the study of carbon partitioning by utilizing multiscale phenotypes.

59 BASIC BIOLOGICAL SCIENCES↗

DeepM&Mnet: Inferring the electroconvection multiphysics fields based on operator approximation by neural networks

Electroconvection is a multiphysics problem involving coupling of the flow field with the electric field as well as the cation and anion concentration fields. Here, we use electroconvection as a benchmark problem to put forward a new data assimilation framework, the DeepM&Mnet, for simulating multiphysics and multiscale problems at speeds much faster than standard numerical methods using pre-trained neural networks. We first pre-train DeepONets that can predict independently each field, given general inputs from the rest of the fields of the coupled system. DeepONets can approximate nonlinear operators and are composed of two sub-networks, a branch net for the input fields and a trunk net for the locations of the output field. DeepONets, which are extremely fast, are used as building blocks in the DeepM&Mnet and form constraints for the multiphysics solution along with some sparse available measurements of any of the fields. We demonstrate the new methodology and document the accuracy of each individual DeepONet, and subsequently we present two different DeepM&Mnet architectures that infer accurately and efficiently 2D electroconvection fields for unseen electric potentials. Furthermore, the DeepM&Mnet framework is general and can be applied for building any complex multiphysics and multiscale models based on very few measurements using pre-trained DeepONets in a “plug-and-play” mode.

97 MATHEMATICS AND COMPUTING↗