Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “L CODES”

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 487 records · Page 27

NCSU RDMA data

NCSU RDMA Data Level b1: QC checks applied to measurements Data Format: CSV Description: See "instrument description Site: Houston, TX; Tracking Aerosol Convection interactions ExpeRiment (HOU) Location: Houston, TX; AMF1 (main site for TRACER) Facility Code: M1 Category: Aerosol Properties Data Type: PI Data Source Instrument/Data: Nano-Scanning Mobility Particle (SMPS); A Radial Differential Mobility Analyzer (RDMA) coupled with 1 Condensation Particle Counter Start Date: 2022-06-01 End Date: 2022-9-26 Contact PI: Markus Petters (mdpetter@ncsu.edu) Funding Source: DOE ASR award US Department of Energy, Office of Science, Biological and Environment Research (grant no. DE-SC 0021074) Instrument description The NCSU RDMA was operated at a sheath-to-sample flow ratio of 5:1.5 L min−1. The RDMA was configured to scan from 5 to 55 nm. It was located into the temperature-controlled trailer adjacent to the NCSU Flux Tower to observe size distributions of the aerosols.The sample line was dried with three silica-gel driers in series, and then neutralized with X-ray neutralizer. Please contact mdpetter@ncsu.edu for further information.

54 ENVIRONMENTAL SCIENCES↗

Accelerating nuclear fuel development and qualification: Modeling and simulation integrated with separate-effects testing

In this work, an approach to transform and accelerate nuclear fuel development and qualification is outlined. The approach exploits advanced modeling and simulation at the outset to inform constituent and system selection and to enable integral fuel performance analyses. Analyses using these tools identify and prioritize the most important fuel performance parameters and phenomena for subsequent targeted characterization with separate-effects tests. Separate-effects testing spans out-of-pile and in-pile tests and is meant to iterate with and inform engineering-scale integral fuel performance analyses throughout the development process. Exercising this cycle in an agile fashion will increase confidence in the integral fuel performance predictions while reducing uncertainties. This process sets the stage for executing a much more limited set of well-defined integral irradiation tests designed to validate engineering-scale fuel performance codes and to confirm the performance and safety of the fuel system under prototypic conditions. This approach will reduce the time for development and qualification of a new fuel system, and it will also reduce associated costs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Adaptive finite differencing in high accuracy electronic structure calculations

Abstract A multi-order Adaptive Finite Differencing (AFD) method is developed for the kinetic energy operator in real-space, grid-based electronic structure codes. It uses atomic pseudo orbitals produced by the corresponding pseudopotential codes to optimize the standard finite difference (SFD) operators for improved precision. Results are presented for a variety of test systems and Bravais lattice types, including the well-known Δ test for 71 elements in the periodic table, the Mott insulator NiO, and borax decahydrate, which contains covalent, ionic, and hydrogen bonds. The tests show that an 8th-order AFD operator leads to the same average Δ value as that achieved by plane-wave codes and is typically far more accurate and has a much lower computational cost than a 12th-order SFD operator. The scalability of real-space electronic calculations is demonstrated for a 2016-atom NiO cell, for which the computational time decreases nearly linearly when scaled from 18 to 144 CPU-GPU nodes.

Briggs, E. L. (ORCID:0000000343983492)↗

Design study of a combined interferometer and polarimeter for a high-field, compact tokamak

This article is the first design study of a combined interferometer and polarimeter on a compact, high-field, high-density, net-energy tokamak. Recent advances in superconducting technology have made possible designs for compact, high magnetic field fusion power plants, such as ARC [Sorbom et al., Fusion Eng. Des. 100, 378 (2015)], and experiments, such as SPARC [Greenwald et al., PSFC Report No. RR-18-2 (2018)]. These new designs create both challenges and opportunities for plasma diagnostics. The diagnostic proposed in this work, called InterPol, takes advantage of unique opportunities provided by high magnetic field and density to measure both line-averaged density and poloidal magnetic field with a single set of CO 2 and quantum cascade lasers. These measurements will be used for fast density feedback control, constraint of density and safety factor profiles, and density fluctuation measurements. Synthetic diagnostic testing using a model machine geometry, called MQ1 (Mission Q ≥ 1), and profiles simulated with Tokamak Simulation Code indicate that InterPol will be able to measure steady state density and poloidal magnetic field, as well as fluctuations caused by toroidal Alfven eigenmodes and other phenomena on a high-field compact tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cross-code comparison of the impact of the fill tube on high yield implosions on the National Ignition Facility

