Search NASA⌕ Search

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.

At least 1,081 records · Page 60

Native Natural circulation in VERA

This work seeks to introduce a native natural circulation solver to VERA's thermohydraulics capabilities. This work begins with initial investigations into the viability of such a solver by determining the fundamental equations as well as necessary inputs. The method was then implemented in the form of a Python script coupled to CTF which was compared to TRACE. Once viability of the natural circulation solver was established, it was then implemented natively in MPACT for use with previously created simplified thermohydraulics solvers. This new implementation of the native natural circulation solver in MPACT is then tested for single physics problems for both a BWR and PWR single channel model and compared to TRACE where it shows an error of about 30% for BWRs and about 15% for PWRs. The solver is additionally tested for NuScale based models where it demonstrates an error to the reported flow rate of about 10%, resulting in an eigenvalue error of only 138 pcm when comparing fixed flow rate multiphysics runs. Finally, the solver is tested in multiphysics systems for a small PWR assembly, a small BWR assembly, and a full core approximate model of the NuScale design. In these multiphysics runs, it demonstrates no impact on multiphysics convergence stability, as well as negligible impact on runtimes. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Integral Full Core Multi-Physics PWR Benchmark with Measured Data

In recent years, the importance of modeling and simulation has been highlighted extensively in the DOE research portfolio with concrete examples in nuclear engineering with the CASL and NEAMS programs. These research efforts and similar efforts worldwide aim at the development of high-fidelity multi-physics analysis tools for the simulation of current and next-generation nuclear power reactors. Like all analysis tools, verification and validation is essential to guarantee proper functioning of the software and methods employed. The current approach relies mainly on the validation of single physic phenomena (e.g. critical experiment, flow loops, etc.) and there is a lack of relevant multiphysics benchmark measurements that are necessary to validate high-fidelity methods being developed today. This work introduces a new multi-cycle full-core Pressurized Water Reactor (PWR) depletion benchmark based on two operational cycles of a commercial nuclear power plant that provides a detailed description of fuel assemblies, burnable absorbers, in-core fission detectors, core loading and re-loading patterns. This benchmark enables analysts to develop extremely detailed reactor core models that can be used for testing and validation of coupled neutron transport, thermal-hydraulics, and fuel isotopic depletion. The benchmark also provides measured reactor data for Hot Zero Power (HZP) physics tests, boron letdown curves, and three-dimensional in-core flux maps from 58 instrumented assemblies. The benchmark description is now available online and has been used by many groups. However, much work remains to be done on the quantification of uncertainties and modeling sensitivities. This work aims to address these deficiencies and make this benchmark a true non-proprietary international benchmark for the validation of high-fidelity tools. This report details the BEAVRS uncertainty quantification for the first two cycle of operations and serves as the final report of the project.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Solving macroscopic and microscopic pin-fin heat sink problems by adapted spectral method

In this paper, simulation of heat transfer in a heat sink with macroscopic and microscopic scales when one pin-fin is added to the system is investigated by the proposed spectral method. In the microscopic problems, heat transfer model uses dual-phase lag formulas in contrast with macroscopic problems when Fourier law is used to formulate the governing equation. In macroscopic problem, the results are compared with COMSOL multiphysics software results and a good agreement between the results are shown. In microscopic problems, 3D Gaussian heat source is used and boundary conditions obey the Newton law. Comparisons show the efficiency of the current method, while the results are compared with existed literature. It is shown that: by adding one pin fin to the heat sink problem, temperature of the system will decrease. Numerical results are given and graphs of the temperature for various kinds of pin-fin heat sink problems are depicted.

97 MATHEMATICS AND COMPUTING↗

Bypass Flow Model Implementation for VERA BWR

The Virtual Environment for Reactor Applications (VERA) has been developed by the Consortium for Advanced Simulation of LWRs (CASL) over the past 10 years to address difficult problems facing commercial light-water reactors (LWRs). VERA comprises several high-fidelity physics packages, each of which models a different aspect of the reactor. These individual pieces are then linked to provide a fully coupled multiphysics solution for various LWR analyses.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Advances in Geologic Disposal Safety Assessment and an Unsaturated Alluvium Reference Case

The Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Depat ment of Energy (DOE) Office of Nuclear Energy (NE), Office of Fuel Cycle Technology (OFCT) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and high level nuclear waste (HLW). Two high priorities for SFWST disposal R&D are design concept development and disposal system modeling (DOE 2011, Table 6). These priorities are directly addressed in the SFWST Geologic Disposal Safety Assessment (GDSA) work package, which is charged with developing a disposal system modeling and analysis capability for evaluating disposal system performance for nuclear waste in geologic media. This report describes specific GDSA activities in fiscal year 2018 (FY 2018) toward the development of GDSA Framework, an enhanced disposal system modeling and analysis capability for geologic disposal of nuclear waste. GDSA Framework employs the PFLOTRAN thermal-hydrologic-chemical multiphysics code (Hammond et al. 2011a; Lichtner and Hammond 2012) and the Dakota uncertainty sampling and propagation code (Adams et al. 2012; Adams et al. 2013). Each code is designed for massivelyparallel processing in a high-performance computing (HPC) environment. Multi-physics representations in PFLOTRAN are used to simulate various coupled processes including heat flow, fluid flow, waste dissolution, radionuclide release, radionuclide decay and ingrowth, precipitation and dissolution of secondary phases, and radionuclide transport through engineered barriers and natural geologic barriers to the biosphere. Dakota is used to generate sets of representative realizations and to analyze parameter sensitivity.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Investigation of Ring-heterogeneous Geometry Approximation for Efficient and Practical SFR Analysis

A novel geometry approximation scheme named ring-heterogeneous approximation for intermediate resolution multiphysics analysis of SFRs is presented. This approach particularly targets neutronics coupled with thermo-mechanics feedback, which dominate the passive safety of SFR designs. The ring-heterogeneous approximation casts an array of pins into volume-preserving bands of materials that constitute the pins. It allows to efficiently simulate SFRs with an unstructured mesh and also to explicitly consider the deformation of each material by retaining their heterogeneity. Verification results with single assembly and core problems of the ABTR benchmark demonstrated that the ring-heterogeneous approximation is a reasonable approximation for SFRs. For a 2D full core problem, the ring-heterogeneous approximation significantly reduced the number of elements by at least a factor of 70 from a fully-heterogeneous mesh while almost exactly preserving the reference fully-heterogeneous solutions with the RMS error of only 0.02% in the ring-wise fission reaction rate.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Computational multiphysics modeling of radioactive aerosol deposition in diverse human respiratory tract geometries

The evaluation of aerosol exposure relies on generic mathematical models that assume uniform particle deposition profiles over the human respiratory tract and do not account for subject-specific characteristics. Here we introduce a hybrid-automated computational workflow that generates personalized particle deposition profiles in 3D reconstructed human airways from computed tomography scans using Computational Fluid and Particle Dynamics simulations. This is the first large-scale study to consider realistic airways variability, where 380 lower and 40 upper human respiratory tract 3D geometries are reconstructed and parameterized. The data is clustered into nine groups using random forest regression. Computational fluid and particle dynamics simulations are conducted on these representative geometries using a realistic heavy-breathing respiratory cycle and radioactive iodine-131 as a source term. Monte Carlo radiation transport simulations are performed to obtain detailed energy deposition maps. Our findings emphasize the importance of personalized studies, as minor respiratory tract variations notably influence deposition patterns rather than global parameters of the lower airways, observing more than 30% variance in the mass deposition fraction.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Release of a High Temperature Engineering Test Reactor (HTTR) Steady State Multiphysics Model to the Virtual Test Bed

