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 145 records · Page 8

Constraining nucleon effective masses with flow and stopping observables from the S π RIT experiment

Properties of the nuclear equation of state (EoS) can be probed by measuring the dynamical properties of nucleus-nucleus collisions. In this study, we present the directed flow (υ 1 ), elliptic flow (υ 2 ) and stopping (VarXZ) measured in fixed target Sn + Sn collisions at 270 AMeV with the SπRIT Time Projection Chamber. We perform Bayesian analyses in which EoS parameters are varied simultaneously within the Improved Quantum Molecular Dynamics-Skyrme (ImQMD-Sky) transport code to obtain a multivariate correlated constraint. The varied parameters include symmetry energy, S 0 , and slope of the symmetry energy, L, at saturation density, isoscalar effective mass, $m^{⁎}_{S}$/$m_{N}$ , isovector effective mass, $m^{⁎}_{υ}$/$m_{N}$ and the in-medium cross-section enhancement factor η. We find that the flow and VarXZ observables are sensitive to the splitting of proton and neutron effective masses and the in-medium cross-section. Comparisons of ImQMD-Sky predictions to the SπRIT data suggest a narrow range of preferred values for $m^{⁎}_{S}$/$m_{N}$, $m^{⁎}_{υ}$/$m_{N}$ and η.

Astronomy & Astrophysics↗

SNEWPY: A Data Pipeline from Supernova Simulations to Neutrino Signals

Current neutrino detectors will observe hundreds to thousands of neutrinos from a Galactic supernovae, and future detectors will increase this yield by an order of magnitude or more. With such a data set comes the potential for a huge increase in our understanding of the explosions of massive stars, nuclear physics under extreme conditions, and the properties of the neutrino. However, there is currently a large gap between supernova simulations and the corresponding signals in neutrino detectors, which will make any comparison between theory and observation very difficult. SNEWPY is an open-source software package which bridges this gap. The SNEWPY code can interface with supernova simulation data to generate from the model either a time series of neutrino spectral fluences at Earth, or the total time-integrated spectral fluence. Data from several hundred simulations of core-collapse, thermonuclear, and pair-instability supernovae is included in the package. This output may then be used by an event generator such as sntools or an event rate calculator such as SNOwGLoBES. Additional routines in the SNEWPY package automate the processing of the generated data through the SNOwGLoBES software and collate its output into the observable channels of each detector. In this paper we describe the contents of the package, the physics behind SNEWPY, the organization of the code, and provide examples of how to make use of its capabilities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Thermal Hydraulic Experimental Test Article - Fiscal Year 2022 Final Report

The Thermal Hydraulic Experimental Test Article (THETA) is a METL vessel experiment designed for testing and validating sodium fast reactor components and phenomena. THETA has been scaled using a non-dimensional Richardson number approach to represent temperature distributions during nominal and loss of flow conditions in a sodium fast reactor (SFR), this analysis was detailed in the THETA FY19 report. The facility is being constructed with versatility in mind, allowing for the installation of various immersion heaters, heat pipes, and heat exchangers without significant facility modification. THETA was designed in collaboration with systems code experts to inform the geometry and sensor placement to acquire the highest value code validation data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SNEWPY: A Data Pipeline from Supernova Simulations to Neutrino Signals

Current neutrino detectors will observe hundreds to thousands of neutrinos from Galactic supernovae, and future detectors will increase this yield by an order of magnitude or more. With such a data set comes the potential for a huge increase in our understanding of the explosions of massive stars, nuclear physics under extreme conditions, and the properties of the neutrino. However, there is currently a large gap between supernova simulations and the corresponding signals in neutrino detectors, which will make any comparison between theory and observation very difficult. SNEWPY is an open-source software package that bridges this gap. The SNEWPY code can interface with supernova simulation data to generate from the model either a time series of neutrino spectral fluences at Earth, or the total time-integrated spectral fluence. Data from several hundred simulations of core-collapse, thermonuclear, and pair-instability supernovae is included in the package. This output may then be used by an event generator such as sntools or an event rate calculator such as the SuperNova Observatories with General Long Baseline Experiment Simulator (SNOwGLoBES). Additional routines in the SNEWPY package automate the processing of the generated data through the SNOwGLoBES software and collate its output into the observable channels of each detector. In this paper we describe the contents of the package, the physics behind SNEWPY, the organization of the code, and provide examples of how to make use of its capabilities.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of fast-ion losses measured in MAST-U via infrared thermography and a Fast Ion Loss Detector