Fill tubes are used to inject deuterium and tritium fuel into inertial confinement fusion capsules fielded on the National Ignition Facility. These fill tubes have been shown to have a detrimental effect on capsule performance, primarily by introducing a low-density pathway into the central fuel region that enables the jetting of ablation material into the hot spot. Due to the complexity of the highly nonlinear flow associated with the fill tube and the challenge of diagnosing the evolution of the fill tube jet late in the implosion experiments, the uncertainty in how this perturbation source evolves is great. In this work, we report on the results of a detailed code comparison performed to understand uncertainties in computational modeling of the impact of fill tubes on implosion performance. The study employed two radiation-hydrodynamics codes, HYDRA and xRAGE, which employ very different meshing strategies and hydrodynamics solvers, as well as two radiation transport methodologies, discrete ordinates and multi-group diffusion. Our results demonstrate generally good agreement between codes through most of the implosion although they indicate sensitivity to opacity averaging methods. Late in the implosion, differences arise in the distribution and amount of contaminant although these differences have a remarkably small impact on the amount of yield reduction due to the fill tube. While these results demonstrate sensitivity in fill tube modeling to algorithmic choices, the observed differences between codes are small relative to known sensitivities due to expected variations in the fill tube geometry. Finally, we have developed a methodology for performing multi-group diffusion simulations that show good agreement with the more accurate discrete ordinates method

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Crown Pearl: a draft genome assembly of the European freshwater pearl mussel Margaritifera margaritifera (Linnaeus, 1758)

Abstract Since historical times, the inherent human fascination with pearls turned the freshwater pearl mussel Margaritifera margaritifera (Linnaeus, 1758) into a highly valuable cultural and economic resource. Although pearl harvesting in M. margaritifera is nowadays residual, other human threats have aggravated the species conservation status, especially in Europe. This mussel presents a myriad of rare biological features, e.g. high longevity coupled with low senescence and Doubly Uniparental Inheritance of mitochondrial DNA, for which the underlying molecular mechanisms are poorly known. Here, the first draft genome assembly of M. margaritifera was produced using a combination of Illumina Paired-end and Mate-pair approaches. The genome assembly was 2.4 Gb long, possessing 105,185 scaffolds and a scaffold N50 length of 288,726 bp. The ab initio gene prediction allowed the identification of 35,119 protein-coding genes. This genome represents an essential resource for studying this species’ unique biological and evolutionary features and ultimately will help to develop new tools to promote its conservation.

Gomes-dos-Santos, André (ORCID:0000000199734861)↗

A machine learning approach to galaxy properties: joint redshift–stellar mass probability distributions with Random Forest

We demonstrate that highly accurate joint redshift–stellar mass probability distribution functions (PDFs) can be obtained using the Random Forest (RF) machine learning (ML) algorithm, even with few photometric bands available. As an example, we use the Dark Energy Survey (DES), combined with the COSMOS2015 catalogue for redshifts and stellar masses. We build two ML models: one containing deep photometry in the griz bands, and the second reflecting the photometric scatter present in the main DES survey, with carefully constructed representative training data in each case. We validate our joint PDFs for 10 699 test galaxies by utilizing the copula probability integral transform and the Kendall distribution function, and their univariate counterparts to validate the marginals. Benchmarked against a basic set-up of the template-fitting code bagpipes, our ML-based method outperforms template fitting on all of our predefined performance metrics. In addition to accuracy, the RF is extremely fast, able to compute joint PDFs for a million galaxies in just under 6 min with consumer computer hardware. Such speed enables PDFs to be derived in real time within analysis codes, solving potential storage issues. As part of this work we have developed galpro 1, a highly intuitive and efficient python package to rapidly generate multivariate PDFs on-the-fly. galpro is documented and available for researchers to use in their cosmology and galaxy evolution studies.

