Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional performance testing”

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,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

STAT7 v1.2 User Guide: The STAT7 Code for Statistical Propagation of Uncertainties in Steady-State Thermal Hydraulics Analysis of Plate-Fueled Reactors

The STAT7 software was developed to perform steady-state, single-phase thermal hydraulics analysis of plate-fueled reactors based on statistical propagation of uncertainties. Application of the software is for non-power research and test reactors, including conversion to low-enriched uranium fuel of U.S. High-Performance Research Reactors such as Massachusetts Institute of Technology Research Reactor. Since it can be necessary to repeat analysis during fuel reloading, STAT7 accommodates flexibility in analyzing many realistic aspects of reactor fuel management. STAT7 uses a Monte Carlo approach to model uncertainty in common fuel fabrication parameters and other key reactor operating parameters required for thermal hydraulics analyses of research and test reactors. These safety calculations are ultimately intended to protect against high fuel plate temperatures due to critical heat flux, or onset of flow instability. STAT7 supports water properties based on the IAPWS-IF97 functions (The International Association for the Properties of Water and Steam Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam) in addition to the fit functions. STAT7 predicts axial profiles of fuel, cladding, and coolant temperature along a lateral stripe that runs the full length of the fuel plate from the bottom to the top. STAT7 can simultaneously analyze all of the axial nodes of all of the fuel plates and all of the coolant channels for one latera stripe of a fuel element. Power splits are calculated for each axial node of each plate to determine how much of the power goes out each face of the plate. By running STAT7 multiple times, full core analysis can be performed by analyzing the margin to onset of nucleate boiling and onset of flow instability for each axial node of each stripe of each plate of each fuel element in the core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers

We report accurate simulations of flows in stellar interiors are crucial to improving our understanding of stellar structure and evolution. Because the typically slow flows are merely tiny perturbations on top of a close balance between gravity and the pressure gradient, such simulations place heavy demands on numerical hydrodynamics schemes. We demonstrate how discretization errors on grids of reasonable size can lead to spurious flows orders of magnitude faster than the physical flow. Well-balanced numerical schemes can deal with this problem. Three such schemes were applied in the implicit, finite-volume SEVEN-LEAGUE HYDRO code in combination with a low-Mach-number numerical flux function. We compare how the schemes perform in four numerical experiments addressing some of the challenges imposed by typical problems in stellar hydrodynamics. We find that the α-β and deviation well-balancing methods can accurately maintain hydrostatic solutions provided that gravitational potential energy is included in the total energy balance. They accurately conserve minuscule entropy fluctuations advected in an isentropic stratification, which enables the methods to reproduce the expected scaling of convective flow speed with the heating rate. The deviation method also substantially increases accuracy of maintaining stationary orbital motions in a Keplerian disk on long timescales. The Cargo–LeRoux method fares substantially worse in our tests, although its simplicity may still offer some merits in certain situations. Overall, we find the well-balanced treatment of gravity in combination with low Mach number flux functions essential to reproducing correct physical solutions to challenging stellar slow-flow problems on affordable collocated grids.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Status on the Development of a Modeling and Simulation Framework for the Economic Assessment of Nuclear Hybrid Energy Systems

An effort to design and build a modeling and simulation framework to assess the economic viability of Nuclear Hybrid Energy Systems (NHES) was undertaken in fiscal year 2015 (FY15). The purpose of this report is to document the various tasks associated with the development of such a framework and to provide a status on its progress. Several tasks have been accomplished. First, starting from a simulation strategy, a rigorous mathematical formulation has been achieved in which the economic optimization of a Nuclear Hybrid Energy System is presented as a constrained robust (under uncertainty) optimization problem. Some possible algorithms for the solution of the optimization problem are presented. A variation of the Simultaneous Perturbation Stochastic Approximation algorithm has been implemented in RAVEN and preliminary tests have been performed. The development of the software infrastructure to support the simulation of the whole NHES has also moved forward. The coupling between RAVEN and an implementation of the Modelica language (OpenModelica) has been implemented, migrated under several operating systems and tested using an adapted model of a desalination plant. In particular, this exercise was focused on testing the coupling of the different code systems; testing parallel, computationally expensive simulations on the INL cluster; and providing a proof of concept for the possibility of using surrogate models to represent the different NHES subsystems. Another important step was the porting of the RAVEN code under the Windows™ operating system. This accomplishment makes RAVEN compatible with the development environment that is being used for dynamic simulation of NHES components. A very simplified model of a NHES on the electric market has been built in RAVEN to confirm expectations on the analysis capability of RAVEN to provide insight into system economics and to test the capability of RAVEN to identify limit surfaces even for stochastic constraints. This capability will be needed in the future to enforce the stochastic constraints on the electric demand coverage from the NHES. The development team gained experience with many of the tools that are currently envisioned for use in the economic analysis of NHES and completed several important steps. Given the complexity of the project, preference has been given to a structural approach in which several independent efforts have been used to build the cornerstone of the simulation framework. While this is good approach in establishing such a complex framework, it may delay reaching more complete results on the performance of analyzed system configurations. The integration of the previously reported exergy analysis approach was initially proposed as part of this milestone. However, in reality, the exergy-based apportioning of cost will take place only in a second stage of the implementation since it will be used to properly allocate cost among the different NHES subsystems. Therefore, exergy does not appear at the level of the main drivers in the analysis framework; the latter development of the base framework is the focus of this report.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

