Enabling multiphysics simulation of fusion blankets on exascale architecture using OpenFOAM
Explore the source record for details and available documents.
Sircar, Arpan↗
SEARCH · Search NASA
Results for “Multiphysics simulation”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Multiphysics Simulation Of 2nd Generation 238Pu Production Designs Using COMSOL®
Explore the source record for details and available documents.
Hurt, Christopher↗
High-fidelity multiphysics simulation of BWR assembly with coupled TORT-TD/CTF
This paper describes the application of the coupled codes TORT-TD and CTF to the pin-by-pin modeling of a BWR fuel assembly with thermal-hydraulic feedback. TORT-TD, developed at GRS, is a time-dependent three dimensional discrete ordinates code. CTF is the PSU's improved version of the subchannel code COBRA-TF, which uses a two-fluid, three-field model to represent two-phase flow with entrained droplets, and is commonly applied to evaluate LWR safety margins. The coupled codes system TORT-TD/CTF, already applied to several PWR cases involving MOX, was adapted from PWR to BWR applications. The purpose of the research described in this paper is to verify the coupling for modeling two-phase flow at the pin cell level. Using an ATRIUM-10 assembly, the system's steady-state capabilities were tested on two cases: one without control blade insertion and another with partially inserted blades. The influence of the neutron absorber on local axial and radial parameters is presented. The description of an inlet flow reduction transient is an example for the time-dependent capability of the coupled system. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Multiphysics Simulations of CFRP Degradation in Fire
Explore the source record for details and available documents.
Haluza, Rudy Thomas↗
Physics and Components syntax to enable a systems-based approach to multiphysics
Simulations in MOOSE have traditionally used kernel and boundary condition classes to describe the equations. Downstream applications leveraged a system called Actions to define a pre-packaged discretization of the equations they solve. Unfortunately, the Action base class was very limited, and most applications implemented the same concepts in their Actions. This led to an increased maintenance burden and a reduction in coupling opportunities, save for the use of MultiApps which renders each input mostly independent. With the introduction of multi-system capabilities in MOOSE, there is growing interest in defining entire simulations of complex multiphysics systems in a single input file. By introducing a new Physics system, with its dedicated syntax and a new base class providing wide-ranging capabilities, we are now able to define multiple equations in a single input file in a compact and user-friendly way. With new interactions between Physics and the Component system, these equations can be defined on each component of a complex system. In this talk, we will present the capabilities of these new systems, their interactions, and how to define complex systems multiphysics simulations with Physics and Components.
22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
Multiphysics modeling of SOFC performance degradation caused by interface delamination and active layer cracking
In solid oxide device operation, delamination and cracking cause more catastrophic damage to cell life than other degradation mechanisms, leading to more rapid cell failure. The electrochemical performance degradation caused by electrode/electrolyte interface delamination and active layer cracking were investigated through in-house developed three-dimensional multiphysics simulations. Here the multiphysics simulations, fully calibrated against experimental datasets for button cells in a previous study, were extended to quantitatively investigate the effects of electrode/electrolyte interface delamination and active layer cracking. The simulated polarization curves and impedance behavior with various levels of delamination and cracking provided insights of the performance degradation. Furthermore, the local current distributions within the cell were analyzed for three different operating modes: constant voltage, constant current, and constant power. Finally, a practical strategy was proposed to mitigate the performance degradation in the realistic cell operations, extending overall cell lifetime.
08 HYDROGEN↗
SHARP Assembly-Scale Multiphysics Demonstration Simulations
The NEAMS Reactor Product Line effort aims to develop an integrated multiphysics simulation capability for the design and analysis of future generations of nuclear power plants. The Reactor Product Line code suite’s multi-resolution hierarchy is being designed to ultimately span the full range of length and time scales present in relevant reactor design and safety analyses, as well as scale from desktop to petaflop computing platforms. In this report, we describe our efforts to integrate thermal/hydraulics, neutronics, and structural mechanics modeling codes to perform coupled analysis of a representative fast sodium-cooled reactor fuel assembly. Over the past five years, the Reactor Product Line effort has developed high-fidelity single-physics codes for neutron transport modeling, in the PROTEUS code, and Large Eddy Simulation-based thermal/fluid modeling in the Nek5000 code. Both these codes have been exercised on over 100,000 processors of the IBM Blue Gene/P. The Diablo code has been used to perform structural mechanics and thermo-mechanical modeling. MOAB, the Reactor Geometry Generator (RGG), and MeshKit have been developed to generate and manipulate mesh and mesh-based data, in both serial and parallel environments. These tools together form a strong basis on which to build a multi-physics modeling capability. The goal of developing such a tool is to perform multi-physics neutronics, thermal/fluid, and structural mechanics modeling of the components inside a reactor core. While the focus of this report is on modeling a fast sodium-cooled reactor, it is also the goal that this simulation tool be useful for most reactor types. Here we report on the integration of PROTEUS and Nek5000 into the NEAMS framework. Both codes have been modified to read and write data from and to the MOAB library, in both serial and parallel environments. We focused on broad support for as many original physics code options as possible, rather than focusing only on near-term needs. Both codes were demonstrated through this interface on a relatively simple, seven-pin hexagonal fuel assembly. To support coupled analysis, the new multi-physics driver library COUPE was developed; this library provides flexible support for both loose and tight coupling, and can interface flexibly with codes in multiple programming languages. To couple results between the separate meshes used in each physics, the driver, PROTEUS, and Nek5000 all interface with the MBCoupler solution transfer tool provided with MOAB. Together, these tools form the basis of a new, coupled neutronics and thermal/fluid modeling code, that works in both serial and parallel. In a related effort, a driver in the Diablo code was used to access Nek5000-computed results written to a MOAB database, as well as the temperature data used to drive structural mechanics computations.
Tautges, Timothy J.↗
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↗
Design, fabrication, and testing of a water-cooled tungsten-foam high-heat flux plasma-facing module
Here, we present a comprehensive study detailing the design, fabrication and testing a new modular “finger” component, featuring a hexagonal tungsten foam armor tile bonded to the top surface of the water-cooled finger unit. The mechanical, thermal, and particle transport properties of the tungsten foam armor were determined for a wide range of porosity and cell size, and then combined in an equivalent solid material with “effective” properties suitable for large-scale simulations. A series of 3-D multiphysics simulations for coupled thermal and mechanical behavior of the “finger” module are used to optimize the design before fabrication. The simulation results show that, for the expected experimental conditions, a safe range of heat flux and foam porosity is possible for both water-cooling and helium-cooling. Based on these multiphysics simulations, an optimum finger cooling module has been fabricated at Ultramet, Inc. A series of experiments were designed to study the thermomechanical performance of the “finger” unit under steady-state and transient thermal loading at surface heat flux up to 16 MW/m 2 in the High Energy Flux Test facility (HEFTY) at UCLA. Testing conditions included low-power, steady-state plasma exposure at 6 MW/m 2 , followed by severe high heat flux cycles at 16 MW/m 2 . Results of these tests are explained through comparisons with multiphysics simulations. To eliminate the effects of thermomechanical damage carried out on a module on its thermo-mechanical response on a subsequent test, three different yet identically designed and fabricated finger modules were tested. The results show that the current design is capable of cooling the walls of fusion devices at surface heat flux up to 6 MW/m 2 without failure, and that the failure heat flux limit is below 16 MW/m 2 .
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Combining simulations and data with deep learning and uncertainty quantification for advanced energy modeling
A novel and modern framework for energy modeling is developed in this paper with a focus on nuclear energy modeling and simulation. The framework combines multiphysics simulations and real data, with validation by uncertainty quantification tasks and facilitation by machine and deep learning methods. The hybrid framework is built on the basis of a wide range of physical models, real data, mathematical and statistical methods, and artificial intelligence techniques. The framework is demonstrated in different applications, including quantifying uncertainties in computer simulations, multiphysics coupling, analysis of variance using machine learning surrogate models, deep learning of time series phenomena, and propagating parametric uncertainties of nuclear data. The applications demonstrated are oriented to nuclear engineering simulations, even though majority of the methods are applicable to other energy sources (eg, renewable). Efficient utilization of this framework is expected to yield a much better understanding of the physical phenomena analyzed as well as an improvement in the performance of the energy design/model under construction.
42 ENGINEERING↗
Integration of microseismic monitoring data into coupled flow and geomechanical models with ensemble Kalman filter
Abstract Hydraulic stimulation of low‐permeability rocks in enhanced geothermal systems, shale resources, and CO 2 storage aquifers can trigger microseismic events, also known as microearthquakes (MEQs). The distribution of microseismic source locations in the reservoir may reveal important information about the distribution of hydraulic and geomechanical rock properties. In this paper, we present a framework for conditioning heterogeneous rock permeability and geomechanical property distributions on microseismic data. To simulate the multiphysics processes in these systems, we combine a fully coupled flow and geomechanical model with the Mohr‐Coulomb type rock failure criterion. The resulting multiphysics simulation constitutes the forecast model that relates microseismic source locations to reservoir rock properties. We adopt this forward model in an ensemble Kalman filter (EnKF) data assimilation framework to jointly estimate reservoir permeability and geomechanical property distributions from injection‐induced microseismic response measurements. We show that integration of a large number of spatially correlated microseismic data with practical ensemble sizes can lead to severe underestimation of ensemble spread, and eventually ensemble collapse. To mitigate the variance underestimation issue, two low‐rank data representation schemes are presented and discussed. In the first approach, microseismic data are projected onto a low‐dimensional subspace defined by the left singular vectors of the perturbed observation matrix. The second method uses a coarser grid for representing the microseismic data. A series of numerical experiments is presented to evaluate the performance of the proposed methods and to illustrate their applicability for assimilating microseismic data into coupled flow and geomechanical forward models to estimate multiphysics rock properties.
Tarrahi, Mohammadali↗
Predicting PEMFC performance from a volumetric image of catalyst layer structure using pore network modeling
Not Available
Sadeghi, Mohammad Amin↗
On the use of reverse Brownian motion to accelerate hybrid simulations
Multiscale and multiphysics simulations are two rapidly developing fields of scientific computing. Efficient coupling of continuum (deterministic or stochastic) constitutive solvers with their discrete (stochastic, particle-based) counterparts is a common challenge in both kinds of simulations. We focus on interfacial, tightly coupled simulations of diffusion that combine continuum and particle-based solvers. The latter employs the reverse Brownian motion (rBm), a Monte Carlo approach that allows one to enforce inhomogeneous Dirichlet, Neumann, or Robin boundary conditions and is trivially parallelizable. We discuss numerical approaches for improving the accuracy of rBm in the presence of inhomogeneous Neumann boundary conditions and alternative strategies for coupling the rBm solver with its continuum counterpart. Numerical experiments are used to investigate the convergence, stability, and computational efficiency of the proposed hybrid algorithm.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
An approach for coupled-code multiphysics core simulations from a common input
This study describes an approach for coupled-code multiphysics reactor core simulations that is being developed by the Virtual Environment for Reactor Applications (VERA) project in the Consortium for Advanced Simulation of Light-Water Reactors (CASL). In this approach a user creates a single problem description, called the “VERAIn” common input file, to define and setup the desired coupled-code reactor core simulation. A preprocessing step accepts the VERAIn file and generates a set of fully consistent input files for the different physics codes being coupled. The problem is then solved using a single-executable coupled-code simulation tool applicable to the problem, which is built using VERA infrastructure software tools and the set of physics codes required for the problem of interest. The approach is demonstrated by performing an eigenvalue and power distribution calculation of a typical three-dimensional 17 × 17 assembly with thermal–hydraulic and fuel temperature feedback. All neutronics aspects of the problem (cross-section calculation, neutron transport, power release) are solved using the Insilico code suite and are fully coupled to a thermal–hydraulic analysis calculated by the Cobra-TF (CTF) code. The single-executable coupled-code (Insilico-CTF) simulation tool is created using several VERA tools, including LIME (Lightweight Integrating Multiphysics Environment for coupling codes), DTK (Data Transfer Kit), Trilinos, and TriBITS. Parallel calculations are performed on the Titan supercomputer at Oak Ridge National Laboratory using 1156 cores, and a synopsis of the solution results and code performance is presented. Finally, ongoing development of this approach is also briefly described.
97 MATHEMATICS AND COMPUTING↗
Accelerating Thermochemical Equilibrium Calculations for Nuclear Reactor Applications
Thermochemical properties play a key role in modeling and simulation of several key phenomena in nuclear reactors. There has been an increasing interest in incorporating CALPHAD-based formulations in multiphysics simulations including for Molten Salt Reactors where knowledge of phase evolution of the salt and the chemical potentials of various elements are of utmost importance in source term analyses and redox control. However, the size of such simulations is often limited by the high computational cost of full thermodynamic equilibrium calculations. This work discusses the current efforts aimed at accelerating thermochemical equilibrium calculations for multiphysics simulations performed using the open-source finite element / finite volume code Multiphysics Object Oriented Simulation Environment (MOOSE) [1]. While several methods have been proposed for accelerating phase equilibrium calculations [2], most focus on relatively small systems and often rely on a- priori knowledge of the state-space of the system. Nuclear materials, however, are often multi-component systems owing to the evolution of composition under irradiation and an approach based on a-priori mapping of phase diagram is often not enough. This work is aimed at demonstrating an on-the-fly surrogate modeling framework that uses active learning to reduce the number of full equilibrium calculations that must be performed. By combining with efficient coupling approaches, the surrogate framework helps in reducing the computational cost of thermodynamic equilibrium informed multiphysics simulations of nuclear materials. The performance is benchmarked against full coupling with the thermochemistry library Thermochimica [3]. This work uses a machine learning based approach for constructing surrogate models to predict the stable phases in a multicomponent system. The surrogates were constructed using neural networks and Gaussian process classification. In this work, we compare the relative performance of the two methods. We also demonstrate the use of caching previous calculations by interpolating the values from nearest neighbors. References [1] Lindsay, A.D., et al. "2.0 – MOOSE: Enabling massively parallel multiphysics simulation", SoftwareX, 20 (2022): 101202. [2] Roos, W.A. and Zietsman J.H. "Accelerating complex chemical equilibrium calculations – A Review", Calphad, 77 (2022): 102380. [3] Piro, M.H.A., et al. "The thermochemistry library Thermochimica", Computational Materials Science, 67 (2013): 266-272.
36 MATERIALS SCIENCE↗
Multiphysics Object Oriented Simulation Environment
The Multiphysics Object Oriented Simulation Environment (MOOSE) software library developed at Idaho National Laboratory is a tool. MOOSE, like other tools, doesn't actually complete a task. Instead, MOOSE seeks to reduce the effort required to create engineering simulation applications. MOOSE itself is a software library: a blank canvas upon which you write equations and then MOOSE can help you solve them. MOOSE is comparable to a spreadsheet application. A spreadsheet, by itself, doesn't do anything. Only once equations are entered into it will a spreadsheet application compute anything. Such is the same for MOOSE. An engineer or scientist can utilize the equation solvers within MOOSE to solve equations related to their area of study. For instance, a geomechanical scientist can input equations related to water flow in underground reservoirs and MOOSE can solve those equations to give the scientist an idea of how water could move over time. An engineer might input equations related to the forces in steel beams in order to understand the load bearing capacity of a bridge. Because MOOSE is a blank canvas it can be useful in many scientific and engineering pursuits.
Source record↗
Project Title: Demonstration High Temperature Superconducting NonPlanar Stellarator Magnet with Advanced Manufactured Assemblies
This is the final report for the project “Demonstration High Temperature Superconducting Non- Planar Stellarator Magnet with Advanced Manufactured Assemblies”, funded by DOE, and performed by Type One Energy from September, 2020, to March 2024 involving the Fusion Technology Institute at the University of Wisconsin–Madison, the Plasma Science and Fusion Center (PSFC) at the Massachusetts Institute of Technology (MIT) and Commonwealth Fusion Systems (CFS) to design and fabricate the first non-planar HTS (REBCO) coil for a high-field stellarator based on the SPARC tokamak’s VIPER cable concept. Stellarators at high fields make high-temperature superconducting magnets necessary for a compact fusion device. But the asymmetric and non-planar nature of its components, especially the magnets make it difficult for scalable producibility. To address these challenges, two promising technologies have emerged: advanced manufacturing (AM) for the supporting plates for forming the magnets, and high-temperature superconducting (HTS) cables inside the plates. AM has advanced enough to produce stellarator components with the necessary geometric complexity, size, and the precision, leading to potentially significant reductions in production time, cost, and waste. The cost of HTS tape has decreased dramatically, and progress in HTS planar magnet development has reached a point where it can be proposed for application to complex 3D non-planar magnets. The main objective of this project is to develop, demonstrate and pre-commercialize a novel, non-planar HTS coil shape that remains superconducting to achieve production scalable reductions in time and cost and performance. The proposed technology is based on the novel concept of a precision sub-scale HTS nonplanar coil assembly. This project focuses on the design, fabrication, material optimization of cable design, and validation and demonstration of the high current carrying capability of superconducting magnets and their support in a complex 3D shape needed for application to stellarator magnetic plasma confinement. The specific objectives of this research program include: (1) The successful application of metal AM to build a precision sub-scale HTS nonplanar coil, (2) An HTS cable and cross-section design that can conform to the required nonplanar coil shape (bend radii as tight as 10-cm) and remains superconducting at an engineering current density of 1.35 kA/cm 2 at 77 K and 1 tesla at the conductor (5 kA in the cable). To achieve the above challenging goals, we have formed a multidisciplinary research team consisting of members from Type One Energy and UW-Madison, MIT PSFC and CFS with complementary skills and strong facilities. The team worked collaboratively on fundamental and applied research on the following three major technical areas: (1) Design, fabrication, and optimization of non-planar HTS Cable The ultimate goal of the project is to determine if commercial REBCO tapes and additive manufacturing can be used to fabricate high field (≥ 10T) non-planar coils with tight bending radii (≃ 100mm) and with a degradation of the critical current (Ic) smaller than 20% with respect to the expected performance. We started with shorter length cable to evaluate the scalability of the production process and eventually reached multiple turns for higher magnetic fields. Our findings suggest that a stellarator coil system of a relevant size, characterized by its asymmetric and non-planar components, can be fabricated using a formed cable in plate method. This system can be simulated using a large-scale modeling approach. The use of hybrid modeling 3 techniques will be pivotal in reducing the complexity of the model and in assessing expected performance in designs. (2) Modeling and simulation of the non-planar HTS Cable Multiphysics simulations are performed using the commercial software and are carried out in self-field conditions, involving 2D and 3D models and twisted around one slot of twist-pitched VIPER cable. Multiphysics simulations are mainly focused on the coil for the critical current evaluation, the magnetic field map, self-Lorentz forces and mechanical, and magnetothermal behavior and the quench dynamics. The detailed model and prediction of the superconducting performance of a stellarator-relevant demonstration cable from numerical simulations supports the results from the actual testing backing the results. A detailed description and results are provided in the later sections. (3) Design, fabrication, and optimization of support for the non-planar HTS Cable The team developed an additive manufactured (AM) coil positioning plate that formed into the required non-planar geometry (with bend radii as tight as 10-cm) and to acceptable tolerances required for a stellarator magnet: (+0.25-mm from ideal on dimensions of coil positioning plates and up to +1-mm from ideal for position of wound coil). The plate materials is also included in this selection process from fabrication and 3D printing perspective and commensurate with eventual application to a fusion reactor. From the cost effectiveness point of view, the HTS coil and plate has the potential to cost less than that made in conventional methods with less waste (<75% waste) reducing time (<50%) and cost (<50%), especially as the AM field matures. The application of advanced manufacturing in the construction of the support plates will also lead to cost reduction, as the cables can be easily replaced, thereby making the assembly modular. With such high primary cost and time savings, high current densities and magnetic field, the funded R&D work has validated the designs, proven the feasibility, and characterized the performance of the HTS coil and plate assembly, paving the way for a relevant-size stellarator coil system.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Assembling Multiphysics Nuclear Reactor Simulations Using the MOOSE Framework
The Multiphysics Object Oriented Simulation Environment (MOOSE) [1] is an open-source, parallel finite element framework which provides the foundation for many advanced modeling and simulation tools developed under the Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program [2] for the analysis of advanced reactors. The MOOSE framework provides the common foundational capability on which many NEAMS codes for reactor analysis are built. The MOOSE framework also includes several systems to assemble unique workflows and couplingamong MOOSE-based applications. In particular, the MultiApp and Transfer Systems are widely used to assemble different MOOSE-based or MOOSE-wrapped physics applications together to perform loosely or tightly coupled multiphysics simulations. The National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB) [3] hosts publicly available nuclear reactor multiphysics simulation examples which leverage MOOSE’s MultiApp System to meet the modeling needs of different reactor types. The flexibility and robustness of coupling provided by MOOSE permits rapid development of coupled physics models for a wide range of reactor types and events
22 GENERAL STUDIES OF NUCLEAR REACTORS↗