基于fluent阻力计算

1、百度文库-好好学习，天天向上 基于fluent的兴波阻力计算本文主要研究内容本文的工作主要涉及小型航行器在近水面航行时的绕流场及兴波模拟和阻力的数值模拟两个方面。在阅读大量文献资料的基础上，通过分析、比较上述领域所采用的理论和方法，针对目前需要解决的问题，选择合理的方法加以有机地综合运用。具体工作体现在以下几个方面：1本人利用FLUENT软件的前处理软件GAMBIT自主建立简单回转体潜器模型，利用FLUENT求解器进行计算，得出在不同潜深下潜器直线航行的绕流场、自由面形状及阻力系数的变化情况。2通过对比潜器在不同潜深情况下的阻力系数，论证了增加近水面小型航行器的深度可以有效降低阻力。通过对模型

2、型线的改动，为近水面小型航行器的型线设计提供了一定的参考。通过改变附体形状和位置计算了附体对阻力的影响程度，为附体的优化设计提供了一定的依据。计算模型航行器粘性流场的数值计算理论水动力计算数学模型的建立根据流体运动时所遵循的物理定律，基于合理假设（连续介质假设）用定量的数学关系式表达其运动规律，这些表达式成为流体运动的数学模型，它们是对流体运动的一种定量模型化，称为流体运动控制方程组。根据控制方程组，结合预先给定的初始条件和边界条件，就可以求解反映流体运动的变量值，从而实现对流体运动的数值模拟预报，形成分析报告。基于连续介质假设的流体力学中流体运动必须满足要遵循的物理定律:1）质量守恒定律2）

3、动量守恒定律3）能量守恒定律4）组分质量守恒方程针对具体研究的问题，有选择的满足上述四个定律。船体的粘性不可压缩绕流运动，如果不考虑水温对水物理性质的影响，水的密度和分子粘性系数都是常数，同时没有能量的转换，就仅仅需要满足质量守恒定律、动量守恒定律。在满足这些定律下所建立的数学模型称为Navier-Stokes方程。另外，自由液面的存在也需要建立合适的数学模型。本文是利用FLUENT进行数值模拟，而软件里面关于自由液面模拟是用界面追踪方法的一种一流体体积法（VOF），基于该方法所建立的数学模型称为流体体积分数方程。另外，高雷诺数下的水动力问题还需要考虑粘性不可压缩流体的湍流运动。对于湍流运动的

4、数值模拟一直是流体力学数值计算的一个难点。直接数值模拟（DNS）目前还仅仅在院校中研究，而且也仅限于二维流体问题。大涡模拟（LES）向工程应用的过渡似乎还没有完成，并且就高雷诺数问题而言，对计算机硬件要求很苛刻。目前，从算法的可行性、硬件要求的可实现性、完成任务所消耗时间和人力等方面看，基于湍流模型的数值计算更为工程实际所接受。本章将会对各种湍流模型加以介绍。百度文库-好好学习，天天向上百度文库-好好学习，天天向上 - 粘性不可压缩流体流动数学模型连续方程任何流动问题都必须满足质量守恒方程即：连续方程。根据连续介质假设，单位时间内流体微团的质量变化等于同时间间隔内进入微团的总净质量。按照这一定

5、律，连续方程数学表达式写为：4_匸忙I+鼠“_田卢门=|（）以上是在笛卡尔直角坐标系下表示，上面给出的是瞬态可压流体连续方程。由于对于潜艇粘性流场介质的不可压缩，密度P为常数，引入散度算子，则方程（）变成为：（）式中：速度矢量V=u,v,w。上式为粘性不可压缩流体运动的连续方程。动量方程动量守恒方程也是任何流动系统都必须满足的定律。根据牛顿第二定律，流体微团中流体的动量对时间的变化等于微团所受外力之和，即：()鹫&cexdyCZ式中，p代表流体微团所受的压力；Txx、Txy、Txz等是因分子粘性作用而产生的作用在流xxxyxz体微团表面上的粘性应力t的分量；Fx、Fy、Fz表示直角坐标系下三个

