Search NASA⌕ Search

SEARCH · Search NASA

Results for “Code Robustness”

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

DSD2D-FLS 2010: Bdzil's 2010 DSD Code Base; Computing tb and Dn with Edits to Reduce the Noise in the Dn Field Near HE Boundaries

The full level-set function code, DSD3D, is fully described in LA-14336 (2007) [1]. This ASCI-supported, DSD code project was the last such LANL DSD code project that I was involved with before my retirement in 2007. My part in the project was to design and build the core DSD3D solver, which was to include a robust DSD boundary condition treatment. A robust boundary condition treatment was required, since for an important local “customer,” the only description of the explosives’ boundary was through volume fraction data. Given this requirement, the accuracy issues I had encountered with our “fast-tube,” narrowband, DSD2D solver, and the difficulty we had building an efficient MPI-parallel version of the narrowband DSD2D, I decided DSD3D should be built as a full level-set function code, using a totally local DSD boundary condition algorithm for the level-­set function, phi, which did not rely on the gradient of the level-­set function being one, |grad(phi)| = 1. The narrowband DSD2D solver was built on the assumption that |grad(phi)| could be driven to one, and near the boundaries of the explosive this condition was not being satisfied. Since the narrowband is typically no more than10*dx wide, narrowband methods are discrete methods with a fixed, non-­resolvable error, where the error is related to the thickness of the band: the narrower the band the larger the errors. Such a solution represents a discrete approximation to the true solution and does not limit to the solution of the underlying PDEs under grid resolution.The full level-­set function code, DSD3D, is fully described in LA-14336 (2007) [1]. This ASCI-­supported, DSD code project was the last such LANL DSD code project that I was involved with before my retirement in 2007. My part in the project was to design and build the core DSD3D solver, which was to include a robust DSD boundary condition treatment. A robust boundary condition treatment was required, since for an important local “customer,” the only description of the explosives’ boundary was through volume fraction data. Given this requirement, the accuracy issues I had encountered with our “fast-­tube,” narrowband, DSD2D solver, and the difficulty we had building an efficient MPI-parallel version of the narrowband DSD2D, I decided DSD3D should be built as a full level-­set function code, using a totally local DSD boundary condition algorithm for the level-­set function, phi, which did not rely on the gradient of the level-­set function being one, |grad(phi)| = 1. The narrowband DSD2D solver was built on the assumption that |grad(phi)| could be driven to one, and near the boundaries of the explosive this condition was not being satisfied. Since the narrowband is typically no more than10*dx wide, narrowband methods are discrete methods with a fixed, non-resolvable error, where the error is related to the thickness of the band: the narrower the band the larger the errors. Such a solution represents a discrete approximation to the true solution and does not limit to the solution of the underlying PDEs under grid resolution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Localizing resonant magnetic perturbations for edge localized mode control in KSTAR

An external 3D magnetic perturbation typically drives a resonant response at the rational surfaces from the core to the edge of tokamak plasmas, due to strong mode coupling and amplification. In this paper, we present a method to isolate the edge from core resonant fields using the ideal perturbed equilibrium code and to design an edge-localized resonant magnetic perturbation (RMP) for effective edge localized mode (ELM) control. A robust feature of the edge-localized RMP is the curtailed response to the field at the low-field-side (LFS) midplane, as opposed to typical RMPs which strongly resonate with the LFS fields. This emphasizes the importance of off-midplane coils to improve ELM control without provoking a large core response that could lead to devastating instabilities. The conceptual design of new ELM control coils based on the edge-localized RMP in KSTAR shows how this new insight can be utilized to enhance the efficiency of our ELM suppression capabilities. Simple window-pane coils matching the edge-localized resonant mode structure substantially expand in the ELM suppression window beyond the existing coil. Further optimization using the flexible optimized coils using space-curves code leads to additional enhancement in the edge-localized control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Robust Mercury Methylation across Diverse Methanogenic Archaea

