SEARCH · Engineering Papers

Results for “Multiphysics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,081 records · Page 60

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Modeling physical vapor deposition of energetic materials

Morphology and microstructure of organic explosive films formed using physical vapor deposition (PVD) processes strongly depends on local surface temperature during deposition. Currently, there is no accurate means of quantifying the local surface temperature during PVD processes in the deposition chambers. This study focuses on using a multiphysics computational fluid dynamics tool, STARCCM+, to simulate pentaerythritol tetranitrate (PETN) deposition. The PETN vapor and solid phase were simulated using the volume of fluid method and its deposition in the vacuum chamber on spinning silicon wafers was modeled. The model also included the spinning copper cooling block where the wafers are placed along with the chiller operating with forced convection refrigerant. Implicit time-dependent simulations in two- and three-dimensional were performed to derive insights in the governing physics for PETN thin film formation. PETN is deposited at the rate of 14 nm/s at 142.9 °C on a wafer with an initial temperature of 22 °C. The deposition of PETN on the wafers was calculated at an assumed heat transfer coefficient (HTC) of 400 W/m 2 K. This HTC proved to be the most sensitive parameter in determining the local surface temperature during deposition. Previous experimental work found noticeable microstructural changes with 0.5 mm fused silica wafers in place of silicon during the PETN deposition. This work showed that fused silica slows initial wafer cool down and results in ~10 °C difference for the surface temperature at 500 μm PETN film thickness. It was also found that the deposition surface temperature is insensitive to the cooling power of the copper block due to the copper block's very large heat capacity and thermal conductivity relative to the heat input from the PVD process. Future work should incorporate the addition of local stress during PETN deposition. Lastly, based on simulation results, it is also recommended to investigate the impact of wafer surface energy on the PETN microstructure and morphology formation.

36 MATERIALS SCIENCE↗

DOE OSTI2020

O · , H · , and · OH radical etching probability of polystyrene obtained for a radio frequency driven atmospheric pressure plasma jet

Atmospheric pressure plasma jets have great potential for the surface modification of polymers. In this work, the authors report on polystyrene etching by a radio frequency driven atmospheric pressure plasma jet with a focus on the role of H · , O · , and · OH radicals in this process. The absolute flux of H · , O · , and · OH radicals reaching the surface of the polymer was determined by a COMSOL MULTIPHYSICS reacting fluid dynamics model incorporating detailed transport phenomena in the boundary layer near the substrate. The simulated results of H · and · OH densities in the jet effluent were experimentally verified by two-photon absorption laser induced fluorescence and laser induced fluorescence, respectively. The carbon atom removal flux from the polystyrene surface was taken from previously reported measurements using the same plasma source. The authors show that the boundary layer effects in the interfacial region above the substrate can have a significant impact on the calculated etching probabilities. The reaction probability (β) has a significant uncertainty although a variation of 2 orders of magnitude in β leads to uncertainties of approximately 1 order of magnitude variation in the determined etching probability. The etching probability of polystyrene by · OH radicals was confirmed to be at least an order of magnitude larger than the polystyrene etching probability by O · radicals. The authors also confirmed the weak polystyrene etching probability by H · radicals. The model suggests that the presence of a 30 ppm O 2 impurity can lead to the production of · OH radicals in the far effluent of the Ar + 1% H 2 plasma jet close to the substrate at sufficient densities to enable effective etching.

36 MATERIALS SCIENCE↗

DOE OSTI2014

SU-F-19A-11: Retrospective Evaluation of Thermal Coverage by Thermobrachytherapy Seed Arrangements of Clinical LDR Prostate Implants

