Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “VISCOSITY”

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 1,675 records · Page 93

SIMULATIONS OF ACCRETION POWERED SUPERNOVAE IN THE PROGENITORS OF GAMMA-RAY BURSTS

Observational evidence suggests a link between long-duration gamma-ray bursts (LGRBs) and Type Ic supernovae. Here, we propose a potential mechanism for Type Ic supernovae in LGRB progenitors powered solely by accretion energy. We present spherically symmetric hydrodynamic simulations of the long-term accretion of a rotating gamma-ray burst progenitor star, a 'collapsar', onto the central compact object, which we take to be a black hole. The simulations were carried out with the adaptive mesh refinement code FLASH in one spatial dimension and with rotation, an explicit shear viscosity, and convection in the mixing length theory approximation. Once the accretion flow becomes rotationally supported outside of the black hole, an accretion shock forms and traverses the stellar envelope. Energy is carried from the central geometrically thick accretion disk to the stellar envelope by convection. Energy losses through neutrino emission and nuclear photodisintegration are calculated but do not seem important following the rapid early drop of the accretion rate following circularization. We find that the shock velocity, energy, and unbound mass are sensitive to convective efficiency, effective viscosity, and initial stellar angular momentum. Our simulations show that given the appropriate combinations of stellar and physical parameters, explosions with energies {approx}5 Multiplication-Sign 10{sup 50} erg, velocities {approx}3000 km s{sup -1}, and unbound material masses {approx}> 6 M{sub Sun} are possible in a rapidly rotating 16 M{sub Sun} main-sequence progenitor star. Further work is needed to constrain the values of these parameters, to identify the likely outcomes in more plausible and massive LRGB progenitors, and to explore nucleosynthetic implications.

79 ASTRONOMY AND ASTROPHYSICS↗

ON HYDROMAGNETIC STRESSES IN ACCRETION DISK BOUNDARY LAYERS

Detailed calculations of the physical structure of accretion disk boundary layers, and thus their inferred observational properties, rely on the assumption that angular momentum transport is opposite to the radial angular frequency gradient of the disk. The standard model for turbulent shear viscosity satisfies this assumption by construction. However, this behavior is not supported by numerical simulations of turbulent magnetohydrodynamic (MHD) accretion disks, which show that angular momentum transport driven by the magnetorotational instability (MRI) is inefficient in disk regions where, as expected in boundary layers, the angular frequency increases with radius. In order to shed light on physically viable mechanisms for angular momentum transport in this inner disk region, we examine the generation of hydromagnetic stresses and energy density in differentially rotating backgrounds with angular frequencies that increase outward in the shearing-sheet framework. We isolate the modes that are unrelated to the standard MRI and provide analytic solutions for the long-term evolution of the resulting shearing MHD waves. We show that, although the energy density of these waves can be amplified significantly, their associated stresses oscillate around zero, rendering them an inefficient mechanism to transport significant angular momentum (inward). These findings are consistent with the results obtained in numerical simulations of MHD accretion disk boundary layers and challenge the standard assumption of efficient angular momentum transport in the inner disk regions. This suggests that the detailed structure of turbulent MHD accretion disk boundary layers could differ appreciably from those derived within the standard framework of turbulent shear viscosity.

79 ASTRONOMY AND ASTROPHYSICS↗

OUTWARD MIGRATION OF JUPITER AND SATURN IN EVOLVED GASEOUS DISKS

The outward migration of a pair of resonant-orbit planets, driven by tidal interactions with a gas-dominated disk, is studied in the context of evolved solar nebula models. The planets' masses, M{sub 1} and M{sub 2}, correspond to those of Jupiter and Saturn. Hydrodynamical calculations in two and three dimensions are used to quantify the migration rates and analyze the conditions under which the outward migration mechanism may operate. The planets are taken to be fully formed after 10{sup 6} and before 3 Multiplication-Sign 10{sup 6} years. The orbital evolution of the planets in an evolving disk is then calculated until the disk's gas is completely dissipated. Orbital locking in the 3:2 mean motion resonance may lead to outward migration under appropriate conditions of disk viscosity and temperature. However, resonance locking does not necessarily result in outward migration. This is the case, for example, if convergent migration leads to locking in the 2:1 mean motion resonance, as post-formation disk conditions seem to suggest. Accretion of gas on the planets may deactivate the outward migration mechanism by raising the mass ratio M{sub 2}/M{sub 1} and/or by reducing the accretion rate toward the star, and hence depleting the inner disk. For migrating planets locked in the 3:2 mean motion resonance, there are stalling radii that depend on disk viscosity and on stellar irradiation, when it determines the disk's thermal balance. Planets locked in the 3:2 orbital resonance that start moving outward from within 1-2 AU may reach beyond Almost-Equal-To 5 AU only under favorable conditions. However, within the explored space of disk parameters, only a small fraction-less than a few percent-of the models predict that the interior planet reaches beyond Almost-Equal-To 4 AU.