ABSTRACT Methylmercury (MeHg) production was compared among nine cultured methanogenic archaea that contain hgcAB , a gene pair that codes for mercury (Hg) methylation. The methanogens tested produced MeHg at inherently different rates, even when normalized to growth rate and Hg availability. Eight of the nine tested were capable of MeHg production greater than that of spent- and uninoculated-medium controls during batch culture growth. Methanococcoides methylutens , an hgcAB + strain with a fused gene pair, was unable to produce more MeHg than controls. Maximal conversion of Hg to MeHg through a full batch culture growth cycle for each species (except M. methylutens ) ranged from 2 to >50% of the added Hg(II) or between 0.2 and 17 pmol of MeHg/mg of protein. Three of the species produced >10% MeHg. The ability to produce MeHg was confirmed in several hgcAB + methanogens that had not previously been tested ( Methanocella paludicola SANAE, Methanocorpusculum bavaricum , Methanofollis liminatans GKZPZ, and Methanosphaerula palustris E1-9c). Maximal methylation was observed at low sulfide concentrations (<100 μM) and in the presence of 0.5 to 5 mM cysteine. For M. hollandica , the addition of up to 5 mM cysteine enhanced MeHg production and cell growth in a concentration-dependent manner. As observed for bacterial Hg methylators, sulfide inhibited MeHg production. An initial evaluation of sulfide and thiol impacts on bioavailability showed methanogens responding to Hg complexation in the same way as do Deltaproteobacteria . The mercury methylation rates of several methanogens rival those of the better-studied Hg-methylating sulfate- and iron-reducing Deltaproteobacteria . IMPORTANCE Archaea , specifically methanogenic organisms, play a role in mercury methylation in nature, but their global importance to MeHg production and the subsequent risk to ecosystems are not known. Methanogenesis has been linked to Hg methylation in several natural habitats where methylmercury production incurs risk to people and ecosystems, including rice paddies and permafrost. In this study, we confirm that most methanogens carrying the hgcAB gene pair are capable of Hg methylation. We found that methylation rates vary inherently among hgcAB + methanogens but that several species are capable of MeHg production at rates that rival those of the better-know Hg-methylating sulfate- and iron-reducing bacteria. Methanogens may need to be considered equally with sulfate and iron reducers in evaluations of MeHg production in nature.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Computing the Properties of Matter with Leadership Computing Resources (Closeout Report for DE-SC0018121)

