Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale analysis”

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 127 records · Page 7

Initial Development of Fusion Magnet Simulation Capabilities for Performance and Safety Evaluation Using the MOOSE Framework

Fusion energy holds the promise of being a transformative technology as a carbon-neutral, sustainable source of energy. Whole device modeling and the development of fusion digital twins will be increasingly important for emerging fusion device concepts at both national laboratories and within the commercial fusion industry. However, meeting the challenge of whole device modeling of fusion energy devices requires robust, multiphysics, multiscale modeling and simulation technologies capable of running on large-scale supercomputers. Detailed analysis of individual systems at-scale is also required to ensure safe and efficient operation as well as provide the safety basis for future device designs and licensing activities. In a tokamak, toroidal and poloidal magnets confine and shape the fusion plasma to promote the fusion reaction. High plasma temperatures and high magnetic field requirements in modern design concepts (leading to high amounts of energy stored within each magnet) impose electrical, thermal, and mechanical loads on the magnet components, which in turn impacts the safety considerations of the magnet and their supporting systems. Idaho National Laboratory (INL) has a history of working in this space, including development and benchmarking of the Magnetic System Circuitry Analysis Program (MSCAP) and Magnet Arcing (MAGARC) codes to study magnet quench events; notably, MAGARC was used to study quenching during the ITER Engineering Design Activity. However, these legacy codes and capabilities are not parallel and scalable, and new tools are required for future advances in this area, which leads to the INL-developed Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Developed originally for fission reactor systems under United States Department of Energy, Office of Nuclear Energy modeling and simulation programs, the MOOSE framework has been well-suited to multiscale, multiphysics modeling and simulation needs for nuclear systems. The framework is open-source, well-tested, under continuous development and deployment, and developed to a Nuclear Quality Assurance, Level 1 software quality standard. MOOSE has also been used in the fusion space previously in several projects: INL’s Tritium Migration Analysis Program, Version 8 (TMAP8) for tritium migration, UK Atomic Energy Authority’s A Unified Resource for OpenMC (fusion) Reactor Applications (AURORA) code for fusion thermo-mechanical and neutronics analysis, and Argonne National Laboratory’s Cardinal for high-fidelity computational fluid dynamics and neutronics. However, to model superconducting magnets, several MOOSE enhancements are required: additions to the current MOOSE electromagnetic capabilities, new material libraries for superconductors of interest (such as YBCO), as well as fusion-specific models for thermo-mechanics. This talk will discuss initial development activities to build these capabilities in MOOSE, focusing on initial validation and benchmarking activities. Proposed coupling workflows and future work to support the simulation of fusion magnets and magnet structural assemblies for performance and safety evaluation in MOOSE will also be discussed.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Application of NEAMS Multiphysics Framework for Species Tracking in Molten Salt Reactors

This report from Idaho National Laboratory (INL) summarizes the key modeling and simulation activities conducted under the Department of Energy (DOE) Molten Salt Reactor (MSR) Campaign during the Fiscal Year 2023 (FY23). The focus of the work was to leverage state-of-the-art modeling capabilities from the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) codes to enable novel multiphysics and multiscale modeling and simulation of MSRs. Through collaboration with NEAMS code developers, advanced multiphysics analysis capabilities for MSR systems were demonstrated by coupling depletion, thermal-hydraulics, and thermochemistry into an innovative framework for chemical species transport in MSRs. As a result, the framework can track nuclides throughout their lifetimes in the core, from production (depletion) to advection throughout the salt volume (thermal-hydraulics) and off-gassing or precipitation outside of the salt (thermochemistry). This work supports the near-term deployment of MSRs by integrating the synergistic efforts between the DOE’s MSR Campaign and NEAMS program. The resulting framework will help better connect system design modelers with experimentalists to better understand and predict complex physical behaviors in MSRs. Researchers and MSR developers alike can now leverage these new modeling and simulation capabilities to perform novel analyses with applications including: • MSR dynamics during normal operational transients and accident scenarios • Off-gas system design and performance for fuel cycle and depletion analysis • Corrosion and active chemistry control for reactor component health and lifetime determination • Source term, decay heat and activity determination in accident scenarios • Special nuclear material accountancy and chemical forensic analysis for safeguards • Digital twin development of experiments and experimental reactor demonstrations • Measurement requirements for instrumentation and control design • Uncertainty and sensitivity analysis of missing data to inform future experimental data collection.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Improved multiphysics model of the High Temperature Engineering Test Reactor for the simulation of loss-of-forced-cooling experiments

