Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamic interactions”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Physics of Hard Spheres Experiment: Significant and Quantitative Findings Made

Direct examination of atomic interactions is difficult. One powerful approach to visualizing atomic interactions is to study near-index-matched colloidal dispersions of microscopic plastic spheres, which can be probed by visible light. Such spheres interact through hydrodynamic and Brownian forces, but they feel no direct force before an infinite repulsion at contact. Through the microgravity flight of the Physics of Hard Spheres Experiment (PHaSE), researchers have sought a more complete understanding of the entropically driven disorder-order transition in hard-sphere colloidal dispersions. The experiment was conceived by Professors Paul M. Chaikin and William B. Russel of Princeton University. Microgravity was required because, on Earth, index-matched colloidal dispersions often cannot be density matched, resulting in significant settling over the crystallization period. This settling makes them a poor model of the equilibrium atomic system, where the effect of gravity is truly negligible. For this purpose, a customized light-scattering instrument was designed, built, and flown by the NASA Glenn Research Center at Lewis Field on the space shuttle (shuttle missions STS 83 and STS 94). This instrument performed both static and dynamic light scattering, with sample oscillation for determining rheological properties. Scattered light from a 532- nm laser was recorded either by a 10-bit charge-coupled discharge (CCD) camera from a concentric screen covering angles of 0 to 60 or by sensitive avalanche photodiode detectors, which convert the photons into binary data from which two correlators compute autocorrelation functions. The sample cell was driven by a direct-current servomotor to allow sinusoidal oscillation for the measurement of rheological properties. Significant microgravity research findings include the observation of beautiful dendritic crystals, the crystallization of a "glassy phase" sample in microgravity that did not crystallize for over 1 year in 1g (Earth's gravity), and the emergence of face-centered-cubic (FCC) crystals late in the coarsening process (as small crystallites lost particles to the slow ripening of large crystallites). Significant quantitative findings from the microgravity experiments have been developed describing complex interactions among crystallites during the growth process, as concentration fields overlap in the surrounding disordered phase. Time-resolved Bragg scattering under microgravity captures one effect of these interactions quite conclusively for the sample at a volume fraction of 0.528. From the earliest time until the sample is almost fully crystalline, the size and overall crystallinity grow monotonically, but the number of crystallites per unit volume (number density) falls. Apparently nucleation is slower than the loss of crystallites because of the transfer of particles from small to large crystals. Thus, coarsening occurs simultaneously with growth, rather than following the completion of nucleation and growth as is generally assumed. In the same sample, an interesting signature appears in the apparent number density of crystallites and the volume fraction within the crystallites shortly before full crystallinity is reached. A brief upturn in both indicates the creation of more domains of the size of the average crystallite simultaneous with the compression of the crystallites. Only the emergence of dendritic arms offers a reasonable explanation. The arms would be "seen" by the light scattering as separate domains whose smaller radii of curvature would compress the interior phase. In fiscal year 1999, numerous papers, a doctoral dissertation, and the PHaSE final report were produced. Although this flight project has been completed, plans are in place for a follow-on colloid experiment by Chaikin and Russel that employs a light microscope within Glenn's Fluids and Combustion Facility on the International Space Station. PHaSE is providing us with a deeper understanding of the nure of phase transitions. The knowledge derived has added to the understanding of condensed matter. In addition, the burgeoning study of the dynamics of colloidal self-assembly may lead to the development of a range of photonic materials that control the desirable properties of light. Thus, applications of ordered colloidal structures include not only ultrastructure ceramics, but also photonic crystals and photothermal nanosecond light-switching devices. Industries dealing with semiconductors, electro-optics, ceramics, and composites stand to benefit from such advancements.

Doherty, Michael P.↗
DOE OSTI2021

Achieving record hot spot energies with large HDC implosions on NIF in HYBRID-E

HYBRID-E is an inertial confinement fusion implosion design that increases energy coupled to the hot spot by increasing the capsule scale in cylindrical hohlraums while operating within the current experimental limits of the National Ignition Facility. HYBRID-E reduces the hohlraum scale at a fixed capsule size compared to previous HYBRID designs, thereby increasing the hohlraum efficiency and energy coupled to the capsule, and uses the cross-beam energy transfer (CBET) to control the implosion symmetry by operating the inner (23° and 30°) and outer (44° and 50°) laser beams at different wavelengths (Δλ> 0). Small case to capsule ratio designs can suffer from insufficient drive at the waist of the hohlraum. We show that only a small amount of wavelength separation between the inner and outer beams (Δλ1-2 Å) is required to control the symmetry in low-gas-filled hohlraums (0.3 mg/cm3 He) with enough drive at the waist of the hohlraum to symmetrically drive capsules 1180 μm in outer radius. This campaign is the first to use the CBET to control the symmetry in 0.3 mg/cm3 He-filled hohlraums, the lowest gas fill density yet fielded with Δλ> 0. We find a stronger sensitivity of hot spot P2 in μm per Angstrom (40–50 μm/Å wavelength separation) than observed in high-gas-filled hohlraums and previous longer pulse designs that used a hohlraum gas fill density of 0.6 mg/cm3. There is currently no indication of transfer roll-off with increasing Δλ, indicating that even longer pulses or larger capsules could be driven using the CBET in cylindrical hohlraums. We show that the radiation flux symmetry is well controlled during the foot of the pulse, and that the entire implosion can be tuned symmetrically in the presence of the CBET in this system, with low levels of laser backscatter out of the hohlraum and low levels of hot electron production from intense laser–plasma interactions. Radiation hydrodynamic simulations can accurately represent the early shock symmetry and be used as a design tool, but cannot predict the late-time radiation flux symmetry during the peak of the pulse, and semi-empirical models are used to design the experiments. Deuterium–tritium (DT)-layered tests of 1100 μm inner radius implosions showed performance close to expectations from simulations at velocities up to ~360 km/s, and record yields at this velocity, when increasing the DT fuel layer thickness to mitigate hydrodynamic mixing of the ablator into the hot spot as a result of defects in the ablator. However, when the implosion velocity was increased, mixing due to these defects impacted performance. The ratio of measured to simulated yield for these experiments was directly correlated with the level of observed mixing. These simulations suggest that reducing the mixing, e.g., by improving the capsule defects, could result in higher performance. In addition, future experiments are planned to reduce the coast time at this scale, delay between the peak compression and the end of the laser, to increase the hot spot convergence and pressure. To reduce the coast time by several hundred ps compared to the 1100 μm inner radius implosions, HYBRID-E has also fielded 1050 μm inner radius capsules, which resulted in higher hot spot pressure and a fusion energy yield of ~170 kJ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Strangeness enhancement at its extremes: multiple (multi-)strange hadron production in pp collisions at \(\sqrt{s}=5.02\) TeV

The probability to observe a specific number of strange and multi-strange hadrons (nS), denoted as P(nS), is measured by ALICE at midrapidity (|y| < 0.5) in $$\sqrt{s}=5.02$$ TeV proton-proton (pp) collisions, dividing events into several multiplicity-density classes. Exploiting, for the first time, a technique based on counting the number of strange-particle candidates event-by-event, this measurement allows one to extend the study of strangeness production beyond the mean of the distribution. This constitutes a new test bench for production mechanisms, probing events with a large imbalance between strange and non-strange content. The analysis of a large-statistics data sample makes it possible to extract P(nS) up to a maximum nS of 7 for $${\text{K}}_{\text{S}}^{0}$$, 5 for Λ and $$\overline{\Lambda }$$, 4 for Ξ− and $${\overline{\Xi } }^{+}$$, and 2 for Ω− and $${\overline{\Omega } }^{+}$$. From this, the probability of producing strange hadron multiplets per event is calculated, thereby enabling the extension of the study of strangeness enhancement to extreme situations where several strange quarks hadronize in a single event at midrapidity. Moreover, comparing hadron combinations with different u and d quark compositions and equal overall s quark content, the contribution to the enhancement pattern coming from non-strangeness related mechanisms is isolated. The results are compared with state-of-the-art phenomenological models implemented in commonly used Monte Carlo event generators, including PYTHIA 8 Monash 2013, PYTHIA 8 with QCD-based Color Reconnection and Rope Hadronization (QCD-CR + Ropes), and EPOS LHC, which incorporates both partonic interactions and hydrodynamic evolution. These comparisons show that the new approach dramatically enhances the sensitivity to the different underlying physics mechanisms modeled by each generator.

Abualrob, I J↗
DOE OSTI2020

Nature of Non-Adiabatic Electron-Ion Forces in Liquid Metals

We report that an accurate description of electron-ion interactions in materials is crucial for our understanding of their equilibrium and non-equilibrium properties. Here, we assess the properties of frictional forces experienced by ions in non-crystalline metallic systems, including liquid metals and warm dense plasmas, that arise from electronic excitations driven by the nuclear motion due to the presence of a continuum of low-lying electronic states. To this end, we perform detailed ab-initio calculations of the full friction tensor that characterizes the set of friction forces. The nonadiabatic electron-ion interactions introduce hydrodynamic couplings between the ionic degrees of freedom, which are sizeable between nearest neighbors. The friction tensor is generally inhomogeneous, anisotropic and non-diagonal, especially at lower densities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Antielectrophoretic Response-Driven Bending–Tilting Deformation of Cationic Polyelectrolyte Brushes Drives Nonlinear Electroosmotic Transport in Brush-Grafted Nanochannels

In this paper, we use all-atom molecular dynamics (MD) simulations to describe a non-linearly enhanced electroosmotic (EOS) flow, where, in a nanochannel grafted with cationic PMETAC ([Poly(2-(Methacryloyloxy)Ethyl) Trimethylammonium Chloride]) brushes, a two-fold increase in the electric field strength leads to a several-fold (more than two-fold) increase in the EOS flow strength and volume flow rate. The electric field enforces the PMETAC brushes to undergo a bending-tilting driven deformation with a significant portion of the brush layer becoming parallel to the grafting surface. In response, a substantial fraction of the counterions leave the brush layer (hence become more mobile), but instead of going into the bulk, accumulate at the brush-bulk interface, i.e., stay in proximity of the brush segments aligned parallel to the grafting surface. This creates an interesting situation, where the counterions are not completely within the brush layer, yet they fully screen the brush charges. Such “freer” conditions enable the counterions to achieve very high velocity, thereby ensuring that the water solvating the counterions themselves move very fast triggering the significantly augmented EOS transport. Probing deeper we can identify that the bending-tilting driven brush deformation, enforcing the brushes to align parallel to the substrate, results from the anti-electrophoretic behavior of the brushes, where despite being positively charged, the brushes move against the electric field direction. Such an anti-electrophoretic behavior of the PE brushes, which has not been reported before, can be associated with the very fast velocities of the negatively charged counterions and the electrostatic and hydrodynamic coupling of the counterions with the positive functional groups of the brushes. Here, we anticipate that the findings of this paper will shed light on strategies for nanochannel flow, the anti-electrophoretic response of charged polymer chains, and the significance of capturing the detailed chemical architecture of polyelectrolytes in nanoscale science and engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2024

On the statistical theory of self-gravitating collisionless dark matter flow: Scale and redshift variation of velocity and density distributions

The statistics of velocity and density fields are crucial for cosmic structure formation and evolution. Here, this paper extends our previous work on the two-point second-order statistics for the velocity field [Phys. Fluids 35, 077105 (2023)] to one-point probability distributions for both density and velocity fields. The scale and redshift variation of density and velocity distributions are studied by a halo-based non-projection approach. First, all particles are divided into halo and out-of-halo particles so that the redshift variation can be studied via generalized kurtosis of distributions for halo and out-of-halo particles, respectively. Second, without projecting particle fields onto a structured grid, the scale variation is analyzed by identifying all particle pairs on different scales $r$. We demonstrate that: (i) Delaunay tessellation can be used to reconstruct the density field. The density correlation, spectrum, and dispersion functions were obtained, modeled, and compared with the N-body simulation; (ii) the velocity distributions are symmetric on both small and large scales and are non-symmetric with a negative skewness on intermediate scales due to the inverse energy cascade on small scales with a constant rate $\varepsilon_u$; (iii) On small scales, the even order moments of pairwise velocity $\Delta u_L$ follow a two-thirds law $\propto{(-\varepsilon_ur)}^{2/3}$, while the odd order moments follow a linear scaling $\langle(\Delta u_L)^{2n+1}\rangle=(2n+1)\langle(\Delta u_L)^{2n}\rangle\langle\Delta u_L\rangle\propto{r}$; (iv) The scale variation of the velocity distributions was studied for longitudinal velocities $u_L$ or $u_L^{'}$, pairwise velocity (velocity difference) $\Delta u_L$=$u_L^{'}$-$u_L$ and velocity sum $\Sigma u_L$=$u^{'}_L$+$u_L$. Fully developed velocity fields are never Gaussian on any scale, despite that they can initially be Gaussian; (v) On small scales, $u_L$ and $\Sigma u_L$ can be modeled by a $X$ distribution to maximize the entropy of the system. The distribution of $\Delta u_L$ can be different; (vi) On large scales, $\Delta u_L$ and $\Sigma u_L$ can be modeled by a logistic or a $X$ distribution, while $u_L$ has a different distribution; (vii) the redshift variation of the velocity distributions follows the evolution of the $X$ distribution involving a shape parameter $\alpha(z)$ decreasing with time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A review of five decades of strategic decisions and key results on the path to ignition on the National Ignition Facility

Successful demonstration of an implosion that meets the physics metrics for ignition on the National Ignition Facility (NIF) in August 2021 and demonstration of G > 1 in December 2022 represented the culmination of five decades of research. This review covers the key strategic decisions and physics results from the 20 J Janus laser to the 2 MJ NIF laser. NIF's specifications were set in order to have a margin in performance to account for uncertainty in the physics challenges leading to ignition, particularly for laser–plasma interaction and hydrodynamic instabilities including long-wavelength modes responsible for implosion symmetry. The ignition experiments utilized all that margin. Achieving the laser and target performance utilized in the Hybrid-E experiments that achieved ignition and G > 1 proved to be a challenge that consumed a decade.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Checking nonflow assumptions and results via PHENIX published correlations in p+p, p+Au, d+Au, and 3 He+Au at √ s NN=200 GeV

Recently the PHENIX Collaboration has made available two-particle correlation Fourier coefficients for multiple detector combinations in minimum bias p+p and 0–5% central p+Au, d+Au, and 3 He+Au collisions at √ s NN = 200 GeV [Phys. Rev. C 105, 024901 (2022)]. Using these coefficients for three sets of two-particle correlations, azimuthal anisotropy coefficients v 2 and v 3 are extracted for midrapidity charged hadrons as a function of transverse momentum. In this paper, we use the available coefficients to explore various nonflow hypotheses as well as to compare the results with theoretical model calculations. The nonflow methods fail basic closure tests with ampt and pythia/angantyr, particularly when including correlations with particles in the low multiplicity light-projectile going direction. In data, the nonflow adjusted v 2 results are modestly lower in p+Au and the adjusted v 3 results are more significantly higher in p+Au and d+Au. However, the resulting higher values for the ratio v 3 /v 2 in p+Au at RHIC compared to p+Pb at the LHC is additional evidence for a significant overcorrection. Incorporating these additional checks, the conclusion that these flow coefficients are dominated by initial geometry coupled with final-state interactions (e.g., hydrodynamic expansion of quark-gluon plasma) remains true, and explanations based solely on initial-state glasma are ruled out. The detailed balance between intrinsic and fluctuation-driven geometry and the exact role of weakly versus strongly coupled prehydrodynamic evolution remains an open question for triangular flow, requiring further theoretical and experimental investigation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

MFIX DEM Enhancement for Industry-Relevant Flows (Final Report)

The overall goal of this two-phase project is to implement performance improvements of the Multiphase Flow with Interphase Exchanges (MFIX) Discrete Element Model (DEM) code that enable a transformative shift for industrial use. Prior to this effort, the largest simulations performed using MFIX are O(10 7 ) particles. This falls short of the O(10 9 ) particle simulations that must be completed on a timescale of days or weeks (vs. months or years) to enable simulations with physically-relevant domain sizes to be incorporated into industrial design cycles within five years. This was accomplished by tailoring best-in-class practices to bear on the unique challenges posed by the MFIX-DEM algorithm and code base. Scientific simulations (e.g., in cosmology, turbulent combustion) routinely use massively parallel computing to update far more particles in short wall clock times. Results from Phase 1 (1.5 years in duration) indicated significant gains in speed were possible for a wide range of benchmark cases. Moreover, a survey sent to >35 companies indicates that the timing is ideal for such an enhanced tool, with >80% of the respondents indicating that DEM is already value-added or will be within the next 5 years, and >70% of the respondents indicating that improved speed is the top computational priority. In Phase 2 (3.5 years in duration), the two major barriers that hinder industry from effectively using multiphase Computational Fluid Dynamics (CFD) to cut costs and improve performance, namely computational overhead and confidence in predictions, continued to be addressed. Regarding the former, the results from Phase 1 to guide the effort, with enhancements focused on an improved time-stepping algorithm and particle sorting. Four target problems of 1 billion particles each and increasing complexity were identified: homogeneous cooling, tumbler with continuous particle size distribution, discharge from a rectangular hopper and a cylindrical riser. Each of these were successfully simulated for relevant time scales (on order of seconds) using less than 24 hours of wall clock time. These represent the first 1-billion particle DEM simulations performed with MFIX, namely using the MFIX-Exa code. This code is currently under development at NETL in collaboration with Lawrence Berkeley National Laboratory. Regarding the second barrier on predictive uncertainty, experiments from Phase 1 (interacting nozzles - hydrodynamics only) and Phase 2 (very small-scale segregation experiments) were used to demonstrate the ability of two simplified approaches to uncertainty quantification (UQ). By limiting the number of particles, UQ based on the simplified treatment was compared to standard UQ, which was shown to have much higher computational demands. Experiments were also performed on a pilot-scale stripper unit to provide validation data for future CFD-DEM simulations and UQ.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Underground Hydrogen Storage: Transforming Subsurface Science into Sustainable Energy Solutions

As the global economy moves toward net-zero carbon emissions, large-scale energy storage becomes essential to tackle the seasonal nature of renewable sources. Underground hydrogen storage (UHS) offers a feasible solution by allowing surplus renewable energy to be transformed into hydrogen and stored in deep geological formations such as aquifers, salt caverns, or depleted reservoirs, making it available for use on demand. This study thoroughly evaluates UHS concepts, procedures, and challenges. This paper analyzes the most recent breakthroughs in UHS technology and identifies special conditions needed for its successful application, including site selection guidelines, technical and geological factors, and the significance of storage characteristics. The integrity of wells and caprock, which is important for safe and efficient storage, can be affected by the operating dynamics of the hydrogen cycle, notably the fluctuations in pressure and stress within storage formations. To evaluate its potential for broader adoption, we also examined economic elements such as cost-effectiveness and the technical practicality of large-scale storage. We also reviewed current UHS efforts and identified key knowledge gaps, primarily in the areas of hydrogen–rock interactions, geochemistry, gas migration control, microbial activities, and geomechanical stability. Resolving these technological challenges, regulatory frameworks, and environmental sustainability are essential to UHS’s long-term and extensive integration into the energy industry. This article provides a roadmap for UHS research and development, emphasizing the need for further research to fully realize the technology’s promise as a pillar of the hydrogen economy.

25 ENERGY STORAGE↗
DOE OSTI2020

Star Formation in Splash Bridges

Splash bridges are created from the direct collision of two gas-rich disk galaxies. These direct collisions can eject gas masses on the order of 10 10 M ⊙ stripped from the stellar disks of each galaxy. The Taffy Galaxy system (UGC 1294/5) is a prototypical example of a splash bridge system. CO observations of the Taffy revealed that its splash bridge contains a mass of H 2 equal to the Milky Way's H 2 mass. However, the little visible star formation occurring within the bridge highlights the need for models of direct gas-rich disk collisions. The Arp 194 system displays what may be another splash bridge resulting from the collision between two disk galaxies. The region between the two stellar disks contains bright clumps of active star formation. We aim to better understand the conditions for star formation in splash bridges by employing a Jeans criterion to determine where gravitational instabilities occur in the shocked and cooling gas of gas-rich disk collisions. The splash bridge results are obtained from our previous work using a sticky particle code and post-processed. We find that the inclination between the gas disks and the collision velocity with which the gas collides strongly affects the fraction of gas that will become gravitationally unstable. Low inclinations between gas disks produce starbursts whereas high inclinations result in steady-rate star formation. The offset of the gas disks at impact will determine how many gas elements directly collide but does not strongly affect the resulting star formation.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Theoretical spectra of circumstellar dust shells around carbon-rich asymptotic giant branch stars

Realistic modeling of circumstellar dust shells around evolved stars has to include a physical treatment of the interaction among hydrodynamics, thermodynamics, radiative transfer, chemistry and dust formation and -growth. A self-consistent solution of this problem is presented in the case of stationary, spherical symmetric dust-driven winds. The resulting shell structure and the mass-loss rate are completely determined by the three fundamental stellar parameters stellar mass M(stellar), stellar luminosity L(stellar) and effective temperature T(sub eff) and by the element abundances epsilon(sub i). A detailed calculation of the transport coefficients of the dust component by means of the particle size distribution function and the solution of the non-grey radiative transfer problem provide realistic synthetic spectra of the dust shell models. We discuss the dependence of the resulting spectra on the stellar parameters in terms of infrared two color diagrams. Application of these model calculations to the prominent infrared object IRC +10216 yields best agreement with the observed spectrum and the visibility data at maximum light for the stellar parameters M(stellar) = 0.7 solar mass, L(stellar) = 2.4 x 10(exp 4) solar luminosity, T(stellar) = 2010K and a carbon to oxygen ratio of epsilon(sub c)/epsilon(sub o) = 1.40, which corresponds to a mass-loss rate of M-dot = 8 x 10(exp -5) solar mass/yr. In this model only amorphous carbon grains are considered as the main opacity source. From this model a distance to IRC +10216 of d = 170pc is deduced. The total mass contained in the circumstellar dust shell implies and initial main sequence mass of M(sub ZAMS) greater than or = 1.3 solar mass for IRC +10216.

Winters, J. M.↗
NASA NTRS1993

Study of bubble behavior in weightlessness (effects of thermal gradient and acoustic stationary wave) (M-16)

The aim of this experiment is to understand how bubbles behave in a thermal gradient and acoustic stationary wave under microgravity. In microgravity, bubble or bubbles in a liquid will not rise upward as they do on Earth but will rest where they are formed because there exists no gravity-induced buoyancy. We are interested in how bubbles move and in the mechanisms which support the movement. We will try two ways to make bubbles migrate. The first experiment concerns behavior of bubbles in a thermal gradient. It is well known than an effect of surface tension which is masked by gravity on the ground becomes dominant in microgravity. The surface tension on the side of the bubble at a lower temperature is stronger than at a higher temperature. The bubble migrates toward the higher temperature side due to the surface tension difference. The migration speed depends on the so-called Marangoni number, which is a function of the temperature difference, the bubble diameter, liquid viscosity, and thermal diffusivity. At present, some experimental data about migration speeds in liquids with very small Marangoni numbers were obtained in space experiments, but cases of large Marangoni number are rarely obtained. In our experiment a couple of bubbles are to be injected into a cell filled with silicon oil, and the temperature gradient is to be made gradually in the cell by a heater and a cooler. We will be able to determine migration speeds in a very wide range of Marangoni numbers, as well as study interactions between the bubbles. We will observe bubble movements affected by hydrodynamical and thermal interactions, the two kinds of interactions which occur simultaneously. These observation data will be useful for analyzing the interactions as well as understanding the behavior of particles or drops in materials processing. The second experiment concerns bubble movement in an acoustic stationary wave. It is known that a bubble in a stationary wave moves toward the node or the loop according to whether its diameter is larger or smaller than that of the main resonant radius. In our experiment fine bubbles will be observed to move according to an acoustic field formed in a cylindrical cell. The existence of bubbles varies the acoustic speed, and the interactive force between bubbles will make the bubble behavior collective and complicated. This experiment will be very useful to development of bubble removable technology as well as to the understanding of bubble behavior.

Azuma, H.↗
NASA NTRS1993

Hydrodynamical simulations of star-gas interactions in the interstellar medium with an external gravitational potential

We have calculated 2D models for the global structure of the interstellar medium in a disk galaxy where stars and gas are coupled through star formation, mass loss, stellar heating of the gas, and optically thin radiative cooling. In one set of models, the two dimensions are in the plane of a disk, and for one case, the total size of the region is 1 kpc square. Using a more accurate numerical approach, we confirm the results of Chiang and Bregman (1988) that H I forms into filaments or sheets separated by hotter, more diffuse gas. Although the filaments often have denser subregions, individual clouds of neutral gas are rare. We calculated the vertical distribution of the gas in another set of models where one dimension was in the plane and the other was perpendicular to it. A filament and hot gas network still occurs, although dense neutral complexes are more common and occur near the midplane. The hot gas often cools as it rises, leading to the highest velocity material having both positive and negative velocities. We find other features that are qualitatively similar to observations; for instance, H I 'worms' are produced.

Rosen, Alexander↗
NASA NTRS1994

Hydrodynamics of rotating stars and close binary interactions: Compressible ellipsoid models

We develop a new formalism to study the dynamics of fluid polytropes in three dimensions. The stars are modeled as compressible ellipsoids, and the hydrodynamic equations are reduced to a set of ordinary differential equations for the evolution of the principal axes and other global quantities. Both viscous dissipation and the gravitational radiation reaction are incorporated. We establish the validity of our approximations and demonstrate the simplicity and power of the method by rederiving a number of known results concerning the stability and dynamical oscillations of rapidly rotating polytropes. In particular, we present a generalization to compressible fluids of Chandrasekhar's classical results for the secular and dynamical instabilities of incompressible Maclaurin spheroids. We also present several applications of our method to astrophysical problems of great current interest, such as the tidal disruption of a star by a massive black hole, the coalescence of compact binaries driven by the emission of gravitational waves, and the development of instabilities in close binary systems.

Lai, Dong↗
NASA NTRS1966

Gravitational collapse and weak interactions.

Massive star behavior during final catastrophic evolution stages, stressing effect of electron type neutrino interactions, using numerical hydrodynamics coupled with energy transfer methods

HYDRODYNAMICS↗
DOE OSTI2020

C HIMERA : A Massively Parallel Code for Core-collapse Supernova Simulations

Herein we provide a detailed description of the C HIMERA code, a code developed to model core collapse supernovae (CCSNe) in multiple spatial dimensions. The CCSN explosion mechanism remains the subject of intense research. Progress to date demonstrates that it involves a complex interplay of neutrino production, transport, and interaction in the stellar core, three-dimensional stellar core fluid dynamics and its associated instabilities, nuclear burning, and the fundamental physics of the neutrino–stellar core weak interactions and the equations of state of all stellar core constituents—particularly, the nuclear equation of state associated with core nucleons, both free and bound in nuclei. C HIMERA , by incorporating detailed neutrino transport, realistic neutrino–matter interactions, three-dimensional hydrodynamics, realistic nuclear, leptonic, and photonic equations of state, and a nuclear reaction network, along with other refinements, can be used to study the role of neutrino radiation, hydrodynamic instabilities, and a variety of input physics in the explosion mechanism itself. It can also be used to compute observables such as neutrino signatures, gravitational radiation, and the products of nucleosynthesis associated with CCSNe. The code contains modules for neutrino transport, multidimensional compressible hydrodynamics, nuclear reactions, a variety of neutrino interactions, equations of state, and modules to provide data for post-processing observables such as the products of nucleosynthesis, and gravitational radiation. C HIMERA is an evolving code, being updated periodically with improved input physics and numerical refinements. We detail here the current version of the code, from which future improvements will stem, which can in turn be described as needed in future publications.

79 ASTRONOMY AND ASTROPHYSICS↗