Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale modeling and 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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quantum-Accurate Multiscale Modeling of Shock Hugoniots, Ramp Compression Paths, Structural and Magnetic Phase Transitions, and Transport Properties in Highly Compressed Metals

Fully characterizing high energy density (HED) phenomena using pulsed power facilities (Z machine) and coherent light sources is possible only with complementary numerical modeling for design, diagnostic development, and data interpretation. The exercise of creating numerical tests, that match experimental conditions, builds critical insight that is crucial for the development of a strong fundamental understanding of the physics behind HED phenomena and for the design of next generation pulsed power facilities. The persistence of electron correlation in HED materials arising from Coulomb interactions and the Pauli exclusion principle is one of the greatest challenges for accurate numerical modeling and has hitherto impeded our ability to model HED phenomena across multiple length and time scales at sufficient accuracy. An exemplar is a ferromagnetic material like iron, while familiar and widely used, we lack a simulation capability to characterize the interplay of structure and magnetic effects that govern material strength, kinetics of phase transitions and other transport properties. Herein we construct and demonstrate the Molecular-Spin Dynamics (MSD) simulation capability for iron from ambient to earth core conditions, all software advances are open source and presently available for broad usage. These methods are multi-scale in nature, direct comparisons between high fidelity density functional theory (DFT) and linear-scaling MSD simulations is done throughout this work, with advancements made to MSD allowing for electronic structure changes being reflected in classical dynamics. Main takeaways for the project include insight into the role of magnetic spins on mechanical properties and thermal conductivity, development of accurate interatomic potentials paired with spin Hamiltonians, and characterization of the high pressure melt boundary that is of critical importance to planetary modeling efforts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Examining the effects of soil entrainment during nuclear cloud rise on fallout predictions using a multiscale atmospheric modeling framework

Current operational models for nuclear cloud rise over land were developed and validated using observations from shallow-buried or surface detonations, where lofted soil quickly mixed with fission products from the detonation. These models poorly predict fallout from elevated detonations near the fallout-free height of burst (FFHOB), where interactions with the ground are limited and the mixing of fission products and lofted soil is incomplete. Fallout-free is a misnomer at this HOB, as fallout was observed in these cases, but was below the levels of concern, especially off-grounds of the nuclear test site. To correctly characterize and model fallout from detonations near the FFHOB, models must be developed which can capture the stratified nature of the particle and activity-size distributions within the cloud. Previously, it was shown that the Weather Research and Forecasting (WRF) model can accurately simulate nuclear cloud rise for airbursts with little to no ground interactions (Arthur et al., 2021). That work is expanded here by (1) using a radiation-hydrodynamics code to improve the fireball initialization in WRF, (2) further developing an aerosol package from WRF-Chem to simulate lofted soil, and (3) combining the WRF cloud rise simulations with the operational models used at the National Atmospheric Release Advisory Center (NARAC) for fallout modeling. Using this combination of codes, the Upshot-Knothole Grable detonation, which was just below the FFHOB, is simulated from seconds after detonation through cloud rise and fallout, and results are compared to historical test data. Here, the results show improved prediction of dose rate and highlight the need to correctly characterize the entrainment of material into the cloud and the subsequent mixing of fission products with entrained material.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

Multiscale-Informed Modeling of High Temperature Component Response with Uncertainty Quantification

