Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation tools”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

SU-G-TeP1-12: Random Repainting as Mitigation for Scanned Ion Beam Interplay Effects

Purpose: Interference of dose application in scanned beam particle therapy and organ motion may lead to interplay effects with distorted dose to target volumes. Interplay effects depend on the speed and direction of the scanning beam, leading to fringed field edges (scanning parallel to organ motion direction) or over- and under-dosed regions (both directions are orthogonal). Current repainting methods can mitigate interplay effects, but are susceptible to artefacts when only a limited number of repaints are applied. In this study a random layered-repainting strategy was investigated. Methods: Mono-energetic proton beams were irradiated to a 10 ×10 cm{sup 2} scanned field at a Varian ProBeam facility. Applied dose was measured with a 2D amorphous silicon detector mounted on a motion platform (CIRS dynamic platform). Motion was considered with different cycles, directions and translations up to ±8 mm. Dose distributions were measured for a static case, regular repainting (repeated meander-like path) and random repainting. Latter was realized by randomly distributing single spot locations during irradiation for a given number of repaints. Efficiency of repainting was analyzed by comparison to the static case. A simulation tool based on treatment logs and motion information was developed to compare measurement results to expected dose distributions. Results: Regular repainting could reduce motion artefacts, but dose distortion was strongly dependent on motion direction. Random repainting with same number of repaints (N=4) showed superior results, independent of target movement direction, while introducing slight penalty on delivery times, caused by an increase of overall scanning travel distance. The simulation tool showed good agreement to measured results. Conclusion: The results demonstrate significant improvement in terms of dose conformity when layered repainting is applied in a randomized fashion. This allows for reduced target margins during treatment planning and limited number of repaints. A combination with e.g. respiratory gating is straight-forward. Authors are employees of Varian Medical Systems Particle Therapy GmbH.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

New Virtual Test Bed Capabilities: Virtual DOME Model and New Updates to Repository

The Department of Energy (DOE) Office of Nuclear Energy National Reactor Innovation Center accelerates the deployment of novel reactor concepts by establishing both physical and virtual spaces for building and testing various components, systems, and complete pilot plants. The Virtual Test Bed represents the virtual arm of the National Reactor Innovation Center and is a joint effort with the DOE Nuclear Energy Advanced Modeling and Simulation Program. The Virtual Test Bed mission is to accelerate the deployment of advanced reactors by facilitating the adoption of cutting-edge DOE advanced modeling and simulation tools to design, evaluate, and license reactors. This is primarily achieved by storing example challenge problems in an externally available repository and by developing models to fill the M&S gaps needed for potential demonstrators. Activities conducted this fiscal year focused on developing of a Demonstration of Microreactor Experiments shield model to help accelerate the confirmatory analysis required for the reactor demonstration. This model and workflow will allow developers to leverage advanced modeling and simulation tools to ensure their reactor demonstration concept will meet dose requirements and that the surrounding shield will stay within concrete temperature limits during steady-state and transient operation conditions. An initial model has been developed to evaluate the temperature distribution in the concrete shield during steady-state operation, including neutron and gamma heating effects. Various modeling strategies have been examined to understand their applicability and limitations with different reactor designs to make the workflow as reactor-agnostic as possible and computationally effective to maximize its usability. In addition to describing the Demonstration of Microreactor Experiments shield model and associated results, this report summarizes other accomplishments regarding repository maintenance and improvement and new external models hosted on the repository.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

The Next-Generation Ecosystem Experiment Arctic Rainfall Simulator: a tool to understand the effects of changing rainfall patterns in the Arctic

Rainfall frequency and intensity are expected to increase in the Arctic, with potential detrimental impacts on permafrost, leading to enhanced thawing and carbon release to the atmosphere. However, there have been very few studies on the effect of discrete rain events on permafrost in the Arctic and sub-Arctic. Conducting controlled rainfall experiments within permafrost landscapes can provide an improved understanding of the effect of changing intensity, duration, and timing of rain events on permafrost tundra ecosystems. Here, we describe the design and implementation of the Next-Generation Ecosystem Experiment Arctic Rainfall Simulator (NARS), a variable intensity (4–82 mm/h) rainfall simulator that can be used to study the effects of rainfall on permafrost stability. The NARS design includes a 3D-printed 4 cm H-flume and uses an eTape resistivity sensor that was calibrated (R 2 = 0.9–0.96) to measure discharge from the system. NARS is designed to be lightweight, simple to construct, and can be easily deployed in remote locations. As a field validation of updated rainfall simulator design and modernized controls, NARS was tested on the Seward Peninsula, AK. Because of its portability, versatility in deployment, dimensions, and rainfall intensity, NARS represents a methodological innovation for researching the impacts of rainfall on permafrost environments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

