Free Access
Issue
A&A
Volume 631, November 2019
Article Number A66
Number of page(s) 6
Section Astrophysical processes
DOI https://doi.org/10.1051/0004-6361/201936215
Published online 22 October 2019

© ESO 2019

1. Introduction

In Marchand et al. (2018) (hereafter, Paper I), we presented the numerical implementation of the Hall effect in the AMR code RAMSES (Teyssier 2002), aimed for application in protostellar collapse simulations. While the implementation successfully passes several tests, showing the second-order convergence in space, the gas angular momentum is not conserved in star formation simulations. As expected, the Hall effect generates rotation in an initially non-rotating cloud and counter-rotating envelopes form on both sides of the mid-plane (Krasnopolsky et al. 2011; Tsukamoto et al. 2015, 2017; Wurster et al. 2017). However, shortly after the formation of the first Larson core, a large amount of rotation is generated in the first core and violates the conservation of the total angular momentum, which is a purely numerical issue. This problem arises in every simulation with the Hall effect, severely limiting the validity of our results. We could not find the origin of the problem. Krasnopolsky et al. (2011) encountered a similar issue and assumed it was due to their boundary conditions. We however did not find significant angular momentum transport through the box boundaries. In this work, as well as in Paper I, we only consider the angular momentum of the fluid, not the magnetic field component, because there is no transfer between both in our framework.

In this paper, we present a minor modification of our numerical scheme that significantly decreases the spurious generation of angular momentum. This method is presented in Sect. 2, then its impact on the test case of Paper I in Sect. 3. In Sect. 4, we make a comparison with a previous study of a more realistic scenario (initial rotation and consistent resistivities), and Sects. 5 and 6 are dedicated to discussion and conclusions.

2. Methods

2.1. The Hall effect implementation

The implementation of the Hall effect in the eulerian RAMSES code, as described in details in Paper I, has been inspired by Lesur et al. (2014). We briefly summarize it here.

The magnetic field is updated at every time-step using the Constrained Transport scheme (Evans & Hawley 1988) on cell interfaces. For the x-component, the integration reads

B x , i 1 2 , j , k n + 1 B x , i 1 2 , j , k n Δ t = E z , i 1 2 , j + 1 2 , k n + 1 2 E z , i 1 2 , j 1 2 , k n + 1 2 Δ x E y , i 1 2 , j , k + 1 2 n + 1 2 E y , i 1 2 , j , k 1 2 n + 1 2 Δ x . $$ \begin{aligned} \frac{B^{n+1}_{x,i-\frac{1}{2},j,k}-B^{n}_{x,i-\frac{1}{2},j,k}}{\Delta t}&= \frac{E^{n+\frac{1}{2}}_{z,i-\frac{1}{2},j+\frac{1}{2},k}-E^{n+\frac{1}{2}}_{z,i-\frac{1}{2},j-\frac{1}{2},k}}{\Delta x} \nonumber \\&\quad - \frac{E^{n+\frac{1}{2}}_{y,i-\frac{1}{2},j,k+\frac{1}{2}}-E^{n+\frac{1}{2}}_{y,i-\frac{1}{2},j,k-\frac{1}{2}}}{\Delta x}. \end{aligned} $$(1)

B x , i 1 2 , j , k n + 1 $ B^{n+1}_{x,i-\frac{1}{2},j,k} $ is the x-component of the magnetic field at time-step n + 1 on cell face (i − 1/2, j, k), with (i, j, k)i, j, k ∈ ℕ being the indexes of the center of the cell. E , i 1 2 , , n + 1 2 $ E^{n+\frac{1}{2}}_{*,i-\frac{1}{2},*,*} $ are the electric fields on cell edges drawing the contour of the cell interface. These electric fields are computed after the prediction step of the Monotone Upstream-centered Scheme for Conservative Laws (MUSCL) scheme (van Leer 1976) using the Harten Lax-Van Leer (HLL) 2D-Riemann solver (Londrillo & del Zanna 2004)

