Many important problems in astrophysics, space physics, and geophysics involve flows of (possibly ionized) gases in the vicinity of a spherical object, such as a star or planet. The geometry of such a system naturally favors numerical schemes based on a spherical mesh. Despite its orthogonality property, the polar (latitude-longitude) mesh is ill suited for computation because of the singularity on the polar axis, leading to a highly non-uniform distribution of zone sizes. The consequences are (a) loss of accuracy due to large variations in zone aspect ratios, and (b) poor computational efficiency from a severe limitations on the time stepping. Geodesic meshes, based on a central projection using a Platonic solid as a template, solve the anisotropy problem, but increase the complexity of the resulting computer code. We describe a new finite volume implementation of Euler and MHD systems of equations on a triangular geodesic mesh (TGM) that is accurate up to fourth order in space and time and conserves the divergence of magnetic field to machine precision. The paper discusses in detail the generation of a TGM, the domain decomposition techniques, three-dimensional conservative reconstruction, and time stepping.

Magnetic fields advected along with low angular momentum accretion flows predominantly become toroidal due to the strong azimuthal velocity close to a black hole. We study self-consistently the movements of these flux tubes inside an advective disc and how they dynamically influence the flow. We find that the centrifugal barrier slows down the radial motion of the flux tubes. In this case, the large magnetic flux tubes with a significant drag force escape along the vertical axis due to buoyancy. Magnetic pressure rises close to the black hole and together with the centrifugal force, it combats gravity. The tug-of-war among these forces causes the centrifugal pressure supported shock to oscillate radially. We study the effects of successive injection of flux tubes and find how the flux tube could be trapped inside the disc in regions of highest entropy. Most interestingly, the shock wave remains at its average location and is not destroyed. We show that the toroidal field loops contribute significantly to collimate and accelerate the outflows from the centrifugal barrier and suggest this mechanism to be a way to collimate and accelerate jets.

There is a great need in several areas of astrophysics and space-physics to carry out high order of accuracy, divergence-free MHD simulations on spherical meshes. This requires us to pay careful attention to the interplay between mesh quality and numerical algorithms. Methods have been designed that fundamentally integrate high order isoparametric mappings with the other high accuracy algorithms that are needed for divergence-free MHD simulations on geodesic meshes. The goal of this paper is to document such algorithms that are implemented in the geodesic mesh version of the RIEMANN code. The fluid variables are reconstructed using a special kind of WENO-AO algorithm that integrates the mesh geometry into the reconstruction process from the ground-up. A novel divergence-free reconstruction strategy for the magnetic field that performs efficiently at all orders, even on isoparametrically mapped meshes, is then presented. The MHD equations are evolved in space and time using a novel ADER predictor algorithm that is efficiently adapted to the isoparametrically mapped geometry. The application of one-dimensional and multidimensional Riemann solvers at suitable locations on the mesh then provides the corrector step. The corrector step for the magnetic field uses a Yee-type staggering of magnetic fields. This results in a scheme with divergence-free update for the magnetic field. The use of ADER enables a one-step update which only requires one messaging operation per complete timestep. This is very beneficial for parallel processing. Several accuracy tests are presented as are stringent test problems. PetaScale performance is also demonstrated on the largest available supercomputers.

The small-scale turbulent dynamo is an important process contributing to the cosmic magnetization. In partially ionized astrophysical plasmas, the dynamo growth of magnetic energy strongly depends on the coupling state between ions and neutrals and the ion-neutral collisional damping effect. A new damping stage of turbulent dynamo in a weakly ionized medium was theoretically predicted by Xu \& Lazarian (2016). By carrying out a 3D two-fluid dynamo simulation, here we for the first time numerically confirmed the physical conditions and the linear-in-time growth of magnetic field strength of the damping stage of dynamo. The dynamo-amplified magnetic field has a characteristic length as the damping scale, which increases with time and can reach the injection scale of turbulence after around eight largest eddy-turnover times given sufficiently low ionization fraction and weak initial magnetic field. Due to the weak coupling between ions and neutrals, most turbulent energy carried by neutrals cannot be converted to the magnetic energy, resulting in a relatively weak magnetic field at the end of dynamo. This result has important implications for the growth of magnetic fields in the partially ionized interstellar medium and shock acceleration of Galactic cosmic rays.

