Search NASA⌕ Search

SEARCH · Search NASA

Results for “EQUATIONS - BOUNDARY LAYER - PRANDTL”

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 235 records · Page 13

Exploring the Neutron Substructure with Advanced Polarized Helium-3 Targets (Or: How I Learned to Stop Worrying and Love Spectroscopy)

As we seek to understand the smallest, physical aspects of our universe, we cannot simply rely on our senses to probe the world around us as we did in the past. The smallest physical elements of our universe behave in strange, probabilistic ways and are completely invisible to the naked eye/ear/etc. So, we design clever experiments (such as scattering experiments) to probe these minute realms. Then, just as with the larger, observable world, we devise models and equations to describe what we think is happening. Due to the nature of the physical universe at the quantum scale and with the aid of symmetries such as Lorentz invariance, we can write down equations that describe the scattering, but the expressions contain functions, which we call ?form factors? and ?structure functions?, that we cannot compute from first principles. We can, however, formulate models that make predictions for these functions. By comparing our predictions with the observed data, we can gain insight into the validity of our models and thus a better physical understanding of what is happening at these minuscule scales. Studying the constituents inside of the nucleus of an atom adds another layer of difficulty if we can?t remove those components from the nucleus. This is the case with the neutron. When not bound in the nucleus with protons and other neutrons, the neutron will decay into a proton after about 15 minutes. So, we?re forced to study the neutron while it is still bound in the nucleus of an atom such as helium-3 (3He). For the last 1,000 years (rounding up), our group has developed high quality, polarized 3He targets made of an aluminosilicate glass. These targets are made in order to perform experiments at Jefferson Lab (JLab), experiments which let us determine the form factors and structure functions of the neutron by scattering polarized electrons from polarized neutrons (or rather polarized 3He). The specific experiments reported on in this thesis push the bounds of our understanding of the internal structure of the neutron. Good science is often about pushing experimental techniques to a new level. Toward that goal we study our polarized 3He targets both to advance the technology and to choose the best ones for our experiments. We do this using a process called nuclear magnetic resonance (NMR) to gauge the maximum polarization of a target and how fast the polarization decays with time. While these tests primarily provide us information that make analysis of our experimental scattering data possible, they also let us determine whether or not a target-cell is useful or even, dare I say, of spectacular quality. Our latest targets utilize a novel convection design allowing 3He to be polarized and quickly moved in front of the electron-beam, making it possible to use larger targets with higher electron-beam currents than ever before. This means more electrons scatter and we get more data. And by studying our targets in detail prior to using them in our experiments, we have found techniques to take effects which could have been detrimental to target quality and turn them to our advantage! It?s a real case of making lemonade out of lemons. We also use laser spectroscopy to study the absorption lines of alkali-metals in the target (potassium and rubidium, specifically). We add these alkali-metals to our target to facilitate polarizing the 3He. We can use the measurement of these pressure broadened absorption lines to determine the 3He density inside of the target with great precision. Historically, we understood the width of these lines would be dependent on the temperature of the target. Specifically, if I raise the temperature, the width should get bigger. I found that was not the case, which was very confusing at first, though very exciting now that I realize the data are self-consistent and suggestive of unexpected behavior. This thesis details the development of high quality, glass, polarized 3He targets for the 2020 An 1 /dn 2 and 2023 Gn E experiments, which utilized the first 3He convection targets and broke records in target quality. This thesis also covers the initial development of metal windows for the next-generation of 3He target-cells. Finally, this thesis documents the temperature dependence of the width of potassium (K) and rubidium (Rb) absorption lines as measured with laser spectroscopy.

Jantzi, Christopher↗

Asymmetric crossing of the attractive and repulsive magnetic potential by Abrikosov vortices

