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 343 records · Page 19

A Massively Parallel Hybrid Dusty-Gasdynamics and Kinetic Direct Simulation Monte Carlo Model for Planetary Applications

In order to understand the global structure, dynamics, and physical and chemical processes occurring in the upper atmospheres, exospheres, and ionospheres of the Earth, the other planets, comets and planetary satellites and their interactions with their outer particles and fields environs, it is often necessary to address the fundamentally non-equilibrium aspects of the physical environment. These are regions where complex chemistry, energetics, and electromagnetic field influences are important. Traditional approaches are based largely on hydrodynamic or magnetohydrodynamic (MHD) formulations and are very important and highly useful. However, these methods often have limitations in rarefied physical regimes where the molecular collision rates and ion gyrofrequencies are small and where interactions with ionospheres and upper neutral atmospheres are important. At the University of Michigan we have an established base of experience and expertise in numerical simulations based on particle codes which address these physical regimes. The Principal Investigator, Dr. Michael Combi, has over 20 years of experience in the development of particle-kinetic and hybrid kinetichydrodynamics models and their direct use in data analysis. He has also worked in ground-based and space-based remote observational work and on spacecraft instrument teams. His research has involved studies of cometary atmospheres and ionospheres and their interaction with the solar wind, the neutral gas clouds escaping from Jupiter s moon Io, the interaction of the atmospheres/ionospheres of Io and Europa with Jupiter s corotating magnetosphere, as well as Earth s ionosphere. This report describes our progress during the year. The contained in section 2 of this report will serve as the basis of a paper describing the method and its application to the cometary coma that will be continued under a research and analysis grant that supports various applications of theoretical comet models to understanding the inner comae of comets (grant NAGS- 13239 from the Planetary Atmospheres program).

Combi, Michael R.↗

60 years of science in ICF: from conception to scientific breakeven on the National Ignition Facility

The recent achievements of a burning plasma, fusion ignition, and scientific energy gain with deuterium-tritium (DT) fuel at Lawrence Livermore National Laboratory’s National Ignition Facility (NIF) represents a major milestone in the development of inertial confinement fusion (ICF) and all of fusion research. In these experiments, fuel pressures well in excess of hundreds of GBars were achieved in the compressed fuel, and robust alpha heating of the fuel, far in excess of the energy provided by the implosion, were demonstrated for the first time. These achievements occurred 60 years after the inception of ICF and the first laser demonstration, and were made possible by more than five decades of research at laser facilities around the world. Advances in laser technology both in wavelength and precision, motivated by improved understanding of laser-plasma interaction physics and the demands of targets; improvements in target fabrication inspired by the need to control and minimize hydrodynamic instabilities in the implosion; and multi-dimensional simulations and diagnostics have been critical to this achievement. This paper will summarize the scientific and technical advances, the surprises, and the challenges that had to be overcome to achieve these goals.

fusion↗

A statistical mechanics approach to macroscopic limits of car-following traffic dynamics

Herein we study the derivation of macroscopic traffic models from car-following vehicle dynamics by means of hydrodynamic limits of an Enskog-type kinetic description. We consider the superposition of Follow-the-Leader (FTL) interactions and relaxation towards a traffic-dependent Optimal Velocity (OV) and we show that the resulting macroscopic models depend on the relative frequency between these two microscopic processes. If FTL interactions dominate then one gets an inhomogeneous Aw–Rascle–Zhang model, whose (pseudo) pressure and stability of the uniform flow are precisely defined by some features of the microscopic FTL and OV dynamics. Conversely, if the rate of OV relaxation is comparable to that of FTL interactions then one gets a Lighthill–Whitham–Richards model ruled only by the OV function. We further confirm these findings by means of numerical simulations of the particle system and the macroscopic models. Unlike other formally analogous results, our approach builds the macroscopic models as physical limits of particle dynamics rather than assessing the convergence of microscopic to macroscopic solutions under suitable numerical discretisations.

97 MATHEMATICS AND COMPUTING↗

Dark matter cores and cusps in spiral galaxies and their explanations