E z n + 1 2 = S L S B E z n + 1 2 , LB S L S T E z n + 1 2 , LT S R S B E z n + 1 2 , RB + S R S T E z n + 1 2 , RT ( S R S L ) ( S T S B ) S T S B S T S B ( B x T B x B ) + S R S L S R S L ( B y R B y L ) . $$ \begin{aligned}&E^{n+\frac{1}{2}}_{z} \nonumber \\&\quad =\frac{S_{\rm L}S_{\rm B}E^{n+\frac{1}{2},\mathrm{LB}}_{z}-S_{\rm L}S_{\rm T}E^{n+\frac{1}{2},\mathrm{LT}}_{z}-S_{\rm R}S_{\rm B}E^{n+\frac{1}{2},\mathrm{RB}}_{z}+S_{\rm R}S_{\rm T}E^{n+\frac{1}{2},\mathrm{RT}}_{z}}{(S_{\rm R}-S_{\rm L})(S_{\rm T}-S_{\rm B})} \nonumber \\&\qquad -\frac{S_{\rm T}S_{\rm B}}{S_{\rm T}-S_{\rm B}}(B^\mathrm{T}_{x} - B^\mathrm{B}_{x})+\frac{S_{\rm R}S_{\rm L}}{S_{\rm R}-S_{\rm L}}(B^\mathrm{R}_{y} - B^\mathrm{L}_{y}). \end{aligned} $$(2)

Here, E n + 1 2 , ( LB , LT , RB , RT ) $ E^{n+\frac{1}{2},(\mathrm{LB,LT,RB,RT})} $ represents the electric fields at the corners of the 4 cells adjacent to the edge1, B(L, B, R, T) are the magnetic fields at the cell corners averaged on the interfaces, and S(L, B, R, T) are the minimum and maximum characteristic wave speeds of the system in the 2 directions perpendicular to the edge. Figure 3 of Paper I summarizes all the notations.

The electric fields are computed by adding the flux-freezing electric field of ideal magnetohydrodynamics (MHD) u × B and the Hall electric field uH × B. uH = −ηHJ/||B|| is the Hall speed and is averaged over the 4 cells, ηH is the Hall resistivity and J = ∇ × B is the electric field. Hence, for, e.g., corner LB,

E z n + 1 2 , LB = ( u x n + 1 2 , LB + u H , x ) B y L ( u y n + 1 2 , LB + u H , y ) B x B . $$ \begin{aligned} E^{n+\frac{1}{2},\mathrm{LB}}_z = \left(u^{n+\frac{1}{2},\mathrm{LB}}_x + u_{\mathrm{H} ,x}\right) B^\mathrm{L}_y - (u^{n+\frac{1}{2},\mathrm{LB}}_y + u_{\mathrm{H} ,y}) B^\mathrm{B}_x. \end{aligned} $$(3)

The very last sentence of Sect. 3.5 in Paper I states that the whistler waves speeds are accounted for in the 1D-Riemann problems at cell interfaces during the prediction step of the MUSCL scheme. Actually, that sentence is incorrect, as the prediction step does not use the Riemann solvers. The intended meaning is that whistler speeds are added to the wave fan of other variables (density, momentum ...), that use a 1D Riemann solver to compute the fluxes at cell interfaces in the correction step. The wave fan should be the same for every variable because the (characteristic) waves speeds correspond to the eigenvalues of the system of equation. The 1D HLL (Harten et al. 1983) flux for the non-magnetic variables U reads

F HLL n + 1 2 = S R F L n + 1 2 S L F R n + 1 2 + S L S R ( U R n + 1 2 U L n + 1 2 ) S R S L , $$ \begin{aligned} F_{\rm HLL}^{n+\frac{1}{2}} = \frac{S_{\rm R} F_{\rm L}^{n+\frac{1}{2}} - S_{\rm L} F_{\rm R}^{n+\frac{1}{2}} + S_{\rm L} S_{\rm R} \left(U_{\rm R}^{n+\frac{1}{2}}-U_{\rm L}^{n+\frac{1}{2}}\right)}{S_{\rm R} S_{\rm L}}, \end{aligned} $$(4)