Here, we present studies of Abrikosov vortex motion across the magnetically charged edges of long thin ferromagnetic stripes placed above or under a thin superconducting film. The magnetic charges at the stripe edges form attractive or repulsive potential wells for vortices. Using a relatively small in-plane magnetic field to polarize the stripe edges and normal-to-plane magnetic field to induce up or down polarized vortices, we tune the attractive or repulsive stripe-vortex interactions. Imaging of the vortex dynamics reveals that the repulsive magnetic potential U m + acts as a robust vortex pinning barrier, while the attractive U m - has practically no effect on the vortex motion irrespective of the position of the stripes located above or underneath the superconducting film. We analyze the observed asymmetry using equations of the overdamped vortex motion. The formal analytical solution yields an asymmetry, but the numerical modeling with and without noise terms does not confirm it. Instead, we find that the asymmetry is caused by the creation of spontaneous vortices at the maxima of U m , which depends on the edge polarity and suppress U m - . Furthermore, our experiment and time-dependent Ginzburg-Landau simulations demonstrate that magnetic pinning dominates over vortex pinning due to corrugations of the superconducting layer deposited on top of the ferromagnetic stripes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quantum oscillations with angular dependence in PdTe 2 single crystals

The layered transition-metal dichalcogenide PdTe 2 has been discovered to possess bulk Dirac points as well as topological surface states. By measuring the magnetization (up to 7 T) and magnetic torque (up to 35 T) in single crystalline PdTe2, we observe distinct de Haas–van Alphen (dHvA) oscillations. Eight frequencies are identified with H||c, with two low frequencies (F α = 8 T and F β = 117 T) dominating the spectrum. The effective masses obtained by fitting the Lifshitz–Kosevich (LK) equation to the data are $m_α^* =0.059m_0$ and $m_β^* =0.067m_0$ where $m_0$ is the free electron mass. The corresponding Landau fan diagrams allow the determination of the Berry phase for these oscillations resulting in values of ~0.67π for the 3D α band (hole-type) (down to the 1st Landau level) and ~0.23π–0.73π for the 3D β band (electron-type) (down to the 3rd Landau level). By investigating the angular dependence of the dHvA oscillations, we find that the frequencies and the corresponding Berry phase (Φ B ) vary with the field direction, with a Φ B ~ 0 when H is 10°–30° away from the ab plane for both α and β bands. Finally, the multiple band nature of PdTe 2 is further confirmed from Hall effect measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Spectral broadening from turbulence in multiscale lower hybrid current drive simulations

Abstract The scattering of lower hybrid (LH) waves due to scrape-off layer filaments is investigated. It is revealed that scattering can account for the LH spectral gap without any ad hoc modification to the wave-spectrum. This is shown using a multiscale simulation approach which allows, for the first time, the inclusion of full-wave scattering physics in ray-tracing/Fokker–Planck calculations. In this approach, full-wave scattering probabilities are calculated for a wave interacting with a statistical ensemble of filaments. These probabilities are coupled to ray-tracing equations using radiative transfer theory. This allows the modeling of scattering along the entire ray-trajectory, which can be important in the multi-pass regime. Simulations are conducted for LH current drive (LHCD) in Alcator C-Mod, resulting in excellent agreement with experimental current and hard x-ray profiles. A region in filament parameter space is identified in which the impact of scattering on LHCD is saturated. Such a state coincides with experimental LHCD measurements, suggesting saturation indeed occurs in C-Mod, and therefore the exact statistical properties of the filaments are not important.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron-mediated entanglement of two distant macroscopic ferromagnets within a nonequilibrium spintronic device

