The Smokewise simulation framework: Integrating multiple models for strategic wildfire smoke modeling
Explore the source record for details and available documents.
Crockett, Joseph Lafayette [Sandia National Labora↗
SEARCH · Search NASA
Results for “Integrated 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.
INtegrated Fluid & paRticles simulatioN cOde (INF&RNO) v1.0
INF&RNO (INtegrated Fluid & paRticle simulation cOde) is a 2D cylindrical (r-z) code to model the interaction of short laser pulse(s) with an underdense plasma. The code is based on an envelope model for the laser field, and the action of the laser on the plasma is modeled with the time-averaged ponderomotive force approximation. The code features an improved laser envelope solver which enables an accurate description of the laser pulse evolution deep into depletion even at a reasonably low resolution.
Benedetti, Carlo [Lawrence Berkeley National Labor↗
Prediction of transport in the JET DTE2 discharges with TGLF and NEO models using the TGYRO transport code
Abstract The JET Deuterium-Tritium-Experiment Campaign 2 (DTE2) has demonstrated the highest-ever fusion energy production. To forecast the transport dynamics within these discharges, the TGLF and NEO models within the TGYRO transport code were employed. A critical development in this study is the new quasilinear transport model, TGLF-SAT2, specifically designed to resolve discrepancies identified in JET deuterium discharges. This model accurately describes the saturated three-dimensional (3D) fluctuation spectrum, aligning closely with a database of nonlinear CGYRO turbulence simulations, thereby enhancing the predictive accuracy of TGYRO simulations. In validating against the JET DTE2 discharges across two primary operating scenarios, TGYRO effectively predicted the temperature profiles within a broad radial window ( ρ ∼ 0.2–0.85), though with minor ion temperature discrepancies near the core. However, a consistent underprediction of electron density profiles by 20% across the simulation domain was noted, indicating areas for future refinement. To achieve a self-consistent steady-state solution based on the JET DTE2 discharges, an integrated modeling workflow TGYRO-STEP within the OMFIT framework was introduced. This workflow iterates among the core transport, the pedestal pressure and the MHD equilibrium, ultimately yielding a converged solution that significantly reduces dependence on experimental boundary conditions for temperature and density profiles. The integrated simulation results show negligible differences in electron density and temperature profiles compared to standalone TGYRO modeling, while the ion temperature profile is lower due to the updated boundary condition in TGYRO-STEP. The application of the TGYRO-STEP workflow to JET DTE2 discharges serves as a crucial test to validate its robustness and highlights its limitations, providing valuable insights for its potential future application in ITER and Fusion Power Plant deuterium and tritium prediction modeling.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Influence of simple terrain on the spatial variability of a low-level jet and wind farm performance in the AWAKEN field campaign
In wind energy research, scientific challenges are often associated with complex terrain sites, where orography, vegetation, and buildings disrupt flow uniformity. However, even sites characterized as simple terrain can exhibit significant spatial variability in wind speed, particularly during stable boundary layers (SBLs) and low-level jets (LLJs). This study investigates these terrain interactions using both simulations and observations from the American WAKe ExperimeNt (AWAKEN). We employ a multiscale Weather Research and Forecasting (WRF) model simulation, integrating mesoscale forcing in the coarse domains and representing three rows of turbines from the King Plains wind farm as generalized actuator disks (GAD) in the large-eddy simulation (LES) domains. During a nocturnal LLJ event on 3 April 2023, the downstream, wake-affected turbine rows outperformed the upstream, unwaked row by 25 %–51 %. This counterintuitive result arises from terrain-induced streamwise variations in hub-height wind speed of approximately 4 m s −1 over 5 km – equivalent to ∼50 % of the upstream reference speed. This enhancement outweighs the wake-induced reduction in mean wind speed (∼12 %) and global blockage effects reported in the literature (∼1 %–3.4 %). The multiscale simulations capture the intra-farm spatial variability in power performance observed in SCADA data. Terrain-induced vertical displacement of the LLJ, coupled with large wind shear below the jet maximum, drives the substantial streamwise acceleration within the wind farm. These findings underscore the importance of accounting for spatial variability related to terrain, even in simple landscapes, particularly during LLJ conditions. Incorporating such effects into reduced-order modeling frameworks for wind farm design and control could significantly enhance their effectiveness.
17 WIND ENERGY↗
Field Emission Analysis in SRF Cavities for PIP-II Using GEANT4
Field emission (FE) remains a key limitation to achieving optimal performance and reliability in superconducting radiofrequency (SRF) cavities used in accelerator cryomodules. Addressing this challenge requires a detailed understanding of the generation and transport of FE-induced radiation; however, the absence of standardized measurement protocols complicates comparisons across different testing stages and facilities. This work presents initial steps toward a quantitative, systematic framework for analyzing FE-induced radiation using Geant4 simulations. We model four cavity types to be deployed in Fermilab’s Proton Improvement Plan-II (PIP-II) linac - HB650, LB650, SSR2, and SSR1 - to investigate radiation signatures and their dependence on emission conditions. By integrating simulation results with radiation diagnostics, our goal is to identify probable emission sites, enhance diagnostic accuracy, and optimize detector placement. Ultimately, this approach seeks to inform cavity preparation, assembly, and testing procedures, enabling the delivery of FE-free SRF cryomodules for PIP-II.
Wijethunga, Sajini [Fermilab] (ORCID:0000000287659↗
Field Emission Analysis in SRF Cavities for PIP-II Using GEANT4
Field emission (FE) remains a key limitation to achieving optimal performance and reliability in superconducting radiofrequency (SRF) cavities used in accelerator cryomodules. Addressing this challenge requires a detailed understanding of the generation and transport of FE-induced radiation; however, the absence of standardized measurement protocols complicates comparisons across different testing stages and facilities. This work presents initial steps toward a quantitative, systematic framework for analyzing FE-induced radiation using Geant4 simulations. We model four cavity types to be deployed in Fermilab’s Proton Improvement Plan-II (PIP-II) linac - HB650, LB650, SSR2, and SSR1 - to investigate radiation signatures and their dependence on emission conditions. By integrating simulation results with radiation diagnostics, our goal is to identify probable emission sites, enhance diagnostic accuracy, and optimize detector placement. Ultimately, this approach seeks to inform cavity preparation, assembly, and testing procedures, enabling the delivery of FE-free SRF cryomodules for PIP-II.
Wijethunga, Sajini [Fermilab] (ORCID:0000000287659↗
Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)
Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Inertial Confinement Fusion Design Search Using Bayesian Optimization
Inertial confinement fusion (ICF) experiments rely on complex multi-physics simulation codes such as the Lawrence Livermore National Laboratory-developed HYDRA to guide design work. However, these simulations have several dozen tunable parameters and can be computationally expensive. This makes searching the parameter space challenging and time-consuming. Recently developed automated tools utilize Bayesian optimization to search these high-dimensional parameter spaces for optimal designs. The optimization tools run 2D integrated simulations in HYDRA to converge on a design that produces specified scalar or vector outputs. In this paper, we apply the Bayesian optimization tools to two common tuning scenarios. First, we tune simulation inputs to match measurements of a well-characterized experiment at the National Ignition Facility. This type of tuning is commonly performed to compensate for the use of simplified simulation settings (e.g. reduced resolution) or to account for missing physics in the simulations. Second, we search for an ICF simulation design that has a particular radiation drive profile. These optimizations replicate the kinds of tuning researchers routinely perform, but do so with significantly reduced manual effort. This approach demonstrates a powerful and efficient pathway toward autonomous, high-fidelity design optimization for future ICF experiments.
Bayesian optimization↗
Computational Methods for Multi-Physics Simulation of Melting in Steelmaking
Iron and steel production accounts for approximately 8% of global carbon dioxide (CO) emissions. Pathways to decarbonize include replacing fossil fuels in iron ore reduction and electrifying other steelmaking processes. Iron pellets produced by hydrogen, called Hydrogen Direct Reduced Iron (HDRI), have property differences from those produced using conventional DRI processes. These differences may impact melting in electric arc furnaces (EAF) and other downstream processes. The physical properties of iron pellets vary significantly with temperature during heating, complicating predictions of their behavior. In this project, we seek to develop an integrated simulation, including the fluid flow and convective thermal transport around the pellet particle. We also examine conduction and phase changes within the particle as they impact the melting process. We use adaptive mesh refinement (AMR) to resolve both the changing size of the particle and the complex physics of the interaction between the pellet and the surrounding fluid. We base our simulations on the AMReX-incflo module, which allows large-scale Navier-Stokes simulation while resolving the changing particle size during melting. As we advance our numerical tools, we anticipate an improved understanding of the dynamics of HDRI melting, which will, in turn, accelerate the adoption of low-carbon technologies in the steelmaking industry.
AMReX↗
Simulation results for "Integration of Plant and Microbial Oil Processing at Oilcane Biorefineries for More Sustainable Biofuel Production" publication
Simulation results for "Integration of Plant and Microbial Oil Processing at Oilcane Biorefineries for More Sustainable Biofuel Production" publication.
09 BIOMASS FUELS↗
Simulation results for "Integration of Plant and Microbial Oil Processing at Oilcane Biorefineries for More Sustainable Biofuel Production" publication
Simulation results for "Integration of Plant and Microbial Oil Processing at Oilcane Biorefineries for More Sustainable Biofuel Production" publication.
BioSTEAM↗
Mechanism of Tyrosine-Driven Deprotonation in Photosystem II Revealed by Multiscale Simulations
Photosystem II (PSII) drives light-induced water oxidation via stepwise redox transitions of its oxygen-evolving complex (OEC), a Mn 4 Ca cluster advancing through five intermediate S-states (S 0 –S 4 ). The S 2 → S 3 transition involves a redox event in which a Mn ion donates an electron to the redox-active tyrosine YZ, coupled to deprotonation of an OEC-bound water ligand─yet the underlying coupling mechanism remains unresolved. Time-resolved serial femtosecond crystallography (TR-SFX) has revealed transient electron density shifts near the redox-active tyrosine Y Z , interpreted as sequential oxidation and reduction, with reduction initiating ∼1 μs after excitation and substantially progressed by 30 μs. However, this interpretation conflicts with kinetics from photothermal beam deflection (PBD), time-resolved X-ray absorption spectroscopy (TR-XAS), and electron paramagnetic resonance (EPR), which place electron transfer at 190–400 μs and proton transfer around 30 μs. Here, we reconcile these discrepancies using quantum mechanics/molecular mechanics (QM/MM) and molecular dynamics (MD) simulations. We show that oxidation of P680 and Y Z breaks the symmetry of the nearby hydrogen bonds involving water molecule W4, displacing Y Z and replicating the TR-SFX features of Y Z and Q165 observed at 1 μs. This local perturbation propagates through a hydrogen-bond network, transmitting the electrostatic signal from Y Z to the E65-E312 dyad and triggering redox-coupled deprotonation via the Cl1 channel. By 30 μs, the hydrogen-bond symmetry is restored through deprotonation of W2 (or alternatively W1), reproducing the disappearance of TR-SFX density differences around Y Z and Q165 without requiring Y Z reduction. Our proposed mechanism also gives molecular insights into the O6* density, assigning it to water reorganization rather than a discrete Ca-bound hydroxide species. Here, our results reveal a detailed atomistic mechanism linking Y Z oxidation to long-range proton release and suggest a functional role for the nearby Cl – ion in proton transfer. More broadly, this study underscores the importance of hydrogen-bond dynamics in mediating redox-driven proton transport and demonstrates how integrative simulations can resolve mechanistic ambiguities.
Deprotonation↗
Fleet Algorithm Design for Pooled Rideshare: Integrating Human Factors, Simulation, and Optimization
This dissertation explores the study the integration of human factors modeling and rideshare fleet control algorithms. Pooled rideshare is a unique transportation mode offering that allows riders increased flexibility and accessibility over public transportation, and decreased cost relative to personal vehicles or traditional rideshare. Additionally, relative to personal vehicles, pooled rideshare offers reduced costs and options for those with difficulty obtaining transportation. Prior research in the space typically focused on modeling human behavior, or optimizing system performance, but a lack of integration of the concepts leads to unrealistic or underutilized outcomes. To tackle this problem, novel rideshare assignment, and repositioning strategies were designed and implemented in a simulation environment. Through a series of successive studies, improvements to current rideshare processes were identified, and beneficial outcomes for profitability, accessibility, and traffic were explored. Further, improved metrics to assess rideshare performance were designed and analyzed in the context of improved rideshare offerings. This research contributes to the field of transportation by tackling novel but pragmatic approaches to challenges facing the rideshare industry. Through the course of this dissertation, rideshares impacts on users, operators, and even regulators will be explored in detail. The justification behind the use of a simulation environment, a set of simulated regional models for testing, and the focus on realism and deployability is illustrated. The research identifies holes in potential markets for the use of both private, and public rideshare systems.
Paul, Joseph↗
Toward a Generalizable Prediction Model of Molten Salt Mixture Density with Chemistry‐Informed Transfer Learning
Optimally designing applications of molten salts requires knowledge of their thermophysical properties over a wide range of temperatures and compositions. There exist significant gaps in existing databases and this data can be challenging to experimentally measure due to high temperatures, salt corrosivity, and salt hygroscopicity. Existing databases have been used to create Redlich–Kister (RK) models for mixture density showing improved accuracy with respect to ideal mixing assumptions, but these models require subcomponent data measurements for each new system, therefore lacking generality. In order to address generalizability and data sparsity, a transfer learning procedure is proposed to train deep neural networks (DNNs) using a combination of semi‐empirical relationships (RK), data from the thermophysical arm of the molten salt thermal properties database and universal ab initio properties of component mixtures taken from the joint automated repository for various integrated simulations (JARVIS) classical force‐field inspired descriptors database to predict density in molten salts. Herein, it is shown that DNNs predict molten salt density with an r 2 over 0.99 and a mean absolute percentage error under 1%, outperforming alternative methods.
inorganic materials↗
Design and Analysis of Short Period 2G-HTS Undulators
Recent advancements in second-generation high-temperature superconducting (2G-HTS) materials have opened new possibilities for superconducting undulator (SCU) technology, particularly in achieving shorter periods and higher magnetic fields. Here, this paper presents the design and analysis of a short-period, as small as 10 mm, 2G-HTS undulator, focusing on magnetic performance, mechanical force minimization, and quench protection. A force-balanced coil pack design is introduced to significantly reduce mechanical stresses on the conductors, eliminating the need for heavy support structures. Rapid prototyping techniques such as3D printing is employed to refine the design for scalability and manufacturability. Correction schemes for end-field effects are developed using a combination of 2G-HTS tapes and NbTi conductors to optimize field integrals. Simulation results demonstrated that optimized magnetic performance and enhanced mechanical stability, positioning 2G-HTS undulators as promising candidates for storage ring and free-electron based light sources.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Machine learning-accelerated path integral molecular dynamics simulations of reactive organic electrolytes
Hydrogen bonded electrolytes that exhibit accelerated proton transport via sequential reactive hops have drawn interest for their promise in clean energy applications. Molecular dynamics simulations of these electrolytes offer the opportunity to uncover microscopic mechanistic details that could be used to design and tune the properties of candidate electrolyte technologies. However, accurately modeling the proton transfer reactions and transport properties that give rise to high charge conductivites in these electrolytes proves computationally challenging because of the need to perform lengthy condensed phase simulations, treating both the electronic and nuclear degrees of freedom quantum mechanically. In this paper, we demonstrate that such a modeling task can be efficiently achieved with the use of density functional theory (DFT)-trained machine learning potentials (MLP) to accelerate path integral molecular dynamics (PIMD) simulations. We highlight the practical utility of this approach by using it to benchmark how closely PIMD simulations employing different DFT exchange–correlation functionals reproduce the composition-dependent densities, diffusion coefficients, and electrical conductivities of mixtures consisting of imidazole and levulinic acid. Even with the speedup afforded by our MLPs, PIMD simulations remain quite expensive. Furthermore, in order to render PIMD more computationally tractable, we introduce and benchmark the accuracy of a ring polymer contraction approach that leverages a computationally efficient short-range MLP to accelerate our PIMD simulations by an additional factor of four.
Chemical bonding↗
Estimating ionization states and continuum lowering from ab initio path integral Monte Carlo simulations for warm dense hydrogen
Warm dense matter (WDM) is an active field of research, with applications ranging from astrophysics to inertial confinement fusion. Ionization degree and continuum lowering are important quantities to understand how materials behave under these conditions, but can be difficult to diagnose since experimental campaigns are limited and often require model-dependent analysis. This is especially true for hydrogen, which has a comparably low scattering cross section, making high-quality data particularly difficult to obtain. Consequently, building equation of state tables often relies on simulations in combination with untested approximations to extract properties from experiments. Here, we investigate an approach for extracting the ionization potential depression and ionization degree—quantities which are otherwise not directly accessible from the physical model—from first-principles path integral Monte Carlo (PIMC) simulations utilizing a chemical model. In contrast to experimental measurements, where noise and nonequilibrium effects add to the uncertainty of the inferred parameters, PIMC simulations provide a clean signal with well-defined thermodynamic conditions. Comparisons against commonly used models show a qualitative agreement, but we find deviations primarily for the high-density and high-temperature cases. We also demonstrate the decreasing sensitivity of the dynamic structure factor with respect to both ionization and continuum lowering for increasing scattering angles in x-ray Thomson scattering experiments. Our work has important implications for the design of future experiments, but also offers qualitative understanding of structure factors and the imaginary-time correlation function obtained from first-principles quantum Monte Carlo simulations.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Integrating Experiments and Simulations to Reveal Anisotropic Growth Mechanisms and Interfaces of a One-Dimensional Zeolite
Zeolites are nanoporous crystalline materials critical for diverse industrial applications, yet their growth mechanisms are poorly understood. Here, this study presents a novel integrated framework combining experimental synthesis, high-resolution imaging, coarse-grained molecular dynamics simulations, and computer vision to uncover the mechanisms of growth of SSZ-24, a 1D channel zeolite. We demonstrate how synthesis conditions, such as temperature and reactant concentration, govern crystal anisotropy and surface roughness with growth dynamics differing markedly by crystallographic orientation. Along the channels, growth involves minimal energy barriers and rapid nucleation, resulting in rough surfaces. In contrast, growth perpendicular to the channels requires cooperative molecular organization and is highly sensitive to thermodynamic and kinetic conditions, yielding smooth anisotropic surfaces under low driving forces. By simulating transmission electron microscopy (TEM) images, we bridge molecular-scale simulations with experimental observations, identifying distinct growth mechanisms along different crystal planes. This work offers molecular-level insights into zeolite crystallization, advancing the rational design of nanoporous materials. The integration of cross-disciplinary methodologies establishes a transformative framework for optimizing zeolite synthesis, with implications for broader classes of materials.
Bertolazzo, Andressa A. [Univ. of Utah, Salt Lake ↗