with subscripts L and R indicating the left and right side of the interface. The flux is then used in a second-order Godunov scheme to update the flow variables.

The speed of the whistler waves is

c w = η H π 2 Δ x + ( η H π 2 Δ x ) 2 + c A 2 , $$ \begin{aligned} c_\mathrm{w} = \frac{\eta _\mathrm{H} \pi }{2\Delta x} + \sqrt{\left(\frac{\eta _\mathrm{H} \pi }{2\Delta x}\right)^2 + c_\mathrm{A} ^2}, \end{aligned} $$(5)

with c A = B / ρ $ c_{\mathrm{A}}=B/\sqrt{\rho} $ the Alfven speed, ρ the density and Δx the cell-size. For all variables, the characteristic speeds used in Eqs. (2) and (4) account for the whistler speeds. Since they are usually the fastest waves at high resolution, we have SL = ux − cw, SR = ux + cw, SB = uy − cw and ST = uy + cw.

2.2. Modification of the scheme

The modification we propose consists in not accounting for whistler waves in the Riemann problems of variables other than the magnetic field. In other words, we do not use them to compute SL and SR in Eq. (4). Two reasons motivate this modification. Firstly, the truncation error increases with the characteristic speeds. At the center of a protostellar collapse simulation, the whistler speed can reach several hundred times the value of the second fastest wave, the fast magneto-sonic wave, which is usually the fastest wave in the absence of the Hall effect. In the test case presented in Paper I, there is a factor ∼300 between both speeds at the end of the protostellar collapse simulation. Truncation errors increase then significantly in the first Larson core even for the purely hydro variable, the momentum in particular. The second reason is that the Hall effect does not directly affect the fluid motion but only indirectly through the Lorentz force (by changing the magnetic fields). Removing the whistler speed in the magnetic Riemann problems as well leads to a magnetic field instability.

In Appendix A, we show that our new scheme still propagates the whistler waves at the correct frequencies and keeps its second order convergence in space.

3. Test case: non-rotating sphere

3.1. Initial conditions

We use the same initial condition and numerical parameters as in Sect. 5 of Paper I, a uniform 1.5 M non-rotating sphere of radius 3712 au at T = 10 K with a uniform magnetic field B = 90.3 μG. The temperature is given by the following equation of state (EOS)

T = 10 [ 1 + ( ρ 10 13 g cm 3 ) 2 3 ] . $$ \begin{aligned} T=10 \left[ 1+\left(\frac{\rho }{10^{-13}\,\mathrm{g} \,\mathrm{cm} ^{-3}} \right)^{\frac{2}{3}} \right] . \end{aligned} $$(6)

We include only the Hall effect with a constant resistivity ηH = 1020 cm2 s−1. We use the generalized monotonized central (moncen) slope limiter with a coefficient 1.5 as in Paper I, with two different refinement criteria, 8 and 16 points per Jeans length respectively. We also perform one simulation with a moncen coefficient of 1.05 and 8 points per Jeans length for consistency with Sect. 4.

3.2. Results

Our results are qualitatively similar to Paper I. The Hall effect generates rotation in the cloud, especially in the mid-plane, and counter-rotating envelopes develop above and below the mid-plane to compensate the generation of angular momentum. Figure 1 represents a slice at x = 0 of the azimuthal velocity with magnetic field vectors. The scale and the shape of the counter-rotating envelopes are similar to the collapse in Paper I.

thumbnail Fig. 1.

Density slice at x = 0 of the azimuthal velocity at t = tC + 700 yr. Yellow-red colors represent rotation in the positive azimuthal direction, while blue-purple colors represent the negative direction. Gray vectors indicate the direction of the magnetic field.

