Search NASA⌕ Search

SEARCH · Search NASA

Results for “multiphysics coupling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 145 records · Page 8

Verification methods for drift–diffusion reaction models for plasma simulations

Abstract Compared to other computational physics areas such as codes for general computational fluid dynamics, the documentation of verification methods for plasma fluid codes remains under developed. Current analytical solutions for plasma are often highly limited in terms of testing highly coupled physics, due to the harsh assumptions needed to derive even simple plasma equations. This work highlights these limitations, suggesting the method of manufactured solutions (MMSs) as a potential option for future verification efforts. To demonstrate the flexibility of MMS in verifying these highly coupled systems, the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework was utilized. Thanks to the MOOSE framework’s robustness and modularity, as well as to its physics module capabilities and ecosystem applications (i.e. Zapdos and the chemical reaction network) developed for plasma physics modeling and simulation, this report lays the groundwork for a structured method of conducting plasma fluid code verification.

DeChant, Corey (ORCID:0000000293779732)↗

New developments and verification of fusion blanket simulation capabilities in the MOOSE framework

Multiphysics modeling capabilities have a crucial role to play in the accelerated deployment of fusion energy. To that end, we developed new multiphysics fusion blanket simulation capabilities in the Multiphysics Object-Oriented Simulation Environment (MOOSE). Firstly, we expanded on the existing capabilities of the previously published work, by coupling 3D tritium transport modeling capabilities using the Tritium Migration Analysis Program, version 8 (TMAP8) to an existing tool including thermal hydraulics, fully three-dimensional (3D) heat transfer, and loosely coupled neutronics analysis. Secondly, we performed a thorough verification of the new capabilities and increased testing code coverage to meet MOOSE’s software quality standards. The MOOSE framework follows a strict software quality assurance plan to be Nuclear Quality Assurance, Level 1 compliant. The new multiphysics fusion blanket simulation capabilities are now held to the same standard. Thirdly, to demonstrate MOOSE’s new fusion blanket modeling capabilities, we performed a fully integrated, multiphysics simulation of a 3D solid ceramic breeder blanket design. This proof-of-concept simulation provides the temperature and tritium distribution across the blanket. In conclusion, the combined efforts towards software quality and the development of multiphysics coupling capabilities provide an effective and reliable framework for modeling solid ceramic fusion blankets using MOOSE.

modeling and simulation↗

Development of a BlueCRAB/MELCOR Framework for Supporting Realistic Mechanistic Source Term Calculations in Microreactors

Efforts are currently underway to deploy microreactor modeling and simulation tools to better support vendors and regulatory authorities in submitting and reviewing licensing applications. In particular, the Nuclear Regulatory Commission is expected to rely on the Comprehensive Reactor Analysis Bundle (BlueCRAB) multiphysics toolset in performing design- and beyond-design-basis accident analyses. In addition, the Nuclear Regulatory Commission has been using the MELCOR code to estimate mechanistic source terms during accidents. As MELCOR relies on isotopic inventory and reactor temperature/power evolution profiles during accident conditions—all of which can theoretically be obtained from BlueCRAB—the ultimate goal of this activity is to establish a common BlueCRAB-MELCOR framework. However, prior to the present research, BlueCRAB had never been used to calculate such quantities of interest at the full-core level. While there are many Monte Carlo (MC) codes capable of computing such quantities of interest, they are unable to readily account for multiphysics feedback. BlueCRAB allows for the coupling of different physics codes together to perform multiphysics-informed calculations. Therefore, the purpose of this fiscal year 2023 work is to investigate the feasibility and challenges of performing such calculations within BlueCRAB so as to generate the data that MELCOR relies on. To demonstrate the methodology, the proposed workflow was applied to a prototypical heat pipe-cooled microreactor model. To predict isotopic concentrations (taking into account the ac- cumulation of fission products during operation), the necessary microscopic cross sections were generated via OpenMC and tabulated with respect to temperature and burnup. Next, a recently developed capability in Griffin (the reactor physics application in BlueCRAB) was used to convert the OpenMC output format into the ISOXML format used by Griffin. A multiphysics microscopic depletion calculation that involved performing a coupled full-core, heterogeneous neutron trans- port and thermal calculation at each depletion step was conducted to deplete the core to end of life (EOL) conditions so as to provide both isotopics and the initial condition for the transient calculation. Following a brief null-transient to verify that the initial condition had been properly restarted and was indeed in thermal equilibrium, a heat pipe failure transient was simulated. Thus, the entire workflow of using BlueCRAB to generate MELCOR inputs, from cross-section generation to producing isotopic inventory and power/temperature evolution profiles during transients, is demonstrated. This report also details the identified gaps in the workflow and how they were (for the most part) addressed. Future work should focus on directly including MEL- COR into the workflow by performing a MELCOR calculation using the BlueCRAB-generated input data. In addition, the heat pipe reactor design should be improved so as to reflect more prototypical burnup characteristics at EOL.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗

