Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mass density ratio”

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.

882 records · Page 49

Detonation to Deflagration-Mode Transitions in Pulsed Plasma Accelerators

PULSED plasma accelerators typically operate by storing energy in a capacitor bank and then discharging this energy through a gas. The current in such discharges will ionize the gas and produce a strong magnetic field, which interacts with the flowing current to accelerate the plasma through the Lorentz body force. For the present work, two plasma accelerator types employing this general scheme are of interest: the gas-fed pulsed plasma thruster (PPT) 1 and the quasi-steady magnetoplasmadynamic (MPD) accelerator. 2 The gas-fed pulsed plasma accelerator is generally understood as a completely transient device discharging in ∼1-10 μs. When the capacitor bank is discharged through the gas, a current sheet forms at the breech of the thruster and propagates forward under a j × B body force, entraining and accelerating propellant it encounters. This process is sometimes referred to in literature as 'snowplowing' the propellant or accelerating the gas in a detonation-mode because the current sheet representation approximates that of a strong detonation shockwave propagating through the gas. For these devices, acceleration of the initial current sheet ceases when either the current sheet reaches the end of the device and is ejected or when the current in the circuit reverses, striking a new 'crowbar' discharge at the breech and depriving the initial sheet of additional acceleration. In general, PPTs typically claim thrust efficiencies (ratio of jet kinetic energy to input electrical energy) registering in the teens or lower. 3 In the quasi-steady MPD accelerator, the pulse is lengthened to ∼1 ms or longer and maintained at an approximately constant level during discharge through the use of a pulse-forming network (PFN) of capacitors. After an initial transient discharge, which is typically short relative to the overall discharge period, the plasma assumes a relatively steady-state configuration, known as 'quasi-steady' MPD operation. 2 In this state, ionized gas flows through a stationary current channel in a manner that is sometimes referred to as deflagration-mode operation owing to the similarities to deflagration waves in gases. The plasma experiences electromagnetic acceleration as the plasma flows through the current channel towards the exit of the device. Quasi-steady MPD thrusters claim efficiencies up to 50% for certain propellants.4 There has been significant and sustained research over several decades on both gas-fed PPTs and quasi-steady MPD thrusters, however there have been pulsed thrusters that do not appear to exactly fit either classification, instead possessing a mixture of operational qualities characteristic of both thruster variants. The Coaxial High ENerGy (CHENG) thruster by Cheng, et al.5 operated on the short 10 μs timescales characteristic of PPTs, but claimed the high thrust densities, high efficiencies, and low electrode erosion rates that are more consistent with the MPD/deflagration mode of plasma acceleration. Gas-fed PPT research by Ziemer, et al. 3, 6 identified two separate regimes of performance in those thrusters. The regime at higher mass bits (termed Mode I in that work) possessed relatively constant thrust efficiency as a function of mass bit, while the second regime at very low mass bits (termed Mode II) exhibited an increase in efficiency with decreasing mass bit. Work by Poehlmann et al.7 and by Sitaraman and Raja8 sought to understand the performance of the CHENG thruster and the Mode I/Mode II performance in PPTs by modeling the acceleration using the Hugoniot Relation, with the detonation and deflagration modes of plasma acceleration representing two distinct sets of solutions to the relevant conservation laws. In these works, it was proposed that the values of the various controllable parameters determined whether the accelerator would operate in detonation or deflagration mode. Our hypothesized view of the acceleration process in the CHENG thruster and in PPTs experiencing a transition from Mode I to Mode II is inspired by observations of the transition from the PPT mode of operation to the quasi-steady MPD mode. Specifically, the quasi-steady MPD was discovered by driving a PPT to extended pulse lengths. Above a certain pulse length threshold the transient plasma current sheet transitions into a stable plasma acceleration mode that 'replicates in every observable detail steady flow self-field magnetoplasmadynamic acceleration.' 9 In the present work, instead of treating the accelerator as if it were only operating in a single mode during a pulse, we consider the initial stage of the discharge in all cases as a current sheet forming at the breach of the accelerator and moving towards the exit as a detonation wave. If the current sheet reaches the exit of the accelerator before the discharge is completed, the view of the acceleration mode transitions to the deflagration mode-type found in quasi-steady MPD thrusters. In previous work10 we presented a modeling framework that first captured the time-evolution of the current sheet (detonation) mode of the thruster and then transitioned into the quasi-steady MPD (deflagration) mode of plasma acceleration. In the present work, variations of the controllable parameters - specifically the pulsed circuit properties, the amount of mass injected into the thruster, and the relative timing between the initial gas injection and the initiation of the plasma current sheet - will be used to explore the thruster performance. A range of parameters are explored to demonstrate that standard gas-fed pulsed plasma accelerators, the CHENG thruster, and the quasi-steady MPD accelerator are variations of the same device, with the overall acceleration of the plasma depending upon the behavior of the plasma discharge during initial transient phase and the relative lengths of the detonation and deflagration modes of operation.

Polzin, Kurt A.↗

Development of a Testbed for Pipeline Flushing - 20425

