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,657 records · Page 92

Lattice Boltzmann simulation of immiscible fluid displacement in porous media: Homogeneous versus heterogeneous pore network

Injection of anthropogenic carbon dioxide (CO{sub 2}) into geological formations is a promising approach to reduce greenhouse gas emissions into the atmosphere. Predicting the amount of CO{sub 2} that can be captured and its long-term storage stability in subsurface requires a fundamental understanding of multiphase displacement phenomena at the pore scale. In this paper, the lattice Boltzmann method is employed to simulate the immiscible displacement of a wetting fluid by a non-wetting one in two microfluidic flow cells, one with a homogeneous pore network and the other with a randomly heterogeneous pore network. We have identified three different displacement patterns, namely, stable displacement, capillary fingering, and viscous fingering, all of which are strongly dependent upon the capillary number (Ca), viscosity ratio (M), and the media heterogeneity. The non-wetting fluid saturation (S{sub nw}) is found to increase nearly linearly with logCa for each constant M. Increasing M (viscosity ratio of non-wetting fluid to wetting fluid) or decreasing the media heterogeneity can enhance the stability of the displacement process, resulting in an increase in S{sub nw}. In either pore networks, the specific interfacial length is linearly proportional to S{sub nw} during drainage with equal proportionality constant for all cases excluding those revealing considerable viscous fingering. Our numerical results confirm the previous experimental finding that the steady state specific interfacial length exhibits a linear dependence on S{sub nw} for either favorable (M ≥ 1) or unfavorable (M < 1) displacement, and the slope is slightly higher for the unfavorable displacement.

54 ENVIRONMENTAL SCIENCES↗

Periodically sheared 2D Yukawa systems

We present non-equilibrium molecular dynamics simulation studies on the dynamic (complex) shear viscosity of a 2D Yukawa system. We have identified a non-monotonic frequency dependence of the viscosity at high frequencies and shear rates, an energy absorption maximum (local resonance) at the Einstein frequency of the system at medium shear rates, an enhanced collective wave activity, when the excitation is near the plateau frequency of the longitudinal wave dispersion, and the emergence of significant configurational anisotropy at small frequencies and high shear rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of post-treatment on the characteristics of electrospun poly (vinyl alcohol)/chitosan nanofibers

Electrospun nanofibers have many advantages such as high porosity, easy to be fabricated in various size and high ratio of surface area to volume. This paper presents the preparation of electrospun PVA/Chitosan nanofibers and more specifically focuses on the effect of post-treatment on the permeability and morphology of electrospun PVA/chitosan nanofibers. The mixtures of various concentrations of PVA (6,7,8 wt%)and 2 wt%.chitosan solution (with the ratio of 3:1)were used in electrospun with a constant rate of 0.7 ml/hour. The post-treatment was conducted by immersing in a ethanol or glutaraldehyde solution to performed crosslink structure. The electrospun PVA/Chitosan nanofiber was characterized by scanning electron microscopy (SEM) and fourier transform infrared (FTIR) spectroscopy. The results revealed that the viscosity of the mixture solution is directly proportional to its concentration. Increasing the viscosity increased the diameter of fiber but also made the larger beads formation. FTIR measurement exhibited the existence of relevant functional groups of both PVA and chitosan in the composites.The crosslinked structure was found for the electrospun PVA/Chitosan nanofibers treated with glutaraldehyde solution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Catalytic cracking of the top phase fraction of bio-oil into upgraded liquid oil

The energy consumption is increasing, while oil reserves as a primary energy resource are decreasing, so that is the reason seeking alternative energy source is inevitable. Biomass especially oil palm empty fruit bunches (EFB) which is abundant in Indonesia can be processed into bio-oil by pyrolysis process. The potential for direct substitution of bio-oil for petroleum may be limited due to the high viscosity, high oxygen content, low heating value, and corrosiveness. Consequently, upgrading of the bio-oil before use is inevitable to give a wider variety of applications of its liquid product. Furthermore, upgrading process to improve the quality of bio-oil by reduction of oxygenates involves process such as catalytic cracking. The objective of this research is to study the effect of operation temperature on yield and composition of upgraded liquid oil and to determine physical properties. Bio-oil derived from EFB was upgraded through catalytic cracking using series tubular reactor under atmospheric pressure on a silica-alumina catalyst. Results show that increasing temperature from 450 to 600 °C, resulting in decreasing of upgraded liquid oil (ULO) yield, decreasing viscosity and density of ULO, but increasing in calorimetric value of ULO. The increasing temperature of cracking also will increase the concentration of gasoline and kerosene in ULO.

09 BIOMASS FUELS↗

Structure of cyano-anion ionic liquids: X-ray scattering and simulations

