We propose an FFT-based algorithm for computing fundamental solutions of difference operators with constant coefficients. Our main contribution is to handle cases where the symbol has zeros.
The iterative solution of systems of equations arising from partial differential equations (PDEs) governing boundary layer flow for large Reynolds numbers is studied. The PDEs are discretized using finite volume or finite difference approximations on tensor product grids. We consider a convergence acceleration technique, where a semicirculant approximation of the spatial difference operator is employed as preconditioner. A relevant model problem is derived, and the spectrum of the preconditioned coefficient matrix is analyzed. It is proved that, asymptotically, the time step for the forward Euler method could be chosen as a constant, which is independent of the number of gridpoints and the Reynolds number Re. The same type of result is also derived for finite size grids, where the solution fulfills a given accuracy requirement. By linearizing the Navier-Stokes equations around an approximate solution, we form a system of linear PDEs with variable coefficients. When utilizing the semicirculant preconditioner for this problem, which has properties very similar to the full nonlinear equations, the results show that the favorable convergence properties hold also here. We compare the semicirculant preconditioner to a multigrid scheme. The number of iterations and the arithmetic complexities are considered, and it clear that semicirculant method is much more efficient for problems where the Reynolds number is large. The number of iterations for the multigrid method grows like √ Re, while the convergence rate for the scheme using semicirculant preconditioning is independent of Re. Also, the multigrid scheme is very sensitive to the level of artificial dissipation, while the method using semicirculant preconditioning is not. Introduction.The following strategy for solving the Navier-Stokes equations governing compressible, stationary flow is frequently employed in computer codes used by the computational fluid dynamics (CFD) community:(i) Generate a tensor product grid or a multiblock grid where the partition grids are tensor product grids.(ii) Discretize the equations in space utilizing a finite volume or finite difference scheme.(iii) Integrate the corresponding time-dependent problem in time by using an explicit Runge-Kutta time-marching method.(iv) Improve the convergence properties by employing some convergence acceleration technique.In section 2 of this paper, we review the semicirculant convergence acceleration, or preconditioning, technique, which is applicable to the type of computations described above. Then in section 3, we describe the boundary layer equations, which yield an approximate solution to the compressible Navier-Stokes equations for the 2D flow over a flat plate. We employ the boundary layer equations for formulating a relevant scalar model problem, and for linearizing the Navier-Stokes equations, resulting in a
Abstract. We present a new preconditioner for the iterative solution of linear systems of equations arising from discretizations of systems of first order partial differential equations (PDEs) on structured grids. Such systems occur in many important applications, including compressible fluid flow and electormagnetic wave propagation. The preconditioner is a truncated convolution operator, with a kernel that is a fundamental solution of a difference operator closely related to the original discretization.Analysis of a relevant scalar model problem in two spatial dimensions shows that grid independent convergence is obtained using a simple one-stage iterative method. As an example of a more involved problem, we consider the steady state solution of the non-linear Euler equations in a two dimensional, non-axisymmetric duct. We present results from numerical experiments, verifying that the preconditioning technique again achieves grid independent convergence, both for an upwind discretization and for a centered second order discretization with fourth order artificial viscosity.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
customersupport@researchsolutions.com
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.