Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BUBBLE”

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 811 records · Page 45

High-speed imaging of transition from fluid breakup to phase explosion in electric explosion of tungsten wires in air

High-speed visible imaging of sub-microsecond electric explosion of wires at the low specific energy deposition threshold reveals three distinct modes of wire failure as capacitor charge voltage and energy deposition are increased. For 100 micron diameter gold-plated tungsten wires of 2 cm length, deposited energies of 1.9 eV/atom produces a liquid column that undergoes hydrodynamic breakup into droplets with radii of order of wire diameter on timescales of 200 μs. Instability growth, column breakup, and droplet coalescence follow classic Rayleigh-Plateau predictions for instability of viscous fluid column. Above 3.2 eV/atom of deposited energy, wires are seen to abruptly transition to an expanding mixture of micron scale liquid-droplets and vapor within one frame (less than 3.33 μs), which has been termed ‘phase explosion’ in literature. Between these two limits, at 2.5 eV/atom of deposited energy, wire radius is unchanged for the first 10 μs before the onset of a rapid expansion and disintegration that resembles homogenous nucleation of mechanically unstable bubbles. Here, thermodynamic calculations are presented that separate cases by temperature obtained during heating: below boiling point, near boiling point, and exceeding boiling point.

42 ENGINEERING↗

Non-dilute helium-related defect interactions in the near-surface region of plasma-exposed tungsten

We report a systematic energetic analysis of helium-related defect interactions that mediate helium (He) segregation on surfaces of plasma-exposed tungsten at different levels of He ion implantation. We focus on high He fluences that increase the He content in the plasma-exposed material well beyond the dilute limit of He concentration and employ atomic configurations generated by large-scale molecular dynamics simulations of He-implanted tungsten. We perform systematic molecular statics computations of cluster– defect interaction energetics in the highly defect-rich near-surface region of plasma-exposed tungsten for small mobile helium clusters as a function of the cluster distance from the surface. In this region, mobile clusters are also subjected to the stress fields generated by defects such as helium bubbles and other clusters, which govern cluster–defect interactions in addition to the cluster–surface interaction. Based on systematic investigation of individual cluster–defect interactions, we develop a mathematical framework to describe the interaction energy landscapes consisting of elastic interaction potential perturbations to finite-width square-well potentials, where the potential well accounts for cluster trapping by the defect at close range and subsequent coalescence and the perturbation potential is parameterized according to elastic inclusion theory. Superposition of all the relevant interaction potentials provides a comprehensive description of the interaction energy landscape that would be experienced by a small mobile cluster along its migration path toward the plasma-exposed surface at high He fluence. Such descriptions are particularly important for developing atomistically-informed, hierarchical multi-scale models of helium cluster dynamics in plasma-facing materials.

36 MATERIALS SCIENCE↗

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↗

Structural organization of space polymers

Extraterrestrial polymers of glycine with iron have been characterized by mass spectrometry to have a core mass of 1494 Da with dominant rod-like variants at mass to charge ratios of 1567 and 1639 [McGeoch et al., “Meteoritic proteins with glycine, iron and lithium,” arXiv:2102.10700 (2021)]. Several principal macro-structural morphologies are observed in solvent extracts from a chondritic Vigarano class alteration type 3 meteoritic material. The first is an extended sheet of linked (three-legged) triskelia containing the 1494 Da core entity that encloses gas bubbles in the solvent. A second is of fiber-like crystals found here, via x-ray diffraction, to be multiple-walled nanotubes made from a square lattice of the 1494 Da polymer. A third is a dispersion of floating phantom-like short tubes of up to 100 μm length with characteristic angled bends that suggest the influence of a specific underlying protein structure. Here, it is proposed that the angled tubes are the observable result of a space-filling superpolymerization of 1638 Da polymer subunits guided by the tetragonal symmetry of linking silicon bonds. Distorted hexagonal sheets are linked by perpendicular subunits in a three-dimensional hexagonal diamond structure to fill the largest possible volume. This extended very low-density structure is conjectured to have dominated in a process of chemical selection because it captured a maximum amount of molecular raw material in the ultra-low density of molecular clouds or of the proto-solar nebula. This could have led ultimately to the accretion of the earliest planetary bodies.

59 BASIC BIOLOGICAL SCIENCES↗

Mechanisms of shape transfer and preheating in indirect-drive double shell collisions