79 ASTRONOMY AND ASTROPHYSICS↗

Metastable states from multinucleon excitations in 202 Tl and 203 Pb

The excited level structures of 202 Tl and 203 Pb, above the 7 + and 29/2 – isomers, respectively, have been studied. An isomer with I π = 20 + and T 1/2 = 215(10) μs has been established in 202 Tl, and the level scheme extended from I = 10 to 20ℏ with the placement of fifteen new transitions. In 203 Pb, the I π = 37/2 + state is established to be metastable, with T 1/2 = 2.5(3) ns. Levels in both nuclei arise from intrinsic excitations, with likely particle-hole character for the higher-lying states in 203 Pb. The 20 + isomer in 202 Tl is most likely associated with a $πh$ $^{–1}_{11/2}$ Ⓧ ν($i$ $^{–2}_{13/2}$, $f$ $^{–1}_{5/2}$) configuration, while the 37/2 + state in 203 Pb results from the excitation of five neutrons. Furthermore, calculations, using both an empirical approach and the oxbash code, have been performed to aid in the description of the excited level structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Montepy

Montepy is an object-oriented Python interface for MCNP input files. It excels at reading, editing, creating, and writing MCNP input files. The Monte Carlo N Particle (MCNP) code is a Monte Carlo transport solver for simulating many different types of radiation problems and is developed by Los Alamos National Laboratory. Montepy creates an object-oriented interface for the MCNP input file to allow for efficient automated modifications to large, complex simulations.

