Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parabolic”

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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Multigrid Reduction in Time for Chaotic Dynamical Systems

As CPU clock speeds have stagnated and high performance computers continue to have ever higher core counts, increased parallelism is needed to take advantage of these new architectures. Traditional serial time-marching schemes can be a significant bottleneck, as many types of simulations require large numbers of time-steps which must be computed sequentially. Parallel-in-time schemes, such as the Multigrid Reduction in Time (MGRIT) method, remedy this by parallelizing across time-steps and have shown promising results for parabolic problems. However, chaotic problems have proved more difficult, since chaotic initial value problems (IVPs) are inherently ill-conditioned. MGRIT relies on a hierarchy of successively coarser time-grids to iteratively correct the solution on the finest time-grid, but due to the nature of chaotic systems, small inaccuracies on the coarser levels can be greatly magnified and lead to poor coarse-grid corrections. Here we introduce a modified MGRIT algorithm based on an existing quadratically converging nonlinear extension to the multigrid Full Approximation Scheme (FAS), as well as a novel time-coarsening scheme. Together, these approaches better capture long-term chaotic behavior on coarse-grids and greatly improve convergence of MGRIT for chaotic IVPs. Further, we introduce a novel low-memory variant of the algorithm for solving chaotic PDEs with MGRIT which not only solves the IVP, but also provides estimates for the unstable Lyapunov vectors of the system. Finally, we provide supporting numerical results for the Lorenz system and demonstrate parallel speedup for the chaotic Kuramoto–Sivashinsky PDE over a significantly longer time-domain than in previous works.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Explicit Monotone Stable Super-Time-stepping Methods for Finite Time Singularities

We explore a novel way to numerically resolve the scaling behavior of finite-time singularities in solutions of nonlinear parabolic PDEs. The Runge–Kutta–Legendre (RKL) and Runge–Kutta–Gegenbauer (RKG) super-time-stepping methods were originally developed for nonlinear complex physics problems with diffusion. These are multistage single step second-order, forward-in-time methods with no implicit solves. The advantage is that the time-step size for stability scales with stage number 𝑠 as $\mathcal{O}$⁡(𝑠 2 ). Many interesting nonlinear PDEs have finite-time singularities, and the presence of diffusion often limits one to using implicit or semi-implicit time-step methods for stability constraints. Finite-time singularities are particularly challenging due to the large range of scales that one desires to resolve, often with adaptive spatial grids and adaptive time steps. Here, in this study, we show two examples of nonlinear PDEs for which the self-similar singularity structure has time and space scales that are resolvable using the RKL and RKG methods, without forcing even smaller time steps. Compared to commonly used implicit numerical methods, we achieve a significantly smaller run time while maintaining comparable accuracy. We also prove numerical monotonicity for both the RKL and RKG methods under their linear stability conditions for the constant coefficient heat equation, in the case of infinite domain and periodic boundary condition, leading to a theoretical guarantee of the superiority of the RKL and RKG methods over traditional super-time-stepping methods, such as the Runge-Kutta-Chebyshev and the orthogonal Runge-Kutta-Chebyshev methods. Code can be found at https://github.com/ZT220501/SRK-Singularity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

NSO Processing Scripts (NSO wind loading data set processing scripts) [SWR-23-86]

NREL conducted comprehensive field measurements of the atmospheric turbulent wind conditions and the resulting structural wind loads on parabolic troughs at the Nevada Solar One plant from October 2021 to June, 2023. The measurement set-up included meteorological masts and structural load sensors on four trough rows. Additionally, NREL commissioned a lidar scanning the horizontal plane over the trough field. The high-resolution data set, characterizing the complex flow field and resulting structural loads, is available at https://data.openei.org/submissions/5938. The processing routines published here were used to create the data from the instruments' raw data.

Egerer, Ulrike↗
DOE OSTI2026

SANDWake3D