Fast-ion losses need to be monitored to avoid damage to plasma facing components. In existing experimental devices, the scintillator-based fast-ion loss detector (FILD) is the most advanced diagnostic for measuring fast-ion losses. However, FILDs provide only local information about the losses. Infrared (IR) thermography can be used as a complementary tool for more global monitoring of the deposition of fast-ion losses on the wall, at the expense of no velocity-space resolution. IR cameras measure the temperature of the plasma facing components. This measurement, determined by a combined effect of the thermal plasma, radiation, neutrons and fast-ion losses, can be decomposed to infer the fast-ion load on the tokamak wall. In this manuscript, a workflow to estimate fast-ion losses via IR thermography is applied to the MAST-U spherical tokamak, using a 1D approximation to extract the experimental heat flux on the FILD front face from IR data. To numerically estimate the different contributions to this total heat flux, the field-line tracing environment SMITER is used to calculate the thermal plasma contribution, the orbit-following Monte-Carlo code ASCOT to estimate the fast-ion losses, and bolometry measurements for the radiation. To validate the workflow, two discharges, L-mode plasmas with low MHD activity, were executed using on and off-axis beams, respectively. The experimentally and numerically estimated heat flux are of the same order of magnitude for the on-axis heated scenario, with a strong dependence of the estimated fast-ion losses contribution on the fit to the kinetic profiles used as input. This is also true for the off-axis heated scenario, where the total numerically estimated heat flux is 2.1 or 1.3 times higher than the maximum experimentally estimated heat flux, depending on the ASCOT input used.

FILD↗

BEYONDPLANCK VIII. Efficient sidelobe convolution and corrections through spin harmonics

We introduce a new formulation of the Conviqt convolution algorithm in terms of spin harmonics, and apply this to the problem of sidelobe correction for BEYONDPLANCK, the first end-to-end Bayesian Gibbs sampling framework for CMB analysis. We compare our implementation to the previous Planck LevelS implementation, and find good agreement between the two codes in terms of accuracy, but with a speed-up reaching a factor of 3–10, depending on the frequency bandlimits, l max and m max . The new algorithm is significantly simpler to implement and maintain, since all low-level calculations are handled through an external spherical harmonic transform library. We find that our mean sidelobe estimates for Planck LFI are in good agreement with previous efforts. Additionally, we present novel sidelobe rms maps that quantify the uncertainty in the sidelobe corrections due to variations in the sky model.

79 ASTRONOMY AND ASTROPHYSICS↗

Error-Detectable Bosonic Entangling Gates with a Noisy Ancilla

Bosonic quantum error correction has proven to be a successful approach for extending the coherence of quantum memories, but to execute deep quantum circuits, high-fidelity gates between encoded qubits are needed. To that end, we present a family of error-detectable two-qubit gates for a variety of bosonic encodings. From a new geometric framework based on a “Bloch sphere” of bosonic operators, we construct ZZ L ⁡(θ) and exponential-swap(θ) gates for the binomial, four-legged cat, dual-rail, and several other bosonic codes. The gate Hamiltonian is simple to engineer, requiring only a programmable beam splitter between two bosonic qubits and an ancilla dispersively coupled to one qubit. This Hamiltonian can be realized in circuit QED hardware with ancilla transmons and microwave cavities. The proposed theoretical framework was developed for circuit QED but is generalizable to any platform that can effectively generate this Hamiltonian. Crucially, one can also detect first-order errors in the ancilla and the bosonic qubits during the gates. We show that this allows one to reach error-detected gate fidelities at the 0.01% level with today’s hardware, limited only by second-order hardware errors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

2D radial-azimuthal particle-in-cell benchmark for E × B discharges

In this paper we propose a representative simulation test-case of E × B discharges accounting for plasma wall interactions with the presence of both the electron cyclotron drift instability and the modified-two-stream-instability. Seven independently developed particle-in-cell (PIC) codes have simulated this benchmark case, with the same specified conditions. The characteristics of the different codes and computing times are given. Here, results show that both instabilities were captured in a similar fashion and good agreement between the different PIC codes is reported as main plasma parameters were closely related within a 5% interval. The number of macroparticles per cell was also varied and statistical convergence was reached. Detailed outputs are given in the supplementary data, to be used by other similar groups in the perspective of code verification.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Machine learning enabled discovery of application dependent design principles for two-dimensional materials