We study time evolution of sub-Keplerian transonic accretion flows onto black holes using a general relativistic numerical simulation code. We perform simulations around the black holes having non-zero rotation. We first compare one-dimensional simulation results with theoretical results and validate the performance of our code. Next, we present results of axisymmetric, two-dimensional simulation of advective flows. In the literature, there is no solution which describes steady shock solutions in two dimensions. However, our simulations produce these centrifugal force supported steady shock waves even in presence of strong dragging of inertial frames. Since the post-shock region could be hot and upscatter photons through Comptonization, these shock would put imprints on the spectra. Thus, our solutions, which represent truly new results, could be useful to measure spins through radiation spectrum of accreting Kerr black holes.

Two Component Advective Flow (TCAF) successfully explains the spectral and tem- poral properties of outbursting or persistent sources. Images of static TCAF with Compton cloud or CENtrifugal pressure supported Boundary Layer (CENBOL) due to gravitational bending of photons have been studied before. In this paper, we study time dependent images of advective flows around a Schwarzschild black hole which include cooling effects due to Comptonization of soft photons from a Keplerian disks well as the self-consistently produced jets and outflows. We show the overall image of the disk-jet system after convolving with a typical beamwidth. A long exposure image with time dependent system need not show the black hole horizon conspicuously, un- less one is looking at a soft state with no jet or the system along the jet axis. Assuming these disk-jet configurations are relevant to radio emitting systems also, our results would be useful to look for event horizons in high accretion rate Supermassive Black Holes in Seyfert galaxies, RL Quasars.

We study time evolution of sub-Keplerian transonic accretion flows onto black holes using a general relativistic numerical simulation code. We perform simulations in Schwarzschild spacetime. We first compare one-dimensional simulation results with theoretical results and validate the performance of our code. Next, we present results of axisymmetric, two-dimensional simulation of advective flows. We find that even in this case, for which no complete theoretical analysis is present in the literature, steady state shock formation is possible.

In the magnetosphere of a rotating black hole, an inner Alfven critical surface (IACS) must be crossed by inflowing plasma. Inside the IACS, Alfven waves are inward directed toward the black hole. The majority of the proper volume of the active region of spacetime (the ergosphere) is inside of the IACS. The charge and the totally transverse momentum flux (the momentum flux transverse to both the wave normal and the unperturbed magnetic field) are both determined exclusively by the Alfven polarization. Thus, it is important for numerical simulations of black hole magnetospheres to minimize the dissipation of Alfven waves. Elements of the dissipated wave emerge in adjacent cells regardless of the IACS, there is no mechanism to prevent Alfvenic information from crossing outward. Thus, numerical dissipation can affect how simulated magnetospheres attain the substantial Goldreich-Julian charge density associated with the rotating magnetic field. In order to help minimize dissipation of Alfven waves in relativistic numerical simulations we have formulated a one-dimensional Riemann solver, called HLLI, which incorporates the Alfven discontinuity and the contact discontinuity. We have also formulated a multidimensional Riemann solver, called MuSIC, that enables low dissipation propagation of Alfven waves in multiple dimensions. The importance of higher order schemes in lowering the numerical dissipation of Alfven waves is also catalogued.

In various astrophysics settings it is common to have a two-fluid relativistic plasma that interacts with the electromagnetic field. While it is common to ignore the displacement current in the ideal, classical magnetohydrodynamic limit, when the flows become relativistic this approximation is less than absolutely well-justified. In such a situation, it is more natural to consider a positively charged fluid made up of positrons or protons interacting with a negatively charged fluid made up of electrons. The two fluids interact collectively with the full set of Maxwell's equations. As a result, a solution strategy for that coupled system of equations is sought and found here. Our strategy extends to higher orders, providing increasing accuracy. Three important innovations are reported here. In our first innovation, the magnetic field within each zone is reconstructed in a divergence-free fashion while the electric field within each zone is reconstructed in a form that is consistent with Gauss' law. In our second innovation, a multidimensionally upwinded strategy is presented which ensures that the magnetic field can be updated via a discrete interpretation of Faraday's law and the electric field can be updated via a discrete interpretation of the generalized Ampere's law. Our third innovation consists of an efficient design for several popular IMEX schemes so that they provide strong coupling between the finite-volume-based fluid solver and the electromagnetic fields at high order. Several accuracy analyses are presented showing that our method meets its design accuracy in the MHD limit as well as in the limit of electromagnetic wave propagation. Several stringent test problems are also presented. We also present a relativistic version of the GEM problem, which shows that our algorithm can successfully adapt to challenging problems in high energy astrophysics.

