Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “H 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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

PETSc Users Manual (Rev. 3.13)

This manual describes the use of PETSc for the numerical solution of partial differential equations and related problems on high-performance computers. The Portable, Extensible Toolkit for Scientific Computation (PETSc) is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all message-passing communication. PETSc includes an expanding suite of parallel linear solvers, nonlinear solvers, and time integrators that may be used in application codes written in Fortran, C, C++, and Python. PETSc provides many of the mechanisms needed within parallel application codes, such as parallel matrix and vector assembly routines. The library is organized hierarchically, enabling users to employ the level of abstraction that is most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users. PETSc is a sophisticated set of software tools; as such, for some users it initially has a much steeper learning curve than a simple subroutine library. In particular, for individuals without some computer science background, experience programming in C, C++, python, or Fortran and experience using a debugger such as gdb or dbx, it may require a significant amount of time to take full advantage of the features that enable efficient software use. However, the power of the PETSc design and the algorithms it incorporates may make the efficient implementation of many application codes simpler than \rolling them" yourself.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Quantifying Moss Response to Metal Contaminant Exposure Using Laser-Induced Fluorescence

Tracing sources of contamination, including potentially toxic elements (PTEs), has historically been achieved through sampling and analysis of soil or biota, which are labor-intensive, costly, and destructive methods. Thus, availability of a non-destructive in situ remote sensing method for monitoring metals deposited in biota is of great interest. Laser-induced fluorescence (LIF) is an emerging spectroscopic and imaging technique that documents changes in molecular energy level in plants as a biological response to metal contamination. For a proof-of-concept study and preliminary experiment, moss was selected for experimentation due to its long history of use in tracing atmospheric deposition of PTEs. Consecutive treatments of copper chloride (CuCl2) were administered to three moss samples, simulating wet deposition every 48 h over 10 days until reaching cumulative Cu concentrations of 2.690 to 8.075 μmol/cm2. While these Cu amounts are above environmentally relevant concentrations, they allowed the best conditions for testing and fine tuning of the imaging and data processing protocols presented in this paper. Moss fluorescence was induced using both 532 nm green and 355 nm UV lasers. A CMOS camera captured images of the LIF response, and red–green–blue (RGB) decimal code values were extracted for each pixel in the images, and pixel densities of color channels from treated and untreated moss samples were compared. Results show a shift towards lower color decimal codes corresponding to increased Cu concentration. We developed and contrasted multiple quantitative analyses of color distributions and demonstrated that LIF shows great promise for remote sensing of Cu accumulation in moss at μmol/cm2 levels. Though currently, the method would be limited to highly toxic sites, it illustrates the possibility and provides a framework for development of higher-sensitivity methods to detect nmol/cm2 that are viable for urban contamination level monitoring.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2021

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

JETSCAPE framework: p+p results

The JETSCAPE framework is a modular and versatile Monte Carlo software package for the simulation of high energy nuclear collisions. In this work we present a new tune of JETSCAPE, called PP19, and validate it by comparison to jet-based measurements in p + p collisions, including inclusive single jet cross sections, jet shape observables, fragmentation functions, charged hadron cross sections, and dijet mass cross sections. These observables in p+p collisions provide the baseline for their counterparts in nuclear collisions. Quantifying the level of agreement of JETSCAPE results with p+p data is thus necessary for meaningful applications of JETSCAPE to A+A collisions. The calculations use the JETSCAPE PP19 tune, defined in this paper, based on version 1.0 of the JETSCAPE framework. For the observables discussed in this work calculations using JETSCAPE PP19 agree with data over a wide range of collision energies at a level comparable to standard Monte Carlo codes. These results demonstrate the physics capabilities of the JETSCAPE framework and provide benchmarks for JETSCAPE users.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Addition of tabulated equation of state and neutrino leakage support to illinoisgrmhd