Pipelines which carry high-level waste within the DOE complex should be properly flushed in the cases where stationary or moving beds of solid sediment occur, or lines are prone to hydrogen gas buildup. Present guidelines establish a minimum for the amount of water (flush volume) and flush velocity values used for post-transfer flushing operations to achieve satisfactory cleaning of pipelines. However, further studies are needed to find optimal operating velocity modes/values in order to minimizing the flush volume and consequent downstream waste. These efforts are significantly helpful to DOE waste remediation sites by preserving tank storage, preventing additional waste processing, and minimizing dilution and changes in waste chemistry. An experimental test loop was recently developed for study of non-Newtonian slurry flushing at the Florida International University (FIU). This loop was designed to create sediment beds of various materials and bed heights and investigate parameters that effect the efficiency of flushing operations. The objective is to find flush velocity values/modes which lead to satisfactory cleaning of transport lines with a minimum amount of water usage. A 165 ft experimental test loop made of 3-inch carbon steel pipes was constructed and used in initial testing campaigns. This pipe loop was equipped with real-time monitoring instruments to record variation of pressure, density, mass flow rate, and temperature during testing. Post-flush in-situ no-flow evaluations were performed using ultrasonic, endoscopy, and visual inspection. Additional post-flush evaluations were conducted using filtration and density monitoring of residuals within a designated circulations loop. These provisions enable comparison of pipeline cleanness before the start and after termination of flushing operations for efficiency evaluation. This paper presents efforts associated with flushing of kaolin-water mixtures at different concentrations and with various flush volume and flush velocity values/modes. In initial testing, efforts focused on creating repeatable sediment beds inside the pipeline in both fully-flooded (no pre-flush drainage) and gravity-drained (pre-flush drained) systems. Bed characterization (i.e., measurement of sediment height and solids concentration) was conducted to ensure about consistency of initial conditions between tests. Preliminary results for flushing study of 10% vol. kaolin-water mixtures starting from a fully-flooded initial condition showed that flushing with 1-line volume (65 gallons) could remove up to 72 percent of solids from the pipeline. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Galaxy cluster matter profiles - I. Self-similarity, mass calibration, and observable-mass relation validation employing cluster mass posteriors

We present a study of the weak lensing inferred matter profiles ΔΣ(R) of 698 South Pole Telescope (SPT) thermal Sunyaev-Zel’dovich effect (tSZE) selected and MCMF optically confirmed galaxy clusters in the redshift range 0.25 < z < 0.94 that have associated weak gravitational lensing shear profiles from the Dark Energy Survey (DES). Rescaling these profiles to account for the mass dependent size and the redshift dependent density produces average rescaled matter profiles ΔΣ(R/R200c)/(ρcritR200c) with a lower dispersion than the unscaled ΔΣ(R) versions, indicating a significant degree of self-similarity. Galaxy clusters from hydrodynamical simulations also exhibit matter profiles that suggest a high degree of self-similarity, with RMS variation among the average rescaled matter profiles with redshift and mass falling by a factor of approximately six and 23, respectively, compared to the unscaled average matter profiles. We employed this regularity in a new Bayesian method for weak lensing mass calibration that employs the so-called cluster mass posterior P(M200|ζ̂, λ̂, z), which describes the individual cluster masses given their tSZE (ζ̂) and optical (λ̂, z) observables. This method enables simultaneous constraints on richness λ-mass and tSZE detection significance ζ-mass relations using average rescaled cluster matter profiles. We validated the method using realistic mock datasets and present observable-mass relation constraints for the SPT×DES sample, where we constrained the amplitude, mass trend, redshift trend, and intrinsic scatter. Our observable-mass relation results are in agreement with the mass calibration derived from the recent cosmological analysis of the SPT×DES data based on a cluster-by-cluster lensing calibration. Our new mass calibration technique offers a higher efficiency when compared to the single cluster calibration technique. We present new validation tests of the observable-mass relation that indicate the underlying power-law form and scatter are adequate to describe the real cluster sample but that also suggest a redshift variation in the intrinsic scatter of the λ-mass relation may offer a better description. In addition, the average rescaled matter profiles offer high signal-to-noise ratio (S/N) constraints on the shape of real cluster matter profiles, which are in good agreement with available hydrodynamical ΛCDM simulations. This high S/N profile contains information about baryon feedback, the collisional nature of dark matter, and potential deviations from general relativity.Key words: gravitational lensing: weak / galaxies: clusters: general / large-scale structure of Universe

79 ASTRONOMY AND ASTROPHYSICS↗

Population of Merging Compact Binaries Inferred Using Gravitational Waves through GWTC-3

We report on the population properties of compact binary mergers inferred from gravitational-wave observations of these systems during the first three LIGO-Virgo observing runs. The Gravitational-Wave Transient Catalog 3 (GWTC-3) contains signals consistent with three classes of binary mergers: binary black hole, binary neutron star, and neutron star–black hole mergers. We infer the binary neutron star merger rate to be between 10 and 1700 Gpc − 3 yr − 1 and the neutron star–black hole merger rate to be between 7.8 and 140 Gpc − 3 yr − 1 , assuming a constant rate density in the comoving frame and taking the union of 90% credible intervals for methods used in this work. We infer the binary black hole merger rate, allowing for evolution with redshift, to be between 17.9 and 44 Gpc − 3 yr − 1 at a fiducial redshift ( z = 0.2 ). The rate of binary black hole mergers is observed to increase with redshift at a rate proportional to ( 1 + z ) κ with κ = 2. 9 − 1.8 + 1.7 for z ≲ 1 . Using both binary neutron star and neutron star–black hole binaries, we obtain a broad, relatively flat neutron star mass distribution extending from 1.2 − 0.2 + 0.1 to 2.0 − 0.3 + 0.3 M ⊙ . We confidently determine that the merger rate as a function of mass sharply declines after the expected maximum neutron star mass, but cannot yet confirm or rule out the existence of a lower mass gap between neutron stars and black holes. We also find the binary black hole mass distribution has localized over- and underdensities relative to a power-law distribution, with peaks emerging at chirp masses of 8.3 − 0.5 + 0.3 and 27.9 − 1.8 + 1.9 M ⊙ . While we continue to find that the mass distribution of a binary’s more massive component strongly decreases as a function of primary mass, we observe no evidence of a strongly suppressed merger rate above approximately 60 M ⊙ , which would indicate the presence of a upper mass gap. Observed black hole spins are small, with half of spin magnitudes below χ i ≈ 0.25 . While the majority of spins are preferentially aligned with the orbital angular momentum, we infer evidence of antialigned spins among the binary population. We observe an increase in spin magnitude for systems with more unequal-mass ratio. We also observe evidence of misalignment of spins relative to the orbital angular momentum. Published by the American Physical Society 2023

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Asymptotic Giant Branch stars as a source of short-lived radioactive nuclei in the solar nebula

