An official website of the United States government
ABSTRACT Recent high angular resolution ALMA observations have revealed rich information about protoplanetary discs, including ubiquitous substructures and three-dimensional gas kinematics at different emission layers. One interpretation of these observations is embedded planets. Previous 3D planet–disc interaction studies are either based on viscous simulations or non-ideal magnetohydrodynamics (MHD) simulations with simple prescribed magnetic diffusivities. This study investigates the dynamics of gap formation in 3D non-ideal MHD discs using non-ideal MHD coefficients from the look-up table that is self-consistently calculated based on the thermochemical code. We find a concentration of the poloidal magnetic flux in the planet-opened gap (in agreement with previous work) and enhanced field-matter coupling due to gas depletion, which together enable efficient magnetic braking of the gap material, driving a fast accretion layer significantly displaced from the disc mid-plane. The fast accretion helps deplete the gap further and is expected to negatively impact the planet growth. It also affects the corotation torque by shrinking the region of horseshoe orbits on the trailing side of the planet. Together with the magnetically driven disc wind, the fast accretion layer generates a large, persistent meridional vortex in the gap, which breaks the mirror symmetry of gas kinematics between the top and bottom disc surfaces. Finally, by studying the kinematics at the emission surfaces, we discuss the implications of planets in realistic non-ideal MHD discs on kinematics observations.
more »
« less
Foundations of magnetohydrodynamics
In this tutorial, a derivation of magnetohydrodynamics (MHD) valid beyond the usual ideal gas approximation is presented. Non-equilibrium thermodynamics is used to obtain conservation equations and linear constitutive relations. When coupled with Maxwell's equations, this provides closed fluid equations in terms of material properties of the plasma, described by the equation of state and transport coefficients. These properties are connected to microscopic dynamics using the Irving–Kirkwood procedure and Green–Kubo relations. Symmetry arguments and the Onsager–Casimir relations allow one to vastly simplify the number of independent coefficients. Importantly, expressions for current density, heat flux, and stress (conventionally Ohm's law, Fourier's law, and Newton's law) take different forms in systems with a non-ideal equation of state. The traditional form of the MHD equations, which is usually obtained from a Chapman–Enskog solution of the Boltzmann equation, corresponds to the ideal gas limit of the general equations.
more »
« less
- Award ID(s):
- 2205506
- PAR ID:
- 10612244
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- Physics of Plasmas
- Volume:
- 32
- Issue:
- 7
- ISSN:
- 1070-664X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Abstract We develop a two-dimensional model for the transient diffusion of gas from the cavities in ridge-type structured surfaces to a quiescent liquid suspended above them in the Cassie state to predict the location of the liquid vapor-interface (meniscus) as a function of time. The transient diffusion equation is numerically solved by a Chebyshev collocation (spectral) method coupled to the Young-Laplace equation and the ideal gas law. We capture the effects of variable meniscus curvature and, subsequently, when applicable, movement of triple contact lines. Results are presented for the evolution of the dissolved gas concentration field in the liquid and, when applicable, the time it takes for a meniscus to depin and that for longevity, i.e., the onset of the Cassie to Wenzel state transition. Two configurations are examined; viz., one where an impermeable membrane pressurizes the liquid above the ridges and one where hydrostatic pressure is considered and the top of the liquid is exposed to non-condensable gas.more » « less
-
null (Ed.)Abstract We establish several finiteness properties of groups defined by algebraic difference equations. One of our main results is that a subgroup of the general linear group defined by possibly infinitely many algebraic difference equations in the matrix entries can indeed be defined by finitely many such equations. As an application, we show that the difference ideal of all difference algebraic relations among the solutions of a linear differential equation is finitely generated.more » « less
-
We present a reduced magnetohydrodynamic (MHD) mathematical model describing the dynamical behavior of highly conducting plasmas with frozen-in magnetic fields, constrained by the assumption that there exists a frame of reference, where the magnetic field vector, B, is aligned with the plasma velocity vector, u, at each point. We call this solution “stream-aligned MHD” (SA-MHD). Within the framework of this model, the electric field, E = − u × B ≡ 0, in the induction equation vanishes identically and so does the electromagnetic energy flux (Poynting flux), E × B ≡ 0, in the energy equation. At the same time, the force effect from the magnetic field on the plasma motion (the Ampère force) is fully taken into account in the momentum equation. Any steady-state solution of the proposed model is a legitimate solution of the full MHD system of equations. However, the converse statement is not true: in an arbitrary steady-state magnetic field, the electric field does not have to vanish identically (its curl has to, though). Specifically, realistic three-dimensional solutions for the steady-state (“ambient”) solar atmosphere in the form of so-called Parker spirals can be efficiently generated within the stream aligned MHD (SA-MHD) with no loss in generality.more » « less
-
ABSTRACT Rings and gaps are commonly observed in the dust continuum emission of young stellar discs. Previous studies have shown that substructures naturally develop in the weakly ionized gas of magnetized, non-ideal MHD discs. The gas rings are expected to trap large mm/cm-sized grains through pressure gradient-induced radial dust–gas drift. Using 2D (axisymmetric) MHD simulations that include ambipolar diffusion and dust grains of three representative sizes (1 mm, 3.3 mm, and 1 cm), we show that the grains indeed tend to drift radially relative to the gas towards the centres of the gas rings, at speeds much higher than in a smooth disc because of steeper pressure gradients. However, their spatial distribution is primarily controlled by meridional gas motions, which are typically much faster than the dust–gas drift. In particular, the grains that have settled near the mid-plane are carried rapidly inwards by a fast accretion stream to the inner edges of the gas rings, where they are lifted up by the gas flows diverted away from the mid-plane by a strong poloidal magnetic field. The flow pattern in our simulation provides an attractive explanation for the meridional flows recently inferred in HD 163296 and other discs, including both ‘collapsing’ regions where the gas near the disc surface converges towards the mid-plane and a disc wind. Our study highlights the prevalence of the potentially observable meridional flows associated with the gas substructure formation in non-ideal MHD discs and their crucial role in generating rings and gaps in dust.more » « less