The “positive” and “negative” angular momenta in the simulations ( L z + $ L_z^ + $ and L z $ L_z^ - $ as defined in Paper I) are plotted in Fig. 2 alongside the reference case of Paper I. There is an obvious improvement of the conservation of the total angular momentum with the present method. The divergence starts at tC + 700 yr instead of tC + 300 yr, and the increase of L z + $ L_z^ + $ is one order of magnitude slower than in the previous case. The conservation of the lower resolution case is even better than the higher resolution case of Paper I for both slope limiters. If we compare the excess of angular momentum (difference between solid and dashed lines) to the total angular momentum of the disk (dotted lines), there is also a significant improvement. In the previous method, 85% of the disk’s angular momentum was due to the numerical error. With the new scheme, this fraction decreases to less than 50%. Additionally, both resolutions and both slope limiters show the same disk’s angular momentum.

thumbnail Fig. 2.

Evolution of the positive and negative angular momentum ( L z + $ L_z^ + $ and L z $ L_z^ - $) in solid and dashed line respectively. Dotted lines indicate the angular momentum of the disk. The blue lines represent the reference case of Paper I, while the red and green lines show the evolution for the modified scheme, for a moncen coefficient of 1.5 and 1.05 respectively. The lower resolution simulations with 8 points per Jeans length are in light colors and higher resolution simulations with 16 points per Jeans length are in darker colors.

In this test-case, the initial angular momentum is zero, while in reality dense cores exhibit rotational motions (Goodman et al. 1993). Moreover, a realistic Hall resistivity, as computed in Marchand et al. (2016), is likely to be one order of magnitude lower than in this case. The creation of spurious angular momentum would be then damped by this factor approximatively. For these reasons, and the addition of ambipolar and Ohmic diffusion, we can expect the angular momentum increase to be negligible compared to the total angular momentum, and even the accretion of angular momentum in the disk, in a more realistic setup. This point is tested in next section.

4. Models with initial rotation

4.1. Initial conditions

In the following simulations, we use the same initial conditions as Tsukamoto et al. (2015; hereafter T15). The purpose is to assess the validity of our numerical methods in realistic conditions by comparison with an independent implementation. The initial cloud is a uniform sphere of M = 1 M, with a radius R = 2957 au and a temperature of 10 K (thermal-to-gravitational energy ratio of α = 0.3). The sphere undergoes a solid rotation characterized by a rotational-to-gravitational energy ratio of β = 0.011. The initial magnetic field is uniform, either parallel (θ = 0°) or anti-parallel (θ = 180°) to the rotation axis, with a mass-to-flux ratio of μ = 4. All three non-ideal MHD effects are included. We also perform another simulation without the Hall effect for comparison. The three cases are summarized in Table 1. Contrary to T15, we do not include the effects of radiation-hydrodynamics, and instead use the barotropic EOS (6) to compute the temperature. While the first core is not able to do its second collapse with such EOS, the temperature rise prevents high densities to be reached quickly, which would slow the simulation. We are here interested in the formation of structures in and around the first core and the disk rather than the formation of a protostar. We define the disk as a rotationally-supported structure, with the same criteria as in Joos et al. (2012) and Paper I. Contrary to T15, the gas outside the sphere is at rest, with a density 30 times lower than the sphere, and boundary conditions are periodic. The refinement criterion is 8 points per Jeans length.

Table 1.

Simulations.

Simulations have been performed with the generalized moncen slope limiter with a coefficient 1.05 for the magnetic field and minmod for the other variables. We use the shallow slope limiters to prevent any overshooting of magnetic field while reconstructing states at cell interfaces. In the following sections, tC denotes the formation time of the first Larson core, that is when the maximum density reaches 10−13 g cm−3.

4.2. Magnetic resistivities

To compute the resistivities of non-ideal MHD effects, we use the table of Marchand et al. (2016), which contains the equilibrium abundances of a reduced chemical network across a wide range of density and temperature. The network includes species relevant to the star formation environment and grains following the Mathis, Rumpl and Nordseick (MRN) size distribution (Mathis et al. 1977). They take into account thermal ionisations, the thermionic emission of grains (Desch & Turner 2015) and the grain evaporation. During simulations, the non-ideal MHD resistivities are calculated for each cell using the local state variables.