SAND2026-23958O SANDWake3D models atmospheric surface layers and turbine wakes using a 3D parabolic RANS solver. This Python-based tool uses an efficient method to simulate wakes, capturing complex wake behaviors such as wind shear, veer, and atmospheric stratification, naturally. It has been calibrated and validated against higher-fidelity large eddy simulations (LES), providing reliable results for wind energy and atmospheric research applications. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗
DOE OSTI2021

Non-intrusive nonlinear model reduction via machine learning approximations to low-dimensional operators

Abstract Although projection-based reduced-order models (ROMs) for parameterized nonlinear dynamical systems have demonstrated exciting results across a range of applications, their broad adoption has been limited by their intrusivity: implementing such a reduced-order model typically requires significant modifications to the underlying simulation code. To address this, we propose a method that enables traditionally intrusive reduced-order models to be accurately approximated in a non-intrusive manner. Specifically, the approach approximates the low-dimensional operators associated with projection-based reduced-order models (ROMs) using modern machine-learning regression techniques. The only requirement of the simulation code is the ability to export the velocity given the state and parameters; this functionality is used to train the approximated low-dimensional operators. In addition to enabling nonintrusivity, we demonstrate that the approach also leads to very low computational complexity, achieving up to $$10^3{\times }$$ 10 3 × in run time. We demonstrate the effectiveness of the proposed technique on two types of PDEs. The domain of applications include both parabolic and hyperbolic PDEs, regardless of the dimension of full-order models (FOMs).

42 ENGINEERING↗
DOE OSTI2022

Modeling the hydrodynamic impact on the tool influence function during hemispherical subaperture optical polishing

To fabricate high-precision and accurate optics relative to the optical design surface, a high level of deterministic control of material removal (i.e., the tool influence function, TIF) during subaperture tool polishing is required. In this study, a detailed analysis of the pressure distribution, which is a key component of the TIF, has been performed using finite element analysis to couple together solid mechanics and fluid dynamics. Modeling experimental parameters of recently published work reveals that, when considering tool deformation, which in turn influences the fluid film thickness between the tool and workpiece, the effective pressure profile has a flat-top distribution. This flat-top pressure profile differs from the parabolic pressure distributions predicted by Hertzian mechanics. Furthermore, the shear contribution is shown here to be a key contributor to material removal, inducing the removal at the periphery of the contact edge and even outside the generally accepted contact area. Finally, the simulated fluid velocities provide evidence of mixed-mode contact polishing, supporting recent experimental findings that also suggest that onset of hydroplaning contributions lead to material removal drop-off.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Temporal reflection of an optical pulse from a short soliton: impact of Raman scattering

We study temporal reflection of an optical pulse from the refractive-index barrier created by a short pump soliton inside a nonlinear dispersive medium such as an optical fiber. One feature is that the soliton’s speed changes continuously as its spectrum redshifts because of intrapulse Raman scattering. We use the generalized nonlinear Schrödinger equation to find the shape and spectrum of the reflected pulse. Both are affected considerably by the soliton’s trajectory. The reflected pulse can become considerably narrower compared to the incident pulse under conditions that involve a type of temporal focusing. This phenomenon is explained through space–time duality by showing that the temporal situation is analogous to an optical beam incident obliquely on a parabolic mirror. We obtain an approximate analytic expression for the reflected pulse’s spectrum and use it to derive the temporal version of the transformation law for the

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Numerical investigation of gas-filled multipass cells in the enhanced dispersion regime for clean spectral broadening and pulse compression

We show via numerical simulations that the regime of enhanced frequency chirp can be achieved in gas-filled multipass cells. Our results demonstrate that there exists a region of pulse and cell parameters for which a broad and flat spectrum with a smooth parabolic-like phase can be generated. This spectrum is compatible with clean ultrashort pulses, whose secondary structures are always below the 0.5% of its peak intensity such that the energy ratio (the energy contained within the main peak of the pulse) is above 98%. This regime makes multipass cell post-compression one of the most versatile schemes to sculpt a clean intense ultrashort optical pulse.

Segundo Staels, V. W. (ORCID:0000000200056157)↗
DOE OSTI2021