A black hole accretion is necessarily transonic. In presence of sufficiently high viscosity and cooling effects, a low-angular momentum transonic flow can become a standard Keplerian disc except close to the where hole where it must pass through the inner sonic point. However, if the viscosity is not high everywhere and cooling is not efficient everywhere, the flow cannot completely become a Keplerian disc. In this paper, we show results of rigorous numerical simulations of a transonic flow having vertically varying viscosity parameter (being highest on the equatorial plane) and optical depth dependent cooling processes to show that the flow indeed segregates into two distinct components as it approaches a black hole. The component on the equatorial plane has properties of a standard Keplerian disc, though the flow is not truncated at the inner- most stable circular orbit. This component extends till the horizon as a sub-Keplerian flow. This standard disc is found to be surrounded by a hot, low angular momentum component forming a centrifugal barrier dominated oscillating shock wave, consistent with the Chakrabarti-Titarchuk two component advective flow configuration.

A black hole accretion may have both the Keplerian and the sub-Keplerian components. We consider the most general accretion flow configuration, namely, two-component advective flow (TCAF) in which the Keplerian disk is immersed inside a low angular momentum, accreting sub-Keplerian halo component around a black hole. Low energy (soft) photons from the Keplerian component and hot electrons in the sub-Keplerian component exchange their energy through Comptonization or inverse-Comptonization processes. In the sub-Keplerian component, a shock is generally formed due to the centrifugal force. The post-shock region is known as the CENtrifugal pressure dominated BOundary Layer (CENBOL). The spectral and the timing properties of TCAF have been extensively studied using mostly analytical and some time dependent numerical simulations since the model was proposed by Chakrabarti & Titarchuk in 1995. The findings are the key inputs of understanding several observed features of black hole candidates. In this thesis, using numerical simulation, we rigorously prove some of the conjectures of the TCAF model. In the work presented in this thesis, we have considered for the first time the presence of both the Keplerian and the sub-Keplerian flow in a single simulation. The Keplerian disk resides on the equatorial plane and is the standard disk from which low energy photons having multi-color blackbody spectrum is emitted. The hydrodynamics as well as the thermal properties of the sub-Keplerian halo are simulated using a finite difference code which uses the principle of total variation diminishing (TVD). The Comptonization between the photons and the hot electrons is simulated using a Monte Carlo code. These two codes are then coupled and the resulting localized heating and cooling are included in the coupled code. Using this code, we study the spectral and timing properties of the TCAF.

We study the spectral and timing properties of a two component advective flow (TCAF) around a black hole by numerical simulation. Several cases have been simulated by varying the Keplerian disk rate and the resulting spectra and lightcurves have been produced for all the cases. The dependence of the spectral states and quasi-periodic oscillation (QPO) frequencies on the flow parameters is discussed. We also find the earlier explanation of arising of QPOs as the resonance between infall time scale and cooling time scale remain valid even for Compton cooling.

Low and intermediate frequency quasi-periodic oscillations (QPOs) in black hole candidates are believed to be due to oscillations of the Comptonizing regions in an accretion flow. Assuming that the general structure of an accretion disk is a Two Component Advective Flow (TCAF), we numerically simulate the light curves emitted from an accretion disk for different accretion rates and find how the QPO frequencies vary. We use a standard Keplerian disk residing at the equatorial plane as a source of soft photons. These soft photons, after suffering multiple scattering with the hot electrons of the low angular momentum, sub-Keplerian, flow emerge out as hard radiation. The hydrodynamic and thermal properties of the electron cloud is simulated using a Total Variation Diminishing (TVD) code. The TVD code is then coupled with a radiative transfer code which simulates the energy exchange between the electron and radiation using Monte Carlo technique. The resulting localized heating and cooling are included also. We find that the QPO frequency increases and the spectrum becomes softer as we increase the Keplerian disk rate. However, the spectrum becomes harder if we increase the sub-Keplerian accretion rate. We find that an earlier prediction that QPOs occur when the infall time scale roughly matches with the cooling time scale, originally obtained using a power-law cooling, remains valid even for Compton cooling. Our findings agree with the general observations of low frequency QPOs in black hole candidates.