We present a multiphysics model of the High Temperature Engineering Test Reactor for comparison with past and predict future loss-of-forced-cooling (LOFC) experiments. The approach selected combines (1) 3-D full-core superhomogenization-corrected neutronics, (2) 3-D full-core homogenized or semi-heterogeneous heat transfer (macroscale), (3) 2-D axisymmetric fuel rod heat transfer (mesoscale), and (4) 1-D thermal-hydraulics channels. Although large sensitivities remain and the available experimental data are insufficient for use in a rigorous validation exercise, the overall transient behavior is successfully reproduced. The novelty of our work includes (1) a new macroscale-mesoscale heat transfer coupling approach relying on gap conductance to drastically speed up numerical convergence by two orders of magnitude, (2) determination of a radial effective thermal conductivity that yields reasonably close results between the homogenized and semi-heterogeneous macroscale thermal models, and (3) a preliminary study of the reactor's early behavior following a LOFC event, thus enabling further assessment of numerical models against fission power measurements.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Continuous Emulation and Multiscale Visualization of Traffic Flow Using Stationary Roadside Sensor Data

With the advent of the next-generation traffic monitoring systems, there has been a significant increase in the spatial-temporal resolution of vehicle mobility data in many cities. Effective analysis and visualization of such data can provide transportation planners with data-driven insights, which can facilitate the understanding of multiscale traffic dynamics. In this paper, we present a web-based traffic emulator for emulating and visualizing near-real-time and historical traffic flows on highways using data from road-side sensors. To construct a continuous traffic flow, the emulator adopts an analytical pipeline that can (a) integrate traffic data collected from discrete road-side radar detection sensors, (b) interpolate traffic conditions (vehicle speed and volume) on unmeasured road segments based on traffic flow theory, and (c) generate lane-specific vehicle trajectories and movements using a mathematically optimized representation of the road network. Our app also provides an integrated visual workflow that allows users to explore the interconnected traffic dynamics using an appropriate traffic flow visualization selected based on the level of detail. We devise two innovative geo-visualization techniques that utilize an animated strips-network representation and a lane usage matrix to visualize lane performances. To ensure a smooth emulation of large-scale traffic flow in an easy-to-access web environment, we implement the emulator using client-side GPU-accelerated techniques. Lastly, we close with a case study that visualizes traffic dynamics of two scenarios - an afternoon peak hour and a traffic accident - in Chattanooga, Tennessee. Our app visualizes the responses of traffic dynamics during different traffic conditions, and to the presence of the traffic accident at different spatial scales.

42 ENGINEERING↗

Toward predictive permeabilities: Experimental measurements and multiscale simulation of methanol transport in Nafion

A polymer membrane's permeability to solutes determines its suitability for various applications: a permeability value is essential for predicting performance in diverse contexts. Using aqueous methanol permeation through Nafion as an example, we describe a methodology for determining membrane permeability that accounts for boundary layer effects and the possibility of swelling. For the materials and apparatus used herein, analysis of a permeance measurement and computational fluid dynamics simulations show that the mass transfer boundary layer is on the order of ones to tens of microns. Additionally, the data are used to develop and validate a multiscale model describing solute permeation through a hydrated membrane as a series of physical mechanistic steps: reversible adsorption from solution at the membrane interface, diffusion driven by a concentration gradient within the membrane, and reversible desorption into solution at the opposite membrane interface. The validated model is used to predict methanol transport across a solar-driven CO 2 reduction device and to assess the impact of polymer changes on the measured value. The approach of combining experimental data, computational fluid dynamics, and the mechanistic multiscale model is expected to provide more accurate analysis of membrane permeation data in cases with polymer swelling or unusual device geometries, among others.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An efficient protocol for excited states of large biochromophores