Abstract The unique electronic and mechanical properties of two-dimensional (2D) materials make them promising next-generation candidates for a variety of applications. Large-scale searches for high-performing 2D materials are limited to calculating descriptors with computationally demanding first-principles density functional theory. In this work, we alleviate this issue by extending and generalizing crystal graph convolutional neural networks to systems with planar periodicity and train an ensemble of models to predict thermodynamic, mechanical and electronic properties. We carry out a screening of nearly 45,000 structures for two separate applications: mechanical strength and photovoltaics. By collecting statistics of the screened candidates, we investigate structural and compositional design principles that impact the properties of the structures surveyed. Our approach recovers some well-accepted design rules: hybrid organic-inorganic perovskites with lead and tin tend to be good candidates for solar cell applications and titanium based MXenes usually have high stiffness coefficients. Interestingly, other members of the group 4 elements also contribute to increasing the mechanical strength of MXenes. For all-inorganic perovskites, we discover some compositions that have not been deeply studied in the field of photovoltaics and thus open up paths for further investigation. We open-source the code-base to spur further development in this space.

Venturi, Victor (ORCID:0000000156727549)↗

Internal Radiation Dose Evaluation For An Unruptured Post Release Tristructural Istropic Fuel Particle for Advanced and Micro-reactor Applications

There are unique benefits from advanced/micro reactor designs and fuel types that offer safety features in the case of an accident that may reduce environmental consequences compared to conventional reactors and fuels. Several micro reactor designs, for both civilian and military applications, are proposed using tristructural-isotropic (TRISO) fuel. TRISO particles are a robust advanced nuclear fuel type that leads to the unique question of how unruptured, activated TRISO particles will interact with humans. Human exposure to nuclear fuel, or radionuclides, includes all pathways of exposure. However, TRISO particles are 0.9 mm in size, and that particle size restricts internal dose assessment to the ingestion pathway; eliminating the inhalation and dermal absorption pathways. Furthermore, activity of the TRISO particle was established by High Temperature Engineering Test Reactor simulations assuming the TRISO particle was one day post-release from the reactor. International Commission of Radiological Protection publications 100 and 133 were used to develop the internal dosimetry model for an unruptured TRISO particle for the International Commission of Radiological Protection reference male and female dose assessment. The rectosigmoid colon had the highest internal organ dose for both reference male (2.1 Sv) and reference female (2.3 Sv). The internal dose from ingestion of the scenario-defined TRISO particle was 0.25 Sv for the reference male and 0.29 Sv for the refence female, which exceeds the annual occupational effective dose limit of 0.05 Sv in the code of federal regulations, 10.CFR.20 Subpart C.

61 RADIATION PROTECTION AND DOSIMETRY↗

Hot-electron generation at direct-drive ignition-relevant plasma conditions at the National Ignition Facility

Laser–plasma interaction instabilities can be detrimental for direct-drive inertial confinement fusion by generating high-energy electrons that preheat the target. Moreover, an experimental platform has been developed and fielded on the National Ignition Facility to investigate hot-electron production from laser–plasma instabilities at direct-drive ignition-relevant conditions. The radiation-hydrodynamic code DRACO has been used to design planar-target experiments that generate plasma and interaction conditions comparable to direct-drive ignition designs: I L ~10 15 W/cm 2 , T e > 3 keV, and density-gradient scale lengths of L n ~ 600 μ m in the quarter-critical density region. The hot-electron properties were inferred by comparing the experimentally observed hard x-ray spectra to Monte Carlo simulations of hard x-ray emission from hot electrons depositing energy in the target. Hot-electron temperatures of ~40 keV to 60 keV and the fraction of laser energy converted to hot electrons of ~0.5% to 5% were inferred in plastic targets for laser intensities at the quarter-critical density surface of (~4 to 14) × 10 14 W/cm 2 . The use of silicon ablators was found to mitigate the hot-electron preheat by increasing the threshold laser intensity for hot-electron generation from ~3.5 × 10 14 W/cm 2 in plastic to ~6 × 10 14 W/cm 2 in silicon. The overall hot-electron production is further reduced in silicon ablators when the intensity threshold is exceeded.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Deep learning for NLTE spectral opacities