We compare proposed solutions to the core vs cusp issue of spiral galaxies, which has also been framed as a diversity problem, and demonstrate that the cuspiness of dark matter halos is correlated with the stellar surface brightness. We compare the rotation curve fits to the SPARC sample from a self-interacting dark matter (SIDM) model, which self-consistently includes the impact of baryons on the halo profile, and hydrodynamical N-body simulations with cold dark matter (CDM). The SIDM model predicts a strong correlation between the core size and the stellar surface density, and it provides the best global fit to the data. The CDM simulations without strong baryonic feedback effects fail to explain the large dark matter cores seen in low surface brightness galaxies. On the other hand, with strong feedback, CDM simulations do not produce galaxy analogs with high stellar and dark matter densities, and therefore they have trouble in explaining the rotation curves of high surface brightness galaxies. This implies that current feedback implementations need to be modified. We also explicitly show how the concentration-mass and stellar-to-halo mass relations together lead to a radial acceleration relation (RAR) in an averaged sense, and reiterate the point that the RAR does not capture the diversity of galaxy rotation curves in the inner regions. These results make a strong case for SIDM as the explanation for the cores and cusps of field galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

Neutral hydrogen in cometary comas.

The strong L alpha radiation observed recently in comets Tago-Sato-Kosaka and Bennett can be explained in terms of the resonant scattering of solar L alpha radiation on neutral hydrogen formed by the photodissociation of H2O which is vaporized from a nucleus having an ice core. A complete hydrodynamic description of an atmosphere composed of H2O and its daughter products OH, H and O coupled through frictional interaction as well as production and loss processes is given. Numerical results are computed in a typical case, and it is found that a temperature of about 3000 K for the cometary atmosphere provides the best fit with observation.

Mendis, D. A.↗

Initiating solar system formation through stellar shock waves

Isotopic anomalies in presolar grains and other meteoritical components require nucleosynthesis in stellar interiors, condensation into dust grains in stellar envelopes, transport of the grains through the interstellar medium by stellar outflows, and finally injection of the grains into the presolar nebula. The proximity of the presolar cloud to these energetic stellar events suggests that a shock wave from a stellar outflow might have initiated the collapse of an otherwise stable presolar cloud. We have begun to study the interactions of stellar shock waves with thermally supported, dense molecular cloud cores, using a three spatial dimension (3D) radiative hydrodynamics code. Supernova shock waves have been shown by others to destroy quiescent clouds, so we are trying to determine if the much smaller shock speeds found in, e.g., asymptotic giant branch (AGB) star winds, are strong enough to initiate collapse in an otherwise stable, rotating, solar-mass cloud core, without leading to destruction of the cloud.

Boss, A. P.↗

Experimental study of unsteady hydrodynamic force matrices on whirling centrifugal pump impellers

An experimental facility was constructed and instrumented. A set of centrifugal flow pumps whose impellers were made to follow a controlled circular whirl motion were studied. The aim was to characterize the steady and unsteady fluid forces measured on the impeller under various pump operating conditions. The postulation was that the unsteady lateral forces result from interactions between the impeller and the surrounding diffuser and/or volute (via the working fluid), and that under certain flow regimes these forces can drive unstable lateral motions of the pump rotor. The lateral hydrodynamic forces were decomposed into their steady and unsteady parts, the latter being further expressed in terms of a generalized fluid stiffness matrix. Conclusions regarding the effect of impeller geometry could not be reached given the similarity of the tested designs. However, other results on phenomena such as skin friction and leakage flow are presented. Some of the findings are compared to experimental and theoretical data from other sources. Finally, the rotordynamic consequences of the results are discussed as the present data were applied to the case of the Space Shuttle Main Engine (SSME) High Pressure Oxidizer Turbopump (HPOTP).

Belgacem, Jery↗

Collinear dual-pulse laser optical breakdown and energy deposition

A non-equilibrium model for laser-generated plasmas is used to represent collinear dual nano-second-pulse interactions. In the specific case considered, the breakdown is initiated with an ultraviolet (UV) laser pulse at 266 nm, which pre-ionizes the gas, and it is followed by a near-infrared (NIR) pulse at 1064 nm, which deposits significant energy into the ionized mixture. The model is validated against corresponding experiments, and simulation results are interrogated to understand key features of the plasma-kernel dynamics and the post-discharge hydrodynamics. The hydrodynamics of the non-equilibrium plasma is governed by the two-temperature Navier–Stokes equations accounting for both multiphoton ionization and inverse bremsstrahlung. The interaction between the laser beam and the plasma is modeled based on the radiative transfer equation. The temporal and spatial offsets of the two pulses can generate ionization kernels with different topology and dynamics. It is shown that the UV pre-ionization pulse can tailor the plasma region, leading to a larger ionized volume than would occur for a single-pulse breakdown, and it increases the efficiency of the energy deposition for the following NIR discharge. Vorticity in the early post-discharge phase is generated via baroclinic torque from the misalignment of the radial gradient of density (sudden gas expansion) and the strong pressure gradient (initiated by energy deposition from the laser beam). It is predicted that the collinear dual pulse can be designed to specify the plasma kernel size, shape, and maximum temperature in the breakdown phase, and the initial post-breakdown vorticity and plasma-core decay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the X-ray spectrum of Kepler's supernova remnant

