Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LLNL”

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 577 records · Page 32

Ceramic High Efficiency Particulate Air (HEPA) Filter Research and Development at Lawrence Livermore National Laboratory

Ceramic high efficiency particulate air (HEPA) filter development program began at Lawrence Livermore National Laboratory (LLNL) over 20 years ago. Historical incidents motivated current safety systems at nuclear facilities. A ceramic filter that can survive a fire can provide robust, passive safety protection. Research has been conducted to improve the safety of nuclear facilities and to reduce operational and lifecycle costs, through ceramic filter technology that can survive fire conditions. Here, the research focused on applications in both new facilities as well as meeting operational requirements necessary to retrofit filters into existing DOE facilities. The research has developed multiple filter technologies spanning traditional HEPA filter materials to advanced manufacturing technologies (e.g., electrospinning, additive manufacturing, etc.). This communication will present highlights of selected development efforts for ceramic HEPA filter research, current state-of-the-art for ceramic filters, and future needs including technical, regulatory, and commercial efforts.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Leveraging Radiofrequency Identification Success Beyond Hazardous Material Inventory Management at a National Laboratory

Effective inventory management can be overshadowed by conflicting priorities in organizational procedures, particularly in research-focused institutions such as national laboratories that handle expensive, delicate, and hazardous materials. Here, this study investigated the potential of radiofrequency identification (RFID) technology, currently used for hazardous chemical inventory, in applications with higher metal interference and absorption, specifically pressure release device (PRD) compliance and nuclear container management, at Lawrence Livermore National Laboratory (LLNL). This study was done to document best practices to enhance inventory identification speeds for inventory reconciliation and inventory recall and to explore optimal configurations for RFID implementation compared to traditional manual methods of equipment management. Tests were conducted to determine the ideal RFID tag orientation (read at angles of 0°, 90°, and 270°), various container layouts (linear, separated, curved, operational), and ID methods such as manual, barcode, and RFID performing three trials per method per orientation. Results indicated that 0° was the optimal read angle for minimizing metallic interference, and the operational and curved arrangements significantly outperformed the linear and separated configurations in read speed. 3D printed mounts were developed and tested, increasing the read range of the RFID reader by up to 235% in cases of high metallic interference. The RFID technology demonstrated an average speed increase of 65% over a simplified manual identification, which supports the conclusion that RFID is a more efficient method for large hazardous inventory management and equipment reconciliation. Additionally, capturing meta-data, such as location and date, can be used to query for inventory recall and automated updating of record information.

42 ENGINEERING↗

Nucleation and phase transformation pathways in electrolyte solutions investigated by in situ microscopy techniques

Identification of crystal nucleation and growth pathways is of fundamental importance for synthesis of functional materials, which requires control over size, orientation, polymorph, and hierarchical structure, often in the presence of additives used to tune the energy landscape defining these pathways. Furthermore we summarize the recent progress in application of in situ TEM and AFM techniques to monitor or even tune the pathway of crystal nucleation and growth.

36 MATERIALS SCIENCE↗

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗

Energy Dependence of Fission Product Yields from {sup 235}U, {sup 238}U and {sup 239}Pu for Incident Neutron Energies Between 0.5 and 14.8 MeV