This report summarizes a joint effort between Argonne National Laboratory, Idaho National Laboratory, and Los Alamos National Laboratory to develop and deploy constitutive models targeted at predicting the life of Grade 91 alloy components subjected to high temperature environments typical of those that structural components in advanced nuclear reactors would experience. Two distinct, but complementary constitutive modeling approaches have been taken here. The first employs a phenomenological viscoplastic model for which parameters have been calibrated based on experimental data for a wide range of Grade 91 alloy that has undergone a variety of processing. A Bayesian approach was used to derive distributions of uncertain parameters for this model based on this data set. The second approach is a reduced order model suitable for engineering-scale analysis that is based on the results of a large set of mesoscale simulations. Mesoscale models allow for the microstructure and composition of a particular alloy to be directly taken into account in the computation of the viscoplastic response, but are computationally expensive, which makes it impractical to directly call those models for the material constitutive response in an engineering-scale simulation. The reduced-order representation of the response of the underlying model used here allows for an engineering-scale model to take into account the characteristics of the underlying microstructure, while only incurring a reasonable computational expense. Both of these approaches have different strengths, and are applicable for different parts of the design/analysis process. The phenomenological models can be readily parameterized based on a set of experimental data for a given class of materials and used for scoping calculations. Once a specific material is chosen and adequately characterized, the reduced order models can accurately predict the response of that specific alloy, and because the models are based on predictive models of the underlying microstructure, they can be used to more confidently predict the response under conditions in regions where there is limited experimental data. Both of these models have been integrated in the Grizzly code, which is used here to perform proof-of-concept uncertainty quantification analyses of a simple component under prototypical conditions. The built- in stochastic analysis capabilities in the MOOSE framework that Grizzly is built on are used here to run large sets of simulations for this uncertainty quantification analysis. As would be expected, because the reduced order models are developed for a much more tightly defined alloy, they predict tighter distributions of the time to failure than the phenomenological models, which are calibrated to a broader set of data. Also important is that these simulations demonstrate that a reduced order modeling approach can be successfully deployed to propagate uncertainties from the material scale to practical engineering-scale component simulations.

42 ENGINEERING↗
DOE OSTI2025

Learning Physically Interpretable Atmospheric Models From Data With WSINDy

The multiscale and turbulent nature of Earth's atmosphere has historically rendered accurate weather modeling a hard problem. Recently, there has been an explosion of interest surrounding data-driven approaches to weather modeling, which in many cases show improved forecasting accuracy and computational efficiency when compared to traditional methods. However, many of the current data-driven approaches employ highly parameterized neural networks, often resulting in uninterpretable models and limited gains in scientific understanding. In this work, we address the interpretability problem by explicitly discovering partial differential equations governing atmospheric phenomena, identifying symbolic mathematical models with direct physical interpretations. The purpose of this paper is to demonstrate that, in particular, the weak-form sparse identification of nonlinear dynamics (WSINDy) algorithm can learn effective atmospheric models from both simulated and assimilated data. Our approach adapts the standard WSINDy algorithm to work with high-dimensional fluid data of arbitrary spatial dimension.

58 GEOSCIENCES↗
DOE OSTI2025

COCOMO2: A Coarse-Grained Model for Interacting Folded and Disordered Proteins

Biomolecular interactions are essential in many biological processes, including complex formation and phase separation processes. Coarse-grained computational models are especially valuable for studying such processes via simulation. Here, we present COCOMO2, an updated residue-based coarse-grained model that extends its applicability from intrinsically disordered peptides to folded proteins. This is accomplished with the introduction of a surface exposure scaling factor, which adjusts interaction strengths based on solvent accessibility, to enable the more realistic modeling of interactions involving folded domains without additional computational costs. COCOMO2 was parametrized directly with solubility and phase separation data to improve its performance on predicting concentration-dependent phase separation for a broader range of biomolecular systems compared to the original version. COCOMO2 enables new applications including the study of condensates that involve IDPs together with folded domains and the study of complex assembly processes. COCOMO2 also provides an expanded foundation for the development of multiscale approaches for modeling biomolecular interactions that span from residue-level to atomistic resolution.

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

Correction to “COCOMO2: A Coarse-Grained Model for Interacting Folded and Disordered Proteins”

Biomolecular interactions are essential in many biological processes, including complex formation and phase separation processes. Coarse-grained computational models are especially valuable for studying such processes via simulation. Here, we present COCOMO2, an updated residue-based coarse-grained model that extends its applicability from intrinsically disordered peptides to folded proteins. This is accomplished with the introduction of a surface exposure scaling factor, which adjusts interaction strengths based on solvent accessibility, to enable the more realistic modeling of interactions involving folded domains without additional computational costs. COCOMO2 was parametrized directly with solubility and phase separation data to improve its performance on predicting concentration-dependent phase separation for a broader range of biomolecular systems compared to the original version. COCOMO2 enables new applications including the study of condensates that involve IDPs together with folded domains and the study of complex assembly processes. COCOMO2 also provides an expanded foundation for the development of multiscale approaches for modeling biomolecular interactions that span from residue-level to atomistic resolution.

