Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hydrodynamic instability”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The Evolution of Inclined Binary Black Holes in the Disks of Active Galactic Nuclei

The accretion disks that fuel active galactic nuclei (AGNs) may house numerous stars and compact objects, formed in situ or captured from nearby star clusters. Embedded neutron stars and black holes may form binaries and eventually merge, emitting gravitational waves detectable by LIGO/VIRGO. AGN disks are a particularly promising environment for the production of high-mass gravitational-wave events involving black holes in the pair-instability mass gap, and may facilitate electromagnetic counterparts to black hole binary mergers. However, many orders of magnitude separate the typical length scales of binary formation and those on which gravitational waves can drive binary inspirals, making binary mergers inside the disk uncertain. Previous hydrodynamical simulations of binaries have either been restricted to two dimensions entirely, or focused on binaries aligned with the midplane of the disk. Herein we present the first three-dimensional, high-resolution, local-shearing-box, inviscid hydrodynamical simulations of disk-embedded binaries over a range of orbital inclinations. We find that retrograde binaries can shrink up to 4 times as quickly as prograde binaries, and that all binaries not perfectly aligned (or anti-aligned) with the AGN disk are driven into alignment. An important consequence of this is that initially retrograde binaries will traverse the inclinations where von Zeipel–Lidov–Kozai oscillations can drive binary eccentricities to large values, potentially facilitating mergers. We also find that interactions with the AGN disk may excite eccentricities in retrograde binaries and cause the orbits of embedded binaries to precess.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

RANS Simulation of Variable Density Turbulent Round Jets with Coflow using xRAGE Hydrodynamic Code and BHR Turbulence Models

This work for the fiscal year 2023 (FY23) is a continuation of previous efforts to evaluate the BHR turbulence models for their ability to accurately simulate variable density turbulent round jets with coflow. As before, RANS simulations are carried out using the xRAGE hydrodynamic code. The following are some of the previous findings. Israel showed that i) the symmetry boundary conditions for the BHR models in axisymmetric simulations were in error, ii) three grids of different resolutions did not lead to converging solutions, and iii) BHR 3.1 simulation exhibited instabilities and did not reach a steady state. Saenz and Rauenzahn derived and implemented into xRAGE the correct BHR boundary conditions at the symmetry axis. Cline conducted sensitivity studies with various parameters including the BHR models (versions 2, 3.1 and 4), gravity, material pressure, specific heat, initial turbulent kinetic energy and initial turbulent length scale, and found that the largest factor impacting on simulation results was the BHR model version, followed by the initial turbulent length scale. In addition, freeze boundary conditions at the exit and the side wall of the computational domain were used to remove anomalous flow behavior. Cline adjusted the inlet jet velocity, initial turbulent kinetic energy and initial turbulent length scale to obtain the best reasonable match with the experimental data of Charonko and Prestridge. It was found that the BHR 2 and 3.1 models performed in a similar manner, but the BHR 2 model produced much lower levels of density-specific-volume covariance and turbulent kinetic energy. The main focus for the FY23 is to study the effects of computational parameters related to the boundary conditions, mesh, domain size and timestep size. The reasoning behind this is that, unless simulation results can be shown to be reasonably independent from the aforementioned computational parameters, it would be difficult to attribute any discrepancies between simulation and experimental results to turbulence models. This important aspect has largely been overlooked in the previous years, and therefore will be studied comprehensively here. Additionally, effects of varying the initial turbulent length scale will be examined because it was previously identified as a major factor affecting the flow fields.

42 ENGINEERING↗
DOE OSTI2025

A fluctuating hydrodynamics model for nanoscale surfactant-laden interfaces

