沉浸边界法的改进及其在鱼游推进机理研究中的应用研究.pdfVIP

第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会文集 沉浸边界法的改进及其在鱼游推进机理 研究中的应用· 潘定一邓见邵雪明 浙江大学航空航天学院力学系，流体传动及控制图家重点实验室，杭州，310027 摘要：沉浸边界法是一种基于非贴体笛卡尔网格的计算方法，在处理含有运动物体流场的流 动问题时，有着较高的计算效率且前处理所耗费的时间较少．本研究给出了一种沉浸边界法 的改进，对于流场中的运动物体边界，通过一系列的离散点来描述，同时在这些离散点和网 格点之间建立了一种较简洁的变量分布过程．在保证计算准确性的同时，简化了对沉浸边界 的处理过程．文中同时介绍了此方法在鱼游推进机理研究中的应用，主要包括：尾鳍拍动、 双鱼游动、鱼群巡游和整鱼游动等。 关键词：沉浸边界法；鱼游推进；计算流体力学；仿生力学 1引言 在自然界和工程技术领域，复杂边界和运动边界的问题普遍存在，对复杂动边界处理方 法的研究是目前计算流体力学的一大热点。对于复杂动边界问题，通常有基于坐标变换和网 格映射的有限差分法和基于非结构网格的有限体积法。近年来，非贴体网格技术被越米越广 Domain 泛的应用。这其中主要有：虚拟区域法(Fictitious Interface Grid BoundaryMethods)、浸没界面法(Immersed Methods)等。相比于传统的贴体网格方法，非贴体网格方法通常采用结构化网格，针对动边界 问题，在每一时间步，不需要重新生成网格。这大大地减少了前处理的工作量，也提高了每 一时间步的计算效率。 沉浸边界法的主要思想是：认为流体充满整个计算区域，即浸没物体的区域也被流体所 占据，通过在动量方程中加入一个虚拟附加力来满足流体在浸没物体上的边界条件，一般情 况下，只有在浸没物体边界附近，虚拟力不为零。沉浸边界法最早由PeskinIlJ在1972年提出， 用于模拟血液在心脏瓣膜附近的流动。对于沉浸边界法的研究主要集中于两个方面：①虚拟 附加力的求解形式；②插值方法(即如何使物理量在浸没边界和计算网格之间建立联系)。 Goldstein12J、Mohd．Yusof【3J、Fadlun I7 I引、Kim【引、Balar丛lol和MittalJ分别基丁二以上两方面对沉 ‘基金项目：高等学校博七学科点专项科研堆金(20070335066) 一58— 第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会文集 浸边界法进行了发展改进，提出了反馈力(FeedbackForce)、直接力(Di∞ctForce)和Ghos沁e¨ 等处理方法，并将沉浸边界法与有限差分法和有限体积法结合使用。 本研究作者对沉浸边界法中处理边界和网格点之间的插值及虚拟附加力的求解进行了一 些改进，在一定程度上简化了求解过程，并保证了相应的精度。同时将改进后的方法应用于 鱼游运动的模拟研究。本研究将分别对不同模型的数值模拟结果进行简要介绍。 2改进的沉浸边界法 2．1控制方程 带虚拟附加力的黏性不可压N·S方稗如卜．(，为附加力，在沉浸边界附近时，不为零)： 、 (1) 一十II，·V-l，=一V尸+一V-，+一 詈邶∽y一尸+去V2y+F8t Re V．y=0 (2) 2．2沉浸边界的处理 对于沉浸边界，通常选取一系列的离散控制点x。(即拉格朗日点)，用丁．描述边界的形状。 这样的优点是，可以描述任意的边界形状，特别是无法用解析公式表示的不规则边界。在这 些控制点上计算相应的附加力，在此，称之为拉格朗日力，(一)，然后通过特定的分布函数 将拉格朗日力分布剑附近的网格