We have devised a method to do nonequilibrium ionization calculations on the results of two-dimensional hydrodynamical simulations, based on the algorithm of Hughes & Helfand (1985). We have calculated the ionization structure and X-ray emission for a two-dimensional numerical hydrodynamical simulation for the remnant of Kepler's supernova (SN); the hydrodynamical model was presented in a previous paper. In this model, the progenitor of Kepler's SN is assumed to have been a massive runaway star ejected from the Galactic plane. In its red supergiant stage, its dense stellar wind was distorted and compressed into a bow shock by the ram pressure of the tenuous interstellar medium. The subsequent interaction of the supernova ejecta with this asymmetric circumstellar matter produced a strongly asymmetric supernova remnant (SNR). In this paper, we present calculated X-ray spectra for this hydrodynamical model. A comparison with observations implies a moderate overabundance of Fe in Kepler's SNR (only 50% larger than its cosmic value), in contrast to a large (6 to 15) Fe overabundance derived previously. However, we confirm earlier conclusions that Si and S abundances are 2 to 3 times solar. These modest enhancements of Si, S, and Fe may be attributed either to heavy-element enriched SN ejecta or to the initial chemical abundances of the SN progenitor, which originated in the metal-rich inner Galaxy. The comparison of our models with the observed spectra confirm theoretical predictions that moderate electron heating occurs at strong collisionless shock fronts, with the implied electron/mean temperature ratio of approximately 0.5.

Borkowski, Kazimierz J.↗

The role of cosmic rays in magnetic hydrodynamics of interstellar medium

Cosmic ray (CR) propagation in the Galaxy and generally in the cosmic plasma is usually considered in the diffusion approximation. The diffusion is regarded to result from CR scattering due to their interaction with a magnetic and an electric field. In most cases the fields are assumed to be given. Meanwhile, in the Galaxy the CR energy density w sub cr is similar to I eV/cm, i.e., it is comparable with the energy densities of the magnetic field and turbulent motions in the interstellar gas. Therefore, for the Galaxy it becomes necessary to take into account the influence of CR on the gas dynamics and on the magnetic fields in this gas. The simplest way to this is to use the hydrodynamic approximation, but this is possible only on scales greatly exceeding the CR free path lambda before scattering and only for times larger than lambda/v approx. equals lambda/c. One should thus obtain corresponding MHD equations and establish the limits of their applicability.

Ginzburg, V. L.↗

Disentangling Planets from Photoelectric Instability in Gas-rich Optically Thin Dusty Disks

Structures in circumstellar disks such as gaps and rings are often attributed to planets. This connection has been difficult to show unequivocally, as other processes may also produce these features. In particular, a photoelectric instability (PEI) has been proposed, operating in gas-rich optically thin disks, that generates structures predicted by planet–disk interactions. We examine the question of how to disentangle the planetary effects on disk structure from the effects of the PEI. We use the Pencil Code to perform 2D global hydrodynamical models of the dynamics of gas and dust in a thin disk with and without planetary perturbers. Photoelectric heating is modeled with an equation of state where pressure is proportional to dust surface density. The drag force on grains and its backreaction on the gas are included. Analyzing the situation without PEI, we find that gas–dust interactions alter the shape of the planetary gap from the dust-free case when the local dust-to-gas ratio ε approaches unity. This result also applies to primordial disks, because dust drifting inward accumulates at the edge of the planetary gap, and any initial dust-to-gas ratio eventually achieves ε = 1 if the dust reservoir is sufficient. We find a result particular to high dust-to-gas ratio disks as well: as dust drifts inward, the dust front becomes a sharp transition, and the backreaction triggers the Rossby wave instability. When PEI is included, we find that it obscures structures induced by planets unless the planet's mass is sufficiently large to carve a noticeable gap. Specifically, the instability generates arcs and rings of regular spacing: a planet is discernible when it carves a dust gap wider than the wavelength of the PEI.

Areli Castrejon↗