6、方向上流体微团的体积力分量，如果体积力只有重力，且Z竖直向上，则Fx=Fy=O,Fz=-pg。式（）（）是对任何类型的流体（包括非牛顿流体）均成立的动量守恒方程。本文研究的范围属于牛顿流体，故粘性应力t与流体的变形率成比例，有：式中，卩是动力粘度系数，入是第二粘度，一般可取入=-2/3,将（）代入式（）（）式（）就是动量守恒方程。方程（）和（）组成了控制粘性不可压缩流体运动的基本数学模型。对于低雷诺数的层流运动，上述方程组已经可以确切描述流体运动。但湍流流动以脉动的速度场为基本特征，各速度在时间和空间上变化很快，给流场的数值模拟带来很大困难再则，湍流是一种极度复杂的物理现象，包含无规律性，扩散

7、性，三维涡旋波动及耗散。在实际工程计算中要对湍流进行数值模拟代价十分高昂。然而研究表明，大尺度涡在流体运动中起主要作用。由此可见，若采用时间平均、集合平均或者其他人工处理方法略去小尺度运动，将小尺度运动模型化后代入大尺度中，从而替代求解原有瞬时控制方程，就会花费较小的计算代价获得较高精度的数值解。以此为出发点，提出了将速度分解成平均值和脉动值，则瞬时速度分量u可以表达为：将式（）代入（）和（）再对时间积分就会得到下面的平均流方程。()()迪+他妆一璽+2（赵竺+漏H耳dtdx.fixcXjf方程（）是时均形式的连续方程，方程（）是时均形式的Navier-Stokes方程。方程（）为Reynol

8、ds应力。由于式（）采用的是Reynolds平均法，因此方程（）被成为Reynolds平均Navier-Stokes方程（Reynolds-AveragedNavier-Stokes简称RANS方程）。有式（）和（）组成的方程组共有五个方程（RANS方程实际是3个）现在新增了6个Reynolds应力，再加上原来4个时均未知量，总共9个未知量，因此，方程组不封闭，必须引入新的湍流模型（方程）才能使方程组（）和（）封闭。湍流模型为了使雷诺平均N-S方程（RANS方程）封闭可解，要根据湍流的运动规律来寻求附加的条件和关系式，这就形成了不同的湍流模型。在FLUENT计算软件中可以供选择的湍流模型有：一

9、方程模型Spalart-Allmaras（S-A）、两方程模型k-（包括Sk-、RNGk-）和k-s（包括Sk-s和SSTk-s）以及雷诺应力模型（RSM），下面将本文所用到的四种湍流模型加以介绍。标准k-模型（Sk）标准k-模型是典型的两方程模型，该模型是目前应用最广泛的湍流模型。k和是两个基本未知量，与之相对应的输运方程为：()()()湍流动能k方程为：湍流耗散率&方程为：式中的湍流涡粘度卩t可表示为：（其中=，为一常数）式中：Gk式由于平均速度梯度引起的湍动能k的产生项，Gb是由于浮力引起湍动能k的产生项，Ym代表可压缩湍流中脉动扩张的贡献，C花、C2和C3为经验常数，ok和叭分别是与湍

10、动能k和耗散率&对应的湍流普朗特数，Sk和S是用户定义的源项。模型常数C比、C2厂C口、ok、叭的取值为：C1尸，C2=，5=,ok=,oe=RNGk-模型RNGk-湍流模型是由Yakhot及Orzag提出的，该模型中的RNG是英文“renormalizationgroup”的缩写。在RNGk-湍流模型中，通过在大尺度运动和修正后的粘度项体现小尺度的影响，而使这些尺度运动有系统地从控制方程中去除。所得到的k方程和&方程，与标准k-模型非常相似:湍流动能k方程：()湍流耗散率&方程:()与标准k-模型相比较发现，RNGk-模型的主要变化：通过修正湍动粘度，考虑了平均流动中的旋转流动情况；在&方程

11、中增加了一项，从而反映了主流的时均应变率E耳，这样，RNGk-模型中产生项不仅与流动情况有关，而且在同一问题中也还是空间坐标的函数。模型常数C1，C2由RNG理论：C1尸，C2=，其他常数：C卩=,ok=,og=Sk-们模型本文采用的Sk-3湍流模型是基于湍流动能k和特殊湍流动能耗散率3的输运方程建立起来的经验公式。是由Wilcox在1998年提出的对原k-3模型的改进模型。湍流动能k方程为：百度文库-好好学习，天天向上百度文库-好好学习，天天向上- rat二3联_务曲)=三厂上竺CtOXf禺GXj()特殊耗散率3方程为:at()rkru表示k、3的有效扩散率，表示为:()上式中a*=Bt/3