Using the nascent concept of quantum spin-transfer torque, we demonstrate that a current pulse can be harnessed to entangle quantum localized spins of two spatially separated ferromagnets (FMs) which are initially unentangled. The envisaged setup is composed of a spin-polarizer (FM p ) and a spin-analyzer (FM a ) FM layers separated by a normal metal (NM) spacer. The injection of a current pulse into the device leads to a time-dependent superposition of many-body states characterized by a high degree of entanglement between the spin degrees of freedom of the two distant FM layers. The nonequilibrium dynamics are due to the transfer of spin angular momentum from itinerant electrons to the localized spins via a quantum spin-torque mechanism that remains active even for collinear but antiparallel arrangements of the FM p and FM a magnetizations (a situation in which the conventional spin torque is absent). We quantify the mixed-state entanglement generated between the FM layers by tracking the time evolution of the full density matrix and analyzing the build-up of the mutual logarithmic negativity over time. The effect of decoherence and dissipation in the FM layers due to coupling to bosonic baths at finite temperature, the use of multielectron current pulses, and the dependence on the number of spins are also considered in an effort to ascertain the robustness of our predictions under realistic conditions. Lastly, we propose a “current-pump–x-ray-probe” scheme, utilizing ultrafast x-ray spectroscopy, that can witness nonequilibrium and transient entanglement of the FM layers by extracting its time-dependent quantum Fisher information.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Surrogate Model Development of Spent Fuel Degradation for Repository Performance Assessment

In model simulations of deep geologic repositories, UO 2 fuel matrix degradation typically begins as soon as the waste package breaches and groundwater contacts the fuel surface. The initial degradation rate depends on the timing of these events, burnup of the fuel, temperature, and concentrations of dissolved reactants. Estimating the initial rate of degradation is fairly straightforward, but as UO 2 corrosion products precipitate on the fuel surface and the movement of dissolved species between the fuel surface and environment is impeded by the precipitated solids, the rate is more difficult to quantify. At that point, calculating the degradation rate becomes a reactive-transport problem in which a large number of equations must be solved by iteration for a large number of grid cells at each time step. The consequence is that repository simulations, which are already expensive, become much more expensive, especially when hundreds or thousands of waste packages breach. The Fuel Matrix Degradation (FMD) model is the process model of the Spent Fuel and Waste Science and Technology (SFWST) campaign of the US Department of Energy (DOE). It calculates spent fuel degradation rates as a function of radiolysis, redox reactions, electrochemical reactions, alteration layer growth, and diffusion of reactants through the alteration layer. Like other similar fuel degradation process models, it is a complicated model requiring a large number of calculations and iterations at each time step.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Spatiotemporal parallelization of an analytical heat conduction model for additive manufacturing via a hybrid OpenMP + MPI approach

The ability to do thermal simulations for entire additive manufacturing builds is a key computational problem facing the additive manufacturing community; however, complex numerical models considering multiple physical phenomena currently do not have the capacity for simulations at this scale. To this end, conduction only analytic models offer a viable approach due to the massive drop in computational expense. In this work, we extend an existing implementation which uses a governing equation which can be evaluated at any point in space and time. This implementation already utilizes OpenMP with a spatial decompositions scheme stemming from a melt pool tracking algorithm. Furthermore, we then combine this with a parallel in time (PinT) approach to make the problem highly parallelizable. The new scheme, which uses MPI for internode communication and OpenMP for intranode communication, is shown to scale very well across multiple computational nodes. This approach results in the ability to simulate the 3D solidification conditions for entire layers of additively manufactured parts in minutes making part scale thermal simulations more practical.

36 MATERIALS SCIENCE↗

A Free-Standing Boron-Doped Diamond Grid Electrode for Fundamental Spectroelectrochemistry