4.3. Results

We choose to stop the slower simulation, Th0, at t = tf = tc + 1400 yr, and the two faster simulations at t = tc + 3000 yr. At the 1400 years mark, the maximum density is ρmax = 10−10 g cm−3 for Th0, ρmax = 10−11 g cm−3 for Th180 and ρmax = 6.6 × 10−11 g cm−3 for NoHall. Figure 3 displays density maps of the mid-plane for the three simulations at this time. As in T15, the Hall effect reduces the magnetic braking in the anti-parallel case, allowing the formation of a large disk, and enhances it in the parallel case, speeding up the collapse. The Hall effect modifies the size of the disk by up to 50% in this setup, resulting in a factor 2 between the models with parallel and anti-parallel magnetic fields. The disks in Th180 and Nohall develop a m = 2 instability at t = tc + 2000 yr and t = tc + 3000 yr, respectively.

thumbnail Fig. 3.

Density slices of the mid-plane for (from left to right) Th0, NoHall, and Th180. Black vectors represent the direction of the gas velocity.

Figure 4 represents the azimuthal velocity maps with the disk seen from the edge for Th180, as Fig. 5 of T15, at t = tC + 1400 yr and t = tC + 300 yr. The global rotation (in the negative direction in this simulation) is reinforced in the mid-plane by the Hall effect, leading to a high rotation velocity, up to uϕ ≈ 2 km s−1. Counter-rotating envelopes form in this case, with velocity of uϕ ≈ 1 km s−1. The scale and geometry of the various rotating regions match almost perfectly the Fig. 5 of T15. This geometry evolves over time. The infalling negatively-rotating gas eventually mixes with the positively-rotating envelopes, that almost disappear at the end of the simulation, 1.5 years after this snapshot (see bottom panel of Fig. 4), for a total lifetime of ≲3000 years in this case.

thumbnail Fig. 4.

Azimuthal velocity slice in the plane y = 0 (edge-on) for Th180. Blue-purple colors represent the rotation in the negative direction, and red-yellow colors the rotation in the positive direction. Top panel: t − tC = 1.4 kyr, bottom panel: t − tC = 2.9 kyr.

We now compare the angular momentum in the various cases. Figure 5 shows the angular momentum contained in the core (ρ >  10−13 g cm−3) and the disk. Unsurprisingly, the run in which the Hall effect enhances (reduces) the magnetic braking has the lowest (highest) angular momentum, with a factor two between both, and NoHall being an intermediate case. The solid line is the total angular momentum in the simulation box, shifted such as it equals zero at t = tc. It then represents the accumulated error on the total angular momentum as function of the time after the first core formation. In all three cases, this quantity increases at similar rates. After 1400 years, the excess reaches ≈2−3 × 1051 g cm2 s−1, which is less than 1% of the total angular momentum Ltot ≈ 3.7 × 1053 g cm2 s−1. At t = tc + 3000 yr, this error represents ≈2.7% of the total angular momentum. The relative value between the dashed and solid lines gives the upper-limit of the fraction of the disk’s angular momentum that can be attributed to the spurious error linked to the Hall effect. This fraction is 30% in Th180 and 80% in Th0. However, given that the “excess” of angular momentum in NoHall is similar to Th0 and Th180 (and even higher in this case), the spurious fraction due to the Hall effect is most likely minor.

thumbnail Fig. 5.

Angular momentum evolution after the first core formation for Th0 (red curves), Th180 (blue curves) and Nohall (green curves). The solid lines represent the “excess” of angular momentum in the whole box compared to t = tC.

5. Discussions

In the simulations with initial rotation, the angular momentum excess due to the Hall effect is now a minor factor with respect to the total angular momentum, but our results show that conservation is still not perfect, independently of the Hall effect. This can however be significantly improved by simply increasing the resolution to 16 points per Jeans length, which remains computationally affordable. Also, the additional thermal support provided by the stiff EOS keeps the gas away from the center. Any numerical error on the linear momentum therefore results in a larger error on the angular momentum compared to a case with a more compact core.