Alquimia v1.0: a generic interface to biogeochemical codes – a tool for interoperable development, prototyping and benchmarking for multiphysics simulators

Alquimia v1.0 is a generic interface to geochemical solvers that facilitates development of multiphysics simulators by enabling code coupling, prototyping and benchmarking. The interface enforces the function arguments and their types for setting up, solving, serving up output data and carrying out other common auxiliary tasks while providing a set of structures for data transfer between the multiphysics code driving the simulation and the geochemical solver. Alquimia relies on a single-cell approach that permits operator splitting coupling and parallel computation. We describe the implementation in Alquimia of two widely used open-source codes that perform geochemical calculations: PFLOTRAN and CrunchFlow. We then exemplify its use for the implementation and simulation of reactive transport in porous media by two open-source flow and transport simulators: Amanzi and ParFlow. We also demonstrate its use for the simulation of coupled processes in novel multiphysics applications including the effect of multiphase flow on reaction rates at the pore scale with OpenFOAM, the role of complex biogeochemical processes in land surface models such as the E3SM Land Model (ELM) and the impact of surface–subsurface hydrological interactions on hydrogeochemical export from watersheds with the Advanced Terrestrial Simulator (ATS). These applications make it apparent that the availability of a well-defined yet flexible interface has the potential to improve the software development workflow, freeing up resources to focus on advances in process models and mechanistic understanding of coupled problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantification of electrochemical-mechanical coupling in lithium-ion batteries

Lithium-ion battery safety and durability by nature are dependent on electrochemical and mechanical coupling. Interdisciplinary efforts are required to understand and quantify coupling behaviors. Here we design and conduct mechanically constrained charge and discharge characterizations with efforts supported by multiphysics modeling to unravel the coupling mechanisms of solid-liquid electrode-electrolyte and solid-solid active materials in lithium-ion batteries. We demonstrate that a lithium-ion battery cell under mechanical constraint exhibits a higher voltage during charging and a shorter charging time because of increased electrolyte resistance and decreased diffusivity caused by decreased electrode porosity. The reaction force response of the cell is a combined result of the cell structural response mechanically and lithium-ion intercalation/de-intercalation-induced volume variation electrochemically. Under mechanical constraint, cell capacity is significantly reduced in fast-charge scenarios; however, it can be recovered by a constant-voltage charge protocol. The results highlight the promise of multiphysics approaches to unravel the electrochemical-mechanical coupling mechanisms to direct battery system design and management.

25 ENERGY STORAGE↗

Toward a high-fidelity tritium transport modeling for retention and permeation experiments

