Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model code”

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 739 records · Page 41

Data Evaluation of Actinide Cross Sections: 230 Th and 231 Th

A new evaluation of the ENDL cross section set for Thorium (Z=90) is developed using the TALYS statistical model cross section code. The primary goal of this effort is to produce an evaluation that attempts to match as closely as possible fission cross sections developed through surrogate reaction techniques on actinide targets 230 Th and 231 Th. This evaluation effort and the processing needed to render its results into data libraries is a necessary step in making the efforts of nuclear experimentalists useful to the broad community of researchers engaged in simulations of nuclear fusion for basic and applied science. Another aspect, verification and validation against various AGEX experiments, is also presented. The end-product is an updated library that includes the latest measurements of fission cross sections for comparison against those measured via traditional techniques. All the steps in the evaluation, processing, validation and verification, and library release are described in the following sections. For completeness, the appendix contains all the parameters used in the TALYS cross section evaluation for neutrons incident on 230 Th and 231 Th.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Materials response under hypersonic flow conditions probed by multi-modal diagnostics in a benchtop wind tunnel

In this LDRD “Materials response under hypersonic flow conditions probed by multi-modal diagnostics in a benchtop wind tunnel” a peer-reviewed lab-scale hypersonic wind tunnel (WT) was planned to be assembled in-house, along with time and space resolved and related flow and materials diagnostics to derive high quality measurements and the science that is able to inform the development of our modeling and simulation codes for the hypersonic regime. The baseline capability to achieve Mach 3 flow was scheduled to be completed in Q4 2020, in the first year of the LDRD. However due to the COVID-19 delays and travel restrictions starting in March, 2020 this task could not be completed to carry out the materials response studies in-house, or those planned with external partners that were engaged during this LDRD (NASA Langley, NASA Ames, Texas A&M university, and Virginia Tech). Still, using our hypersonic computational fluid simulation codes, many of the hypersonic components have been designed (for Mach 3, 5, and 7) and procured (or in procurement phase), and a new dedicated hypersonic laboratory was facilitized, 2 WT test sections have been manufactured, and a Mach 3 nozzle was 3D printed for prototyping. In parallel, unique materials diagnostics have been tested under static conditions and fully specified and awaiting validation in our partner facilities. Due to the new hypersonic initiative announced in January, 2020 we proposed a related scope increase in hypersonics research requiring energy interactions and new diagnostics well beyond what could be supported by an LDRD ER. A decision was made to end the current LDRD ER and formally propose a new LDRD SI that leverages the current effort, which will be starting in FY21. As a result of the uniqueness of this LDRD work, our early efforts have been encouraged and well received in the hypersonic research community, and is now part of multiple unsolicited internal and external project proposals expected to be funded in the coming years, which is encouraging to make LLNL a leader in experimentally probing extreme physics of local energetics critical to hypersonics.

36 MATERIALS SCIENCE↗

Nemchinov-Dyson Solutions of the Two-Dimensional Axisymmetric Inviscid Compressible Flow Equations [Review]

We investigate the two-dimensional (2D) inviscid compressible flow equations in axisymmetric coordinates, constrained by an ideal gas equation of state (EOS). Beginning with the assumption that the 2D velocity field is space-time separable and linearly variable in each corresponding spatial coordinate, we proceed to derive an in nite family of elliptic or hyperbolic, uniformly expanding or contracting \gas cloud" solutions. Construction of specific example solutions belonging to this family is dependent on the solution of a system of nonlinear, coupled, second-order ordinary differential equations, and the prescription of an additional physical process of interest (e.g., uniform temperature or uniform entropy flow). The physical and computational implications of these solutions as pertaining to quantitative code verification or model qualification studies are discussed in some detail.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Confronting Dark Matter with the Multiwavelength Sky (Final Technical Report)

The existence of dark matter is a key clue to new fundamental physics, yet its nature and non-gravitational interactions remain mysterious. In the course of this five-year project (spanning July 2015 - July 2020), we have developed novel analysis techniques to separate dark matter signals from astrophysical backgrounds in the inner Galaxy, and measured hitherto-unknown properties of the Fermi Bubbles and the Galactic Center excess in gamma-rays from the inner Milky Way. We have improved predictions for the spatial distribution and energy spectra of photon signals from dark matter, and set new limits from astrophysical and cosmological datasets on dark matter inhabiting complex dark sectors, multi-particle dark matter annihilation, and tiny primordial black holes. We have successfully adapted effective field theory techniques from collider physics to make high-precision predictions for signals from heavy weakly-interacting DM, allowing us to essentially exclude the possibility that DM is a thermal wino, and have discovered novel theoretical and observational implications for dark matter bound states and self-interactions. We have developed a comprehensive public code package for modeling the effects of exotic energy injections, such as from dark matter annihilation and decay, on the ionization and temperature history of the early cosmos, and have used these tools to set robust, powerful and broadly-applicable constraints on models of light dark matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Saline: An API for Thermophysical Properties

