IB-LBM浸入边界法_第1页
IB-LBM浸入边界法_第2页
IB-LBM浸入边界法_第3页
IB-LBM浸入边界法_第4页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、精选优质文档-倾情为你奉上浸入边界法1 IB-LBM的基本数学构造Hofler和Schwarzer在2000年提出了一种在对NS方程求解的有限差分方程中增加一项体积力的IBM来模拟充满粒子流场的固体粒子的运动。其粘性不可压流的控制方程表示为： （1） （2）f：作用在流体上的外力上述两个公式在IB-LBM中可有格子Boltzmann方程所代替，即：（3）（4）（5）F表示作用在物体表面上的力，其物理过程包括：碰撞过程：（6）迁移过程：（7）由此可知，碰撞仅仅是局部的，但迁移和邻近点皆有关系。平衡态分布函数为：（8）其中为声速，式中表示各离散速度方向的加权系数，其值分别为： （9）该模型中，流体

2、的松弛时间和运动粘性系数的关系可表示为：（10）1.1 IB-LBM速度修正法在LBM中，同时满足NS方程的宏观量是通过Chapman-Enskog展开得到的。在速度修正法里其表达式方程有如下表示：（11） （12）其中，为中间层速度，为修正速度。结合方程（11），（12），方程（4）可变为： （13）在传统的IB-LBM中，f由胡克定理、直接力法或者动量交换法来求得。其中胡克定理表达式为： （14） （15）其中称为狄拉克插值函数，其作用是将作用在物面上的拉格朗日坐标系下的外力插值到欧拉坐标系下对流体点的作用力，描述了流体和固体之间相互作用的关系。如果单从上述方程来求外力f的话，流体会穿越物

3、面，导致质量不守恒。为了解决流体穿透的问题，该方法在求解f时做了相应处理，即假设f为非预知力，由满足无滑移边界条件的速度修正场求得。方程（13）中，设为欧拉坐标系下的修正速度，有：（16）其中，可由方程（11）求得。依此，由方程（14）、（15）类似推导出的求法。假设表示拉格朗日坐标系下的修正速度，相应地，则有：（17）其中狄拉克插值函数可以光滑地近似为： (18)Peskin提出表达式为：（19）因此可应用方程（18）将方程（17）离散：（20）由此类推可得：？（21）将方程（20）代入方程（16）中，则有（22）将方程（22）代入方程（21）可得（23）将方程（23）改写成：（24）其中，于是可将方程（24）简化为：宏观量压强和密度分别为：， （25）作用在流体上的外力可由方程（12），求得（26）计算步骤如下：（1）给定初始状态值，计算出矩阵A的值,并求得矩阵A的逆A-1;(2) 给定的初始值为0，即，应用式(3)求时的密度分布函数值，通过式(11)和式(25)，计算该时刻的流场宏观量，如宏观量、等;（3）求解方程组，获得物面边界点的修正速度，再通过式(20)将求得的物面边界点的修正速度插值到欧拉点上，也即求得欧拉点的修正速度；（4）利用式（16）修正流场中欧拉点的速度，同时应用式(2

人人文库> 全部分类> 教育资料 > 辅导培训

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论