We carried out a theoretical evaluation of the contribution of Asymptotic Giant Branch (AGB) stars to some short-lived (10(exp 6) less than or equal to Tau-bar less than or equal to 2 x 10(exp 7) yr) isotopes in the Interstellar Medium (ISM) and in the early solar system using stellar model calculations for thermally pulsing evolutionary phases of low-mass stars. The yields of s-process nuclei in the convective He-shell for different neutron exposures tau(sub 0) were obtained, and AGB stars were shown to produce several radioactive nuclei (especially Pd-107, Pb-205, Fe-60, Zr-93, Tc-99, Cs-135, and Hf-182) in diferent amounts. Assuming either contamination of the solar nebula from a single AGB star or models for continuous injection and mixing from many stars into the ISM, we calculate the ratios of radioactive to stable nuclei at the epoch of the Sun's formation. The dilution factor between the AGB ejecta and the early solar system matter is obtained by matching the observed Pd-107/Pd-108 and depends on the value of tau(sub 0). It is found that small masses M(sub He) of He-shell material (10(exp -4)-10(exp -7) solar mass) enriched in s-process nuclei are sufficient to contaminate 1 solar mass of the ISM to produce the Pd-107 found in the early solar system. Predictions are made for all of the other radioactive isotopes. The optimal model to explain several observed radioactive species at different states of the proto-solar nebula involves a single AGB star with a low neutron exposure (tau(sub 0) = 0.03 mbarn(sup -1)) which contaminated the cloud with a dilution factor of M(sub He)/solar mass approximately 1.5 x 10(exp -4). This will also contribute newly synthesized stable s-process nuclei in the amount of approximately 10(exp -4) of their abundances already present in the proto-solar cloud. Variations in the degree of homogenization (approximately 30%) of the injected material may account for some of the small general isotopic anomalies found in meteorites. It is also found that Fe-60 is produced in small but significant quantities that may be sufficient to explain the observations if the time elapsed delta from the contamination of the ISM to the formation of protoplanetary bodies is not higher than delta = 5 x 10(exp 6) yr. If delta is longer, up to 10 x 10(exp 6) yr, this would require the single AGB star to experience enhanced neutron densities (n(sub n) approximately 3 x 10(exp 9)n/cu cm) in the s-processing zone in order to compensate for the branching at Fe-59. The alternative model of long-term continuous ejection of matter from many AGB stars does not appear to match the observations. We also estimate the Al-26 production from the H-shell and find that the Al-26 abundance in the early solar system may be readily explained in a self-consistent manner. Moreover, Al-26 from AGB stars may contribute substantially to the galactic Al-26 gamma-source, while no significant gamma-flux from Co-60 (deriving from Fe-60 decay) is to be expected.

Wasserburg, G. J.↗

Boson localization and universality in YBa2Cu(3-x)M(x)O(7-delta)

We consider a two component mixture of charged fermions on neutralizing background with all sign combinations and arbitrarily small mass ratios. In the two impurity limit for the heavier component we show that the pair forms a bound state for all charge combinations. In the lowest order approximation we derive a closed form expression Veff(r) for the binding potential which has short-range repulsion followed by attraction. In the classical limit, when the mass of embedded particles is large m2 much greater than m, we can calculate from Veff(r) also the cohesive energy E and the bond length R of a metallic crystal such as lithium. The lowest order result is R = 3.1 A, E = -0.9 eV, not entirely different from the experimental result for lithium metal. The same interaction for two holes on a parabolic band with m2 greater than m gives the quantum mechanical bound state which one may interpret as a boson or local pair in the case of high-Te and heavy fermion superconductors. We also show that for compounds of the type YBa2Cu(3 - x)M(x)O(7 - delta) one can understand most of the experimental results for the superconducting and normal states with a single temperature dependent boson breaking function f(T) for each impurity content x governing the decay of bosons into pairing fermions. In the normal state f(T) turns out to be a linear, universal function, independent of the impurity content I and the oxygen content delta. We predict with universality a depression in Tc(x) with slight down bending in agreement with experiment. As a natural consequence of the model the bosons become localized slightly above Tc due to the Wigner crystallization, enhanced with lattice local field minima. The holes remain delocalized with a linearly increasing concentration in the normal state, thus explaining the rising Hall density. The boson localization temperature T(sub BL) shows up as a minimum in the Hall density R(sub ab)(exp -1). We also give explanation for very recently observed scaling of temperature dependent Hall effect in La(2 - x)Sr(x)CuO4.

Kallio, A.↗

The ISM in the M82 starburst

We have observed (O 1) (63 microns) and (Si 2) (35 microns) in the central 700 pc of the starburst galaxy M82. The luminosities in these transitions are 7.1 x 10(exp 7) solar luminosity and 6.2 x 10(exp 7) solar luminosity, respectively, which are each approx. 0.15% of the bolometric luminosity from this region. The ratios of (O 1) line luminosity to (O 3), (Si 2) (35 microns), and to bolometric luminosities in M82 are similar to those in M42, M17, and Sgr A. These similarities, and the association of the bulk of the (O 1) and (Si 2) emission with the ionized emission, suggest that the dominant emission mechanism for (O 1) and (Si 2) in M82 is the same as in these Galactic regions, namely warm gas photodissociated by UV flux from the OB stars responsible for the nearby H 2 regions. We argue that shock or x ray heated gas or H 2 plasma is a minor contributor to the intensities of these fine structure lines. Both the (O 1) (53 microns) and the (Si 2) (35 microns) spectrum show an asymmetric line profile indistinguishable in shape from those of the (O 3) (52 and 88 microns) and (N 3) (57 microns) lines and similar to that of the more extended (C 2) 158 micron line measured previously in M82. We detect two distinct velocity components, which we attribute to emission from two regions at either end of the central bar, where the bar connects to an orbiting torus of neutral gas seen in H 1 and CO J = 1-0. We model separately the two velocity components and derive the physical conditions in these two regions. The clouds in these regions are small, R approx. 1-2 pc, have warm neutral gas surfaces, T approx. 200 K, and are concentrated with volume filling factors of approx. 0.02 and area filling factors of 1-5. The entire central region (R approx. 700 pc) is characterized by a large number, approx. 5 x 10(exp 4), of 2 x 10(exp 3) solar mass clouds with surface densities of approx. 3 x 10(exp 4) cm(exp -3), illuminated by FUV fluxes 10(exp 4) times the average local interstellar value for the Milky Way. These clouds reside in the harsh conditions of a starburst nucleus, with photoevaporation times of 10(exp 6) yr, and collision timescales only about an order of magnitude longer.

