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 667 records · Page 37

Optimizing hybrid rotational femtosecond/picosecond coherent anti-Stokes Raman spectroscopy in nitrogen at high pressures and temperatures

We demonstrate the use of hybrid rotational femtosecond/picosecond (fs/ps) coherent anti-Stokes Raman spectroscopy (HR-CARS) as a technique for temperature measurements in nitrogen gas at high pressures and temperatures. A broadband pulse shaper-adjusted 42 fs pulse interacts with a narrow-bandwidth, frequency-upconverted 5.5 ps pulse in a cell containing N 2 at pressures of 1–70 atm and temperatures of 300–1000 K. A computational code is used to model spectra and fit experimental results to obtain best-fit temperatures. We demonstrate good qualitative fits as well as good accuracy and precision between thermocouple measured and best-fit temperatures over the explored pressure and temperature regimes. The overall average percentage temperature difference between thermocouple measurements and best-fit temperatures is − <#comment/> 0 . 3 % <#comment/> with a standard deviation of 7.1%, showing the suitability of HR-CARS for characterizing high-pressure and -temperature environments.

Torge Mecker, Nils (ORCID:0000000238977170)↗

Sierra/SolidMechanics 4.56 Theory Manual

Presented in this document are the theoretical aspects of capabilities contained in the Sierra/SM code. This manuscript serves as an ideal starting point for understanding the theoretical foundations of the code. For a comprehensive study of these capabilities, the reader is encouraged to explore the many references to scientific articles and textbooks contained in this manual. It is important to point out that some capabilities are still in development and may not be presented in this document. Further updates to this manuscript will be made as these capabilities come closer to production level.

97 MATHEMATICS AND COMPUTING↗

Towards Predictive Plasma Science and Engineering through Revolutionary Multi-Scale Algorithms and Models (Final Report)

This report describes the high-level accomplishments from the Plasma Science and Engineering Grand Challenge LDRD at Sandia National Laboratories. The Laboratory has a need to demonstrate predictive capabilities to model plasma phenomena in order to rapidly accelerate engineering development in several mission areas. The purpose of this Grand Challenge LDRD was to advance the fundamental models, methods, and algorithms along with supporting electrode science foundation to enable a revolutionary shift towards predictive plasma engineering design principles. This project integrated the SNL knowledge base in computer science, plasma physics, materials science, applied mathematics, and relevant application engineering to establish new cross-laboratory collaborations on these topics. As an initial exemplar, this project focused efforts on improving multi-scale modeling capabilities that are utilized to predict the electrical power delivery on large-scale pulsed power accelerators. Specifically, this LDRD was structured into three primary research thrusts that, when integrated, enable complex simulations of these devices: (1) the exploration of multi-scale models describing the desorption of contaminants from pulsed power electrodes, (2) the development of improved algorithms and code technologies to treat the multi-physics phenomena required to predict device performance, and (3) the creation of a rigorous verification and validation infrastructure to evaluate the codes and models across a range of challenge problems. These components were integrated into initial demonstrations of the largest simulations of multi-level vacuum power flow completed to-date, executed on the leading HPC computing machines available in the NNSA complex today. These preliminary studies indicate relevant pulsed power engineering design simulations can now be completed in (of order) several days, a significant improvement over pre-LDRD levels of performance.

42 ENGINEERING↗

TOUGH3-FLAC3D: a modeling approach for parallel computing of fluid flow and geomechanics

The recent development of the TOUGH3 code allows for a faster and more reliable fluid flow simulator. At the same time, new versions of FLAC3D are released periodically, allowing for new features and faster execution. In this paper, we present the first implementation of the coupling between TOUGH3 and FLAC3Dv6/7, maintaining parallel computing capabilities for the coupled fluid flow and geomechanical codes. We compare the newly developed version with analytical solutions and with the previous approach, and provide some performance analysis on different meshes and varying the number of running processors. Finally, we present two case studies related to fault reactivation during CO 2 sequestration and nuclear waste disposal. The use of parallel computing allows for meshes with a larger number of elements, and hence more detailed understanding of thermo-hydro-mechanical processes occurring at depth.

58 GEOSCIENCES↗

ExaCA: A performance portable exascale cellular automata application for alloy solidification modeling

