A finite-element coarse-grid projection method for incompressible flow simulations

被引:13
|
作者
Kashefi, Ali [1 ,2 ]
Staples, Anne E. [2 ]
机构
[1] Stanford Univ, Inst Computat & Math Engn, Dept Geophys, Stanford, CA 94305 USA
[2] Virginia Tech, Engn Sci & Mech Program, Dept Biomed Engn & Mech, Blacksburg, VA 24061 USA
关键词
Pressure-correction schemes; Coarse-grid projection; Unstructured grids; Finite elements; Semi-implicit time integration; NAVIER-STOKES EQUATIONS; BACKWARD-FACING STEP; MULTIGRID METHOD; CIRCULAR-CYLINDER; LAMINAR-FLOW; SCHEME; DISCRETIZATIONS; APPROXIMATION; ALGORITHM; PRESSURE;
D O I
10.1007/s10444-017-9573-5
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
Coarse grid projection (CGP) methodology is a novel multigrid method for systems involving decoupled nonlinear evolution equations and linear elliptic Poisson equations. The nonlinear equations are solved on a fine grid and the linear equations are solved on a corresponding coarsened grid. Mapping operators execute data transfer between the grids. The CGP framework is constructed upon spatial and temporal discretization schemes. This framework has been established for finite volume/difference discretizations as well as explicit time integration methods. In this article we present for the first time a version of CGP for finite element discretizations, which uses a semi-implicit time integration scheme. The mapping functions correspond to the finite-element shape functions. With the novel data structure introduced, the mapping computational cost becomes insignificant. We apply CGP to pressure-correction schemes used for the incompressible Navier-Stokes flow computations. This version is validated on standard test cases with realistic boundary conditions using unstructured triangular meshes. We also pioneer investigations of the effects of CGP on the accuracy of the pressure field. It is found that although CGP reduces the pressure field accuracy, it preserves the accuracy of the pressure gradient and thus the velocity field, while achieving speedup factors ranging from approximately 2 to 30. The minimum speedup occurs for velocity Dirichlet boundary conditions, while the maximum speedup occurs for open boundary conditions.
引用
收藏
页码:1063 / 1090
页数:28
相关论文
共 50 条