We report Tritium Migration Analysis Program version 8 (TMAP8), the latest version of TMAP, was developed within the framework of the Multiphysics Object-Oriented Simulation Environment (MOOSE). Created at Idaho National Laboratory (INL), MOOSE is an open-source, dimension-agnostic, fully coupled, and fully implicit multiphysics platform featuring massively parallel computation capabilities. Using TMAP8, tritium transport in a divertor monoblock was analyzed to elucidate the effects of pulsed operation (up to fifty 1,600 s plasma discharge and cool-down cycles) on the tritium in-vessel inventory source term and ex-vessel release term (i.e., tritium retention and permeation) for safety analysis. With its built-in Message Passing Interface capability, TMAP8 can, in under 2 h, simulate tritium transport in three different layered materials (i.e., tungsten, copper, and copper-chromium-zirconium alloy) in 2D geometry, using a single device/computer with 10 cores. The MOOSE-based TMAP8 code can leverage other MOOSE tools developed under the Nuclear Energy Advanced Modeling and Simulation program to perform tritium and thermal transport in complex geometries and multiphysics environments. And via its massively parallel computation, MOOSE will enable the fusion pilot plant designers to conduct high-fidelity multiphysics modeling for the design of the divertor and blanket systems as well as for the safety analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiphysics modeling of accelerators through code integration

This work aims to improve the ability of particle accelerator researchers to develop high-performance accelerator cavity designs by creating an overall multiphysics framework that integrates and couples existing application codes. This framework will allow accelerator researchers to build multiphysics models that will optimize cavity design, improve understanding of whole-device performance, and reduce the development and fabrication costs of accelerator research. We utilize the open-source VizSchema data standard as an intermediate data structure interface layer to standardize interfaces between individual application codes. VizScema is extensively documented online, and plugins for VizSchema are available for popular visualization packages, including VisIt and ParaView. Currently, the work focuses on coupling the EM field solver COMSOL and the electron gun code MICHELLE to allow COMSOL field-solve results to be seamlessly used by MICHELLE for particle-solve. Later work will extend this integration to include other fields, particles, and thermodynamics simulation codes.

43 PARTICLE ACCELERATORS↗

Data transfers for full core heterogeneous reactor high- fidelity multiphysics studies

Multiphysics simulations for nuclear reactor analysis are usually performed by resorting to operator splitting and fixed point iterations between single-physics solvers. This enables the separate solution of each physics, such as neutronics, fuel performance, and thermal hydraulics, on meshes tailored to the requirements of the respective numerical discretizations of the equations. As the equations are coupled, several fields must be transferred between single-physics solves. Projecting fields between meshes while preserving order of accuracy, conservation properties, and mapping non-overlapping geometries is a complex endeavor. This conference paper will present the transfers as implemented in MOOSE, which can handle arbitrary meshes, arbitrary mappings, conservation of integral quantities, and are made to scale with distributed simulations on both ends of the transfers. Their adequacy for advanced nuclear reactor multiphysics coupling is shown through examples and numerical studies.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

High-fidelity multiphysics load following and accidental transient modeling of microreactors using NEAMS tools: Application of NEAMS codes to perform multiphysics modeling analyses of micro-reactor concepts