In response to climate change, global governments and private industry have established a common goal of achieving net-zero emissions by 2050 \cite{osti_1865910}. This goal requires a reassessment of current energy demands and production methods. Reducing emissions at an affordable cost while maintaining grid reliability requires a nationwide collaborative effort among government and industry in the United States. Nuclear power is the leading low-carbon electricity generation method. In the past 50 years, the use of nuclear power has reduced carbon dioxide emissions by over 60 gigatons and has played a crucial role in the security of energy supply~\cite{IEA}. In the U.S., nuclear power accounts for 20\% of the electrical supply and provides energy reliably. Advanced reactors will operate at higher temperatures, operate more efficiently, utilize more energy stored within fuel, and reduce the amount of waste produced \cite{osti_1616270}. To face these challenges and goals, the U.S. Department of Energy has created an initiative to focus on the modeling and simulation tools to support future nuclear power plant design, licensing, and operations. The Virtual Test Bed (VTB)~\cite{vtb2023} was launched by the National Reactor Innovation Center (NRIC) in collaboration with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to support the advanced nuclear reactor community. The VTB involves teams from both Idaho National Laboratory and Argonne National Laboratory and aims to provide example models for a broad range of both current and future advanced reactor designs. A feature of the VTB is the automatic testing of these models to ensure continued functionality as simulation tools are further developed. The VTB and the advanced reactor models documented there are important resources for this initiative. This work describes the inclusion of a new model on the VTB---a High Temperature Engineering Test Reactor (HTTR) steady-state model \cite{LABOURE2023109838}.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SiC cladding R&D execution plan

To advance the development of accident-tolerant fuel cladding for light-water reactors, the Advanced Fuels Campaign actively addresses the R&D gap between the specialized technological progress in SiC-based composite materials, led by industrial partners, and the fundamental scientific understanding of material behavior under reactor operating conditions. The laboratory activities are systematically organized based on the Development Strategy document, which outlines critical barriers to the use of SiC composites in nuclear energy applications for the coming decade and beyond. To complement the strategy, this execution plan specifies a detailed year-by-year work package and establishes key milestones for assessing progress in technological advancement. This plan has been formulated in collaboration with principal investigators at national laboratories, stakeholders, and industry partners with the goal of closing technology gaps and readying SiC-based cladding technologies for full-scale commercial testing. The research conducted at the national laboratories is categorized into three primary subjects: (1) addressing challenges associated with the performance of SiC cladding under normal operating environments, (2) identifying the failure limits of SiC cladding under accident conditions, and (3) modeling fuel performance for SiC cladding systems. Comprehensive laboratory activities integrate the development of advanced modeling tools for SiC-based cladding systems with experimental validation via integrated testing of multiphysics phenomena. This document lays out the laboratory plan to support the development of SiC-based technologies, which have the potential for significant operational and safety benefits for the U.S. nuclear industry. This plan aligns with industry developer needs, relevant Nuclear Energy University Programs, and Nuclear Energy Advanced Modeling and Simulation Programs to ensure SiC cladding technologies are matured as efficiently as possible and will be maintained and updated to accommodate the dynamically evolving progress of development, capability needs, and resource availability.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Laura Users Manual: 5.1-41601

This users manual provides in-depth information concerning installation and execution of LAURA, version 5. LAURA is a structured, multi-block, computational aerothermodynamic simulation code. Version 5 represents a major refactoring of the original Fortran 77 LAURA code toward a modular structure afforded by Fortran 95. The refactoring improved usability and maintainability by eliminating the requirement for problem-dependent re-compilations, providing more intuitive distribution of functionality, and simplifying interfaces required for multiphysics coupling. As a result, LAURA now shares gas-physics modules, MPI modules, and other low-level modules with the FUN3D unstructured-grid code. In addition to internal refactoring, several new features and capabilities have been added, e.g., a GNU-standard installation process, parallel load balancing, automatic trajectory point sequencing, free-energy minimization, and coupled ablation and flowfield radiation.

Mazaheri, Alireza↗

LAURA Users Manual: 5.2-43231

This users manual provides in-depth information concerning installation and execution of LAURA, version 5. LAURA is a structured, multi-block, computational aerothermodynamic simulation code. Version 5 represents a major refactoring of the original Fortran 77 LAURA code toward a modular structure afforded by Fortran 95. The refactoring improved usability and maintainability by eliminating the requirement for problem-dependent re-compilations, providing more intuitive distribution of functionality, and simplifying interfaces required for multiphysics coupling. As a result, LAURA now shares gas-physics modules, MPI modules, and other low-level modules with the FUN3D unstructured-grid code. In addition to internal refactoring, several new features and capabilities have been added, e.g., a GNU-standard installation process, parallel load balancing, automatic trajectory point sequencing, free-energy minimization, and coupled ablation and flowfield radiation.

Mazaheri, Alireza↗

