跳到主页内容
美国国旗

美国政府的官方网站

Dot政府

gov意味着它是官方的。
联邦政府网站通常以.gov或.mil结尾。之前分享敏感信息,确保你在联邦政府政府网站。

Https系统

该站点是安全的。
这个https://确保您连接到官方网站,并且您提供的任何信息都是加密的并安全传输。

访问密钥 NCBI主页 MyNCBI主页 主要内容 主导航
.2008年8月10日;227(16):7587-7620.
doi:10.1016/j.jcp.2008.04.028。

复杂三维刚体流固耦合模拟的曲线浸没边界法

附属公司

复杂三维刚体流固耦合模拟的曲线浸没边界法

伊曼·博拉兹贾尼等。 计算机物理杂志. .

摘要

Ge和Sotiropoulos的sharp-interface CURVIB方法[L.Ge,F.Sotiropoolos,《求解复杂浸没边界曲线域中三维非定常不可压Navier-Stokes方程的数值方法》,计算物理杂志225(2007)1782-1809]被扩展到模拟流体-结构相互作用(FSI)涉及复杂三维刚体承受大结构位移的问题。FSI求解器采用分区FSI求解方法,实现了松耦合和强耦合策略。浸没体和流体之间的界面用拉格朗日网格离散,并用显式前跟踪方法跟踪。为了快速识别背景网格和运动物体之间的关系,提出了一种高效的光线追踪算法。对两个FSI问题进行了数值实验:弹性安装圆柱体的涡激振动和生理条件下双叶机械心脏瓣膜的流动。对于这两种情况,计算结果与基准模拟和实验测量结果都非常一致。数值实验表明,结构的特性(质量、几何形状)和局部流动条件在决定FSI算法的稳定性方面起着重要作用。在一定条件下,即使采用强耦合FSI,也会产生无条件不稳定的迭代格式。然而,对于这种情况,将强耦合迭代与欠松弛结合起来,再加上艾特肯加速技术,可以有效地解决稳定性问题。给出了一个理论分析来解释数值实验的结果。结果表明,附加质量与结构质量的比值以及流体对结构施加的力或力矩的局部时间变化率的符号决定了FSI算法的稳定性和收敛性。还阐明了欠松弛的稳定作用,并导出了稳定所需欠松弛系数的上限。

PubMed免责声明

数字

图1
图1
3D BMHV几何形状包括直主动脉中的外壳和小叶(左)以及小叶角度的定义(右)。传单铰链轴平行于x个-位置处全局参考系的轴(小时,z(z)小时)在YZ平面中。
图2
图2
显示sharp-interface浸没边界法思想的草图。
图3
图3
显示二维点多边形测试的光线投射方法的草图。
图4
图4
显示3D传单的光线投射方法的草图。
图5
图5
(a) 为双叶机械阀创建的边界框。阀门用2048个三角形离散化。(b) 边界框被划分为较小的控制单元。
图6
图6
自由流中弹性安装的圆柱体
图7
图7
网格用于弹性安装圆柱周围的二维涡激振动(VIV)研究。
图8
图8
作为函数的最大气缸位移U型红色来自当前(圆圈)和Ahn和Kallinderis[28](填充正方形)模拟。流动条件:重新= 150,M(M)红色= 2.
图9
图9
比较LC-FSI和SC-FSI获得的气缸运动学。流动条件:重新= 150,U型红色= 4,M(M)红色= 2.
图10
图10
圆柱体附近的瞬时涡度轮廓。流动条件:重新= 150,U型红色= 4,M(M)红色= 2.
图11
图11
六个弹性安装气缸的VIV。两张非平面涡度等值线快照。流动条件:重新= 100,U型红色= 4,M(M)红色= 2.
图12
图12
BMHV模拟的计算网格。a) 侧视图显示用于离散空主动脉的背景曲线网格,以及用于离散传单和外壳的非结构化网格。b) 背景网格的横截面视图。为了清晰起见,只绘制了三分之一网格。
图13
图13
BMHV模拟。入口处规定的生理输入流量波形。
图14
图14
BMHV模拟,在一个时间步长内,作用在具有恒定系数欠松弛和Aitken加速度的传单上的计算力矩的SC-FSI收敛历史的比较。
图15
图15
BMHV模拟(精细网格)。计算的小叶运动学(实线)与实验观察结果的比较[1](圆)。
图16
图16
精细网格上的BMHV模拟(左)与Dasi等人的PIV测量(右)进行了比较[1]。心脏周期内四个不同时间点阀门中平面上的瞬时非平面涡度轮廓。每个图右侧所示流入波形上的点表示循环期间的时间瞬间。轮廓级别相同。
图17
图17
BMHV模拟(精细网格)。如图16所示,在一个心动周期内的相同四个时间瞬间,通过q标准的等值面可视化的瞬时涡结构。
图18
图18
细网格(蓝色实线)和粗网格(红色虚线)上的BMHV模拟结果。阀门中间平面上的瞬时非平面涡度等值线,涡度大小等于1。每个图列底部显示的流入波形上的圆点表示心脏周期中的时间瞬间。
图19
图19
BMHV模拟。SC-FSI在细(黑色)和粗(红色)网格上获得的开放阶段上下叶运动学的比较。fg和cg分别指细网格和粗网格。
图20
图20
BMHV模拟。通过SC-FSI和LC-FSI比较闭合阶段的上下叶运动学。
图21
图21
BMHV模拟。根据小叶角度(顶部)、角速度(中部)和力矩系数(底部),比较闭合阶段其中一个小叶的LC-FSI和SC-FSI。
图22
图22
欠松弛系数的比较α通过Aitken的SC-FSI加速技术在早期打开和关闭阶段进行计算。
图23
图23
在静止流体(A)和(b)、(c)的自由流中绕铰链旋转的板。板通过(a)恒定外部力矩和(b)、(c)流体施加的力矩旋转。(b)中的结构和流向类似于开启阶段,(c)类似于关闭阶段。

类似文章

引用人

工具书类

    1. Dasi LP、Ge L、Simon HA、Sotiropoulos F、Yoganathan AP。轴对称主动脉中双叶机械瓣的涡动力学。流体物理学。2007;19:067105.
    1. Ge L,Sotiropoulos F.求解复杂浸没边界曲线域中三维非定常不可压Navier-Stokes方程的数值方法。计算物理杂志。2007;225:1782–1809.-项目管理咨询公司-公共医学
    1. Donea J、Giuliani S、Halleux J。瞬态动态流体-结构相互作用的任意拉格朗日-欧拉有限元方法。《应用力学与工程计算机方法》1982;33(1–3):689–723.
    1. 佩斯金CS。心脏瓣膜周围的流动模式:一种数值方法。计算物理杂志。1972;10:252–271.
    1. Farhat C,Lesoinne M,LeTallec P.非匹配离散界面流体/结构相互作用问题的载荷和运动传递算法:动量和能量守恒,最优离散化以及在气动弹性中的应用。应用力学和工程中的计算机方法。1998年;157(1):95–114.

LinkOut-更多资源