A multispecies diffuse interface model is formulated in a fluctuating hydrodynamics framework for the purpose of simulating surfactant interfaces at the nanoscale. The model generalizes previous work to ternary mixtures, employing a Cahn-Hilliard free energy density combined with incompressible, isothermal fluctuating hydrodynamics where dissipative fluxes include both deterministic and stochastic terms. The intermolecular parameters in the free energy are chosen such that one species acts as a partially miscible surfactant. From Laplace pressure measurements, we show that in this model the surface tension decreases linearly with surfactant concentration, leading to Marangoni convection for interfaces with concentration gradients. In the capillary wave spectrum for interfaces with and without surfactant, we find that for the former, the spectrum deviates significantly from classical capillary wave theory, presumably due to Gibbs elasticity. In non-equilibrium simulations of the Rayleigh-Plateau instability, deterministic simulations showed that the surfactant delays pinching of a fluid cylinder into droplets. However, stochastic simulations indicate that thermal fluctuations disrupt the surfactant's stabilizing effect. Similarly, the spreading of a patch of surfactant, driven by Marangoni convection, was found to be partially suppressed by thermal fluctuations.

Capillary waves↗
DOE OSTI2024

Identification of Kelvin-Helmholtz generated vortices in magnetised fluids

The Kelvin-Helmholtz Instability (KHI), arising from velocity shear across the magnetopause, plays a significant role in the viscous-like transfer of mass, momentum, and energy from the shocked solar wind into the magnetosphere. While the KHI leads to growth of surface waves and vortices, suitable detection methods for these applicable to magnetohydrodynamics (MHD) are currently lacking. A novel method is derived based on the well-established λ-family of hydrodynamic vortex identification techniques, which define a vortex as a local minimum in an adapted pressure field. The J × B Lorentz force is incorporated into this method by using an effective total pressure in MHD, including both magnetic pressure and a pressure-like part of the magnetic tension derived from a Helmholtz decomposition. The λ MHD method is shown to comprise of four physical effects: vortical momentum, density gradients, fluid compressibility, and the rotational part of the magnetic tension. A local three-dimensional MHD simulation representative of near-flank magnetopause conditions (plasma β’s 0.5 – 5 and convective Mach numbers M f ∼ 0.4) under northward interplanetary magnetic field (IMF) is used to validate λ MHD . Analysis shows it correlates well with hydrodynamic vortex definitions, though the level of correlation decreases with vortex evolution. Overall, vortical momentum dominates λ MHD at all times. During the linear growth phase, density gradients act to oppose vortex formation. By the highly nonlinear stage, the formation of small-scale structures leads to a rising importance of the magnetic tension. Compressibility was found to be insignificant throughout. Finally, a demonstration of this method adapted to tetrahedral spacecraft observations is performed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Three-dimensional signatures of self-similarity in a high-energy-density plasma shear-driven mixing layer

A hydrodynamic shear mixing layer experiment at the National Ignition Facility had previously demonstrated Eulerian scaling of integrated, late-time quantities, including turbulent kinetic energy. In this manuscript, the experiment is repeated with new materials. In this work, using the new dataset, we demonstrate that Euler-number scalings hold not just for late time, but dynamically throughout the experiment, for measurements in all three spatial dimensions. Incorporating the dynamic scaling leads to an enhanced calculation that the heavier of the two scaled experiments has approached three generations of mergers of its primary instability's structures and a consistent observation of such a merger in action in the lighter of the two scaled experiments. Furthermore, the improved scrutiny of the time evolution of instability structures leads to sharper estimates of turbulent kinetic energy, including a demonstration of different behaviors correlating with surface roughness (quantitatively consistent with transitions between laminar and turbulent initial states), as predicted by a Reynolds-averaged turbulent model, which evidently correctly handles the differing shock-roughness interactions to drive its internal state of the model into different regimes. Altogether, a picture arises of the analytical improvements in treating these variations (of times, densities, and roughnesses) as a unified whole and of multiple ways by which deviations from the scaling could indicate an onset of non-hydrodynamic behavior. Such deviations were not expected for these experiments (which models correctly indicated would remain hydrodynamic) but could be introduced by, for example, imposing external fields or increasing drive energy to test conditions relevant to inertial confinement fusion or other high-energy-density experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Nonlinear second order electromagnetic gyrokinetic theory for a tokamak plasma