In order to add more capabilities to Halide, we have designed a new framework called Tiramisu and integrated this framework into Halide. Since Tiramisu enables Halide to target heterogeneous architectures, our development efforts have been refocused on Tiramisu. Most high-performance computer systems today are complex and increasingly heterogeneous; they may have CPUs, GPUs and FPGAs. Achieving best performance requires taking full advantage of all these different architectures. To address this issue, we have designed Tiramisu, an optimization framework that enables Halide (and other DSLs) to target heterogeneous architectures. Tiramisu is an optimization framework that takes as input a high level, architecture-independent representation of code and a set of scheduling and data mapping commands that guide code transformation. The input can either be generated by a domain-specific language (DSL) compiler such as Halide or directly written by a programmer. Tiramisu then applies the user-specified code and data-layout transformations and generates an architecture-specific, low-level intermediate representation (IR) that takes advantage of modern architectural features such as multicore parallelism, non-uniform memory (NUMA) hierarchies, clusters, and accelerators like GPUs and FPGAs. We integrated Tiramisu within Halide and implemented a representative set of benchmarks to evaluate this integration. Tiramisu is now open source and is available for public use (http://tiramisu-compiler.org/). A paper about Tiramisu was published, it shows that Tiramisu extends Halide with many new capabilities and that Tiramisu can generate efficient code for multicores, GPUs, FPGAs and distributed heterogeneous systems. The performance of code generated by the Tiramisu backends matches or exceeds hand optimized reference implementations. For example, the multicore backend matches the highly optimized Intel MKL library on many kernels and shows speedups reaching 4x over the original Halide. In addition to making Tiramisu more robust, we have used Tiramisu to implement a set of representative tensor operation for constructing baryon building blocks required for multi baryon contractions in LQCD. In order to implement this code, we needed to generalize Tiramisu in two ways: first we needed to support indirect array accesses, and second, we needed to add support for complex numbers to Tiramisu. The code generated by Tiramisu is 6x faster than the reference code. Our efforts towards an MPI based multi-node version of tiramisu have matured and the resulting code scales well on multiple nodes (tests up to 512 KNL nodes have been undertaken).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

The DESI N -body simulation project – I. Testing the robustness of simulations for the DESI dark time survey

ABSTRACT Analysis of large galaxy surveys requires confidence in the robustness of numerical simulation methods. The simulations are used to construct mock galaxy catalogues to validate data analysis pipelines and identify potential systematics. We compare three N-body simulation codes, abacus, gadget-2, and swift, to investigate the regimes in which their results agree. We run N-body simulations at three different mass resolutions, 6.25 × 108, 2.11 × 109, and 5.00 × 109 h−1 M⊙, matching phases to reduce the noise within the comparisons. We find systematic errors in the halo clustering between different codes are smaller than the Dark Energy Spectroscopic Instrument (DESI) statistical error for $s\ \gt\ 20\ h^{-1}$ Mpc in the correlation function in redshift space. Through the resolution comparison we find that simulations run with a mass resolution of 2.1 × 109 h−1 M⊙ are sufficiently converged for systematic effects in the halo clustering to be smaller than the DESI statistical error at scales larger than $20\ h^{-1}$ Mpc. These findings show that the simulations are robust for extracting cosmological information from large scales which is the key goal of the DESI survey. Comparing matter power spectra, we find the codes agree to within 1 per cent for k ≤ 10 h Mpc−1. We also run a comparison of three initial condition generation codes and find good agreement. In addition, we include a quasi-N-body code, FastPM, since we plan use it for certain DESI analyses. The impact of the halo definition and galaxy–halo relation will be presented in a follow-up study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

From plane waves to local Gaussians for the simulation of correlated periodic systems

Here, we present a simple, robust, and black-box approach to the implementation and use of local, periodic, atom-centered Gaussian basis functions within a plane wave code, in a computationally efficient manner. The procedure outlined is based on the representation of the Gaussians within a finite bandwidth by their underlying plane wave coefficients. The core region is handled within the projected augment wave framework, by pseudizing the Gaussian functions within a cutoff radius around each nucleus, smoothing the functions so that they are faithfully represented by a plane wave basis with only moderate kinetic energy cutoff. To mitigate the effects of the basis set superposition error and incompleteness at the mean-field level introduced by the Gaussian basis, we also propose a hybrid approach, whereby the complete occupied space is first converged within a large plane wave basis, and the Gaussian basis used to construct a complementary virtual space for the application of correlated methods. We demonstrate that these pseudized Gaussians yield compact and systematically improvable spaces with an accuracy comparable to their non-pseudized Gaussian counterparts. A key advantage of the described method is its ability to efficiently capture and describe electronic correlation effects of weakly bound and low-dimensional systems, where plane waves are not sufficiently compact or able to be truncated without unphysical artifacts. We investigate the accuracy of the pseudized Gaussians for the water dimer interaction, neon solid, and water adsorption on a LiH surface, at the level of second-order Møller–Plesset perturbation theory.

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

Scaling Studies for Advanced High Temperature Reactor Concepts, Final Technical Report: October 2014—December 2017

Computer simulations of nuclear reactor thermal-hydraulic phenomena are often used in the design and licensing of nuclear reactor systems. In order to assess the accuracy of these computer simulations, computer codes and methods are often validated against experimental data. This experimental data must be of sufficiently high quality in order to conduct a robust validation exercise. In addition, this experimental data is generally collected at experimental facilities that are of a smaller scale than the reactor systems that are being simulated due to cost considerations. Therefore, smaller scale test facilities must be designed and constructed in such a fashion to ensure that the prototypical behavior of a particular nuclear reactor system is preserved. The work completed through this project has resulted in scaling analyses and conceptual design development for a test facility capable of collecting code validation data for the following high temperature gas reactor systems and events— 1. Passive natural circulation core cooling system, 2. pebble bed gas reactor concept, 3. General Atomics Energy Multiplier Module reactor, and 4. prismatic block design steam-water ingress event. In the event that code validation data for these systems or events is needed in the future, significant progress in the design of an appropriate integral-type test facility has already been completed as a result of this project. Where applicable, the next step would be to begin the detailed design development and material procurement. As part of this project applicable scaling analyses were completed and test facility design requirements developed. Conceptual designs were developed for the implementation of these design requirements at the Oregon State University (OSU) High Temperature Test Facility (HTTF). The original HTTF is based on a ¼-scale model of a high temperature gas reactor concept with the capability for both forced and natural circulation flow through a prismatic core with an electrical heat source. The peak core region temperature capability is 1400°C. As part of this project, an inventory of test facilities that could be used for these experimental programs was completed. Several of these facilities showed some promise, however, upon further investigation it became clear that only the OSU HTTF had the power and/or peak temperature limits that would allow for the experimental programs envisioned herein. Thus the conceptual design and feasibility study development focused on examining the feasibility of configuring the current HTTF to collect validation data for these experimental programs. In addition to the scaling analyses and conceptual design development, a test plan was developed for the envisioned modified test facility. This test plan included a discussion on an appropriate shakedown test program as well as the specific matrix tests. Finally, a feasibility study was completed to determine the cost and schedule considerations that would be important to any test program developed to investigate these designs and events.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

MDSplus quality improvement project

MDSplus is a data acquisition and analysis system used worldwide predominantly in the fusion research community. Development began 29 years ago on the OpenVMS operating system. Since that time there have been many new features added and the code has been ported to many different operating systems. There have been contributions to the MDSplus development from the fusion community in the way of feature suggestions, feature implementations, documentation and porting to different operating systems. The bulk of the development and support of MDSplus, however, has been provided by a relatively small core developer group of three or four members. Given the size of the development team and the large number of users much more effort was focused on providing new features for the community than on keeping the underlying code and documentation up to date with the evolving software development standards. To ensure that MDSplus will continue to provide the needs of the community in the future, the MDSplus development team along with other members of the MDSplus user community has commenced on a major quality improvement project. The planned improvements include changes to software build scripts to better use GNU Autoconf and Automake tools, refactoring many of the source code modules using new language features available in modern compilers, using GNU MinGW-w64 to create MS Windows distributions, migrating to a more modern source code management system, improvement of source documentation as well as improvements to the www.mdsplus.org web site documentation and layout, and the addition of more comprehensive test suites to apply to MDSplus code builds prior to releasing installation kits to the community. This paper should lead to a much more robust product and establish a framework to maintain stability as more enhancements and features are added. Finally, this paper will describe these efforts that are either in progress or planned for the near future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

UQpy v4.1: Uncertainty quantification with Python

This paper presents the latest improvements introduced in Version 4 of the UQpy, Uncertainty Quantification with Python, library. In the latest version, the code was restructured to conform with the latest Python coding conventions, refactored to simplify previous tightly coupled features, and improve its extensibility and modularity. To improve the robustness of UQpy, software engineering best practices were adopted. A new software development workflow significantly improved collaboration between team members, and continuous integration and automated testing ensured the robustness and reliability of software performance. Continuous deployment of UQpy allowed its automated packaging and distribution in system agnostic format via multiple channels, while a Docker image enables the use of the toolbox regardless of operating system limitations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Demonstration of x-ray fluorescence spectroscopy as a sensitive temperature diagnostic for high-energy-density physics experiments

We present the use of x-ray fluorescence spectroscopy (XFS) to a sensitive temperature diagnostic in shocked foams at temperatures of 30–75 eV. Cobalt-doped foams were shock compressed using a planar drive at the OMEGA laser facility and photo-pumped with a Zn He⁢𝛼 x-ray source. Analysis of the resulting cobalt 𝐾⁢𝛽 x-ray fluorescence spectra using collisional radiative codes allows the temperature to be determined in the shocked foams. Furthermore, this method provides a sensitive and robust technique to determine temperatures in high-energy-density physics experiments in the tens of electronvolts temperature range. In these experiments, we find that radiation hydrodynamic simulations predict a lower temperature in the shocked foams compared to analysis of the XFS data using collisional radiative models. Although additional experiments with an independent temperature diagnostic to absolutely calibrate XFS spectra for these conditions will be required to resolve this discrepancy, these results demonstrate the excellent temperature sensitivity of XFS spectra for high-energy-density physics experiments.

Atomic spectra↗
DOE OSTI2015

Exploration of direct neutron capture with covariant density functional theory inputs

Predictions of direct neutron capture are of vital importance for simulations of nucleosynthesis in supernovae, merging neutron stars, and other astrophysical environments. Here, we calculated direct capture cross sections using nuclear structure information obtained from a covariant density functional theory as input for the fresco coupled reaction channels code. We investigated the impact of pairing, spectroscopic factors, and optical potentials on our results to determine a robust method to calculate cross sections of direct neutron capture on exotic nuclei. Our predictions agree reasonably well with experimental cross section data for the closed shell nuclei 16 O and 48 Ca, and for the exotic nucleus 36 S. We then used this approach to calculate the direct neutron capture cross section on the doubly magic unstable nucleus 132 Sn which is of interest for the astrophysical r-process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Shocks in the Scrape-off Layer: Hydrodynamical Aspects (Summary) [Slides]

Hydrodynamical model provides a robust description of mass, momentum, and energy transfer. It can be used for benchmarking hydrodynamical modules of SOL codes. It clarifies the role of boundary conditions at the end-plates. It provides a “geometrical” model for partition of power flows between two divertor legs, both for the steady-state and pulsed flows. It indicates possible presence of shocks in the divertor legs. A framework for finding the shock location and strength has been developed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

The Fungal Cell Death Regulator czt-1 Is Allelic to acr-3

Fungal infections have far-reaching implications that range from severe human disease to a panoply of disruptive agricultural and ecological effects, making it imperative to identify and understand the molecular pathways governing the response to antifungal compounds. In this context, CZT-1 (cell death-activated zinc cluster transcription factor) functions as a master regulator of cell death and drug susceptibility in Neurospora crassa. Here we provide evidence indicating that czt-1 is allelic to acr-3, a previously described locus that we now found to harbor a point mutation in its coding sequence. This nonsynonymous amino acid substitution in a low complexity region of CZT-1/ACR-3 caused a robust gain-of-function that led to reduced sensitivity to acriflavine and staurosporine, and increased expression of the drug efflux pump abc-3. Thus, accumulating evidence shows that CZT-1 is an important broad regulator of the cellular response to various antifungal compounds that appear to share common molecular targets.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Evaluation of Shot-to-Shot In-Nozzle Flow Variations in a Heavy-Duty Diesel Injector Using Real Nozzle Geometry

Cyclic variability in internal combustion engines (ICEs) arises from multiple concurrent sources, many of which remain to be fully understood and controlled. This variability can, in turn, affect the behavior of the engine resulting in undesirable deviations from the expected operating conditions and performance. Shot-to-shot variation during the fuel injection process is strongly suspected of being a source of cyclic variability. This study focuses on the shot-to-shot variability of injector needle motion and its influence on the internal nozzle flow behavior using diesel fuel. High-speed x-ray imaging techniques have been used to extract high-resolution injector geometry images of the sac, orifices, and needle tip that allowed the true dynamics of the needle motion to emerge. These measurements showed high repeatability in the needle lift profile across multiple injection events, while the needle radial displacement was characterized by a much higher degree of randomness. A robust and previously validated computational setup from the authors’ research group using a commercial computational fluid dynamics (CFD) code has been adapted for the eight-hole heavy-duty common-rail diesel injector used for the measurements. The simulation results obtained using the x-ray-scanned geometry have been validated against available experimental data of mass flow rate at the nozzle exit. The fuel mass flow rate has then been analyzed at three different injection pressures that cover the conditions at which the injector typically operates during normal engine operation. Finally, the average off-axis motion of the needle has been perturbed (based on the variability found in the experimental measurements) to generate three new cases that present different amplitude and phasing of the radial displacement with respect to the baseline average motion. Furthermore, this revealed the effects of off-axis motion on shot-to-shot and orifice-to-orifice variations.

42 ENGINEERING↗
DOE OSTI2016

THE AGORA HIGH-RESOLUTION GALAXY SIMULATIONS COMPARISON PROJECT. II. ISOLATED DISK TEST

Using an isolated Milky Way-mass galaxy simulation, we compare results from nine state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt–Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly formed stellar clump mass functions show more significant variation (difference by up to a factor of ∼3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low-density region, and between more diffusive and less diffusive schemes in the high-density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Warthog: A MOOSE-Based Application for the Direct Code Coupling of BISON and PROTEUS

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program from the Department of Energy's Office of Nuclear Energy provides a robust toolkit for the modeling and simulation of current and future advanced nuclear reactor designs. This toolkit provides these technologies organized across product lines: two divisions targeted at fuels and end-to-end reactor modeling, and a third for integration, coupling, and high-level workflow management. The Fuels Product Line and the Reactor Product line provide advanced computational technologies that serve each respective field well, however, their current lack of integration presents a major impediment to future improvements of simulation solution fidelity. There is a desire for the capability to mix and match tools across Product Lines in an effort to utilize the best from both to improve NEAMS modeling and simulation technologies. This report details a new effort to provide this Product Line interoperability through the development of a new application called Warthog. This application couples the BISON Fuel Performance application from the Fuels Product Line and the PROTEUS Core Neutronics application from the Reactors Product Line in an effort to utilize the best from all parts of the NEAMS toolkit and improve overall solution fidelity of nuclear fuel simulations. To achieve this, Warthog leverages as much prior work from the NEAMS program as possible, and in doing so, enables interoperability between the disparate MOOSE and SHARP frameworks, and the libMesh and MOAB mesh data formats. This report describes this work in full. We begin with a detailed look at the individual NEAMS framework technologies used and developed in the various Product Lines, and the current status of their interoperability. We then introduce the Warthog application: its overall architecture and the ways it leverages the best existing tools from across the NEAMS toolkit to enable BISON-PROTEUS integration. Furthermore, we show how Warthog leverages a tool known as DataTransferKit to seamlessly enable the transfer for solution data between disparate frameworks and mesh formats. To end, we demonstrate tests for the direct software coupling of BISON and PROTEUS using Warthog, and discuss current impediments and solutions to the construction of physically realistic input models for this coupled BISON-PROTEUS system.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

A NEW GENERATION OF LOS ALAMOS OPACITY TABLES

We present a new, publicly available set of Los Alamos OPLIB opacity tables for the elements hydrogen through zinc. Our tables are computed using the Los Alamos ATOMIC opacity and plasma modeling code, and make use of atomic structure calculations that use fine-structure detail for all the elements considered. Our equation of state model, known as ChemEOS, is based on the minimization of free energy in a chemical picture and appears to be a reasonable and robust approach to determining atomic state populations over a wide range of temperatures and densities. In this paper we discuss in detail the calculations that we have performed for the 30 elements considered, and present some comparisons of our monochromatic opacities with measurements and other opacity codes. We also use our new opacity tables in solar modeling calculations and compare and contrast such modeling with previous work.

79 ASTRONOMY AND ASTROPHYSICS↗