A quantitative comparison to T15 is difficult here because we use a stiff barotropic EOS instead of radiation-hydrodynamics, and our chemical models to compute the resistivities are different. We also did not evolve the simulation for such a long time. Though, we reproduce very well their results, especially concerning the size of the disk and the geometry of the various rotating regions.

Counter-rotating envelopes of 200 au scale develop on each side of the disk, but eventually disappear in 3 kyr. While such structures could be observed by modern instruments, as was claimed by Takakuwa et al. (2018), their short lifetime makes their discovery extremely unlikely. Different magnetic field inclinations can produce more prominent envelopes, especially θ = 90° and θ = 135° (Tsukamoto et al. 2017). We however do not expect them to live significantly longer. Should one be detected, we would expect to find it in the vicinity of a first hydrostatic core.

6. Conclusions

The angular momentum conservation has clearly been improved by the minor modification of our numerical scheme. While the issue has not completely disappeared, the spurious generation of angular momentum has been reduced by one order of magnitude. It is now negligible in regard to the disk’s angular momentum of a simulation with initial rotation and with realistic resistivities. The origin of the problem is still unknown, but we can now be confident in the results of our simulations with the Hall effect if we have sufficiently high resolution.


1

L,B,R,T stand for Left, Bottom, Right, Top.

References

  1. Desch, S. J., & Turner, N. J. 2015, ApJ, 811, 156 [NASA ADS] [CrossRef] [Google Scholar]
  2. Evans, C. R., & Hawley, J. F. 1988, ApJ, 332, 659 [Google Scholar]
  3. Goodman, A. A., Benson, P. J., Fuller, G. A., & Myers, P. C. 1993, ApJ, 406, 528 [NASA ADS] [CrossRef] [Google Scholar]
  4. Harten, A., Lax, P., & Van Leer, B. 1983, SIAM Rev., 25, 35 [Google Scholar]
  5. Joos, M., Hennebelle, P., & Ciardi, A. 2012, A&A, 543, A128 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Krasnopolsky, R., Li, Z.-Y., & Shang, H. 2011, ApJ, 733, 54 [NASA ADS] [CrossRef] [Google Scholar]
  7. Lesur, G., Kunz, M. W., & Fromang, S. 2014, A&A, 566, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Londrillo, P., & del Zanna, L. 2004, J. Comput. Phys., 195, 17 [NASA ADS] [CrossRef] [Google Scholar]
  9. Marchand, P., Masson, J., Chabrier, G., et al. 2016, A&A, 592, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Marchand, P., Commeron, B., & Chabrier, G. 2018, A&A, 619, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425 [NASA ADS] [CrossRef] [Google Scholar]
  12. Sano, T., & Stone, J. M. 2002, ApJ, 570, 314 [NASA ADS] [CrossRef] [Google Scholar]
  13. Takakuwa, S., Tsukamoto, Y., Saigo, K., & Saito, M. 2018, ApJ, 865, 51 [NASA ADS] [CrossRef] [Google Scholar]
  14. Teyssier, R. 2002, A&A, 385, 337 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Tsukamoto, Y., Iwasaki, K., Okuzumi, S., Machida, M. N., & Inutsuka, S. 2015, ApJ, 810, L26 [NASA ADS] [CrossRef] [Google Scholar]
  16. Tsukamoto, Y., Okuzumi, S., Iwasaki, K., Machida, M. N., & Inutsuka, S.-I. 2017, PASJ, 69, 95 [CrossRef] [Google Scholar]
  17. van Leer, B. 1976, in Computing in Plasma Physics and Astrophysics, ed. D. Biskamp (Amsterdam: Elsevier) [Google Scholar]
  18. Wurster, J., Price, D. J., & Bate, M. R. 2017, MNRAS, 466, 1788 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Scheme convergence