On the relation between 'mixing length' and 'direct interaction approximation' theories of turbulence

The capabilities of direct interaction approximations (DIA) and weak-coupling approximations (WCA) in modeling turbulence in plasma physics studies are examined. DIA and WCA formalisms are reviewed and a DIA expression for homogeneous incompressible hydrodynamics is simplified to derive an energy spectrum in an arbitrary wavenumber space. A solution is defined for the expression when applied to long wavelength phenomena generated by an instability. The equation is further generalized to yield a transport equation for modeling the spectrum generated by homogeneous turbulence in a linearly unstable system undergoing nonlinear interaction. Finally, the technique is applied to two-dimensional turbulence in a low pressure, weakly ionized plasma embedded in a homogeneous magnetic field and experiencing instabilities caused by electron density and electric potential gradients, such as found in the ionospheric E region.

Sudan, R. N.↗

Quantum simulations of hydrodynamics via the Madelung transformation

Developing numerical methods to simulate efficiently nonlinear fluid dynamics on universal quantum computers is a challenging problem. In this paper, a generalization of the Madelung transform is defined to solve quantum relativistic charged fluid equations interacting with external electromagnetic forces via the Dirac equation. The Dirac equation is discretized into discrete-time quantum walks which can be efficiently implemented on universal quantum computers. A variant of this algorithm is proposed to implement simulations using current noisy intermediate scale quantum (NISQ) devices in the case of homogeneous external forces. High resolution (up to N=2 17 grid points) numerical simulations of relativistic and nonrelativistic hydrodynamical shocks on current IBM NISQs are performed with this algorithm. Here, this paper demonstrates that fluid dynamics can be simulated on NISQs, and opens the door to simulating other fluids, including plasmas, with more general quantum walks and quantum automata.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Diffusivity and Structure of Room Temperature Ionic Liquid in Various Organic Solvents

Room-temperature ionic liquids (RTILs) hold promise for applications in electric double layer capacitors (EDLCs), owing to a much wider potential window, lower vapor pressure, and better thermal and chemical stabilities compared to conventional aqueous and organic electrolytes. However, because the low diffusivity of ions in neat RTILs negates the EDLCs’ advantage of high power density, the ionic liquids are often used in mixture with organic solvents. In this study, we measured the diffusivity of cations and anions in RTIL, 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl) ([BMIM + ][TFSI – ]), mixed with 10 organic solvents, by using the pulsed-field gradient NMR method. The ion diffusivity was found to follow that of neat solvents and in most studied solvents showed an excellent agreement with the predicted values reported in the recent molecular dynamics (MD) study [Thompson, M. W.; J. Phys. Chem. B 2019, 123, 1340-1347]. In two solvents consisting of long-chain molecules, however, the MD simulations predictions slightly underestimated the ionic diffusivities. The degree of ion dissociation was also estimated for each solvent by comparing the ionic conductivity with the molar conductivity derived from the diffusion measurements. The degree of ion dissociation and the hydrodynamic radius of ions suggest that the ions are coordinated by ~1 solvent molecule. Finally, the scarcity of solvent–ion interactions explains the fact that the diffusivity of ions in the mixture significantly depends on the viscosity of the solvent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The turbulence driving parameter of molecular clouds in disc galaxies

ABSTRACT Supersonic turbulence plays a pivotal role during the formation of molecular clouds and stars in galaxies. However, little is known about how the fraction of compressive and solenoidal modes in the velocity field evolves over time and how it depends on properties of the molecular cloud or the galactic environment. In this work, we carry out magnetohydrodynamical simulations of disc galaxies and study the time evolution of the turbulence driving parameter for an ensemble of clouds. We find that the time-averaged turbulence driving parameter is insensitive to the position of the cloud within the galaxy. The ensemble-averaged driving parameter is found to be rather compressive with b ∼ 0.5–0.7, indicating almost time-independent global star formation properties. However, each individual cloud shows a highly fluctuating driving parameter, which would strongly affect the cloud’s star formation rate. We find that the mode of turbulence driving can rapidly change within only a few Myr, both from solenoidal to compressive and vice versa. We attribute these changes to cloud collisions and to tidal interactions with clouds or overdensities in the environment. Last, we find no significant differences in the average driving parameter between hydrodynamic and initially strongly magnetized galaxies. However, the magnetic field tends to reduce the overall fluctuation of the driving parameter. The average driving and its uncertainty are seen to be in agreement with recent constraints on the turbulence driving mode for solar neighbourhood clouds.