Purpose: The superiority of a properly-administered combination of radiation therapy and hyperthermia over radiation alone in treatment of human cancers has been demonstrated in multiple studies examining radiobiology, local control, and survival. Unfortunately, hyperthermia is not yet a common modality in oncology practice, due in part to the technical difficulty of heating a deep-seated target volume to sufficient temperature. To address this problem, our group has invented a thermobrachytherapy (TB) seed based on a commonly-used low dose-rate permanent brachytherapy seed for implant in solid tumors. Instead of the tungsten radiographic marker of the standard seed, the TB seed contains one of a self-regulating ferromagnetic alloy. Placement of a patient implanted with such seeds in an oscillating magnetic field generates heat via induction of eddy currents. We present the results of studies of the capability of clinically-realistic TB seed arrangements to adequately heat defined target volumes. Methods: Seed distributions for several past LDR prostate permanent implant brachytherapy patients were reproduced in the finite element analysis software package COMSOL Multiphysics 4.4, with the difference that TB seeds were modelled, rather than the radiation-only seeds actually used for their treatments. The implant geometries were mainly of the modified peripheral loading type; a range of prostatic volumes and blood perfusion rates likely to be seen in a clinical setting were examined. Results: According to the simulations, when distributed to optimize radiation dose, TB seeds also produce sufficient heat to provide thermal coverage of the target given proper selection of the magnetic field strength. However, the thermal distributions may be improved by additional use of hyperthermia-only seeds. Conclusion: A dual-modality seed intended as an alternative to and using the same implantation apparatus and technique as the standard LDR permanent implant seed has been successfully evaluated for its ability to provide sufficient hyperthermia in clinically-realistic implants. This work was partially supported by the National Institutes of Health (NIH) STTR Grant No. R41 CA153631-01A1.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

WE-G-17A-08: Electron Gun Operation for in Line MRI-Linac Configurations: An Assessment of Beam Fidelity and Recovery Techniques for Different SIDs and Magnetic Field Strengths

Purpose: To test the functionality of medical electron guns within the fringe field of a purpose built superconducting MRI magnet, and to test different recovery techniques for a variety of imaging field strengths and SIDs. Methods: Three different electron guns were simulated using Finite Element Modelling; a standard diode gun, a standard triode gun, and a novel diode gun designed to operate within parallel magnetic fields. The approximate working regime of each gun was established by assessing exit current in constant magnetic fields of varying strength and defining ‘working’ as less than 10% change in injection current. Next, the 1.0T MRI magnet was simulated within Comsol Multiphysics. The coil currents in this model were also scaled to produce field strengths of .5, 1, 1.5 and 3T. Various magnetic shield configurations were simulated, varying the SID from 800 to 1300mm. The average magnetic field within the gun region was assessed together with the distortion in the imaging volume - greater than 150uT distortion was considered unacceptable. Results: The conventional guns functioned in fields of less than 7.5mT. Conversely, the redesigned diode required fields greater than .1T to function correctly. Magnetic shielding was feasible for SIDS of greater than 1000mm for field strengths of .5T and 1T, and 1100mm for 1.5 and 3.0T. Beyond these limits shielding resulted in unacceptable MRI distortion. In contrast, the redesigned diode could perform acceptably for SIDs of less than 812, 896, 931, and 974mm for imaging strengths of 0.5, 1.0, 1.5, 3.0T. Conclusions: For in-line MRIlinac configurations where the electron gun is operating in low field regions, shielding is a straight forward option. However, as magnetic field strength increases and the SID is reduced, shielding results in too great a distortion in the MRI and redesigning the electron optics is the preferable solution. The authors would like to acknowledge funding from the National Health and Research Council (AUS), National Institute of Health (NIH), and Cancer Institute NSW.

43 PARTICLE ACCELERATORS↗

DOE OSTI2014

TH-C-12A-08: New Compact 10 MV S-Band Linear Accelerator: 3D Finite-Element Design and Monte Carlo Dose Simulations

Purpose: To design a new compact S-band linac waveguide capable of producing a 10 MV x-ray beam, while maintaining the length (27.5 cm) of current 6 MV waveguides. This will allow higher x-ray energies to be used in our linac-MRI systems with the same footprint. Methods: Finite element software COMSOL Multiphysics was used to design an accelerator cavity matching one published in an experiment breakdown study, to ensure that our modeled cavities do not exceed the threshold electric fields published. This cavity was used as the basis for designing an accelerator waveguide, where each cavity of the full waveguide was tuned to resonate at 2.997 GHz by adjusting the cavity diameter. The RF field solution within the waveguide was calculated, and together with an electron-gun phase space generated using Opera3D/SCALA, were input into electron tracking software PARMELA to compute the electron phase space striking the x-ray target. This target phase space was then used in BEAM Monte Carlo simulations to generate percent depth doses curves for this new linac, which were then used to re-optimize the waveguide geometry. Results: The shunt impedance, Q-factor, and peak-to-mean electric field ratio were matched to those published for the breakdown study to within 0.1% error. After tuning the full waveguide, the peak surface fields are calculated to be 207 MV/m, 13% below the breakdown threshold, and a d-max depth of 2.42 cm, a D10/20 value of 1.59, compared to 2.45 cm and 1.59, respectively, for the simulated Varian 10 MV linac and brehmsstrahlung production efficiency 20% lower than a simulated Varian 10 MV linac. Conclusion: This work demonstrates the design of a functional 27.5 cm waveguide producing 10 MV photons with characteristics similar to a Varian 10 MV linac.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2015

