A method for enhancing the stability and robustness of explicit schemes in astrophysical fluid dynamics

A method for enhancing the stability and robustness of explicit schemes in computational fluid dynamics is presented. The method is based in reformulating explicit schemes in matrix form, which cane modified gradually into semi or strongly-implicit schemes. From the point of view of matrix-algebra, explicit numerical methods are special cases in which the global matrix of coefficients is reduced to the identity matrix $I$. This extreme simplification leads to severer stability range, hence of their robustness. In this paper it is shown that a condition, which is similar to the Courant-Friedrich-Levy (CFL) condition can be obtained from the stability requirement of inversion of the coefficient matrix. This condition is shown to be relax-able, and that a class of methods that range from explicit to strongly implicit methods can be constructed, whose degree of implicitness depends on the number of coefficients used in constructing the corresponding coefficient-matrices. Special attention is given to a simple and tractable semi-explicit method, which is obtained by modifying the coefficient matrix from the identity matrix $I$ into a diagonal-matrix $D$. This method is shown to be stable, robust and it can be applied to search for stationary solutions using large CFL-numbers, though it converges slower than its implicit counterpart. Moreover, the method can be applied to follow the evolution of strongly time-dependent flows, though it is not as efficient as normal explicit methods. In addition, we find that the residual smoothing method accelerates convergene toward steady state solutions considerably and improves the efficiency of the solution procedure.Comment: 33 pages, 15 figure

Problem-orientable numerical algorithm for modelling multi-dimensional radiative MHD flows in astrophysics -- the hierarchical solution scenario

We present a hierarchical approach for enhancing the robustness of numerical solvers for modelling radiative MHD flows in multi-dimensions. This approach is based on clustering the entries of the global Jacobian in a hierarchical manner that enables employing a variety of solution procedures ranging from a purely explicit time-stepping up to fully implicit schemes. A gradual coupling of the radiative MHD equation with the radiative transfer equation in higher dimensions is possible. Using this approach, it is possible to follow the evolution of strongly time-dependent flows with low/high accuracies and with efficiency comparable to explicit methods, as well as searching quasi-stationary solutions for highly viscous flows. In particular, it is shown that the hierarchical approach is capable of modelling the formation of jets in active galactic nuclei and reproduce the corresponding spectral energy distribution with a reasonable accuracy.Comment: 28 pages, 9 figure

A Meshfree Generalized Finite Difference Method for Surface PDEs

In this paper, we propose a novel meshfree Generalized Finite Difference Method (GFDM) approach to discretize PDEs defined on manifolds. Derivative approximations for the same are done directly on the tangent space, in a manner that mimics the procedure followed in volume-based meshfree GFDMs. As a result, the proposed method not only does not require a mesh, it also does not require an explicit reconstruction of the manifold. In contrast to existing methods, it avoids the complexities of dealing with a manifold metric, while also avoiding the need to solve a PDE in the embedding space. A major advantage of this method is that all developments in usual volume-based numerical methods can be directly ported over to surfaces using this framework. We propose discretizations of the surface gradient operator, the surface Laplacian and surface Diffusion operators. Possibilities to deal with anisotropic and discontinous surface properties (with large jumps) are also introduced, and a few practical applications are presented

Iterative spectral methods and spectral solutions to compressible flows

A spectral multigrid scheme is described which can solve pseudospectral discretizations of self-adjoint elliptic problems in O(N log N) operations. An iterative technique for efficiently implementing semi-implicit time-stepping for pseudospectral discretizations of Navier-Stokes equations is discussed. This approach can handle variable coefficient terms in an effective manner. Pseudospectral solutions of compressible flow problems are presented. These include one dimensional problems and two dimensional Euler solutions. Results are given both for shock-capturing approaches and for shock-fitting ones

Approximate tensor-product preconditioners for very high order discontinuous Galerkin methods