The feasibility of modeling microreactors using high-fidelity models with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) tools is investigated in this report. Three overarching questions guided this research: can NEAMS tools readily be applied for high-fidelity multiphysics modeling of different types of transients in microreactor designs; how accurate are the results obtained; and are improvements needed in accuracy or user experience of NEAMS tools, especially considering newly developed capabilities? This work builds upon FY-2022 work, and two microreactor concepts considering heat pipe (HP-MR) and gas-cooled (GC-MR) technologies were further analyzed using high-fidelity multiphysics simulations. The NEAMS tools considered and coupled within the MultiApp environment are Griffin for neutronics, BISON for thermo-mechanics, Sockeye for heat pipe modeling (in HP-MR), SAM for 1D Fluid – 3D solid modeling of coolant channels and system modeling of balance of plant components (in GC-MR), and the SWIFT code for hydrogen redistribution in hydride moderator. The Heat Pipe MicroReactor (HP-MR) concept was further analyzed in FY-2023 to demonstrate the stochastic TRISO failure modeling capability in BISON to check operational limits of the TRISO fuel. A new full-core Gas-Cooled MicroReactor (GC-MR) model was developed based on the initial assembly-model used in Y-2022 and used for steady-state and accidental depressurization transient simulations. Accuracy of the simulations performed was assessed through 1) verification analyses completed on the different physics with code-to-code comparison, and 2) validation of the multiphysics simulations based on modeling of the Kilopower Reactor Using Stirling Technology (KRUSTY) experiment. In FY-2023, the mesh and model of KRUSTY was updated to closely match publicly available data, and the neutronic model was verified and validated against experimental control rod worth measurements. The multiphysics model of KRUSTY was developed and used for steady-state analysis and for modeling reactivity insertion transient. The calculated power increase and stabilization agrees well with experimental data following adjustment in fuel thermal expansion coefficient. As an important component of this project, the ANL team gathered experience with a wide range of NEAMS tools: the MOOSE Mesh System, Griffin, BISON, SWIFT, Sockeye, SAM, Workbench, and the MOOSE MultiApp System, and provided assessment of new capabilities. Noteworthy are the user assessment of the “vapor-only” flow model in Sockeye and development of a multiphysics startup transient in HP-MR unit cell for use as tutorial in Sockeye. The full-core GC-MR model was used for assessment of SAM for balance of plant modeling and for demonstrating the SWIFT code capability for hydrogen redistribution modeling in multiphysics transient analyses. In this process, several bugs/issues were identified and reported to developers. Finally, the assembly GC-MR model developed in FY-2022 coupling Griffin, BISON and SAM through flow blockage and rod ejection transients was published to the National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB). The Heat Pipe MicroReactor (HP-MR) concept high-fidelity multiphysics coupling of Griffin/BISON/Sockeye in load-following and heat pipe failure transients was also published on the VTB. Those submissions are enabling thorough review of these models as well as wide distribution to industry, regulator, and university users. In this analysis, several new research questions were uncovered, and follow-up analyses are recommended to further improve some models, consider additional transients, and continue development of VTB models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Coupling of Spark Plasma Sintering with Advanced Modeling to Enable Process Scale-Up: Presentation to DOE-NE [Slides]

The research goal of this project is to develop at Idaho National Laboratory (INL) a first-of-its-kind Multiphysics Object-Oriented Simulation Environment (MOOSE)-based, multiscale, multiphysics spark plasma sintering (SPS) modeling and simulation code application, termed “Freya.” Freya will simulate the thermo-mechanical-electrical aspects of the SPS fabrication process and will be paired with lower length scale sub-models, such as phase-field, to predict the resulting microstructure. SPS is an advanced manufacturing process that can be used to solve a variety of material manufacturing challenges; however, this process is an extremely challenging problem for modeling and simulation. The SPS process is inherently multiphysics and multi-scale, with the macroscale electro-thermo-mechanical behavior linked intricately to the microstructure evolution of the part being sintered. Accurate modeling and simulation tools, specifically geared towards the SPS process, are needed to predict the influence of the multiple variables involved in the manufacturing process. Modeling and simulation accuracy is achieved and demonstrated through comparison to multiple validation experiments. The validation efforts for Freya include both separate effects and complete multiphysics SPS process experiments. One of the key benefits this Laboratory Directed Research & Development (LDRD) project offers stems from the emphasis placed on experimental validation of the Freya models, both on the individual length scales and of the final coupled multiscale multiphysics simulations. Experimental validation of Freya’s multiscale coupling capability provides the technical credibility necessary for potential future industry and research partners to accept the simulation predictions.

36 MATERIALS SCIENCE↗

Modeling kinetic effects of charged vacancies on electromechanical responses of ferroelectrics: Rayleighian approach

