全欧拉法流固耦合问题

摘要

在这项工作中，我们提出一个关于流固耦合的一个完全欧拉结构(FSI）问题即耦合不可压缩Navier-Stokes超弹性固体方程。完全欧拉结构是一个单一变量结构耦合问题。与此相反的任意拉格朗日欧拉（ALE）坐标也是行之有效的，但是完全欧拉结构，包含两个子问题，液体和固体问题。在欧拉坐标，这个概念绕过与ALE各种困难坐标的联系，因为没有其他人的领域成果可以使用，所以该结构是一种首创全新的方法。本方法主要研究该固体的变形，它作为一种固体变形研究的的扩展，建立在初始点的设置，以此检测连接点的位置。由于涉及到大变形，所以本文尽可能的利用了像固体接触边界的变化或其他固体领域的研究。1介绍

我们为流固耦合问题提供了一种完全统一的可变有限单元法。重点强调于大变形结构领域的应用，处于流域和该流域边界接触的结构的自由运动，和其他自接触结构的应用。这项工作中出现的方程是Eulerian-Eulerian型的，并且这种奇特的方法第一次出现是由Dunne引入的[14,15]。

现在有数不清的不同方法来建模和模拟流固耦合问题。在这些方法中我们专注于单片模型，这种模型整个问题被描述为一种包含固体和流体表面的双系统。单片模型允许含糊的描述方案，大的时间步长并且提供使用基于错误假设的可能和最优化方法。他们已经很好的适用于大流体密度问题中如血液动力学。当流体问题被本质的描述为一种混合欧拉或朗格朗日的结构，材料的描述通常是固体问题的基础。所有的描述流固耦合问题的单片模型不管怎么样都要符合这两种结构。

在拉格朗日法或专一的拉格朗日法中，流体问题被定位于一种涉及域的匹配的位置。经典的方法是ALE法，见[29.4.35]或变空间域/稳态空间法（DSD/SST），见实例[53,51]。这些方程有相同之处，运动学和动力学结合很容易嵌入路径空间并有多种技术支持。拉格朗日法的缺点是流体问题的转化可在大变形或大的固体运动中破坏。

欧拉-拉格朗日法为解决流体问题使用一种欧拉混合计算单元，为解决固体

问题使用朗格朗日单元。通过使用力密度法结合这两种结构的应用于大边界[43]或大表面[39]。其它介绍的方法在于额外的表面不同。这些例子都是虚构的域法，其他最近的方法都是基于扩展有限单元法的。表面混合法的调查由Felippa[17]

给出。欧拉-拉格朗日法是一种表面捕捉的方法。这种表面不是修复的欧拉单元

的一部分。一种早期的表面捕捉是流体体积法[30]，其中的轨迹函数在流域中取值1在固体中取值0。轨迹函数与表面的速度有关。找到一种合适的轨迹函数的数值近似以此达到从流域到固体的转换是这种方法的一个困难。另一种可能的表面捕捉方法首先被用于构造多相流，其中的功能函数表示的是域中的任意一点到表面的距离。界面本身被设置为零等级的轮廓函数。基于等级设置方案，存在一种非常有效和简单的数字方案可以描述自由边界流固耦合。通常表示边界的等级设置方案缺乏表示尖角的能力。欧拉拉格朗日方法的普遍问题是靠近接口的近似性。由于流体单元被界面分隔，那么解可以是不连续的（或至少不可微）跨越界面，与标准近似的有限单元是困难的。在此，扩展有限元方法有助于提高精度。

最后，对于大变形问题，结构问题的描述欧拉法是可取。使用混合欧拉边界网格，液体和固体之间的界面可以自由地在域内移动。欧拉法一直是界面捕捉的类型。邓恩[14]和[15]首次使用了一个具有欧拉-欧拉模型的不可压缩流体与超弹性固体相互作用的方法。这里，表面捕捉由初始点集（IPS）完成，矢量场ΦIPS 用于传输完整的参考坐标系。IPS法能够捕捉尖锐界面。让全欧拉流体有限差分方法基于流体体积的方法[30]是由Sugiyama等人介绍[50]。这种方法中不可压缩的流体和结构之间的耦合通过校正压力迭代得以解决。全欧拉流体结构基于水平集功能的交互方式是由他和乔介绍的。等级设置函数代表所述移动结构域需要四种。其中两种所扮演的角色类似于由Dune介绍的IPS法。

在这项工作中引入的全欧拉法是基于IPS法的延伸。流体和固体问题在欧拉方程中给出，该边界条件在函数空间通过可变载荷的平衡实现。用于捕捉的界面处的固体变形将扩展到只有一小层的流体域中。流动问题没有进行任何改造的建模，并且没有额外的变量引入使得该方法很有效。

在第二部分中，我们简短提出要求的符号。第三部分致力于在欧拉坐标方程下解决流固耦合问题。细节描述和解决方案在第四部分给出。最后，在第五部分我们介绍展示这个解决流固耦合问题的全欧拉法的不同数值计算例子。在第六部分我们给出结论。

符号规定

2 符号规定

令Ω为二维空间向量。在初始时刻t=0，这个空间划分为非重叠分区，分别为流体部分f Ω和固体部分s Ω，并且还有一个共同的接触面i Γ。我们会考虑的问题，在这里域分区随时间变换，但总的区域Ω是固定的：

在t=0时刻我们称f Ω和s Ω为初始值。在本文的流固耦合问题中s Ω为朗格朗日观测点，而()s t Ω指的是欧拉观测点。为简单起见，我们只考虑狄氏边界条件问题。在功能区我们在Ω域中使用勒贝格空间2L ，在运动的子域和索博列夫空间中2L 是弱衍生设置。我们定义2L 在两个域空间边界上为线性的。通过1H 表

示运动的子域和索博列夫空间中2L 的弱衍生设置。并且用1

0H 表示

1H 功能的0号路径。

3. 流固耦合问题的欧拉方程

在本节中，我们得出方程的耦合系统描述一个结合用斯托克斯相互作用或纳维 - 斯托克斯流体与圣维南基尔霍夫类型描述的的弹性结构。所有方程给出配方，这使得欧拉方程和拉格朗日方程之间容易转化为协调结构。

3.1 流体问题的欧拉描述

和拉各朗日方法不同，欧拉方法不同，欧拉方法的着眼点不是流体质点而是空间点。设法在空间中的每一点上描述出流体运动随时间的变化状况。如果，每一点的流体运动都已知道，则整个流体的运动状况也就清楚了。因为不同时刻将有不同流体质点经过空间某固定点，所以站在固定点上就无法观测和记录掠过的流体质点以前和以后的详细历史。也就是说我们无法象拉各朗日方法那样直接测量出每个质点的位置随时间的变化情况。虽然如此，不同时刻经过固定点的流体质点的速度是可以测出的，这样采用速度矢量来描写固定点上流体运动的变化状况就是十分自然的了。考虑到上面所说的情形，欧拉方法中流体质点的运动规律数学上可表示为下列矢量形式：

()

v=v r,t （1）

在直角坐标系中有： (),,,u u x y z t = (),,,v v x y z t = (),,,w w x y z t = （2） 要完全描述运动流体的状况还需要给定状态函数压力、密度、温度等 (),,,p p x y z t = (),,,x y z t ρρ= (),,,T T x y z t = （3）