Lord, Steven↗

First study of Hohlraum x-ray preheat asymmetry inside an ICF capsule

In indirect drive inertial confinement fusion (ICF), laser induced Hohlraum preheat radiation (so-called M-band, >1.8 keV) asymmetry will lead to asymmetric ablation front and ablator–fuel interface hydrodynamic instability growth in an imploding capsule. First experiments to infer the M-band asymmetries at the capsule were performed on the National Ignition Facility for high density carbon (HDC) ICF capsules in low density fill (0.3 mg/cc 4 He) Au Hohlraums by time resolved imaging of 2.3 keV fluorescence emission of a smaller Mo sphere placed inside the capsule. In this study, the measured Mo emission is pole hot (P 2 > 0) since M-band is generated mainly by the outer laser beams as their irradiance at the Hohlraum wall is 5× higher than for the inner beams. P 2 has a greater negative than positive swing vs time [Δ(P 2 /P 0 )/Δt ~ 0.2/ns], giving insight into laser heated Hohlraum dynamics. P 4 asymmetry is small at the sphere due to efficient geometric smoothing of Hohlraum asymmetries at large Hohlraum -to-capsule ratios. The M-band P 2 history is qualitatively reproduced by radiation hydrodynamic HYDRA simulations. The smaller P 2 than that calculated earlier suggests either less outer beam spot motion and/or preheat emission. At late times, the observed P 2 swing is larger and P 4 is more negative than simulated, which could be due to inner beams being stopped more in the outer beams wall plasma bubble than simulated. Asymmetry at the HDC capsule inner surface (“ice–ablator interface”) is also inferred from the Mo emission asymmetry by an analytic viewfactor model, accounting for the Mo/HDC radius difference and HDC capsule opacity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of Lunar-Inspired Geopolymer Concrete Formulations Mixed and Cured in Microgravity on the International Space Station (ISS)

The research outlined in this presentation investigates the use of various lunar regolith simulants in geopolymer lunar concrete mixes mixed and cured on the International Space Station (ISS). The motivation for this work is to study the effects of gravity on the microstructure of alkali-activated materials cured with heat, and to develop materials for the construction of long-term infrastructure on the lunar surface with in-situ resource utilization (ISRU). ISRU for construction materials reduces the cost and mass of payloads related to lunar construction. The advantage of geopolymer concrete as opposed to traditional portland cement concrete is that water acts as a medium for the polymerization reaction and leaves the system throughout the process, reducing its demand. Twelve samples of lunar regolith simulant and a solution composed of sodium hydroxide and sodium silicate were sent to the ISS. The three simulants were OPRH2N, OPRL2N, and JSC-1AF, using only particles less than 53 µm in diameter to increase reactivity of the simulant. Simulant to solution ratios were determined by workability while mixing. The simulant and solution were sealed Burst Pouches® along with 2 other sealed bags to prevent material from leaking. Crew member F-14 conducted testing on the ISS by introducing the solution to the simulant in the Burst Pouch®, mixing the sample with a spatula, and then clamping the specimen in the fresh state to prevent flow inside the Burst Pouch®. These specimens were then put in a thermos heated to 80C via sealed drinking water bags to cure for 24 hours with a temperature logger. The cured specimens remained in microgravity for at least 28 days and were returned from the ISS in February 2025. The specimens were then brought to the NASA Marshall Space Flight Center (MSFC) to analyze. Material characterization consisted of conducting Micro-CT tests of entire samples in their sealed apparatus to a resolution of 25µm. 2D image slices were saved in each orthogonal direction of each specimen at a 0.03 mm step size from the 3D model to conduct analytical porosity calculations. Representative samples from each specimen were sampled to perform helium gas pycnometery and were then mounted in resin for SEM imaging, EDS, and nanoindentation. Porosity was analyzed analytically using micromechanics modelling with the assistance of the NASA Multiscale Analysis Tool (NASMAT), as well as the NASA Advanced Supercomputing (NAS) servers (V. Saseendran & N. Yamamoto, 2024). Density was measured using helium gas pycnometery and was then compared to the theoretical density for experimental porosity calculation. Due to the samples’ non-uniform shape being cured in a pouch, traditional compression and tensile strength testing could not be performed. Nanoindentation was conducted at Clarkson University to determine the microhardness and reduced modulus of elasticity. Results from flight samples can be compared to ground samples currently in DLR’s possession to determine the effect on microstructure from being mixed and cured in microgravity. This study gives further insight and understanding of geopolymer lunar concrete and its viability as a lunar construction material with ISRU.

Adam Johnson↗

Infrared coronal emission lines and the possibility of their maser emission in Seyfert nuclei