Understanding the time-dependent effects of charged vacancies on the electromechanical responses of materials is at the forefront of research for designing materials exhibiting metal-insulator transitions and memristive behavior. A Rayleighian approach is used to develop a model for studying the nonlinear kinetics of the reaction leading to generation of vacancies and electrons via the dissociation of vacancy-electron pairs. Also, diffusion and elastic effects of charged vacancies are considered to model polarization-electric potential and strain-electric potential hysteresis loops. The model captures multiphysics phenomena by introducing couplings among polarization, the electric potential, stress, strain, and concentrations of charged (multivalent) vacancies and electrons (treated as classical negatively charged particles), where the concentrations can vary due to association-dissociation reactions. A derivation of coupled time-dependent equations based on the Rayleighian approach is presented. Three limiting cases of the governing equations are considered, highlighting the effects of (1) nonlinear reaction kinetics on the generation of charged vacancies and electrons, (2) Vegard's law (i.e., the concentration-dependent local strain) on asymmetric strain-electric potential relations, and (3) coupling between a fast component and the slow component of the net polarization on the polarization-electric-field relations. The Rayleighian approach discussed in this work should pave the way for developing a multiscale modeling framework in a thermodynamically consistent manner while capturing multiphysics phenomena in ferroelectric materials. Published by the American Physical Society 2025

Kumar, Rajeev (ORCID:0000000194943488)↗

Coupled Monte Carlo Transport and Conjugate Heat Transfer for Wire-Wrapped Bundles Within the MOOSE Framework

Cardinal is an open-source application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics (CFD) to the Multiphysics Object-Oriented Simulation Environment (MOOSE), closing neutronics and thermal-fluid gaps in conducting high-resolution multiscale and multiphysics analyses of nuclear systems. We first provide a brief introduction to Cardinal's software design, data mapping, and coupling strategy to highlight our approach to overcoming common challenges in high-fidelity multiphysics simulations. Here we then present two Cardinal simulations for hexagonal pin bundles. The first is a validation of Cardinal's conjugate heat transfer coupling of NekRS's Reynolds-Averaged Navier Stokes model with MOOSE's heat conduction physics for a bare seven-pin Freon-12 bundle flow experiment. Predictions for pin surface temperatures under three different heating modes agree reasonably well with experimental data and similar CFD modeling from the literature. The second simulation is a multiphysics coupling of OpenMC, NekRS, and BISON for a reduced-scale, seven-pin wire-wrapped version of an Advanced Burner Reactor bundle. Wire wraps are approximated using a momentum source model, and coupled predictions are provided for velocity, temperature, and power distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multiphysics Analysis of the MSRE Experiment Using Griffin-SAM Coupled Code System

Molten salt reactors (MSRs) with flowing fuel have the unique feature of utilizing the fuel salt for heat generation and extraction at the same time since the fuel salt is circulating through the whole primary loop of the reactor. This movement of the fuel salt results in a partial decay of the delayed neutron precursors (DNPs) outside the core and corresponding redistribution in the active core region. To capture this phenomenon accurately, neutronics and thermal hydraulics computational tools need to be modified to handle the movement of the DNPs and their decay. Idaho and Argonne National Laboratories are actively working on developing neutronics and thermal hydraulics tools to model MSRs with flowing fuel. The neutronics code Griffin and the thermal hydraulics code SAM were extended to handle flowing fuel with the drift of the DNPs under the MOOSE framework. In this work, a simplified Multiphysics model of the Molten Salt Reactor Experiment (MSRE) is developed and utilized to perform steady-state and transient analyses. The MSRE experiments were designed to show the impact of the DNPs losses on core reactivity. The following section provides a description of the MSRE experiment with the core main parameters used to develop the model. Then, the developed Multiphysics model of the MSRE is presented. Finally, steady-state verification tests and unprotected loss of flow transient tests are discussed followed by a summary and conclusions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A new coupling of a GPU-resident large-eddy simulation code with a multiphysics wind turbine simulation tool

The development of new wind farm control strategies can benefit from combined analysis of flow dynamics in the farm and the behavior of individual turbines within one simulation environment. In this work, we present such an environment by developing a new coupling between the large-eddy simulation (LES) code GRASP and the multiphysics wind turbine simulation tool OpenFAST via an actuator line model (ALM). In addition, the implementation of the recently proposed filtered actuator line model (FALM) within the coupling is described. The new ALM implementation is cross-verified with results from four other commonly used research LES codes. The results for the blade loads and the near wake obtained with the new coupling are consistent with the other codes. Deviations are observed in the far wake. The results further indicate that the FALM is able to reduce the lift and power overprediction from which the traditional ALM suffers on coarse LES grids. This new simulation environment paves the way for future wind farm simulations under realistic weather conditions by leveraging GRASP's ability to impose data from large-scale meteorological models as boundary conditions.