Modeling the as-solidified grain structures that form during alloy processing is a critical component in understanding process-property relationships, particularly for additive manufacturing (AM) where grain structure is very sensitive to processing conditions. While cellular automata (CA)-based models have proven able to predict aspects of microstructure for several alloys and AM process conditions, long run times and large resource sets required limit the utility and the problem size to which existing CA models can be applied. As part of the ExaAM project, an initiative within the Exascale Computing Project (ECP) to develop, test, and optimize an exascale-capable coupled and self-consistent model of AM parts, we developed ExaCA (https://github.com/LLNL/ExaCA) for the liquid–solid phase transformation in the wake of AM melt pools. The CA-based code is parallelized using MPI and the Kokkos programming model, the latter enabling simulation on both CPUs and GPUs within a single-source implementation. Here, we detail the steps taken to transform a baseline, MPI-based CA code into one that is performant on CPUs and GPUs. Performance testing of ExaCA on Summit (a pre-exascale machine at Oak Ridge National Laboratory) was used to quantify CPU–GPU speedup comparing with equal numbers of nodes. Testing showed comparable CPU performance to the MPI-only CA code and a 5-20x speedup when running AM-based test problems using GPUs. The improved performance of CA through GPU utilization and the performance portable nature of ExaCA will enable accurate part-scale modeling by harnessing the power of current and future generations of high performance computing resources. Future work will include improving the strong scaling of ExaCA on GPUs by reducing load imbalance associated with the locality of the problem, and continuing performance optimization across exascale hardware.

36 MATERIALS SCIENCE↗

Simulating Magnetized Electron Cooling for EIC with JSPEC

We present a possible electron cooling configuration for the proposed Electron Ion Collider (EIC) facility, developed using a Nelder-Mead Simplex optimization procedure built into JSPEC, an electron cooling code developed at Jefferson Lab. We show the time evolution of the emittance of the ion beam in the presence of this cooler evaluated assuming the ion distribution remains Gaussian. We also show that bi-gaussian distributions emerge in simulations of ion macro-particles. We show how intra-beam scattering can be treated with a core-tail model in simulations of ion macro-particles. The Sirepo/JSPEC* and Sirepo/Jupyter** apps will be presented, with instructions enabling the community to reproduce our simulations.

Coleman, S. J.↗

Particle-in-cell Simulations of Mildly Relativistic Outflows in Kilonova Emissions

Abstract The electromagnetic emission from neutron star mergers is comprised of multiple components. Synchrotron emission from the disk-powered jet and thermal emission from the merger ejecta (powered by a variety of sources) are among the most studied sources. The low masses and high velocities of the merger ejecta quickly develop conditions where emission from collisionless shocks becomes critical and synchrotron emission from the merger ejecta constitutes a third component to the observed signal. The aim of this project is to examine shock development, magnetic field generation, and particle acceleration in the case of mildly relativistic shocks, which are expected when the tidal ejecta of neutron star mergers drive a shock into the external medium. Using LANL’s vector particle-in-cell (VPIC) code, we have run simulations of such mildly relativistic, collisionless, weakly magnetized plasmas and computed the resultant magnetic fields and particle energy spectra. We show the effects of varying plasma conditions, as well as explore the validity of using different proton-to-electron mass ratios in VPIC. Our results have implications for observing late-time electromagnetic counterparts to gravitational wave detections of neutron star mergers.

79 ASTRONOMY AND ASTROPHYSICS↗

Shifting institutional culture to develop climate solutions with Open Science

To address our climate emergency, “we must rapidly, radically reshape society”—Johnson & Wilkinson, All We Can Save. In science, reshaping requires formidable technical (cloud, coding, reproducibility) and cultural shifts (mindsets, hybrid collaboration, inclusion). We are a group of cross-government and academic scientists that are exploring better ways of working and not being too entrenched in our bureaucracies to do better science, support colleagues, and change the culture at our organizations. We share much-needed success stories and action for what we can all do to reshape science as part of the Open Science movement and 2023 Year of Open Science.

54 ENVIRONMENTAL SCIENCES↗

Measuring the oscillator strength of intercombination lines of helium-like V ions in a laser-produced-plasma

In this work, we present results of measurements of the oscillator strength of intercombination lines of He-like Vanadium ions in high energy density (HED) laser-produced-plasmas and compare them with the simulations from commonly used codes and data from the NIST database. Whilst not yet sufficiently accurate to constrain different trusted atomic-physics models for the particular system studied, our results are in agreement with the available data within experimental error bars, yet differ from cruder approximations of the oscillator strength used in certain atomic-kinetics packages, suggesting that this general method could be further extended to be used as a measurement of the oscillator strength of additional atomic transitions under the extreme conditions that are achieved in HED experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical sciences by providing tools to support and guide the experimental effort and for prediction of chemical and materials properties. In this regard, a special role has been played by electronic structure packages where complex chemical and materials processes can be modeled using first-principle-driven methodologies. Over the last few decades, the rapid development of computing technologies and a tremendous increase in computational power has offered a unique chance to study complex chemical transformations using sophisticated and predictive many-body techniques to describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, a critical role has been played by novel parallel algorithms capable of taking advantage of computational resources to address polynomial scaling of electronic structure methods. NWChem was among the first electronic structure codes that focused on delivering scalable performance for electronic structure simulations. Herein, we briefly review the NWChem suite of computational codes including its history, design principles, parallel tools, current capabilities, outreach and outlook.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reactor Performance and Safety Characteristics of Beryllium-Based Composite Moderators as Replacements for Graphite in mHTGRs

Here, this study evaluates beryllium-based two-phase composite moderators as an alternative to graphite in an evaluation of reactor performance and safety characteristics. Historically, modular high-temperature gas-cooled reactors (mHTGRs) use graphite as a moderator because of its high moderating ratio and reasonable thermal properties; however, graphite has unfavorable properties under irradiation, which can require component replacement and a significant radioactive waste burden. In this assessment, we explore advanced moderators comprised of magnesium oxide (MgO) as the host matrix and beryllium metal and/or beryllium oxide (Be and/or BeO) as the entrained moderating phase. For the reactor performance and thermal-hydraulic safety analysis, the core design model of the General Atomics mHTGR-350 was used to demonstrate the feasibility of a “drop-in” replacement of graphite using the beryllium-based moderators. We employed the neutronics code Serpent to analyze the moderating behavior of the composite moderators with comparisons drawn to graphite. We performed a scoping analysis of accidents for mHTGRs using RELAP to show that these moderators do not present impediments to safety and are expected to stay within temperature limits. Measured thermophysical properties of the composite moderators are used in the thermal-hydraulic assessments. Our analysis reveals that the two-phase composite MgO-matrix beryllium-based moderators are a suitable replacement for graphite.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Not just for programmers: How GitHub can accelerate collaborative and reproducible research in ecology and evolution

Abstract Researchers in ecology and evolutionary biology are increasingly dependent on computational code to conduct research. Hence, the use of efficient methods to share, reproduce, and collaborate on code as well as document research is fundamental. GitHub is an online, cloud‐based service that can help researchers track, organize, discuss, share, and collaborate on software and other materials related to research production, including data, code for analyses, and protocols. Despite these benefits, the use of GitHub in ecology and evolution is not widespread. To help researchers in ecology and evolution adopt useful features from GitHub to improve their research workflows, we review 12 practical ways to use the platform. We outline features ranging from low to high technical difficulty, including storing code, managing projects, coding collaboratively, conducting peer review, writing a manuscript, and using automated and continuous integration to streamline analyses. Given that members of a research team may have different technical skills and responsibilities, we describe how the optimal use of GitHub features may vary among members of a research collaboration. As more ecologists and evolutionary biologists establish their workflows using GitHub, the field can continue to push the boundaries of collaborative, transparent, and open research.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Hybrid-gyrokinetic simulations of low- n toroidal Alfvén eigenmodes using gKPSP

Here, we report a benchmark study of toroidal Alfvén eigenmode (TAE) simulation using the hybrid-gyrokinetic code GyroKinetic Plasma Simulation Program (gKPSP). A simulation capability for energetic particles based on the gyrokinetic δf method has been newly implemented in the gKPSP code. Benchmark simulations have been performed in both circular and realistic tokamak geometries. Good agreement has been found with previously reported results, demonstrating the new capability of the gKPSP code. We have investigated the effects of the distribution function on TAE stability by examining both isotropic and anisotropic slowing-down distributions of energetic particles. The slowing-down distribution produces a higher linear growth rate than a Maxwellian distribution, while a growth rate scan with the anisotropy shows an opposite trend. This can be attributed to competition between Landau damping and the linear drive, which are correlated with the fraction of resonant passing particles and their distribution in phase space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DRACO: An Overview [Slides]

DRACO (Diffusion ReACtiOn) is a diffusion and chemistry code designed to: 1) Operate on 3D with an unstructured grid defining an arbitrary geometry of interacting parts. 2) Generate its own meshes and use meshes created by other software. 3)Model the transport of any number of diffusing quantities: Concentrations, pressures, temperature, etc. 4) Allow diffusion coefficients to depend in an arbitrary way on concentration, temperature, position, time, etc. 5) Model general chemistry between concentrations with arbitrary reaction rates. 6) Allow arbitrary initial conditions, boundary conditions, and sources/sinks. 7) Allow all of the above to be specified by the user.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Omics-to-Reactive-Transport (ORT): A workflow linking genome-scale metabolic models with reactive transport codes