The steep plasma pressure gradient that forms at the edge of the high confinement, H-mode regime of tokamak operation provides free energy to drive electromagnetic micro-instabilities that are widely believed to influence the transport processes in this so-called pedestal region. This high pressure gradient also provides a high current density (bootstrap current), known to influence ballooning mode stability and to be important for driving kink modes in the ideal magneto-hydrodynamic plasma model (so-called peeling-ballooning modes). Furthermore, efficient, steady state future tokamak power plants must operate with a large bootstrap current in the core and especially concerning spherical tokamaks, confinement will be influenced by electromagnetic turbulence. To accommodate these important situations, conventional electromagnetic gyrokinetic theory is extended to incorporate neoclassical effects in the equilibrium drives, allowing $B_{\vartheta} ~ B_0$ (B0 is the confining magnetic field, and $B_{\vartheta}$ is its poloidal component). This provides a global gyrokinetic model that self-consistently captures the consequences of large bootstrap current fractions on the equilibrium distribution functions.

kink modes↗
DOE OSTI2021

Shock Ignition Laser-Plasma Interactions in Ignition-Scale Plasmas

Here, we use a subignition scale laser, the 30 kJ Omega, and a novel shallow-cone target to study laser-plasma interactions at the ablation-plasma density scale lengths and laser intensities anticipated for direct drive shock-ignition implosions at National Ignition Facility scale. Our results show that, under these conditions, the dominant instability is convective stimulated Raman scatter with experimental evidence of two plasmon decay (TPD) only when the density scale length is reduced. Particle-in-cell simulations indicate this is due to TPD being shifted to lower densities, removing the experimental back-scatter signature and reducing the hot-electron temperature. The experimental laser energy-coupling to hot electrons was found to be 1%–2.5%, with electron temperatures between 35 and 45 keV. Radiation-hydrodynamics simulations employing these hot-electron characteristics indicate that they should not preheat the fuel in MJ-scale shock ignition experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Shock-driven three-fluid mixing with various chevron interface configurations

When a shock wave crosses a density interface, the Richtmyer–Meshkov instability causes perturbations to grow. Richtmyer–Meshkov instabilities arise from the deposition of vorticity from the misaligned density and pressure gradients at the shock front. In many engineering applications, microscopic surface roughness will grow into multi-mode perturbations, inducing mixing between the fluid on either side of an initial interface. Applications often have multiple interfaces, some of which are close enough to interact in the later stages of instability growth. In this study, we numerically investigate the mixing of a three-layer system with periodic zigzag (or chevron) interfaces, calculating the dependence of the width and mass of mixed material on properties such as the shock timing, chevron amplitude, multi-mode perturbation spectrum, density ratio, and shock mach number. The multi-mode case is also compared with a single-mode perturbation. The Flash hydrodynamic code is used to solve the Euler equations in three dimensions with adaptive grid refinement. Key results include a significant increase in mixed mass when changing from a single-mode to a multi-mode perturbation on one of the interfaces. The mixed width is mainly sensitive to the density ratio and chevron amplitude, whereas the mixed mass also depends on the multi-mode spectrum. In conclusion, steeper initial perturbation spectra have lower mixed mass at early times but a greater mixed mass after the reflected shock transits back across the layer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Fast flavor instability in hypermassive neutron star disk outflows