SU-E-J-111: Finite Element-Based Deformable Image Registration of Pleural Cavity for Photodynamic Therapy

Purpose: The pleural volumes will deform during surgery portion of the pleural photodynamic therapy (PDT) of lung cancer when the pleural cavity is opened. This impact the delivered dose when using highly conformal treatment techniques. In this study, a finite element-based (FEM) deformable image registration is used to quantify the anatomical variation between the contours for the pleural cavities obtained in the operating room and those determined from pre-surgery computed tomography (CT) scans. Methods: An infrared camera-based navigation system (NDI) is used during PDT to track the anatomical changes and contour the lung and chest cavity. A series of CTs of the lungs, in the same patient, are also acquired before the surgery. The structure contour of lung and the CTs are processed and contoured in Matlab and MeshLab. Then, the contours are imported into COMSOL Multiphysics 5.0, where the FEM-based deformable image registration is obtained using the deformed mesh - moving mesh (ALE) model. The NDI acquired lung contour is considered as the reference contour, and the CT contour is used as the target one, which will be deformed. Results: The reconstructed three-dimensional contours from both NDI and CT can be converted to COMSOL so that a three-dimensional ALE model can be developed. The contours can be registered using COMSOL ALE moving mesh model, which takes into account the deformation along x, y and z-axes. The deformed contour has good matches to the reference contour after the dynamic matching process. The resulting 3D deformation map can be used to obtain the locations of other critical anatomic structures, e.g., heart, during surgery. Conclusion: Deformable image registration can fuse images acquired by different modalities. It provides insights into the development of phenomenon and variation in normal anatomical structures over time. The initial assessments of three-dimensional registration show good agreement.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

Polarity effects and apparent ion recombination in microionization chambers

Purpose: Microchambers demonstrate anomalous voltage-dependent polarity effects. Existing polarity and ion recombination correction factors do not account for these effects. As a result, many commercial microchamber models do not meet the specification of a reference-class ionization chamber as defined by the American Association of Physicists in Medicine. The purpose of this investigation is to determine the cause of these voltage-dependent polarity effects. Methods: A series of microchamber prototypes were produced to isolate the source of the voltage-dependent polarity effects. Parameters including ionization-chamber collecting-volume size, stem and cable irradiation, chamber assembly, contaminants, high-Z materials, and individual chamber components were investigated. Measurements were performed with electrodes coated with graphite to isolate electrode conductivity. Chamber response was measured as the potential bias of the guard electrode was altered with respect to the collecting electrode, through the integration of additional power supplies. Ionization chamber models were also simulated using COMSOL Multiphysics software to investigate the effect of a potential difference between electrodes on electric field lines and collecting volume definition. Results: Investigations with microchamber prototypes demonstrated that the significant source of the voltage-dependent polarity effects was a potential difference between the guard and collecting electrodes of the chambers. The voltage-dependent polarity effects for each prototype were primarily isolated to either the guard or collecting electrode. Polarity effects were reduced by coating the isolated electrode with a conductive layer of graphite. Polarity effects were increased by introducing a potential difference between the electrodes. COMSOL simulations further demonstrated that for a given potential difference between electrodes, the collecting volume of the chamber changed as the applied voltage was altered, producing voltage-dependent polarity effects in the chamber response. Ionization chamber measurements and COMSOL simulations demonstrated an inverse relationship between the chamber collecting volume size and the severity of voltage-dependent polarity effects on chamber response. The effect of a given potential difference on chamber polarity effects was roughly ten times greater for microchambers as compared to Farmer-type chambers. Stem and cable irradiations, chamber assembly, contaminants, and high-Z materials were not found to be a significant source of the voltage-dependent polarity effects. Conclusions: A potential difference between the guard and collecting electrodes was found to be the primary source of the voltage-dependent polarity effects demonstrated by microchambers. For a given potential difference between electrodes, the relative change in the collecting volume is smaller for larger-volume chambers, illustrating why these polarity effects are not seen in larger-volume chambers with similar guard and collecting electrode designs. Thus, for small-volume chambers, it is necessary to reduce the potential difference between the guard and collecting electrodes in order to reduce polarity effects for reference dosimetry measurements.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-F-I-50: Finite Element-Based Deformable Image Registration of Lung and Heart