Fission Product Yields (FPY) have historically been one of the most observable features of the fission process. They are known to have strong variations that are dependent on the fissioning species, the excitation energy, and the angular momentum of the compound system. However, consistent and systematic studies of the variation of these FPY with energy have proved challenging. This is caused primarily by the nature of the experiments that have traditionally relied on radiochemical procedures to isolate specific fission products. Although radiochemical procedures exist that can isolate all products, each element presents specific challenges and introduces varying degrees of systematic errors that can make inter-comparison of FPY uncertain. Although of high importance in fields such as nuclear forensics and Stockpile Stewardship, accurate information about the energy dependence of neutron induced FPY are sparse, due primarily to the lack of suitable monoenergetic neutron sources. There is a clear need for improved data, and to address this issue, a collaboration was formed between Los Alamos National Laboratory (LANL), Lawrence Livermore National Laboratory (LLNL) and the Triangle Universities Nuclear Laboratory (TUNL) to measure the energy dependence of FPY for {sup 235}U, {sup 238}U and {sup 239}Pu. The measurements have been performed at TUNL, using a 10 MV Tandem Van de Graaff accelerator to produce monoenergetic neutrons at energies between 0.6 MeV to 14.8 MeV through a variety of reactions. The measurements have utilized a dual-fission chamber, with thin (10-100 μg/cm2) reference foils of similar material to a thick (100-400 mg) activation target held in the center between the chambers. This method allows for the accurate determination of the number of fissions that occurred in the thick target without requiring knowledge of the fission cross section or neutron fluence on target. Following activation, the thick target was removed from the dual-fission chamber and gamma-ray counted using shielded HPGe detectors for a period of 1-2 months to determine the yield of various fission products. To the extent possible all irradiation and counting procedures were kept the same to minimize sources of systematic errors. FPY have been determined at incident neutron energies of 0.6, 1.4, 2.4, 3.5, 4.6, 5.5, 8.9 and 14.8 MeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Exploring utilization options of heterogeneous architectures for multi-physics simulations

Top supercomputers in the TOP500 list have transitioned from homogeneous node architectures toward heterogeneous manycore nodes with accelerators and CPUs. These new architectures present significant challenges to developers of large-scale multiphysics applications, especially at DOE laboratories that have invested heavily in scalable MPI codes over decades. Much of these scientific application porting efforts for the new heterogeneous architectures are focused on running the computation on the accelerators, which usually comprise > 90% of the FLOPS of the system, typically using one core per GPU to offload work. Furthermore, we describe an approach to utilizing the remaining FLOPS on a heterogeneous machine by running a portion of the computation on the CPUs cooperatively with the GPU computation. Orthogonally, we evaluate using the additional cores to also offload work to GPUs, effectively oversubscribing the GPUs. We present a proof-of-concept implementation in ARES, a multiphysics ALE-AMR code at LLNL. ARES uses a portability layer, RAJA, which enables us to utilize the same source code for both the CPU and the GPU. We develop an approach to utilize both types of processors cooperatively in a heterogeneous system. Our implementation divides the work between the computing resources via domain decomposition, and utilizes all cores of the CPU and all of the GPUs on the node for computation. Load balancing is necessary to use the heterogeneous resources effectively. We present preliminary results on early delivery pre-Sierra machines at LLNL, showing up to 40% performance benefit of oversubscribing the GPUs (with caveats), and up to an 18% performance benefit of using the CPUs on the heterogeneous nodes for computing in addition to using the GPUs.

97 MATHEMATICS AND COMPUTING↗

Rarefaction-driven Rayleigh–Taylor instability. Part 2. Experiments and simulations in the nonlinear regime

Experiments and large eddy simulation (LES) were performed to study the development of the Rayleigh–Taylor instability into the saturated, nonlinear regime, produced between two gases accelerated by a rarefaction wave. Single-mode two-dimensional, and single-mode three-dimensional initial perturbations were introduced on the diffuse interface between the two gases prior to acceleration. The rarefaction wave imparts a non-constant acceleration, and a time decreasing Atwood number, $A=(\unicode[STIX]{x1D70C}_{2}-\unicode[STIX]{x1D70C}_{1})/(\unicode[STIX]{x1D70C}_{2}+\unicode[STIX]{x1D70C}_{1})$ , where $\unicode[STIX]{x1D70C}_{2}$ and $\unicode[STIX]{x1D70C}_{1}$ are the densities of the heavy and light gas, respectively. Experiments and simulations are presented for initial Atwood numbers of $A=0.49$ , $A=0.63$ , $A=0.82$ and $A=0.94$ . Nominally two-dimensional (2-D) experiments (initiated with nearly 2-D perturbations) and 2-D simulations are observed to approach an intermediate-time velocity plateau that is in disagreement with the late-time velocity obtained from the incompressible model of Goncharov ( Phys. Rev. Lett. , vol. 88, 2002, 134502). Reacceleration from an intermediate velocity is observed for 2-D bubbles in large wavenumber, $k=2\unicode[STIX]{x03C0}/\unicode[STIX]{x1D706}=0.247~\text{mm}^{-1}$ , experiments and simulations, where $\unicode[STIX]{x1D706}$ is the wavelength of the initial perturbation. At moderate Atwood numbers, the bubble and spike velocities approach larger values than those predicted by Goncharov’s model. These late-time velocity trends are predicted well by numerical simulations using the LLNL Miranda code, and by the 2009 model of Mikaelian ( Phys. Fluids. , vol. 21, 2009, 024103) that extends Layzer type models to variable acceleration and density. Large Atwood number experiments show a delayed roll up, and exhibit a free-fall like behaviour. Finally, experiments initiated with three-dimensional perturbations tend to agree better with models and a simulation using the LLNL Ares code initiated with an axisymmetric rather than Cartesian symmetry.