Energetic emitting regions have traditionally been studied via x-ray, UV and optical emission lines of highly ionized intermediate mass elements. Such lines are often referred to as 'coronal lines' since the ions, when produced by collisional ionization, reach maximum abundance at electron temperatures of approx. 10(exp 5) - 10(exp 6) K typical of the sun's upper atmosphere. However, optical and UV coronal lines are also observed in a wide variety of Galactic and extragalactic sources including the Galactic interstellar medium, nova shells, supernova remnants, galaxies and QSOs. Infrared coronal lines are providing a new window for observation of energetic emitting regions in heavily dust obscured sources such as infrared bright merging galaxies and Seyfert nuclei and new opportunities for model constraints on physical conditions in these sources. Unlike their UV and optical counterparts, infrared coronal lines can be primary coolants of collisionally ionized plasmas with 10(exp 4) less than T(sub e)(K) less than 10(exp 6) which produce little or no optical or shorter wavelength coronal line emission. In addition, they provide a means to probe heavily dust obscured emitting regions which are often inaccessible to optical or UV line studies. In this poster, we provide results from new model calculations to support upcoming Infrared Space Observatory (ISO) and current ground-based observing programs involving infrared coronal emission lines in AGN. We present a complete list of infrared (lambda greater than 1 micron) lines due to transitions within the ground configurations 2s(2)2p(k) and 3s(2)3p(k) (k = 1 to 5) or the first excited configurations 2s2p and 3s3p of highly ionized (x greater than or equal to 100 eV) astrophysically abundant (n(X)/n(H) greater than or equal to 10(exp -6)) elements. Included are approximately 74 lines in ions of O, Ne, Na, Mg, Al, Si, S, Ar, Ca, Fe, and Ni spanning a wavelength range of approximately 1 - 280 microns. We present new results from detailed balance calculations, new critical densities for collisional de-excitation, intrinsic photon rates, branching ratios, and excitation temperatures for the majority of the compiled transitions. The temperature and density parameter space for dominant cooling via infrared coronal lines is presented, and the relationship of infrared to optical coronal lines is discussed.

Greenhouse, Matthew A.↗

Preliminary modeling of triply periodic minimal surface (TPMS) structures using RELAP5-3D

With the United States Department of Energy (DOE)’s goal of quadrupling the nation’s nuclear energy supply by 2050, and with the Advanced Fuels Campaign pushing for new types of advanced reactor fuels and geometries, the need has arisen for new nuclear fuel designs. One such design is to swap out current nuclear fuel geometries in exchange for another type of geometry, called a Triply Periodic Minimal Surface (TPMS). TPMSs are self-supporting, infinitely repeating lattices—attributes that lend themselves well to additive manufacturing. These surfaces also possess enhanced heat transfer properties thanks to their internal area changes and large surface-area-to-volume ratios. Their drawback, however, is an increased pressure drop. Given the small amount of correlations and data (Reynolds numbers in the 2,000–8,000 range), and the minimal amount of experience so far obtained by modeling TPMS structures using 1D systems codes such as the Reactor Excursion and Leak Analysis Program (RELAP5-3D), further research into this topic was needed. Using data from the University of Wisconsin - Madison (UW), curve fits were created for both a Heat Transfer Coefficient (HTC) correlation and a Darcy friction factor empirical coefficient correlation. The curves’ coefficients and multipliers were then output and utilized in RELAP5-3D models of two upcoming experiments—Flow Loop for INFLUX Pressure drop (FLIP) and Microreactor Agile Non-nuclear Experimental Test (MAGNET)—aimed at increasing the available data for Reynolds numbers to the 16,000–36,000 range for TPMS structures. The models were run under the conditions utilized by a Computational Fluid Dynamics (CFD) analysis performed by another group at Idaho National Laboratory. Only CFD pressure drop values were obtained from the FLIP test, and those values showed that the RELAP5-3D models had a lower rate of pressure increase in comparison to the CFD values. In addition, there seemed to be a vertical shift upward in the pressure drop for both models whenever the TPMS porosity decreased, and the RELAP5-3D models showed a higher vertical shift in comparison to the CFD values. The MAGNET results did not correspond to any CFD or experimental results against which they could be compared, so they were instead compared against the proposed CFD input conditions. These values were then compared with each other to make sure the model seemed to be performing as expected, paving the way for future tests that can be run for the purpose of further analyses and comparisons. The pressure drop increased with temperature and mass flow rate independently. The temperature change would decrease with increasing mass flow rate and temperature, which was just as we expected based on the fact that the lower viscosity and decreased density would result in higher friction and churning losses. The last metric that was assessed was the enthalpy flow change, which increased with increasing mass flow rate and decreasing temperature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Ocean Turbulence, III: New GISS Vertical Mixing Scheme

We have found a new way to express the solutions of the RSM (Reynolds Stress Model) equations that allows us to present the turbulent diffusivities for heat, salt and momentum in a way that is considerably simpler and thus easier to implement than in previous work. The RSM provides the dimensionless mixing efficiencies Gamma-alpha (alpha stands for heat, salt and momentum). However, to compute the diffusivities, one needs additional information, specifically, the dissipation Epsilon. Since a dynamic equation for the latter that includes the physical processes relevant to the ocean is still not available, one must resort to different sources of information outside the RSM to obtain a complete Mixing Scheme usable in OGCMs. As for the RSM results, we show that the Gamma-alpha s are functions of both Ri and Rq (Richardson number and density ratio representing double diffusion, DD); the Gamma-alpha are different for heat, salt and momentum; in the case of heat, the traditional value Gamma-h = 0.2 is valid only in the presence of strong shear (when DD is inoperative) while when shear subsides, NATRE data show that Gamma-h can be three times as large, a result that we reproduce. The salt Gamma-s is given in terms of Gamma-h. The momentum Gamma-m has thus far been guessed with different prescriptions while the RSM provides a well defined expression for Gamma-m(Ri,R-rho). Having tested Gamma-h, we then test the momentum Gamma-m by showing that the turbulent Prandtl number Gamma-m/Gamma-h vs. Ri reproduces the available data quite well. As for the dissipation epsilon, we use different representations, one for the mixed layer (ML), one for the thermocline and one for the ocean;s bottom. For the ML, we adopt a procedure analogous to the one successfully used in PB (planetary boundary layer) studies; for the thermocline, we employ an expression for the variable epsilon/N(exp 2) from studies of the internal gravity waves spectra which includes a latitude dependence; for the ocean bottom, we adopt the enhanced bottom diffusivity expression used by previous authors but with a state of the art internal tidal energy formulation and replace the fixed Gamma-alpha = 0.2 with the RSM result that brings into the problem the Ri, R-rho dependence of the Gamma-alpha; the unresolved bottom drag, which has thus far been either ignored or modeled with heuristic relations, is modeled using a formalism we previously developed and tested in PBL studies. We carried out several tests without an OGCM. Prandtl and flux Richardson numbers vs. Ri. The RSM model reproduces both types of data satisfactorily. DD and Mixing efficiency Gamma-h(Ri,Rq). The RSM model reproduces well the NATRE data. Bimodal epsilon-distribution. NATRE data show that epsilon (Ri < 1) approximately equals 10epsilon(Ri > 1), which our model reproduces. Heat to salt flux ratio. In the Ri much greater than 1 regime, the RSM predictions reproduce the data satisfactorily. NATRE mass diffusivity. The z-profile of the mass diffusivity reproduces well the measurements at NATRE. The local form of the mixing scheme is algebraic with one cubic equation to solve.