An initial interface for thermophysical data was generalized for broader incorporation into modeling and simulation codes used for molten salt reactor analysis. A quality assurance review of the thermophysical properties data was conducted; errors were fixed, and improvements were made with consistency and accuracy to references. The Saline code application programming interface was created to facilitate access and integration of thermophysical properties. This report describes the thermophysical data, the quality assurance improvements, and the Saline interface.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Summary of Thermodynamic Data Integration (M4SF-23LL010301052)

This progress report (Level 4 Milestone Number M4SF-23LL010301052) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Argillite Activity Number M4SF-23LL010301052. SUPCRTNE is being developed as the primary engine for thermodynamic database development in support of geologic disposal of high-level nuclear waste. As used here, “SUPCRT” refers to both a computer program and its supporting database. Additional letters or numbers refer to specific variants (SUPCRT92, Johnson et al., 1992; SUPCRTBL, Zimmer et al., 2016; and SUPCRTNE). A SUPCRT database contains the data required to calculate the thermodynamic properties of solid, gas, and aqueous species over a wide range of temperature and pressure. Normally SUPCRT is used to create a higher-level data base that directly supports modeling and simulation codes such as EQ3/6, GWB, PHREEQC, and PFLOTRAN.

58 GEOSCIENCES↗

MCNP Simulation of the SFRC and Minimum Detection Threshold Determination (Rev. 1)

In this report the Minimum Detection Threshold (MDT) for the Spent Fuel Rod Counter (SFRC) was determined using Monte Carlo N-Particle transport computer code (MCNP). A model of the SFRC was simulated in addition to a fuel transport pipe and spent Magnox fuel at descending positions down the pipe. The neutron MDT for Magnox fuel free-falling along a discharge pipe was calculated to be equal to 4080 n/s. The gamma MDT for the same free falling fuel as measured by the SFRC’s ionization chamber was determined to be 455 γ/s. These values were determined to be a conservative value to use when the instrument is deployed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

LYNM-PE1 Seismic Parameters from Borehole Log, Laboratory, and Tabletop Measurements

The goal of this work is to provide a database of quality-checked seismic parameters that can be integrated with the Geologic Framework Model (GFM) for the LYNM-PE1 (Low Yield Nuclear Monitoring – Physical Experiment 1) testbed. We integrated data from geophysical borehole logs, tabletop measurements on collected core, and laboratory measurements. We reviewed for internal consistency among each measurement type, documented the caveats of measurement conditions, and integrated lithologic logs to check the validity of outlier values. The resulting consolidated parameter tables can be used as inputs for modeling and analysis codes and are designed to interface with the GFM, which is being actively developed.

58 GEOSCIENCES↗

Faster Tensor Network Decoding for Topological Quantum Codes

We present a fast and Bayes-optimal-approximating tensor network decoder for planar quantum LDPC codes based on the tensor renormalization group algorithm, originally proposed by Levin, and Nave. By precomputing the renormalization group flow for the null syndrome, we need only recompute tensor contractions in the causal cone of the measured syndrome at the time of decoding. This allows us to achieve an overall runtime complexity of ($pnχ^6$) where p is the depolarizing noise rate, and χ is the cutoff value used to control singular value decomposition approximations used in the algorithm. We apply our decoder to the surface code in the code capacity noise model and compare its performance to the original matrix product state (MPS) tensor network decoder introduced by Bravyi, Suchara, and Vargo. The MPS decoder has a p-independent runtime complexity of $\mathcal{O}(nχ^3)$ resulting in significantly slower decoding times compared to our algorithm in the low-p regime.

97 MATHEMATICS AND COMPUTING↗

M3SF-24LL010301052-Summary of SUPCRTNE testing and Rev0 database

This progress report (Level 3 Milestone Number M3SF-24LL010301052) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Argillite Disposal R&D work package SF-24LL01030105. SUPCRTNE is being developed as the primary engine for thermodynamic database development in support of geologic disposal of high-level nuclear waste. As used here, “SUPCRT” refers to both a computer program and its supporting database. Additional letters or numbers refer to specific variants (SUPCRT92, (Johnson et al., 1992); SUPCRTBL, (Zimmer et al., 2016); and SUPCRTNE). A SUPCRT database contains the data required to calculate the thermodynamic properties of solid, gas, and aqueous species over a wide range of temperature and pressure. Normally, SUPCRT is used to create a higher-level database that directly supports modeling and simulation codes such as EQ3/6, GWB, PHREEQC, and PFLOTRAN.

58 GEOSCIENCES↗

ATOMIC Simulations and Experimental Data for Sodium and Copper Mixtures