79 ASTRONOMY AND ASTROPHYSICS↗

DISK CORONA INTERACTION: MECHANISM FOR THE DISK TRUNCATION AND SPECTRUM CHANGE IN LOW-LUMINOSITY ACTIVE GALACTIC NUCLEI

The truncation of an optically thick, geometrically thin accretion disk is investigated in the context of low-luminosity AGNs (LLAGNs). We generalize the disk evaporation model used in the interpretative framework of black hole X-ray binaries by including the effect of a magnetic field in accretion disks surrounding supermassive black holes. The critical transition mass accretion rate for which the disk is truncated is found to be insensitive to magnetic effects, but its inclusion leads to a smaller truncation radius in comparison to a model without its consideration. That is, a thin viscous disk is truncated for LLAGNs at an Eddington ratio less than 0.03 for a standard viscosity parameter ({alpha} = 0.3). An increase of the viscosity parameter results in a higher critical transition mass accretion rate and a correspondingly smaller truncation distance, the latter accentuated by greater magnetic energy densities in the disk. Based on these results, the truncation radii inferred from spectral fits of LLAGNs published in the literature are consistent with the disk evaporation model. The infrared emission arising from the truncated geometrically thin accretion disks may be responsible for the red bump seen in such LLAGNs.

79 ASTRONOMY AND ASTROPHYSICS↗

VERTICAL STRUCTURE OF NEUTRINO-DOMINATED ACCRETION DISKS AND NEUTRINO TRANSPORT IN THE DISKS

We investigate the vertical structure of neutrino-dominated accretion disks by self-consistently considering the detailed microphysics, such as the neutrino transport, the vertical hydrostatic equilibrium, the conservation of lepton number, as well as the balance between neutrino cooling, advection cooling, and viscosity heating. After obtaining the emitting spectra of neutrinos and antineutrinos by solving the one-dimensional Boltzmann equation of neutrino and antineutrino transport in the disk, we calculate the neutrino/antineutrino luminosity and their annihilation luminosity. We find that the total neutrino and antineutrino luminosity is about 10{sup 54} erg s{sup -1} and their annihilation luminosity is about 5 Multiplication-Sign 10{sup 51} erg s{sup -1} with an extreme accretion rate 10 M{sub Sun} s{sup -1} and an alpha viscosity {alpha} = 0.1. In addition, we find that the annihilation luminosity is sensitive to the accretion rate and will not exceed 10{sup 50} erg s{sup -1}, which is not sufficient to power the fireball of most energetic gamma-ray bursts (GRBs) if the accretion rate is lower than 1 M{sub Sun} s{sup -1}. Therefore, the effects of the spin of the black hole or/and the magnetic field in the accretion flow might play a role in powering the central engine of GRBs.

79 ASTRONOMY AND ASTROPHYSICS↗

ON THE UNUSUAL GAS COMPOSITION IN THE {beta} PICTORIS DEBRIS DISK