Here, we examine the effect of neutrino flavor transformation by the fast flavor instability (FFI) on long-term mass ejection from accretion disks formed after neutron star mergers. Neutrino emission and absorption in the disk set the composition of the disk ejecta, which subsequently undergoes r-process nucleosynthesis upon expansion and cooling. Here we perform 28 time-dependent, axisymmetric, viscous-hydrodynamic simulations of accretion disks around hypermassive neutron stars (HMNSs) of variable lifetime, using a 3-species neutrino leakage scheme for emission and an annular-lightbulb scheme for absorption. We include neutrino flavor transformation due the FFI in a parametric way, by modifying the absorbed neutrino fluxes and temperatures, allowing for flavor mixing at various levels of flavor equilibration, and also in a way that aims to respect the lepton-number preserving symmetry of the neutrino self-interaction Hamiltonian. We find that for a promptly-formed black hole (BH), the FFI lowers the average electron fraction of the disk outflow due to a decrease in neutrino absorption, driven primarily by a drop in electron neutrino/antineutrino flux upon flavor mixing. For a long-lived HMNS, the disk emits more heavy lepton neutrinos and reabsorbs more electron neutrinos than for a BH, with a smaller drop in flux compensated by a higher neutrino temperature upon flavor mixing. The resulting outflow has a broader electron fraction distribution, a more proton-rich peak, and undergoes stronger radiative driving. Disks with intermediate HMNS lifetimes show results that fall in between these two limits. In most cases, the impact of the FFI on the outflow is moderate, with changes in mass ejection, average velocity, and average electron fraction of order ~10%, and changes in the lanthanide/actinide mass fraction of up to a factor ~2.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The Effect of Spatially Varying Collision Frequency on the Development of the Rayleigh–Taylor Instability

Abstract The Rayleigh–Taylor (RT) instability is ubiquitously observed, yet has traditionally been studied using ideal fluid models. Collisionality can vary strongly across the fluid interface, and previous work demonstrates the necessity of kinetic models to completely capture dynamics in certain collisional regimes. Where previous kinetic simulations used spatially and temporally constant collision frequency, this work presents five-dimensional (two spatial, three velocity dimensions) continuum-kinetic simulations of the RT instability using a more realistic spatially varying collision frequency. Three cases of collisional variation are explored for two Atwood numbers: low to intermediate, intermediate to high, and low to high. The low-to-intermediate case exhibits no RT instability growth, while the intermediate-to-high case is similar to a fluid-limit kinetic case with interface widening biased toward the lower-collisionality region. A novel contribution of this work is the low-to-high collisionality case that shows significantly altered instability growth through an upward movement of the interface and damped spike growth due to increased free-streaming particle diffusion in the lower region. Contributions to the energy flux from the non-Maxwellian portions of the distribution function are not accessible to fluid models and are greatest in magnitude in the spike and regions of low collisionality. Increasing the Atwood number results in greater RT instability growth and reduced upward interface movement. Deviation of the distribution function from Maxwellian is inversely proportional to collision frequency and concentrated around the fluid interface. The linear phase of RT instability growth is well described by theoretical linear growth rates accounting for viscosity and diffusion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Exploration of Rayleigh-Taylor seeding mechanisms in laser- and pulsed-power-driven inertial fusion

The ubiquitous Rayleigh-Taylor instability (RTI) has been observed in fluids and plasmas in a wide range of parameter regimes from astrophysical to laboratory plasmas. In high-energy-density (HED) laboratory plasmas, such as laser-driven inertial confinement fusion (ICF) plasmas, the RTI can have a detrimental impact on achieving fusion ignition through the generation of hydrodynamic mix. Hence, an understanding of the seeding mechanisms that produce the RTI and identifying ways to mitigate hydrodynamic mix is of critical importance. Appropriately aligned magnetic fields have been known to stabilize short-wavelength RTI. A number of numerical and experimental studies have demonstrated the benefits of using imposed magnetic fields in laser-driven ICF to achieve higher ion temeratures, higher neutron yields, and a relaxation of the ignition criteria. This work addresses seeding mechanisms in laser-driven implosions to understand critical early-stage physics that ultimately leads to substantial growth of the RTI along with mechanisms for mitigation of this growth. Surface perturbations due to machining tolerances and single-feature seeds (for example due to fill tubes) can produce substantial RTI growth in the ignition-relevant high-convergence ratio targets for laser- driven ICF implosions. RTI growth from single-feature seeding in laser-based implosions has the potential to be mitigated through appropriately aligned externally applied magnetic fields and the goal of this work has been to quantify that numerically and experimentally. Resistive-magnetohydrodynamic (MHD) simu- lations are used to study the seeding and evolution of the RTI leveraging previous support from the DOE HEDLP program. Furthermore, experimental data has been obtained and applied towards code validation of unmagnetized and magnetized evolution of single-feature seeded RTI growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Turbulent dynamo action and its effects on the mixing at the convective boundary of an idealized oxygen-burning shell