DESIGN AND EXPERIMENTAL STUDY OF A HIGH PRESSURE AND SUPERCRITICAL METHANE-OXYGEN BURNER

Directly heated supercritical oxy-fuel power cycles have potential to offer a higher thermal efficiency and lower pollutant emissions compared to existing power cycles. Recent thermodynamic analysis of the cycle performed by several groups including the UTEP-Air Liquide research team show that combustion in the vicinity of 300 bar pressure and 1000-1400 K temperature allows for relatively high system efficiencies while operating within the limit of practical combustor materials. However, the realization of directly heated supercritical power cycle requires combustion systems be designed to operate in supercritical conditions and at temperature far below the blowout limit of conventional flames (above 1500 K), where not only the thermodynamic properties but also the combustion properties and kinetics are unexplored. To minimize these knowledge gaps, some intermediate pressure ranges (up to 20 bar) are experimented and modeled using a CFD simulation tool. The knowledge obtained from the high pressure test will assist in understanding the combustion chamber pressurization mechanism, ignition and flame behavior at the elevated pressure. This is a systematic first step in testing at higher pressures of 100 and 300 bar pressures. The primary objective of this dissertation is to perform qualitative analysis on oxy-methane combustion at high pressure (< 20 bar) and compare to CFD model for future scale up to supercritical condition. For modeling of the system, a commercial computational fluid dynamics simulation tool, ANSYS Fluent, is used. The inlet conditions for the CFD analysis are obtained from the experiments. The geometry used for the study is same as the test apparatus and operates at the same power ratings. The expertise gained from this experimental study is important to accurately and safely design combustors at elevated pressures up to 300bar. A high pressure methane-oxygen fuel burner and combustor have been developed to accommodate oxycombustion environment. A detailed CFD analysis is conducted to understand the flame length and cooling phenomenon inside the combustor first at pressures up to 20 bar to compare to experiments then at supercritical pressure. The first part of the dissertation compares CFD results with high pressure burner experiments at pressures up to 20 bar. The test is conducted for 30 s within 150 kW – 250 kW thermal power inputs. Afterwards, based on the simulation results a cooling system is proposed for steady state high pressure combustion experiments. It is observed from the experimental study that the proposed pressurization mechanism able to pressurize the vessel up to 20 bar. It is also determined that the current combustor can operate up to 20 bar for short term. The steady state CFD simulation demonstrates that a cooling system must be incorporated for continuous operation. Additionally, the study focuses on the ignition delay due to added diluent during the combustion process. It is observed that as the carbon dioxide recirculation ratio increases the ignition delay time increases. Later part of the dissertation provides a preliminary guideline for developing a laboratory scale supercritical oxy combustor. The employment of supercritical fluid in gas turbine is fairly a new concept. Hence, there are many questions that need to be answered. During this dissertation study an investigation is performed to understand the impact of equation of state in the supercritical combustion simulation. Finally, a CFD model for developing supercritical oxy-methane burner and the combustor, incorporating real gas model, is presented. Although Lee-Kesler equation of state provides better accuracy than Peng Robinson, due to the computational time Peng Robinson equation of state is used for the simulation. It is found that existing knowledge is not enough to simulate combustion at supercritical phase. Detailed combustion chemistry at the supercritical condition must be developed and should be incorporated to accurately replicate the combustion phenomenon.

Chowdhury, ASM Arifur Rahim↗
DOE OSTI2014

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↗
DOE OSTI2018

Phenomena Important in Molten Salt Reactor Simulations

The U.S. Nuclear Regulatory Commission (NRC) is preparing for the future licensing of advanced reactors that will be very different from current light water reactors. Part of the NRC preparation strategy is to identify the simulation tools that will be used for confirmatory safety analysis of normal operation and abnormal situations in those reactors. This report advances that strategy for reactors that will use molten salts (MSRs). This includes reactors with the fuel within the salt as well as reactors using solid fuel. Although both types are discussed in this report, the emphasis is on those reactors with liquid fuel because of the perception that solid-fuel MSRs will be significantly easier to simulate. These liquid-fuel reactors include thermal and fast neutron spectrum alternatives. The specific designs discussed in the report are a subset of many designs being considered in the U.S. and elsewhere but they are considered the most likely to submit information to the NRC in the near future. The objective herein, is to understand the design of proposed molten salt reactors, how they will operate under normal or transient/accident conditions, and what will be the corresponding modeling needs of simulation tools that consider neutronics, heat transfer, fluid dynamics, and material composition changes in the molten salt. These tools will enable the NRC to eventually carry out confirmatory analyses that examine the validity and accuracy of applicant’s calculations and help determine the margin of safety in plant design.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Efficiency-Optimized Dual Fuel Engine with In-Cylinder Gasoline/CNG Blending