Molecular interactions↗
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 OSTI2023

Model America - Arizona extract from ORNL's AutoBEM v1.1

Oak Ridge National Laboratory (ORNL) has developed the Automatic Building Energy Modeling (AutoBEM) software suite to process multiple types of data, extract building-specific descriptors, generate building energy models, and simulate them on High Performance Computing (HPC) resources. For more information, see AutoBEM-related publications (bit.ly/AutoBEM).Two sets of sample data are provided for 2,555,152 buildings located within the boundary of Arizona in the United States:Data (846.3MB *.csv) - minimalist list of each building (rows) for the following fields (columns) • ID - unique building ID • Centroid - building center location in latitude/longitude (from Footprint2D) • Footprint2D - building polygon of 2D footprint (lat1/lon1_lat2/lon2_...) • State_abbr - state name • Area - estimate of total conditioned floor area (ft2) • Area2D - footprint area (ft2) • Height - building height (ft) • NumFloors - number of floors (above-grade) • WWR_surfaces - percent of each facade (pair of points from Footprint2D) covered by fenestration/windows (average 14.5% for residential, 40% for commercial buildings) • CZ - ASHRAE Climate Zone designation • BuildingType - DOE prototype building designation (IECC=residential) as implemented by OpenStudio-standards • Standard - building vintage • Sample Models (114GB*.zip by county) - OpenStudio and EnergyPlus building energy models named according to IDThis data is made free and openly available in hopes of stimulating any simulation-informed use case. Data is provided as-is with no warranties, express or implied, regarding fitness for a particular purpose. We wish to thank our sponsors which include Oak Ridge National Laboratory (ORNL), U.S. Dept. of Energy’s (DOE) Building Technologies Office (BTO), Office of Electricity (OE), and Biological and Environmental Research (BER).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Computer Simulation of Proton Transport in Fuel Cell Membranes (Final Report)

This DOE-supported research grant focused on understanding the nature of proton transport in complex systems such as proton exchange membranes (PEMs). The most unique aspect of the research was the development and implementation of a novel multiscale reactive molecular dynamics (MS-RMD) methodology. In this approach, covalent bonds can dynamically break and form, allowing one to accurately treat the proton hopping process essential to capturing the physics of proton transport. The Voth group applied this method to proton exchange membrane systems, providing insight into their proton transport mechanism. They found that protons can diffuse most rapidly in the water-rich regions, but that protons actually spend so little time in such regions that transport along the hydrophobic – hydrophilic interface controls the membrane performance. In addition, the group worked to increase understanding of acidic solutions, developing methods for simulating and interpreting experimental infrared vibrational spectroscopy for excess protons (acidic solutions). The group also implemented novel tools for developing proton transport reactive MD models using a relative entropy minimization scheme.

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

State of the art, gaps, and prospects in fusion materials theory and modelling

Advancing the theory and simulation of materials for fusion applications remains a key component of global roadmaps aimed at delivering much-needed fusion power. Especially as the drive for commercial application increases, prototypes must be designed against radiation damage before the relevant experimental data can be collected and cost reductions that are possible by testing materials in silico become even more important. Here, we summarise the state of the art as it emerged during the 7 th Fusion Materials Theory & Modelling Workshop that took place in 2024, with the aim to highlight present gaps and future directions for the fusion materials modelling community. Of particular interest were the effects of transmutations, chemical complexity with the development of novel alloys and interatomic potentials, advancements in modelling high-dose microstructures, comparison with experimental data and multiscale models for structural assessment relying on high-performance computing and virtual reality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Multi-GPU immersed boundary method hemodynamics simulations

