Results for “Simulation Cases” — Page 3

DOE OSTI2021

3D simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK

In recent years, the nonlinear 3D magnetohydrodynamic codes JOREK, M3D-C1, and NIMROD developed the capability of modeling realistic 3D vertical displacement events (VDEs) including resistive walls. In this paper, a comprehensive 3D VDE benchmark is presented between these state-of-the-art codes. The simulated case is based on an experimental NSTX plasma but with a simplified rectangular wall. There are differences between the physics models and numerical methods, and the VDE evolution leads to sensitivities on the initial conditions that cannot be avoided as can be done in edge localized modes (ELM) and sawtooth simulations (due to the non-cyclical nature of VDEs). Nonetheless, the comparison serves to quantify the level of agreement in the relevant quantities used to characterize disruptions, such as the 3D wall forces and energy decay. The results bring confidence regarding the use of the mentioned codes for disruption studies, and they distinguish aspects that are specific to the models used (e.g., reduced vs full MHD models). Here, the simulations show important 3D fea- tures for a NSTX plasma, such as the self-consistent evolution of the halo current and the origin of the wall forces. In contrast to other reduced MHD models based on an ordering in the aspect ratio, the ansatz-based JOREK reduced MHD model allows capturing many aspects of the 3D dynamics even in the spherical tokamak limit considered here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Dynamic mode decomposition of nonequilibrium electron-phonon dynamics: accelerating the first-principles real-time Boltzmann equation

Abstract Nonequilibrium dynamics governed by electron–phonon ( e -ph) interactions plays a key role in electronic devices and spectroscopies and is central to understanding electronic excitations in materials. The real-time Boltzmann transport equation (rt-BTE) with collision processes computed from first principles can describe the coupled dynamics of electrons and atomic vibrations (phonons). Yet, a bottleneck of these simulations is the calculation of e –ph scattering integrals on dense momentum grids at each time step. Here we show a data-driven approach based on dynamic mode decomposition (DMD) that can accelerate the time propagation of the rt-BTE and identify dominant electronic processes. We apply this approach to two case studies, high-field charge transport and ultrafast excited electron relaxation. In both cases, simulating only a short time window of ~10% of the dynamics suffices to predict the dynamics from initial excitation to steady state using DMD extrapolation. Analysis of the momentum-space modes extracted from DMD sheds light on the microscopic mechanisms governing electron relaxation to a steady state or equilibrium. The combination of accuracy and efficiency makes our DMD-based method a valuable tool for investigating ultrafast dynamics in a wide range of materials.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Vidyut3d: A GPU accelerated fluid solver for non-equilibrium plasmas on adaptive grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure three-electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate ~ 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Vidyut3d: A Gpu Accelerated Fluid Solver for Non-Equilibrium Plasmas on Adaptive Grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure twin electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate approximately 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

Sitaraman, Hariswaran↗

DOE OSTI2023

Fuel Cell Stack Model for Real-Time Simulation of Grid-Connected Applications

Fuel cell stacks coupled with electrolyzers and hydrogen storage sites can be a promising category of distributed energy resources for both grid-connected and stand-alone power systems. However, because of high costs, at-scale hardware testing of fuel cell stacks for grid-connected applications is not economically viable at present. A model-based system that can accurately captures the steady-state and dynamic response of fuel cell stacks over long time periods (hours), is needed. This paper demonstrate a real-time electromagnetic transient model of a megawatt-scale, grid-connected proton exchange membrane fuel cell stack, coupled with a mass-based hydrogen storage system. This model can emulate the electrical steady-state and dynamic response of grid-connected fuel cell stacks. We validate the model using the response of commercial hardware fuel cell stacks and analytical models in the literature - using a digital real-time simulator (RSCAD). The proposed real-time model is then used to simulate cases spanning different time horizons and to design controller-hardware-in-the-loop experiments to evaluate controllers for hydrogen stations.

activation potential↗

DOE OSTI2022

Using Energy Storage-Based Grid Forming Inverters for Operational Reserve in Hybrid Diesel Microgrids