Canuto, V. M.↗

Low-temperature etching of silicon oxide and silicon nitride with hydrogen fluoride

Etching of high aspect ratio features into alternating SiO 2 and SiN layers is an enabling technology for the manufacturing of 3D NAND flash memories. In this paper, we study a low-temperature or cryo plasma etch process, which utilizes HF gas together with other gas additives. Compared with a low-temperature process that uses separate fluorine and hydrogen gases, the etching rate of the SiO 2 /SiN stack doubles. Both materials etch faster with this so-called second generation cryo etch process. Pure HF plasma enhances the SiN etching rate, while SiO 2 requires an additional fluorine source such as PF 3 to etch meaningfully. Further, the insertion of H 2 O plasma steps into the second generation cryo etch process boosts the SiN etching rate by a factor of 2.4, while SiO 2 etches only 1.3 times faster. We observe a rate enhancing effect of H 2 O coadsorption in thermal etching experiments of SiN with HF. Ammonium fluorosilicate (AFS) plays a salient role in etching of SiN with HF with and without plasma. AFS appears weakened in the presence of H 2 O. Density functional theory calculations confirm the reduction of the bonding energy when NH 4 F in AFS is replaced by H 2 O.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spectral and Temporal Properties of Galactic Black Hole Systems

Kusunose, Mineshige & Yamada (1996; hereafter KMY) extended the model of Kusunose & Mineshige (1995) to the Galactic black hole candidates by considering nonthermal electron injection with gamma(EQ\0(,\s\up2(less than),\s\do-l(_))) 10. The effects of pair escape and advection on the disk structure and general relativistic effects on the emission spectrum were also examined. They found that the energy spectral index (alpha)(sub x) of the power law X-rays is about-0.8 and-2.0 when 1(sub soft)/1 = 0.2 and 2, respectively, where 1(sub soft)/1 is the ratio of the compactness of the injected soft photons to that of the gravitational energy. The power law index was found to be nearly independent of the mass accretion which is consistent with the observed luminosity independence. The model with small 1(sub soft)/1 (less than 1) shows promise for explaining the low state observed in Galactic black hole candidates. Model fits were provided for GX339-4 and Cyg X- 1 data from COMPTEL and OSSE on the Compton Gamma Ray Observatory. The difference in emission spectra between thermal disks and the model of KMY appears only in the energy range greater than 100 keV. Li, Kusunose and Liang (1996) studied stochastic particle acceleration to produce nonthermal particle distributions which then were used in the model of Kusunose & Mineshige (1995) to model the spectrum above 1 Mev from GBHC's. Under certain conditions, stochastic electron acceleration overcame Coulomb and Compton losses resulting in a suprathermal electron population. Good fits were obtained by COMPTEL and OSSE observations of Cyg X-1 and GRO J0422+23. Kusunose & Mineshige (1996a) examined the role of electron-positron pairs in advection-dominated disks. They found that the results for advection-dominated disks without pairs are not qualitatively changed by including pairs. Summaries of work sponsored by this grant are given in Wheeler, Kim, Moscoso, Kusunose & Mineshige (1996) and Kusunose (1996) Work was also done on developing a model for an e(+-) pair wind from the inner disk region of a black hole. The model consists of three zones: a pair production/annihilation zone at the base of the wind, a pair annihilation zone slightly further out from the disk and a wind zone extending to infinity where no annihilation occurs. The model assumes an input X-ray / gamma-ray power-law spectrum as a function of photon energy and angular distribution. Pairs are created in the pair production/annihilation zone via photon-photon, photonparticle and particle-particle collisions. The bulk velocity of the pairs is obtained from the momentum component perpendicular to the disk taking into account the radiation pressure acceleration of the pairs. Energy balance in the pair production/annihilation zone is calculated by assuming that the momentum component parallel to the disk is thermalized and by taking into account heating/cooling via anisotropic thermal Compton scattering along with the cooling processes of bremsstrahlung and synchrotron radiation. The pair density is also calculated. The outer boundary of the pair production/annihilation zone is where the optical depth for photon-photon pair produciton is unity. Beyond this point only pair annihilation will occur in the pair annihilation zone. When the pair density becomes small, pair annihilation will be negligible and the pairs will flow freely to infinity. In this model we have found that the X-ray power-law spectral index, alpha(sub x) is the primary parameter which determines the density, temperature and velocity of the pair production/annihilation zone (Moscoso, Kusunose & Wheeler 1996).

Wheeler, J. Craig↗

A Warm Garage for a Lunar Rover