12、,Ret为雷诺数:Gk、03分别为湍流动能k和湍流耗散率3的普朗特数，湍流涡粘度卩t可表示为:()(&a*为低湍流雷诺数修正系数:,J広“+Rr/&()气14阿绮丿以上各式中的常数取值为：口；二1，P=0J72,crt=2.Dtcra-2.0剪切应力输运k-们模型(SSTk-m)SSTk-e湍流模型由Menter提出，该模式的湍流动能方程和湍流耗散率方程与标准Sk-e模型的形式相似：湍流动能k方程为：()特殊耗散率3方程为:rkru和卩t见式，常数ok、G3表示为:I丽齐口57石式中F1是混合函数：耳=ranht()()()（）Z；=max2/?()式中：Gk式由于平均速度梯度引起的湍动能k的

13、产生项，Gu是由于特殊湍流动能耗散率s的产生，D为横向扩散项，Yk、Ys表示湍流k、s的消耗，Sk和Se是用户定义项。边界条件边界条件类型简介流体在运动的过程中会受到边界的限制，反映到物理模型上，就是要给控制方程加一些关于变量q、P、k、&相应的边界条件。最常见的线性边界条件有两大类：第一类边界条件（Dirichlet条件）和第二类边界条件（Neumann条件）。前者描述的是计算区域的边界或部分边界上变量的值，后者则描述边界上变量梯度的法向分量值，即：Dirichlet条件：=b在边界上bNeumann条件：nV=n在边界上式中为任意的物理量，n表示物体表面的单位外法线矢量，b为给定的边界上的

14、数值，n为给定的V巾在边界上的法向分量。对于潜器粘性绕流，入流边界是一种人工边界，它不由物体的性质决定，因而不是固定不变的，它需要取得离潜器表面足够远才能尽量地反映真实情况。入口处边界条件属于Dirichlet条件：其速度是预先给定的，一般是均匀来流条件，湍动能k和耗散率s也是预先给定的。出流边界条件则是虚拟的，出流边界到艇尾的距离也要合理确定以消除对流场计算的影响。对于粘性流动，在固壁边界（如艇体表面）须满足对速度和湍动能k的无滑移边界条件，即：u=v=w=0，k=0然而在靠近壁面的区域，由于湍动能被强烈地耗减，耗散率达到最大值，在固壁上不易给出s的边界条件，因为它在壁面上不等于零。在与壁面

15、相邻的粘性子层中，由于粘性的影响，局部雷诺数变得很小。由于前述k-模型是一种高雷诺数模型，因而对粘性子层不再适用，一般采取Launder和Spalding提出的壁函数方法来处理。使用边界条件的注意事项1边界条件的组合在CFD计算域内的流动是由边界条件驱动的。从某种意义上说，求解实际问题的过程，就是将边界线或边界面上的数据，外推扩展到计算域内部的过程。因此，提供符合物理实际且适合的边界条件是极其重要的，否则，求解过程将很难进行。CFD模拟过程中迅速发散的一个最常见的原因就是边界条件选择的不合理。例如，只给定进口边界和壁面边界，而没有给定出口边界，那么，将不可能得到计算域的稳定解，CFD将越计算越

16、发散。这样的边界条件组合显然是不合理的。在使用出口边界时需要特别注意，该边界只在进入计算域的流动是以进口边界条件给定（如在进口给定速度和标量）时才使用，而且仅推荐在只有一个出口的计算域中使用。物理上，出口压力控制着流体在多出口间的分流情况，因此，在出口给定压力值要比给定出口条件（梯度为零）合理。将出口条件和一个或多个恒压边界结合使用是不允许的，因为零梯度的出口条件不能指定出口的流量，也不能指定出口的压力，这样将使问题不可解。2流动出口边界的位置选取如果流动出口边界太靠近固体障碍物，流动可能尚未达到充分发展的状态（在流动方向上梯度为零），这将导致相当大的误差。一般来讲，为了得到准确的结果，出口边

17、界必须位于最后一个障碍物后10倍于障碍高度或更远的位置。对于更高的精度要求，还要研究模拟结果对出口位于不同距离时的影响的敏感程度，以保证内部模拟不受出口位置选取的影响3近壁面网格在CFD模拟时，为了获得较高的精度，常需要加密计算网格，而另一方面，在近壁面处为快速得到解，就必须将k-&模型与结合了准确经验数据的壁面函数法一起使用。要保证壁面函数法有效，就需要使离壁面最近的以内节点位于湍流的对数律层之中，即Y+必须大于（最好是在30500之间）。这就相当于给最靠近壁面的网格到壁面的距离设定了一个下限。但是，在流动的任意位置都使上述要求得到保证常常不太可能，典型的例子就是包含回流的流动。4随时间变化

18、的边界条件这类边界条件是针对非稳定问题而言的。就是说边界上的有关流动变量并不是一成不变，而是随着时间变化的。对于这类边界条件，需要将边界条件离散成与时间步长相应的离散结果，然后存储起来，供计算到相应的时间步时调用。这类边界条件一般是与初始条件一同给定的。k和的计算公式在入口、出口或远场边界流入流域的流动，FLUENT需要指定输运标量的值。在某些情况下流动流入开始时，将边界处的所有湍流量指定为统一值是适当的。在大多数湍流流动中，湍流的更高层次产生于边界层而不是流动边界进入流域的地方，因此这就导致了计算结果对流入边界值相对来说不敏感。然而必须注意的是要保证边界值不是非物理边界。非物理边界会导致你的

19、解不准确或者不收敛。对于外部流来说这一特点尤其突出如果自由流的有效粘性系数具有非物理性的大值，边界层就会找不到了。湍流强度I定义为相对于平均速度uavg的脉动速度的均方根。_avg小于或等于1%的湍流强度通常被认为低强度湍流，大于10%被认为是高强度湍流。从外界测量数据的入口边界，你可以很好的估计湍流强度。例如：如果你模拟风洞试验，自由流湍流强度通常可以从风洞指标中得到。在现代低湍流风洞中自由流湍流强度通常低到对于内部流动，入口的湍流强度完全依赖于上游流动的历史，如果上游流动没有完全发展者没有被扰动，你就可以使用低湍流强度。如果流动完全发展，湍流强度可能就达到了百之几。完全发展的管流的核心的湍

20、流强度可以用下面的经验公式计算：I=X（雷诺数Re）-1/8雷诺数Re=速度X当量直径X（密度三粘度）湍流尺度l是和携带湍流能量的大涡的尺度有关的物理量。例如在完全发展的管流中，l被管道的尺寸所限制，因为大涡不能大于管道的尺寸。L和管的物理尺寸之间的计算关系如下：L=Xl百度文库-好好学习，天天向上百度文库-好好学习，天天向上- - 其中L为管道的相关尺寸。因子是基于完全发展湍流流动混合长度的最大值的。在速度入流，压力边界条件中，需要输入如下两个物理量，其计算公式如下：k=(IW)2wall0.093/4k3/2=Lw其中是湍流模型中指定的经验常数(近似为)。wa为最大速度，I为湍流强度，L为

21、相当尺寸。数值计算方法网格生成用CFD方法进行流场计算时，首先要将计算区域离散化，即划分网格。网格是CFD模型的几何表达形式，也是模拟与分析的载体。计算网格的好坏直接影响到数值计算的可行性、收敛性以及计算精度。对于复杂的CFD问题，网格生成极为耗时，且极易出错，生成网格所需时间常常大于实际CFD计算的时间。因此，有必要对网格生成方式给以足够的关注。网格(grid)分为结构网格和非结构网格两大类。把节点看成是控制体积的代表。在离散过程中，将一个控制体积上的物理量定义并存储在该节点处。若节点排列有序，即当给出了一个节点的编号后，立即可以得出其相邻节点的编号。这种网格称之为结构网格(structur

22、edgrid)。结构网格是一种传统的网格形式，网格自身利用了几何体的规则形状。近几年来，还出现了非结构网格(unstructuredgrid)。非结构网格的节点以一种不规则的方式布置在流场中。这种网格虽然生成过程比较复杂，但却有着极大的适应性，尤其对具有复杂边界的流场计算问题特别有效。无论是结构网格还是非结构网格，都需要按下列过程生成网格：1建立几何模型。几何模型是网格和边界的载体。对于二维问题，几何模型是二维面对于三维问题，几何模型是三维实体。2划分网格。在所生成的几何模型上应用特定的网格类型、网格单元和网格密度对面或体进行划分，获得网格。3指定边界区域。为模型的每个区域指定名称和类型，为后

23、续给定模型的物理属性、边界条件和初始条件做好准备。网格生成是一个“漫长而枯燥”的工作过程，经常需要进行大量的试验才能取得成功。因此，出现了许多商品化的专业网格生成软件。如GAMBIT、TGrid、GeoMesh、preBFC和ICEMCFD等。此外，一些CFD或有限元结构分析软件，如ANSYS、一DEAS、NASTRAN、PATRAN和ARIES等，也提供了专业化的网格生成工具。这些软件或工具的使用方法大同小异，且各软件之间往往能够共享所生成的网格文件，例如FLUENT就可读取上述各软件所生成的网格。有一点需要说明，由于网格生成涉及几何造型，特别是3D实体造型，因此，许多网格生成软件除自己提供

24、几何建模功能外，还允许用户利用CAD软件(AutoCAD、Pro/ENGINEER)先生成几何模型，然后再导入到网格软件中进行网格划分。因此，使用前处理软件，往往需要涉及CAD软件的造型功能。方程离散数值计算是将描述物理现象的偏微分方程在一定的网格系统内离散，用网格节点处的场变量值近似描述微分方程中各项所表示的数学关系，按一定的物理定律或数学原理构造与微分方程相关的离散代数方程组。引入边界条件后求解离散代数方程组，得到各网格节点的场变量分布，用这一离散的场变量分布近似代替原微分方程的解析解。当前求解流体流动和传热方程的数值计算方法比较多，如有限差分法(FiniteDifferenceMetho

25、d),有限元法(FiniteElementMethod)、有限体积法(FiniteVolumeMethod)、边界元法、特征线法、谱方法、有限分析法和格子类方法等。每种数值计算方法有各自的特点和各自的适用范围，其中通用性比较好、应用比较广泛的是前3种。有限差分法有限差分法(FiniteDifferenceMethod,简称FDM)是数值解法中最经典的方法。它是将求解域划分成差分格式，用有限个网格节点代替连续的求解域，然后将偏微分方程(控制方程)的导数用差商代替，推导出含有离散点上有限个未知数的差分方程组。求差分方程组(代数方程组)的解，就是微分方程定解问题的数值近似解，这是一种直接将微分问题变

26、成代数问题的近似数值解法。这种方法发展较早，比较成熟，较多的用于求解双曲型和抛物型问题。用它求解边界条件复杂、尤其是椭圆问题不如有限元或有限体积法方便。有限元法有限元法(FiniteElementMethod，简称FEM)与有限差分法都是广泛应用的流体动力学数值计算方法。有限元是将一个连续的求解域任意分成适当形状的许多微小单元，并与各小单元分片构造插值函数，然后根据极值原理(变分或加权余量法)，将问题的控制方程转化为所有单元上的有限元方程，把总体的极值作为各单元极值之和，即将局部单元总体合成，形成嵌入了指定边界条件的代数方程组，求解该方程就得到各节点上待求的函数值。有限元法的基础是极值原理和划

27、分插值，它吸收了有限差分法中离散处理的内核，又采用了变分计算中选择逼近函数并对区域进行积分的合理方法，是这两类方法互相结合、取长补短发展的结果。它具有很广泛的适应性，特别适用于几何及物理条件比较复杂的问题，而且便于程序的标准化。对椭圆型方程问题又更好的适应性。有限元法因求解速度较有限差分法和有限体积法慢，因此，在商用CFD软件中应用不普遍。有限体积法有限体积法(FiniteVolumeMethod)又称控制体积法(ControlVolumeMethodCVM)。其基本思想是：将计算区域划分为网格，并使每个网格点周围有一个互不重复的控制体积；将待解微分方程(控制方程)对每个控制体积积分，从而得出

28、一组离散方程。其中的未知量是网格点上因变量。为了求出控制体积的积分，必须假定因变量在网格点之间的变化规律。从积分区域的选取方式来看，有限体积法属于加权余量法中的子域法，从未知解的近似方法看来，有限体积法属于采用局部近似的离散方法。简而之，子域法加离散，就是有限体积法的基本方法。就离散方法而言，有限体积法可视作有限元法和有限差分法的中间物。有限元法必须假定值在网格节点之间的变化规律(即插值函数)，并将其作为近似解。有限差分法只考虑网格点上的数值而不考虑值在网格点之间如何变化。有限体积法只寻求值节点值，这与有限差分法相类似；但有限体积法在寻求控制体积的积分时，必须假定值在网格点之间的分布，这又与有

29、限单元法相类似。在有限体积法中，插值函数只用于计算控制体积的积分，得出离散方程之后，便可忘掉插值函数；如果需要的话，可以对微分方程中不同的项采取不同的插值函数。综上所述，有限体积法是目前在流体流动和传热问题求解中最有效的数值计算方法。有限体积法也称为控制容积积分法，是20世纪六七十年代逐步发展起来的一种主要用于求解流体流动和传热问题的数值计算方法。有限体积法是在有限差分法的基础上发展起来的，同时它又吸收了有限元法的一些优点。有限体积法与有限元法和有限差分法一样，要对求解域进行离散，将其分割成有限大小的离散格式。在有限体积法中每一网格点按一定的方式形成一个包围该节点的控制容积矿，有限体积法的关键

30、步骤是将控制微分方程式在控制容积内进行积分。有限体积法获得的离散方程，物理上表示的是控制容积的通量平衡，方程中各项有明确的物理意义。有限体积法区域离散的节点网格与进行积分的控制容积分立。由于有限体积法的诸多优点，当今大多数CFD软件都采用它来离散求解，如PHOENICS、FLUENT、STAR-CD、NUMECA等。SIMPLE算法流场计算方法的本质是对离散后的控制方程组的求解。目前各种商用CFD软件普遍采纳的算法是压力耦合方程组的半隐式方法（SIMPLE算法），它属于压力修正法的一种。SIMPLE是英文SemiImplicitMethodforPressure-LinkedEquations

31、的缩写。SIMPLE方法由Patankar与Spalding于1972年提出，是种主要用于求解不可压流场的数值方法，也可用于求解可压流动。它的核心是采用“猜测一修正”的过程，在交错网格的基础上来计算压力场，从而达到求解动量方程（NavierStokes方程）的目的。SIMPLE方法的基本思想是对于给定的压力场（它可以是假定的值或者是上一次迭代计算所得到的结果），求解离散形式的动量方程，得出速度场。因为压力场是假定的或不精确的，由此得到的速度场一般不满足连续方程，所以，必须对给定的压力场加以修正。修正的原则是与修正后的压力场相对应的速度场能满足这一迭代层次上的连续方程。据此原则，把由动量方程的离

32、散形式所规定的压力与速度的关系代入连续方程的离散形式，从而得到压力修正项，由压力修正方程得出压力修正值。接着，根据修正后的压力场，求得新的速度场。然后检查速度场是否收敛，若不收敛，用修正后的压力值作为给定的压力场，开始下一层次的计算。如此反复，直到获得收敛的解。VOF方法由于本文研究的是近水面潜体周围的绕流场，所以必须考虑自由面问题。20世纪50、60年代，自由面的数值模拟方法有了实质的进展，其中美国LosAlamos的科学家们提出和发展的格子类方法，应用效果较好。著名的格子类方法主要有PIC方法、FLIC方法和MAC方法等。从20世纪70年代开始，自由面追踪的数值方法有了进一步发展，比较著名

33、的有标高法、线段法、VOF方法等。VOF（VolumeofFluid）方法是使用固定网格系统捕捉两种或两种以上不相掺混流体交界面的一种数值方法。在运动界面追踪问题的数值模拟方法中，VOF方法是最为重要的方法之一，它的特点是将运动界面在空间网格内定义成一种流体体积函数，并构造这种流体体积函数的发展方程，从而界面追踪问题的目的就是如何随着主场的模拟过程，通过流体输运，精细地确定该运动界面的位置、形状和变形方向，达到追踪的目的。VOF方法是美国学者Hirt和Nichols等人在MAC方法基础上提出的，它是一种可以处理任意自由面的方法。其基本原理是利用计算网格单元中流体体积量的变化和网格单元本身体积的

34、比值函数，来确定自由面的位置和形状。VOF方法追踪的是网格单元中流体体积的变化，而非追踪自由液面流体质点的运动，这与Harlow和疆welch提出的MAC方法不同，后者则是从流体质点入手。相对于MAC方法，VOF法可以处理自由面重入等强非线性现象，所需计算时间更短、存储量更少，但在处理网格单元中体积比函数，的变化时，稍显繁琐，而且有一定的人为因素。VOF法同MAC法一样，以压力P和速度“，V作为独立原始变量，边界条件易处理，为计算程序的编制提供了很大方便，对于研究多相流体交界面的运动变化有着非常大的吸引力。在VOF方法中，所有流体满足同一组动量方程，在整个计算区域上跟踪每一种流体在每个计算单元

35、中的体积分数，根据各个时刻流体在网格单元中所占体积分数来构造和追踪自由面。假设第q种流体在单元中的体积分数为Fq,第q种流体标量函数定义为fq,存在第q种流体空间点的fq值等于1,其它不被第q种流体占据点的fq值为0。在各网格单元中对fq值积分，并把此积分值除以单元的体积，得到单元的fq平均值，即网格单元中第q种流体所占据的体积分数Fq。若在某时刻网格单元中Fq=1，则说明单元中充满第q种流体；若Fq=0,则单元中不含第q种流体；当0Fq2Re=td/u图二维光滑圆柱体绕流的CD-Re关系曲线ReVV1,称为低雷诺数流动或蠕动流。几乎无流动分离，阻力以摩擦阻力为主，且与速度一次方成比例。lWR

36、eW500,有流动分离。当Re=100,圆柱后部有一对驻涡。当Re100时从圆柱后部交替释放出旋涡，组成卡门涡街。阻力由摩擦阻力和压差阻力两部分组成，且大致与速度的次方成比例。500WRe2X105,流动分离严重，大约从Re=104起，边界层甚至从圆柱的前部就开始分离，涡街破裂成为湍流,形成很宽的分离区。阻力以压差阻力为主，且与速度的二次方成比例，即CD几乎不随Re数变化。2X105Re3X106,CD与Re无关，称为自模区。本文中模型特征尺寸采用3m,速度为1m/s时雷诺数为X106，速度为3m/s时雷诺数约为X106。计算区域基本全部处于自模区。所以阻力系数大体不变。粘压阻力粘压阻力主要是

37、由于界层分离，其次由于边界层排挤厚度的存在使流线受挤压造成压力降低而产生。影响粘压阻力的因素有：物体后端形状，后端收缩缓和者粘压阻力较低；前端型线变化亦不宜太陡；物面型线应光顺；分离点之前的流态对分离点位置影响很大；层流中法向流速分布曲线较瘦，易发生分离，百度文库-好好学习，天天向上百度文库-好好学习，天天向上 - 分离点偏前，导致分离区域扩大，R增大；当层流转变为紊流时，分离点突然pv后移，将使阻力突然下降。在同一流态下；若界层未发生分离，则随速度的增加，C逐渐下降；在超过pv界分离的临界雷诺数之后，C与速度的增加或Re的增加基本无关，即C=常数。pp减小R的重要方法是改进船型，使之流线型化

38、，避免界层分离。p小结本章模拟了不同潜深的近水面潜器直航时的绕流情况，得到了它的兴波图形和阻力系数并与相关资料对比，发现计算结果正确的可能性很大。本章还介绍了本人使用成功的Gambit建模、选择控制区域及网格划分，边界条件的选择和初始条件的计算。得到了许多直观的数据图像，如速度云图，总压力云图等，为分析阻力提供了直观的依据。由本节的计算结果可以得到如下结论：1近水面航行器在水下航行时兴波阻力并不占绝对优势，即使是Fr=(v=2.7m/s)附近时阻力系数也没有明显表现出兴波阻力的变化趋势。由此可知完全在水下航行的近水面航行器阻力的主要成分依然是粘性阻力(摩擦阻力和粘压阻力)。当潜深大于3倍于潜器

39、直径时，兴波阻力的影响将很小。增加潜深可以减小近水面航行器的阻力系数。由压力云图可以看出圆柱体的前端与半圆球体接触的地方压力变化明显，这一点对压阻力的影响很大。由此可以推断更加光滑的流线型潜器的阻力情况可以有所改善。本文将继续研究型线方面的优化。对潜器型线的优化和对附体形状的优化型线优化流线型的计算模型及边界条件和计算条件在上一章发现圆柱形潜器首部压力变化明显的研究基础上，本文又建立了如下流线型的3D计算模型。该模型针对上述研究做出很大优化。模型及网格划分情况如图。图模型对称面附近的六面体网格该模型型线大体为流线型，总长仍然是3米，最宽处仍然是米。模型的控制区域仍然同前一章一样，边界条件方面来

40、流采用速度入口，出流采用自由出流（outflow）,其余均采用了壁面条件，计算介质为水，用来模拟深水直航时的情况。与之形成对比的前一章所提到的模型采用同样的边界条件和初始条件进行精确计算。计算条件大体与前一章一样，采用unsteady状况k-&模型并考虑了重力影响，用SIMPLE算法。壁面粗糙度设置仍然是默认的。其余参数由计算公式计算得到。计算结果本章将在这一节里对两种型线的潜器的计算结果进行对比，得到阻力系数更小的型线并分析原因，为设计更加快速的潜器提供依据。首先对比潜器的压力情况，分析两种型线的潜器表面压力有什么不同，从而得到阻力不同的原因。然后对比计算所得的阻力系数，确认结果与事实相符。

41、两种模型的压力云图如图。图v=lm/s时的动压力云图对比图v=s时的动压力云图对比图v=lm/s时的总压力云图对比图v=s时的总压力云图对比由两个模型的动压力云图对比可知，流线型的潜器并没有像圆柱形的潜器那样在首部有很明显的压力变化，它的压力分布比较均匀，压力变化较为缓和，因此粘压阻力减小了。由总压云图的对比可知，流线型的潜器只有在首部那一小部分压力比较大，尾部的一小部分压力比较小，其它地方总压力分布比较均匀，变化不大。首尾压差也不大。圆柱形的潜器首部较大面积里压力很大，而且压力较流线型潜器更趋向于指向后方。所以压力的水平分量比较大。过了半球体以后，圆柱体前半截的压力变小，这样就形成了更大的压

42、差，压阻力显然比较大些。所以总阻力也比流线型大一些。图v=s时垂直面上的动压力云图图是垂直面上的动压力图。由于重力的影响，云图并不是对称的。这张图反映了模型周围流体的压力变化情况。由图可知压力流畅内压力分布均匀且变化（形成压差）的区域不大。计算结果表明，当速度为lm/s时，圆柱形潜器的阻力系数为，而流线型潜器的阻力系数仅为，是圆柱形潜器阻力系数的%，还不及圆柱形潜器阻力系数的一半；当速度为2m/s时，圆柱形潜器的阻力系数为，流线型潜器的阻力系数仅为，是圆柱形潜器阻力系数的%；当速度为s时，圆柱形潜器的阻力系数为，而流线型潜器的阻力系数仅为，是圆柱形潜器阻力系数的%，同样不及圆柱形潜器阻力系数的

43、一半。由此可见型线对潜器阻力系数的影响之大。这并不奇怪，因为潜器的阻力系数主要是粘压阻力系数，摩擦阻力系数所占比例不是很大。另外由上述数据可以看出潜器速度越高，流线型潜器阻力系数减小越明显。不同附体对潜器阻力系数的影响计算模型、边界条件和计算条件为了使潜器能够在水中按照人们的意愿行驶，潜器要有完成控制方向等功能的附体。本文仅对两种不同形状的“舵”对阻力系数的影响做了简单的计算研究。一种是方形舵，另一种是较方形舵超过潜器直径较少的三角形舵。虽然两种舵形状不同，但面积相差不大。三角形舵沿着潜器轴线布置得更长一些。两种模型的网格划分情况与前面几个雷同，这里不再给出。其形状将在计算结果的压力云图中看到

44、。边界条件及计算条件和上一节中完全相同。计算结果依然按照前面所述的方法，首先对比压力云图。图。图v=1m/s时的动压力云图对比图v=s时的动压力云图对比百度文库-好好学习，天天向上百度文库-好好学习，天天向上 #- 图v=1m/s时的总压力云图对比图v=s时的总压力云图对比由压力云图可知，方形附体的潜器压力分布均匀，说明附体的压力与潜器主体的压力比较接近。而三角形附体附体部分的压力较主体小，主体部分压力图类似于无附体潜器。附体部分不难看出附体前端的压力大，后端的压力小，附体前后存在压力差。可见附体同样存在压阻力。但由于附体厚度很小，压阻力不大。图v=s时水平面上的总压力云图对比图是潜器外部流场的压力情况。计算结果表明，当速