PDMS ‐silica composite gas separation membranes by direct ink writing

Abstract Polydimethylsiloxane (PDMS)‐based membranes containing amine‐functionalized and unfunctionalized silica particles were fabricated via direct ink writing for CO 2 /N 2 gas separation. The printability of the inks was evaluated by rheological measurements, while spectroscopy, microscopy, thermal measurements, and mechanical testing were employed to characterize the printed membranes. The surface morphology of the membranes revealed the absence of voids, demonstrating their suitability for gas separation. The printed membranes also exhibited desirable thermal and mechanical properties (i.e., thermal degradation temperature of 518 °C and tensile strength as high as 1.178 MPa with 529% elongation). The PDMS‐based membranes generally displayed high permeability for CO 2 but slightly low selectivity for the CO 2 /N 2 gas pair. The best‐combined permeability‐selectivity performance of 8794 barrer and selectivity of 11.64 was demonstrated by the printed PDMS membrane containing no SiO 2 fillers. The inclusion of unfunctionalized SiO 2 particles generally increased the membrane's gas permeability but compromised the selectivity. In contrast, membranes with amine‐functionalized silica showed improved selectivity compared to membranes containing unfunctionalized silica. Overall, the performance and characteristics offered by the PDMS/silica composite membranes demonstrated the potential of 3D printing as an economical and sustainable fabrication approach to developing materials for carbon capture applications.

Gutierrez, Dianne B.↗
DOE OSTI2016

On the applicability of one- and many-electron quantum chemistry models for hydrated electron clusters

We evaluate the applicability of a hierarchy of quantum models in characterizing the binding energy of excess electrons to water clusters. In particular, we calculate the vertical detachment energy of an excess electron from water cluster anions with methods that include one-electron pseudopotential calculations, density functional theory (DFT) based calculations, and ab initio quantum chemistry using MP2 and eom-EA-CCSD levels of theory. The examined clusters range from the smallest cluster size (n = 2) up to nearly nanosize clusters with n = 1000 molecules. The examined cluster configurations are extracted from mixed quantum-classical molecular dynamics trajectories of cluster anions with n = 1000 water molecules using two different one-electron pseudopotenial models. We find that while MP2 calculations with large diffuse basis set provide a reasonable description for the hydrated electron system, DFT methods should be used with precaution and only after careful benchmarking. Strictly tested one-electron psudopotentials can still be considered as reasonable alternatives to DFT methods, especially in large systems. The results of quantum chemistry calculations performed on configurations, that represent possible excess electron binding motifs in the clusters, appear to be consistent with the results using a cavity structure preferring one-electron pseudopotential for the hydrated electron, while they are in sharp disagreement with the structural predictions of a non-cavity model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Structure/property relationships in copolymers comprising renewable isosorbide, glucarodilactone, and 2,5-bis(hydroxymethyl)furan subunits