Implosions of Hohlraum-driven double shell targets as an alternative inertial confinement fusion concept are underway at the National Ignition Facility. The double shell system relies on a series of energy transfer processes starting from thermal x-ray absorption by the outer shell, followed by collisional transfer of kinetic energy to a heavy metal inner shell, and finally, conversion to the internal energy of the deuterium-tritium fuel. During each of these energy transfer stages, low-mode asymmetries can act to reduce the ideal transfer efficiency degrading double shell performance. Mechanisms, such as hard x-ray preheat from the Hohlraum, not only decrease the efficiency of kinetic energy transfer but may also be a source of low-mode asymmetry. In this article, we evaluate the shape transfer processes through the time of shell collision using two-dimensional integrated Hohlraum and capsule computations. We find that the dominant mode of the shape transfer is well described using a “radial impulse” model from the shape of the foam pressure reservoir. To evaluate the importance of preheat on inner shell shape, we also report on first measurements of Au L-shell preheat asymmetry in a double shell with a tungsten pusher. These measurements showed a 65% higher preheat velocity at the pole of the capsule relative to the equator. We also found that the experiments provided rigorous constraints by which to test the Hohlraum model settings that impact the amount and symmetry of Au L-shell preheat via the plasma conditions inside the outer cone Au bubble.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Local measurement of bulk thermal diffusivity using photothermal radiometry

Here, an experimental methodology using photothermal radiometry is developed for the accurate measurement of bulk thermal diffusivity of nuclear fuels and materials irradiated to high doses. Under these conditions, nuclear fuels, such as uranium oxide, and moderator materials, such as graphite, become friable, which requires characterization techniques that can accommodate irregularly shaped fragments. Photothermal radiometry, a good candidate for this application, involves locally heating a sample by using a laser and measuring the temperature field by monitoring blackbody radiation. The interaction volume for this study, less than a millimeter, is carefully chosen to sample a statistically significant number of large-scale structural features, such as pores and gas filled bubbles, and is small enough that the sample fragments can be treated as a thermal half-space. The thermal diffusivity standards considered in this study cover a range of thermal diffusivities representative of both fresh and spent nuclear fuels. We also consider a sample having a porous microstructure representative of large-scale structures found in materials irradiated to high doses. Our measurement methodology circumvents complex thermal wave models that address optical diffraction, nonlinear transfer function associated with blackbody radiation, and finite sample size effects. Consequently, the large measurement uncertainty associated with modeling these effects can be avoided. While the emphasis here is on nuclear fuels and materials, this measurement approach is well suited to measure thermal transport in a variety of technologically important materials associated with advanced synthesis techniques. Examples range from small, exotic single crystals grown using hydrothermal growth techniques to additively manufactured components having complex geometries.

36 MATERIALS SCIENCE↗

Constraining the release of Sn to the ambient melting point following shock loading using time-resolved x-ray diffraction

The formation of liquid following release from a shocked state governs the transition from spall to cavitation and the formation of ejecta in metals. In order to build physics-based models of these processes, it is necessary to critically evaluate the relative importance of kinetics and entropy generation during the release along with the accuracy of multiphase equations of state. Tin (Sn) has served as a testbed for a variety of experiments examining strength and ejecta due to its accessible melt boundary and solid–solid phase transitions. This work presents experiments examining the phase evolution of high purity Sn following the shock and release to ambient stress near the melting point. Sn is found to release to states between its ambient solidus and liquidus from approximately 19 to 33 GPa under uniaxial loading, with the two-phase region being characterized by a reduction in the intensity of the (220), (211) β -Sn doublet. Jetting experiments performed at 27–28 GPa exhibit comparable diffraction patterns with what is observed following the uniaxial release. The solid fractions of β-Sn in the ambient mixed phase region are found to decrease linearly with increasing shock stress as increasing liquid Sn is formed. Here, the results provide much needed information for interpreting measurements of dynamic strength at a high strain rate and experiments examining cavitation and shallow bubble collapse in Sn.

, spall↗

Mixture Model for Refrigerant Pairs R-32/1234yf, R-32/1234ze(E), R-1234ze(E)/227ea, R-1234yf/152a, and R-125/1234yf