In remote arctic communities, where access to a bulk electrical grid interconnection is not available, the implementation of islanded microgrids has been the most viable way to produce and distribute electricity services to their inhabitants. Historically, these islanded grids have relied primarily on diesel generators or hydropower resources to supply the baseload. However, this practice can result in increased expense due to the high costs associated with fuel transportation and the significant amounts of on-site storage necessary when fuel transportation is unavailable during winter months. In order to mitigate this problem, arctic microgrids have started to transition to a hybrid-source operational mode by incorporating renewable energy sources that are inherently variable in nature, such as wind or solar. Due to their highly stochastic behavior, these hybrid-source islanded microgrids can pose potential issues related to power quality due to introduction of rapid net load fluctuations and inability of diesel generators to respond rapidly. In addition, non-firm stochastic sources can require significant idling diesel generator resources to serve as spinning reserves, which is inefficient and wasteful. This work studies the problems that may arise in the transient dynamics of a real-world hybrid diesel microgrid when subjected to a loss of wind generation. Moreover, this work proposes a transition from a diesel spinning reserve to a battery energy-storage system (BESS) operating reserve scheme. The study of the proposed transition is important in establishing the fundamental implication of transient dynamics and the potential benefits of integrating a BESS as a spinning reserve in terms of stability, frequency nadir, and transient voltage deviation. The methods to investigate and validate the transient dynamics relied on both electromagnetic simulation models of GFMIs and a commercially available GFMI in an experimental power hardware-in-the-loop setup. The simulation results showed that the proposed operating reserve scheme improves the power quality of the system in terms of voltage deviation and frequency nadir when the microgrid is subjected to a loss of wind generation scenario. Depending on the simulation cases, adding a GFMI reduced the frequency nadir between 65.3% and 86.7%. Moreover, the reductions in the voltage deviations improved between 3.6% and 23.0%. From these results it can be concluded that the integration of a GFMI can reduce the frequency nadir in a hybrid diesel microgrid, and in turn, reduce diesel consumption, which improves system reliability while reducing fuel expenses. Furthermore, the novelty of this work relies on the fact that the offline simulation results were validated using a power hardware-in-the-loop platform that incorporated a 100 kVA commercially available GFMI as the device under test.

Hernandez-Alvidrez, Javier↗

DOE OSTI2025

Portable Software Environment for Ultrahigh-Resolution ELM Development on GPUs

This paper presents our endeavors in developing the large-scale, ultra-high-resolution E3SM Land Model (uELM), specifically designed for exascale computers furnished with accelerators such as Nvidia GPUs. The uELM is a sophisticated code that substantially relies on High-Performance Computing (HPC) environments, necessitating particular machine and software configurations. To facilitate community-based uELM developments employing GPUs, we have created a portable, standalone software environment preconfigured with uELM input datasets, simulation cases, and source code. This environment, utilizing Docker, encompasses all essential code, libraries, and system software for uELM development on GPUs. It also features a functional unit test framework and an offline model testbed for comprehensive numerical experiments. From a technical perspective, the paper discusses GPU-ready container generations, uELM code management, and input data distribution across computational platforms. Lastly, the paper demonstrates the use of environment for functional unit testing, end-to-end simulation on CPUs and GPUs, and collaborative code development.

E3SM Land Model↗

DOE OSTI2024

Comprehensive parametric study of CO 2 sequestration in deep saline aquifers

Carbon dioxide injection in deep saline aquifers is a key method for permanently sequestering anthropogenic CO 2 . Here this study employs a reactive transport model to explore mineral precipitation/dissolution and its impact on reservoir properties in deep saline aquifers. We also assess capillary pressure and relative permeability hysteresis on various CO 2 trapping mechanisms. Results from this study reveal the significant influence of initial brine composition on mineral precipitation/dissolution. The dissolution and precipitation of minerals have different effects around the wellbores compared to the overall reservoir. Additionally, salt concentration (Ca ++ and Mg ++ ) and quartz surface area affect CO 2 mineralization, while Na + impacts halite precipitation, altering flow properties. The effect of capillary pressure is significant, as including the capillary pressure in the simulation case resulted in significantly improved CO 2 trapping, achieving almost total dissolution of the injected CO 2 in around 300 years. This study offers novel insights into the interactions of reservoir minerals, brine properties, and the injected CO 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