42 ENGINEERING↗

Subtropical Eastern North Pacific SST bias in earth system models

Karthik Balaguru (PNNL), Luke Van Roekel (LLNL), L. Ruby Leung (PNNL) and Milena Veneziani (LLNL) Abstract: This study systematically evaluates the warm SST bias in the Subtropical Eastern North Pacific (SENP), a problem plaguing most CMIP6 models, using the Energy Exascale Earth System Model version 1 (E3SM). Compared to the model at its standard resolution (1o atmosphere, 30-60 km ocean), the SST bias improves in several regions at high resolution (0.25o atmosphere, 18-6 km ocean) except in subtropical eastern boundary regions, such as the SENP. A mixed layer heat budget analysis reveals that errors in the surface radiative fluxes and horizontal heat advection are primarily responsible for the bias. Analysis of HighResMIP models indicates that the lack of improvement in the SENP SST bias with resolution is not native to E3SM alone. A larger CMIP6 inter-model correlation of SST bias with biases in surface winds than radiation highlights the important role of large-scale circulation biases in the SENP SST bias in coupled models.

Balaguru, Karthik↗

Understanding the impact of nuclear-data covariances on various integral responses using adjustment

The EUCLID (Experiments Underpinned by Computational Learning for Improvements in Nuclear Data) project created a library of sensitivities for nine different integral responses with respect to nuclear data. These integral responses were obtained from measurements at LLNL (Lawrence Livermore National Laboratory) pulsed spheres, critical and sub-critical assemblies. At the same time, covariances for ENDF/B-VIII.0 were processed at LANL (Los Alamos National Laboratory). The combination of these data allow us to study the impact of nuclear-data covariances on various integral responses, either by forward-propagating covariances via sensitivities, or by using nuclear data, integral responses, and sensitivities for adjustment. Here, we will present: the impact of 1 H, 9 Be, 12 C, 27 Al, 56 Fe, 235,238 U, and 239,240 Pu ENDF/B-VIII.0 covariances on simulated bounds of the following integral responses: LLNL pulsed-spheres neutron-leakage spectra, the effective neutron multiplication factor, reaction rates, and reactivity coefficients of ICSBEP critical assemblies. Also, adjustment results with the same nuclear-data covariances and responses will be discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Recent results and future plans for a 45 actuator adaptive x-ray optics experiment at the advanced light source

We report on the current status of the Adaptive X-ray Optics project run by Lawrence Livermore National Laboratory (LLNL). LLNL is collaborating with the Advanced Light Source (ALS) to demonstrate a near real-time adaptive X-ray optic. To this end, a custom-built 45 cm long deformable mirror has been installed at ALS beamline 5.3.1 (end station 2) for a two-year period that started in September 2014. We will outline general aspects of the instrument, present results from a recent experimental campaign and touch on future plans for the project.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗

The long road to ignition: An eyewitness account

