Convergence of the Critical Cooling Rate for Protoplanetary Disk Fragmentation Achieved: The Key Role of Numerical Dissipation of Angular Momentum
- Center for Theoretical Astrophysics and Cosmology, Institute for Computational Science, University of Zurich, Winterthurerstrasse 190, 8057 Zurich (Switzerland)
- Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CB3 0HA (United Kingdom)
We carry out simulations of gravitationally unstable disks using smoothed particle hydrodynamics (SPH) and the novel Lagrangian meshless finite mass (MFM) scheme in the GIZMO code. Our aim is to understand the cause of the nonconvergence of the cooling boundary for fragmentation reported in the literature. We run SPH simulations with two different artificial viscosity implementations and compare them with MFM, which does not employ any artificial viscosity. With MFM we demonstrate convergence of the critical cooling timescale for fragmentation at β{sub crit}≈3. Nonconvergence persists in SPH codes. We show how the nonconvergence problem is caused by artificial fragmentation triggered by excessive dissipation of angular momentum in domains with large velocity derivatives. With increased resolution, such domains become more prominent. Vorticity lags behind density, due to numerical viscous dissipation in these regions, promoting collapse with longer cooling times. Such effect is shown to be dominant over the competing tendency of artificial viscosity to diminish with increasing resolution. When the initial conditions are first relaxed for several orbits, the flow is more regular, with lower shear and vorticity in nonaxisymmetric regions, aiding convergence. Yet MFM is the only method that converges exactly. Our findings are of general interest, as numerical dissipation via artificial viscosity or advection errors can also occur in grid-based codes. Indeed, for the FARGO code values of β{sub crit} significantly higher than our converged estimate have been reported in the literature. Finally, we discuss implications for giant planet formation via disk instability.
- OSTI ID:
- 22875780
- Journal Information:
- Astrophysical Journal, Vol. 847, Issue 1; Other Information: Country of input: International Atomic Energy Agency (IAEA); ISSN 0004-637X
- Country of Publication:
- United States
- Language:
- English
Similar Records
Consistent SPH Simulations of Protostellar Collapse and Fragmentation
The Fragmentation Criteria in Local Vertically Stratified Self-gravitating Disk Simulations