ADVANCED THERMAL-HYDRAULIC MODEL OF HEAT RECOVERY STEAM GENERATORS

Steam-side oxidation at elevated temperatures in heat recovery steam generators (HRSGs) can lead to overheating of local tubes, potential oxide scale exfoliation, and damage to critical equipment. Predicting the oxide growth distribution requires detailed knowledge of the distribution of tube metal temperatures. To address these challenges, a computational fluid dynamics (CFD) model of a high-pressure superheater (HPSH) was developed within the Siemens’ STAR-CCM+ CFD commercial software. The HPSH section of the steam generator was selected for modeling because of its high oxidation rates.The model development included several steps. In the first step, an oxidation model based on a parabolic oxidation relationship was implemented. The oxidation growth predicted by the model was then verified for a full-load boiler operation case. In the second step, a three-dimensional CFD model was developed to account for the thermal and flow physics involved in HRSGs. The CFD model was then verified and validated on a simplified geometry that contained six row-tubes which mimicked the tube arrangements of an HPSH section. Numerical results compared well with published experimental results. In the third step, the oxidation model was coupled to the HRSG’s CFD model, and a comparison study was performed for two cases: (a) regular finless tubes and (b) finned tubes. Preliminary results showed that the oxide thickness predicted by the case with tube fins was higher compared to that of the finless tube case, thus demonstrating the need to model the effects of the tube fins on oxide thickness.

Panicker, Nithin↗
DOE OSTI2021

Renewable Thermal Energy Systems Designed for Industrial Process Solutions in Multiple Industries

The need for renewable heat in industry is vital for the next decade and beyond. Industrial decarbonization is a key area that must be accelerated, to foster the removal of fossil fuels from the provision of heat, especially at low temperatures. This paper looks at the development and results of case studies for understanding the economics and potential for renewable thermal energy systems (RTES), particularly in hybrid configurations to provide industrial process heat (IPH). For the case studies, these include non-concentrating e.g., heat pumps, and concentrating collectors e.g., parabolic trough collectors and direct steam generation (DSG)-linear Fresnel collectors (LFCs). The results show that a levelized cost of heat (LCOH) of $6-$8 per million British Thermal Units (MMBTU) is possible, depending on the direct normal irradiance (DNI) and the system sizing e.g., to increase the solar fractions. In Arizona for example, with a DNI of 7.36 kWh/m2/day, the base case for the DSG-LFC system with 6hrs of thermal energy storage could potentially meet a 1 MWth load 80% of the year.

concentrated solar thermal↗
DOE OSTI2024

Single-gap isotropic s- wave superconductivity in single crystals AuSn 4

In this article, London, λ L (T), and Campbell, λ C (T), penetration depths were measured in single crystals of a topological superconductor candidate AuSn 4 . At low temperatures, λ L (T) is exponentially attenuated and, if fitted with the power law, λ(T) ~ T n , gives exponents n > 4, indistinguishable from the isotropic single s- wave gap Bardeen-Cooper-Schrieffer (BCS) asymptotic. The superfluid density fits perfectly in the entire temperature range to the BCS theory. The superconducting transition temperature, T c = 2.40 ± 0.05 K, does not change after 2.5 MeV electron irradiation, indicating the validity of the Anderson theorem for isotropic s- wave superconductors. Campbell penetration depth before and after electron irradiation shows no hysteresis between the zero-field cooling (ZFC) and field cooling (FC) protocols, consistent with the parabolic pinning potential. Interestingly, the critical current density estimated from the original Campbell theory decreases after irradiation, implying that a more sophisticated theory involving collective effects is needed to describe vortex pinning in this system. In general, our thermodynamic measurements strongly suggest that the bulk response of the AuSn 4 crystals is fully consistent with the isotropic s- wave weak-coupling BCS superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

3D Printing Prototypes for a Reflector-Based Broadband Axion Search