The objective of the proposed project is to demonstrate targeted blending of gasoline and natural gas as a means to reduce petroleum consumption of light-duty spark-ignition engine vehicles by at least 50%. This goal will be achieved through petroleum displacement with natural gas (NG) combined with a relative engine efficiency improvement in excess of 10% while meeting retrofitability requirements with an anticipated payback period of less than 36 months. Further through NG direct injection (DI) a 10% improvement in power density over the PFI baseline will be demonstrated. The proposed project employs a sound R&D approach and combines experimental work with state-of-the-art computational modeling to be able to demonstrate the ambitious project goals. The project spans a 30 month period and is divided in 3 phases with the initial phase focused on experimental characterization of the gaseous jet structure of the high-pressure NG DI injector which is used to validate the 3D-CFD simulation tool. The phase also includes an experimental demonstration of the fundamental feasibility of the proposed dual-fuel approach on an automotive single-cylinder research engine. The scope of the second phase includes optimization of the combustion system by means of virtual design combined with 3D-CFD analysis. The resulting optimized combustion system will be built and experimentally evaluated by the end of the 12 month period. The third and final phase expands the scope by integrating the engine results in a vehicle level simulation tool with the objective of developing an optimized engine control strategy that maximizes the efficiency benefits attainable with the limited on-board NG supplies. A detailed cost analysis to verify that the 36 month amortization period requirement is met will conclude the proposed effort.

42 ENGINEERING↗
NASA NTRSDate not supplied

The Influence of Computer Architecture on Performance and Scaling for Hypersonic Flow Simulations

It is critical to understand how hypersonic simulation tools perform on a range of computational platforms. This information will aid in the acquisition of appropriate hardware and the potential refactoring of hypersonic codes to run on different systems. In this paper, we consider two representative high-speed reacting flow cases: a model Mach 8 hypersonic waverider glide vehicle and a model hydrocarbon-fueled hypersonic ramjet propulsion system. In both scenarios, the flow fields are in chemical non-equilibrium and are modeled by the multi-species reacting Navier-Stokes equations. For these simulations we use several hypersonic simulation tools, including US3D, Kestrel, FUN3D, and the JENRER© flow solver. We explore several high performance computing systems containing IntelR© XeonR© Platinum processors, AMD EPYCTM7702 processors, and NVIDIAR© Tesla V100 devices. We compare performance and strong scaling between the different systems.

CPU↗
DOE OSTI2023

Real-Sim Flexible Interface for X-in-the-loop Simulation (FIXS)

Real-Sim is a multi-resolution X-in-the-loop simulation framework to support inclusive testing and evaluation of CAVs before actual implementation in the real world. Real-Sim framework connects different simulation components and integrates various traffic and vehicle simulation tools, virtual environments, and XIL systems through a Flexible Interface for XIL Simulation (FIXS). Connections among different tools are transparent to users and handled by the interface behind the scenes. Co-simulation of various vehicle and traffic simulation tools with different XIL systems can be easily achieved and become a transparent "plug-and-play" process to users.

Shao, Yunli (0000000185287925)↗
NASA NTRS2010

Microstructure Modeling of 3rd Generation Disk Alloys

The objective of this program is to model, validate, and predict the precipitation microstructure evolution, using PrecipiCalc (QuesTek Innovations LLC) software, for 3rd generation Ni-based gas turbine disc superalloys during processing and service, with a set of logical and consistent experiments and characterizations. Furthermore, within this program, the originally research-oriented microstructure simulation tool will be further improved and implemented to be a useful and user-friendly engineering tool. In this report, the key accomplishment achieved during the second year (2008) of the program is summarized. The activities of this year include final selection of multicomponent thermodynamics and mobility databases, precipitate surface energy determination from nucleation experiment, multiscale comparison of predicted versus measured intragrain precipitation microstructure in quench samples showing good agreement, isothermal coarsening experiment and interaction of grain boundary and intergrain precipitates, primary microstructure of subsolvus treatment, and finally the software implementation plan for the third year of the project. In the following year, the calibrated models and simulation tools will be validated against an independently developed experimental data set, with actual disc heat treatment process conditions. Furthermore, software integration and implementation will be developed to provide material engineers valuable information in order to optimize the processing of the 3rd generation gas turbine disc alloys.

