Light curves of SN 1998A and comparison with similar unusual SNe

Following the proposal by Damineli that the central object of Eta Carinae may be an early-type binary, we perform numerical simulations of the X-ray emission from colliding stellar winds. A synthetic light curve has been generated which qualitatively agrees with the recent X-ray variability, and provides further support for the binary model. In particular, the model predicts a rise in the observed X-ray emission towards periastron, followed by a sharp drop and subsequent recovery. This is indeed what is seen in the RXTE light curve, although some problems concerning the X-ray spectrum at periastron still need to be explained. The simulations suggest that the width of the periastron dip will provide strong constraints on the binary and stellar wind properties of the components of Eta Car.     

Non-axisymmetric relativistic Bondi–Hoyle accretion on to a Schwarzschild black hole

We present the results of an exhaustive numerical study of fully relativistic non-axisymmetric Bondi–Hoyle accretion on to a moving Schwarzschild black hole. We have solved the equations of general relativistic hydrodynamics with a high-resolution shock-capturing numerical scheme based on a linearized Riemann solver. The numerical code was previously used to study axisymmetric flow configurations past a Schwarzschild black hole. We have analysed and discussed the flow morphology for a sample of asymptotically high Mach number models. The results of this work reveal that initially asymptotic uniform flows always accrete on to the hole in a stationary way, which closely resembles the previous axisymmetric patterns. This is in contrast with some Newtonian numerical studies where violent flip-flop instabilities were found. As discussed in the text, the reason can be found in the initial conditions used in the relativistic regime, as they cannot exactly duplicate the previous Newtonian setups where the instability appeared. The dependence of the final solution on the inner boundary condition as well as on the grid resolution has also been studied. Finally, we have computed the accretion rates of mass and linear and angular momentum.     

Radiatively driven mass accretion on to galactic nuclei by circumnuclear starbursts

We examine the physical processes of radiatively driven mass accretion on to galactic nuclei, owing to intensive radiation from circumnuclear starbursts. The radiation from a starburst not only causes the inner gas disc to contract via radition flux force, but also extracts angular momentum owing to relativistic radiation drag, thereby inducing an avalanche of the surface layer of the disc. To analyse such a mechanism, the radiation–hydrodynamical equations are solved, including the effects of the radiation drag force as well as the radiation flux force. As a result, it is found that the mass accretion rate owing to the radiative avalanche is given by M ˙ ( r )= η ( L _*/ c ²)( r / R )² (Δ R / R )(1 −  e ^−τ) at radius r , where the efficiency η ranges from 0.2 up to 1, L _* and R are respectively the bolometric luminosity and the radius of the starburst ring, Δ R is the extent of the emission regions, and τ is the face-on optical depth of the disc. In an optically thick regime, the rate depends upon neither the optical depth nor the surface mass density distribution of the disc. The present radiatively driven mass accretion may provide a physical mechanism which enables mass accretion from 100-pc scales down to ∼ parsec scales, and it may eventually be linked to advection-dominated viscous accretion on to a massive black hole. The radiation–hydrodynamical and self-gravitational instabilities of the disc are briefly discussed. In particular, the radiative acceleration possibly builds up a dusty wall, which 'shades' the nucleus in edge-on views. This provides another version of the model for the formation of an obscuring torus.     

Lagrangian meshfree particles method (SPH) for large deformation and failure flows of geomaterial using elastic–plastic soil constitutive model

Simulation of large deformation and post‐failure of geomaterial in the framework of smoothed particle hydrodynamics (SPH) are presented in this study. The Drucker–Prager model with associated and non‐associated plastic flow rules is implemented into the SPH code to describe elastic–plastic soil behavior. In contrast to previous work on SPH for solids, where the hydrostatic pressure is often estimated from density by an equation of state, this study proposes to calculate the hydrostatic pressure of soil directly from constitutive models. Results obtained in this paper show that the original SPH method, which has been successfully applied to a vast range of problems, is unable to directly solve elastic–plastic flows of soil because of the so‐called SPH tensile instability. This numerical instability may result in unrealistic fracture and particles clustering in SPH simulation. For non‐cohesive soil, the instability is not serious and can be completely removed by using a tension cracking treatment from soil constitutive model and thereby give realistic soil behavior. However, the serious tensile instability that is found in SPH application for cohesive soil requires a special treatment to overcome this problem. In this paper, an artificial stress method is applied to remove the SPH numerical instability in cohesive soil. A number of numerical tests are carried out to check the capability of SPH in the current application. Numerical results are then compared with experimental and finite element method solutions. The good agreement obtained from these comparisons suggests that SPH can be extended to general geotechnical problems. Copyright © 2008 John Wiley & Sons, Ltd.     

AMRVAC and relativistic hydrodynamic simulations for gamma-ray burst afterglow phases

We apply a novel adaptive mesh refinement (AMR) code, AMRVAC (Adaptive Mesh Refinement version of the Versatile Advection Code), to numerically investigate the various evolutionary phases in the interaction of a relativistic shell with its surrounding cold interstellar medium (ISM). We do this for both 1D isotropic and full 2D jet-like fireball models. This is relevant for gamma-ray bursts (GRBs), and we demonstrate that, thanks to the AMR strategy, we resolve the internal structure of the shocked shell–ISM matter, which will leave its imprint on the GRB afterglow. We determine the deceleration from an initial Lorentz factor  γ= 100  up to the almost Newtonian     phase of the flow. We present axisymmetric 2D shell evolutions, with the 2D extent characterized by their initial opening angle. In such jet-like GRB models, we discuss the differences with the 1D isotropic GRB equivalents. These are mainly due to thermally induced sideways expansions of both the shocked shell and shocked ISM regions. We found that the propagating 2D ultrarelativistic shell does not accrete all the surrounding medium located within its initial opening angle. Part of this ISM matter gets pushed away laterally and forms a wide bow-shock configuration with swirling flow patterns trailing the thin shell. The resulting shell deceleration is quite different from that found in isotropic GRB models. As long as the lateral shell expansion is merely due to ballistic spreading of the shell, isotropic and 2D models agree perfectly. As thermally induced expansions eventually lead to significantly higher lateral speeds, the 2D shell interacts with comparably more ISM matter and decelerates earlier than its isotropic counterpart.     

Shock oscillation model for quasi-periodic oscillations in stellar mass and supermassive black holes

We numerically examine centrifugally supported shock waves in 2D rotating accretion flows around a stellar mass  (10 M_⊙)  and a supermassive  (10⁶ M_⊙)  black holes over a wide range of input accretion rates of     . The resultant 2D shocks are unstable with time and the luminosities show quasi-periodic oscillations (QPOs) with modulations of a factor of 2–3 and with periods of a tenth of a second to several hours, depending on the black hole masses. The shock oscillation model may explain the intermediate frequency QPOs with 1–10 Hz observed in the stellar mass black hole candidates and also suggest the existence of QPOs with the period of hours in active galactic nuclei. When the accretion rate     is low, the luminosity increases in proportion to the accretion rate. However, when     greatly exceeds the Eddington critical rate     , the luminosity is insensitive to the accretion rate and is kept constantly around  ∼3 L _E  . On the other hand, the mass-outflow rate     increases in proportion to     and it amounts to about a few per cent of the input mass-flow rate.