ITER cold VDEs in the limit of perfectly conducting walls

Recently, it has been shown that a vertical displacement event (VDE) can occur in ITER even when the walls are perfect conductors, as a consequence of the current quench (CQ) [A. H. Boozer, Phys. Plasmas 26, 114501 (2019)]. We used the extended-MHD code M3D-C1 with an ITER-like equilibrium and induced a CQ to explore cold VDEs in the limit of perfectly conducting walls, using different wall geometries. In the case of a rectangular first wall with the side walls far away from the plasma, we obtained very good agreement with the analytical model developed by Boozer that considers a top/bottom flat-plates wall. We show that the solution in which the plasma remains at the initial equilibrium position is improved when bringing the side walls closer to the plasma. When approximating the ITER first wall as a perfect conductor, the plasma remains stable at the initial equilibrium position far beyond the value predicted by the flat-plates wall limit. When considering an opposite limit in which only the inner shell of the ITER vacuum vessel acts as a perfect conductor, the plasma is displaced during the CQ, but the edge safety factor remains above 2 longer in the current decay compared to the flat-plates wall limit. In all the simulated cases, the VDE is found to be strongly dependent on the plasma current, in agreement with a similar finding in the flat-plates wall limit, showing an important difference with hot VDEs in which the CQ is not a necessary condition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Machine Learning Surrogates of a Fuel Matrix Degradation Process Model for Performance Assessment of a Nuclear Waste Repository

Spent nuclear fuel repository simulations are currently not able to incorporate detailed fuel matrix degradation (FMD) process models due to their computational cost, especially when large numbers of waste packages breach. The current paper uses machine learning to develop artificial neural network and k-nearest neighbor regression surrogate models that approximate the detailed FMD process model while being computationally much faster to evaluate. Further, using fuel cask temperature, dose rate, and the environmental concentrations of CO 3 2- , O 2 , Fe 2+ , and H 2 as inputs, these surrogates show good agreement with the FMD process model predictions of the UO 2 degradation rate for conditions within the range of the training data. A demonstration in a full-scale shale repository reference case simulation shows that the incorporation of the surrogate models captures local and temporal environmental effects on fuel degradation rates while retaining good computational efficiency.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2025

An update to the Sandia method for creating Typical Meteorological Years from a limited pool of calendar years

Typical Meteorological Years (TMYs) are essential for the efficient evaluation of energy system performance. Ideally, 30 years of weather data are required to generate TMYs, but significantly fewer years are typically available due to practical limitations. To address this issue, an update to the Sandia method was developed, referred to as the Argonne method, to create TMYs from a limited number of years. Furthermore, this method enhances candidate diversity by systematically shifting original candidate months forward or backward by specific days, creating an expanded pool of candidates. The effectiveness of the Argonne method was validated through statistical testing, comparison of monthly average weather parameters, and numerical simulations. The results demonstrate a high probability of identifying at least one shifted month whose cumulative distribution functions of weather parameters closely align with long-term distributions. In 67 % of all comparisons, the monthly average weather parameters in TMYs generated using the Argonne method exhibit better agreement with long-term averages than TMY3. Moreover, in 74 % of the 318 building simulation cases, the Argonne method outperforms TMY3 in estimating long-term average building heating and cooling demands. Therefore, the Argonne method effectively diversifies the candidate pool and produces typical years that provide more accurate estimations of long-term averages compared to TMY3 when only a limited pool of calendar years (10 years or fewer) is available.

Building energy modeling↗

DOE OSTI2023

The influence of physical and algorithmic factors on simulated far-field waveforms and source–time functions of underground explosions using unsupervised machine learning