Purpose: Photodynamic therapy (PDT) is used after surgical resection to treat the microscopic disease for malignant pleural mesothelioma and to increase survival rates. Although accurate light delivery is imperative to PDT efficacy, the deformation of the pleural volume during the surgery impacts the delivered light dose. To facilitate treatment planning, we use a finite-element-based (FEM) deformable image registration to quantify the anatomical variation of lung and heart volumes between CT pre-(or post-) surgery and surface contours obtained during PDT using an infrared camera-based navigation system (NDI). Methods: NDI is used during PDT to obtain the information of the cumulative light fluence on every cavity surface point that is being treated. A wand, comprised of a modified endotrachial tube filled with Intralipid and an optical fiber inside the tube, is used to deliver the light during PDT. The position of the treatment is tracked using an attachment with nine reflective passive markers that are seen by the NDI system. Then, the position points are plotted as three-dimensional volume of the pleural cavity using Matlab and Meshlab. A series of computed tomography (CT) scans of the lungs and heart, in the same patient, are also acquired before and after the surgery. The NDI and CT contours are imported into COMSOL Multiphysics, where the FEM-based deformable image registration is obtained. The NDI and CT contours acquired during and post-PDT are considered as the reference, and the Pre-PDT CT contours are used as the target, which will be deformed. Results: Anatomical variation of the lung and heart volumes, taken at different times from different imaging devices, was determined by using our model. The resulting three-dimensional deformation map along x, y and z-axes was obtained. Conclusion: Our model fuses images acquired by different modalities and provides insights into the variation in anatomical structures over time.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

Volcanic arc rigidity variations illuminated by coseismic deformation of the 2011 Tohoku-oki M9

Rock strength has long been linked to lithospheric deformation and seismicity. However, independent constraints on the related elastic heterogeneity are missing, yet could provide key information for solid Earth dynamics. Using coseismic Global Navigation Satellite Systems (GNSS) data for the 2011 M9 Tohoku-oki earthquake in Japan, we apply an inverse method to infer elastic structure and fault slip simultaneously. We find compliant material beneath the volcanic arc and in the mantle wedge within the partial melt generation zone inferred to lie above ~100 km slab depth. We also identify low-rigidity material closer to the trench matching seismicity patterns, likely associated with accretionary wedge structure. Along with traditional seismic and electromagnetic methods, our approach opens up avenues for multiphysics inversions. Those have the potential to advance earthquake and volcano science, and in particular once expanded to InSAR type constraints, may lead to a better understanding of transient lithospheric deformation across scales.

58 GEOSCIENCES↗

DOE OSTI2023

Machine learning–aided real-time detection of keyhole pore generation in laser powder bed fusion

Porosity defects are currently a major factor that hinders the widespread adoption of laser-based metal additive manufacturing technologies. One common porosity occurs when an unstable vapor depression zone (keyhole) forms because of excess laser energy input. With simultaneous high-speed synchrotron x-ray imaging and thermal imaging, coupled with multiphysics simulations, we discovered two types of keyhole oscillation in laser powder bed fusion of Ti-6Al-4V. Amplifying this understanding with machine learning, we developed an approach for detecting the stochastic keyhole porosity generation events with submillisecond temporal resolution and near-perfect prediction rate. Finally, the highly accurate data labeling enabled by operando x-ray imaging allowed us to demonstrate a facile and practical way to adopt our approach in commercial systems.

42 ENGINEERING↗

DOE OSTI2020

A Scalable Multigrid Reduction Framework for Multiphase Poromechanics of Heterogeneous Media