This work sought to begin research and development on a horn antenna which can be used to transfer the energy of microwave photons produced by axions in a strong magnetic field to a coaxial cable for measurement. 3D prototypes of pyramidal horn antennas and a parabolic reflector were printed and then coated with conductive spray paint. The S 11 parameter was measured and simulated for both 3D printed horn antenna and the conical horn antenna both with and without the reflector present to determine how much of the radiated power was reflected back. Ultimately, for both the conical and pyramidal antennas, less than one percent of the energy emitted by the antennas was reflected back. When the pyramidal antenna was tilted by about 45 degrees, however, around 30-40 percent of the microwaves emitted by the antennas was reflected back.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Sodium Ion Expansion Power Block for Distributed CSP

The Sodium Ion Expansion Power Block for Distributed CSP was a three-plus-one-year effort under the Concentrating Solar Power: Advanced Projects Offering Low LCOE Opportunities (CSP: APOLLO) funding program within the U.S. Department of Energy Solar Energy Technologies Office. The primary objective of this project is to develop a dual-stage modular sodium thermal electrochemical converter (Na-TEC) heat engine power block, which can be potentially integrated with either a small-scale dish solar or large-scale heliostats and parabolic trough CSP. Na-TEC is a heat engine that generates electricity through the isothermal expansion of sodium ions. The Na-TEC is a closed system that can theoretically achieve conversion efficiencies above 45% when operating between thermal reservoirs at 1150 K and 550 K. However, thermal designs have confined previous single-stage devices to thermal efficiencies below 20%. To mitigate some of these limitations, we consider dividing the isothermal expansion into two stages; one at the evaporator temperature (1150 K) and another at an intermediate temperature (650 K –1050 K). This dual-stage Na-TEC takes advantage of regeneration and reheating, and could be amenable to better thermal management. In light of this, we first designed and developed a thermo-electrochemical model, and thermodynamically demonstrated how the dual-stage device can improve the efficiency by up to 8% points over the best performing single-stage device. We also established an application regime map for the single- and dual-stage Na-TEC in terms of the power density and the total thermal parasitic loss. Moreover, a thermal design of an axisymmetric dual-stage Na-TEC is developed to guide the scale-up and fabrication of sub-components of prototype module. A reduced-order finite-element model is used in conjunction with a Na-TEC thermodynamic model that was developed to determine the total parasitic heat loss of this dual-stage design. A number of simplifications are applied in the reduced-order model to decrease the computational time while maintaining acceptable accuracy. According to this analysis, a maximum efficiency of 29% and a maximum power output of 125 W can be achieved. Ultimately, we were able to demonstrate thermal efficiency improvements of the Na-TEC heat engine from 19% up to 40.3%, in a dual-stage (non-optimized) prototype module that we designed, fabricated, and tested with high temperature stage at 923 K. Furthermore, a cost-performance analysis for this improved dual-stage design was carried out for distributed-CSP systems. A high-level techno-economic analysis (TEA) explores four scenarios where a Na-TEC is used as the heat engine for a distributed-CSP system. Overnight capital cost and levelized cost of electricity (LCOE) are estimated for a system lifetime of 30 years, revealing that overnight capital costs in a range from $3.57 to $17.71 per We are feasible, which equate to LCOEs from 6.9 to 17.2 cents/kWh e -1 . This analysis makes a significant contribution by concurrently quantifying the efficiency and unit costs for a range of multistage configurations, and demonstrating that a Na-TEC may be a promising alternative to Stirling engines for distributed-CSP systems at residential scale of 1–5 kW e .

14 SOLAR ENERGY↗
DOE OSTI2020

Enhancement of Optical Efficiency of CSP Mirrors for Reducing O&M Cost via Near-Continuous Operation of Self-Cleaning Electrodynamic Screens (EDS). Final Report