Spectroelectrochemistry (SEC) is a powerful technique that enables a variety of redox properties to be studied, including formal potential (E o ), thermodynamic values (ΔG, ΔH, ΔS), diffusion coefficient (D), electron transfer stoichiometry (n), and others. SEC requires an electrode which light can pass through while maintaining sufficient electrical conductivity. This has been traditionally composed of metal or metal oxide films atop transparent substrates like glass, quartz, or metallic mesh. Robust electrode materials like boron-doped diamond (BDD) could help expand the environments in which SEC can be performed, but most designs are limited to thin films (~100–200 nm) on transparent substrates less resilient than free-standing BDD. Here, this work presents a free-standing BDD grid electrode (G-BDD) for fundamental SEC measurements, using the well-characterized Fe(CN) 6 3–/4– redox couple as proof-of-concept. With a combination of cyclic voltammetry (CV), thin-layer SEC, and chronoabsorptometry, several of the redox properties mentioned above were calculated and compared. For E o' , n, and D, similar results were obtained when comparing the CV [E o' = +0.279 (±0.002) V vs Ag/AgCl; n = 0.97; D = 4.1 × 10 –6 cm 2 ·s –1 ] and SEC [E o' = +0.278 (±0.001) V vs Ag/AgCl; n = 0.91; D = 5.2 × 10 –6 cm 2 ·s –1 ] techniques. Both values align with what has been previously reported. To calculate D from the SEC data, modification of the classical equation used in chronoabsorptometry was required to accommodate the G-BDD electrode geometry. Overall, this work expands on the applicability of SEC techniques and BDD as a versatile electrode material.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Radiation GRMHD Models of Accretion onto Stellar-mass Black Holes. I. Survey of Eddington Ratios

We summarize results from a survey of radiation-dominated black hole accretion flows across a wide range of mass accretion rates, as well as two values of black hole spin and initial magnetic field geometry. These models apply an algorithm targeting direct solutions to the radiation transport equation in full general relativity and have been enabled by access to modern exascale computing systems. Super-Eddington accretion flows form geometrically thick radiation-pressure-supported disks that drive powerful equatorial outflows. A narrow funnel-shaped photosphere in the inner region results in very low radiative efficiencies in this regime. The structure of near- and sub-Eddington accretion depends on whether there is net vertical magnetic flux at the midplane of the disk. With net flux, the disk forms a thin, dense layer at the midplane surrounded by a magnetically dominated corona, whereas without net flux, the disk remains magnetically dominated everywhere. Although none of our models achieve the magnetically arrested disk regime, those with net vertical flux and a rapidly spinning black hole still produce powerful relativistic jets. Our calculations adopt simple opacity models (with scalings appropriate to stellar-mass black hole accretion). We discuss the application of our results to observations of X-ray binaries and ultraluminous X-ray sources such as Cyg X-3 and SS433. We also speculate on the application of our super-Eddington models to the interpretation of little red dots recently discovered by JWST.

79 ASTRONOMY AND ASTROPHYSICS↗

Facile and quantitative estimation of strain in nanobubbles with arbitrary symmetry in 2D semiconductors verified using hyperspectral nano-optical imaging

When layers of van der Waals materials are deposited via exfoliation or viscoelastic stamping, nanobubbles are sometimes created from aggregated trapped fluids. Though they can be considered a nuisance, nanobubbles have attracted scientific interest in their own right owing to their ability to generate large in-plane strain gradients that lead to rich optoelectronic phenomena, especially in the semiconducting transition metal dichalcogenides. Determination of the strain within the nanobubbles, which is crucial to understanding these effects, can be approximated using elasticity theory. However, the Föppl–von Kármán equations that describe strain in a distorted thin plate are highly nonlinear and often necessitate assuming circular symmetry to achieve an analytical solution. Here, we present an easily implemented numerical method to solve for strain tensors of nanobubbles with arbitrary symmetry in 2D crystals. The method only requires topographic information from atomic force microscopy and the Poisson ratio of the 2D material. We verify that this method reproduces the strain for circularly symmetric nanobubbles that have known analytical solutions. Finally, we use the method to reproduce the Grüneisen parameter of the E' mode for 1L-WS2 nanobubbles on template-stripped Au by comparing the derived strain with measured Raman shifts from tip-enhanced Raman spectroscopy, demonstrating the utility of our method for estimating localized strain in 2D crystals.

2D materials↗

Theoretical investigation of braking of tearing mode rotation by resistive walls in ITER

The locking of the 2/1 tearing mode to the resistive wall in the ITER tokamak (15 MA inductive scenario 2) is investigated theoretically using a cylindrical asymptotic matching model. The model takes into account the fact that ITER plasmas will effectively be surrounded by two walls; the inner blanket module layer with a time constant of about 23 ms, and the outer vacuum vessel with a time constant of about 380 ms. The model also takes cognizance of the fact that neither the blanket module layer nor the vacuum vessel can be accurately described as “thin” walls (in the ordinarily accepted sense). The model incorporates changes in both the plasma poloidal and the toroidal angular velocity profiles, in response to the electromagnetic braking torque that develops at the rational surface, because it turns out that neoclassical poloidal flow-damping is not strong enough to completely suppress changes in the poloidal velocity. Finally, the model accurately calculates changes in the poloidal and toroidal plasma angular velocity profiles by evolving the full angular equations of motion, taking the electromagnetic braking torque, plasma inertia, plasma viscosity, and poloidal flow-damping into account. The time required for the 2/1 tearing mode to grow from a small amplitude to a sufficient one to lock to the walls is found to be about 3.5 s. As a result, the critical full radial island width at which wall locking is triggered is found to be about 9% of the plasma minor radius.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Particle-in-cell Monte Carlo-collision modeling of non-ideal effects in wave-heated dense microplasmas

A computational model for non-ideal plasma effects during the time evolution of a second-stage laser-heated discharge at high pressures is presented. The model extends a classical one-dimensional particle-in-cell Monte Carlo-collision (PIC-MCC) approach coupled with Maxwell's equations for the laser-heating process of a xenon plasma at 300 K temperature and 10–100 bar pressure. Plasma non-ideality resulting from Coulomb coupling at high plasma densities is manifested as a depression in the effective ionization potential of atoms and enhanced collision cross sections. These non-ideal effects are represented using the Ecker–Kröll model in the context of the PIC-MCC approach. We find that full ionization of the plasma is obtained on the picosecond timescale, starting from the skin layer and quickly expanding throughout the domain through an anomalous extension of the skin depth. More critically, we show that the inclusion of the non-ideal plasma effects results in more rapid ionization compared to an ideal plasma, especially at higher pressures. The ionization delay reduction is of the order of a fraction of a picosecond, corresponding to a 16% decrease at 100 bar. As the amplitude of the wave field is lowered, the ionization rate is lowered, making the plasma non-ideality effects more prominent.

Solmaz, Evrim (ORCID:0000000211830881)↗

Radar Retrieval Evaluation and Investigation of Dendritic Growth Layer Polarimetric Signatures in a Winter Storm

Abstract This study evaluates ice particle size distribution and aspect ratio φ Multi-Radar Multi-Sensor (MRMS) dual-polarization radar retrievals through a direct comparison with two legs of observational aircraft data obtained during a winter storm case from the Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) campaign. In situ cloud probes, satellite, and MRMS observations illustrate that the often-observed K dp and Z DR enhancement regions in the dendritic growth layer can either indicate a local number concentration increase of dry ice particles or the presence of ice particles mixed with a significant number of supercooled liquid droplets. Relative to in situ measurements, MRMS retrievals on average underestimated mean volume diameters by 50% and overestimated number concentrations by over 100%. IWC retrievals using Z DR and K dp within the dendritic growth layer were minimally biased relative to in situ calculations where retrievals yielded −2% median relative error for the entire aircraft leg. Incorporating φ retrievals decreased both the magnitude and spread of polarimetric retrievals below the dendritic growth layer. While φ radar retrievals suggest that observed dendritic growth layer particles were nonspherical (0.1 ≤ φ ≤ 0.2), in situ projected aspect ratios, idealized numerical simulations, and habit classifications from cloud probe images suggest that the population mean φ was generally much higher. Coordinated aircraft radar reflectivity with in situ observations suggests that the MRMS systematically underestimated reflectivity and could not resolve local peaks in mean volume diameter sizes. These results highlight the need to consider particle assumptions and radar limitations when performing retrievals. significance statement Developing snow is often detectable using weather radars. Meteorologists combine these radar measurements with mathematical equations to study how snow forms in order to determine how much snow will fall. This study evaluates current methods for estimating the total number and mass, sizes, and shapes of snowflakes from radar using images of individual snowflakes taken during two aircraft legs. Radar estimates of snowflake properties were most consistent with aircraft data inside regions with prominent radar signatures. However, radar estimates of snowflake shapes were not consistent with observed shapes estimated from the snowflake images. Although additional research is needed, these results bolster understanding of snow-growth physics and uncertainties between radar measurements and snow production that can improve future snowfall forecasting.

