Artificial viscosity in simulation of shock waves by smoothed particle hydrodynamics

This paper reconsiders artificial viscosity in smoothed particle hydrodynamics in order to prevent interparticle penetration, unwanted heating, and unphysical solutions. The coefficients in the Monaghan’s standard artificial viscosity are considered as a time variable, and we propose a restriction on them such that we avoid the undesired effects in the subsonic regions. We use the shock formation in adiabatic and isothermal cases to study the ability of this modified artificial viscosity recipe. Our computer experiments show that the proposal appears to work and the accuracy of this restriction is acceptable.     

Astrophysical Radiation Dynamics: The Prospects for Scaling

The general principles of scaling are discussed, followed by a survey of the important dimensionless parameters of fluid dynamics including radiation and magnetic fields, and of non-LTE spectroscopy. The values of the parameters are reviewed for a variety of astronomical and laboratory environments. It is found that parameters involving transport coefficients – the fluid and magnetic Reynolds numbers – have enormous values for the astronomical problems that are not reached in the lab. The parameters that measure the importance of radiation are also scarcely reached in the lab. This also means that the lab environments are much closer to LTE than the majority of astronomical examples. Some of the astronomical environments are more magnetically dominated than anything in the lab. The conclusion is that a good astronomical environment for simulation in a given lab experiment can be found, but that the reverse is much more difficult. PACS NOS: 95.30.Jx, 95.30.Lz, 97.10.Ex, 97.10.Gz, 98.62.Mw The U.S. Government’s right to retain a non-exclusive, royalty-free license in and to any copyright is acknowledged.     

Magnetic-Tower Jet Solution for Launching Astrophysical Jets

In spite of the large number of global three-dimensional (3-D) magnetohydrodynamic (MHD) simulations of accretion disks and astrophysical jets, which have been developed since 2000, the launching mechanisms of jets is somewhat controversial. Previous studies of jets have concentrated on the effect of the large-scale magnetic fields permeating accretion disks. However, the existence of such global magnetic fields is not evident in various astrophysical objects, and their origin is not well understood. Thus, we study the effect of small-scale magnetic fields confined within the accretion disk. We review our recent findings on the formation of jets in dynamo-active accretion disks by using 3-D MHD simulations. In our simulations, we found the emergence of accumulated azimuthal magnetic fields from the inner region of the disk (the so-called magnetic tower) and also the formation of a jet accelerated by the magnetic pressure of the tower. Our results indicate that the magnetic tower jet is one of the most promising mechanisms for launching jets from the magnetized accretion disk in various astrophysical objects. We will discuss the formation of cosmic jets in the context of the magnetic tower model.     

Numerical Simulations and Astrophysical Applications of Laboratory Jets at Omega

We have conducted experiments on the Omega laser at the University of Rochester that have produced jets of supersonic Ti impacting and being deflected by a ball of high density plastic. These mm-sized jets of dense plasma are highly complex, have large Reynolds numbers, and, given sufficient time and shear, should produce a fully turbulent flow. The experiments are diagnosed with a point-projection backlighter, resulting in a single image per shot. Simulations of the 3D hydrodynamics capture the large-scale features of the experimental data fairly well while missing some of the smaller scale turbulent-like phenomena. This is to be expected given the limited characterization of the targets as well as the finite resolution of the 3D simulations. If Euler scaling holds, these experiments should model larger astrophysical jets in objects such as HH 110 where an outflow can be seen colliding with a molecular cloud. However, Euler scaling demands that not only the isothermal internal Mach numbers of the two systems be similar but also that any dissipative mechanisms, such as radiative cooling or viscous dissipation, be of equal importance relative to each other. Similar equations of state are required as well. We discuss such issues in the context of these experiments and simulations.     

The importance of continuum radiation for the stellar wind hydrodynamics of hot stars

The contribution of bound-free and free-free processes to the outward acceleration of Pup is studied and is found to be negligible.     

Outflow Propagation in Collapsars: Collimated Jets and Expanding Outflows

We investigate the outflow propagation in the collapsar in the context of gamma-ray bursts (GRBs) with 2D relativistic hydrodynamic simulations. We vary the specific internal energy and bulk Lorentz factor of the injected outflow from non-relativistic regime to relativistic one, fixing the power of the outflow to be 10⁵¹erg s⁻¹. We observed the collimated outflow, when the Lorentz factor of the injected outflow is roughly greater than 2. To the contrary, when the velocity of the injected outflow is slower, the expanding outflow is observed. The transition from collimated jet to expanding outflow continuously occurs by decreasing the injected velocity. Different features of the dynamics of the outflows would cause the difference between the GRBs and similar phenomena, such as, X-ray flashes.     

