Aims. The importance of radiation to the physical structure of protoplanetary disks cannot be understated. However, protoplanetary disks evolve with time, and so to understand disk evolution and by association, disk structure, one should solve the co
mbined and time-dependent equations of radiation hydrodynamics. Methods. We implement a new implicit radiation solver in the AZEuS adaptive mesh refinement magnetohydrodynamics fluid code. Based on a hybrid approach that combines frequency-dependent ray-tracing for stellar irradiation with non-equilibrium flux limited diffusion, we solve the equations of radiation hydrodynamics while preserving the directionality of the stellar irradiation. The implementation permits simulations in Cartesian, cylindrical, and spherical coordinates, on both uniform and adaptive grids. Results. We present several hydrostatic and hydrodynamic radiation tests which validate our implementation on uniform and adaptive grids as appropriate, including benchmarks specifically designed for protoplanetary disks. Our results demonstrate that the combination of a hybrid radiation algorithm with AZEuS is an effective tool for radiation hydrodynamics studies, and produces results which are competitive with other astrophysical radiation hydrodynamics codes.
Context. The growth process of dust particles in protoplanetary disks can be modeled via numerical dust coagulation codes. In this approach, physical effects that dominate the dust growth process often must be implemented in a parameterized form. Due
to a lack of these parameterizations, existing studies of dust coagulation have ignored the effects a hydrodynamical gas flow can have on grain growth, even though it is often argued that the flow could significantly contribute either positively or negatively to the growth process. Aims. We intend to provide a quantification of hydrodynamical effects on the growth of dust particles, such that these effects can be parameterized and implemented in a dust coagulation code. Methods. We numerically integrate the trajectories of small dust particles in the flow of disk gas around a proto-planetesimal, sampling a large parameter space in proto-planetesimal radii, headwind velocities, and dust stopping times. Results. The gas flow deflects most particles away from the proto-planetesimal, such that its effective collisional cross section, and therefore the mass accretion rate, is reduced. The gas flow however also reduces the impact velocity of small dust particles onto a proto-planetesimal. This can be beneficial for its growth, since large impact velocities are known to lead to erosion. We also demonstrate why such a gas flow does not return collisional debris to the surface of a proto-planetesimal. Conclusions. We predict that a laminar hydrodynamical flow around a proto-planetesimal will have a significant effect on its growth. However, we cannot easily predict which result, the reduction of the impact velocity or the sweep-up cross section, will be more important. Therefore, we provide parameterizations ready for implementation into a dust coagulation code.
A new adaptive mesh refinement (AMR) version of the ZEUS-3D astrophysical magnetohydrodynamical (MHD) fluid code, AZEuS, is described. The AMR module in AZEuS has been completely adapted to the staggered mesh that characterises the ZEUS family of cod
es, on which scalar quantities are zone-centred and vector components are face-centred. In addition, for applications using static grids, it is necessary to use higher-order interpolations for prolongation to minimise the errors caused by waves crossing from a grid of one resolution to another. Finally, solutions to test problems in 1-, 2-, and 3-dimensions in both Cartesian and spherical coordinates are presented.