Convection is one of the most important mixing processes in stellar interiors. Hydrodynamic mass entrainment can bring fresh fuel from neighboring stable layers into a convection zone, modifying the structure and evolution of the star. Because flows in stellar convection zones are highly turbulent, multidimensional hydrodynamic simulations are fundamental to accurately capture the physics of mixing processes. Under some conditions, strong magnetic fields can be sustained by the action of a turbulent dynamo, adding another layer of complexity and possibly altering the dynamics in the convection zone and at its boundaries. In this study, we used our fully compressible S EVEN -L EAGUE H YDRO code to run detailed and highly resolved three-dimensional magnetohydrodynamic simulations of turbulent convection, dynamo amplification, and convective boundary mixing in a simplified setup whose stratification is similar to that of an oxygen-burning shell in a star with an initial mass of 25 M ⊙ . We find that the random stretching of magnetic field lines by fluid motions in the inertial range of the turbulent spectrum (i.e., a small-scale dynamo) naturally amplifies the seed field by several orders of magnitude in a few convective turnover timescales. During the subsequent saturated regime, the magnetic-to-kinetic energy ratio inside the convective shell reaches values as high as 0.33, and the average magnetic field strength is ∼10 10 G. Such strong fields efficiently suppress shear instabilities, which feed the turbulent cascade of kinetic energy, on a wide range of spatial scales. The resulting convective flows are characterized by thread-like structures that extend over a large fraction of the convective shell. The reduced flow speeds and the presence of magnetic fields with strengths up to 60% of the equipartition value at the upper convective boundary diminish the rate of mass entrainment from the stable layer by ≈20% as compared to the purely hydrodynamic case.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Three-dimensional Core-collapse Supernova Simulations with 160 Isotopic Species Evolved to Shock Breakout

We present three-dimensional simulations of core-collapse supernovae using the FLASH code that follow the progression of the explosion to the stellar surface, starting from neutrino radiation hydrodynamic simulations of the neutrino-driven phase performed with the Chimera code. We consider a 9.6 M ⊙ zero-metallicity progenitor starting from both 2D and 3D Chimera models and a 10 M ⊙ solar-metallicity progenitor starting from a 2D Chimera model, all simulated until shock breakout in 3D while tracking 160 nuclear species. The relative velocity difference between the supernova shock and the metal-rich Rayleigh–Taylor (R-T) "bullets" determines how the metal-rich ejecta evolves as it propagates through the density profile of the progenitor and dictates the final morphology of the explosion. We find maximum 56 Ni velocities of ~1950 and ~1750 km s -1 at shock breakout from 2D and 3D 9.6 M ⊙ Chimera models, respectively, due to the bullets' ability to penetrate the He/H shell. When mapping from 2D, we find that the development of higher-velocity structures is suppressed when the 2D Chimera model and 3D FLASH model meshes are aligned. The development of faster-growing spherical-bubble structures, as opposed to the slower-growing toroidal structure imposed by axisymmetry, allows for interaction of the bullets with the shock and seeds further R-T instabilities at the He/H interface. We see similar effects in the 10 M ⊙ model, which achieves maximum 56 Ni velocities of ~2500 km s -1 at shock breakout.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Magnetorotational instability in a swirling partially ionized gas