Motivation: Nutrient and contaminant behavior in the subsurface are governed by multiple coupled hydrobiogeochemical processes which occur across different temporal and spatial scales. Accurate description of macroscopic system behavior requires accounting for the effects of microscopic and especially microbial processes. Microbial processes mediate precipitation and dissolution and change aqueous geochemistry, all of which impacts macroscopic system behavior. As `omics data describing microbial processes is increasingly affordable and available, novel methods for using this data quickly and effectively for improved ecosystem models are needed. Results: We propose a workflow (`Omics to Reactive Transport – ORT) for utilizing metagenomic and environmental data to describe the effect of microbiological processes in macroscopic reactive transport models. This workflow utilizes and couples two open-source software packages: KBase (a software platform for systems biology) and PFLOTRAN (a reactive transport modeling code). We describe the architecture of ORT and demonstrate an implementation using metagenomic and geochemical data from a river system. Our demonstration uses microbiological drivers of nitrification and denitrification to predict nitrogen cycling patterns which agree with those provided with generalized stoichiometries. While our example uses data from a single measurement, our workflow can be applied to spatiotemporal metagenomic datasets to allow for iterative coupling between KBASE and PFLOTRAN. Live, interactive models, which incorporate the results from this narrative into a PFLOTRAN simulation, are available (without login) at https://pflotranmodeling.paf.subsurfaceinsights.com/pflotran-simple-model/.