Efficient energy transport in photosynthetic antenna is a long-standing source of inspiration for artificial light harvesting materials. However, characterizing the excited states of the constituent chromophores poses a considerable challenge to mainstream quantum chemical and semiempirical excited state methods due to their size and complexity and the accuracy required to describe small but functionally important changes in their properties. In this paper, we explore an alternative approach to calculating the excited states of large biochromophores, exemplified by a specific method for calculating the Q y transition of bacteriochlorophyll a, which we name Chl-xTB. Using a diagonally dominant approximation to the Casida equation and a bespoke parameterization scheme, Chl-xTB can match time-dependent density functional theory’s accuracy and semiempirical speed for calculating the potential energy surfaces and absorption spectra of chlorophylls. We demonstrate that Chl-xTB (and other prospective realizations of our protocol) can be integrated into multiscale models, including concurrent excitonic and point-charge embedding frameworks, enabling the analysis of biochromophore networks in a native environment. We exploit this capability to probe the low-frequency spectral densities of excitonic energies and interchromophore interactions in the light harvesting antenna protein LH2 (light harvesting complex 2). The impact of low-frequency protein motion on interchromophore coupling and exciton transport has routinely been ignored due to the prohibitive costs of including it in simulations. Our results provide a more rigorous basis for continued use of this approximation by demonstrating that exciton transition energies are unaffected by low-frequency vibrational coupling to exciton interaction energies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiscale Reduced Order Modeling and Parameter Estimation for Climate Sciences

Several problems in earth system modeling are dependent on highly multiscale phenomena, such as turbulence, where computational modeling is challenging and expensive. This issue is exacerbated in atmospheric and oceanic domains, due to inherent high-dimensionality of the problem. One approach to this problem has been reduced order modeling (ROM); which aims to represent the key physics of the phenomena as a low-dimensional system. AI methods have huge potential in building accurate, stable ROMs and parameter estimation for these ROMs, as it requires extracting nonlinearities and patterns from simulation and/or observational data. Developing physics-based AI approaches specialized for the complexities of multiscale data, along with strategies to account for uncertainties, will revolutionize rapid modeling, analysis and decision making for earth system problems of practical interest.

58 GEOSCIENCES↗

Multidimensional Diagnostics of the Interface Evolutions in Solid-State Lithium Batteries

All-solid-state lithium metal batteries hold tremendous potential as next-generation energy storage solutions. However, their performance is hindered by complex chemo-mechanical instabilities at the interface, which encompass electrolyte decomposition, loss of contact due to volume changes, and the growth of lithium filaments. Given the complexity of these failure mechanisms, it is imperative to employ multidimensional diagnostics that leverage a combination of tools capable of quantifying formed voids and filaments, identifying the chemical and mechanical properties of lithium dendrites and electrolyte decomposition products, as well as monitoring the evolution of these processes in situ. The primary objective of this project is to develop a comprehensive characterization platform for solid-state battery interfaces, integrating FIB-SEM tomography, ToF-SIMS, in-SEM nanoindentation, and in-SIMS AFM-based stiffness mapping to enable structural, chemical, and mechanical analyses. The project team has accomplished significant milestones, including: 1) The successful development of a characterization platform that combines FIB-SEM, ToF-SIMS, nanoindentation, and AFM. 2) The creation of solid-state cells at various scales. 3) Real-time visualization of the plating and stripping of lithium. 4) Operando analysis of the chemical composition and evolution of electrolyte decomposition products, encompassing intermediate and metastable species, at both the cathode and anode interfaces. 5) In situ analysis of mechanical properties coupled with chemical mapping. 6) The implementation of the multiscale characterization technique has yielded valuable insights into cell failure mechanisms and has led to the development of a solid-state lithium cell that exhibits exceptional durability, surpassing 2500 cycles. This project has resulted in 10 publications, 1 invention disclosure, 8 manuscripts under review or in preparation.