ABSTRACT The magnetorotational instability (MRI) has been proposed as the method of angular momentum transport that enables accretion in astrophysical discs. However, for weakly ionized discs, such as protoplanetary discs, it remains unclear whether the combined non-ideal magnetohydrodynamic (MHD) effects of Ohmic resistivity, ambipolar diffusion, and the Hall effect make these discs MRI stable. While much effort has been made to simulate non-ideal MHD MRI, these simulations make simplifying assumptions and are not always in agreement with each other. Furthermore, it is difficult to directly observe the MRI astrophysically because it occurs on small scales. Here, we propose the concept of a swirling gas experiment of weakly ionized argon gas between two concentric cylinders threaded with an axial magnetic field that can be used to study non-ideal MHD MRI. For our proposed experiment, we derive the hydrodynamic equilibrium flow and a dispersion relation for MRI that includes the three non-ideal effects. We solve this dispersion relation numerically for the parameters of our proposed experiment. We find it should be possible to produce a non-ideal MRI in such an experiment because of the Hall effect, which increases the MRI growth rate when the vertical magnetic field is anti-aligned with the rotation axis. As a proof of concept, we also present experimental results for a hydrodynamic flow in an unmagnetized prototype. We find that our prototype has a small, but non-negligible, α-parameter that could serve as a baseline for comparison to our proposed magnetized experiment, which could be subject to additional turbulence from the MRI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Exploring implosion designs for increased compression on the National Ignition Facility using high density carbon ablators

It has long been recognized that high compression, and hence good confinement, is essential to achieving high yields in inertial confinement fusion implosions. In pursuit of multi-megajoule yields on the National Ignition Facility (NIF), a new campaign has begun aimed at testing the hypothesis that controlling hydrodynamic stability is key to achieving effective higher compression with the high density carbon ablators currently fielded on NIF. This campaign is built around a new implosion design, called SQ-n, that is derived from the uniquely stable Bigfoot design tested on NIF in 2016–2019. While very stable and with performance that was quite close to one-dimensional expectations, Bigfoot was a relatively high adiabat, and consequently lower compression design. The goal of SQ-n is then to evolve Bigfoot toward a higher compression design but without compromising its unique stability characteristics. Specifically, SQ-n adopts a ramped foot pulse shape to minimize early time Richtmyer–Meshkov instability growth and uses an ablator dopant distribution extending all of the way to the fuel–ablator interface that simulations suggest further reduces perturbation growth. In conclusion, this paper describes the design philosophy pursued with SQ-n, the results of instability modeling of the candidate design, and the experimental campaign planned to test these ideas in the near future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Density variance dynamics in disparate shock tubes

This report discusses two experiments which investigate thin layer, heavy curtain fragmentation from the perspective of a Reynolds-Averaged mix model, in drastically disparate experimental regimes. The first, the centimeter/millisecond-scale “Horizontal Shock Tube” (HST) is a compressed-gas piston-driven shock tube experiment. The second, the “Multishock thin layer” (Mshock) experiment performed at the National Ignition Facility, is a micrometer/nanosecond-scale laser-driven shock tube experiment. Both are situations in which a heavy plane layer (a ‘curtain’) is initially suspended in a lighter medium. After being shocked from at least one side, the layer translates while its interfaces evolve due to the excitation of the Richtmyer-Meshkov instability at its surfaces. The evolution of density variance, which initially exists only on the surface of the layer, as it comes to encompass the whole layer interior is used as a description of layer fragmentation and dissolution. These experiments have each been simulated in the Los Alamos National Laboratory multi-physics code xRAGE, which includes fundamental hydrodynamics, extended plasma physics and radiation effects which are important to drive the high-energy density experiment, and the Besnard-Harlow-Rauenzahn (BHR) turbulence model. In each, the principal diagnostic for comparison is an experimental metric for the density (co)variance, b, which tracks the moments of the density field at the curtain interfaces and body. Due to experimental constraints in different regimes (i.e. optical diagnostics can be deployed on conventional shock tubes, while the plasma shock tubes must be imaged by x-rays; interfaces can be imposed to specification on laser-driven experiments, which are stored in the solid phase, while conventional experiments have imperfect control of the flow fields which separate the layer, etc.) the experiments are not designed to be perfect scaled cognates of one another. However, despite the separation of six orders of magnitude of scaling in time, and four in space, we are able to demonstrate that the same turbulence model, operating in the same fashion in the same computer code, is able to reproduce results in each experiment, by tracking evolution due to common relevant physics. Additionally, we will present preliminary work toward density variance comparisons in a single-interface Richtmyer-Meshkov configuration, the conventional fluid “Vertical Shock Tube” (VST) experiment, and the Modal Initial Conditions (ModCons) campaign fielded at the OMEGA-EP laser facility.

