SEARCH · Engineering Papers
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.
In continuous flow electrophoresis the axial flow structure changes from a fully developed rectilinear form to one characterized by meandering as power levels are increased. The origin of this meandering is postulated to lie in a hydrodynamic instability driven by axial (and possibly lateral) temperature gradients. Experiments done at MSFC show agreement with the theory.
This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main-sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses(MM0.1510)and six black hole masses( MM105 105BH7). We present here the results most relevant to observations. The effects of internal stellar structure and relativity decouple for both the disruption cross section and the characteristic energy width of the debris. Moreover, the full disruption cross section is almost independent of Må for Må/Me3. Independent of Må, relativistic effects increase the critical pericenter distance for full disruption events by up to a factor of∼3 relative to the Newtonian prediction. The probability of a direct capture is also independent of Må; at MBH/Me;5×106 this probability is equal to the probability of a complete disruption. The breadth of the debris energy distribution ΔE can differ from the standard estimate by factors of0.35−2, depending on Må and MBH, implying a corresponding change(∝(ΔE)−3/2)in the characteristic mass-return timescale. We provide analytic forms, suitable for use in both event rate estimates and parameter inference, to describe all these trends. For partial disruptions, we find a nearly universal relation between the star’s angular momentum and the fraction of Må remaining. Within the “empty loss-cone” regime, partial disruptions must precede full disruptions. These partial disruptions can drastically affect the rate and appearance of subsequent total disruptions.
This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main-sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses (0.15 < M_/M<10)and six black hole masses( ́10<MM105 105BH7). We present here the results most relevant to observations. The effects of internal stellar structure and relativity decouple for both the disruption cross section and the characteristic energy width of the debris. Moreover, the full disruption cross section is almost independent of Må for Må/Me3. Independent of Må, relativistic effects increase the critical pericenter distance for full disruption events by up to a factor of ∼3 relative to the Newtonian prediction. The probability of a direct capture is also independent of Må; at MBH/Me5×106this probability is equal to the probability of a complete disruption. The breadth of the debris energy distribution ΔE can differ from the standard estimate by factors of 0.35−2, depending on Må and MBH, implying a corresponding change(∝(ΔE)−3/2)in the characteristic mass-return timescale. We provide analytic forms, suitable for use in both event rate estimates and parameter inference, to describe all these trends. For partial disruptions, we find a nearly universal relation between the star’s angular momentum and the fraction of Må remaining. Within the “empty loss-cone” regime, partial disruptions must precede full disruptions. These partial disruptions can drastically affect the rate and appearance of subsequent total disruptions
In this LDRD we investigated the application of machine learning methods to understand dimensionality reduction and evolution of the Rayleigh-Taylor instability (RTI). As part of the project, we undertook a significant literature review to understand current analytical theory and machine learning based methods to treat evolution of this instability. We note that we chose to refocus on assessing the hydrodynamic RTI as opposed to the magneto-Rayleigh-Taylor instability originally proposed. This choice enabled utilizing a wealth of analytic test cases and working with relatively fast running open-source simulations of single-mode RTI. This greatly facilitated external collaboration with URA summer fellowship student, Theodore Broeren. In this project we studied the application of methods from dynamical systems learning and traditional regression methods to recover behavior of RTI ranging from the fully nonlinear to weakly nonlinear (wNL) regimes. Here we report on two of the tested methods SINDy and a more traditional regression-based approach inspired by analytic wNL theory with which we had the most success. We conclude with a discussion of potential future extensions to this work that may improve our understanding from both theoretical and phenomenological perspectives.
Whether or not femto-scale droplets of quark-gluon plasma (QGP) are formed in so-called small systems at high-energy colliders is a pressing question in the phenomenology of the strong interaction. For proton-proton or proton-nucleus collisions the answer is inconclusive due to the large theoretical uncertainties plaguing the description of these processes. While upcoming data on collisions of 16 O nuclei may mitigate these uncertainties in the near future, here we demonstrate the unique possibilities offered by complementing 16 O + 16 O data with collisions of 20 Ne ions. We couple both nuclear lattice effective field theory (NLEFT) and projected generator coordinate method (PGCM) ab initio descriptions of the structure of 20 Ne and 16 O to hydrodynamic simulations of 16 O + 16 O and 20 Ne + 20 Ne collisions at high energy. We isolate the imprints of the bowling-pin shape of 20 Ne on the collective flow of hadrons, which can be used to perform quantitative tests of the hydrodynamic QGP paradigm. In particular, we predict that the elliptic flow of 20 Ne + 20 Ne collisions is enhanced by as much as 1.174(8) stat (31) syst for NLEFT and 1.139(6) stat (39) syst for PGCM relative to 16 O + 16 O collisions for the 1% most central events. At the same time, theoretical uncertainties largely cancel when studying relative variations of observables between two systems. This demonstrates a method based on experiments with two light-ion species for precision characterizations of the collective dynamics and its emergence in a small system.
Comparison of hydrodynamic calculations with experimental data inevitably requires a model for converting the fluid to particles. For shear viscous fluids, several different conversion models have been proposed, such as the quadratic Grad corrections, the Strickland-Romatschke (SR) ansatz, self-consistent shear corrections from linearized kinetic theory, or corrections from the relaxation time approach. We show, via comparison to nonlinear 2 → 2 kinetic theory evolution, that even for low specific shear viscosities 0.03 ≲ $\textit{η/s}$ ≲ 0.2 , these four models have varying accuracy and still significant errors in reproducing particle phase space densities from hydrodynamic fields. Moreover, we demonstrate that the overall reconstruction error of additive shear viscous $f = f_{eq} + δf$ models dramatically reduces by up to one order of magnitude, if one ensures through exponentiation that $\textit{f}$ is always positive. Further, we also illustrate how even more accurate viscous $\textit{δf}$ models can be constructed by incorporating the first time derivatives of hydrodynamic fields to account for the past evolution of the system. Such time derivatives are readily available in hydrodynamic simulations, though usually not included in the output. Though our comparisons are limited to a one-component massless system undergoing a (0 + 1)-dimensional (0 + 1D) longitudinal boost-invariant expansion, we expect that the improved models will be useful in 2 + 1D and 3 + 1D hydrodynamic studies as well.
Recent experimental, numerical and analytical work has shown that the stoichiometric mixture fraction (Z(sub st)) can have a profound effect on soot formation in diffusion flames. These findings were obtained at constant flame temperature (T(sub ad)), employing the approach described in Du and Axelbaum (1995, 1996). For example, a fuel mixture containing 1 mole of ethylene and 11.28 moles of nitrogen burning in pure oxygen ((Z(sub st)) = 0.78) has the same adiabatic flame temperature (2370 K) as that of pure ethylene burning in air ((Z(sub st)) = 0.064). An important finding of these works was that at sufficiently high (Z(sub st)), flames remain blue as strain rate approaches zero in counterflow flames, or as flame height and residence time approach infinity in coflowing flames. Lin and Faeth (1996a) coined the term permanently blue to describe such flames. Two theories have been proposed to explain the appearance of permanently-blue flames at high (Z(sub st)). They are based on (1) hydrodynamics and (2) flame structure. Previous experimental studies in normal gravity are not definitive as to which, if either, mechanism is dominant because both hydrodynamics and structure suppress soot formation at high (Z(sub st)) in coflowing and counterflowing diffusion flames. In counterflow flames with (Z(sub st)) < 0.5 streamlines at the flame sheet are directed toward the fuel. Newly formed soot is convected into richer regions, favoring soot growth over oxidation. For (Z(sub st)) > 0.5, convection at the flame is toward the oxidizer, thus enhancing soot oxidization. Thus, in counterflow flames, hydrodynamics causes soot to be convected towards the oxidizer at high (Z(sub st)) which suppresses soot formation. Axelbaum and co-workers maintain that while the direction of convection can impact soot growth and oxidation, these processes alone cannot cause permanently-blue flames. Soot growth and oxidation are dependent on the existence of soot particles and the presence of soot is invariably accompanied by yellow luminosity. Soot-particle inception, on the other hand, arises from gas-phase reactions and its dependence on flow direction is weak, similar to that of other gas-phase reactions in flames. For example, when the flame moves across the stagnation plane no significant changes in flame chemistry are observed. Furthermore, since the soot-inception zone has a finite thickness, soot has been produced in counterflow flames with (Z(sub st)) > 0.5. For large (Z(sub st)) the fuel concentration decreases and oxygen concentration increases in the soot forming regions of the flame. This yields a shift in the OH profile toward the fuel side of the flame, and this shift can dramatically influence soot inception because it essentially narrows the soot inception zone. Soot-free (permanently-blue) conditions can be realized when the structure of the flame is adjusted to the extent that significant oxidizing species exist on the fuel side of the flame at temperatures above the critical temperature for soot inception, ca. 1250 K. In previously considered flames it was impossible to independently vary flame structure and convection direction. In contrast, spherical diffusion flames (which generally require microgravity) allow both properties to be varied independently. We altered structure (Z(sub st)) by exchanging inert between the oxidizer and the fuel and we independently varied convection direction at the flame sheet by interchanging the injected and ambient gases. In this work we established four flames: (a) ethylene issuing into air, (b) diluted ethylene issuing into oxygen, (c) air issuing into ethylene, and (d) oxygen issuing into diluted ethylene. (Z(sub st)) is 0.064 in flames (a) and (c) and 0.78 in flames (b) and (d). The convection direction is from fuel to oxidizer in flames (a) and (b) and from oxidizer to fuel in flames (c) and (d). Under the assumption of equal diffusivities of all species and heat, the stoichiometric contours of these flames have identical temperatures and nitrogen concentrations.
ABSTRACT We investigate the 3D spin alignment of galaxies with respect to the large-scale filaments using the MaNGA survey. The cosmic web is reconstructed from the Sloan Digital Sky Survey using disperse and the 3D spins of MaNGA galaxies are estimated using the thin disc approximation with integral field spectroscopy kinematics. Late-type spiral galaxies are found to have their spins parallel to the closest filament’s axis. The alignment signal is found to be dominated by low-mass spirals. Spins of S0-type galaxies tend to be oriented preferentially in perpendicular direction with respect to the filament’s axis. This orthogonal orientation is found to be dominated by S0s that show a notable misalignment between their kinematic components of stellar and ionized gas velocity fields and/or by low-mass S0s with lower rotation support compared to their high-mass counterparts. Qualitatively similar results are obtained when splitting galaxies based on the degree of ordered stellar rotation, such that galaxies with high spin magnitude have their spin aligned, and those with low spin magnitude in perpendicular direction to the filaments. In the context of conditional tidal torque theory, these findings suggest that galaxies’ spins retain memory of their larger scale environment. In agreement with measurements from hydrodynamical cosmological simulations, the measured signal at low redshift is weak, yet statistically significant. The dependence of the spin-filament orientation of galaxies on their stellar mass, morphology, and kinematics highlights the importance of sample selection to detect the signal.
A new theory is proposed to explain the heating of the solar chromosphere, and possibly the corona, by the dissipation of hydrodynamic compression waves. The basis of the dissipative mechanism, here referred to as ionization pumping, is hysteresis caused by irreversible relaxation of the chromospheric medium to ionization equilibrium following pressure perturbations. In the middle chromosphere, where hydrogen is partially ionized, it is shown that ionization pumping will cause strong dissipation of waves whose periods are 200s or less. This could cause heating of the chromosphere sufficient to compensate for the radiative losses. The mechanism retains a high efficiency for waves of arbitrarily small amplitude and, thus, can be more efficient than shock dissipation for small perturbations in pressure. The formation of shocks therefore is not required for the dissipation of waves whose periods are several minutes or less.
The quark–gluon plasma (QGP) produced by collisions between ultrarelativistic heavy nuclei is well described in the language of hydrodynamics. Noncentral collisions are characterized by very large angular momentum, which in a fluid system manifests as flow vorticity. This rotational structure can lead to a spin polarization of the hadrons that eventually emerge from the plasma, and thus these collisions provide experimental access to flow substructure at unprecedented detail. Recently, the first observations of Λ hyperon polarization along the direction of collisional angular momentum were reported. These measurements are in broad agreement with hydrodynamic and transport-based calculations and reveal that the QGP is the most vortical fluid ever observed. However, there remain important tensions between theory and observation that might be fundamental in nature. In the relatively mature field of heavy-ion physics, the discovery of global hyperon polarization and 3D simulations of the collision have opened an entirely new direction of research. We discuss the current status of this rapidly developing area and directions for future research.
We investigate the transport properties of a kinetic theory model that is tuned to describe the thermodynamic properties of QCD at zero chemical potential using a new formulation of the relaxation time approximation. In contrast to previous approaches, the latter is constructed to preserve the fundamental properties of the collision term of the Boltzmann equation for any energy dependence of the relaxation time. A novel choice of matching conditions is implemented to ensure that the background mean-field depends only on the temperature even when the system is out of equilibrium. We provide a consistent analysis of how the transport coefficients of relativistic Navier-Stokes theory vary with the energy dependence of the relaxation time. We also show that the entropy production of this theory is consistent with the second law of thermodynamics and verify that it is independent of the matching conditions employed. We used this fact to calculate the matching independent combination of transport coefficients.
Recent data from SAS-2 on the galactic gamma-ray line flux as a function of longitude reveal a broad maximum in the region below 30 deg. These data, as unfolded here, imply that the low-energy (1-10 GeV) galactic cosmic-ray flux varies with the radial distance from the galactic center and is about an order of magnitude higher than the local value in a toroidal region for radial distances between 4 and 5 kpc. We further show that this enhancement can be plausibly accounted for by Fermi acceleration and compression caused by a hydrodynamic shock driven by the expanding gas in the 3-kpc arm and invoked in some versions of galactic structure theory.
The present paper provides the results of molecular dynamics calculations on a Na(+) ion in aqueous solution. Attention is given to the sodium-oxygen and sodium-hydrogen radial distribution functions, the velocity autocorrelation function for the Na(+) ion, the autocorrelation function of the force on the stationary ion, and the accuracy of Brownian motion assumptions which are basic to hydrodynamic models of ion dyanmics in solution. It is pointed out that the presented calculations provide accurate data for testing theories of ion dynamics in solution. The conducted tests show that it is feasible to calculate Brownian friction constants for ions in aqueous solutions. It is found that for Na(+) under the considered conditions the Brownian mobility is in error by only 60 percent.
Based on the theory of the thermodynamic equilibrium in a system of quantum vortices in superfluid liquids and in a Bose – Einstein condensate in the presence of a counterflow of normal and superfluid components, we investigate the structure of a chaotic tangle of quantum vortices in turbulent superfluid liquids. Using the characteristic functional method, the properties of hydrodynamic vortex filaments are examined. It is shown that the average curvature of the vortex lines is on the order of the intervortex distance, with the proportionality coefficient being independent of the counterflow velocity. It is found that the degree of anisotropy of the vortex loops does not depend on the applied counterflow velocity. The obtained results explain the origin of the anisotropy as well as the relationship between the curvature of the lines and the intervortex space and their dependence on the parameters of the problem. (paper)
Here, in this study, we introduce and utilize GEOSPH, a Lagrangian particle-based continuum Smoothed Particle Hydrodynamics (SPH) code to simulate the deformation of geologic materials in sandbox-style experiments using an elastoplastic constitutive model. Due to the meshfree nature of SPH, our proposed framework is capable of handling the large deformations characteristic of sandbox-style experiments. SPH is also able to capture strain localization and discontinuities, allowing us to resolve faulting patterns emerging in the sandbox experiments. We successfully validate our SPH implementation against the benchmark experiments of Buiter and co-workers, and then apply SPH to understand the deformational processes occurring in doubly vergent orogens resultant from the S-point experimental setup and from finitely high wedge-like indenter backstops of different inclinations. We find that the topographic slope of the prowedge results from the competing processes of vertical uplift due to the indenter and forward thrusting, while the topographic slope of the retrowedge is a product of repeated shallow slope failure. Our findings show that SPH is a promising numerical method which should be added to the current palette of computational tools used by structural geologists for solving problems relating to tectonics and crustal deformation.
We present precision measurements of the bias parameters of the one-loop power spectrum model of the Lyman- α (Ly- α ) forest, derived within the effective field theory (EFT) of large-scale structure. We fit our model to the three-dimensional flux power spectrum measured from the ACCEL 2 hydrodynamic simulations. The EFT model fits the data with an accuracy of below 2% up to k = 2 h Mpc − 1 . Further, we analytically derive how nonlinearities in the three-dimensional clustering of the Ly- α forest introduce biases in measurements of the baryon acoustic oscillations (BAOs) scaling parameters in radial and transverse directions. From our EFT parameter measurements, we obtain a theoretical error budget of Δ α ∥ = − 0.2 % ( Δ α ⊥ = − 0.3 % ) for the radial (transverse) parameters at redshift z = 2.0 . This corresponds to a shift of − 0.3 % (0.1%) for the isotropic (anisotropic) distance measurements. We provide an estimate for the shift of the BAO peak for Ly- α -quasar cross-correlation measurements assuming analytical and simulation-based scaling relations for the nonlinear quasar bias parameters resulting in a shift of − 0.2 % ( − 0.1 % ) for the radial (transverse) dilation parameters, respectively. This analysis emphasizes the robustness of Ly- α forest BAO measurements to the theory modeling. We provide informative priors and an error budget for measuring the BAO feature—a key science driver of the currently observing Dark Energy Spectroscopic Instrument (DESI). Our work paves the way for full-shape cosmological analyses of Ly- α forest data from DESI and upcoming surveys such as the Prime Focus Spectrograph, WEAVE-QSO, and 4MOST. Published by the American Physical Society 2025
Correlations between mean transverse momentum [pT] and anisotropic flow coefficients v 2 or v 3 are measured as a function of centrality in Pb–Pb and Xe–Xe collisions at $\sqrt{s{NN}}$ = 5.02 TeV and 5.44 TeV, respectively, with ALICE. In addition, the recently proposed higher-order correlation between [p T ], v 2 , and v 3 is measured for the first time, which shows an anticorrelation for the presented centrality ranges. These measurements are compared with hydrodynamic calculations using IP-Glasma and T R ENTo initialstate shapes, the former based on the Color Glass Condensate effective theory with gluon saturation, and the latter a parameterized model with nucleons as the relevant degrees of freedom. The data are better described by the IP-Glasma rather than the T R ENTo based calculations. In particular, Trajectum and JETSCAPE predictions, both based on the T R ENTo initial state model but with different parameter settings, fail to describe the measurements. As the correlations between [p T ] and v n are mainly driven by the correlations of the size and the shape of the system in the initial state, these new studies pave a novel way to characterize the initial state and help pin down the uncertainty of the extracted properties of the quark–gluon plasma recreated in relativistic heavy-ion collisions
Two-fluid theory applied to flow stability of helium II down inclined plane