SUMMARY Characterizing explosion sources and differentiating between earthquake and underground explosions using distributed seismic networks becomes non-trivial when explosions are detonated in cavities or heterogeneous ground material. Moreover, there is little understanding of how changes in subsurface physical properties affect the far-field waveforms we record and use to infer information about the source. Simulations of underground explosions and the resultant ground motions can be a powerful tool to systematically explore how different subsurface properties affect far-field waveform features, but there are added variables that arise from how we choose to model the explosions that can confound interpretation. To assess how both subsurface properties and algorithmic choices affect the seismic wavefield and the estimated source functions, we ran a series of 2-D axisymmetric non-linear numerical explosion experiments and wave propagation simulations that explore a wide array of parameters. We then inverted the synthetic far-field waveform data using a linear inversion scheme to estimate source–time functions (STFs) for each simulation case. We applied principal component analysis (PCA), an unsupervised machine learning method, to both the far-field waveforms and STFs to identify the most important factors that control variance in the waveform data and differences between cases. For the far-field waveforms, the largest variance occurs in the shallower radial receiver channels in the 0–50 Hz frequency band. For the STFs, both peak amplitude and rise times across different frequencies contribute to the variance. We find that the ground equation of state (i.e. lithology and rheology) and the explosion emplacement conditions (i.e. tamped versus cavity) have the greatest effect on the variance of the far-field waveforms and STFs, with the ground yield strength and fracture pressure being secondary factors. Differences in the PCA results between the far-field waveforms and STFs could possibly be due to near-field non-linearities of the source that are not accounted for in the estimation of STFs and could be associated with yield strength, fracture pressure, cavity radius and cavity shape parameters. Other algorithmic parameters are found to be less important and cause less variance in both the far-field waveforms and STFs, meaning algorithmic choices in how we model explosions are less important, which is encouraging for the further use of explosion simulations to study how physical Earth properties affect seismic waveform features and estimated STFs.

58 GEOSCIENCES↗

DOE OSTI2020

The First Habitable-zone Earth-sized Planet from TESS. III. Climate States and Characterization Prospects for TOI-700 d

We present self-consistent three-dimensional climate simulations of possible habitable states for the newly discovered habitable-zone Earth-sized planet TOI-700 d. We explore a variety of atmospheric compositions, pressures, and rotation states for both ocean-covered and completely desiccated planets in order to assess the planet’s potential for habitability. For all 20 of our simulated cases, we use our climate model outputs to synthesize transmission spectra, combined-light spectra, and integrated broadband phase curves. These climatologically informed observables will help the community assess the technological capabilities necessary for future characterization of this planet—as well as similar transiting planets discovered in the future—and will provide a guide for distinguishing possible climate states if one day we do obtain sensitive spectral observations of a habitable planet around an M star. We find that TOI-700 d is a strong candidate for a habitable world and can potentially maintain temperate surface conditions under a wide variety of atmospheric compositions. Unfortunately, the spectral feature depths from the resulting transmission spectra and the peak flux and variations from our synthesized phase curves for TOI-700 d do not exceed 10 ppm. This will likely prohibit the James Webb Space Telescope from characterizing its atmosphere; however, this motivates the community to invest in future instrumentation that perhaps can one day reveal the true nature of TOI-700 d and to continue to search for similar planets around less distant stars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Large eddy simulation of lean blow-off in swirl-stabilized flame with the subgrid dissipation concept