Computer simulations of high energy density science experiments are computationally challenging, consisting of multiple physics calculations including radiation transport, hydrodynamics, atomic physics, nuclear reactions, laser–plasma interactions, and more. To simulate inertial confinement fusion (ICF) experiments at high fidelity, each of these physics calculations should be as detailed as possible. However, this quickly becomes too computationally expensive even for modern supercomputers, and thus many simplifying assumptions are made to reduce the required computational time. Much of the research has focused on acceleration techniques for the various packages in multiphysics codes. In this work, we explore a novel method for accelerating physics packages via machine learning. The non-local thermodynamic equilibrium (NLTE) package is one of the most expensive calculations in the simulations of indirect drive inertial confinement fusion, taking several tens of percent of the total wall clock time. We explore the use of machine learning to accelerate this package, by essentially replacing the physics calculation with a deep neural network that has been trained to emulate the physics code. Overall, we demonstrate the feasibility of this approach on a simple problem and perform a side-by-side comparison of the physics calculation and the neural network inline in an ICF Hohlraum simulation. We show that the neural network achieves a 10× speed up in NLTE computational time while achieving good agreement with the physics code for several quantities of interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effective field theory for dark matter absorption on single phonons

Single phonon excitations, with energies in the 1–100 meV range, are a powerful probe of light dark matter (DM). Utilizing effective field theory, we derive a framework to compute DM absorption rates into single phonons starting from general DM-electron, proton, and neutron interactions. We apply the framework to a variety of DM models: Yukawa coupled scalars, axionlike particles with derivative interactions, and vector DM coupling via gauge interactions or Standard Model electric and magnetic dipole moments. We find that GaAs or Al 2 O 3 targets can set powerful constraints on a U ( 1 ) B − L model, and targets with electronic spin ordering are similarly sensitive to DM coupling to the electron magnetic dipole moment. Lastly, we make the code, (an extension of the existing code which computes general DM–single phonon scattering rates), publicly available. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hydrodynamic instabilities and heat transfer characteristics in the duct flow of a fluid in the supercritical thermodynamic regime

The behavior of fluids at supercritical thermodynamic conditions is inherently complex due to large variations in thermodynamic and transport properties. Recent numerical and experimental investigations illustrate ongoing interest for these fluids, especially supercritical CO 2 and supercritical water, for a variety of applications. For example, supercritical water reactors (SCWR) operate in this extreme condition of high-pressure and temperature, resulting in highly dynamic flow fields and unexpected heat transfer regimes. The potential heat transfer benefits in this regime are directly associated with the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. This work characterizes the hydrodynamic instabilities that arise for fluids at supercritical thermodynamic conditions when buoyancy forces are significant. Two specific configurations are considered, a natural convection cavity flow, and a mixed convection, heated, horizontal channel flow. Natural convection flow in a cavity is a classical configuration with expected behavior below the critical point. This configuration aids in characterizing the effect of the variable properties in the supercritical thermodynamic regime. Further, limited studies in the existing literature have been conducted for low-Reynolds and intermediate-Rayleigh numbers, mixed-convection channel flows for supercritical water, which is the focus of the channel flow configuration. To investigate the thermally driven hydrodynamic instabilities in this regime, a high-order fully-implicit numerical method is used. Such strong variations in thermophysical properties (in particular, density) are difficult to simulate and an altogether compressible framework is needed. Therefore, the compressible Navier-Stokes equations are solved without any additional assumptions. The fully implicit, high-order in space and time, reconstructed discontinuous Galerkin method as implemented within the multi-physics code called ALE3D (Arbitrary Lagrangian and Eulerian in 2D and 3D), developed at Lawrence Livermore National Laboratory (LLNL), is used. This fully implicit, L-stable method accurately captures the compressible nature of the ow in the limit of very low Mach number. It has been widely accepted that above the critical point, only one phase is observed. However, recent research has indicated the existence of the distinct gas-like and liquid-like regions separated by the Widom line, the locus of the maxima of the specific heat. Along the Widom line, density decreases 6-fold, viscosity drops by a factor of 2, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film boiling. A pseudo-film at the heated wall of the cavity and the horizontal channel is observed where buoyancy forces induce mixing through the specific configurations. Further the local Rayleigh and Richardson numbers provide maps of the flow field and the buoyancy forces driving the microscopic mixing. In the first chapter, I describe a background of supercritical fluid and the various applications. The second chapter focuses on the mathematical model and numerical method used for simulations, where a description of the equation of state for supercritical water is described. The third chapter focuses on the natural convection cavity with a heated bottom wall. In this cavity a gas-like and a liquid-like flow within the supercritical thermodynamic regime are observed. The fourth chapter focuses on a forced convection, horizontal channel, distinguishing between the gas-like, liquid-like, and mixed flow regimes. Mixed convection flow, with the addition of gravitational forces in the horizontal channel show the influence of variable properties on the hydrodynamic development, heat transfer, and rising instabilities. The last chapter of this research focuses on characterizing the unstable hydrodynamics through time-averaging processes and analysis of the movement of energy through the developing plumes.

42 ENGINEERING↗