25 ENERGY STORAGE↗

Validation of a Global Geospace Model With a Systems Science Approach Based on Canonical Correlation Analysis

A systems science approach based on canonical correlation analysis (CCA) is applied as a new, behavioral way to validate global geospace models. The biggest novelty of the technique is that it validates models at a system level, whereby a side‐by‐side comparison is performed of CCA applied to a 30‐day observational and the corresponding simulation data sets comprising quiet, moderate and active times. The simulation used the Multiscale Atmosphere‐Geospace Environment (MAGE) model. It is shown that (a) CCA must be combined with sensitivity analysis to be effective, (b) the MAGE model generally reproduces the observed behavior (more so for quieter time intervals), quantified by the intercorrelations between different variables and (c) the technique identifies the SuperMAG SML index as a quantity for which refinements of the model are needed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Separability of microtearing mode and electron temperature gradient turbulence regimes

The separability of microtearing mode (MTM)-dominated and electron temperature gradient (ETG)-driven turbulence regimes is studied with multiscale nonlinear gyrokinetic simulations. The simulations are based on National Spherical Torus Experiment-like, high-confinement mode pedestal parameters, where electromagnetic perturbations are large. Linear analysis indicates a wide scale-separation between the MTM and ETG modes in binormal wavenumber space (perpendicular to the magnetic field line), with no unstable modes at intermediate scales. Likewise, single-scale nonlinear analyses, retaining ion-only or electron-only spatio-temporal scales, produce seemingly well-converged transport states. Surprisingly, the multiscale simulation, which contains both the ion and electron scales, closely follows the transport from the electron-scale simulation. This trend is robust over a wide range of electron temperature gradient. Remarkably, compared to ion-scale simulations, MTM turbulence is significantly reduced at multiscale resolution even when ETG turbulence is low. In this case, traditional ion-scale resolution overestimates the electron energy flux, and it is not possible to accurately simulate the MTM turbulence with separable ion-scale simulations. Here, while the analysis confirms the validity of electron-scale simulations for predicting the electron transport, it also indicates that multiscale simulation may be required for reproducing the turbulence spectrum for systems with coupled MTM-ETG turbulence.

H-mode↗

Beyond Price-Taker: Multiscale Optimization of Wind and Battery Integrated Energy Systems

Integrating renewable energy into the electric grid is challenging due to the intermittency and variability of wind and other non-dispatchable resources. Integrated energy systems (IESs) combine multiple energy technologies (e.g., fossil, nuclear, renewables, storage) to reduce costs and improve flexibility and reliability. However, standard techno-economic analysis (TEA) methods often overestimate the benefits of IESs because they fail to account for energy market adjustments. This paper systematically studies the limitations of the prevailing price-taker assumption for TEA and optimization of hybrid energy systems. As an illustrative case study, we retrofit an existing wind farm in the RTS-GMLC test system (which loosely mimics the Southwest U.S.) with battery energy storage to form an IES. We show that the standard price-taker model overestimates the electricity revenue and the net present value (NPV) of the IES up to 178% and 30.4%, respectively, compared to our more rigorous multiscale optimization. These differences arise because introducing storage creates a more flexible resource that impacts the larger wholesale electricity market. Moreover, this work highlights the impact of the IES has on the market via various strategic bidding, and underscores the importance of moving beyond price-taker for optimal storage sizing and TEA of IESs. We conclude by discussing opportunities to generalize the proposed framework to other IESs, and highlight emerging research questions regarding the complex interactions between IESs and markets.

25 ENERGY STORAGE↗

tell: a Python package to model future total electricity loads in the United States