Jou, Herng-Jeng↗
DOE OSTI2023

Validating Connected, Automated, and Electric Vehicle Models and Simulation - Research Performance Progress Report

The objective of this project is to test connected and automated vehicles with both electrified and internal combustion engine powertrains to support updates and validation of modeling and simulation tools. This includes the development of the components and network architecture to execute and collect empirical data for multiple scenarios and traffic interactions. Specific program objectives include: • Translate Lab algorithms into vehicle and infrastructure controls • Conduct physical testing at realistic scale • Evaluate system performance • Improve models using empirical data • Improve control algorithms from lessons learned • Identify system and algorithm assumptions which need refinement It is important to note that the objective of this project was not to demonstrate the efficacy of the selected algorithms to improve energy efficiency but rather to validate and improve modeling and simulation tools using empirical data. While it is a desirable outcome to concurrently demonstrate improved energy efficiency through use of these algorithms, and in most cases that was the outcome, the success of this project was not predicated on the performance of the algorithm towards improving energy efficiency across all scenarios and test matrices.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

A Blockchain Based Co-Simulation Framework for Integrating DERs to Wholesale Electricity Markets

As the number of distributed energy resources (DERs) continue to increase across energy-delivery systems, there remains a need for integrating their capabilities into traditional grid operations. In this paper, a blockchain-based solution is proposed to facilitate FERC's Order No. 2222 implementations. The presented use-case enables small-scale DERs to participate in wholesale market operations through DER aggregators, while also enabling local distribution system operators to enforce distribution system constraints in a secure and traceable manner. The presented use case is built around the Unified Testing Platform (UTP) being developed as a part of the Blockchain for Optimized Security and Energy Management (BLOSEM) project. This is a multi-lab effort intended to simplify the deployment of blockchain-powered grid solutions by enabling the integration of simulation tools, and blockchain technologies through the use of system-agnostic interfaces that provide a modular, interoperable, and reusable connectivity layer.

blockchain interoperability↗
DOE OSTI2017

Two-Dimensional Neutronic and Fuel Cycle Analysis of the Transatomic Power Molten Salt Reactor

This status report presents the results from the first phase of the collaboration between Transatomic Power Corporation (TAP) and Oak Ridge National Laboratory (ORNL) to provide neutronic and fuel cycle analysis of the TAP core design through the Department of Energy Gateway for Accelerated Innovation in Nuclear, Nuclear Energy Voucher program. The TAP design is a molten salt reactor using movable moderator rods to shift the neutron spectrum in the core from mostly epithermal at beginning of life to thermal at end of life. Additional developments in the ChemTriton modeling and simulation tool provide the critical moderator-to-fuel ratio searches and time-dependent parameters necessary to simulate the continuously changing physics in this complex system. Results from simulations with these tools show agreement with TAP-calculated performance metrics for core lifetime, discharge burnup, and salt volume fraction, verifying the viability of reducing actinide waste production with this design. Additional analyses of time step sizes, mass feed rates and enrichments, and isotopic removals provide additional information to make informed design decisions. This work further demonstrates capabilities of ORNL modeling and simulation tools for analysis of molten salt reactor designs and strongly positions this effort for the upcoming three-dimensional core analysis.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Performance of Coupled Physics Solvers for Multidisciplinary Hypersonic Flow Simulations on Several Classes of Computer Architectures

The application of hypersonic flow simulation tools to realistic flight scenarios will require the coupling of multiple physical effects to the baseline fluid dynamics. Such multiphysics effects can include the aerooelastic response of the airframe or engine components, dynamic transport of atmospheric particles, the deformation of solid-fluid interfaces that can ablate, pyrolyze, or erode, as well as a host of other processes, all of which are governed by unique sets of physical equations and models. Coupling multiple (and potentially disparate) physics solvers to a robust compressible flow solver poses additional challenges related to the stability, performance and scalability of the combined solver. The choices made during the software design process can therefore lead to a variation in simulation efficiency across different computer architectures. In this paper, we will consider two representative multiphysics hypersonic flow scenarios: the interaction of solid particulates with the flow field created by a hypersonic lifting body and the aerooelastic deformation of a model airframe under high-Mach-number flow conditions. For these simulations we explore the behavior of several hypersonic simulation tools, including Kestrel, FUN3D, US3D, and JENRE multiphysics framework, on several high performance computing systems containing various CPU and GPU architectures.