Numerical Improvements in Magnetohydrodynamic, Pulsed Power Simulations of Near-Target Plasmas

Magnetohydrodynamic (MHD) simulations of pulsed power experiments frequently result in unphysical runaway heating when extended to low values of mass density. Traditionally, this has been addressed by the use of conductivity floors below which the plasma is given an arbitrarily low conductivity value. Here, a low-density treatment is presented that allows for low-density material to carry current while maintaining a stable temperature. This treatment is implemented in the Ares multiphysics code. Numerical modifications to the conventional vacuum treatment include an energy-conserving density floor, a modified averaging procedure to determine the thermal conductivity at the edge of adjacent zones, and a modified averaging procedure for determining the thermal conductivity in mixed zones that contain material from differing regions. Additionally, anomalous resistivity and Bohm diffusion models are implemented as simplified models for microphysics-induced enhancement of collisional transport. Importantly, the advantage of these various improvements are illustrated in a simple 1-D pulsed power target where the combined changes result in stable temperatures within the lower density regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment (Final Report)

This report documents the FY25 Theory and Simulation Performance Target (TSPT) of developing an integrated modeling framework for fusion reactor design and assessment (FREDA). Over Q1-Q4, new capabilities were developed across both plasma and engineering domains and demonstrated on an example representation of a Compact Advanced Tokamak with a Dual Cooled Lead Lithium blanket. This represents a first-of-a-kind demonstration of coupled core-to-wall-to-engineering for a reactor. Self-consistent CESOL workflows were applied to provide core, pedestal, and SOL prediction; new modules were developed for energetic particle stability (FAR3D) and transport (TGLF-EP) analysis; and boundary plasma modeling (SOLPS-ITER, BOUT++/Hermes-3) was expanded to evaluate wall and divertor heat fluxes and interface with engineering thermal analysis. A parameterized CAD tool, TRACER, was expanded to generate medium-fidelity divertor, blanket, and coil geometries; OpenFOAM and Diablo workflows were applied for first-wall and divertor thermal analyses with helium cooling; and reduced-order models were created for high-mass-flux divertor cooling. Magnet multiphysics capabilities were verified between Elmer, Diablo, and a new MFEM-based solver, and workflows enable stress, thermal, and neutron-fluence analysis of TF coils with neutronics-driven heating. Nuclear and blanket analysis workflows were demonstrated, including tritium breeding, transport, and CFD-informed thermo-mechanical assessment. Preliminary multi-fidelity uncertainty quantification workflows were applied to boundary modeling codes and shown to achieve variance reductions with fewer high-fidelity boundary simulations. Key findings highlight the challenges of resolving the ITEP gap to find suitable balance between wall and divertor loads, neutron heating, and practical limits of PFC cooling. Next step priorities are to develop automated workflows to check boundary code convergence and detachment, implement tighter physics-engineering CAD provenance tracking, and inclusion of plasma-material interface models for SLAG and tungsten cracking behavior. Collectively, these developments establish sophisticated capabilities for predictive, multi-fidelity, whole-device modeling that integrates plasma physics, materials, magnets, and nuclear engineering to guide pathways to viable Fusion Pilot Plant design points.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhancing EnergyPlus capabilities to model dynamic building envelopes using python plugin

Nearly half of the energy consumption in the United States is related to buildings, resulting in an urgent need to develop innovative technologies to improve building energy efficiency. Dynamic building envelopes, comprising switchable insulation and thermal energy storage materials, have been proposed recently as a promising solution to reduce buildings' heating and cooling loads by thermally coupling the indoor environment with the ambient environment when beneficial while decoupling them when outdoor conditions are not favorable. Although various related technologies are still underway, the whole-building energy modeling tools, like EnergyPlus, do not have the capability to simulate the transient and dynamic nature of dynamic envelope materials and components to accurately predict their impact on building energy use. The objective of this study is to formulate a method in EnergyPlus simulation engine to model multilayer envelopes, comprising dynamic building materials with variable thermophysical properties, and discuss the changes made to the program using a Python plugin. Furthermore, the thermal performance of the dynamic envelopes using the proposed method is compared and verified with the results from a well-established commercial code, COMSOL Multiphysics. A parametric assessment is also conducted to evaluate the energy efficiency benefits of dynamic envelopes in a single-family residential building, demonstrating total annual energy savings up to 11.6 %, when a dynamic envelope operates alone, and up to 18.2 % when it is combined with a thin layer of phase change material as a thermal storage medium. Finally, a United States wide energy efficiency assessment is presented to showcase the geographical spread of the energy savings. The method designed and implemented in this study provides the researchers with the ability to implement their dynamic insulation methods in EnergyPlus and evaluate the whole building energy impact.