The lean blow-off mechanism of the premixed swirl flame is numerically investigated by large eddy simulation (LES) with the subgrid dissipation concept (SDC) combustion model. Three simulated cases cover stable, near blow-off, and transient conditions. Compared with the experiment, the LES-SDC approach captures the flow and combustion features for stable and near blow-off conditions. More importantly, the predictions of the blowoff procedure and duration agree satisfactorily with the experiment, indicating that the LES-SDC approach is a promising tool for predicting strong, unsteady turbulent combustion processes. Further, the numerical results are used to investigate the blow-off mechanism. Two stages in the blow-off procedure are specified. The first is the necking and extinction of the downstream flame surface, and the second is the shrinking of the upstream flame surface. The blow-off mechanism is well explained by the theory of stretched flame extinction. At the end of the recirculation zone, the large negative radial velocity pushes the flame to the central line. The combustion process here can be abstracted as the stretched counter-flame of the reactant-reactant configuration. The excessive flame stretch dominates the flame extinction and triggers the blow-off event. The upstream flame resists the intense stretch with the help of hot product recirculation, and the combustion here can be idealized as the counter-flame of the reactant-product configuration. The alignment of the temperature gradient and flow velocity, together with the excessive stretch, clearly indicates the tendency of flame local extinction. A Damköhler number-based blow-off criterion is raised from the mean flow strain rate and laminar flame bulk extinction strain rate.

42 ENGINEERING↗

DOE OSTI2025

Modeling inter‐reader variability in clinical target volume delineation for soft tissue sarcomas using diffusion model

Abstract Background Accurate delineation of the clinical target volume (CTV) is essential in the radiotherapy treatment of soft tissue sarcomas. However, this process is subject to inter‐reader variability due to the need for clinical assessment of risk and extent of potential microscopic spread. This can lead to inconsistencies in treatment planning, potentially impacting treatment outcomes. Most existing automatic CTV delineation methods do not account for this variability and can only generate a single CTV for each case. Purpose This study aims to develop a deep learning‐based technique to generate multiple CTV contours for each case, simulating the inter‐reader variability in the clinical practice. Methods We employed a publicly available dataset consisting of fluorodeoxyglucose positron emission tomography (FDG‐PET), x‐ray computed tomography (CT), and pre‐contrast T1‐weighted magnetic resonance imaging (MRI) scans from 51 patients with soft tissue sarcoma, along with an independent validation set containing five additional patients. An experienced reader drew a contour of the gross tumor volume (GTV) for each patient based on multi‐modality images. Subsequently, two additional readers, together with the first one, were responsible for contouring three CTVs in total based on the GTV. We developed a diffusion model‐based deep learning method that is capable of generating arbitrary number of different and plausible CTVs to mimic the inter‐reader variability in CTV delineation. The proposed model incorporates a separate encoder to extract features from the GTV masks, leveraging the critical role of GTV information in accurate CTV delineation. Results The proposed diffusion model demonstrated superior performance with the highest Dice Index (0.902 compared to values below 0.881 for state‐of‐the‐art models) and the best generalized energy distance (GED) (0.209 compared to values exceeding 0.221 for state‐of‐the‐art models). It also achieved the second‐highest recall and precision metrics among the compared ambiguous image segmentation models. Results from both datasets exhibited consistent trends, reinforcing the reliability of our findings. Additionally, ablation studies exploring different model structures and input configurations highlighted the significance of incorporating prior GTV information for accurate CTV delineation. Conclusions The proposed diffusion model successfully generates multiple plausible CTV contours for soft tissue sarcomas, effectively capturing inter‐reader variability in CTV delineation.

Dong, Yafei [Yale Biomedical Imaging Institute Yal↗

DOE OSTI2020

Protection Scheme for Multiterminal HVDC System Based on Wavelet Transform Modulus Maxima

Fault protection is a major issue in the development of multiterminal high‐voltage direct current (HVDC) system based on voltage‐source converter. Owing to the rapid increase in fault current under DC faults, protection requires a fast fault detection method to isolate the fault current. However, the selectivity and reliability of the fast protection scheme should be primarily guaranteed. Therefore, a new protection strategy based on modulus maximum under wavelet transform is proposed in this paper. Fault detection and discrimination are achieved using single‐end signals of fault current. The external and internal faults can be discriminated via local information. The reliability and robustness of the proposed scheme against potential factors are verified. Numerous simulation cases are performed in PSCAD/EMTDC under different fault scenarios, and signal data are analyzed in MATLAB to demonstrate the effectiveness of the proposed fault protection scheme. © 2020 Institute of Electrical Engineers of Japan. Published by Wiley Periodicals LLC.

Huai, Qing↗