Simulation of multiphase poromechanics involves solving a multiphysics problem in which multiphase flow and transport are tightly coupled with the porous medium deformation. To capture this dynamic interplay, fully implicit methods, also known as monolithic approaches, are usually preferred. The main bottleneck of a monolithic approach is that it requires solution of large linear systems that result from the discretization and linearization of the governing balance equations. Because such systems are nonsymmetric, indefinite, and highly ill-conditioned, preconditioning is critical for fast convergence. Recently, most efforts in designing efficient preconditioners for multiphase poromechanics have been dominated by physics-based strategies. Current state-of-the-art “black-box” solvers such as algebraic multigrid (AMG) are ineffective because they cannot effectively capture the strong coupling between the mechanics and the flow subproblems, as well as the coupling inherent in the multiphase flow and transport process. In this work, we develop an algebraic framework based on multigrid reduction (MGR) that is suited for tightly coupled systems of PDEs. Using this framework, the decoupling between the equations is done algebraically through defining appropriate interpolation and restriction operators. One can then employ existing solvers for each of the decoupled blocks or design a new solver based on knowledge of the physics. We demonstrate the applicability of our framework when used as a “black-box” solver for multiphase poromechanics. Here, we show that the framework is flexible to accommodate a wide range of scenarios, as well as efficient and scalable for large problems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Landmark-Warped Emulators for Models with Misaligned Functional Response

Many computer models output functional data, and in some cases, these functional data have similar, but misaligned, shape characteristics. In this paper, we introduce a general approach for building emulators for computer models that output misaligned functional data when key values in the functional response (landmarks) can be easily identified. This approach has two main parts: modeling the aligned (using the landmarks) functional data, and modeling the functions that map the misaligned data to the aligned space (warping functions). As the warping functions are required to be monotonic, we give special attention to modeling monotonic functional response data. We discuss how our approach can be easily applied for a variety of typical emulators, such as Gaussian processes, Bayesian multivariate adaptive regression splines, and Bayesian additive regression trees, and how sensitivity analysis can be performed. We demonstrate our approach by building emulators for two applications: (1) a high-energy-density physics computer model used to simulate inertial confinement fusion ignition experiments, where model outputs are highly misaligned, and (2) a multiphysics continuum hydrocode used to simulate high-velocity impact experiments, where model outputs are only slightly misaligned. In case (1) traditional methods cannot be applied, while in (2) they can be applied, but the proposed method performs significantly better.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

When Do Extended Physics-Informed Neural Networks (XPINNs) Improve Generalization?

Physics-informed neural networks (PINNs) have become a popular choice for solving high-dimensional partial differential equations (PDEs) due to their excellent approximation power and generalization ability. Recently, extended PINNs (XPINNs) based on domain decomposition methods have attracted considerable attention due to their effectiveness in modeling multiscale and multiphysics problems and their parallelization. However, theoretical understanding of their convergence and generalization properties remains unexplored. In this study, we take an initial step towards understanding how and when XPINNs outperform PINNs. Specifically, for general multilayer PINNs and XPINNs, we first provide a prior generalization bound via the complexity of the target functions in the PDE problem and a posterior generalization bound via the posterior matrix norms of the networks after optimization. Moreover, based on our bounds, we analyze the conditions under which XPINNs improve generalization. Concretely, our theory shows that the key building block of XPINN, namely, the domain decomposition, introduces a tradeoff for generalization. On the one hand, XPINNs decompose the complex PDE solution into several simple parts, which decreases the complexity needed to learn each part and boosts generalization. On the other hand, decomposition leads to less training data being available in each subdomain, and hence such a model is typically prone to overfitting and may become less generalizable. Empirically, we choose five PDEs to show when XPINNs perform better than, similar to, or worse than PINNs, hence demonstrating and justifying our new theory.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Enabling Hyper-Differential Sensitivity Analysis for Ill-Posed Inverse Problems

Inverse problems constrained by partial differential equations (PDEs) play a critical role in model development and calibration. In many applications, there are multiple uncertain parameters in a model that must be estimated. However, high dimensionality of the parameters and computational complexity of the PDE solves make such problems challenging. A common approach is to reduce the dimension by fixing some parameters (which we will call auxiliary parameters) to a best estimate and use techniques from PDE-constrained optimization to estimate the other parameters. In this article, hyper-differential sensitivity analysis (HDSA) is used to assess the sensitivity of the solution of the PDE-constrained optimization problem to changes in the auxiliary parameters. Foundational assumptions for HDSA require satisfaction of the optimality conditions which are not always practically feasible as a result of ill-posedness in the inverse problem. Here we introduce novel theoretical and computational approaches to justify and enable HDSA for ill-posed inverse problems by projecting the sensitivities on likelihood informed subspaces and defining a posteriori updates. Our proposed framework is demonstrated on a nonlinear multiphysics inverse problem motivated by estimation of spatially heterogeneous material properties in the presence of spatially distributed parametric modeling uncertainties.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Two-Level Sketching Alternating Anderson Acceleration for Complex Physics Applications

