浸没边界法的高效半隐式格式

浸没边界法的高效半隐式格式

论文摘要

在研究流体和弹性体结构相互作用的问题时,浸没边界法(Immersed BoundaryMethod)是一种非常重要而且常用的计算方法。在很多领域尤其是生物力学领域有着广泛的应用。但是浸没边界法本身也有一些局限,影响了其应用的范围,其中一个就是当使用显式格式时浸没边界法的数值刚性特别大,这就使得在计算时,时间步长必须取得非常小,那么如果想研究系统长时间的演化,计算量就会变得不可忍受,大大限制了浸没边界法的应用范围。为了克服这个困难,很多的半隐式和近似隐式格式被发展了出来,但是迄今为止还没有一种格式有令人满意的表现。在本文中,我们针对二维的定常Stokes流动,非定常Stokes流动以及Navier-Stokes流动,分别提出了几种高效的半隐式格式。这些半隐式格式具有非常好的数值稳定性,可以采用比显式格式大的多的时间步长,同时每一步的计算量与显式格式是相当的,这样在进行长时间的计算时,就可以大大的降低计算量。我们的半隐式格式的构造分为两步:第一步是找到一种无条件稳定的半隐式格式,这种半隐式格式稳定性非常好,但是单步计算量仍嫌太大,在实际计算中相对于显式格式并没有明显的优势。为了降低每一步的计算量,在第二步中我们对这种无条件稳定的半隐式格式进行小尺度分解(Small Scale Decomposition)。然后只对小尺度分解之后得到的首阶项使用隐式格式处理,经过一系列的简化,首阶项可以表示成卷积的形式,这样我们的半隐式格式就可以在傅立叶空间中显式求解,从而获得与显式格式相当的单步计算量。数值结果也证明了我们构造出的高效半隐式格式在大大改善数值稳定性的同时单步的计算量并没有很大的增加。在实际计算中,我们给出的高效半隐式格式是显式格式的一种非常好的替代方法。

论文目录

  • 摘要
  • Abstract
  • 第1章 引言
  • 第2章 浸没边界方法
  • 2.1 什么是浸没界面
  • 2.2 浸没界面的数学描述
  • 2.3 跳跃边界模型
  • 2.4 浸没边界方法的显式计算格式
  • h的选择'>2.5 近似δ函数δh的选择
  • 2.6 一维界面的局部弧长—切向角描述
  • 2.7 空间离散的谱方法描述
  • 第3章 定常stokes流动下浸没边界法的半隐式稳定格式
  • 3.1 数学描述
  • 3.2 绝对稳定的半隐式格式
  • 3.2.1 方法的描述
  • 3.2.2 稳定性分析
  • 3.3 小尺度分解
  • 3.4 定常Stokes流动的高效半隐式格式
  • 3.5 高阶精度的时间离散格式
  • 3.6 计算结果及分析
  • 3.6.1 模型问题及计算参数设定
  • 3.6.2 数值结果及分析
  • 第4章 非定常stokes流动下浸没边界法的半隐式稳定格式
  • 4.1 数学描述
  • 4.2 绝对稳定的半隐式格式
  • 4.2.1 方法的描述
  • 4.2.2 稳定性分析
  • 4.3 小尺度分解
  • 4.4 非定常Stokes流动的高效半隐式格式
  • 4.5 二阶时间精度离散格式
  • 4.5.1 绝对稳定的二阶时间精度半隐式格式
  • 4.5.2 二阶格式的小尺度分解
  • 4.5.3 二阶时间精度的高效半隐式格式
  • 4.6 计算结果及分析
  • 4.6.1 模型问题及计算参数设定
  • 4.6.2 一阶格式的数值结果
  • 4.6.3 二阶格式的数值结果
  • 第5章 Navier-Stokes流动下浸没边界法的半隐式稳定格式
  • 5.1 数学描述
  • 5.2 绝对稳定的半隐式格式
  • 5.2.1 稳定性分析
  • 5.3 小尺度分解
  • 5.4 Navier-Stokes流动的高效半隐式格式
  • 5.5 二阶时间精度离散格式
  • 5.5.1 二阶时间精度半隐式格式
  • 5.5.2 二阶时间精度的高效半隐式格式
  • 5.6 计算结果及分析
  • 5.6.1 模型问题及计算参数设定
  • 5.6.2 一阶格式的数值结果
  • 5.6.3 二阶格式的数值结果
  • 第6章 数值方法的应用:粘性涡面问题的模拟
  • 6.1 数学描述
  • 6.2 数值格式
  • 6.3 计算结果及分析
  • 第7章 结论
  • 参考文献
  • 致谢
  • 个人简历、在学期间发表的学术论文与研究成果
  • 相关论文文献

    • [1].复杂地形上非恒定浅水二维流动数值模拟[J]. 水动力学研究与进展A辑 2013(01)
    • [2].无质量误差的具有复杂地形的浅水问题[J]. 水动力学研究与进展A辑 2008(05)
    • [3].一种基于并行计算的美式期权定价方法[J]. 辽宁石油化工大学学报 2013(02)

    标签:;  ;  ;  ;  

    浸没边界法的高效半隐式格式
    下载Doc文档

    猜你喜欢