In this work, thermodynamic models based on the corresponding states framework with departure terms are developed for the refrigerant pairs R-32/1234yf, R-32/1234ze(E), R-1234ze(E)/227ea, R-1234yf/152a, and R-125/1234yf. These models are based on new measurements of density, speed of sound, and phase equilibria, combined with the data available in the literature. The model for R-32/1234yf is most comprehensive in its data coverage, with speed of sound deviations within 1%, density deviations within 0.1%, and bubble- and dew-point pressure deviations within 1%. In conclusion, the other mixtures have generally more limited data availability but a similar goodness of fit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theory and simulations of linear and nonlinear two-dimensional Rayleigh–Taylor dynamics with variable acceleration

Interfacial Rayleigh–Taylor mixing is crucial to describing important natural and engineering processes, such as exploding supernovae, laser micromachining, hot spots in inertial confinement fusion, and optical telecommunications. These require the characterization of the time dependence of the driving acceleration. We compare our theoretical formulation based on group theory foundations with interface-capturing numerical simulations for linear and nonlinear two-dimensional Rayleigh–Taylor instabilities in a finite-sized domain with time-varying acceleration over broad ranges of Atwood numbers and acceleration exponents. Detailed corroboration between theory and simulations is provided for this foundational case. Both demonstrate the strong interfacial nature of Rayleigh–Taylor instabilities, which suggests that practical flow fields can be reconstructed from the derived fluid potential using the proposed theory. A robust agreement is also obtained for the early and late-time evolution of the amplitudes of the bubble and spike, which demonstrate that the Rayleigh–Taylor flow can transition to the mixing regime even for a single-mode initial perturbation. Corroboration with experiments of high energy density plasmas motivated by studies of supernovae is also achieved. In addition, a long-standing puzzle in Rayleigh–Taylor dynamics on the interplay between the acceleration, the shear, and the interface morphology in the theory and simulations is resolved by accounting for finite viscosity of the fluids. The characterization of Rayleigh–Taylor instabilities as a highly interfacial phenomenon provides valuable insight into its multiscale nature, which enhances the design and understanding of numerous processes of practical interest.

42 ENGINEERING↗

A Deep Potential model for liquid–vapor equilibrium and cavitation rates of water

Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology—a machine learning approach—to study this ubiquitous phase transition, starting from the phase diagram in the liquid–vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. Further, we find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid–vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid–vapor coexistence and water cavitation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Compressible Rayleigh–Taylor instability subject to isochoric initial background state

The effect of compressibility on the single-mode Rayleigh–Taylor instability is examined using two (2D) and three-dimensional (3D) direct numerical simulations. To isolate compressibility from background stratification effects, this work employs a constant density profile on each side of the interface. The numerical simulations are performed at various Reynolds numbers using the gas kinetic method for static Mach numbers up to M = 0.4. The most important finding is that compressibility acting in isolation enhances the instability and perturbations grows faster with increasing Mach number, unlike previous results with background isothermal state, which show suppression of the instability at higher static Mach numbers. In addition, compressibility is also shown to increase the bubble-spike asymmetry. While the instability grows faster for the 3D case, the findings are qualitatively similar in 2D and 3D. Furthermore, the dynamical reasons underlying the effect of compressibility are elucidated by examining the evolution of vorticity and turbulent kinetic energy transport equations.

42 ENGINEERING↗

Simulations and models for the Richtmyer–Meshkov instability with broadband perturbations

The Richtmyer–Meshkov instability (RMI) is shock driven and affects many phenomena from inertial fusion to supernova explosions. The behavior of single-modes in the RMI has been studied extensively but less is known with the broadband perturbations that occur in applications. Here, we describe extensive numerical simulations and modeling of the RMI with broadband perturbations with an initial power spectrum of the form P ∝ km, where k is the wavenumber. The hydrodynamic condition is the same as in the θ-Group Collaboration [Thornber et al., “Late-time growth rate, mixing, and anisotropy in the multimode narrowband Richtmyer–Meshkov instability: The θ-group collaboration,” Phys. Fluids 29, 105107 (2017)] with a Mach ∼1.86 shock and gamma-law = 5/3 fluids with Atwood number A = −0.5. The bubble amplitude hB is found to grow in two stages. Initially, hB undergoes a phase-inversion and grows linearly in time (t) at a rate consistent with a linear theory. Asymptotically, hB grows as a power law ∼tθ when k|hB| > O(1), where k is the average wavenumber for the initial spectrum. The RMI behavior in simulations and modeling agree over a wide range of exponent m, spectral width, initial amplitude, and time. The agreement is quantified objectively using statistical analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