Approach: One approach to heating a rover during the lunar night is the so-called thermal wadis concept [1]. This involves heating the regolith with solar concentrators and placing the rover on the heated surface for the night. Since the regolith is heated by a relatively weak heat flux, a high thermal conductivity is required for heating a sufficiently large mass of regolith. However, lunar regolith has a low thermal conductivity. Therefore, the concept involves increasing the conductivity by sintering the regolith, which requires a significant energy input and complex procedures. Here we propose an alternative approach where the low thermal conductivity of regolith is an advantage. Specifically, we propose to use a highly exothermic combustible mixture for heat generation. The mixture pellets are placed in the surface layer of regolith and ignited. The combustion forms condensed products and releases heat, which then slowly spreads to the surrounding regolith. Heat can also be transferred, for example, by heat pipes, into radiant heating surfaces installed on the ground. A greenhouse that transmits sunlight during the day and decreases the radiative heat losses during the night can also be installed. Selection of the Heat-generating Mixture: The reactive mixture should have a high specific energy and generate only condensed products since gases could disturb the regolith layer, carry enthalpy out of the system, and lead to an explosion. There are mixtures, (sometimes called pyrolants) that possess very high specific energies. One example is magnesium-Teflon-Viton mixtures used in flares. However, they produce gases and may cause explosions. Other mixtures that include magnesium cannot be used either because of the high vapor pressure of Mg at temperatures well below the combustion temperature. Recently, mixtures that involve lithium peroxide (Li2O2) have been proposedfor using in space power systems [2]. However, they produce lithium oxide (Li2O), which boils at 2800 K at 1 atm and hence at a lower temperature in vacuum. Fortunately, there exist many mixtures that release a lot of heat and form only condensed products during the combustion. Many such mixtures have been used for self-propagating high-temperature synthesis (SHS) of various materials [3, 4]. For the application discussed here, t itanium/boron (1:2 mole ratio) mixture appears to be particularly promising. The specific energy is 4.0 MJ/kg (1.1 kWh/kg), the adiabatic flame temperature is about 3200 K, and the reaction forms solid titanium diboride (TiB2, melting point: 3500 K). The mixture can be ignited easily with a heated tungsten wire, and it has been used widely as a booster to ignite the main mixture in the SHS process.Estimates: Assuming that specific heat of regolith is 500 J/(kg∙K) [5] and all generated heat is transferred to regolith, 12.5 kg of the Ti/B mixture would increase the temperature of 1000 kg of regolith by 100 K. To evaluate the rate of heat transfer in the regolith, a spherical model was analyzed where the heat released by a 12.5 kg Ti/B core propagates by thermal conduction through a 1000 kg regolith shell with no heat loss from its outer surface. At a bulk density of 1500 kg/m3 [5], the radius of the shell was 54 cm, while the radius of the core was about 11 cm. The calculations were conducted using Thermal Desktop SINDA/FLUINT (Cullimore and Ring Technologies) software at two constant values of bulk thermal conductivityk of the regolith: 0.001 and 0.01 W/(m∙K). The results show that after 14.5 days the core lost 31% of the released heat at the lower k and 77% at the higher k. At a distance of 20 cm from the core surface, the temperature of the regolith increased by only 1 K at the lower k and by 132 K at the higher k. In reality, the regolith near the heat source will be melted, so its thermal conductivity will increase significantly. Nevertheless, the conducted estimates indicate that combustion-based heat generators, placed directly in the regolith, could provide heat during a rather long period such as the lunar night.Conclusion: Heat generators based on gasless combustion of highly energetic reactive mixtures could be installed directly in the surface layer of lunar regolith. Because of the low thermal conductivity of the regolith, such generators would keep thermal energy for days and gradually supply heat to a rover/lander.Acknowledgment: The material presented in this work is based upon the work supported by National Aeronautics and Space Administration (NASA) under Grant #80NSSC20K0293.References: [1] Balasubramaniam R. et al. (2011) J. Thermophys. Heat Trans., 25,130−139. [2] Blair R.G. and Vasu S.S. (2022) Conf. Advanced Power Systems for Deep Space Exploration. [3] Varma A. et al. (1998) Adv. Chem. Eng., 24,79−226. [4] Levashov E.A. et al. (2017) Int. Mater. Rev., 62,203−239. [5] Wood-Robinson R. et al. (2019) J. Geophys. Res. Planets, 124, 1989−2011.

lunar↗

Performance Testing of an Integrated Magnetic Power Take-Off

A wave energy converter (WEC) and the power take-off (PTO) generator system can be represented by using a mass-spring-dampener model. By incorporating a negative stiffness spring within the PTO, the overall stiffness of the PTO can be lowered allowing the impedance of the PTO to be more closely matched with the WEC. By designing for impedance matching the WEC can greatly enhance its power generation capability. This project has involved the design, fabrication and testing of a new type of linear-stroke length and rotary stroke length adjustable negative stiffness magnetic springs for use within a wave energy converter (WEC). The magnetic springs were studied by using 3-D finite element analysis with the objective of creating high energy density and a long linear stroke length. After the construction and testing of both a proof-of-principle adjustable linear and rotary magnetic spring prototype the rotary (torsional) magnetic spring was selected for scaling up analysis. The selected scaled-up proof-of-principle magnetic spring had a peak torque of 850 Nm with a ±45 degree stroke length. By translating the inner rotor relative to the outer rotor, the stiffness could be adjusted to be either negative or positive stiffness. During this project a magnetic lead screw was also studied, and it was shown that by combining the magnetic lead screw with a linearly translating magnetic spring a very long rotary stroke length could be attained. However, as the magnetic lead screw increased complexity and reduced overall energy density relative to a rotary (torsional) spring, this design approach was not further pursed. Dynamometer testing of the variable stiffness magnetic springs was first completed by Portland State University and following this the magnetic spring performance was independently verified by Sandia National Laboratory (Sandia). A WEC analysis when using a variable stiffness magnetic spring was completed by using the WecOptTool. WecOptTool is an open-source WEC optimization software developed by Sandia that supports efficient power take-off (PTO) and control optimization. Wave condition data from the Oregon PacWave test-site was used in this study. The analysis showed that a WEC with a tunable stiffness value could consistently achieve about 80% of its maximum theoretical power production. The use of a tunable stiffness WEC, rather than zero-stiffness or a constant stiffness WEC was also shown to lead to a smaller maximum PTO force requirement. The variable stiffness magnetic spring was integrated into an experimental WEC developed by Sandia, called a Wave-Bot. Sandia successfully completed water-tank testing of the Wave-Bot at the Navy’s Carderock, Maryland, wave-basin test site. The wave-basin testing helped to experimentally demonstrate the operating capabilities and increased power generation capability of a WEC containing an integrated variable stiffness magnetic spring. A WEC capacity factor analysis was also completed. The capacity factor was defined as the ratio of annual average WEC generator power to the maximum (RMS) generated power. This capacity factor provided a means of identifying the ratio of potential revenue to cost. It was calculated that if a tunable variable stiffness magnetic spring has RMS constraints it could operate with a capacity factor above 30%. This is comparable to a wind turbine’s capacity factor.