Körtgen, Bastian↗

Techno-Economic Modelling of Tidal Energy Converter Arrays in the Tacoma Narrows

Hydrokinetic tidal energy converter (TEC) technology is yet to become cost competitive with other renewable energy sources. Understanding the interaction between energy production and the costs incurred harvesting that energy may unlock the economic potential of this technology. Although hydrodynamic simulation of TEC arrays has matured over time, including demonstration of how small and large arrays affect the resource, integration of cost modelling is often limited. The advanced ocean energy array techno-economic modelling tool ‘DTOcean’ enables designers to calculate and improve the levelised cost of energy (LCOE) of an array through parametric simulation of the energy extraction, design of the electrical network, moorings and foundations, and simulation of the installation and lifetime operations and maintenance of the array. This work presents a verification of DTOcean’s ability to simulate the techno-economic performance of TEC arrays by reproducing the hypothetical RM1 reference model, a semi-analytical model of a TEC array based in the Tacoma Narrows of Washington state, U.S.A. It is demonstrated that DTOcean can produce a reasonable estimate to the LCOE predicted by the reference model, giving (in Euro cents per kiloWatt hour) 36.69 ¢/kWh against the reference model’s 34.612 ¢/kWh for 10 TECs, while for 50 TECs, DTOcean calculated 20.34 ¢/kWh compared to 17.34 ¢/kWh for the reference model.

Topper, Mathew B. R. (ORCID:0000000347324347)↗

Hot Circumsingle Disks Drive Binary Black Hole Mergers in Active Galactic Nucleus Disks

Abstract Binary black hole (BBH) mergers, particularly those with component masses in the pair-instability gap, may be produced by hierarchical mergers in the disks surrounding Active Galactic Nuclei (AGNs). While the interaction of an embedded BBH with an AGN disk is typically assumed to facilitate a merger, recent high-resolution hydrodynamical simulations challenge this assumption. However, these simulations often have simplified treatments for gas thermodynamics. In this work, we model the possible consequence of various feedback from an embedded BBH with a simple model that maintains an enhanced temperature profile around each binary component. We show that when the minidisks around each BH become hotter than the background by a factor of three, the BBH orbital evolution switches from expansion to contraction. By analyzing the gravitational torque profile, we find that this change in direction is driven by a weakening of the minidisk spirals and their positive torque on the binary. Our results highlight the important role of thermodynamics around BBHs and its effect on their orbital evolution, suggesting that AGN disks could be efficient factories for BBH mergers.

79 ASTRONOMY AND ASTROPHYSICS↗

Migration of Accreting Planets in Radiative Discs From Dynamical Torques

We present the results of hydrodynamical simulations of the orbital evolution of planets undergoing runaway gas accretion in radiative discs. We consider accreting disc models with constant mass flux through the disc, and where radiative cooling balances the effect of viscous heating and stellar irradiation. We assume that 20–30 M(sub ⨁) giant planet cores are formed in the region where viscous heating dominates and migrate outward under the action of a strong entropy-related corotation torque. In the case where gas accretion is neglected and for an α viscous stress parameter α = 2 × 10(exp -3), we find evidence for strong dynamical torques in accreting discs with accretion rates M ̇ ≳ 7 × 10(exsp -8)M(sub ⨀) yr(exp -1). Their main effect is to increase outward migration rates by a factor of ∼2 typically. In the presence of gas accretion, however, runaway outward migration is observed with the planet passing through the zero-torque radius and the transition between the viscous heating and stellar heating dominated regimes. The ability for an accreting planet to enter a fast migration regime is found to depend strongly on the planet growth rate but can occur for values of the mass flux through the disc of M ̇ ≳ 5 × 10(exp -8)M(sub ⨀) yr(exp -1).We find that an episode of runaway outward migration can cause an accreting planet formed in the 5–10 au region to temporarily orbit at star–planet separations as large as∼60–70 au. However, increase in the amplitude of the Lindblad torque associated with planet growth plus change in the streamline topology near the planet systematically cause the direction of migration to be reversed. Subsequent evolution corresponds to the planet migrating inward rapidly until it becomes massive enough to open a gap in the disc and migrate in the type II regime. Our results indicate that a planet can reach large orbital distances under the combined effect of dynamical torques and gas accretion, but an alternative mechanism is required to explain the presence of massive planets on wide orbits.

Accretion↗