42 ENGINEERING↗
DOE OSTI2026

Geometric optics model of Thomson scattering enhanced by parametric coupling

We formulate a theory of Thomson scattering (TS) from inhomogeneous plasmas using a ray-tracing approach. The scattered radiation is described within the eikonal approximation, including convective gain arising from scattering instabilities driven by the TS probe laser. In the limit of a weak probe and vanishing gain, the theory recovers the standard TS cross section and scattered power due to volumetric particle noise. The ray-tracing formalism is implemented in a TS fitting framework, enabling accurate analysis of inhomogeneous plasmas whose background profiles can be provided by large-scale hydrodynamic simulations or experimental diagnostics. Here, several examples of synthetic TS spectra are presented, and further extensions of the theory are discussed.

Carleton, D. E. [University of Alberta, Edmonton, ↗
DOE OSTI2020

Gravitational-wave signal of a core-collapse supernova explosion of a 15M ⊙ star

In this paper we report on the gravitational-wave signal computed in the context of a three-dimensional simulation of a core-collapse supernova explosion of a 15 M ⊙ star. The simulation was performed with our neutrino hydrodynamics code chimera . We detail the gravitational wave strains as a function of time, for both polarizations, and discuss their physical origins. We also present the corresponding spectral signatures. Gravitational wave emission in our model has two key features: low-frequency emission (less than 200 Hz) emanates from the gain layer as a result of neutrino-driven convection and the standing accretion shock instability (SASI), and high-frequency emission (greater than 600 Hz) emanates from the proto–neutron star due to Ledoux convection within it. The high-frequency emission dominates the gravitational wave emission in our model and emanates largely from the convective layer itself, not from the convectively stable layer above it, due to convective overshoot. Moreover, the low-frequency emission emanates from the gain layer itself, not from the proto–neutron star, due to accretion onto it. We provide evidence of the SASI in our model and demonstrate that the peak of our low-frequency gravitational wave emission spectrum corresponds to it. Given its origin in the gain layer, we classify the SASI emission in our model as p-mode emission and assign a purely acoustic origin, not a vortical–acoustic origin, to it. We compare the results of our three-dimensional model analysis with those obtained from the model’s two-dimensional counterpart and find a significant reduction in the strain amplitudes in the former case, as well as significant reductions in all related quantities. Our dominant proto–neutron star gravitational wave emission is not well characterized by emission from surface g modes, complicating the relationship between peak frequencies observed and the mass and radius of the proto–neutron star expressed by analytic estimates under the assumption of surface g-mode emission. We present our frequency normalized characteristic strain along with the sensitivity curves of current- and next-generation gravitational wave detectors. This simple analysis indicates that the spectrum of gravitational wave emission between approximately 20 Hz and approximately 1 kHz, stemming from neutrino-driven convection, the SASI, accretion onto the proto–neutron star, and proto–neutron star convection, will be accessible for a Galactic event.

79 ASTRONOMY AND ASTROPHYSICS↗