Large-scale simulations of blood flow that resolve the 3D deformation of each comprising cell are increasingly popular owing to algorithmic developments in conjunction with advances in compute capability. Among different approaches for modeling cell-resolved hemodynamics, fluid structure interaction (FSI) algorithms based on the immersed boundary method are frequently employed for coupling separate solvers for the background fluid and the cells within one framework. GPUs can accelerate these simulations; however, both current pre-exascale and future exascale CPU-GPU heterogeneous systems face communication challenges critical to performance and scalability. In this paper, we describe, to our knowledge, the largest distributed GPU-accelerated FSI simulations of high hematocrit cell-resolved flows with over 17 million red blood cells. We compare scaling on a fat node system with six GPUs per node and on a system with a single GPU per node. Through comparison between the CPU- and GPU-based implementations, we identify the costs of data movement in multiscale multi-grid FSI simulations on heterogeneous systems and show it to be the greatest performance bottleneck on the GPU.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

IM3 Data Center Driven Grid Stress Dataset for the U.S. Western Interconnection

This dataset provides projected grid stress and reliability results (including all model inputs and outputs from an open-source grid operations modeling framework - GO), for the Integrated Multisector, Multiscale Modeling (IM3) project, under varying levels of data center demand growth between 2025 and 2035 in the U.S. Western Interconnection. The scenarios and sensitivity experiments are combinations of different data center demand growth rates and energy, weather, population and economic pathways. Data center demand growth projections were sourced from the Electric Power Research Institute (EPRI). The data center demand growth projection names are: Low (3.71% annual data center demand growth) Moderate (5% annual data center demand growth) High (10% annual data center demand growth) Higher (15% annual data center demand growth) Energy, weather, population and economic pathways are informed by two Shared Socioeconomic Pathways (SSP3 and SSP5) and two Representative Concentration Pathways (RCP4.5 and RCP8.5) following the hotter general circulation model (GCM) forcing group from a set of perturbed thermodynamics simulations. The resulting pathway names are: rcp45hotter_ssp3 rcp45hotter_ssp5 rcp85hotter_ssp3 rcp85hotter_ssp5 The main scenarios and sensitivity experiments are detailed below. Reference scenario: The projected grid stress and reliability results for the U.S. Western Interconnection from a previous study. This scenario does not consider data center demand growth explicitly. Data center scenario: Building on the reference scenario, this scenario considers various data center growth rates and how they impact the U.S. Western Interconnection. Data center loads are modeled as flat 8760-hr profiles. This scenario does not consider new generation and transmission capacities specifically designed to meet the new data center demands. The related folder is named "flat". Delayed generator retirements sensitivity experiment: Building on the data center scenario, this experiment explores the impact of different levels of natural gas and nuclear generator retirement delays. The resulting scenario names are: (1) postponing 100% nuclear retirements; (2) postponing 100% nuclear and 25% natural gas retirements; (3) postponing only 50% natural gas retirements; (4) postponing 100% nuclear and 50% natural gas retirements; (5) postponing 100% nuclear and 75% natural gas retirements; and (6) postponing 100% nuclear and 100% natural gas retirements. The related folder names are: no_gen_retire_0_gas, no_gen_retire_25_gas, no_gen_retire_50_gas, no_gen_retire_50_gas_only, no_gen_retire_75_gas, and no_gen_retire_100_gas. Demand response through curtailment sensitivity experiment: Building on the data center scenario, this experiment explores the impact of different participation and compensation levels of data center demand response. The resulting scenario names are: (1) 5% demand available for curtailment with 750 $/MWh compensation; (2) 5% demand available for curtailment with 500 $/MWh compensation; (3) 5% demand available for curtailment with 250 $/MWh compensation; (4) 15% demand available for curtailment with 750 $/MWh compensation; (5) 15% demand available for curtailment with 500 $/MWh compensation; and (6) 15% demand available for curtailment with 250 $/MWh compensation. The related folder names are: dr_cost_250_drup_0_drdown_5, dr_cost_250_drup_0_drdown_15, dr_cost_500_drup_0_drdown_5, dr_cost_500_drup_0_drdown_15, dr_cost_750_drup_0_drdown_5, and dr_cost_750_drup_0_drdown_15. Combination of delayed generator retirements and demand response through curtailment sensitivity experiment: The impact of combining postponing 100% nuclear and 25% natural gas retirements with 5% demand available for curtailment with 750 $/MWh compensation is simulated. The related folder is named "dr_cost_750_drup_0_drdown_5_nuc_100_gas_25". Please refer to the README file for a detailed description of the dataset including individual files and references.