In this paper, we develop a new tensor-product based preconditioner for discontinuous Galerkin methods with polynomial degrees higher than those typically employed. This preconditioner uses an automatic, purely algebraic method to approximate the exact block Jacobi preconditioner by Kronecker products of several small, one-dimensional matrices. Traditional matrix-based preconditioners require $\mathcal{O}(p^{2d})$ storage and $\mathcal{O}(p^{3d})$ computational work, where $p$ is the degree of basis polynomials used, and $d$ is the spatial dimension. Our SVD-based tensor-product preconditioner requires $\mathcal{O}(p^{d+1})$ storage, $\mathcal{O}(p^{d+1})$ work in two spatial dimensions, and $\mathcal{O}(p^{d+2})$ work in three spatial dimensions. Combined with a matrix-free Newton-Krylov solver, these preconditioners allow for the solution of DG systems in linear time in $p$ per degree of freedom in 2D, and reduce the computational complexity from $\mathcal{O}(p^9)$ to $\mathcal{O}(p^5)$ in 3D. Numerical results are shown in 2D and 3D for the advection and Euler equations, using polynomials of degree up to $p=15$. For many test cases, the preconditioner results in similar iteration counts when compared with the exact block Jacobi preconditioner, and performance is significantly improved for high polynomial degrees $p$.Comment: 40 pages, 15 figure

An implicit numerical algorithm for solving the general relativistic hydrodynamical equations around accreting compact objects

An implicit algorithm for solving the equations of general relativistic hydrodynamics in conservative form in three-dimensional axi-symmetry is presented. This algorithm is a direct extension of the pseudo-Newtonian implicit radiative magnetohydrodynamical solver -IRMHD- into the general relativistic regime. We adopt the Boyer-Lindquist coordinates and formulate the hydrodynamical equations in the fixed background of a Kerr black hole. The set of equations are solved implicitly using the hierarchical solution scenario (HSS). The HSS is efficient, robust and enables the use of a variety of solution procedures that range from a purely explicit up to fully implicit schemes. The discretization of the HD-equations is based on the finite volume formulation and the defect-correction iteration strategy for recovering higher order spatial and temporal accuracies. Depending on the astrophysical problem, a variety of relaxation methods can be applied. In particular the vectorized black-white Line-Gauss-Seidel relaxation method is most suitable for modeling accretion flows with shocks, whereas the Approximate Factorization Method is for weakly compressible flows. The results of several test calculations that verify the accuracy and robustness of the algorithm are shown. Extending the algorithm to enable solving the non-ideal MHD equations in the general relativistic regime is the subject of our ongoing research.Comment: 30 pages, 8 figures, to be published in New Astronom

Spectral/hp element methods: recent developments, applications, and perspectives

The spectral/hp element method combines the geometric flexibility of the classical h-type finite element technique with the desirable numerical properties of spectral methods, employing high-degree piecewise polynomial basis functions on coarse finite element-type meshes. The spatial approximation is based upon orthogonal polynomials, such as Legendre or Chebychev polynomials, modified to accommodate C0-continuous expansions. Computationally and theoretically, by increasing the polynomial order p, high-precision solutions and fast convergence can be obtained and, in particular, under certain regularity assumptions an exponential reduction in approximation error between numerical and exact solutions can be achieved. This method has now been applied in many simulation studies of both fundamental and practical engineering flows. This paper briefly describes the formulation of the spectral/hp element method and provides an overview of its application to computational fluid dynamics. In particular, it focuses on the use the spectral/hp element method in transitional flows and ocean engineering. Finally, some of the major challenges to be overcome in order to use the spectral/hp element method in more complex science and engineering applications are discussed

Spectral multigrid methods for the solution of homogeneous turbulence problems

New three-dimensional spectral multigrid algorithms are analyzed and implemented to solve the variable coefficient Helmholtz equation. Periodicity is assumed in all three directions which leads to a Fourier collocation representation. Convergence rates are theoretically predicted and confirmed through numerical tests. Residual averaging results in a spectral radius of 0.2 for the variable coefficient Poisson equation. In general, non-stationary Richardson must be used for the Helmholtz equation. The algorithms developed are applied to the large-eddy simulation of incompressible isotropic turbulence