Over the past decade, techno-economical advancements in solar energy systems, particularly the improvement in the conversion efficiency of PV modules made with mono-crystalline silicon solar cells from 12% to 20%, as well as the cost reduction in manufacturing by a factor of about 10%, have made it possible to achieve a levelized cost of electricity (LCOE) in PV plants that is comparable to, or less than, the cost of deriving electricity from fossil fuels. Operating PV plants in mid-latitude sun-belt regions, where the solar irradiance level is highest, provides a high annual energy-yield (kWh/kWp) due to two factors: (1) the availability of predictable high solar irradiance throughout the year with the fewest interruptions in solar flux from clouds and rain, and (2) the increased conversion efficiency of crystalline solar cells, as recombination loss has decreased with increased intensity of the sunlight that illuminates the silicon solar cells. Semi-arid and desert regions, however, are plagued by high atmospheric dust concentrations and frequent sand storms. The deposition of a layer of dust on the optical surfaces of solar collectors such as PV modules and concentrating mirrors reduces the transmission efficiency of sunlight that actually reaches the solar cells or receivers, resulting in high energy-yield soiling loss. There are two major cost components to operating a solar plant: (1) installation costs, and (2) operation-and-maintenance (O&M) costs. There is no fuel cost; hence operating a solar plant in a semi-arid or desert region provides high returns on investment if soiling losses are mitigated via efficient cleaning methods and optimized cleaning frequency. If solar collectors are not cleaned, the accumulation of dust layers on solar collectors may cause the operation of such plants in arid regions to become economically unviable. Washing solar collectors with water and detergent, as is most commonly done now, is an efficient method for cleaning. The conventional approach in utility-scale solar plants is to use a large truck with a water tank and pump system for spraying deionized water on the surface of the solar collectors. Robotic cleaning with brushes, used for many solar plants, requires lesser water for cleaning. The water consumed using semi-automated cleaning of PV modules in utility-scale solar plants is approximately 2 liters/m 2 per cleaning cycle. The total optical surface area of the solar collectors in 1 TW-scale solar installation will be more than 3 × 109 m 2 ; hence an enormous amount of water be needed for cleaning. There simply is not enough fresh water in the sun-belt areas of the world for predicted cleaning needs. In solar power plants, the estimated cost of cleaning solar collectors includes expenses related to the equipment used, labor, cost of transportation of water, the energy required for cleaning, plus ancillary costs whereas the cost of water is not considered. The water used is obtained from sources located close to plant sites, unmindful of the environmental and societal impact as the power plants are oftentimes located in regions that face severe drought. This practice is very similar to the cost calculations in deriving the levelized cost of electricity (LCOE) in conventional power plants based on burning fossil fuel such as coal and gas, while disregarding the cost of climate change and health effects. Unless a water-free or low-water cleaning method is established, the expansion of solar plants may lose public support in areas suffering long intervals of drought. The goal of this research project has been the development and application of the Electrodynamic Screen (EDS) as a means for a water-free, scalable cleaning process applicable to solar-power installations, including rooftop applications. We describe here the development of an EDS film-based cleaning process as an emerging method for use on PV modules, parabolic troughs, and heliostats. This report aims to show the feasibility of integrating or retrofitting EDS films onto the optical surfaces of solar collectors (both PV and CSP) while maintaining high transmission or reflection efficiency. The cleaning action provided by the EDS film is an active method to remove dust deposits by electrodynamic force. Current lab-scale prototype EDS films, retrofitted onto solar panels and mirrors, have shown to be capable of maintaining optical transmission or specular-reflection efficiencies higher than 90% of initial values under clean conditions. The optical surfaces of solar collectors laminated with EDS films can remove more than 90% of deposited dust when the EDS is activated for less than two minutes. As an electrodynamic dust removal process, the EDS film-based method is designed primarily for the removal of dust in solar installations located in semi-arid and desert areas, where the atmosphere is often dry and dusty and rainfall is infrequent. While EDS film application minimizes water consumption and facilitates cleaning as frequently as needed, it has limitations in removing contaminants such as soot, organic pollutants deposited as fine films on the surface, and bird droppings. The dust removal efficiency of the EDS is maximum at relative humidity RH is 40 to 50% and decreases at when RH > 65%. Many solar plant sites undergo diurnal and seasonal cycles of high ambient, early afternoon temperatures and an RH that reaches the dew point early in the morning. These variations in atmospheric conditions do not limit the operation of the EDS film. This report presents a brief review of the progress and the potential of EDS film technology for mitigating the impact of dust on solar collectors via water-free cleaning, as well as current technical challenges regarding efficiency and durability. Our experimental data on the performance of EDS films show that: (1) the dust-removal efficiency (DRE) can reach levels higher than 90%, (2) the specular reflectivity (SR) of EDS film-laminated second-surface mirrors reach levels in excess of 90%, (3) the specular reflectivity restoration (SRR) can exceed 90%, (4) the output-power restoration (OPR) of PV modules can exceed 95%, and (5) the optical transmission efficiency (TE) of the EDS films can be greater than 90%. Working with Sandia National Laboratories (NM), Corning Research and Development Corporation (NY), Eastman Kodak (NY), Tomark-Worthen Industries (NH), and EDS Chile SPA (Chile), we have produced EDS film-laminated PV modules and demonstrated their self-cleaning functions without requiring water. We have demonstrated that the operational range of EDS films will cover the expected ambient temperature of solar fields at RH cycling varying from 20 to 95% as long as the EDS films are activated in the RH range 20 to 50%. (Typical solar-field climates in deserts and semi-arid lands often reach near dew point in coastal areas.) Our experiments on the application of hydrophobic-fluorinated nanoparticle coatings on EDS film surfaces show that the EDS operational range can be extended to higher RH levels that approach the dew point. At Eastman Kodak, as one of our industrial partners, we were able to establish a process for manufacturing EDS films using flexographic printing of the electrodes onto transparent polymer films. This process utilizes an existing manufacturing line at Eastman Kodak that allows fabrication of medium-scale EDS films (26 cm × 30 cm). The manufacturing process has the capacity to produce EDS films at high production speeds. The medium-scale EDS films that have been printed at Kodak were evaluated in the lab at Boston University. The EDS films produced at Kodak were laminated at Tomark-Worthen using an industrial scale vacuum laminator to produce EDS film stacks, which can be affixed onto the optical surfaces of PV modules or concentrating mirrors. The EDS film stack consists of the EDS films that have Willow® Glass (WG) which has a thickness of 100 μm as the front surface. The WG sheets obtained from Corning Research and Development Corporation are customized in size and shape to cover the active area of the EDS films. The back surface of the EDS film stack is integrated onto the optical surface of the solar panel or mirror using optically clear adhesive (OCA) films or silicone adhesives. The OCA films (thickness 25 μm) are produced by 3M. The EDS film stacks have the architecture: WG/OCA/EDS Film/OCA/ over PV module or solar mirror. The power-supply units needed for activating the electrodes of the EDS film were designed and produced at Boston University. These power supply units provide three-phase, 1.2 kV voltage pulses at a very low current (micro-ampere) level and at a low frequency (≈ 5 Hz). The voltage pulses are applied to the electrodes in a sequence such that the train of pulses resembles a unidirectional traveling wave of electrical field on the surface of the EDS film. The dust particles on the surface become charged electrostatically and are levitated by the Coulomb force. The lateral sweeping action of the traveling electric field created by the three-phase voltages pulses then sweeps the dust off the surface. The energy consumed by the EDS electrodes is less than 0.2 Wh/m 2 /cleaning cycle, enabling energy-efficient restoration of output power (OPR) of PV module or specular reflectivity restoration (SRR) for solar mirrors. The EDS system consists of (1) an EDS film stack laminated onto the solar collectors, (2) connection of the EDS film to its power supply unit and (3) Interconnection of the power supply to PV modules or solar mirrors. Design, construction and assessment of field-testing units that have EDS stack laminated PV modules for evaluating the performance of the EDS films in solar fields is being carried out at BU. Our progress under this project, aimed at the advancement of EDS film technology, has reached DOE Technology Readiness Level (TRL) 6. Based on the extensive laboratory evaluations and limited field trials, as well as contacts with potential users, we believe that the technology has reached its commercial stage. We are conducting a cost analysis using the National Renewable Energy Laboratory (NREL) System Advisory Model (SAM) and are preparing for field trials of EDS films in different solar fields in the US, Chile, India and in the Middle East. A brief description of EDS film performance, construction and testing of the field-test unit, autonomous operation of the field-test unit for evaluating EDS performance in increasing energy yield, associated revenue savings, and water conservation is presented.

