Massive stars in quasi-static equilibrium.
Quasi-static equilibrium in massive stars during gravitational contraction including discussion of general relativistic and classical hydrodynamic results
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.
Quasi-static equilibrium in massive stars during gravitational contraction including discussion of general relativistic and classical hydrodynamic results
The collisional and dynamical processes in moon and planet formation are discussed. A hydrodynamic code of collision calculations, the orbital element changes due to gravitational scattering, a validation of the mass shifting algorithm, a theory of rotations, and the origin of asteroids are studied. A numerical model of planet growth is discussed and a methodology to evaluate the rate at which megaregolith increases its depth as a function of total accumulate number of impacts on an initially smooth, coherent surface is described.
This is the final progress report for our Astrophysics Theory Program (NRA 97-OSS12) grant NAG5-7723. We have made considerable progress on incorporating photoionization calculations with a 2.5D hydrodynamical code to model disk winds in AGNs. Following up on our simultaneous broad band monitoring campaign of the type I Seyfert galaxy NGC 5548, we have investigated the constraints imposed on models of accretion in Seyfert galaxies by their optical, UV, and X-ray spectral energy distributions (SEDs). Using results from thermal Comptonization models that relate the physical properties of the hot inner accretion flow to the thermal reprocessing that occurs in the surrounding colder thin disk, we find that we can constrain the central black hole mass, accretion rate and size scale of the hot central flow. We have applied our model to observations of Seyfert galaxies NGC 3516, NGC 7469 and NGC 5548. Our mass and accretion rate estimates for these objects roughly agree with those found using other methods.
Correlated many-fermion systems emerge in a broad range of phenomena in warm dense matter, plasmonics, and ultracold atoms. Quantum hydrodynamics (QHD) complements first-principles methods for many-fermion systems at larger scales. We illustrate the failure of the standard Bohm potential central to QHD for strong perturbations when the density perturbation is larger than about 10^{-3} 10 − 3 of the mean density. We then extend QHD to this regime via the many-fermion Bohm potential from first-principles. This may lead to more accurate QHD simulations beyond their common application domain in the presence of strong perturbations at scales unattainable with first-principles methods.
We have developed a method to extract density fluctuation measurements from x-ray radiographs of high-energy density (HED) instability growth and turbulence experiments. We use this information to calculate density fluctuation statistics for constraining the performance of turbulent mix models in HED systems. The density calculation combines image filtering, removal of systemic effects such as backlighter variation, calculation of transmission across multiple materials, and use of tracer materials to generate an approximate single-material density field. From the density map, we calculate both average density and a variance-like moment b (density-specific-volume covariance), which we compare to our models. We infer both quantities from a single image, which is significantly more information than the historic single scalar mix width measurements. We also develop a method of analyzing simulation outputs that incorporate both the density fluctuation metric from a turbulence model and the bulk material maps from the hydrodynamic code. This analysis helps address the question of how to initialize the simulations for best comparison to data from systems with large separations of scale in the mixing perturbation initial condition. We find that our data analysis method yields 1D average density and b curves with similar morphology and amplitudes as those from preliminary simulation comparisons.
A brief status report on bearing technology and present and near-term future problems that warrant research support is presented. For rolling element bearings a material with improved fracture toughness, life data in the low Lambda region, a comprehensive failure theory verified by life data and incorporated into dynamic analyses, and an improved corrosion resistant alloy are perceived as important needs. For hydrodynamic bearings better definition of cavitation boundaries and pressure distributions for squeeze film dampers, and geometry optimization for minimum power loss in turbulent film bearings are needed. For gas film bearings, foil bearing geometries that form more nearly optimum film shapes for maximum load capacity, and more effective surface protective coatings for high temperature operation are needed.
If the mass density of supersonic, collimated material is less than that of the surrounding medium, a so-called light jet will be enveloped by a cocoon of overpressured shocked gas. Hydrodynamical simulations are used to understand the evolution of the cocoon. The cocoon's evolution is also compared to a simple analytic theory. To reconcile the theory with the simulations, the growth of the jet head must be taken into account. The overpressured cocoon stage exists for a relatively short astronomical time, after which only the region of the cocoon near the jet head remains overpressured. The spatial distribution of the optical emission often observed in distant extragalactic jet systems can be explained with this improved understanding of cocoon evolution.
The author's current simulations of Giant Impacts on the protoearth show the development of large hot rock vapor atmospheres. The Balbus-Hawley mechanism will pump mass and angular momentum outwards in the equatorial plane; upon cooling and expansion the rock vapor will condense refractory material beyond the Roche distance, where it is available for lunar formation. During the last seven years, the author together with several colleagues has carried out a series of numerical investigations of the Giant Impact theory for the origin of the Moon. These involved three-dimensional simulations of the impact and its aftermath using Smooth Particle Hydrodynamics (SPH), in which the matter in the system is divided into discrete particles whose motions and internal energies are determined as a result of the imposed initial conditions. Densities and pressures are determined from the combined overlaps of the particles, which have a bell-shaped density distribution characterized by a smoothing length. In the original series of runs all particle masses and smoothing lengths had the same values; the matter in the colliding bodies consisted of initial iron cores and rock (dunite) mantles. Each of 41 runs used 3,008 particles, took several weeks of continuous computation, and gave fairly good representations of the ultimate state of the post-collision body or bodies but at best crude and qualitative information about individual particles in orbit. During the last two years an improved SPH program was used in which the masses and smoothing lengths of the particles are variable, and the intent of the current series of computations is to investigate the behavior of the matter exterior to the main parts of the body or bodies subsequent to the collisions. These runs are taking times comparable to a year of continuous computation in each case; they use 10,000 particles with 5,000 particles in the target and 5,000 in the impactor, and the particles thus have variable masses and smoothing lengths (the latter are dynamically adjusted so that a particle typically overlaps a few tens of its neighbors). Since the matter in the impactor provides the majority of the mass left in orbit after the collision, and since the masses of the particles that originated in the impactor are smaller than those in the target, the mass resolution in the exterior parts of the problem is greatly improved and the exterior particles properly simulate atmospheres in hydrostatic equilibrium.
A NOAA P-3 instrumented aircraft observed an intense, fast-moving narrow cold frontal Farmhand as it approached the Pacific Northwest coast on 19 February 2001 during the Pacific Coastal Jets Experiment. Pseudo-dual-Doppler analyses performed on the airborne Doppler radar data while the frontal system was well offshore indicated that a narrow ribbon of very high radar reflectively convective cores characterized the Farmhand at low levels with echo tops to approximately 4-5 km. The NCFR exhibited gaps in its narrow ribbon of high reflectively, probably as a result of hydrodynamic instability all no its advancing cold pool leading edge. In contrast to some earlier studies of cold frontal rainbands, density current theory described well the motion of the overall front. The character of the updraft structure associated with the heavy rainfall at its leading edge varied across the gap region. The vertical shear of the cross-frontal low-level ambient flow exerted a strong influence on the updraft character, consistent with theoretical arguments developed for squall lines describing the balance of vorticity at the leading edge. In short regions south of the gaps the vertical wind shear was strongest with the updrafts and rain shafts more intense, narrower, and more erect or even downshear tilted. North of the gaps the wind shear weakened with less intense Dihedrals which tilted upshear with a broader band of rainfall. Simulations using a nonhydrostatic mesoscale nested grid model are used to investigate the gap regions, particularly the balance of cold pool induced to pre-frontal ambient shears at the leading edge. Observations confirm the model results that the updraft character depends on the balance of vorticity at the leading edge. Downshear-tilted updrafts imply that convection south of the gap regions would weaken with time relative to the frontal segments north of the gaps since inflow air would be affected by passage through the heavy rain region before ascent, suggesting a mechanism for gap filling.
Recent advances in cosmological observations have provided an unprecedented opportunity to investigate the distribution of baryons relative to the underlying matter. In this work, we show that the gas is more extended than the dark matter, and the amount of baryonic feedback at $z \lesssim 1$ disfavors low-feedback models such as that of state-of-the-art hydrodynamical simulation IllustrisTNG compared with high-feedback models such as that of the original Illustris simulation. This has important implications for bridging the gap between theory and observations and understanding galaxy formation and evolution. Furthermore, a better grasp of the baryon-dark matter link is critical to future cosmological analyses, which are currently impeded by our limited knowledge of baryonic feedback. Here, we measure the kinematic Sunyaev-Zel'dovich (kSZ) effect from the Atacama Cosmology Telescope (ACT), stacked on the luminous red galaxy (LRG) sample of the Dark Energy Spectroscopic Instrument (DESI) imaging survey. This is the first analysis to use photometric redshifts for reconstructing galaxy velocities. Due to the large number of galaxies comprising the DESI imaging survey, this is the highest signal-to-noise stacked kSZ measurement to date: we detect the signal at 13$σ$, finding strong evidence that the gas is more spread out than the dark matter, as well as a preference for larger feedback compared to some commonly used state-of-the-art hydrodynamical simulations. Here, our work opens up the possibility of recalibrating large hydrodynamical simulations using the kSZ effect. In addition, our findings highlight the importance of properly accounting for baryonic feedback with future surveys such as LSST through direct probes such as the kSZ, and shed light on long-standing enigmas in astrophysics, such as the “missing baryon” problem.
This workshop presentation discusses the design and implementation of numerical methods for the quantification of statistical uncertainty, including a-posteriori error bounds, for output quantities computed using CFD methods. Hydrodynamic realizations often contain numerical error arising from finite-dimensional approximation (e.g. numerical methods using grids, basis functions, particles) and statistical uncertainty arising from incomplete information and/or statistical characterization of model parameters and random fields. The first task at hand is to derive formal error bounds for statistics given realizations containing finite-dimensional numerical error [1]. The error in computed output statistics contains contributions from both realization error and the error resulting from the calculation of statistics integrals using a numerical method. A second task is to devise computable a-posteriori error bounds by numerically approximating all terms arising in the error bound estimates. For the same reason that CFD calculations including error bounds but omitting uncertainty modeling are only of limited value, CFD calculations including uncertainty modeling but omitting error bounds are only of limited value. To gain maximum value from CFD calculations, a general software package for uncertainty quantification with quantified error bounds has been developed at NASA. The package provides implementations for a suite of numerical methods used in uncertainty quantification: Dense tensorization basis methods [3] and a subscale recovery variant [1] for non-smooth data, Sparse tensorization methods[2] utilizing node-nested hierarchies, Sampling methods[4] for high-dimensional random variable spaces.
Here we introduce fluctuating hydrodynamics approaches on surfaces for capturing the drift-diffusion dynamics of particles and microstructures immersed within curved fluid interfaces of spherical shape. We take into account the interfacial hydrodynamic coupling, traction coupling with the surrounding bulk fluid, and thermal fluctuations. For fluid-structure interactions, we introduce Immersed Boundary Methods (IBM) and related Stochastic Eulerian-Lagrangian Methods (SELM) for curved surfaces. We use these approaches to investigate the statistics of surface fluctuating hydrodynamics and microstructures. For velocity autocorrelations, we find characteristic power-law scalings $τ^{-1}, τ^{-2}$, and plateaus can emerge. This depends on the physical regime associated with the geometry, surface viscosity, and bulk viscosity. This differs from the characteristic $τ^{-3/2}$ scaling for bulk three dimensional fluids. We develop theory explaining these observed power-laws associated with time-scales for dissipation within the fluid interface and coupling to the surrounding fluid. We then use our introduced methods to investigate a few example systems and roles of hydrodynamic coupling and thermal fluctuations including for the kinetics of passive particles and active microswimmers in curved fluid interfaces.
We numerically investigate the internal evolution of multiphase clouds, which are at rest with respect to an ambient, highly ionized medium (HIM) representing the hot component of the circumgalactic medium. Time-dependent saturated thermal conduction and its implications like condensation rates and mixing efficiency are assessed in multiphase clouds. Our simulations are carried out by using the adaptive mesh refinement code FLASH. The model clouds are initially in both hydrostatic and thermal equilibrium and are in pressure balance with the HIM. Thus, they have steep gradients in both temperature and density at the interface to HIM leading to non-negligible thermal conduction. Several physical processes are considered numerically or semi-analytically: thermal conduction, radiative cooling and external heating of gas, self-gravity, mass diffusion, and dissociation of molecules and ionization of atoms. It turns out that saturated thermal conduction triggers a continuous condensation irrespective of cloud mass. Dynamical interactions with ambient HIM all relate to the radial density gradient in the clouds: (1) mass flux due to condensation is the higher the more homogeneous the clouds are; (2) mixing of condensed gas with cloud gas is easier in low-mass clouds, because of their shallower radial density gradient; and thus (3) accreted gas is distributed more efficiently. A distinct and sub-structured transition zone forms at the interface between cloud and HIM, which starts at smaller radii and is much narrower as deduced from analytical theory.
In this work, we show that the correlation between the elliptic momentum anisotropy v 2 and the average transverse momentum [ p T ] at fixed multiplicity in small system nuclear collisions carries information on the origin of the observed momentum anisotropy. A calculation using a hybrid IP-Glasma+ music +UrQMD model that includes contributions from final state response to the initial geometry as well as initial state momentum anisotropies of the color glass condensate predicts a characteristic sign change of the correlator ρ ^ ( v 2 2 , [ p T ] ) as a function of charged particle multiplicity in p + Au and d + Au collisions at s = 200 GeV , and p + Pb collisions at s = 5.02 TeV . This sign change is absent in calculations without initial state momentum anisotropies. The model further predicts a qualitative difference between the centrality dependence of ρ ^ ( v 2 2 , [ p T ] ) in Au + Au collisions at s = 200 GeV and Pb + Pb collisions at s = 5.02 TeV , with only the latter showing a sign change in peripheral events. Predictions for O + O collisions at different collision energy show a similar behavior. Experimental observation of these distinct qualitative features of ρ ^ ( v 2 2 , [ p T ] ) in small and large systems would constitute strong evidence for the presence and importance of initial state momentum anisotropies predicted by the color glass condensate effective theory.
Here, in a number of reported instances, implosions utilizing fuel mixtures have resulted in anomalously low fusion yields below those predicted by radiation-hydrodynamics simulations. Inter-species ion diffusion has been suggested as a possible cause of the observed yield degradation in fuel mixture implosions. An experimental platform utilizing hydro-equivalent deuterium–tritium (DT), deuterium–tritium–hydrogen (DTH), and deuterium-tritium-helium3 (DT 3 He) capsule implosions was developed to determine whether the inter-species ion diffusion theory may describe the resulting fuel mixture implosion behavior. The implosion experiments were performed at the Omega laser facility. X-ray images and shell areal density diagnostics results show that the hydro-equivalent three capsules (DT, DTH, and DT 3 He) have similar compression behavior. However, nuclear yield deviation was observed from the scaling determined using a fusion yield formula. In the DT 3 He mixture, a reduced yield of a factor of 0.65 ± 0.13 was observed, which is similar to a yield reduction observed in D 3 He mixture by Rygg et al. (i.e., Rygg effect). In contrast, in the DTH mixture, a factor of 1.17 ±0.15 yield increase was observed, which we named the inverse Rygg effect. The yield increase observed in the DTH mixture is consistent with the inter-species ion diffusion theory where lighter H diffuses away from the core and concentrated DT in the core produces higher yield. An inter-species ion diffusion model, the Zimmerman–Paquette–Kagan–Zhdanov model, implemented in a Lagrangian radiation-hydrodynamics fluid code, was also used to analyze the present data, without the need to assume hydrodynamic equivalence of the capsules, but it does not completely explain the DTH or DT 3 He capsules although its effects are in the correct direction. Simulation-based Bayesian inference was used in the latter analysis to quantify the uncertainty in the numerical simulations. The simulation-based analysis resulted in an inferred Rygg-effect yield decrease factor of 0.91 ± 0.02 for the DT 3 He mixture, and an inferred inverse-Rygg yield increase factor of 1.21 ± 0.04 for the DTH mixture, based on simulations ignoring ion diffusion.
I review recent developments in the phenomenological study of the quark-gluon plasma (QGP) transport properties based on a personal selection of results that were presented at Quark Matter 2019. Additionally, the constraints on the temperature dependence of QGP shear and bulk viscosity are summarized. I discuss new theory advancements towards more realistic 3D dynamical simulations of heavy-ion collisions at finite baryon density. The challenges and opportunities of applying hydrodynamics to small collision systems are highlighted.
Pumps used in large liquid-fuel rocket engines are examined. The term preliminary design denotes the initial, creative phases of design, where the general shape and characteristics of the machine are determined. This compendium is intended to provide the design engineer responsible for these initial phases with a physical understanding and background knowledge of the numerous special fields involved in the design process. Primary attention is directed to the pumping part of the turbopump and hence is concerned with essentially incompressible fluids. However, compressible flow principles are developed. As much as possible, the simplicity and reliability of incompressible flow considerations are retained by treating the mechanics of compressible fluids as a departure from the theory of incompressible fluids. Five areas are discussed: a survey of the field of turbomachinery in dimensionless form; the theoretical principles of the hydrodynamic design of turbomachinery; the hydrodynamic and gas dynamic design of axial flow turbomachinery; the hydrodynamic and gas dynamic design of radial and mixed flow turbomachinery; and some mechanical design considerations of turbomachinery. Theoretical considerations are presented with a relatively elementary mathematical treatment.
Planar black holes in AdS have long-lived quasinormal modes which capture the physics of charge and momentum diffusion in the dual field theory. How should we characterize the effective dynamics of a probe system coupled to the conserved currents of the dual field theory? Specifically, how would such a probe record the long-lived memory of the black hole and its Hawking fluctuations? We address this question by exhibiting a universal gauge invariant framework which captures the physics of stochastic diffusion in holography: a designer scalar with a gravitational coupling governed by a single parameter, the Markovianity index. We argue that the physics of gauge and gravitational perturbations of a planar Schwarzschild-AdS black hole can be efficiently captured by such designer scalars. We demonstrate that this framework allows one to decouple, at the quadratic order, the long-lived quasinormal and Hawking modes from the short-lived ones. It furthermore provides a template for analyzing fluctuating open quantum field theories with memory. In particular, we use this set-up to analyze the diffusive Hawking photons and gravitons about a planar Schwarzschild-AdS black hole and derive the quadratic effective action that governs fluctuating hydrodynamics of the dual CFT. Along the way we also derive results relevant for probes of hyperscaling violating backgrounds at finite temperature.