3,476 research outputs found

    Uniform convergence of V-cycle multigrid algorithms for two-dimensional fractional Feynman-Kac equation

    Full text link
    In this paper we derive new uniform convergence estimates for the V-cycle MGM applied to symmetric positive definite Toeplitz block tridiagonal matrices, by also discussing few connections with previous results. More concretely, the contributions of this paper are as follows: (1) It tackles the Toeplitz systems directly for the elliptic PDEs. (2) Simple (traditional) restriction operator and prolongation operator are employed in order to handle general Toeplitz systems at each level of the recursion. Such a technique is then applied to systems of algebraic equations generated by the difference scheme of the two-dimensional fractional Feynman-Kac equation, which describes the joint probability density function of non-Brownian motion. In particular, we consider the two coarsening strategies, i.e., doubling the mesh size (geometric MGM) and Galerkin approach (algebraic MGM), which lead to the distinct coarsening stiffness matrices in the general case: however, several numerical experiments show that the two algorithms produce almost the same error behaviour.Comment: 26 page

    Derivation of new staggered compact schemes with application to navier-stokes equations

    Get PDF
    A method is proposed for the derivation of new classes of staggered compact derivative and interpolation operators. The algorithm has its roots in an implicit interpolation theory consistent with compact schemes and reduces to the computation of the required staggered derivatives and interpolated quantities as a combination of nodal values and collocated compact derivatives. The new approach is cost-effective, simplifies the imposition of boundary conditions, and has improved spectral resolution properties, on equal order of accuracy, with respect to classical schemes. The method is applied to incompressible Navier-Stokes equations through the implementation into a staggered flow solver with a fractional step procedure, and tested on classical benchmarks.Postprint (published version
    • …
    corecore