Ionic liquids with cyano anions have long been used because of their unique combination of low-melting temperatures, reduced viscosities, and increased conductivities. Recently we showed that cyano anions in ionic liquids are particularly interesting for their potential use as electron donors to excited state photo-acceptors [B. Wu et al., J. Phys. Chem. B 119, 14790–14799 (2015)]. Here we report on bulk structural and quantum mechanical results for a series of ionic liquids based on the 1-ethyl-3-methylimidazolium cation, paired with the following five cyano anions: SeCN - , SCN - , N(CN) $-\atop{2}$ , C(CN) $-\atop{3}$ , and B(CN) $-\atop{4}$ . By combining molecular dynamics simulations, high-energy X-ray scattering measurements, and periodic boundary condition DFT calculations, we are able to obtain a comprehensive description of the liquid landscape as well as the nature of the HOMO-LUMO states for these ionic liquids in the condensed phase. Features in the structure functions for these ionic liquids are somewhat different than the commonly observed adjacency, charge-charge, and polarity peaks, especially for the bulkiest B(CN) $-\atop{4}$ anion. While the other four cyano-anion ionic liquids present an anionic HOMO, the one for Im $+\atop{2,1}$ /B(CN) $-\atop{4}$ is cationic.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Detailed high-resolution three-dimensional simulations of OMEGA separated reactants inertial confinement fusion experiments

Here, we present results from the comparison of high-resolution three-dimensional (3D) simulations with data from the implosions of inertial confinement fusion capsules with separated reactants performed on the OMEGA laser facility. Each capsule, referred to as a “CD Mixcap,” is filled with tritium and has a polystyrene (CH) shell with a deuterated polystyrene (CD) layer whose burial depth is varied. In these implosions, fusion reactions between deuterium and tritium ions can occur only in the presence of atomic mix between the gas fill and shell material. The simulations feature accurate models for all known experimental asymmetries and do not employ any adjustable parameters to improve agreement with experimental data. Simulations are performed with the RAGE radiation-hydrodynamics code using an Implicit Large Eddy Simulation (ILES) strategy for the hydrodynamics. We obtain good agreement with the experimental data, including the DT/TT neutron yield ratios used to diagnose mix, for all burial depths of the deuterated shell layer. Additionally, simulations demonstrate good agreement with converged simulations employing explicit models for plasma diffusion and viscosity, suggesting that the implicit sub-grid model used in ILES is sufficient to model these processes in these experiments. In our simulations, mixing is driven by short-wavelength asymmetries and longer-wavelength features are responsible for developing flows that transport mixed material towards the center of the hot spot. Mix material transported by this process is responsible for most of the mix (DT) yield even for the capsule with a CD layer adjacent to the tritium fuel. Consistent with our previous results, mix does not play a significant role in TT neutron yield degradation; instead, this is dominated by the displacement of fuel from the center of the implosion due to the development of turbulent instabilities seeded by long-wavelength asymmetries. Through these processes, the long-wavelength asymmetries degrade TT yield more than the DT yield and thus bring DT/TT neutron yield ratios into agreement with experiment. Finally, we present a detailed comparison of the flows in 2D and 3D simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evaluation of toroidal torque by non-resonant magnetic perturbations in tokamaks for resonant transport regimes using a Hamiltonian approach

Toroidal torque generated by neoclassical viscosity caused by external non-resonant, non-axisymmetric perturbations has a significant influence on toroidal plasma rotation in tokamaks. In this article, a derivation for the expressions of toroidal torque and radial transport in resonant regimes is provided within quasilinear theory in canonical action-angle variables. The proposed approach treats all low-collisional quasilinear resonant neoclassical toroidal viscosity regimes including superbanana-plateau and drift-orbit resonances in a unified way and allows for magnetic drift in all regimes. It is valid for perturbations on toroidally symmetric flux surfaces of the unperturbed equilibrium without specific assumptions on geometry or aspect ratio. The resulting expressions are shown to match the existing analytical results in the large aspect ratio limit. Numerical results from the newly developed code NEO-RT are compared to calculations by the quasilinear version of the code NEO-2 at low collisionalities. The importance of the magnetic shear term in the magnetic drift frequency and a significant effect of the magnetic drift on drift-orbit resonances are demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Star formation through thermal instability of radiative plasma with finite electron inertia and finite Larmor radius corrections