This data consists of simulations and experimental measurements of laser-induced breakdown spectroscopy (LIBS). The simulations are produced by ATOMIC, a general purpose plasma modeling and kinetics code that has been designed to compute emission (or absorption) spectra from plasmas [1]. Our overall suite of simulations contains 12 sets of simulations: training and validation sets of simulations for each of two resolutions of each of sodium, copper, and mixtures of the two. The training data were produced using a 500-run design that varies input parameters temperature (T), electron density (Ne), proportion of sodium (pNa), and proportion of copper (pCu). The latter two variables sum to one and are unused in the single-element simulations. The validation data was produced with a 25-run design. The coarse simulations produce spectra over a range of 240nm - 880nm that roughly mimics the range collected by the ChemCam instrument on the Mars rover Curiosity. The fine resolution simulations produce spectra over a range of 550nm to 600nm and were produced to mimic the included experimental results. The experimental results cover four mixtures of sodium and copper. All data sets are kept in zip files whose names indicate the elemental composition (NaCu, Na, or Cu), resolution (coarse or fine), and purpose (training or validation) with file names numbered to indicate the line in the design files used to produce the simulation. The designs are provided as text files with names indicating their purpose. The experimental data is provided as a CSV file. See attached abstract for figures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

ATOMIC Simulations and Experimental Data for Basalt-like Compounds

This data consists of simulations and experimental measurements of laser-induced breakdown spectroscopy (LIBS). The simulations are produced by ATOMIC, a general purpose plasma modeling and kinetics code that has been designed to compute emission (or absorption) spectra from plasmas [2]. The makeup of the plasma was considered to be divided into some proportion water, some proportion Martian atmosphere (CO2), and some proportion target (from the rock or object impacted by the laser), where these proportions add to 1. Based on expert knowledge, the proportion of water was kept in the range [0.0,0.5] and the proportion of atmosphere was kept in the range [0.02, 0.9]. Our overall suite of simulations contains six sets of simulations that differ in which elements were considered to make up the target. Within each set, we used uniformly drawn temperatures and log mass densities within pre-specified ranges. The temperature range was [0.5,1.5] eV and the log (base 10) mass density range was [-7,-4]. The proportion of water, atmosphere, and target were drawn from a symmetric Dirichlet distribution, but draws in which the propor- tion of water or atmosphere exceeded the pre-specified limits were rejected from the design. Up to eleven constituent elements (Si, Al, Fe, Mg, Ca, O, Ti, Mn, Na, K, P) were considered for the target, as they are the most common elements found in basalt compounds and were used in [1]. For each run, the proportions of the constituent elements making up the target were drawn from a symmetric Dirichlet distribution. We ran 1,350 simulations that included nonzero proportions of all eleven elements. We also ran simulations which excluded some of these elements. In particular, we ran 1,000 simulations that only included nonzero proportions for the six most common elements (Si, Al, Fe, Mg, Ca, O). We also ran five sets, each with 500 simulations, that only included nonzero proportions for five of the six most common elements (but where all sets included O). Thus, we generated a total of 4,850 spectra representing basalt-like compounds in which the target was comprised of oxygen and between four and ten other elements. The ATOMIC code produced spectra over a range of 240nm - 880nm that roughly mimics the range collected by the ChemCam instrument on the Mars rover Curiosity. Each spectra had 32,000 wavelengths split across three spectrometer ranges (to mimic ChemCam). The experimental data, described in [1], measures a prepared basalt sample. All files are kept in directories whose names indicate the set of elements considered for the target with file names numbered to indicate the line in the design files used to produce the simulation. The designs are provided as text files with names indicating their purpose. The experimental data is provided as a CSV file which contains a header with measurement information, followed by a collection of 50 shots across a collection of wavelengths, along with the computed median and mean across shots.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

TRISO-form HALEU-fueled Experiment for Transport Applications (THETA) Serpent 2 Benchmarking

A number of new advanced reactor designs use TRI-structural ISO-topic (TRISO) fuel. The uranium enrichment for the TRISO fuel is in the 5- to 19.75-weight percent range, which is considered High-Assay Low Enriched Uranium (HALEU). The Nuclear Regulatory Commission and vendors both have a strong interest in benchmarking these advanced reactor designs and having confidence in the predictive capabilities of the computer codes used to model them. In the recent years, Serpent 2 gained several capabilities focused on the reactor applications making it attractive for many reactor vendors.

Critical Experiment↗

Stability optimization of energetic particle driven modes in nuclear fusion devices: the FAR3d gyro-fluid code

