期刊
JOURNAL OF COMPUTATIONAL PHYSICS
卷 229, 期 11, 页码 4225-4246出版社
ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2010.02.007
关键词
Explicit jump approximation; Ghost cell approach; Viscous flow; Diffusion equation; Heat equation; Jump correction; The immersed interface method; Irregular solid boundary; The polygonal area mapping method; HyPAM
资金
- Humboldt Foundation
- Germany NSF
- NY Sea Grant Research Institute
The ghost cell approaches (GCA) for handling stationary solid boundaries, regular or irregular, are first investigated theoretically and numerically for the diffusion equation with Dirichlet boundary conditions. The main conclusion of this part of investigation is that the approximation for the diffusion term has to be second-order accurate everywhere in order for the numerical solution to be rigorously second-order accurate. Violating this principle, the linear and quadratic GCAs have the following shortcomings: (1) restrictive constraints on grid size when the viscosity is small; (2) susceptibleness to instability of a time-explicit formulation for strongly transient flows; (3) convergence deterioration to zeroth- or first-order for solutions with high-frequency modes. Therefore, the widely-used linear extrapolation for enforcing no-slip boundary conditions should be avoided, even for regular solid boundaries. As a remedy, a simple method based on explicit jump approximation (EJA) is proposed. EJA hinges on the idea that a velocity-derivative jump at the boundary reduces to the value of the velocity-derivative at the fluid side because the velocity of the stationary boundary is zero. Although the time-marching linear system of EJA is not symmetric, it is strictly diagonal dominant with positive diagonal entries. Numerical results show that, over a large range of viscosity and grid sizes, EJA performs much better than GCAs in terms of stability and accuracy. Furthermore, the second-order convergence of EJA does not depend on viscosity and the spectrum of the solution, as those of GCAs do. This paper is written with enough details so that one can reproduce the numerical results. (C) 2010 Elsevier Inc. All rights reserved.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据