The metallic gas associated with the {beta} Pic debris disk is not believed to be primordial, but arises from the destruction of dust grains. Recent observations have shown that carbon and oxygen in this gas are exceptionally overabundant compared to other elements, by some 400 times. We study the origin of this enrichment under two opposing hypotheses: preferential production, where the gas is produced with the observed unusual abundance (as may happen if gas is produced by photodesorption from C/O-rich icy grains), and preferential depletion, where the gas evolves to the observed state from an original solar abundance (if outgassing occurs under high-speed collisions) under a number of dynamical processes. We include in our study the following processes: radiative blowout of metallic elements, dynamical coupling between different species, and viscous accretion onto the star. We find that, if gas viscosity is sufficiently low (the conventional {alpha} parameter {approx}< 10{sup -3}), differential blowout dominates. While gas accumulates gradually in the disks, metallic elements subject to strong radiation forces, such as Na and Fe, deplete more quickly than C and O, naturally leading to the observed overabundance of C and O. On the other hand, if gas viscosity is high ({alpha} {approx}> 10{sup -1}, as expected for this largely ionized disk), gas is continuously produced and viscously accreted toward the star. This removal process does not discriminate between elements so the observed overabundance of C and O has to be explained by a preferential production that strongly favors C and O to other metallic elements. One such candidate is photodesorption off the grains. We compare our calculation against all observed elements ({approx}10) in the gas disk and find a mild preference for the second scenario, based on the abundance of Si alone. If true, {beta} Pic should still be accreting at an observable rate, well after its primordial disk has disappeared.

79 ASTRONOMY AND ASTROPHYSICS↗

MAGNETOHYDRODYNAMIC SIMULATION OF A DISK SUBJECTED TO LENSE-THIRRING PRECESSION

When matter orbits around a central mass obliquely with respect to the mass's spin axis, the Lense-Thirring effect causes it to precess at a rate declining sharply with radius. Ever since the work of Bardeen and Petterson, it has been expected that when a fluid fills an orbiting disk, the orbital angular momentum at small radii should then align with the mass's spin. Nearly all previous work has studied this alignment under the assumption that a phenomenological 'viscosity' isotropically degrades fluid shears in accretion disks, even though it is now understood that internal stress in flat disks is due to anisotropic MHD turbulence. In this paper we report a pair of matched simulations, one in MHD and one in pure (non-viscous) HD in order to clarify the specific mechanisms of alignment. As in the previous work, we find that disk warps induce radial flows that mix angular momentum of different orientation; however, we also show that the speeds of these flows are generically transonic and are only very weakly influenced by internal stresses other than pressure. In particular, MHD turbulence does not act in a manner consistent with an isotropic viscosity. When MHD effects are present, the disk aligns, first at small radii and then at large; alignment is only partial in the HD case. We identify the specific angular momentum transport mechanisms causing alignment and show how MHD effects permit them to operate more efficiently. Last, we relate the speed at which an alignment front propagates outward (in the MHD case) to the rate at which Lense-Thirring torques deliver angular momentum at smaller radii.

79 ASTRONOMY AND ASTROPHYSICS↗

N-body simulations of terrestrial planet formation under the influence of a hot Jupiter

We investigate the formation of multiple-planet systems in the presence of a hot Jupiter (HJ) using extended N-body simulations that are performed simultaneously with semianalytic calculations. Our primary aims are to describe the planet formation process starting from planetesimals using high-resolution simulations, and to examine the dependences of the architecture of planetary systems on input parameters (e.g., disk mass, disk viscosity). We observe that protoplanets that arise from oligarchic growth and undergo type I migration stop migrating when they join a chain of resonant planets outside the orbit of an HJ. The formation of a resonant chain is almost independent of our model parameters, and is thus a robust process. At the end of our simulations, several terrestrial planets remain at around 0.1 AU. The formed planets are not equal mass; the largest planet constitutes more than 50% of the total mass in the close-in region, which is also less dependent on parameters. In the previous work of this paper, we have found a new physical mechanism of induced migration of the HJ, which is called a crowding-out. If the HJ opens up a wide gap in the disk (e.g., owing to low disk viscosity), crowding-out becomes less efficient and the HJ remains. We also discuss angular momentum transfer between the planets and disk.

79 ASTRONOMY AND ASTROPHYSICS↗

No evidence for Bardeen-Petterson alignment in GRMHD simulations and semi-analytic models of moderately thin, prograde, tilted accretion disks