We present a novel two-level sketching extension of the Alternating Anderson–Picard (AAP) method for accelerating fixed-point iterations in challenging single- and multiphysics simulations governed by discretized PDEs. Our approach combines a static, physics-based projection that reduces the least-squares (LS) problem to the most informative field (e.g., via Schur-complement insight) with a dynamic, algebraic sketching stage driven by a backward stability analysis under Lipschitz continuity. We introduce inexpensive estimators for stability thresholds and cache-aware randomized selection strategies to balance computational cost against memory access overhead. The resulting algorithm solves reduced LS systems in place, minimizes memory footprints, and seamlessly alternates between low-cost Picard updates and Anderson mixing. Implemented in Julia, our two-level sketching AAP achieves up to 50% time-to-solution reductions compared to standard Anderson acceleration—without degrading convergence rates—on benchmark problems including Stokes, 𝑝-Laplacian, bidomain, and Navier–Stokes formulations at varying problem sizes. These results demonstrate the method’s robustness, scalability, and potential for integration into high-performance scientific computing frameworks. Our implementation is available open source in the AAP.jl library.

Barnafi, Nicolas [University of Chile, Santiago]↗

DOE OSTI2026

Discovery of Probabilistic Dirichlet-to-Neumann Maps on Graphs

Dirichlet-to-Neumann maps enable the coupling of multiphysics simulations across computational subdomains by ensuring continuity of state variables and fluxes at artificial interfaces. We present a novel method for learning Dirichlet-to-Neumann maps on graphs using Gaussian processes, specifically for problems where the data obey a conservation law arising from an underlying partial differential equation. Our approach combines discrete exterior calculus and nonlinear optimal recovery to infer relationships between vertex and edge values. This framework yields data-driven predictions with uncertainty quantification across the entire graph, even when observations are limited to a subset of vertices and edges. By minimizing the reproducing kernel Hilbert space norm while penalizing kernel complexity through maximum likelihood estimation, our method ensures that the resulting surrogate strictly enforces conservation laws without overfitting. We demonstrate our method on two representative applications: subsurface flow in fracture networks and arterial blood flow. Finally, the results demonstrate that the method maintains high accuracy and well-calibrated uncertainty estimates even under severe data scarcity, highlighting its potential for scientific applications where limited data and reliable uncertainty quantification are critical.

Dirichlet-to-Neumann map↗

DOE OSTI2023

ARKODE: A Flexible IVP Solver Infrastructure for One-step Methods

We describe the ARKODE library of one-step time integration methods for ordinary differential equation (ODE) initial-value problems (IVPs). In addition to providing standard explicit and diagonally implicit Runge–Kutta methods, ARKODE supports one-step methods designed to treat additive splittings of the IVP, including implicit-explicit (ImEx) additive Runge–Kutta methods and multirate infinitesimal (MRI) methods. We present the role of ARKODE within the SUNDIALS suite of time integration and nonlinear solver libraries, the core ARKODE infrastructure for utilities common to large classes of one-step methods, as well as its use of “time stepper” modules enabling easy incorporation of novel algorithms into the library. Numerical results show example problems of increasing complexity, highlighting the algorithmic flexibility afforded through this infrastructure, and include a larger multiphysics application leveraging multiple algorithmic features from ARKODE and SUNDIALS.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Streaming Data in HPC Workflows Using ADIOS

The “IO Wall” problem, in which the gap between computation rate and data access rate grows continuously, poses significant problems to scientific workflows which have traditionally relied upon using the filesystem for intermediate storage between workflow stages. One way to avoid this problem in scientific workflows is to stream data directly from producers to consumers and avoiding storage entirely. However, the manner in which this is accomplished is key to both performance and usability. This paper presents the Sustainable Staging Transport, an approach which allows direct streaming between traditional file writers and readers with few application changes. SST is an ADIOS “engine”, accessible via standard ADIOS APIs, and because ADIOS allows engines to be chosen at run-time, many existing file-oriented ADIOS workflows can utilize SST for direct application-to-application communication without any source code changes. This paper describes the design of SST and presents performance results from various applications that use SST, for feeding model training with simulation data with substantially higher bandwidth than the theoretical limits of Frontier’s file system, for strong coupling of separately developed applications for multiphysics multiscale simulation, or for in situ analysis and visualization of data to complete all data processing shortly after the simulation finishes.

Podhorszki, Norbert [ORNL] (ORCID:000000019647542X↗