[论文解读] Modeling Elasstic Shells Immersed in Fluid

Edward Givelberg|ArXiv.org|Jan 27, 2003

Lattice Boltzmann Simulation StudiesEngineering参考文献 10被引用 19

一句话总结

本文提出了一种基于浸入边界框架的数值方法，通过基尔霍夫-洛夫理论和平面应力壳体理论将该框架扩展至弹性壳体，其公式推导基于微分几何。该方法成功模拟了基底膜原型中的行进波，再现了耳蜗中从基底到尖端的波传播现象，验证了其在完整3D耳蜗建模中的适用性。

ABSTRACT

We describe a numerical method to simulate an elastic shell immersed in a viscous incompressible fluid. The method is developed as an extension of the immersed boundary method using shell equations based on the Kirchhoff-Love and the planar stress hypotheses. A detailed derivation of the shell equations used in the numerical method is presented. This derivation as well as the numerical method, use techniques of differential geometry in an essential way. Our main motivation for the development of this method is its use in the construction of a comprehensive three-dimensional computational model of the cochlea (the inner ear). The central object of study within the cochlea is the ``basilar membrane'', which is immersed in fluid and whose elastic properties rather resemble those of a shell. We apply the method to a specific example, which is a prototype of a piece of the basilar membrane and study the convergence of the method in this case. Some typical features of cochlear mechanics are already captured in this simple model. In particular, numerical experiments have shown a traveling wave propagating from the base to the apex of the model shell in response to external excitation in the fluid.

研究动机与目标

• 开发一种计算方法，用于模拟浸入于粘性不可压缩流体中的弹性壳体，其动机源于对耳蜗机械行为建模的需求。
• 将浸入边界方法扩展至包含基于基尔霍夫-洛夫和面内应力假设的壳体动力学，其理论基础为微分几何。
• 为完整的3D耳蜗计算模型构建基础，重点研究基底膜在听觉信号处理中的作用。
• 通过在脉冲流体激励下对基底膜原型模型进行收敛性和数值实验，验证该方法的收敛性和准确性。
• 在简化的数值设置中重现耳蜗力学的关键特征，如从基底到尖端的行进波传播。

提出的方法

• 该方法通过引入基于基尔霍夫-洛夫和面内应力假设推导出的壳体方程，扩展了浸入边界框架，以模拟弹性壳体的面内行为和弯曲行为。
• 采用微分几何技术推导壳体方程，从而能够准确表示流体-结构相互作用中曲面薄弹性结构的几何特性。
• 流体由不可压缩的纳维-斯托克斯方程控制，而弹性壳体则被建模为沿其中面分布的力元网络。
• 通过重新映射过程在欧拉网格的流体与拉格朗日节点的壳体之间传递力，确保流体-结构耦合的守恒性和一致性。
• 时间积分采用分数步长方法，时间步长受CFL条件和浸入边界刚度的约束。
• 通过细化流体网格和时间步长分析收敛性，当时间步长与流体网格宽度成比例时，显示出一阶收敛性。

实验结果

研究问题

• RQ1 能否将浸入边界方法扩展至在粘性流体流动中模拟具有曲率和面内应力的弹性壳体？
• RQ2 所得到的数值方法是否能再现实验观察到的基底膜从基底到尖端的行进波？
• RQ3 当流体网格和时间步长被细化时，该方法的收敛行为如何？
• RQ4 该方法能否捕捉到在流体激励下具有空间变化柔度的壳体的机械响应？
• RQ5 该方法在大规模模拟完整3D耳蜗几何结构时是否具备可扩展性和稳定性？

主要发现

• 当时间步长按比例减小至流体网格宽度时，该数值方法表现出一阶收敛性，证实了其稳定性和一致性。
• 在脉冲流体力作用下，行进波从壳体模型的基底向尖端传播，与冯·贝凯希实验观察到的耳蜗波力学特征一致。
• 波的传播由壳体柔度的空间变化驱动，基底附近刚度较高的区域首先引发波，波随后向柔度更大的尖端区域传播。
• 该方法成功捕捉了壳体的位移动力学，可视化结果展示了流体脉冲向下作用后弹性反弹及波的传播过程。
• 该方法在基底膜原型模型上得到验证，并随后用于构建完整的3D耳蜗模型，尽管计算需求较高，但证明了其可行性。
• 模拟所需的时长约为50纳秒，对应100 μm的流体网格，与CFL条件和边界刚度带来的稳定性约束一致。