14 SOLAR ENERGY↗
DOE OSTI2020

Plasma Photonics for the Control of Kilojoule Lasers

The main goal of this project was to develop plasma photonic devices (plasma lenses or ellipsoidal plasma mirrors (EPMs)) for controlling the focus of high-energy (kJ-class) laser systems. In order to defocus a kJclass laser beam, a plasma optic must produce a very specific plasma density profile. For a plasma lens, this profile must be peaked on axis and fall off parabolicly with radius. We developed four plasma-lens target designs that were predicted to yield the ideal plasma profile required for defocusing kJ-class laser beams. Particle-in-cell (PIC) simulations using the code OSIRIS were used to explore the design. The resulting plasma lens must also be underdense, which sets an upper limit on the solid density of the plasma-lens targets to ~ 1 mg/cm 3 or less. This low of a density is highly challenging to fabricate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Updated Estimate of Tritium Permeation from TPBAR Disposal Containers in ILV (U)

A tritium source term analysis was performed for TPBAR disposal in the E-Area Intermediate Level Vaults for the E-Area Low-Level Waste Facility Performance Assessment (PA). This analysis is based on an earlier source term analysis which treated the bulk oft he tritium residual as tightly bound by the TPBAR getter material, with only a small fraction existing as tritiated moisture in the lithium aluminate ceramic pellets. Together with atmospheric moisture trapped in the free volume, the tritiated water vapor is assumed to corrode steel surfaces inside the disposal container, covering them with a magnetite film while generating hydrogen. The carbon steel walls of the disposal container are permeable to hydrogen, providing a path way for tritium to escape containment. The rate of hydrogen generation is assumed to be limited by the rate of corrosion, which is assumed to be governed by parabolic reaction kinetics obtained from the literature. This relies on the further assumption that the water vapor consumed by the corrosion reaction is continually replaced by moisture from the lithium aluminate pellets until all of that moisture is gone. In addition to tritium permeation, the analysis also includes a slow leak through the disposal container walls at the maximum allowable leak rate, 1 x 1CH standard cm3/s. Results were obtained for four different combinations of internal and container wall temperatures, established in an earlier thermal analysis that considered two different vault loadings and two different TPBAR activity levels. Instantaneous release rates for these four cases are plotted in Figure ES-1 below

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Laser Network US Experiments at the University of Texas at Austin