In this paper, we introduce the first results that use data extracted directly from numerical simulations as inputs to the analytic twisted disk model of Zhuravlev and Ivanov. In both numerical and analytic approaches, fully relativistic models of tilted and twisted disks having a moderate effective viscosity around a slowly rotating Kerr black hole are considered. Qualitatively, the analytic model demonstrates the same dynamics as the simulations, although with some quantitative offset. Namely, the general relativistic magnetohydrodynamic simulations typically give smaller variations of tilt and twist across the disk. When the black hole and the disk rotate in the same direction, the simulated tilted disk and analytic model show no sign of Bardeen-Petterson alignment, even in the innermost parts of the disk where the characteristic time for relaxation to a quasi-stationary configuration is of the same order as the computation time. In the opposite case, when the direction of the disk's rotation is opposite to that of the black hole, a partial alignment is observed, in agreement with previous theoretical estimates. Thus, both fully numerical and analytic schemes demonstrate that the Bardeen-Petterson effect may not be possible for the case of prograde rotation provided that disk's effective viscosity is sufficiently small. This may have implications in modeling of different astrophysical phenomena such as disk spectra and jet orientation.

79 ASTRONOMY AND ASTROPHYSICS↗

THE GENERAL RELATIVISTIC EQUATIONS OF RADIATION HYDRODYNAMICS IN THE VISCOUS LIMIT

We present an analysis of the general relativistic Boltzmann equation for radiation, appropriate to the case where particles and photons interact through Thomson scattering, and derive the radiation energy-momentum tensor in the diffusion limit with viscous terms included. Contrary to relativistic generalizations of the viscous stress tensor that appear in the literature, we find that the stress tensor should contain a correction to the comoving energy density proportional to the divergence of the four-velocity, as well as a finite bulk viscosity. These modifications are consistent with the framework of radiation hydrodynamics in the limit of large optical depth, and do not depend on thermodynamic arguments such as the assignment of a temperature to the zeroth-order photon distribution. We perform a perturbation analysis on our equations and demonstrate that as long as the wave numbers do not probe scales smaller than the mean free path of the radiation, the viscosity contributes only decaying, i.e., stable, corrections to the dispersion relations. The astrophysical applications of our equations, including jets launched from super-Eddington tidal disruption events and those from collapsars, are discussed and will be considered further in future papers.

79 ASTRONOMY AND ASTROPHYSICS↗

STRIPPED ELLIPTICAL GALAXIES AS PROBES OF ICM PHYSICS. I. TAILS, WAKES, AND FLOW PATTERNS IN AND AROUND STRIPPED ELLIPTICALS

Elliptical cluster galaxies are progressively stripped of their atmospheres due to their motion through the intracluster medium (ICM). Deep X-ray observations reveal the fine-structure of the galaxy’s remnant atmosphere and its gas tail and wake. This fine-structure depends on dynamic conditions (galaxy potential, initial gas contents, orbit through the host cluster), orbital stage (early infall, pre-/post-pericenter passage), and ICM plasma properties (thermal conductivity, viscosity, magnetic field structure). We aim to disentangle dynamic and plasma effects in order to use stripped ellipticals as probes of ICM plasma properties. This first paper of a series investigates the hydrodynamics of progressive gas stripping by means of inviscid hydrodynamical simulations. We distinguish a long-lasting initial relaxation phase and a quasi-steady stripping phase. During quasi-steady stripping, the ICM flow around the remnant atmosphere resembles the flow around solid bodies, including a “deadwater” region in the near wake. Gas is stripped from the remnant atmosphere predominantly at its sides via Kelvin–Helmholtz instabilities. The downstream atmosphere is largely shielded from the ICM wind and thus shaped into a tail. Observationally, both this “remnant tail” and the stripped gas in the wake can appear as a “tail”, but only in the wake can galactic gas mix with the ambient ICM. While the qualitative results are generic, the simulations presented here are tailored to the Virgo elliptical galaxy M89 (NGC 4552) for the most direct comparison to observations. Papers II and III of this series describe the effect of viscosity and compare to Chandra and XMM-Newton observations, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

Halting migration: Numerical calculations of corotation torques in the weakly nonlinear regime