Efficient laser-driven proton acceleration from a petawatt contrast-enhanced second harmonic mixed-glass laser system

Efficient laser-driven plasma acceleration of ion beams requires precision control of the target–plasma profile, which is crucial to optimize the laser energy transfer. Along the laser propagation direction, this can be achieved by tailoring the temporal structure of the laser pulse. We show for the first time that frequency-doubling of a short pulse (hundreds-femtosecond range) petawatt-class mixed-glass laser system, which results in temporal intensity contrast enhancement, enables surface and volumetric laser–energy coupling, and the acceleration of proton beams from few-nanometer-thick foil targets. Experimentally, maximum ion energies and laser-to-proton energy conversion efficiencies were found to be both maximized at optimum laser and target conditions manifested when the normalized target density nearly equalizes the normalized laser vector potential, which is in agreement with theory and simulations. These signatures are recognized as a unique indication of the interaction between ultra-intense laser pulses with high temporal intensity contrast and ultra-thin nanometer-scale targets. Transverse modulations of accelerated proton beams in the form of bubble- and ring-like structures measured in the thinnest targets provide additional evidence of volumetric laser-driven particle acceleration regimes and transitional features in ultra-thin foil targets specific to laser–plasma interactions characterized by a high temporal intensity contrast. These results open avenues in the generation of high contrast laser pulses from short-pulse-femtosecond petawatt mixed-glass laser systems and demonstrate the feasibility of this technique for applications requiring high laser intensity contrast with high efficiency.

Physics↗

Measuring the density, viscosity, and surface tension of molten titanates using electrostatic levitation in microgravity

Rare earth and barium titanates are useful as ferroelectric, dielectric, and optical materials. Measurements of their thermophysical properties in the liquid state can help guide melt processing technologies for their manufacture and advance understanding of fragile liquids' behavior and glass formation. Here, we report the density, thermal expansion, viscosity, and surface tension of molten BaTi 2 O 5 , BaTi 4 O 9 , and 83TiO 2 -17RE 2 O 3 (RE = La or Nd). Measurements were made using electrostatic levitation and droplet oscillation techniques in microgravity, which provide access to quiescent liquid droplets and deep supercooling of 510–815 K below the equilibrium melting points. Densities were measured over 900–2400 K. Viscosities were similar for all four compositions, increasing from ~10 mPa s near 2100 K to ~30 mPa s near 1750 K. Surface tensions were 450–490 dyn cm -1 for the rare earth titanates and 383–395 dyn cm -1 for the barium titanates; surface tensions of all compositions had small or negligible temperature dependence over 1700–2200 K. In conclusion, for solids recovered after melt quenching, x-ray microtomography revealed the fracture mechanics in crystalline products and minimal internal porosity in glass products, likely arising from entrapped gas bubbles. Internal microstructures were generally similar for products processed either in microgravity or in a terrestrial aerodynamic levitator.

42 ENGINEERING↗

Thermal relaxation of strain and twist in ferroelectric hexagonal boron nitride moiré interfaces

New properties can arise at van der Waals (vdW) interfaces hosting a moiré pattern generated by interlayer twist and strain. However, achieving precise control of interlayer twist/strain remains an ongoing challenge in vdW heterostructure assembly, and even subtle variation in these structural parameters can create significant changes in the moiré period and emergent properties. Characterizing the rate of interlayer twist/strain relaxation during thermal annealing is critical to establish a thermal budget for vdW heterostructure construction and may provide a route to improve the homogeneity of the interface or to control its final state. Here, we characterize the spatial and temporal dependence of interfacial twist and strain relaxation in marginally-twisted hBN/hBN interfaces heated under conditions relevant to vdW heterostructure assembly and typical sample annealing. We find that the ferroelectric hBN/hBN moiré at very small twist angles (θ ≤ 0.1°) relaxes minimally during annealing in air at typical assembly temperatures of 170°C. However, at 400°C, twist angle relaxes significantly, accompanied by a decrease in spatial uniformity. Uniaxial heterostrain initially increases and then decreases over time, becoming increasingly non-uniform in direction. Structural irregularities such as step edges, contamination bubbles, or contact with the underlying substrate result in local inhomogeneity in the rate of relaxation.