Gale, Micah [Idaho National Laboratory (INL), Idah↗

Vadose Zone Model for M Pond Area for Composite Analysis

The objectives of the vadose modeling for the updated Hanford Site composite analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates for the plateau to river (P2R) model, version 8.3 (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The following 16 radionuclides were selected for this modeling effort; carbon-14 (C-14), chlorine-36 (Cl-36), tritium (H-3), iodine-129 (I-129), neptunium-237 (Np-237), rhenium-187 (Re-187), strontium-90 (Sr-90), technetium-99 (Tc-99), uranium-232 (U-232), uranium-233 (U-233), uranium-234 (U-234), uranium-235 (U-235), uranium-236 (U-236), uranium-238 (U-238), radium-226 (Ra-226), and thorium-230 (Th-230). The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The parallel version of the Subsurface Transport Over Multiple Phases (STOMP1) simulator, officially named the exascale Subsurface Transport Over Multiple Phases (eSTOMP), is used to simulate flow and transport for the vadose models. The documentation for the STOMP code is comprehensive. The theoretical and numerical approaches applied in the STOMP code are documented in a published theory guide (PNNL-12030, STOMP Subsurface Transport Over Multiple Phases Version 2.0 Theory Guide). The code has undergone a rigorous verification procedure against analytical solutions, laboratory-scale experiments, and field-scale demonstrations. The application guide (PNNL-11216, STOMP Subsurface Transport Over Multiple Phases Application Guide) provides instructive examples in the application of the code to classical groundwater problems. The user’s guide (PNNL-15782, STOMP: Subsurface Transport Over Multiple Phases Version 4.0: User’s Guide) describes the general use, input file formatting, compilation, and execution of the code. The primary output of the vadose zone modeling is radionuclide transfer rates to the groundwater for input into the P2R model. The rates will be summed over the 100 by 100 m P2R grid cells that fall within the vadose zone model source domain.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Sierra/SolidMechanics 4.58 Verification Tests Manual

Presented in this document is a small portion of the tests that exist in the Sierra/SolidMechanics (Sierra/SM) verification test suite. Most of these tests are run nightly with the Sierra/SM code suite, and the results of the test are checked versus the correct analytical result. For each of the tests presented in this document, the test setup, a description of the analytic solution, and comparison of the Sierra/SM code results to the analytic solution is provided. Mesh convergence is also checked on a nightly basis for several of these tests. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems. Additional example problems are provided in the Sierra/SM Example Problems Manual. Note, many other verification tests exist in the Sierra/SM test suite, but have not yet been included in this manual.

97 MATHEMATICS AND COMPUTING↗

Understanding HPC Benchmark Performance on Intel Broadwell and Cascade Lake Processors

Hardware platforms in high performance computing are constantly getting more complex to handle even when considering multicore CPUs alone. Numerous features and configuration options in the hardware and the software environment that are relevant for performance are not even known to most application users or developers. Microbenchmarks, i.e., simple codes that fathom a particular aspect of the hardware, can help to shed light on such issues, but only if they are well understood and if the results can be reconciled with known facts or performance models. The insight gained from microbenchmarks may then be applied to real applications for performance analysis or optimization. In this paper we investigate two modern Intel x86 server CPU architectures in depth: Broadwell EP and Cascade Lake SP. We highlight relevant hardware configuration settings that can have a decisive impact on code performance and show how to properly measure on-chip and off-chip data transfer bandwidths. The new victim L3 cache of Cascade Lake and its advanced replacement policy receive due attention. Finally we use DGEMM, sparse matrix-vector multiplication, and the HPCG benchmark to make a connection to relevant application scenarios.

97 MATHEMATICS AND COMPUTING↗

Opacity Table Generation

Radiation opacity is a key property of material systems which describes the absorption and scattering of photons. Such processes are an important part of energy transport, diagnostics, and astrophysical phenomena. Recent experiments in iron opacity and continuum lowering have sparked renewed interest in and careful vetting of opacity models. Even within the local thermodynamic equilibrium approximation, opacity varies over a large parameter space: material, temperature, density, and energy. Opacity tables that span wide regimes of this parameter space are of prime importance to radiation-hydrodynamics codes. Accurate opacity calculations cannot be achieved inline during a radiation hydrodynamics simulation since they involve accounting for all possible electron transitions of all important atomic states of a thermal ensemble. Livermore’s current framework utilizes precalculated tables at many densities, temperatures, and energies. Last year, the opacity theory team completed an L2 milestone concerning the development of a new code, Opus, which served the dual purpose of creating a modern code infrastructure and version control as well as the basis for training new members of that team, a critical need as key members near or are passed retirement. For the current L2 milestone, we exercised and validated Opus on three elements: boron, carbon, and nitrogen, which have few enough electrons to allow a careful convergence study over many input parameters. Here, we report on the production of new tables, show their convergence with respect to important parameters, state some of their limitations, and document a mostly automated production process. In this section, we address the milestone completion criteria, found on the title page and the subsection headings. In section 2, we validate Opus against Tycho and compare with TabOp and Atomic. Section 3 details some of the physics implemented in Opus and relevant to this report. We then cover details of the density-temperature grid in section 4 and do a spectral comparison between different parameter settings and different codes in section 5. After summarizing these results in section 6, we give instructions how to run the opacity table scripts in appendix A and document the scripts and input/output files in appendix B.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Radially resolved active charge exchange measurements of the hydrogenic isotope fraction on DIII-D

Radially resolved hydrogenic isotope fraction measurement capabilities have been developed for DIII-D using the main-ion charge exchange recombination (MICER) spectroscopy system in preparation for mixed hydrogen and deuterium experiments. Constraints on the hydrogenic ion temperatures and velocities based on measurements of the impurity ion properties are required to accurately fit the spectrum. Corrections for cross sectional distortions, spatial smearing due to the halo, and a neoclassical offset between the impurity and hydrogenic toroidal rotation are applied to the constraints prior to fitting the MICER spectrum. Extensive atomic physics calculations have been performed using the FIDASIM code, which has recently been improved to allow simulations using mixtures of hydrogenic species. Overall, these results demonstrate that for the same plasma parameters, the Dα emission is 20%–30% brighter than Hα due to differences in rate coefficients associated with the different ion thermal velocities for the same temperature and therefore must be taken into consideration when calculating absolute densities. However, despite these differences, the absolute error when estimating the hydrogen isotope fraction [n H /(n H + n D )] by using the Hα radiance fraction [L Hα /(L Hα + L Dα )] is typically less than 5% due to the way the fraction is formed, making the radiance fraction a reasonably accurate estimate of the isotope fraction for most cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Generation and characterization of ultrabroadband polarization–frequency hyperentangled photons

We generate ultrabroadband photon pairs entangled in both polarization and frequency bins through an all-waveguided Sagnac source covering the entire optical C- and L-bands (1530–1625 nm). We perform comprehensive characterization of high-fidelity states in multiple dense wavelength-division multiplexed channels, achieving full tomography of effective four-qubit systems. Additionally, leveraging the inherent high dimensionality of frequency encoding and our electro-optic measurement approach, we demonstrate the scalability of our system to higher dimensions, reconstructing states in a 36-dimensional Hilbert space consisting of two polarization qubits and two frequency-bin qutrits. Our findings hold potential significance for quantum networking, particularly dense coding and entanglement distillation in wavelength-multiplexed quantum networks.

42 ENGINEERING↗

An experimentally informed design process for future inertial confinement fusion facilities

The achievement of ignition in the laboratory has renewed interest in defining the requirements for a future high-gain inertial confinement fusion (ICF) facility. Our best chance of predicting future ICF performance is with 3-D radiation hydrodynamic simulations that have been benchmarked against experimental data, but their high computational cost is prohibitive for use in practical design studies. We introduce a hierarchical approach where 3-D simulations are tuned to match experimental measurements and used to train 3-D degradation models in 1-D simulations allowing for accurate predictions over the entire OMEGA direct-drive database. A genetic algorithm was used in combination with the trained 1-D simulations to search for optimal direct-drive implosion designs at driver energies ranging from 20 kJ to 10 MJ. As the fidelity of 3-D codes improves, this approach will provide a viable experimentally informed tool for defining the next ICF facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nanosecond isomers and the evolution of collectivity in stable, even- A Hg isotopes

Isomeric states and associated collective structures have been studied up to high spin in 198,200,202 Hg using multinucleon transfer reactions and the Gammasphere array. A coupled rotational band, with possible four-quasiparticle character, is established in 198 Hg. Sequences built on two-quasiparticle, positive- and negative-parity levels are assigned to 202 Hg. New isomers in 202 Hg with I π = (7 – ) and (9 – ), and T 1/2 = 10.4(4) ns and 1.4(3) ns, respectively, have been identified. A half-life of 1.0(3) ns is established for the I π = 12 + state in 200 Hg. B(E2) values deduced from isomeric transitions in Hg isotopes indicate that, while collectivity near the ground state gradually diminishes from N = 112 to N = 124, it is found to increase for the 12 + and 9 – states up to N = 118, followed by a reduction for higher neutron numbers. Calculations using the ultimate cranker code provide insight into the variation of deformation with spin and allow for an understanding of observed band crossings. As a result, the evolution of collectivity with spin, and along the isotopic chain, is described.

190 ≤ A ≤ 219↗

Atmospheric neutrinos in next-generation xenon and argon dark matter experiments

In this work, we study the sensitivity of future xenon- and argon-based dark matter and neutrino detection experiments to low-energy atmospheric neutrinos. Not accounting for experimental backgrounds, the primary obstacle for identifying nuclear recoils induced by atmospheric neutrinos in xenon is the tail of the electron recoil distribution due to pp solar neutrinos. We use the NEST code to model the solar and atmospheric neutrino signals in a xenon detector and find that an exposure of 700 tonne-years will produce a 5σ detection of atmospheric neutrinos. We explore the effect of different detector properties and find that a sufficiently long electron lifetime is essential to the success of such a measurement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