Our new scheme successfully passes all tests described in Paper I. We show here the tests described in Sect. 4.1 (Sano & Stone 2002). In a periodic box with a uniform density and pressure, we set up a magnetic wave B = Bxex + By(x)ez, where Bx = 0.1 G and By(x) = 10−3 cos (2πkx) G, with k the wave number (and the number of periods in the box). Only the Hall effect is present with a constant resistivity ηH. First, we use k and ηH as parameters to find the dispersion relation as in Sect. 4.1.1 of Paper I. k ranges from 5 to 20, and ηH ranges from 5 × 10−3 to 0.1 cm2 s−1. The resolution of the box is uniform with 1283 cells. We use the same slope limiters as in Sect. 4, namely generalized moncen with a coefficient 1.05 for the magnetic field, and minmod for the other hydrodynamic variables. The top panel of Fig. A.1 represents the results, with the normalized frequency of propagation ω/ωH, ω H = c A 2 / η H $ {\omega _{\rm{H}}} = c_{\rm{A}}^2/{\eta _{\rm{H}}} $ as a function of klH, with lH = ηH/cA. We recover the dispersion relation of the Hall effect (Eq. (18) of Paper I) with an error of less than 5% for most points. For the convergence test, we fix k = 5, ηH = 0.005 cm2 s−1 and we use resolutions of 323, 643 and 1283. While propagating, the wave is dissipated by the numerical diffusion (after 5 periods for the low resolution case), and we compute the damping rate for each case, as displayed in the bottom panel of Fig. A.1. It shows that the new method keeps the second order spatial convergence of the numerical damping.

thumbnail Fig. A.1.

Top panel: dispersion relation of the numerical scheme. Red points are the measured oscillation frequencies normalised by ω H = c A 2 / η H $ {\omega _{\rm{H}}} = c_{\rm{A}}^2/{\eta _{\rm{H}}} $ as a function of klH with lH = ηH/cA. Solid black lines represent the theoretical dispersion relation, while blue points show the relative error of each measure. Bottom panel: damping rate of whistler waves as a function of resolution, compared with Δx and Δx2 scaling.

All Tables

Table 1.

Simulations.

All Figures

thumbnail Fig. 1.

Density slice at x = 0 of the azimuthal velocity at t = tC + 700 yr. Yellow-red colors represent rotation in the positive azimuthal direction, while blue-purple colors represent the negative direction. Gray vectors indicate the direction of the magnetic field.

In the text
thumbnail Fig. 2.

Evolution of the positive and negative angular momentum ( L z + $ L_z^ + $ and L z $ L_z^ - $) in solid and dashed line respectively. Dotted lines indicate the angular momentum of the disk. The blue lines represent the reference case of Paper I, while the red and green lines show the evolution for the modified scheme, for a moncen coefficient of 1.5 and 1.05 respectively. The lower resolution simulations with 8 points per Jeans length are in light colors and higher resolution simulations with 16 points per Jeans length are in darker colors.

In the text
thumbnail Fig. 3.

Density slices of the mid-plane for (from left to right) Th0, NoHall, and Th180. Black vectors represent the direction of the gas velocity.

In the text
thumbnail Fig. 4.

Azimuthal velocity slice in the plane y = 0 (edge-on) for Th180. Blue-purple colors represent the rotation in the negative direction, and red-yellow colors the rotation in the positive direction. Top panel: t − tC = 1.4 kyr, bottom panel: t − tC = 2.9 kyr.

In the text
thumbnail Fig. 5.

Angular momentum evolution after the first core formation for Th0 (red curves), Th180 (blue curves) and Nohall (green curves). The solid lines represent the “excess” of angular momentum in the whole box compared to t = tC.

In the text
thumbnail Fig. A.1.

Top panel: dispersion relation of the numerical scheme. Red points are the measured oscillation frequencies normalised by ω H = c A 2 / η H $ {\omega _{\rm{H}}} = c_{\rm{A}}^2/{\eta _{\rm{H}}} $ as a function of klH with lH = ηH/cA. Solid black lines represent the theoretical dispersion relation, while blue points show the relative error of each measure. Bottom panel: damping rate of whistler waves as a function of resolution, compared with Δx and Δx2 scaling.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.