16 TIDAL AND WAVE POWER↗

Impact of mid- Z gas fill on dynamics and performance of shock-driven implosions at the OMEGA laser

Shock-driven implosions with 100% deuterium (D 2 ) gas fill compared to implosions with 50:50 nitrogen-deuterium (N 2 ⁢D 2 ) gas fill have been performed at the OMEGA laser facility to test the impact of the added mid-Ζ fill gas on implosion performance. Ion temperature (Τ ion ) as inferred from the width of measured DD-neutron spectra is seen to be 34%±6% higher for the N 2⁢ D 2 implosions than for the D 2 -only case, while the DD-neutron yield from the D 2 -only implosion is 7.2±0.5 times higher than from the N 2⁢ D 2 gas fill. The T ion enhancement for N 2 ⁢D 2 is observed in spite of the higher Z, which might be expected to lead to higher radiative loss, and higher shock strength for the D 2 -only versus N 2 ⁢D 2 implosions due to lower mass, and is understood in terms of increased shock heating of N compared to D, heat transfer from N to D prior to burn, and limited amount of ion-electron-equilibration-mediated additional radiative loss due to the added higher-Z material. Further, this picture is supported by interspecies equilibration timescales for these implosions, constrained by experimental observables. The one-dimensional (1D) kinetic Vlasov-Fokker-Planck code ifp and the radiation hydrodynamic simulation codes hyades (1D) and xrage [1D, two-dimensional (2D)] are brought to bear to understand the observed yield ratio. Comparing measurements and simulations, the yield loss in the N 2 ⁢D 2 implosions relative to the pure D 2 -fill implosion is determined to result from the reduced amount of D 2 in the fill (fourfold effect on yield) combined with a lower fraction of the D 2 fuel being hot enough to burn in the N 2 ⁢D 2 case. The experimental yield and T ion ratio observations are relatively well matched by the kinetic simulations, which suggest interspecies diffusion is responsible for the lower fraction of hot D 2 in the N 2 ⁢D 2 relative to the D 2 -only case. The simulated absolute yields are higher than measured; a comparison of 1D versus 2D XRAGE simulations suggest that this can be explained by dimensional effects. The hydrodynamic simulations suggest that radiative losses primarily impact the implosion edges, with ion-electron equilibration times being too long in the implosion cores. The observations of increased T ion and limited additional yield loss (on top of the fourfold expected from the difference in D content) for the N 2 ⁢D 2 versus D 2 -only fill suggest it is feasible to develop the platform for studying CNO-cycle-relevant nuclear reactions in a plasma environment.

47 OTHER INSTRUMENTATION↗

Achieving record hot spot energies with large HDC implosions on NIF in HYBRID-E

HYBRID-E is an inertial confinement fusion implosion design that increases energy coupled to the hot spot by increasing the capsule scale in cylindrical hohlraums while operating within the current experimental limits of the National Ignition Facility. HYBRID-E reduces the hohlraum scale at a fixed capsule size compared to previous HYBRID designs, thereby increasing the hohlraum efficiency and energy coupled to the capsule, and uses the cross-beam energy transfer (CBET) to control the implosion symmetry by operating the inner (23° and 30°) and outer (44° and 50°) laser beams at different wavelengths (Δλ> 0). Small case to capsule ratio designs can suffer from insufficient drive at the waist of the hohlraum. We show that only a small amount of wavelength separation between the inner and outer beams (Δλ1-2 Å) is required to control the symmetry in low-gas-filled hohlraums (0.3 mg/cm3 He) with enough drive at the waist of the hohlraum to symmetrically drive capsules 1180 μm in outer radius. This campaign is the first to use the CBET to control the symmetry in 0.3 mg/cm3 He-filled hohlraums, the lowest gas fill density yet fielded with Δλ> 0. We find a stronger sensitivity of hot spot P2 in μm per Angstrom (40–50 μm/Å wavelength separation) than observed in high-gas-filled hohlraums and previous longer pulse designs that used a hohlraum gas fill density of 0.6 mg/cm3. There is currently no indication of transfer roll-off with increasing Δλ, indicating that even longer pulses or larger capsules could be driven using the CBET in cylindrical hohlraums. We show that the radiation flux symmetry is well controlled during the foot of the pulse, and that the entire implosion can be tuned symmetrically in the presence of the CBET in this system, with low levels of laser backscatter out of the hohlraum and low levels of hot electron production from intense laser–plasma interactions. Radiation hydrodynamic simulations can accurately represent the early shock symmetry and be used as a design tool, but cannot predict the late-time radiation flux symmetry during the peak of the pulse, and semi-empirical models are used to design the experiments. Deuterium–tritium (DT)-layered tests of 1100 μm inner radius implosions showed performance close to expectations from simulations at velocities up to ~360 km/s, and record yields at this velocity, when increasing the DT fuel layer thickness to mitigate hydrodynamic mixing of the ablator into the hot spot as a result of defects in the ablator. However, when the implosion velocity was increased, mixing due to these defects impacted performance. The ratio of measured to simulated yield for these experiments was directly correlated with the level of observed mixing. These simulations suggest that reducing the mixing, e.g., by improving the capsule defects, could result in higher performance. In addition, future experiments are planned to reduce the coast time at this scale, delay between the peak compression and the end of the laser, to increase the hot spot convergence and pressure. To reduce the coast time by several hundred ps compared to the 1100 μm inner radius implosions, HYBRID-E has also fielded 1050 μm inner radius capsules, which resulted in higher hot spot pressure and a fusion energy yield of ~170 kJ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