Meteorology & Atmospheric Sciences↗

Implications of Sound Velocities of Natural Topaz on the Seismic L–Discontinuity

Compressional and shear wave velocities of topaz have been measured by ultrasonic interferometry in a multi-anvil apparatus up to 13 GPa at room temperature. By fitting the measured velocities under pressure to finite strain equations, the elastic moduli of topaz and their pressure derivatives were obtained, yielding K S0 = 165.5 (25) GPa, K S ' = 2.9 (1), G 0 = 116.5 (17) GPa, and G' = 1.0 (1). Modeling of the velocity profiles of subducted sediments suggests that at 8 GPa (~240 km depth) the P- and S- wave velocities of the subducted sediments exhibit a first order increase of 10% and 11%, respectively. Furthermore, the comparison between seismic observations and velocity contrasts of sediment + MORB model at 240 km shows dehydration of a 1.5–3.9 km thick sediment layer can potentially contribute to the seismically observed Lehmann discontinuity in Tonga, Japan, Philippine Sea, South America, Indonesia, and South Shetland subduction zones.

36 MATERIALS SCIENCE↗

Size Effect of Local Current-Voltage Characteristics of MX 2 Nanoflakes: Local Density of States Reconstruction from Scanning Tunneling Microscopy Experiments

Local current-voltage characteristics for low-dimensional transition-metal dichalcogenides (LDTMD), as well as the reconstruction of their local density of states (LDOS) from scanning tunneling microscopy (STM) experiments, are of fundamental interest and can be useful for advanced applications. Most of the existing models either have limited applicability for complex-shaped LDTMDs (e.g., those based on the Simmons approach) or require solving of an ill-defined integral equation to deconvolute the unknown LDOS (e.g., those based on the Tersoff approach). Using a serial expansion of the Tersoff formulas, we propose a flexible method to reconstruct the LDOS from local current-voltage characteristics measured in STM experiments. We establish a set of key physical parameters, which characterize the tunneling current of a STM-probe–sample contact and the sample LDOS expanded in Gaussian functions. Using a direct variational method coupled with probabilistic analysis, we determine these parameters from the STM experiments for MoS2 nanoflakes with different numbers of layers. The main result is the reconstruction of the LDOS in a relatively wide energy range around the Fermi level, which allows us to gain insight into the local band structure of LDTMDs. The reconstructed LDOS reveal pronounced size effects for the single-layer, two-layer, and three-layer MoS 2 nanoflakes, which we relate to the low dimensionality and strong bending or corrugation of the nanoflakes. We hope that the proposed elaboration of the Tersoff approach, allowing LDOS reconstruction, will be of critical interest for the quantitative description of STM experiments and also of use to better understand the microscopic physical aspects of the surface, strain, and bending contributions to the LDTMDs’ electronic properties.