Carbohydrates and their derivatives have great potential as building blocks for the development of renewable materials that are cost and performance competitive with conventional petroleum-based materials. To this end, the rigid carbohydrate-derived diols, glucarodilactone and isosorbide were functionalized with a castor oil derivative, 10-undecenoic acid, forming α,ω-dienes suitable for acyclic diene metathesis (ADMET) polymerization. Equimolar copolymerizations of these two monomers, glucarodilactone undecenoate (GDLU) and isosorbide undecenoate (GDLU), transform brittle homopolymers into elastic materials with shape memory capabilities. To understand the source of these properties, a series of copolymers consisting of a range of GDLU and IU compositions (100:0, 77:23, 76:24, 52:48, 48:52, 18:82, 0:100 mol percent GDLU to IU) was synthesized. The material and chemical properties were characterized by uniaxial tensile testing, small-amplitude oscillatory shear rheology, X-ray scattering and hydrolytic degradation testing. Small compositional changes in this family were shown to have a drastic impact on the observed mechanical performance and degradability of these materials. Rheological measurements of GDLU-containing copolymers found evidence for the presence of transient networking within these materials. We posit that the transient network within this material is responsible for the elasticity and shape memory abilities of the GDL-containing materials. In confirmation of these properties, a novel α,ω-diene 2,5-bis(hydroxymethylfuran) undecenoate (BHMFU) was developed and served as a direct replacement for either GDLU or IU in copolymerizations. The replacement of GDLU resulted in a total loss of the elasticity and shape memory of the copolymers, supporting the role of GDLU in the mechanical performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

First simultaneous measurement of the γ-ray and neutron emission probabilities in inverse kinematics at a heavy-ion storage ring