17 WIND ENERGY↗

High-order partitioned spectral deferred correction solvers for multiphysics problems

We present an arbitrarily high-order, conditionally stable, partitioned spectral deferred correction (SDC) method for solving multiphysics problems using a sequence of pre-existing single-physics solvers. This method extends the work in [1], [2], which used implicit-explicit Runge-Kutta methods (IMEX) to build high-order, partitioned multiphysics solvers. We consider a generic multiphysics problem modeled as a system of coupled ordinary differential equations (ODEs), coupled through coupling terms that can depend on the state of each subsystem; therefore the method applies to both a semi-discretized system of partial differential equations (PDEs) or problems naturally modeled as coupled systems of ODEs. The sufficient conditions to build arbitrarily high-order partitioned SDC schemes are derived. Based on these conditions, various of partitioned SDC schemes are designed. The stability of the first-order partitioned SDC scheme is analyzed in detail on a coupled, linear model problem. We show that the scheme is conditionally stable, and under conditions on the coupling strength, the scheme can be unconditionally stable. We demonstrate the performance of the proposed partitioned solvers on several classes of multiphysics problems with moderate coupling strength. They include a stiff linear system of ODEs, advection-diffusion-reaction systems, and fluid-structure interaction problems with both incompressible and compressible flows, where we verify the design order of the SDC schemes and study various stability properties. We also directly compare the accuracy, stability, and cost of the proposed partitioned SDC solver with the partitioned IMEX method in [1], [2] on this suite of test problems. The results suggest that the high-order partitioned SDC solvers are more robust than the partitioned IMEX solvers for the numerical examples considered in this work, while the IMEX methods require fewer implicit solves.

97 MATHEMATICS AND COMPUTING↗

Optical neural engine for solving scientific partial differential equations

Abstract Solving partial differential equations (PDEs) is the cornerstone of scientific research and development. Data-driven machine learning (ML) approaches are emerging to accelerate time-consuming and computation-intensive numerical simulations of PDEs. Although optical systems offer high-throughput and energy-efficient ML hardware, their demonstration for solving PDEs is limited. Here, we present an optical neural engine (ONE) architecture combining diffractive optical neural networks for Fourier space processing and optical crossbar structures for real space processing to solve time-dependent and time-independent PDEs in diverse disciplines, including Darcy flow equation, the magnetostatic Poisson’s equation in demagnetization, the Navier-Stokes equation in incompressible fluid, Maxwell’s equations in nanophotonic metasurfaces, and coupled PDEs in a multiphysics system. We numerically and experimentally demonstrate the capability of the ONE architecture, which not only leverages the advantages of high-performance dual-space processing for outperforming traditional PDE solvers and being comparable with state-of-the-art ML models but also can be implemented using optical computing hardware with unique features of low-energy and highly parallel constant-time processing irrespective of model scales and real-time reconfigurability for tackling multiple tasks with the same architecture. The demonstrated architecture offers a versatile and powerful platform for large-scale scientific and engineering computations.

Tang, Yingheng (ORCID:0009000153622546)↗

Coupling Finite Element and Finite Volume Simulation Within MOOSE

This work focuses on adding capability to couple finite element (FE) variables into finite volume (FV) physics within INL’s Multiphysics Object-Oriented Simulation Environment (MOOSE). This coupling can allow for improvement of multiphysics simulations where one set of physics is best suited for an FE discretization, while another set is best suited for an FV discretization. Electrohydrodynamics, which govern liquid metal reactor concepts and plasma dynamics, is a fitting example where the electromagnetic field equations are solved with FE and the fluid flow is solved with FV. The new FE to FV coupling method can be summarized as taking the element or face average of the FE variable value or gradient and applying that quantity directly in the FV equation objects.

97 MATHEMATICS AND COMPUTING↗