一种基于过滤与反卷积的新型高阶浸没边界法-应用数学和力学.PDF

摇 , 34 8 Applied Mathematics and Mechanics摇 摇 2013年8月15 日出版 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 Vol.34,No.8,Aug.15,2013 文章编号:1000鄄0887(2013)08鄄0807鄄08 訫 应用数学和力学编委会,ISSN 1000鄄0887 一一种基于过滤与反卷积的 新型高阶浸没边界法* 方摇 乐1,2,3 (1. 北京航空航天大学 中法工程师学院,北京 100191; 2. 中国科学院非线性力学国家重点实验室,北京 100090; 3. 西南交通大学 高速铁路线路工程教育部重点实验室,成都610031) 摘要:摇 传统浸没边界法在边界附近只有一阶精度,而高精度的改进方法都需要额外引入跳跃条件, 因此不具备普适性.文中设计了一种基于过滤和反卷积的新型算法,既在一定程度上提高了精度,又 避免了以往方法中引入额外跳跃条件的难题.通过一个简单的一维算例验证了新算法可以达到接近 二阶精度,其具体的精度值与反卷积步骤中选取的逆核函数在积分域边界的连续性有关. 关摇 键摇 词:摇 浸没边界法;摇 反卷积;摇 过滤 中图分类号:摇 O242.1;O242.2;O35摇 摇 摇 文献标志码:摇 A DOI:10.3879/ j.issn.1000鄄0887.2013.08.004 引摇 摇 言 在计算流体力学领域中,浸没边界法(immersedboundary method,IBM)是一种近年来得到 [1鄄3] 很快发展的新型数值方法.浸没边界法于上世纪70 年代由Peskin 在计算血液流动和流固 耦合问题中创立,它允许流体边界与计算网格不一致,将无滑移边界条件通过外力项的形式添 加到求解方程中. 由于外力项一般写成啄 函数的形式,因此如何数值求解这样的含间断的微分 方程也就是随之而来的问题.传统的浸没边界法一般采用所谓的“光滑 啄 函数冶(smoothed del鄄 ta function)方法处理啄 函数,然而通常情况下采用此方法的精度很低,以往的工作已经证明了 [4鄄5] 在间断附近用光滑啄 函数只能达到一阶精度 .为改进此精度问题,在过去的十年内有许多 [6] 学者做过重要工作,比如Lai和Peksin 提出了一种形式上二阶的离散方法,但在实际应用中 [7] 发现此方法仍然只能实现一阶精度 .也有其他学者提出了IIM(immersed interface method)、 MIB(matched interface and boundary)等改进方法[8鄄10] ,但这些方法都需要在间断附近额外引入 若干个跳跃条件(jump condition),而实际问题却并不一定能满足这些跳跃条件,因而这些方 法大多不具备普适性.如何在不引入跳跃条件的前提下实现高阶精度,一直是浸没边界法发展 过程中一个重要而棘手的问题. 本文提出一种新型浸没边界法,首先对含间断的微分方程进行过滤操作,将间断方程变为 连续方程;接着对解得的光滑函数采用反卷积即得到最终结果.新方法既在一定程度上提高了 * 收稿日期:摇 2013鄄03鄄29;修订日期:摇 2013鄄05鄄28 基金项目:摇 国家自然科学基金资助项目;非线性力学国家重点实验室开放基金资助项目 作者简介:摇 方乐(1983—),安徽歙县人,副教授,博士(E鄄mail:le.fang@). 807 808 方摇 摇 乐 精度,又避免了引入跳跃条件的问