Approximate analytical solution for induction heating of solid cylinders
Not Available
97 MATHEMATICS AND COMPUTING↗
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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Combining direct ink writing with reactive melt infiltration to create architectured thermoelectric legs
Here, we present a new additive-reactive synthesis method where inks – cast into molds or 3D-additively extruded into architectured shapes – are reacted into intermetallic thermoelectric compounds. The new method, as demonstrated for equiatomic TiNiSn, combines: (i) extrusion printing (or casting) of inks containing Ni and Ti powders, (ii) debinding and reactive sintering to form a porous NiTi network, (iii) network infiltration with liquid Sn and subsequent reaction to synthesize the TiNiSn phase. Thin plates, created through this method, show high phase purity and low residual porosity. A thermoelectric figure of merit zT = 0.47 ± 0.05 is achieved at 800 K, within the broad range of values (zT= 0.3–1.0) reported for TiNiSn created via traditional powder hot-pressing methods. Layered TiNiSn microlattices – after 3D ink-extrusion printing and infiltration – exhibit high relative densities and minimal undesirable secondary phase content, with a Seebeck coefficient on par with the ink-cast TiNiSn plates, demonstrating that thermoelectric legs with far-ranging architectural freedom can be created additively with this novel reactive manufacturing method. Numerous other thermoelectric (and other) compounds, currently limited to basic geometries due to brittleness, are amenable to this new method based on a sequence of reactive sintering/infiltration of cast or ink-extruded precursors.
36 MATERIALS SCIENCE↗
An overlapping Domain Decomposition preconditioning method for monolithic solution of shear bands
Metals subjected to high strain rate impact often exhibit a sudden and profound drop in the material’s load bearing capability, a ductile failure phenomenon known as shear banding. Shear bands, characterized as material instabilities, are driven by shear heating and described as narrow regions that have sustained intense plastic deformation and high temperature rise. This coupled thermo-mechanical localization problem can be formulated as a nonlinear system with two balance equations, momentum and energy, and two constitutive laws for elasticity and plasticity. Here in our formulation, mixed finite elements are used to discretize the equations in space and an implicit finite difference scheme is used to advance the system in time, where at every time step, a Newton type method is used to solve the system monolithically. To that end, a block Jacobian matrix is formed analytically using Gâteaux derivatives. The resulting Jacobian is sparse, nonsymmetric and its sparsity pattern and eigenvalue content vary with the different stages of shear bands formation, consequently posing a significant challenge to iterative solvers. To address this issue, an overlapping Domain Decomposition preconditioner that takes into account the physics of shear bands and the domain in which it forms, is proposed. The key idea is to resolve the shear band domain accurately while only solving the remaining domain approximately, with selective updates when necessary. The preconditioner is formed on the basis of an additive Schwartz method that is applied to a Schur complement partition of the system. The proposed preconditioner is implemented in parallel and tested on three different benchmark problems as compared against off-the-shelf solvers. We study the effect of h- and k-refinement that are obtained from Isogeometric discretizations of the system, the overlapping strategy and its parallel performance. Excellent performance is demonstrated in serial and parallel on all benchmark examples.
42 ENGINEERING↗
SOMAFOAM: An OpenFOAM based solver for continuum simulations of low-temperature plasmas
Here, we report the development of SOMAFOAM, a finite volume framework for performing continuum simulations of low-temperature plasmas. The primary goal of this work is to discuss the features of SOMAFOAM along with representative results provided as examples for a range of operating conditions and geometries. This includes plasma and plasma–dielectric systems operating in direct current, radio frequency, and microwave regimes from pressures as low as 100 mTorr to atmospheric pressure. The code has several useful features including the ability to run massively parallel simulations using arbitrary geometries, structured/unstructured meshes, choice of models such as drift–diffusion/full-momentum at runtime, and species-dependent timesteps to name a few. The verification/validation studies presented include comparison with previously published continuum simulations (low-pressure direct current and radio frequency plasma), with experiments (Gaseous Electronics Conference Reference Cell and microwave microplasma ignited in a split ring resonator), and previously published kinetic simulations (low-pressure radio frequency plasma). Other examples provided include a direct current atmospheric pressure microplasma bounded by dielectric sidewalls and a helium–nitrogen plasma ignited using a needle electrode facing a dielectric. The performance of the code is also discussed with serial and distributed memory parallel runs demonstrated up to 512 cores. The design and implementation of the code in a modular object-oriented framework allows for easy extension and seamless coupling with other codes and can be expected to play an important role in both academia and industry.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Capturing sub-grid temperature and moisture variations for wildland fire modeling
Not Available
54 ENVIRONMENTAL SCIENCES↗
Multiscale Universal Interface: A concurrent framework for coupling heterogeneous solvers
Graphical abstract: - Abstract: Concurrently coupled numerical simulations using heterogeneous solvers are powerful tools for modeling multiscale phenomena. However, major modifications to existing codes are often required to enable such simulations, posing significant difficulties in practice. In this paper we present a C++ library, i.e. the Multiscale Universal Interface (MUI), which is capable of facilitating the coupling effort for a wide range of multiscale simulations. The library adopts a header-only form with minimal external dependency and hence can be easily dropped into existing codes. A data sampler concept is introduced, combined with a hybrid dynamic/static typing mechanism, to create an easily customizable framework for solver-independent data interpretation. The library integrates MPI MPMD support and an asynchronous communication protocol to handle inter-solver information exchange irrespective of the solvers' own MPI awareness. Template metaprogramming is heavily employed to simultaneously improve runtime performance and code flexibility. We validated the library by solving three different multiscale problems, which also serve to demonstrate the flexibility of the framework in handling heterogeneous models and solvers. In the first example, a Couette flow was simulated using two concurrently coupled Smoothed Particle Hydrodynamics (SPH) simulations of different spatial resolutions. In the second example, we coupled the deterministic SPH method with the stochastic Dissipative Particle Dynamics (DPD) method to study the effect of surface grafting on the hydrodynamics properties on the surface. In the third example, we consider conjugate heat transfer between a solid domain and a fluid domain by coupling the particle-based energy-conserving DPD (eDPD) method with the Finite Element Method (FEM)
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Multiscale modelling of evolving foams
We present a set of multi-scale interlinked algorithms to model the dynamics of evolving foams. These algorithms couple the key effects of macroscopic bubble rearrangement, thin film drainage, and membrane rupture. For each of the mechanisms, we construct consistent and accurate algorithms, and couple them together to work across the wide range of space and time scales that occur in foam dynamics. These algorithms include second order finite difference projection methods for computing incompressible fluid flow on the macroscale, second order finite element methods to solve thin film drainage equations in the lamellae and Plateau borders, multiphase Voronoi Implicit Interface Methods to track interconnected membrane boundaries and capture topological changes, and Lagrangian particle methods for conservative liquid redistribution during rearrangement and rupture. We derive a full set of numerical approximations that are coupled via interface jump conditions and flux boundary conditions, and show convergence for the individual mechanisms. We demonstrate our approach by computing a variety of foam dynamics, including coupled evolution of three-dimensional bubble clusters attached to an anchored membrane and collapse of a foam cluster.
97 MATHEMATICS AND COMPUTING↗
The Virtual Environment for Reactor Applications (VERA): Design and architecture
VERA, the Virtual Environment for Reactor Applications, is the system of physics capabilities being developed and deployed by the Consortium for Advanced Simulation of Light Water Reactors (CASL). CASL was established for the modeling and simulation of commercial nuclear reactors. VERA consists of integrating and interfacing software together with a suite of physics components adapted and/or refactored to simulate relevant physical phenomena in a coupled manner. VERA also includes the software development environment and computational infrastructure needed for these components to be effectively used. We describe the architecture of VERA from both software and numerical perspectives, along with the goals and constraints that drove major design decisions, and their implications. We explain why VERA is an environment rather than a framework or toolkit, why these distinctions are relevant (particularly for coupled physics applications), and provide an overview of results that demonstrate the use of VERA tools for a variety of challenging applications within the nuclear industry.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Ultrafast dynamics of laser-metal interactions in additive manufacturing alloys captured by in situ X-ray imaging
Abstract not provided
Martin, Aiden A.↗
Sensitivity study of parameters important to Molten Salt Reactor Safety
This paper presents a molten salt reactor (MSR) design parameter sensitivity study using a nodal dynamic modelling methodology with explicitly modified point kinetics equation and Mann’s model for heat transfer. Six parameters that can impact MSR safety are evaluated. A MATLAB-Simulink model inspired by Thorcon’s 550MW th MSR is used for parameter evaluations. A safety envelope was formed to encapsulate power, maximum and minimum temperature, and temperature-induced reactivity feedback. The parameters are perturbed by ±30%. The parameters were then ranked by their subsequent impact on the considered safety envelope, which ranks acceptable parameter uncertainty. The model is openly available on GitHub.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Modeling and Characterizing the electron backscatter in a cylindrical anode-based distributed X-ray source
Upcoming advancements in computed tomography architectures warrants the investigation of new X-ray source designs and the impacts that electron backscatter can have on these designs. One such design being investigated is a distributed, cylindrical anode-based X-ray source. For such a distributed X-ray source, we developed a modeling pipeline for simulating electron optics and transport to characterize the quality of the primary X-ray beam and the electron backscatter behavior. We report our results on the energy distributions of the bremsstrahlung spectra; electron backscatter ratio; and spatial, temporal, and energy distributions of backscattered electrons that return to the anode.
46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
Evaluation of coupling approaches for thermomechanical simulations
Many problems of interest, particularly in the nuclear engineering field, involve coupling between the thermal and mechanical response of an engineered system. The strength of the two-way feedback between the thermal and mechanical solution fields can vary significantly depending on the problem. Contact problems exhibit a particularly high degree of two-way feedback between those fields. This paper describes and demonstrates the application of a flexible simulation environment that permits the solution of coupled physics problems using either a tightly coupled approach or a loosely coupled approach. In the tight coupling approach, Newton iterations include the coupling effects between all physics, while in the loosely coupled approach, the individual physics models are solved independently, and fixed-point iterations are performed until the coupled system is converged. These approaches are applied to simple demonstration problems and to realistic nuclear engineering applications. The demonstration problems consist of single and multi-domain thermomechanics with and without thermal and mechanical contact. Simulations of a reactor pressure vessel under pressurized thermal shock conditions and a simulation of light water reactor fuel are also presented. Here, problems that include thermal and mechanical contact, such as the contact between the fuel and cladding in the fuel simulation, exhibit much stronger two-way feedback between the thermal and mechanical solutions, and as a result, are better solved using a tight coupling strategy.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
3D modeling of missing pellet surface defects in BWR fuel
One of the important roles of cladding in light water reactor fuel rods is to prevent the release of fission products. To that end, it is essential that the cladding maintain its integrity under a variety of thermal and mechanical loading conditions. Local geometric irregularities in fuel pellets caused by manufacturing defects known as missing pellet surfaces (MPS) can in some circumstances lead to elevated cladding stresses that are sufficiently high to cause cladding failure. Accurate modeling of these defects can help prevent these types of failures. The BISON nuclear fuel performance code developed at Idaho National Laboratory can be used to simulate the global thermo-mechanical fuel rod behavior, as well as the local response of regions of interest, in either 2D or 3D. In either case, a full set of models to represent the thermal and mechanical properties of the fuel, cladding and plenum gas is employed. A procedure for coupling 2D full-length fuel rod models to detailed 3D models of the region of the rod containing a MPS defect is detailed in this paper. The global and local model each contain appropriate physics and behavior models for nuclear fuel. This procedure is demonstrated on a simulation of a boiling water reactor (BWR) fuel rod containing a pellet with an MPS defect, subjected to a variety of transient events, including a control blade withdrawal and a ramp to high power. The importance of modeling the local defect using a 3D model is highlighted by comparing 3D and 2D representations of the defective pellet region. Finally, parametric studies demonstrate the effects of the choice of gaseous swelling model and of the depth and geometry of the MPS defect on the response of the cladding adjacent to the defect.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Assessment of CASL VERA for BWR analysis and application to accident tolerant SiC/SiC channel box
Application of the Virtual Environment for Reactor Applications (VERA) to BWR analysis is assessed in this paper by comparing results to those calculated using other widely-used modeling tools, namely the U.S. Nuclear Regulatory Commission’s PARCS/PATHS and the Serpent Monte Carlo particle transport code. Additionally, VERA is used to calculate 3-D temperature and fast neutron flux distributions in silicon carbide (SiC) fiber-reinforced, SiC matrix composite (SiC/SiC) BWR channel boxes, which are being studied as an Accident Tolerant Fuel core structural material concept. We found that the code-to-code comparisons were favorable, and the SiC/SiC channel box evaluation demonstrates the many advanced modeling features VERA provides while also highlighting the non-uniformity in fast neutron flux distributions that can play a role in potential SiC/SiC channel box deformation. Traditional BWR analysis tools do not have the calculation fidelity necessary for coupled assessment of flux and temperature gradients in a SiC/SiC channel box. VERA is a state-of-the-art modeling environment that was developed to increase the safety and economic competitiveness of nuclear power through improved modeling accuracy. While VERA has already been deployed in the nuclear industry for PWR applications, the current study is a vital initial step in the extensive development, validation, and verification that VERA must go through to be useful for BWR applications.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Coupled neutronics and species transport simulation of the Molten Salt Reactor Experiment
This paper presents the development of coupling between the molten salt reactor species transport code Mole and the reactor physics code Griffin. Here, in this study, tracking of delayed neutron precursors was investigated in the Molten-Salt Reactor Experiment (MSRE), accounting for changes in fuel flow velocity as a function of position in the primary loop. The neutron transport calculations in Griffin were performed using 11 energy groups, and the species advection calculations in Mole used 6 delayed neutron precursor groups to predict spatial distribution of the neutron flux and neutron precursors in the MSRE. Mole–Griffin was used to calculate $k$ eff and $β$ eff in the reactor as a function of different volumetric flow rates.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Computational modeling of graphite degradation in molten salt reactors: Role of infiltration
Molten salt reactors (MSRs) often employ graphite as a moderator and reflector. An important challenge for deploying graphite in these reactors is that, due to limited experimental data, our understanding of graphite’s structural integrity in molten salt environments remains incomplete. Here, this study addresses heat generation from fuel-bearing salt that has infiltrated open pores in the graphite, driven primarily by pressure differentials. This is one of multiple identified physical and chemical mechanisms through which molten salt could potentially degrade graphite. Thermally driven stresses are quantified using the Molten-Salt Reactor Experiment (MSRE) graphite moderator elements as a case study. Finite element simulations predict stress distributions at varying infiltration levels, indicating that thermal stresses increase with higher infiltration. Rare-event simulations using the parallel subset simulation framework identify the combinations and corresponding ranges of input parameters that lead to stresses above a specified threshold. In particular, combinations involving high infiltration amounts, high power density, and low thermal conductivity tend to induce the highest stresses. Under the inputs and assumptions considered in this work, the magnitudes of the thermally driven stresses are quite low, with a very low likelihood of causing failure due to exceeding the graphite’s tensile strength. Additionally, rare-event simulations were performed for two more scenarios: a scaled-up moderator geometry and a localized hotspot in the original geometry. Both cases resulted in increased susceptibility to failure, though not to a detrimental extent. Furthermore, the combined effects of irradiation and infiltration-induced thermal stresses were evaluated. The results showed that thermal stresses from infiltration were negligible compared to those caused by irradiation. The findings of such a study are inherently component-specific, but the methodology presented here could be used for similar assessments of salt-infiltration effects in other graphite components.
36 - MATERIALS SCIENCE↗
3D probabilistic fracture mechanics / computational fluid dynamics simulation of a reactor pressure vessel under transient conditions
Reactor pressure vessels (RPVs) are safety-critical light-water-reactor components that, under irradiation, experience long-term material degradation in the form of embrittlement. This can increase their susceptibility to fracture under thermal-shock conditions, which could occur during off-normal transients such as loss-of-coolant accidents (LOCAs). During a LOCA, the most severe conditions for the RPV occur when emergency core cooling water is injected through the cold legs into the water-and-steam-filled RPV. The rapid cooling of the downcomer and internal RPV surface causes decreased temperature and elevated thermally driven tensile stresses in the RPV wall. This, combined with long-term material embrittlement, may cause fracture initiation at pre-existing flaws, challenging the integrity of the RPV. Assessing RPV integrity during transients with a large spatial variation in the coolant temperature requires a modeling approach that considers the effects of spatially varying coolant temperature on the fracture probability of a population of flaws distributed throughout the RPV, accounting for spatially varying embrittlement. Here, the present study addresses this need by demonstrating first-of-its-kind coupling of a high-fidelity 3D computational fluid dynamics code with 3D probabilistic fracture mechanics This was accomplished using representative models of a pressurized-water reactor subjected to small- and medium-break LOCA conditions, both of which can result in large spatial temperature variations. While the observed impact of accounting for 3D effects was minimal under the small-break LOCA this study indicates a significant increase in the probability of fracture initiation under the medium-break LOCA when 3D effects are considered, relative to a spatially uniform cooling scenario.
22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
Impact of layer thickness and well orientation on caprock integrity for geologic carbon storage
Economic feasibility of geologic carbon storage demands sustaining large storage rates without damaging caprock seals. Reactivation of pre-existing or newly formed fractures may provide a leakage pathway across caprock layers. In this paper, we apply an equivalent continuum approach within a finite element framework to model the fluid-pressure-induced reactivation of pre-existing fractures within the caprock, during high-rate injection of super-critical CO 2 into a brine-saturated reservoir in a hypothetical system, using realistic geomechanical and fluid properties. We investigate the impact of reservoir to caprock layer thickness, wellbore orientation, and injection rate on overall performance of the system with respect to caprock failure and leakage. We find that vertical wells result in locally higher reservoir pressures relative to horizontal injection wells for the same injection rate, with high pressure inducing caprock leakage along reactivated opening-mode fractures in the caprock. After prolonged injection, leakage along reactivated fractures in the caprock is always higher for vertical than horizontal injection wells. Furthermore, we find that low ratios of reservoir to caprock thickness favor high excess pressure and thus fracture reactivation in the caprock. Finally, injection into thick reservoir units thus lowers the risk associated with CO 2 leakage.
58 GEOSCIENCES↗