architecture↗
NASA NTRS2020

Flow Simulations for Martian and Lunar Lander Plume-Surface Interaction Prediction

Landing vehicle propulsive systems introduce several risks associated with plume/surface interactions (PSI) during the final stages of descent. For example, view obscuration during descent due to the liberation of dust driven by the rocket plume’s interaction with the lunar surface was reported on several Apollo missions. Plume-driven debris impacts on vehicle and nearby assets as well as landing instabilities due to crater formation are also potential risks associated with plume/surface interactions. The importance of studying these plume/surface interactions has been elevated due to NASA’s mission to land humans on the Moon by 2024. Existing and under-development predictive simulation tools are currently being brought to application readiness under a project awarded by the NASA Space Technology Mission Directorate (STMD) Game Changing Development (GCD) program. This tool development is divided into four tasks: (1) simulating plume structure in low-pressure environments, (2) simulating crater development and ejecta, (3) regolith particle phase modeling, and (4) gas-particle interaction modeling. Recent efforts under task 1 include demonstrating existing computational fluid dynamics (CFD) tools for simulating plume structure in low pressure environments similar to those found on Mars. The results of these demonstrative simulations are compared to available experimental data. Emerging capability being developed under task 2 is also presented, demonstrating simulations of crater formation in Lunar conditions.

Thomas P Shurtz↗
DOE OSTI2015

Direct observation of electrothermal instability structures in the skin layer of an intensely Ohmically heated conductor

Magnetically driven implosions (MDIs) on the Z Facility assemble high-energy-density plasmas for radiation effects and ICF experiments. MDIs are hampered by the Magneto-Rayleigh-Taylor (MRT) instability, which can grow to large amplitude from a small seed perturbation, limiting achievable stagnation pressures and temperatures. The metallic liners used in Magnetized Liner Inertial Fusion (MagLIF) experiments include astonishingly small (-10 nm RMS) initial surface roughness perturbations; nevertheless, unexpectedly large MRT amplitudes are observed in experiments. An electrothermal instability (ETI) may provide a perturbation which exceeds the initial surface roughness. For a condensed metal resistivity increases with temperature. Locations of higher resistivity undergo increased Ohmic heating, resulting in locally higher temperature, and thus still higher resistivity. Such unstable temperature (and pressure) growth produces density perturbations when the locally overheated metal changes phase, providing the seed perturbation for MRT growth. ETI seeding of MRT on thick conductors carrying current in a skin layer has thus far only been inferred by evaluating MRT amplitude late in the experiment. A direct observation of ETI is vital to ensure our simulation tools are accurately representing the seed of the deleterious MRT instability. In this LDRD project, ETI growth was directly observed on the surface of 1.0-mm-diameter solid Al rods which were pulsed with 1 MA of current in 100 ns. Fine structures resulting from ETI-driven temperature variations were observed directly through high resolution gated optical imaging. Data from two Aluminum alloys (6061 and 5N) and a variety fabrication techniques (conventional machining, single-point diamond turned, electropolished) enable evaluation of which imperfections provide a seed for ETI growth and subsequent plasma initiation. Data is relevant to the early stages of MagLIF liner implosions, when the ETI seed of MRT may be initiated, and provides a fundamentally new dataset with which to test our state-of-the-art simulation tools.

42 ENGINEERING↗
DOE OSTI2022

Verification of MOOSE/Bison's Heat Conduction Solver Using Combined Spatiotemporal Convergence Analysis

Bison is a computational physics code that uses the finite element method to model the thermo-mechanical response of nuclear fuel. Since Bison is used to inform high-consequence decisions, it is important that its computational results are reliable and predictive. One important step in assessing the reliability and predictive capabilities of a simulation tool is the verification process, which quantifies numerical errors in a discrete solution relative to the exact solution of the mathematical model. One step in the verification process—called code verification—ensures that the implemented numerical algorithm is a faithful representation of the underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. In this paper, the code verification process is applied to spatiotemporal heat conduction problems in Bison. Simultaneous refinement of the discretization in space and time is employed to reveal any potential mistakes in the numerical algorithms for the interactions between the spatial and temporal components of the solution. For each verification problem, the correct spatial and temporal order of accuracy is demonstrated for both first- and second-order accurate finite elements and a variety of time-integration schemes. Furthermore, these results provide strong evidence that the Bison numerical algorithm for solving spatiotemporal problems reliably represents the underlying mathematical model in MOOSE. The selected test problems can also be used in other simulation tools that numerically solve for conduction or diffusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