Ionizing radiation in smoothed particle hydrodynamics

A new method for the inclusion of ionizing radiation from uniform radiation fields into 3D smoothed particle hydrodynamics (SPHI) simulations is presented. We calculate the optical depth for the Lyman continuum radiation from the source towards the SPHI particles by ray-tracing integration. The time-dependent ionization rate equation is then solved locally for the particles within the ionizing radiation field. Using test calculations, we explore the numerical behaviour of the code with respect to the implementation of the time-dependent ionization rate equation. We also test the coupling of the heating caused by the ionization to the hydrodynamical part of the SPHI code.     

Variational principles for relativistic smoothed particle hydrodynamics

In this paper we show how the equations of motion for the smoothed particle hydrodynamics (SPH) method may be derived from a variational principle for both non-relativistic and relativistic motion when there is no dissipation. Because the SPH density is a function of the coordinates the derivation of the equations of motion through variational principles is simpler than in the continuum case where the density is defined through the continuity equation. In particular, the derivation of the general relativistic equations is more direct and simpler than that of Fock. The symmetry properties of the Lagrangian lead immediately to the familiar additive conservation laws of linear and angular momentum and energy. In addition, we show that there is an approximately conserved quantity which, in the continuum limit, is the circulation.     

The Jeans instability in smoothed particle hydrodynamics

We derive analytically the Jeans criterion for a gas simulated using an SPH code in which the number of neighbours N _neighb is held constant (approximately) and the gravity-softening length, ε, equals the smoothing length, h (approximately). We show that the Jeans criterion is reproduced accurately for resolved structures, i.e. those represented by >  N _neighb particles. Unresolved structures are stabilized, as long as (i) the smoothing kernel W ( u ) is sufficiently centrally peaked, and (ii) the Jeans mass is resolved. Provided that these conditions are satisfied, then, in simulations of the formation of stars and galaxies, any fragmentation that occurs should be both physical and resolved. In particular there should be no creation of sub-Jeans condensations owing to numerical instability.     

Introducing a hybrid radiative transfer method for smoothed particle hydrodynamics

A new means of incorporating radiative transfer into smoothed particle hydrodynamics (SPH) is introduced, which builds on the success of two previous methods – the polytropic cooling approximation as devised by Stamatellos et al. and flux-limited diffusion. This hybrid method preserves the strengths of its individual components, while removing the need for atmosphere matching or other boundary conditions to marry optically thick and optically thin regions. The code uses a non-trivial equation of state to calculate temperatures and opacities of SPH particles, which captures the effects of H₂ dissociation, H⁰ ionization, He⁰ and He⁺ ionization, ice evaporation, dust sublimation, molecular absorption, bound-free and free–free transitions and electron scattering. The method is tested in several scenarios, including (i) the evolution of a  0.07 M_⊙  protoplanetary disc surrounding a  0.5 M_⊙  star; (ii) the collapse of a  1 M_⊙  protostellar cloud and (iii) the thermal relaxation of temperature fluctuations in a static homogeneous sphere.     

Three-dimensional smoothed particle hydrodynamics simulations of radiation-driven warped accretion discs

We present three-dimensional smoothed particle hydrodynamics calculations of warped accretion discs in X-ray binary systems. Geometrically thin, optically thick accretion discs are illuminated by a central radiation source. This illumination exerts a non-axisymmetric radiation pressure on the surface of the disc, resulting in a torque that acts on the disc to induce a twist or warp. Initially planar discs are unstable to warping driven by the radiation torque and, in general, the warps also precess in a retrograde direction relative to the orbital flow. We simulate a number of X-ray binary systems which have different mass ratios, using a number of different luminosities for each. Radiation-driven warping occurs for all systems simulated. For mass ratios   q ∼ 0.1  a moderate warp occurs in the inner disc while the outer disc remains in the orbital plane (cf. X 1916−053). For less extreme mass ratios, the entire disc tilts out of the orbital plane (cf. Her X–1). For discs that are tilted out of the orbital plane in which the outer edge material of the disc is precessing in a prograde direction, we obtain both positive and negative superhumps simultaneously in the dissipation light curve (cf. V603 Aql).