Here we have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to illinoisgrmhd, an open-source GRMHD code for dynamical spacetimes in the einstein toolkit. These new features are provided by two new, nrpy+-based codes: nrpyeos, which performs highly efficient EOS table lookups and interpolations, and nrpyleakage, which implements a new, adaptive mesh refinement (AMR)-capable neutrino leakage scheme in the einstein toolkit. We have performed a series of strenuous validation tests that demonstrate the robustness of these new codes, particularly on the Cartesian AMR grids provided by carpet. Furthermore, we show results from fully dynamical GRMHD simulations of single unmagnetized neutron stars, and magnetized binary neutron star mergers. This new version of illinoisgrmhd, as well as nrpyeos and nrpyleakage, is pedagogically documented in jupyter notebooks and fully open source. The codes will be proposed for inclusion in an upcoming version of the einstein toolkit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Mitiq: A software package for error mitigation on noisy quantum computers

We introduce Mitiq, a Python package for error mitigation on noisy quantum computers. Error mitigation techniques can reduce the impact of noise on near-term quantum computers with minimal overhead in quantum resources by relying on a mixture of quantum sampling and classical post-processing techniques. Mitiq is an extensible toolkit of different error mitigation methods, including zero-noise extrapolation, probabilistic error cancellation, and Clifford data regression. The library is designed to be compatible with generic backends and interfaces with different quantum software frameworks. We describe Mitiq using code snippets to demonstrate usage and discuss features and contribution guidelines. We present several examples demonstrating error mitigation on IBM and Rigetti superconducting quantum processors as well as on noisy simulators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A heavy-ion production channel of 149 Tb via 63 Cu bombardment of 89 Y

We report the radionuclide 149 Tb (t 1/2 = 4.1 h) is a potential theranostic isotope which can simultaneously be used for targeted-alpha-particle therapy and positron-emission tomography. Feasibility experiments were performed to test a near-symmetric heavy-ion reaction of 63 Cu bombardment on monoisotopic 89 Y. The indirect reaction was studied to avoid isomer production. Offline gamma spectroscopy was used to quantify thick-target physical yields and experimental results show modest agreement to the fusion-evaporation code PACE4. A near-symmetric fission yield was also observed.

149Tb↗
DOE OSTI2021

SOLPS-ITER validation with TCV L-mode discharges

This work presents a quantitative test of SOLPS-ITER simulations against tokamak a configuration variable (TCV) L-mode experiments. These simulations account for drifts, currents, kinetic neutrals, and carbon impurities providing the most complete edge transport simulations for TCV to date. The comparison is performed on nominally identical discharges carried out to assess the effectiveness of TCV’s divertor baffles in the framework of the European Plasma Exhaust program and employs numerous edge diagnostics providing a detailed code-experiment benchmark for TCV. The simulations show a qualitative consistency, but the quantitative differences remain, which are assessed herein. It is found that, for a given separatrix density, the simulations most notably yield a colder, and denser, divertor state with a higher divertor neutral pressure than measured.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Constraining the Halo Mass of Damped Lyα Absorption Systems (DLAs) at z = 2–3.5 Using the Quasar-CMB Lensing Cross-correlation

We study the cross-correlation of damped Lyα systems (DLAs) and their background quasars, using the most updated DLA catalog and the Planck 2018 CMB lensing convergence field. Our measurement suggests that the DLA bias b DLA is smaller than 3.1, corresponding to log (M/M ⊙ h -1 ) ⩽ 12.3 at a confidence of 90%. These constraints are broadly consistent with Alonso et al. and previous measurements by cross-correlation between DLAs and the Lyα forest (e.g., Font-Ribera et al.; Prez-Rfols et al.). Further, our results demonstrate the potential of obtaining a more precise measurement of the halo mass of the high-redshift sources using next generation CMB experiments with a higher angular resolution. The python-based codes and data products of our analysis are available at https://github.com/LittleLin1999/CMB-lensingxDLA.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Simulation budgeting for hybrid effective field theories