Planets in their formative years can migrate due to the influence of gravitational torques in the protoplanetary disk they inhabit. For low-mass planets in an isothermal disk, it is known that there is a strong negative torque on the planet due to its linear perturbation to the disk, causing fast inward migration. The current investigation demonstrates that in these same isothermal disks, for intermediate-mass planets, there is a strong positive nonlinear corotation torque due to the effects of gas being pulled through a gap on horseshoe orbits. For intermediate-mass planets, this positive torque can partially or completely cancel the linear (Type I) torque, leading to slower or outward migration, even in an isothermal disk. The effect is most significant for super-Earth and sub-Jovian planets, during the transition from a low-mass linear perturber to a nonlinear gap-opening planet, when the planet has opened a so-called “partial gap,” though the precise values of these transition masses depend sensitively on the disk model (density profile, viscosity, and disk aspect ratio). In this study, numerical calculations of planet–disk interactions calculate these torques explicitly, and scalings are empirically constructed for migration rates in this weakly nonlinear regime. These results find outward migration is possible for planets with masses in the range 20–100 M{sub ⊕}, though this range depends on the disk model considered. In the disk models where torque reversal occurs, the critical planet-to-star mass ratio for torque reversal was found to have the robust scaling q{sub crit}∝√α(h/r){sup 3}, where α is the dimensionless viscosity parameter and h/r is the disk aspect ratio.

79 ASTRONOMY AND ASTROPHYSICS↗

ON THE SPATIAL SCALES OF WAVE HEATING IN THE SOLAR CHROMOSPHERE

Dissipation of magnetohydrodynamic (MHD) wave energy has been proposed as a viable heating mechanism in the solar chromospheric plasma. Here, we use a simplified one-dimensional model of the chromosphere to theoretically investigate the physical processes and spatial scales that are required for the efficient dissipation of Alfvén waves and slow magnetoacoustic waves. We consider the governing equations for a partially ionized hydrogen-helium plasma in the single-fluid MHD approximation and include realistic wave damping mechanisms that may operate in the chromosphere, namely, Ohmic and ambipolar magnetic diffusion, viscosity, thermal conduction, and radiative losses. We perform an analytic local study in the limit of small amplitudes to approximately derive the lengthscales for critical damping and efficient dissipation of MHD wave energy. We find that the critical dissipation lengthscale for Alfvén waves depends strongly on the magnetic field strength and ranges from 10 m to 1 km for realistic field strengths. The damping of Alfvén waves is dominated by Ohmic diffusion for weak magnetic field and low heights in the chromosphere, and by ambipolar diffusion for strong magnetic field and medium/large heights in the chromosphere. Conversely, the damping of slow magnetoacoustic waves is less efficient, and spatial scales shorter than 10 m are required for critical damping. Thermal conduction and viscosity govern the damping of slow magnetoacoustic waves and play an equally important role at all heights. These results indicate that the spatial scales at which strong wave heating may work in the chromosphere are currently unresolved by observations.

79 ASTRONOMY AND ASTROPHYSICS↗

An extended magnetohydrodynamics model for relativistic weakly collisional plasmas

Black holes that accrete far below the Eddington limit are believed to accrete through a geometrically thick, optically thin, rotationally supported plasma that we will refer to as a radiatively inefficient accretion flow (RIAF). RIAFs are typically collisionless in the sense that the Coulomb mean free path is large compared to GM/c{sup 2}, and relativistically hot near the event horizon. In this paper we develop a phenomenological model for the plasma in RIAFs, motivated by the application to sources such as Sgr A* and M87. The model is derived using Israel–Stewart theory, which considers deviations up to second order from thermal equilibrium, but modified for a magnetized plasma. This leads to thermal conduction along magnetic field lines and a difference in pressure, parallel and perpendicular to the field lines (which is equivalent to anisotropic viscosity). In the non-relativistic limit, our model reduces to the widely used Braginskii theory of magnetized, weakly collisional plasmas. We compare our model to the existing literature on dissipative relativistic fluids, describe the linear theory of the plasma, and elucidate the physical meaning of the free parameters in the model. We also describe limits of the model when the conduction is saturated and when the viscosity implies a large pressure anisotropy. In future work, the formalism developed in this paper will be used in numerical models of RIAFs to assess the importance of non-ideal processes for the dynamics and radiative properties of slowly accreting black holes.

79 ASTRONOMY AND ASTROPHYSICS↗