It has been recently shown that transonic electron positron fluid is the least relativistic, compared to the fluid containing finite proportion of baryons. We compute spectra from these flows in general relativity (GR) including the effect of light bending. We consider the bremsstrahlung process to supply the seed photons. We choose accretion in the advective domain, and for simplicity the radial accretion or Bondi type accretion. We show that electron positron accreting flow produces the softest spectra and the lowest luminosity.

The thermal state of the fluid is governed by the ratio of the thermal and the rest energy. This brings the composition of the fluid into the picture. Although, fluid composed of lighter particles (e.g: electron-positron pair plasma) at same temperature, is more relativistic compared to fluids with finite baryon loading, but this is not necessarily true when baryon poor transonic fluid are compared with each other. It can be shown that the transonic pair-fluid is the least relativistic. This has far reaching consequences on accreting flows around compact objects and are expected to have similar effect on relativistic outflows and explosive events as well.

The effect of cooling on the outflow rate from an accretion disk around a black hole is investigated using a coupled Monte Carlo Total Variation Diminishing code. A correlation between the spectral states and the outflow rates is found as a consequence.

We wish to investigate the effects of cooling of the Compton cloud on the outflow formation rate in an accretion disk around a black hole. We carry out a time dependent numerical simulation where both the hydrodynamics and the radiative transfer processes are coupled together. We consider a two-component accretion flow in which the Keplerian disk is immersed into an accreting low-angular momentum flow (halo) around a black hole. The soft photons which originate from the Keplerian disk are inverse-Comptonized by the electrons in the halo and the region between the centrifugal pressure supported shocks and the horizon. We run several cases by changing the rate of the Keplerian disk and see the effects on the shock location and properties of the outflow and the spectrum. We show that as a result of Comptonization of the Compton cloud, the cloud becomes cooler with the increase in the Keplerian disk rate. As the resultant thermal pressure is reduced, the post-shock region collapses and the outflow rate is also reduced. Since the hard radiation is produced from the post-shock region, and the spectral slope increases with the reduction of the electron temperature, the cooling produces softer spectrum. We thus find a direct correlation between the spectral states and the outflow rates of an accreting black hole.

We carry out a time dependent numerical simulation where both the hydrodynamics and the radiative transfer are coupled together. We consider a two-component accretion flow in which the Keplerian disk is immersed inside an accreting low angular momentum flow (halo) around a black hole. The injected soft photons from the Keplerian disk are reprocessed by the electrons in the halo. We show that in presence of an axisymmetric soft-photon source, the spherically symmetric Bondi flow losses its symmetry and becomes axisymmetric. The low angular momentum flow was observed to slow down close to the axis and formed a centrifugal barrier which added new features into the spectrum. Using the Monte Carlo method, we generated the radiated spectra as functions of the accretion rates. We find that the transitions from a hard state to a soft state is determined by the mass accretion rates of the disk and the halo. We separate out the signature of the bulk motion Comptonization and discuss its significance. We study how the net spectrum is contributed by photons suffering different number of scatterings and spending different amounts of time inside the Compton cloud. We study the directional dependence of the emitted spectrum as well.

A black hole accretion may have both the Keplerian and the sub-Keplerian component. In the so-called Chakrabarti-Titarchuk scenario, the Keplerian component supplies low energy (soft) photons while the sub-Keplerian component supplies hot electrons which exchange their energy with the soft photons through Comptonization or inverse Comptonization processes. In the sub-Keplerian component, a shock is generally produced due to the centrifugal force. The postshock region is known as the CENtrifugal pressure-supported BOundary Layer (CENBOL). In this paper, we compute the effects of the thermal and the bulk motion Comptonization on the soft photons emitted from a Keplerian disk by the CENBOL, the preshock sub-Keplerian disk and the outflowing jet. We study the emerging spectrum when the converging inflow and the diverging outflow (generated from the CENBOL) are simultaneously present. From the strength of the shock, we calculate the percentage of matter being carried away by the outflow and determine how the emerging spectrum depends on the outflow rate. The preshock sub-Keplerian flow is also found to Comptonize the soft photons significantly. The interplay between the up-scattering and down-scattering effects determines the effective shape of the emerging spectrum. By simulating several cases with various inflow parameters, we conclude that whether the preshock flow, or the postshock CENBOL or the emerging jet is dominant in shaping the emerging spectrum depends strongly on the geometry of the flow and the strength of the shock in the sub-Keplerian flow.