有交界面的椭圆型方程论文：1种求解有交界面的椭圆型方程的间断条件捕捉格式.docVIP

有交界面的椭圆型方程论文：一种求解有交界面的椭圆型方程的间断条件捕捉格式 【中文摘要】浸入界面方法主要用来求解带有界面问题的偏微分方程,目前已经被广泛的运用到计算流体力学领域中。界面问题所导出的偏微分方程的解在穿过界面时是间断的,许多偏微分方程常规的数值计算方法都不能够很好的运用到求解界面问题上。本文用浸入界面方法求解了二维的有交界面的椭圆型方程,并且将这种方法推广到界面及间断条件都隐式地定义在界面附近的网格函数点上,给出了一种间断条件捕捉格式。它特别适合于隐式界面跟踪法如水平集方法。针对这种新的情形,对于原浸入界面方法中的界面关系、不规则点差分格式中的系数代数方程组和修正项Cij都进行了相应的修正。紧邻界面的网格点处差分格式的截断误差是一阶,而界面是低一维的,所以解的整体误差可以达到二阶。该格式利用标准的三阶拉格郎日插值计算间断函数沿界面的导数,易于执行。数值计算验证了该法关于最大模的二阶收敛性。 【英文摘要】The immersed interface method is general technique for solving differential equations with interface problems and now it is applied in various aspects of fluid dynamics. Interface problems usually lead to differential equations whose solutions and its derivatives have discontinuities or non-smoothness across some interface. Many differential equation numerical methods do not work efficiently for interface problems.In this paper the immersed interface method for solving elliptic interface prob-lems in 2D is extended to the situation when interface and jumps are implicitly provided by functions on the grid points in a neighborhood of the interface. The resulting scheme is particularly suitable for interface capturing methods such as the level set method. Formulations of the original immersed interface method are mod-ified for this implicit setting, including the interface jump relations, the algebraic equations for determining the coefficients of the finite difference scheme at irregular grid points and the correction terms. The local truncation error at grid points just adjacent to the interface is of first order. Since the interface is of co-dimension one, the global second order accuracy for the solution in maximum norm can be achieved, as demonstrated by numerical examples. It is easier to implement since standard Lagrange interpolation is used to compute the derivatives of the interface quantities. 【关键词】有交界面的椭圆型方程 间断条件捕捉方法 浸入界面方法 水平集方法 差分方法 【英文关键词】elliptic interface problem jump condition capturing scheme im-mersed interface method level set