The probabilities for γ-ray and particle emission as a function of the excitation energy of a decaying nucleus are valuable observables for constraining the ingredients of the models that describe the deexcitation of nuclei near the particle emission threshold. These models are essential in nuclear astrophysics and applications. In this paper, we have for the first time simultaneously measured the γ-ray and neutron emission probabilities of 208 Pb. The measurement was performed in inverse kinematics at the Experimental Storage Ring (ESR) of the GSI/FAIR facility, where a 208 Pb beam interacted through the 208 Pb (p,p') reaction with a hydrogen gas jet target. Instead of detecting the γ rays and neutrons emitted by 208 Pb, we detected the heavy beamlike residues produced after γ and neutron emission. These heavy residues were fully separated by a dipole magnet of the ESR and were detected with outstanding efficiencies. The comparison of the measured probabilities with model calculations has allowed us to test and select different descriptions of the γ-ray strength function and the nuclear level density available in the literature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Biaxial Strains Mediated Oxygen Reduction Electrocatalysis on Fenton Reaction Resistant L1 0 ‐PtZn Fuel Cell Cathode

Abstract PtM alloy catalysts (e.g., PtFe, PtCo), especially in an intermetallic L1 0 structure, have attracted considerable interest due to their respectable activity and stability for the oxygen reduction reaction (ORR) in proton exchange membrane fuel cells (PEMFCs). However, metal‐catalyzed formation of ·OH from H 2 O 2 (i.e., Fenton reaction) by Fe‐ or Co‐containing catalysts causes severe degradation of PEM/catalyst layers, hindering the prospects of commercial applications. Zinc is known as an antioxidant in Fenton reaction, but is rarely alloyed with Pt owing to its relatively negative redox potential. Here, sub‐4 nm intermetallic L1 0 ‐PtZn nanoparticles (NPs) are synthesized as high‐performance PEMFC cathode catalysts. In PEMFC tests, the L1 0 ‐PtZn cathode achieves outstanding activity (0.52 A mg Pt −1 at 0.9 V iR ‐free , and peak power density of 2.00 W cm −2 ) and stability (only 16.6% loss in mass activity after 30 000 voltage cycles), exceeding the U.S. DOE 2020 targets and most of the reported ORR catalysts. Density function theory calculations reveal that biaxial strains developed upon the disorder‐order (A1L1 0 ) transition of PtZn NPs would modulate the surface PtPt distances and optimize PtO binding for ORR activity enhancement, while the increased vacancy formation energy of Zn atoms in an ordered structure accounts for the improved stability.

Liang, Jiashun↗
DOE OSTI2025

Development and evaluation of a multi-functional heat pump with embedded thermal storage

We developed and tested a novel multi-functional packaged vertical heat pump designed for multi-family buildings, capable of providing space cooling, space heating, water heating, and energy storage integration. The system employs a 3-speed scroll compressor and supports energy-efficient cooling and heating across different ambient conditions, while utilizing both indoor and outdoor air sources for water heating. Key performance metrics include an integrated energy efficiency ratio of 18.0 for cooling, a heating seasonal performance factor of 10.0, and a coefficient of performance (COP) of 2.1 for heating in cold climates down to −15 ˚C. Additionally, the system delivers an annual water heating COP greater than 4.0, with outstanding performance in combined space cooling and water heating modes, achieving a total COP of 8.8 by recovering condenser waste heat. The heat pump also demonstrated thermal energy storage capabilities, integrating with phase change material systems to store heating and cooling energy. Furthermore, this versatile system offers a promising solution for improving energy efficiency and sustainability in multi-family buildings, providing high-performance heating, cooling, and water heating across a range of climates.

COP↗
DOE OSTI2017

The structure and propagation of laminar flames under autoignitive conditions

The laminar flame speed s l is an important reference quantity for characterising and modelling combustion. Experimental measurements of laminar flame speed require the residence time of the fuel/air mixture (τ f ) to be shorter than the autoignition delay time (τ). This presents a considerable challenge for conditions where autoignition occurs rapidly, such as in compression ignition engines. As a result, experimental measurements in typical compression ignition engine conditions do not exist. Simulations of freely propagating premixed flames, where the burning velocity is found as an eigenvalue of the solution, are also not well posed in such conditions, since the mixture ahead of the flame can autoignite, leading to the so called “cold boundary problem”. In this paper, a numerical method for estimating a reference flame speed, s R , is proposed that is valid for laminar flame propagation at autoignitive conditions. Two isomer fuels are considered to test this method: ethanol, which in the considered conditions is a single-stage ignition fuel; and dimethyl ether, which has a temperature-dependent single- or two-stage ignition and a negative temperature coefficient regime for τ. Calculations are performed for the flame position in a one-dimensional computational domain with inflow-outflow boundary conditions, as a function of the inlet velocity U I and for stoichiometric fuel–air premixtures. The response of the flame position, L F , to U I shows distinct stabilisation regimes. For single-stage ignition fuels, at low U I the flame speed exceeds U I and the flame becomes attached to the inlet. Above a critical U I value, the flame detaches from the inlet and L f becomes extremely sensitive to U I until, for sufficiently high U I , the sensitivity decreases and L f corresponds to the location expected from a purely autoignition stabilised flame. The transition from the attached to the autoignition regimes has a corresponding peak dL f /dU I value which is proposed to be a unique reference flame speed s R for single-stage ignition fuels. For two-stage ignition fuels, there is an additional stable regime where a high-temperature flame propagates into a pool of combustion intermediates generated by the first stage of autoignition. This results in two peaks in dL f /dU I and therefore two reference flame speed values. The lower value corresponds to the definition of s R for single-stage ignition fuels, while the higher value exists only for two-stage ignition fuels and corresponds to a high temperature flame propagating into the first stage of autoignition and is denoted s R ' . Finally, a transport budget analysis for low- and high-temperature radical species is also performed, which confirms that the flame structures at U I = s R and U I = s R ' do indeed correspond to premixed flames (deflagrations), as opposed to spontaneous ignition fronts which do not have a unique propagation speed.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2014

A high-order discontinuous Galerkin method for fluid–structure interaction with efficient implicit–explicit time stepping

We present a high-order accurate scheme for coupled fluid–structure interaction problems. The fluid is discretized using a discontinuous Galerkin method on unstructured tetrahedral meshes, and the structure uses a high-order volumetric continuous Galerkin finite element method. Standard radial basis functions are used for the mesh deformation. The time integration is performed using a partitioned approach based on implicit–explicit Runge–Kutta methods. The resulting scheme fully decouples the implicit solution procedures for the fluid and the solid parts, which we perform using two separate efficient parallel solvers. We demonstrate up to fifth order accuracy in time on a non-trivial test problem, on which we also show that additional subiterations are not required. We solve a benchmark problem of a cantilever beam in a shedding flow, and show good agreement with other results in the literature. Finally, we solve for the flow around a thin membrane at a high angle of attack in both 2D and 3D, and compare with the results obtained with a rigid plate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Bayesian analysis for estimating statistical parameter distributions of elasto-viscoplastic material models

High temperature design methods rely on constitutive models for inelastic deformation and failure typically calibrated against the mean of experimental data without considering the associated scatter. Variability may arise from the experimental data acquisition process, from heat-to-heat material property variations, or both and need to be accurately captured to predict parameter bounds leading to efficient component design. Applying the Bayesian Markov Chain Monte Carlo (MCMC) method to produce statistical models capturing the underlying uncertainty in the experimental data is an area of ongoing research interest. This work varies aspects of the Bayesian MCMC method and explores their effect on the posterior parameter distributions for a uniaxial elasto-viscoplastic damage model using synthetically generated reference data. From our analysis with the uniaxial inelastic model we determine that an informed prior distribution including different types of test conditions results in more accurate posterior parameter distributions. The parameter posterior distributions, however, do not improve when increasing the number of similar experimental data. Additionally, changing the amount of scatter in the data affects the quality of the posterior distributions, especially for the less sensitive model parameters. Moreover, we perform a sensitivity study of the model parameters against the likelihood function prior to the Bayesian analysis. The results of the sensitivity analysis help to determine the reliability of the posterior distributions and reduce the dimensionality of the problem by fixing the insensitive parameters. The comprehensive study described in this work demonstrates how to efficiently apply the Bayesian MCMC methodology to capture parameter uncertainties in high temperature inelastic material models. Quantifying these uncertainties in inelastic models will improve high temperature engineering design practices and lead to safer, more effective component designs.

42 ENGINEERING↗
DOE OSTI2017

Generalized fluid theory including non-Maxwellian kinetic effects

The results obtained by the plasma physics community for the validation and the prediction of turbulence and transport in magnetized plasmas come mainly from the use of very central processing unit (CPU)-consuming particle-in-cell or (gyro)kinetic codes which naturally include non-Maxwellian kinetic effects. To date, fluid codes are not considered to be relevant for the description of these kinetic effects. Here, after revisiting the limitations of the current fluid theory developed in the 19th century, we generalize the fluid theory including kinetic effects such as non-Maxwellian super-thermal tails with as few fluid equations as possible. The collisionless and collisional fluid closures from the nonlinear Landau Fokker–Planck collision operator are shown for an arbitrary collisionality. Indeed, the first fluid models associated with two examples of collisionless fluid closures are obtained by assuming an analytic non-Maxwellian distribution function. One of the main differences with the literature is our analytic representation of the distribution function in the velocity phase space with as few hidden variables as possible thanks to the use of non-orthogonal basis sets. These new non-Maxwellian fluid equations could initiate the next generation of fluid codes including kinetic effects and can be expanded to other scientific disciplines such as astrophysics, condensed matter or hydrodynamics. As a validation test, we perform a numerical simulation based on a minimal reduced INMDF fluid model. The result of this test is the discovery of the origin of particle and heat diffusion. The diffusion is due to the competition between a growing INMDF on short time scales due to spatial gradients and the thermalization on longer time scales. Here, the results shown here could provide the insights to break some of the unsolved puzzles of turbulence.

42 ENGINEERING↗
DOE OSTI2025

Criticality analysis of nuclear binding energy neural networks

Machine learning methods, in particular deep learning methods such as artificial neural networks (ANNs) with many layers, have become widespread and useful tools in nuclear physics. However, these ANNs are typically treated as ‘black boxes’, with their architecture (width, depth, and weight/bias initialization) and the training algorithm and parameters chosen empirically by optimizing learning based on limited exploration. We test a non-empirical approach to understanding and optimizing nuclear physics ANNs by adapting a criticality analysis based on renormalization group flows in terms of the hyperparameters for weight/bias initialization, training rates, and the ratio of depth to width. This treatment utilizes the statistical properties of neural network initialization to find a generating functional for network outputs at any layer, allowing for a path integral formulation of the ANN outputs as a Euclidean statistical field theory. We use a prototypical example to test the applicability of this approach: a simple ANN for nuclear binding energies. We find that with training using a stochastic gradient descent optimizer, the predicted criticality behavior is realized, and optimal performance is found with critical tuning. However, the use of an adaptive learning algorithm leads to somewhat superior results without concern for tuning and thus obscures the analysis. Nevertheless, the criticality analysis offers a way to look within the black box of ANNs, which is a first step towards potential improvements in network performance beyond using adaptive optimizers.

artificial neural network↗
DOE OSTI2022

Emulation Framework for Distributed Large-Scale Systems Integration

Recent trends in systems engineering include integration of very large-scale systems, which entails significant challenges when they are geographically dispersed. In these scenarios, intelligent integration of distributed large-scale systems requires significant coordination among hardware elements as well as all software components. The approach of integrated systems (both computing platform and experimental equipment) for end-to-end orchestration is called federation. Virtual frameworks can aid in the testing, assessment, and implementation of a functional system of interconnected resources. We present an emulation framework that replicates the software environments of multi-site federations of computing systems and instruments. Our emulation framework allows systems engineers to reduce developmentcost and avoid disruptions to production infrastructure. Our framework was effectively used to develop and test software modules for various tasks including container orchestration and instrument access. For performance assessment, however, the emulated framework is severely limited in providing accurate network and IO measurements at 10 Gbps and higher data rates. The data transfer performance profiles estimated using these emulated measurements are usually inaccurate for high bandwidth and high latency connections, since emulation does not accurately reflect the critical network transport dynamics.We utilize measurements from a physical testbed with hardware network emulators to obtain data transfer profiles that closely match the expected profiles for the emulated federations. We show the effectiveness of our approach by an illustrative example of integrated (federated) multi-site ultra large-scale systems that are connected via high speed wide area networks.

Imam, Neena↗
DOE OSTI2016

SU-F-T-169: A Periodic Quality Assurance Program for a Spot-Scanning Proton Treatment Facility

Purpose: To develop daily and monthly quality assurance (QA) programs in support of a new spot-scanning proton treatment facility using a combination of commercial and custom equipment and software. Emphasis was placed on efficiency and evaluation of key quality parameters. Methods: The daily QA program was developed to test output, spot size and position, proton beam energy, and image guidance using the Sun Nuclear Corporation rf-DQA™3 device and Atlas QA software. The program utilizes standard Atlas linear accelerator tests repurposed for proton measurements and a custom jig for indexing the device to the treatment couch. The monthly QA program was designed to test mechanical performance, image quality, radiation quality, isocenter coincidence, and safety features. Many of these tests are similar to linear accelerator QA counterparts, but many require customized test design and equipment. Coincidence of imaging, laser marker, mechanical, and radiation isocenters, for instance, is verified using a custom film-based device devised and manufactured at our facility. Proton spot size and position as a function of energy are verified using a custom spot pattern incident on film and analysis software developed in-house. More details concerning the equipment and software developed for monthly QA are included in the supporting document. Thresholds for daily and monthly tests were established via perturbation analysis, early experience, and/or proton system specifications and associated acceptance test results. Results: The periodic QA program described here has been in effect for approximately 9 months and has proven efficient and sensitive to sub-clinical variations in treatment delivery characteristics. Conclusion: Tools and professional guidelines for periodic proton system QA are not as well developed as their photon and electron counterparts. The program described here efficiently evaluates key quality parameters and, while specific to the needs of our facility, could be readily adapted to other proton centers.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Error Analysis of the Thornado Supernova Simulation Code

The iron core collapse of a supernova throws 1.5 solar masses at speeds up to 35% the speed of light into a space with radius 50km, which creates a high density environment capable of generating neutrinos, producing gravitational waves, and driving chemical evolution of matter past iron. By designing a simulation of the process, we can better understand the origins and properties of these rare phenomena which are not seen anywhere else in the universe. Typically, supernovas are simulated with a finite volume method, but we want to incorporate a Nodal Discontinuous Galerkin Method to better utilize high performance computing because the DG method allows us to increase precision without increasing the number of computational elements. This project focuses on the fluid dynamics during collapse and the resulting shockwaves. To understand this, we improved on an existing simulation by comparing benchmark tests against published and theoretical results to identify and correct sources of error. The DG Method creates discontinuous functions across computational elements by solving the Euler fluid equations computed over Runge-Kutta time steps. The DG method allows for discontinuous functions to be built from polynomials defined in individual elements; however, it can lead to unphysical oscillation near discontinuities. We apply the characteristic slope limiter to reduce these affects, but when applied globally it leads to a 1% increase in absolute error. By introducing a Trouble Cell Indicator, which marks certain fluid elements based on the rate of change of observables, we can constrain the limiter to elements near shocks and manage the resulting absolute error. By constraining the limiter to discontinuous regions, we can remove what we know will be unrealistic without impacting naturally occurring oscillations.

79 ASTRONOMY AND ASTROPHYSICS↗