This paper reviews the many twists and turns in the long journey that culminated in ignition in late 2022 using the laser heated indirect-drive approach to imploding DT filled targets at the National Ignition Facility (NIF), located at the Lawrence Livermore National Laboratory (LLNL). We describe the early origins of the Laser Program at LLNL and key developments such as the paradigm shifting birth of high energy density physics (HEDP) studies with lasers, changes in choice of laser wavelength, and the development of key diagnostics and computer codes. Fulfilling the requirements of the multi-faceted Nova Technical Contract was a necessary condition for the approval of the NIF, but more importantly, the end of the Cold War and the cessation of nuclear testing were key catalysts in that approval, along with the ready-and-waiting field of HEDP. The inherent flexibility of the field of laser driven inertial confinement fusion played a fundamental role in achieving success at the NIF. We describe how the ultimately successful ignition target design evolved from the original “point design” target, through the lessons of experiment. All key aspects of that original design changed: The capsule's materials and size were changed; the hohlraum's materials, size, laser entrance hole size, and gas fills were also all changed, as were the laser pulse shapes that go along with all those changes. The philosophy to globally optimize performance for stability (by raising the adiabat and thus lowering the implosion convergence) was also key, as was progress in target fabrication, and in increasing NIF's energy output. The persistence of the research staff and the steadfast backing of our supporters were also necessary elements in this success. We gratefully acknowledge seven decades of researcher endeavors and four decades of the dedicated efforts of many hundreds of personnel across the globe who have participated in NIF construction, operation, target fabrication, diagnostic, and theoretical advances that have culminated in ignition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fully coupled simulation of cosmic reionization. I. numerical methods and tests

Here, we describe an extension of the Enzo code to enable fully coupled radiation hydrodynamical simulation of inhomogeneous reionization in large similar to(100 Mpc)(3) cosmological volumes with thousands to millions of point sources. We solve all dynamical, radiative transfer, thermal, and ionization processes self-consistently on the same mesh, as opposed to a postprocessing approach which coarse-grains the radiative transfer. But, we employ a simple subgrid model for star formation which we calibrate to observations. The numerical method presented is a modification of an earlier method presented in Reynolds et al. differing principally in the operator splitting algorithm we use to advance the system of equations. Radiation transport is done in the gray flux-limited diffusion (FLD) approximation, which is solved by implicit time integration split off from the gas energy and ionization equations, which are solved separately. This results in a faster and more robust scheme for cosmological applications compared to the earlier method. The FLD equation is solved using the hypre optimally scalable geometric multigrid solver from LLNL. By treating the ionizing radiation as a grid field as opposed to rays, our method is scalable with respect to the number of ionizing sources, limited only by the parallel scaling properties of the radiation solver. We test the speed and accuracy of our approach on a number of standard verification and validation tests. We show by direct comparison with Enzo's adaptive ray tracing method Moray that the well-known inability of FLD to cast a shadow behind opaque clouds has a minor effect on the evolution of ionized volume and mass fractions in a reionization simulation validation test. Finally, we illustrate an application of our method to the problem of inhomogeneous reionization in a 80 Mpc comoving box resolved with 3200(3) Eulerian grid cells and dark matter particles.

79 ASTRONOMY AND ASTROPHYSICS↗

Partial defect verification of spent fuel assemblies by PDET: Principle and field testing in Interim Spent fuel Storage Facility (CLAB) in Sweden

Here, the need for the development of a credible method and instrument for partial defect verification of spent fuel has been emphasized over a few decades in the safeguards communities as the diverted spent fuel pins can be the source of nuclear terrorism or devices. The need is increasingly more important and even urgent as many countries have started to transfer spent fuel to so called "difficult-to-access" areas such as dry storage casks, reprocessing or geological repositories. Partial defect verification is required by IAEA before spent fuel is placed into "difficult-to-access" areas. Earlier, Lawrence Livermore National Laboratory (LLNL) has reported the successful development of a new, credible partial defect verification method for pressurized water reactor (PWR) spent fuel assemblies without use of operator data, and further reported the validation experiments using commercial spent fuel assemblies with some missing fuel pins. The method was found to be robust as the method is relatively invariant to the characteristic variations of spent fuel assemblies such as initial fuel enrichment, cooling time, and burn-up. Since then, the PDET system has been designed and prototyped for 17×17 PWR spent fuel assemblies, complete with data acquisition software and acquisition electronics. In this paper, a summary description of the PDET development followed by results of the first successful field testing using the integrated PDET system and actual spent fuel assemblies performed in a commercial spent fuel storage site, known as Central Interim Spent fuel Storage Facility (CLAB) in Sweden will be presented. In addition to partial defect detection initial studies have determined that the tool can be used to verify the operator declared average burnup of the assembly as well as intra-assembly bunrup levels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Experimental Investigation of Subcritical Neutron and Gamma Noise Methods at the Seven Percent Critical Experiment (7uPCX)