Artificial Intelligence↗
DOE OSTI2022

LigninGraphs: lignin structure determination with multiscale graph modeling

Lignin is an aromatic biopolymer found in ubiquitous sources of woody biomass. Designing and optimizing lignin valorization processes requires a fundamental understanding of lignin structures. Experimental characterization techniques, such as 2D-heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR) spectra, could elucidate the global properties of the polymer molecules. Computer models could extend the resolution of experiments by representing structures at the molecular and atomistic scales. We introduce a graph-based multiscale modeling framework for lignin structure generation and visualization. The framework employs accelerated rejection-free polymerization and hierarchical Metropolis Monte Carlo optimization algorithms. We obtain structure libraries for various lignin feedstocks based on literature and new experimental NMR data for poplar wood, pinewood, and herbaceous lignin. The framework could guide researchers towards feasible lignin structures, efficient space exploration, and future kinetics modeling. Its software implementation in Python, LigninGraphs, is open-source and available on GitHub.

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

Toward Addressing the Challenge to Predict the Heat Capacities of RDX and HMX Energetic Materials

Availability of heat capacity as function of pressure and temperature is an essential prerequisite for development of a computational multiscale strategy capable to address the evolution of microstructure and energy release in advanced high energy density materials. In the case of 1,3,5-trinitro-1,3,5-triazinane (RDX) and 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) systems as two of the most studied energetic materials, there are substantial gaps in experimental data, with available heat capacities values distributed only in a region close to standard ambient conditions. In this study we demonstrate how these major experimental limitations can be addressed in the case of the RDX and HMX systems based on the combined use of classical and quantum mechanical calculations. We show that by considering ideal gas properties evaluated using quantum mechanical methods, and residual properties obtained from molecular simulations using fully flexible atomistic force field models, excellent agreement can be obtained for the predicted heat capacities to the most recent experimental values. An important advantage of the current computational methodology is that it allows evaluation of both constant-volume and constant-pressure heat capacities for a broad interval of temperatures and pressures, which encompasses solid and liquid phases conditions. In the case of the solid α and γ phases of RDX and the β phase of HMX, the predicted results follow closely both the available experimental data at standard ambient conditions and the results obtained using density functional theory calculations at high pressures, a regime where experimental data are not available. A perspective to expand the current methodology is also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Coupled Monte Carlo and thermal-hydraulics modeling of a prismatic gas reactor fuel assembly using Cardinal

Cardinal is a MOOSE application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics to the MOOSE framework, closing the neutronics and thermal-fluid gaps in conducting tightly-coupled, high-resolution multiscale and multiphysics analyses. By leveraging MOOSE's interfaces for wrapping external codes, Cardinal overcomes many challenges encountered in earlier multiphysics coupling works, such as file-based I/O or overly-restrictive geometry mapping requirements. In this work, we leverage a subset of the multiphysics interfaces in Cardinal to perform coupling of OpenMC neutron transport, MOOSE heat conduction, and THM thermal-fluids for steady-state modeling of a prismatic gas reactor fuel assembly.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

xSDK: Building an ecosystem of highly efficient math libraries for exascale

