To obtain a simple description of a geometrically thin magnetic accretion disk, we apply the method of asymptotic expansion. For the first time we write a full set of stationary asymptotic approximation equations of a thin magnetic accretion disk. As the obtained equations cannot be solved without knowledge of the solutions at the disk surface, we combine the results from numerical simulations and from analytical equations to find a simple set of functional expressions describing the radial and vertical dependence of physical quantities in the disk. Except very close to the star, the functional form of the disk variables is quite similar in the HD and MHD cases, with the overall scale of density, and the vertical and radial velocity components modified by the stellar magnetic field.

We obtained equations for a thin magnetic accretion disk, using the method of asymptotic approximation. They cannot be solved analytically-without solutions for a magnetic field in the magnetosphere between the star and the disk, only a set of general conditions on the solutions can be derived. To compare the analytical results with numerical solutions, we find expressions for physical quantities in the disk, using our results from resistive and viscous star-disk magnetospheric interaction simulations.

We perform MHD simulations of a thin resistive and viscous accretion disk around a neutron star with the surface dipolar magnetic field of 10$^8$ Gauss. The system evolution is followed during the interval of 500 millisecond pulsar rotations. Matter is accreted through a stable accretion column from the disk onto the star. We also show propagation of the stellar wind through the corona. Analysis of the mass accretion flux and torques on the star shows that the disk reaches the quasi-stationary state.

We performed axisymmetric, grid-based, ideal magnetohydrodynamic (MHD) simulations of oscillating cusp-filling tori orbiting a non-rotating neutron star. A pseudo-Newtonian potential was used to construct the constant angular momentum tori in equilibrium. The inner edge of the torus is terminated by a "cusp" in the effective potential. The initial motion of the model tori was perturbed with uniform sub-sonic vertical and diagonal velocity fields. As the configuration evolved in time, we measured the mass accretion rate on the neutron star surface and obtained the power spectrum. The prominent mode of oscillation in the cusp torus is the radial epicyclic mode. It would appear that vertical oscillations are suppressed by the presence of the cusp. From our analysis it follows that the mass accretion rate carries a modulation imprint of the oscillating torus, and hence so does the boundary layer luminosity.

We have modeled hydrodynamical tori in the Klu{\'z}niak-Lee pseudo-Newtonian potential. The tori in equilibrium were perturbed with uniform sub-sonic velocity fields: vertical, radial and diagonal respectively, and allowed to evolve in time. We identify the eigenmodes corresponding to those of slender tori. The results of our simulations are relevant to the investigation of high-frequency quasi-periodic oscillations observed in stellar-mass black hole binaries.

We performed axisymmetric hydrodynamical simulations of oscillating tori orbiting a non-rotating black hole. The tori in equilibrium were constructed with a constant distribution of angular momentum in a pseudo-Newtonian potential (Klu{\'z}niak-Lee). Motions of the torus were triggered by adding sub-sonic velocity fields: radial, vertical and diagonal to the tori in equilibrium. As the perturbed tori evolved in time, we measured $L_{2}$ norm of density and obtained the power spectrum of $L_{2}$ norm which manifested eigenfrequencies of tori modes. The most prominent modes of oscillation excited in the torus by a quasi-random perturbation are the breathing mode and the radial and vertical epicyclic modes. The radial and the plus modes, as well as the vertical and the breathing modes will have frequencies in an approximate 3:2 ratio if the torus is several Schwarzschild radii away from the innermost stable circular orbit. Results of our simulations may be of interest in the context of high-frequency quasi-periodic oscillations (HF QPOs) observed in stellar-mass black hole binaries, as well as in supermassive black holes.

Our goal is to perform global simulations of thin accretion discs around compact bodies like neutron stars with dipolar magnetic profile and black holes by exploiting the facilities provided by state-of-the-art grid-based, high resolution shock capturing (HRSC) and finite volume codes. We have used the Godunov-type code PLUTO to simulate a thin disc around a compact object prescribed with a pseudo-Newtonian potential in a purely hydrodynamical (HD) regime, with numerical viscosity as a first step towards achieving our goal as mentioned above.

Given a correlation function (or n-point function), can the corresponding nature of space-time be determined ? To answer this question it is required to derive the Ward Identity (WI), analyse the symmetries and arrive at the law of conservation. Modus operandi involves Lie differentiating two-point function considering the symmetry to be non-anomalous. The WI so obtained is shown to form a Lie algebra which determines the nature of space-time. Solving the identity results in a law of conservation, which physically explains the reason for WI to form an algebra and contains in it an equation of motion for four-point function. As a special case, a relation between mass and potential involving the spatial derivatives of four- and five- point function is obtained. Finally, the conservation equation is exploited to get the probability amplitude for the two-point function which shows how correlation functions provide an opportunity to probe the fundamental laws of physics.