As part of a collaborative international effort organized by the Lawrence Livermore National Laboratory (LLNL), with key participants from the Institut de radioprotection et de sûreté nucléaire (IRSN), Los Alamos National Laboratory (LANL), and the Sandia National Laboratories (SNL), a series of high-multiplication subcritical neutron and gamma noise measurements was planned and executed. The primary aim of this article was to advance detector technology, assess the validity of gamma noise for subcriticality measurements, and nuclear criticality safety, focusing on collecting list-mode or time-series data from various reactor configurations with multiplication values ranging from 20 to 310. This comprehensive dataset enabled a detailed comparative analysis of multiple detector systems and the results of both neutron and gamma noise measurements. In this work, we focus on experimentally comparing the results from neutron and gamma noise measurements. We note good agreement between estimations of the prompt neutron decay constant and demonstrate the effects of changing reactor geometry on the efficiency of the differing methods. Our results agree well with independent experimental measurements and simulations performed by LLNL.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

First Experiments and Radiographs on the MegaJOuLe Neutron Imaging Radiography (MJOLNIR) Dense Plasma Focus

We report a dense plasma focus (DPF) is a relatively compact coaxial plasma gun, which completes its discharge as a Z-pinch. These devices are designed to operate at a variety of scales to produce short (<100 ns) pulses of ions, X-rays, and/or neutrons. LLNL recently constructed and brought into operation a new device, the MegaJOuLe Neutron Imaging Radiography (MJOLNIR) DPF, which is designed for radiography and high-yield operations. This device has been commissioned and has achieved neutron yields of up to 3.8E11 neutrons/pulse at 2.5-MA peak current while operating at up to 1 MJ of stored energy in its original pulsed power configuration. MJOLNIR is equipped with a wide range of diagnostics, including nuclear activation detectors, neutron time-of-flight (nToF) detectors, a fast-framing camera, optical light gates, and a time-gated neutron and X-ray imager. LLNL also runs unique particle-in-cell (PIC) simulations of DPF discharges in the Chicago code and has gained significant insight into the various physical factors that influence neutron yield. MJOLNIR is one of the first DPFs whose design and continual upgrades are heavily influenced by model predictions. In this article, we describe insights from modeling, device operation, and recent results. Comparisons between modeling predictions and measurements, as well as X-ray and neutron images are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Sierra Center of Excellence: Lessons Learned

The introduction of heterogeneous computing via GPUs from the Sierra architecture represented a significant shift in direction for computational science at Lawrence Livermore National Laboratory (LLNL), and therefore required significant preparation. Over the last five years, the Sierra Center of Excellence (CoE) has brought employees with specific expertise from IBM and NVIDIA together with LLNL in a concentrated effort to prepare applications, system software, and tools for the Sierra supercomputer. This article shares the process we applied for the CoE and documents lessons learned during the collaboration, with the hope that others will be able to learn from both our success and intermediate setbacks. We describe what we have found to work for the management of such a collaboration and best practices for algorithms and source code, system configuration and software stack, tools, and application performance.

97 MATHEMATICS AND COMPUTING↗

Concepts and Meta-Programming Library

CAMP provides a portable optimized implementation of template meta-programming primitives. All capabilities found within are available in other open source projects, including metal and mpl among others, but no other library we know of supports all compilers required by LLNL codes. CAMP serves as an efficient facility for writing portable compile-time operations in C++ codes. It was developed as part of the RAJA project to facilitate more complex compile-time introspection and better error message delivery to users, but has since proven useful in Umpire, CHAI and other pre-release LLNL software projects and is thus being re-factored to stand alone. The CAMP library uses portable C++ template mechanisms almost exclusively, with the exception of intrinsics provided for performance by some compilers such as clang and gee. There are no private or protected interfaces and no proprietary information included in the source. Any and all future versions will limit themselves to implementing commonly available and well known primitives for compile-time manipulation of types. A new public github repository will be created to host the code independent of its current location in the raja repository, likely at github.com/llnl/camp.

Kunen, AdamJ.↗