AB INITIO SIMULATIONS FOR MATERIAL PROPERTIES ALONG THE JUPITER ADIABAT

We determine basic thermodynamic and transport properties of hydrogen-helium-water mixtures for the extreme conditions along Jupiter's adiabat via ab initio simulations, which are compiled in an accurate and consistent data set. In particular, we calculate the electrical and thermal conductivity, the shear and longitudinal viscosity, and diffusion coefficients of the nuclei. We present results for associated quantities like the magnetic and thermal diffusivity and the kinematic shear viscosity along an adiabat that is taken from a state-of-the-art interior structure model. Furthermore, the heat capacities, the thermal expansion coefficient, the isothermal compressibility, the Grueneisen parameter, and the speed of sound are calculated. We find that the onset of dissociation and ionization of hydrogen at about 0.9 Jupiter radii marks a region where the material properties change drastically. In the deep interior, where the electrons are degenerate, many of the material properties remain relatively constant. Our ab initio data will serve as a robust foundation for applications that require accurate knowledge of the material properties in Jupiter's interior, e.g., models for the dynamo generation.

79 ASTRONOMY AND ASTROPHYSICS↗

Bulk viscous matter-dominated Universes: asymptotic properties

By means of a combined use of the type Ia supernovae and H(z) data tests, together with the study of the asymptotic properties in the equivalent phase space — through the use of the dynamical systems tools — we demonstrate that the bulk viscous matter-dominated scenario is not a good model to explain the accepted cosmological paradigm, at least, under the parametrization of bulk viscosity considered in this paper. The main objection against such scenarios is the absence of conventional radiation and matter-dominated critical points in the phase space of the model. This entails that radiation and matter dominance are not generic solutions of the cosmological equations, so that these stages can be implemented only by means of unique and very specific initial conditions, i. e., of very unstable particular solutions. Such a behavior is in marked contradiction with the accepted cosmological paradigm which requires of an earlier stage dominated by relativistic species, followed by a period of conventional non-relativistic matter domination, during which the cosmic structure we see was formed. Also, we found that the bulk viscosity is positive just until very late times in the cosmic evolution, around z < 1. For earlier epochs it is negative, been in tension with the local second law of thermodynamics.

79 ASTRONOMY AND ASTROPHYSICS↗

The adhesion model as a field theory for cosmological clustering

The adhesion model has been proposed in the past as an improvement of the Zel'dovich approximation, providing a good description of the formation of the cosmic web. We recast the model as a field theory for cosmological large scale structure, adding a stochastic force to account for power generated from very short, highly non-linear scales that is uncorrelated with the initial power spectrum. The dynamics of this Stochastic Adhesion Model (SAM) is reminiscent of the well known Kardar-Parisi-Zhang equation with the difference that the viscosity and the noise spectrum are time dependent. Choosing the viscosity proportional to the growth factor D restricts the form of noise spectrum through a 1-loop renormalization argument. For this choice, the SAM field theory is renormalizable to one loop. We comment on the suitability of this model for describing the non-linear regime of the CDM power spectrum and its utility as a relatively simple approach to cosmological clustering.

79 ASTRONOMY AND ASTROPHYSICS↗

Data-Driven Machine Learning for Wind Plant Flow Modeling

In this study, we introduce a data-driven machine learning framework for improving the accuracy of wind plant flow models by learning turbulence model corrections based on data from higher-fidelity simulations. First, a high-dimensional PDE-constrained optimization problem is solved using gradient-based optimization with adjoints to determine optimal eddy viscosity fields that improve the agreement of a medium-fidelity Reynolds-Averaged Navier Stokes (RANS) model with large eddy simulations (LES). A supervised learning problem is then constructed to find general, predictive representations of the optimal turbulence closure. A machine learning technique using Gaussian process regression is trained to predict the eddy viscosity field based on local RANS flow field information like velocities, pressures, and their gradients. The Gaussian process is trained on LES simulations of a single turbine and implemented in a wind plant simulation with 36 turbines. We show improvement over the baseline RANS model with the machine learning correction, and demonstrate the ability to provide accurate confidence levels for the corrections that enable future uncertainty quantification studies.

17 WIND ENERGY↗