Algebraic Construction of a Strongly Consistent, Permutationally Symmetric and Conservative Difference Scheme for 3D Steady Stokes Flow
暂无分享,去创建一个
By using symbolic algebraic computation, we construct a strongly-consistent second-order finite difference scheme for steady three-dimensional Stokes flow and a Cartesian solution grid. The scheme has the second order of accuracy and incorporates the pressure Poisson equation. This equation is the integrability condition for the discrete momentum and continuity equations. Our algebraic approach to the construction of difference schemes suggested by the second and the third authors combines the finite volume method, numerical integration, and difference elimination. We make use of the techniques of the differential and difference Janet/Grobner bases for performing related computations. To prove the strong consistency of the generated scheme, we use these bases to correlate the differential ideal generated by the polynomials in the Stokes equations with the difference ideal generated by the polynomials in the constructed difference scheme. As this takes place, our difference scheme is conservative and inherits permutation symmetry of the differential Stokes flow. For the obtained scheme, we compute the modified differential system and use it to analyze the scheme’s accuracy.
[1] Vladimir P. Gerdt,et al. Gröbner Bases and Generation of Difference Schemes for Partial Differential Equations , 2006, ArXiv.
[2] Rémi Abgrall,et al. Physics-compatible numerical methods , 2014, J. Comput. Phys..
[3] Vladimir P. Gerdt,et al. Algebraic construction and numerical behavior of a new s-consistent difference scheme for the 2D Navier-Stokes equations , 2017, Appl. Math. Comput..
[4] N. Anders Petersson,et al. Stability of pressure boundary conditions for Stokes and Navier-Stokes equations , 2001 .