2D materials↗

Neutron-producing gas puff Z-pinch experiments on a fast, low-impedance, 0.5 MA linear transformer driver

A study on the neutron production from single and double gas puff Z-pinches on the CESZAR linear transformer driver with ~0.45 MA current and 170 ns rise time is presented. Total neutron yield measurements made with a LaBr activation detector are compared for three configurations, using a double nozzle setup. When a single, hollow, deuterium gas shell was used, reliable implosions could only be attained at higher load mass than the optimal value to match implosion time with the driver rise time, with neutron yields of ~10 6 per pulse. The use of a double gas puff configuration with a deuterium center jet allowed a reduction in the shell density and operation closer to machine-matched conditions, recording up to (4.1 ± 0.3) × 10 7 neutrons/pulse when either Kr or D 2 was used in the shell. For a comparable mass and implosion time, using a higher atomic-number gas in the outer shell results in more unstable plasma surface and smaller plasma radius at the location of instability bubbles, which, however, do not seem to consistently correlate with a higher neutron yield. Comparing implosion dynamics with models and neutron yields with literature scaling suggests that the machine current is not well coupled to the plasma during the final stages of compression. Optimizing current and energy coupling to the pinched plasma is critical to improving performance, particularly in low-impedance drivers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High repetition-rate 0.5 Hz broadband neutron source driven by the Advanced Laser Light Source

Neutron beams are an essential tool to investigate material structure and perform nondestructive analysis, as they give unique access to element composition, thus ideally complementing density analysis allowed by standard x-rays investigation. Laser-driven neutron sources, though compact and cost-effective, currently have lower average flux than conventional neutron sources, due to the limited repetition rate of the lasers used so far. However, advancements in laser technology allow nowadays to address this challenge. Here, we report results obtained at the Advanced Laser Light Source characterizing stable production of broadband (0.1–2 MeV) neutrons produced at a high repetition rate (0.5 Hz). The interaction of laser pulses of 22 fs duration and 3.2 J on-target energy with 2-μm-thick tantalum targets produced protons in the Target Normal Sheath Acceleration (TNSA) regime up to 7.3 MeV. These protons were subsequently converted into neutrons by (p,n) reactions in lithium fluoride (LiF). Activation measurements and bubble detectors were used to characterize neutron emissions, with a neutron fluence of up to ∼1.4×105 neutrons/shot/sr and energies mainly between a few hundred of kilo-electron volt and 2 MeV. The total neutron yield was ∼5×105 neutrons/shot. This paves the way for numerous applications, e.g., in homeland security, materials science, or cultural heritage.

Physics↗

Resolving discrepancies in bang-time predictions for indirect-drive ICF experiments on the NIF: Insights from the Build-A-Hohlraum campaign

This study investigated discrepancies between measured and simulated x-ray drive in Indirect-Drive Inertial Confinement Fusion (ID-ICF) hohlraums at the National Ignition Facility. Despite advances in radiation-hydrodynamic simulations, a consistent “drive deficit” remains. Experimentally measured ID-ICF capsule bang-times are systematically 400–700 ps later than simulations predict. The Build-A-Hohlraum (BAH) campaign explored potential causes for this discrepancy by systematically varying hohlraum features, including laser entrance hole (LEH) windows, capsules, and gas fills. Overall, the agreement between simulated and experimental x-ray drive was found to be largely unaffected by these changes. The data allow us to exclude some hypotheses put forward to potentially explain the discrepancy. Errors in the local thermodynamic equilibrium (LTE) atomic modeling, errors in the modeling of LEH closure, and errors due to a lack of plasma species mix physics in simulations are shown to be inconsistent with our measurements. Instead, the data support the hypothesis that errors in NLTE emission modeling are a significant contributor to the discrepancy. X-ray emission in the 2–4 keV range is found to be approximately 30% lower than in simulations. This is accompanied by higher than predicted electron temperatures in the gold bubble region, pointing to errors in non-LTE modeling. Introducing an opacity multiplier of 0.87 on energy groups above 1.8 keV improves agreement with experimental data, reducing the bang-time discrepancy from 300 to 100 ps. These results underscore the need for refined NLTE opacity models to enhance the predictive power of hohlraum simulations.

Band emission↗