I have studied the effects of finite electron inertia, finite ion Larmor radius (FLR) corrections, and radiative heat-loss function on the thermal instability of an infinite homogeneous, viscous plasma incorporating the effect of thermal conductivity for star formation in interstellar medium (ISM). A general dispersion relation is derived using the normal mode analysis method with the help of relevant linearized perturbation equations of the problem. The wave propagation is discussed for longitudinal and transverse directions to the external magnetic field and the conditions of modified thermal instabilities and stabilities are discussed in different cases. We find that the thermal instability criterion is get modified into radiative instability criterion by inclusion of radiative heat-loss functions with thermal conductivity. The viscosity of medium removes the effect of FLR corrections from the condition of radiative instability. Numerical calculation shows stabilizing effect of heat-loss function, viscosity and FLR corrections, and destabilizing effect of finite electron inertia on the thermal instability. Results carried out in this paper shows that stars are formed in interstellar medium mainly due to thermal instability.

79 ASTRONOMY AND ASTROPHYSICS↗

Collisional transport coefficients of dense high-temperature plasmas within the quantum Landau-Fokker-Planck framework

In this paper, we extend the long-established formulas for the transport coefficients of classical plasmas inside the dense plasma regime for temperatures and densities where the classical Landau equation breaks down but its quantum extension that includes quantum degeneracy effects is valid. To this end, the quantum Landau kinetic equation is solved by the Chapman-Enskog method. The mathematical derivation is done in full generality, i.e., for multicomponent systems and to all orders of the polynomials expansion used to approximate the distribution functions. We apply the general results to two important examples, the electron gas model and an electron-ion plasma model consisting of one type of ions of any charge. We discuss the combined effects of the Pauli exclusion principle, of the electron-electron, and of the electron-ion collisions on the transport coefficients and on the convergence of the Chapman-Enskog method. For the electron gas model, the effect of the Pauli exclusion principle on the transport coefficients rapidly becomes non-negligible outside the domain of validity of the classical Landau equation. For the electron-ion plasmas, the effect of the Pauli exclusion principle depends sensitively on the ion charge Z and varies non-monotonically with Θ. For instance, for ion charge Z = 1, the electrical conductivity is increased by up to ~30% compared to its classical value over the range of degeneracy parameters studied, the thermal conductivity is reduced by up to ~9%, and the shear viscosity coefficient is increased by up to ~13%. Finally, in the Lorentz gas ( Z → ∞ ) limit, the electrical conductivity is reduced by up to ~14% compared to its classical value over the range of degeneracy parameters studied, the thermal conductivity is reduced by up to ~39%, and the shear viscosity coefficient is not affected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Finite boundary effects on the spherical Rayleigh–Taylor instability between viscous fluids

For the Rayleigh–Taylor unstable arrangement of a viscous fluid sphere embedded in a finite viscous fluid spherical shell with a rigid boundary and a radially directed acceleration, a dispersion relation is developed from a linear stability analysis using the method of normal modes. [Formula: see text] is the radially directed acceleration at the interface. ρ i denotes the density, μ i is the viscosity, and R i is the radius, where i = 1 is the inner sphere and i = 2 is the outer sphere. The dispersion relation is a function of the following dimensionless variables: viscosity ratio [Formula: see text], density ratio [Formula: see text], spherical harmonic mode n, [Formula: see text], [Formula: see text], and the dimensionless growth rate [Formula: see text], where σ is the exponential growth rate. We show that the boundedness provided by the outer spherical shell has a strong influence on the instability behavior, which is reflected not only in the modulation of the growth rate but also in the selection of the most unstable modes that are physically possible. This outer boundary effect is quantified by the relative magnitude of the radius ratio H. We find that when H is close to unity, lower order harmonics are excluded from becoming the most unstable within a vast region of the parameter space. In other words, the effect of H has precedence over the other controlling parameters d, B, and a wide range of s in establishing what the lowest most unstable mode can be. When H ~ 1, low order harmonics can become the most unstable only for s >> 1. However, in the limit when s → ∞, we show that the most unstable mode is n = 1 and derive the dispersion relation in this limit. The exclusion of most unstable low order harmonics caused by a finite outer boundary is not realized when the outer boundary extends beyond a certain threshold length-scale in which case all modes are equally possible depending on the value of B.

74 ATOMIC AND MOLECULAR PHYSICS↗

Improved analytic modeling of the linear Rayleigh–Taylor instability with plasma transport

The Rayleigh–Taylor (R–T) instability with transport effects has long been studied analytically in the neutral fluid regime. How this picture may differ in plasmas has been given less consideration, however. Leveraging previous numerical simulations with plasma viscosity and diffusion of R–T at a binary plasma interface, we show how plasma transport and kinetic effects alter the linear R–T dispersion relation. Our insight informs a much improved analytical dispersion relation, which better matches the simulation data than previously established models. Additionally, we demonstrate the importance of employing asymptotically correct viscosity coefficients in fluid codes for accurately capturing plasma R–T growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