42 ENGINEERING↗

Water and Solute Activities Regulate CO 2 Reduction in Gas-Diffusion Electrodes

Electrolysis of CO 2 at gas-diffusion electrodes (GDEs) has typically been limited by the supply of gas to the electrocatalyst, overshadowing the importance of the supply of water. However, at high current densities that approach 1 A cm –2 , where the electrolyte becomes highly concentrated in the catalyst layer of a GDE, the activity of water and solutes deviate from their bulk dilute solution values, potentially slowing reaction rates and changing reaction equilibrium potentials. In addition, as flow plates for the gas stream are introduced to enable larger electrodes and high single pass conversion of CO 2 to product, variations in the gas composition will become important. By drawing upon literature for the oxygen reduction reaction (ORR), here we explain how to account for these effects in future modeling and experimental work, with particular attention to accurate use of the Nernst equation for electrode potentials and the Arrhenius equation for reaction rates. Specifically, using measurements of KOH solvent and solute activity reported in literature, and assuming the second protonation of CO 2 by water as the rate-determining step, we show the Nernst equation dilute-solution approximation of the CO 2 to CO equilibrium potential to be accurate below 5 M KOH, but it has a 74 mV error when increasing the concentration up to 10 M KOH. Finally, a simple one-dimensional model of a serpentine flow-field on a GDE demonstrated that a reactor with constant pressure of 1 bar and 1 A cm –2 at the inlet had only ~0.3 A cm –2 at the outlet for a conversion in CO 2 partial pressure from 0.90 to 0.48 bar, showing the significant practical implications of this work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulation and Reynolds-averaged Navier-Stokes modeling of a three-component Rayleigh-Taylor mixing problem with thermonuclear burn