Under Subcontract B638654, The University of Texas at Austin (UT Austin) Center for High Energy Density Science supported three Laser Network US experiments and conducted upgrade planning for the Texas Petawatt Laser Facility. For the experiment “Investigation of the physics of MeV X-ray and particle generation using Compound Parabolic Concentrators” led by Dr. Andrew Mackinnon of Lawrence Livermore National Laboratory, UT Austin provided postdoctoral associate support in the person of Dr. Ganesh Tiwari. For the experiments “Megatesla magnetic fields and enhanced gamma-ray beam generation” led by Hans Rinderknecht of LLE, University of Rochester (experiment) and Alexey Arefiev of UC San Diego (theory) and “High resolution imaging of resonately accelerated electrons under tight focusing” led by Deepak Kumar of ELI-Beamlines, Czech Republic, UT Austin provided the Texas Petawatt Laser staff and a prorated share of the materials and supplies approximately equal to the average of those required for one 5-week experimental campaign. Summaries of these three experiments are provided in this report. Two staff members were hired—Dr. Sandra Bruce at 40% and Dr. Erhard Gaul at 50% time—to support Texas Petawatt Laser maintenance and operations and to plan for the Texas Relativistic Intensity Ultrafast Multi-beam Plasma Facility (TRIUMPF).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