Current efforts to build increasingly powerful computer architectures are opening up new avenues for more complex and higher fidelity simulations coupled with data analytics and learning, leading to new scientific insights and deeper understanding. At one extreme, exascale computers will be much faster than previous computer generations (performing 10 18 operations per second—that is, 1,000 times faster than petascale). To achieve these performance improvements, computer architectures are becoming increasingly complex, with deep memory hierarchies, very high node and core counts, and heterogeneous features such as graphics processing units (GPUs). Such architectural changes impact the full breadth of computing scales, as heterogeneity pervades even current-generation laptops, workstations, and moderate-sized clusters. While emerging advanced architectures provide unprecedented opportunities, they also present significant challenges for developers of scientific applications, such as multiphysics and multiscale codes, who must adapt their software to handle disruptive changes in architectures and new programming models that have not yet stabilized. Developers must consider increasing concurrency while reducing communication and synchronization, and other complexities such as the potential for using mixed precision to leverage the compute power available in low-precision tensor cores. On one hand, developers must implement new scientific capabilities, which in turn increase code complexity. On the other hand, the codes must be ported to new architectures, requiring the inclusion of new programming models and the restructuring of code to achieve good performance. Addressing these issues is beyond the capability of any single person or team—leading to the need for collaboration among many teams, who encapsulate their expertise in reusable software and work together to create sustainable software ecosystems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Evaluating Mesoscale Convective Systems Over the US in Conventional and Multiscale Modeling Framework Configurations of E3SMv1

Organized mesoscale convective systems (MCSs) contribute a significant amount of precipitation in the Central and Eastern US during spring and summer, which impacts the availability of freshwater and flooding events. However, current global Earth system models cannot capture MCSs well and misrepresent the statistics of precipitation in the region. In this study, we investigate the representation of MCSs in three configurations of the Energy Exascale Earth System Model (E3SMv1) by tracking individual storms based on outgoing longwave radiation using a new application of TempestExtremes. Our results indicate that conventional parameterizations of convection, implemented in both low (LR; ~150 km) and high (HR; ~25 km) resolution configurations, fail to capture almost all MCS-like events, in-part because they underestimate high-level cloud ice associated with deep convection. On the other hand, the multiscale modeling framework (MMF; cloud-resolving models embedded in each grid-column of ~150 km resolution E3SMv1) configuration represents MCSs and their annual cycle better. Nevertheless, relative to observations, the E3SMv1-MMF spatial distribution of MCSs and associated precipitation is shifted eastward, and the diurnal timing is lagged. A comparison between the large-scale environment in E3SMv1-MMF and ERA5 reanalysis suggests that the biases during the summer in E3SMv1-MMF are associated with biases in low-level humidity and meridional moisture transport within the low-level jet. The fact that conventional parameterizations of convection, even with high-resolution, cannot capture MCSs over the US suggests that methods with explicit representation of kilometer-scale convective organization, such as the MMF, may be necessary for improving the simulation of these convective systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Scalable Multiphysics Block Preconditioning for Low Mach Number Compressible Resistive MHD with Application to Magnetic Confinement Fusion

This study investigates multiphysics block preconditioners that are critical in devising scalable Newton–Krylov iterative solvers for longer time-scale fully implicit fluid plasma models. The specific model of interest is the visco-resistive, low Mach number, compressible magnetohydrodynamics (MHD) model. This model describes the dynamics of conducting fluids in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics. The specific application of interest that motivates this study is the macroscopic simulation of longer time-scale stability and disruptions of magnetic confinement fusion devices, specifically the ITER Tokamak. The computational solution of the governing balance equations for mass, momentum, heat transfer, and magnetic induction for resistive MHD systems can be extremely challenging. These difficulties arise from both the strong nonlinear, nonsymmetric coupling of fluid and electromagnetic phenomena as well as the significant range of time and length scales that the interactions of these physical mechanisms produce. To handle the range of time and spatial scales of interest, a fully implicit unstructured variational multiscale finite element formulation is employed. For the scalable solution of the Newton linearized systems, fully coupled block preconditioners are designed to leverage algebraic multigrid subsolves. In conclusion, results are presented for the strong and weak scaling of the method as well as the robustness of these techniques for a large range of Lundquist numbers.

97 MATHEMATICS AND COMPUTING↗