25 ENERGY STORAGE↗

Mitigating Injection‐Induced Seismicity Along Basement Faults by Extraction: Application to 2016–2018 Pohang Earthquakes

Abstract The risk of induced earthquake is considered as the major obstacle for sustainable subsurface energy activities. As a mitigation approach, mechanistic risk assessment can be essential to minimize the potential of seismic hazard. Coupled hydro‐mechanical modeling approach is integrated with geomechanical analysis to better assess the occurrence of seismicity induced by injection and/or extraction, which will provide a basis to update operational conditions (e.g., number and location of mitigating wells, rate and duration of operations) for mitigating earthquake nucleation. The Coulomb stability analysis and seismicity rate estimates from this generic study indicate that simultaneous extraction can relieve injection‐induced pore‐pressure buildup, but also transfer and accumulate elastic energy along the faults. Field‐scale simulation with variation in well design for the Pohang EGS site, South Korea, show that multiple wells aligned across a preexisting fault form gradients of pressure and stress fields normal to the fault, which can enhance the direct impact of diffusion process and poroelastic stressing on the fault stability perturbation. Our mechanistic study suggests multiphysics‐oriented mitigation strategies that consider (1) possible mechanical deformation, caused by hydro‐mechanical coupling effects, and (2) proper configuration of well operations based on site‐specific geological characteristics.

Chang, K. W.↗

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↗

EBS Task Force: Task 9/FEBEX Modeling Final Report: Thermo-Hydrological Modeling with PFLOTRAN

This report outlines Sandia National Laboratories modeling studies applied to Stage 1 and Stage 2 of the Full-scale Engineered Barriers Experiment in Crystalline Host Rock (FEBEX) in situ test for the SKB EBS Task Force Task 9. The FEBEX test was a full-scale test conducted over ~18 years at the Grimsel, Switzerland Underground Research Laboratory (URL) managed by NAGRA. It involved emplacing simulated waste packages, in the form of welded cylindrical heaters, inside a tunnel in crystalline granitic rock and surrounded by a bentonite barrier and cement plug. Sensors emplaced within the bentonite monitored the wetting-up, heating, and drying out of the bentonite barrier, and the large resulting data set provides an excellent opportunity for validation of multiphysics Thermal-Hydrological (TH), Thermal-Hydrologic-Chemical (THC), and Thermal-Hydrological-Mechanical (THM) modeling approaches for underground nuclear waste storage and the performance of engineered bentonite barriers. The present status of the EBS Task Force is finalizing Task 9, which follows years of modeling studies of the FEBEX test, by many notable modeling teams (Gens et al., 2009; Sanchez et al. 2010; 2012; Samper et al., 2018). These modeling studies generally use two-dimensional axisymmetric meshes, ignoring threedimensional effects, gravity and asymmetric wetting and dry out of the bentonite engineered barrier. This study investigates these effects with use of the PFLOTRAN THC code with massively parallel computational methods in modeling FEBEX Stage 1 and Stage 2 results. The PFLOTRAN numerical code is an open source, state-of-the-art, massively parallel subsurface flow and reactive transport code operating in a high-performance computing environment (Hammond et al., 2014). Section 2 describes the applied partial differential equations describing mass, momentum and energy balance used in this study, considerations derived by assuming phase equilibrium between gas and liquid phases, constitutive equations for granite, cement plug, and bentonite domains, and specific approaches for use inthe PFLOTRAN code. Section 3 describes the geometry, meshing, and model set-up. Section 4 describes modeling results, Section 5 compares modeling results to field testing data, and Section 6 gives conclusions. The Appendix provides detailed information required by the EBSTask Force for final reporting.

42 ENGINEERING↗