The development of reduced models provide efficient methods that can be used to perform short term experimental data analysis or narrow down the parametric range of more sophisticated numerical approaches. Reduced models are derived by simplifying the physics description with the goal of retaining only the essential ingredients required to reproduce the phenomena under study. This is the role of the gyro-fluid code FAR3d, dedicated to analyze the linear and nonlinear stability of Alfvén Eigenmodes (AE), Energetic Particle Modes (EPM) and magnetic-hydrodynamic modes as pressure gradient driven mode (PGDM) and current driven modes (CDM) in nuclear fusion devices. Such analysis is valuable for improving the plasma heating efficiency and confinement; this can enhance the overall device performance. The present review is dedicated to a description of the most important contributions of the FAR3d code in the field of energetic particles (EP) and AE/EPM stability. FAR3d is used to model and characterize the AE/EPM activity measured in fusion devices as LHD, JET, DIII-D, EAST, TJ-II and Heliotron J. In addition, the computational efficiency of FAR3d facilitates performing massive parametric studies leading to the identification of optimization trends with respect to the AE/EPM stability. This can aid in identifying operational regimes where AE/EPM activity is avoided or minimized. This technique is applied to the analysis of optimized configurations with respect to the thermal plasma parameters, magnetic field configuration, external actuators and the effect of multiple EP populations. In addition, the AE/EPM saturation phase is analyzed, taking into account both steady-state phases and bursting activity observed in LHD and DIII-D devices. The nonlinear calculations provide: the induced EP transport, the generation of zonal structures as well as the energy transfer towards the thermal plasma and between different toroidal/helical families. Finally, FAR3d is used to forecast the AE/EPM stability in operational scenarios of future devices as ITER, CFETR, JT60SA and CFQS as well as possible approaches to optimization with respect to variations in the most important plasma parameters.

Alfv én Eigenmodes↗

Two-temperature effects in Hall-MHD simulations of the HIT-SI experiment

A two-temperature Hall-magnetohydrodynamic (Hall-MHD) model, which evolves the electron and ion temperatures separately, is implemented in the PSI-Tet 3D MHD code and used to model plasma dynamics in the Helicity Injected Torus–Steady Inductive (HIT-SI) experiment. The two-temperature model is utilized for HIT-SI simulations in both the PSI-Tet and NIMROD codes at a number of different injector frequencies in the 14.5–68.5 kHz range. At all frequencies, the NIMROD two-temperature model results in increased toroidal current, lower chord-averaged density, higher average temperatures, outward radial shift of the current centroid, and axial symmetrization of the current centroid, relative to the single-temperature NIMROD simulations. The two-temperature PSI-Tet model illustrates similar trends, but at high frequency operation, it exhibits lower electron temperature, smaller toroidal current, and decreased axial symmetrization with respect to the single-temperature PSI-Tet model. With all models, average temperatures and toroidal currents increase with the injector frequency. Power balance and heat fluxes to the wall are calculated for the two-temperature PSI-Tet model and illustrate considerable viscous and compressive heating, particularly at high injector frequency. Parameter scans are also presented for artificial diffusivity, wall temperature, and density. Both artificial diffusivity and the density boundary condition significantly modify the plasma density profiles, leading to larger average temperatures, toroidal current, and relative density fluctuations at low densities. As a result, a low density simulation achieves sufficiently high current gain (G > 5) to generate significant volumes of closed flux lasting 1–2 injector periods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multi-species collisions for delta-f gyrokinetic simulations: Implementation and verification with GENE

we report a multi-species linearized collision operator based on the model developed by Sugama et al. has been implemented in the nonlinear gyrokinetic code, GENE. Such a model conserves particles, momentum, and energy to machine precision, and is shown to have negative definite free energy dissipation characteristics, satisfying Boltzmann’s H-theorem, including for realistic mass ratio. Finite Larmor Radius (FLR) effects have also been implemented into the local version of the code. For the global version of the code, the collision operator has been developed to allow for block-structured velocity space grids, allowing for computationally tractable collisional global simulations. The validity of the collision operator has been demonstrated by relaxation and conservation tests, as well as appropriate benchmarks. The newly implemented operator shall be used in future simulations to study magnetically confined fusion plasma turbulence and transport in more extreme regions with higher collisionality.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multi-species collisions for delta-f gyrokinetic simulations: Implementation and verification with GENE

A multi-species linearized collision operator based on the model developed by Sugama et al. has been implemented in the nonlinear gyrokinetic code, GENE. Such a model conserves particles, momentum, and energy to machine precision, and is shown to have negative definite free energy dissipation characteristics, satisfying Boltzmann’s H-theorem, including for realistic mass ratio. Finite Larmor Radius (FLR) effects have also been implemented into the local version of the code. For the global version of the code, the collision operator has been developed to allow for block-structured velocity space grids, allowing for computationally tractable collisional global simulations. The validity of the collision operator has been demonstrated by relaxation and conservation tests, as well as appropriate benchmarks. The newly implemented operator shall be used in future simulations to study magnetically confined fusion plasma turbulence and transport in more extreme regions with higher collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