In this work, we forecast the number of, and requirements on, N-body simulations needed to train hybrid effective field theory (HEFT) emulators for a range of use cases, using a hybrid of HMcode and perturbation theory as a surrogate model. Our accuracy goals, determined with careful consideration of statistical and systematic uncertainties, are 1% accurate in the high-likelihood range of cosmological parameters, and 2% accurate over a broader parameter space volume for k < 1 h Mpc -1 and z < 3. Focusing in part on the 8-parameter w 0 w a CDM+ m ν cosmological model, we find that < 225 simulations are required to meet our error goals over our wide parameter space, including models with rapidly evolving dark energy, given our simulation and emulator recommendations. For a more restricted parameter space volume, as few as 80 simulations are sufficient. We additionally present simulation forecasts for example use cases, and make the code used in our analyses publicly available. These results offer practical guidance for efficient emulator design and simulation budgeting in future cosmological analyses.

cosmological parameters from LSS↗
DOE OSTI2020

Single-Crystal X-Ray Spectropolarimeter

X-ray spectropolarimetry is a diagnostics technique suitable for measuring intense electron beams and strong magnetic fields in hot, dense plasmas. To increase the accuracy achievable, we investigated the possibility of measuring both orthogonal components of linear polarization simultaneously using a single crystal. We demonstrated experimentally for the first time that crystals with cubic symmetry can be used as polarizing beam splitters. We also found that the polarization-splitting planes in cubic crystals can be accessed in two crystal orientations, a property that can be used to decrease the clutter in polarized spectra. In parallel, we developed a ray-tracing code and are working to add the crystal reflectivity to instrument modeling codes, to help us interpret the spectra and to design spectroscopic and spectropolarimetric instruments. This work is motivated by the need to perform measurements of important quantities in difficult experimental conditions, such as the axial magnetic field at stagnation in magnetized liner inertial fusion loads.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

How numerical simulations helped to achieve breakeven on the NIF

The inertial confinement fusion program relies upon detailed simulations with inertial confinement fusion (ICF) codes to design targets and to interpret the experimental results. These simulations treat as much physics from essential principles as is practical, including laser deposition, cross beam energy transfer, x-ray production and transport, nonlocal thermal equilibrium kinetics, thermal transport, hydrodynamic instabilities, thermonuclear burn, and transport of reaction products. Improvements in radiation hydrodynamic code capabilities and vast increases in computing power have enabled more realistic, accurate 3D simulations that treat all known asymmetry sources. We describe how numerical simulations helped to guide the program, assess the impediments to breakeven, and optimize every aspect of target design. A preshot simulation of the first National Ignition Facility experiment that surpassed breakeven predicted an increased yield that matches the experimental result, within the preshot predicted uncertainty, with a target gain of 1.5. We will cover the key developments in Lawrence Livermore National Laboratory ICF codes that enabled these simulations and give specific examples of how they helped to guide the program.

Marinak, M. M. (ORCID:0009000331127696)↗
DOE OSTI2024

Mixed Stochastic-Deterministic Approach for Many-Body Perturbation Theory Calculations

Here, we present an approach for GW calculations of quasiparticle energies with quasi-quadratic scaling by approximating high-energy contributions to the Green's function in its Lehmann representation with effective stochastic vectors. The method is easy to implement without altering the GW code, converges rapidly with stochastic parameters, and treats systems of various dimensionality and screening response. Our calculations on a 5.75° twisted MoS 2 bilayer show how large-scale GW methods include geometry relaxations and electronic correlations on an equal basis in structurally nontrivial materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

IRIS-BLAS: Towards a Performance Portable and Heterogeneous BLAS Library

This paper presents IRIS-BLAS, a novel heterogeneous and performance portable BLAS library. IRIS-BLAS is built on top of the IRIS runtime and multiple vendor and open-source BLAS libraries. It can transparently use all the architectures/devices available in a heterogeneous system, using the appropriate BLAS library based on the task mapping at run time. Thus, IRIS-BLAS is portable across a broad spectrum of architectures and BLAS libraries, alleviating the worry of application developers about modifying the application source code. Even though the emphasis is on portability, IRIS-BLAS provides competitive or even better performance than other state-of-the-art references. Moreover, IRIS-BLAS offers new features such as efficiently using extremely heterogeneous systems composed of multiple GPUs from different hardware vendors.