Rayleigh-Taylor mixing in the presence of a third component with intermediate density is investigated through three-dimensional large-eddy simulation (LES) with a high-order compact finite-difference code. Two configurations are considered: (1) a symmetric configuration in which the Atwood number between the heavy and intermediate components matches the Atwood number between the intermediate and light components and (2) an asymmetric configuration in which the Atwood number between the heavy and intermediate components is an order of magnitude greater than the Atwood number between the intermediate and light components. Mass fraction covariances are extracted, and proposed Reynolds-averaged Navier-Stokes (RANS) closures for density-specific-volume and density-mass-fraction covariances are evaluated in an a priori fashion. Additionally, a multicomponent extension of the k - Φ - L - a - V RANS model [Morgan, Phys. Rev. E 104, 015107 (2021)] is presented which includes model equations for the upper-triangular elements of the mass fraction covariance matrix. This model, referred to as the k - Φ - L - a - C model, is compared against results from LES and against other RANS models. Profiles of average mass fraction, mass-fraction covariance, and density-specific-volume covariance obtained with the k - Φ - L - a - C model are found to agree well with LES data. Finally, the impact of three-component turbulent mixing on average reaction rate is investigated in both premixed and nonpremixed cases by heating the mixing layer and allowing it to undergo thermonuclear (TN) burn. A closure model for average reaction rate is proposed for use with the k - Φ - L - a - C model, and when this model is applied, improved agreement is obtained between LES and RANS in total TN neutron production.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Angular-momentum Transport in Proto-neutron Stars and the Fate of Neutron Star Merger Remnants

Both the core collapse of rotating massive stars, and the coalescence of neutron star (NS) binaries result in the formation of a hot, differentially rotating NS remnant. The timescales over which differential rotation is removed by internal angular-momentum transport processes (viscosity) have key implications for the remnant’s long-term stability and the NS equation of state (EOS). Guided by a nonrotating model of a cooling proto-NS, we estimate the dominant sources of viscosity using an externally imposed angular-velocity profile Ω(r). Although the magneto-rotational instability provides the dominant source of effective viscosity at large radii, convection and/or the Tayler–Spruit dynamo dominate in the core of merger remnants where dΩ/dr ≥ 0. Furthermore, the viscous timescale in the remnant core is sufficiently short that solid-body rotation will be enforced faster than matter is accreted from rotationally supported outer layers. Guided by these results, we develop a toy model for how the merger remnant core grows in mass and angular momentum due to accretion. We find that merger remnants with sufficiently massive and slowly rotating initial cores may collapse to black holes via envelope accretion, even when the total remnant mass is less than the usually considered threshold ≈1.2 M TOV for forming a stable solid-body rotating NS remnant (where M TOV is the maximum nonrotating NS mass supported by the EOS). This qualitatively new picture of the post-merger remnant evolution and stability criterion has important implications for the expected electromagnetic counterparts from binary NS mergers and for multimessenger constraints on the NS EOS

79 ASTRONOMY AND ASTROPHYSICS↗