Rubinstein, Rebecca L↗

First-principles equation of state database for warm dense matter computation

We put together a first-principles equation of state (FPEOS) database for matter at extreme conditions by combining results from path integral Monte Carlo and density functional molecular dynamics simulations of the elements H, He, B, C, N, O, Ne, Na, Mg, Al, and Si as well as the compounds LiF , B 4 C , BN , CH 4 , CH 2 , C 2 H 3 , CH , C 2 H , MgO , and MgSiO 3 . For all these materials, we provide the pressure and internal energy over a density-temperature range from ~ 0.5 to 50 g cm - 3 and from ~ 10 4 to 10 9 K, which are based on ~ 5000 different first-principles simulations. We compute isobars, adiabats, and shock Hugoniot curves in the regime of L - and K -shell ionization. Invoking the linear mixing approximation, we study the properties of mixtures at high density and temperature. Furthermore, we derive the Hugoniot curves for water and alumina as well as for carbon-oxygen, helium-neon, and CH-silicon mixtures. We predict the maximal shock compression ratios of H 2 O , H 2 O 2 , Al 2 O 3 , CO , and CO 2 to be 4.61, 4.64, 4.64, 4.89, and 4.83, respectively. Finally we use the FPEOS database to determine the points of maximum shock compression for all available binary mixtures. We identify mixtures that reach higher shock compression ratios than their end members. We discuss trends common to all mixtures in pressure-temperature and particle-shock velocity spaces. In the Supplemental Material, we provide all FPEOS tables as well as computer codes for interpolation, Hugoniot calculations, and plots of various thermodynamic functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Deformation of the $0^+_{1,2}$ states in 110 Cd from low-energy Coulomb excitation

Electromagnetic properties of 110 Cd were studied via low-energy Coulomb excitation with 32 S and 14 N beams. Magnitudes and relative signs of eight E2 matrix elements, including quadrupole moments of the $2^+_1$ and $2^+_2$ states, were determined using the least-squares code GOSIA. From those, quadrupole deformation parameters of the $0^+_{1,2,3}$ states were inferred, providing for the first time conclusive evidence for the non-axial character of the ground state in 110 Cd. The experimental results were compared with new calculations using the general quadrupole collective Bohr Hamiltonian model with SLy4 and UNEDF0 interactions. The non-axiality of the ground state is reproduced by the present calculations, independently of the interaction used.

Collective models↗

PETSc/TAO Users Manual (Rev. 3.19)

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations and related problems on high-performance computers. PETSc/TAO 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 distributed memory communication. PETSc/TAO includes a large suite of parallel linear solvers, nonlinear solvers, time integrators, and opti mization that may be used in application codes written in Fortran, C, C++, and Python (via petsc4py; see Getting Started). 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 packages such as MATLAB or 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 lldb, 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. For many tasks a package such as MATLAB is often the best tool; PETSc is not intended for the classes of problems for which effective MATLAB code can be written. There are several packages, built on PETSc, that may satisfy your needs without requiring directly using PETSc. We recommend reviewing these packages functionality before starting to code directly with PETSc. PETSc can be used to provide a “MPI parallel linear solver” in an otherwise sequential, or OpenMP parallel code. This approach cannot provide extremely large improvements in the application time by utilizing large numbers of MPI processes but can still improve the performance. Certainly all parts of a previously sequential code need not be parallelized but the matrix generation portion must be parallelized to expect true scalability to large numbers of MPI processes. See PCMPI for details on how to utilize the PETSc MPI linear solver server. Since PETSc is under continued development, small changes in usage and calling sequences of routines will occur. PETSc has been supported for twenty-five years; see mailing list information on our website for information on contacting support.

97 MATHEMATICS AND COMPUTING↗