Miniskar, Narasinga Rao↗
DOE OSTI2024

Design and validation of a fission-product source for radioactive beams

A neutron-generator-based source of fission products (nuCARIBU) is being developed to replace the existing source at Argonne National Laboratory. Neutrons generated by reaction are thermalized in a compact polyethylene [(moderator and impinge on a thin 235 U foil to produce fission products. Experiments have been conducted to verify the nuCARIBU concept and simulation results using a 4 MeV, 0.16 µ⁢A proton beam from the ATLAS accelerator irradiating a LiF neutron production target and two different simplified versions of the nuCARIBU (C 2 ⁢H 4 ) n moderator with thicknesses of 3.81 cm and 1.91 cm separating it from the enriched 235 U fission foil. Here, the emitted γ-ray peaks of 1383.93 keV ( 92 Sr), 358.0 keV ( 104 Tc), 565.992 keV ( 134 Te), and 258.411 keV ( 138 Xe) obtained from fission products are used to determine the neutron-induced fission rates in the fission foil. The measured fission rates are in good agreement with those simulated using the PHITS Monte Carlo code with user-defined neutron production cross sections based on well-known cross section data. These experimental results, consistent with the simulations, show that the current 252 Cf source (CARIBU) can be replaced by the neutron induced fission source (nuCARIBU).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

A combined experimental and modelling study of granite hydrothermal alteration

Geochemical reactions can induce significant changes of rock reservoir porosity and permeability via mineral dissolution and precipitation processes, affecting the long-term fluid behaviour within various geological systems. Here, the understanding and quantification of these reactions rely on field and experimental studies and on the predictions of reactive transport models. The present study was aimed at assessing the extent to which current geochemical models integrating available mineral dissolution/precipitation rate equations can reproduce the experimental data obtained from 4 to 17-day long hydrothermal alteration experiments of a muscovite-biotite granite and, thus, help provide an accurate description of the evolution of geothermal systems within granitic reservoirs. The experiments were conducted at a constant temperature of 180 °C and over an aqueous fluid pH range of 2 to 8.5, using both mixed-flow and static batch reactors. Modelled major element (K, Al, Si, Ca, and Mg) concentrations were generally in satisfactory agreement with the corresponding measured elemental fluxes – the differences between modelled and experimental values were generally within a factor of 5 – and the predicted identity and mass of formed secondary phases were consistent with the microscopic observations of the reacted solids. However, larger differences between measured and modelled element concentrations were observed when significant amounts of secondary phases formed, notably at pH 2 to 3, and for longer-term batch experiments. Much of this concentration difference stems from the underestimation of the amounts of Al-phases formed at acid to near-neutral pH. Although an idealized rock composition was considered, the observed mismatch between model calculations and experimental data can be attributed to inadequate mineral precipitation reaction rates and a poor description of reactive surface areas in existing geochemical modelling codes. More accurate quantification of precipitation kinetics, including nucleation and growth, and improved descriptions of the temporal change of mineral surface area would enhance the predictive capabilities of reactive transport models and benefit, particularly, the efforts aimed at increasing the sustainability of EGS reservoirs.

58 GEOSCIENCES↗
DOE OSTI2020

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

Impurity transport in ion- and electron-root confinement scenarios at Wendelstein 7-X

This paper reports on the observation of enhanced Ar impurity confinement in high density, turbulence reduced ion-root confinement (IRC) scenarios in Wendelstein 7-X (W7-X). Compared to the central electron-root confinement (CERC), Ar densities are substantially increased up to one order of magnitude and Ar flux profiles switch from slightly positive to pronounced negative fluxes all along the plasma radius in IRC. Estimations of the diffusive D and convective V transport parameter profiles using the STRAHL impurity transport code suggest a fundamentally different behavior in CERC and IRC with, compared to neoclassical values, high diffusive D and typical V profiles in CERC versus strongly reduced D and pronounced negative V profiles in IRC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