The purpose of the Total ELectricity Load (tell) model is to generate 21st century profiles of hourly electricity load (demand) across the Conterminous United States (CONUS). tell loads reflect the impact of climate and socioeconomic change at a spatial and temporal resolution adequate for input to an electricity grid operations model. tell uses machine learning to develop profiles that are driven by projections of climate/meteorology and population. tell also harmonizes its results with United States (U.S.) state-level, annual projections from a national- to global-scale energy-economy model. This model accounts for a wide range of other factors affecting electricity demand, including technology change in the building sector, energy prices, and demand elasticities, which stems from model coupling with the U.S. version of the Global Change Analysis Model (GCAM-USA). tell was developed as part of the Integrated Multisector Multiscale Modeling (IM3) project. IM3 explores the vulnerability and resilience of interacting energy, water, land, and urban systems in response to compound stressors, such as climate trends, extreme events, population, urbanization, energy system transitions, and technology change

24 POWER TRANSMISSION AND DISTRIBUTION↗

Probing the Influence of Multiscale Heterogeneity on Effective Properties of Graphite Electrodes

Graphite electrodes in the lithium-ion battery exhibit various particle shapes, including spherical and platelet morphologies, which influence structural and electrochemical characteristics. It is well established that porous structures exhibit spatial heterogeneity, and the particle morphology can influence transport properties. The impact of the particle morphology on the heterogeneity and anisotropy of geometric and transport properties has not been previously studied. Herein, this study characterizes the spatial heterogeneities of 18 graphite electrodes at multiple length scales by calculating and comparing the structural anisotropy, geometric quantities, and transport properties (pore-scale tortuosity and electrical conductivity). We found that the particle morphology and structural anisotropy play an integral role in determining the spatial heterogeneity of directional tortuosity and its dependency on pore-scale heterogeneity. Our analysis reveals that the magnitude of in-plane and through-plane tortuosity difference influences the multiscale heterogeneity in graphite electrodes.

25 ENERGY STORAGE↗

Quantifying the Agyrotropy of Proton and Electron Heating in Turbulent Plasmas

An important aspect of energy dissipation in weakly collisional plasmas is that of energy partitioning between different species (e.g., protons and electrons) and between different energy channels. Here we analyse pressure–strain interaction to quantify the fractions of isotropic compressive, gyrotropic, and nongyrotropic heating for each species. An analysis of kinetic turbulence simulations is compared and contrasted with corresponding observational results from Magnetospheric Multiscale Mission data in the magnetosheath. In assessing how protons and electrons respond to different ingredients of the pressure–strain interaction, we find that compressive heating is stronger than incompressive heating in the magnetosheath for both electrons and protons, while incompressive heating is stronger in kinetic plasma turbulence simulations. Concerning incompressive heating, the gyrotropic contribution for electrons is dominant over the nongyrotropic contribution, while for protons nongyrotropic heating is enhanced in both simulations and observations. Variations with plasma β are also discussed, and protons tend to gain more heating with increasing β.

79 ASTRONOMY AND ASTROPHYSICS↗

Understanding Solution State Conformation and Aggregate Structure of Conjugated Polymers via Small Angle X-ray Scattering

Donor-acceptor (D-A) conjugated polymers are high-performance organic electronic materials that exhibit complex aggregation behavior. Understanding the solution state conformation and aggregation of conjugated polymers is crucial for controlling morphology during thin-film deposition and the subsequent electronic performance. However, a precise multiscale structure of solution state aggregates is lacking. Here, we present an in-depth small-angle X-ray scattering (SAXS) analysis of the solution state structure of an isoindigo-bithiophene-based D-A polymer (PII-2T) in chlorobenzene and decane as our primary system. Modeling the system as a combination of hierarchical fibrillar aggregates mixed with dispersed polymers, we extract information about conformation and multiscale aggregation and also clarify the physical origin of features often observed but unaddressed or misinterpreted in small-angle scattering patterns of conjugated polymers. The persistence length of the D-A polymer extracted from SAXS agrees well with a theoretical model based on the dihedral potentials. Additionally, we show that the broad high q structure factor peak seen in scattering profiles can be attributed to lamellar stacking occurring within the fibril aggregates and that the low q aggregate scattering is strongly influenced by the polymer molecular weight. Overall, the SAXS profiles of D-A polymers in general exhibit a sensitive dependence on the co-existence of fibrillar aggregate and dispersed polymer chain populations. We corroborate our findings from SAXS with electron microscopy of freeze-dried samples for direct imaging of fibrillar aggregates. Finally, we demonstrate the generality of our approach by fitting the scattering profiles of a variety of D-A polymers based on thieno-isoindigo (PTII-2T), diketopyrrolopyrrole (DPP2T-TT, DPP-BTZ, PDPP2FT-C- 16 ), naphthalenediimide (P(NDI2OD-T2)), and a conjugated block copolymer P3HT-b-DPPT-T. The results presented here establish a picture of the D-A polymer solution state structure and provide a general method of interpreting and analyzing their scattering profiles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine learning from RANS and LES to inform coarse grid simulations

Nuclear system thermal hydraulic analysis has historically relied on computationally inexpensive 1D codes. However, such tools are unable to capture multiscale multidimensional effects in large nuclear reactor enclosures. On the other hand, simulations with higher fidelity can be too expensive for such purposes. One of the ways to reduce computational cost is to perform simulations on a coarse grid, which, unfortunately, introduces large discretization errors. In this paper, two high-to-low data-driven approaches are investigated: (1) a coarse grid turbulence model to predict eddy viscosity and (2) correction of errors in coarse grid velocity fields. The approaches aim to reduce grid- and turbulence model-induced errors in coarse grid Reynolds-averaged Navier–Stokes (RANS) simulations. Two sources of high-fidelity data, RANS and large eddy simulations (LES), are explored. To extract the eddy viscosity from the LES data, an inverse optimization problem is solved. However, the LES eddy viscosity is shown to be comparable to the RANS eddy viscosity in terms of error reduction. Therefore, the directly available RANS eddy viscosity was used to develop a coarse grid data-driven turbulence model. Additionally, error correction in velocity is used to reduce the remaining uncertainties and bring the results closer to reality. In conclusion, the performance of the frameworks is demonstrated for a scaled upper plenum of a gas-cooled reactor facility.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

WaterTAP v1.0.0

WaterTAP is a new open source library of water-specific models built on the IDAES Integrated Platform to support the design and optimization of integrated water treatment systems, improve existing systems and enable the analysis of new designs incorporating emerging technologies. WaterTAP enables advanced modeling and optimization for: conceptual and multiscale design, dynamics and model predictive control of water systems and networks.

Beattie, KeithS↗

Predicting Morphological Evolution during Coprecipitation of MnCO 3 Battery Cathode Precursors Using Multiscale Simulations Aided by Targeted Synthesis

The performance of lithium-ion batteries is intimately linked to both the structure and the morphology of the cathode material, which in turn is critically linked to the synthesis conditions. However, few studies focus on understanding synthesis, especially during the coprecipitation of metal oxide precursors, a process that largely determines the final morphology of the material. In this paper, we go beyond the typical equilibrium particle shape analysis conducted in the literature and incorporate kinetic aspects of morphology evolution. We perform these studies using controlled synthesis on a well-defined metal salt system (MnCO 3 ) combined with multiscale simulations and high-resolution microscopy. Results show that with increasing metal concentration, the particles transition from rhombohedral to cubic to spherical shapes. Computational analysis using density functional theory (DFT) reveals that rhombohedral shaped particles evolve under equilibrium conditions. Phase field techniques indicate that at higher metal concentrations, fast growth kinetics of the precipitates result in the transition to cubic and, subsequently, spherical shapes, accompanied by a decrease in particle size. This study, while limited to the one metal salt system, provides an approach to shed light on the synthesis process of mixed transition metal salts, gradient materials, and other cathode materials of interest to the battery community.

25 ENERGY STORAGE↗