时差法超声波流量计的数值模拟。在现场实际应用中， 时差法超声波流量计测量流量是目前广泛应用的一种方法，它具有安装简单、使用方便、测量范围宽等特点。时差法超声波流量计应用于大口径管道水流量的测量中具有明显的优越性。近几年来，随着硬件和数字信号处理技术的改进使流量的测量精度得到了不断提高。
但是，超声波流量计的应用对管道的条件有一定的要求。一般，换能器应安装在上游直管段为10D（D为管径）、下游直管段为5D的位置。因此，选择良好的安装位置和安装条件的管道才能获得满意的应用效果。对于同一个超声波流量计在应用于不同的安装条件时，流量的测量结果应使用不同的修正系数加以修正。通常，修正系数的确定是在实验室中测试得出的，但这只能针对有限的安装条件，而且实验室的建设和测试的过程既费时又费钱，效率很低。
对不同雷诺系数的流体和在不同超声波波束直径的测量进行了数值模拟，计算出流量计的流速修正系数，并通过分析计算的方法计算出修正系数，两者进行比较，比较结果具有很好的一致性。
对于所有的流体流动，计算机数值模拟都是解质量和动量守恒方程。质量守恒方程又称为连续性方程。该方程是质量守恒方程的一般形式，它适用于可压流动和不可压流动。源项Sm是从分散的二级相中加入到连续相的质量，源项也可以是如何自定义源项。
修定系数
由于流体的流速沿管道直径的不均匀分布，致使时差法超声波流量计测得的流速并不等于实际的流速，因此，测量出的结果需乘上一个修正系数k，即
公式中:为超声波流量计测得的平均流速，为实际的流体流速，其中
k=k[u（x，y，z），v（x，y，z），w（x，y，z）] （6）
这里的k系数是基于CFD数值模拟的预测模型模拟出来的，因此需根据超声波流量计的测量原理建立合适的模型，测量的原理图如图1所示。
从图1可以看出，结合时差法超声波流量计的流速计算方法，可以得出测得的平均流速为:
式中:L为超声波传播的距离，α为超声波传播方向与管道轴线的夹角，t1和t2分别是超声波逆流和顺流方向的传播时间。
从图1中分析，由于流体流速的两个分量u和v对超声波传播方向上都有作用，仅垂直于传播方向的流速分量没有作用。因此，要得到k（u，v）则需了解管道内流体的流动状态，这个可以利用CFD软件根据流体在管道入口的初始状态和管道的安装状况模拟得出。
分析计算和CFD建模
当流体在紊流状态时，流速曲线方程为:
式中:r为截面中的任意点到中心轴的距离，umax为中心轴上的流体流速，即截面上的最大流速，p是由雷诺系数决定的参数。
p=0.25-0.023lgRe （9）
当流体在层流状态时，流速曲线方程为:
从图1中可以得知，由超声波测得的流体流速与管径d、夹角α和超声波传播波束的直径b有关，因此利用上述的流速曲线方程和超声波传播的路径可以得出流体的平均流速:
为了计算出流速的修正系数，还要得出流体的实际流速。由紊流理论可知，流体截面上各点的流速不是均匀分布的，从Prandtl定律中的描述可以知道，流体的实际流速可以由雷诺系数得出，即:
式中:utrue为管道内流体的平均流速，d为管道内径，ν为运动粘度（ν=8×10-7m2/s，水温为30℃）。从上面几式中可以得出k系数:
通过改变超声波传播范围的半径b与管道半径的比例，可以得出在不同规格超声波换能器下测量的流速修正系数k。这种分析计算模型可以适用于层流和紊流两种状态下。
这里利用计算流体力学软件Fluent进行管道流体模拟运算，Fluent软件是个工程运用的CFD软件，针对每一种流动的物理问题的特点，采用适合于它的数值解法在计算速度、稳定性和精度方面达到最佳，可以计算流场、传热和化学反应。其思想实际上就是做很多模块，这样只要判断是哪一种流场和边界就可以拿已有的模型来计算[3]。
Fluent软件能推出多种优化的物理模型，如定常和非定常流动、层流（包括各种非牛顿流模型）、紊流（包括最先进的紊流模型）、不可压缩和可压缩流动、传热、化学反应等等。对每一种物理问题的流动特点，有适合它的数值解法，用户可对显式或隐式差分格式进行选择，以期在计算速度、稳定性和精度等方面达到最佳。Fluent将不同领域的计算软件组合起来，成为CFD计算机软件群，软件之间可以方便地进行数值交换，并采用统一的前、后处理工具，这就省却了科研工作者在计算方法、编程、前后处理等方面投入的重复、低效的劳动，而可以将主要精力和智慧用于物理问题本身的探索上。在Fluent软件中，采用Gambit的专用前处理软件，使网格可以有多种形状。对二维流动，可以生成三角形和矩形网格；对于三维流动，则可生成四面体、六面体、三角柱和金字塔等网格；结合具体计算，还可生成混合网格，其自适应功能，能对网格进行细分或粗化，或生成不连续网格、可变网格和滑动网格。
网格划分的好坏将直接影响模拟的精度和收敛性，本文管道截面的网格划分采用的是矩形网格，外层采用的是同心圆的排列方式，里层采用的是交叉排列方式，式样如图2所示，这种混合网格排列的方式有利于管道边界层的模拟效果。
对于管道流速修正系数k（u，v）可以通过Fluent软件模拟的流速数据计算出来，根据网格划分的体积权重来计算超声波传播路径上的平均流速
式中:ui和vi可以从Fluent输出的数据中得到，Vi为体积，如图1中坐标系所标的u和v两个速度分量。根据Fluent的数值模拟可以得出流速修正系数k（u，v）:
式中:是流体的实际流速，是流体在超声波传播路径上的平均流速，α是超声波传播路径和管道轴线的夹角。
试验结果分析
设定的边界条件为:管径D为100mm，长度L为50D，以保证在数值模拟中流体能够得到充分流动，入口流体流速设定为1m/s。
采用三种方式通过分析计算和数值模拟来计算流速修正系数K值:
（1）当流体雷诺系数Re=100，流体流动状态在层流的情况下，通过改变超声波传播波束的直径与管径的比率b/d，计算出流速修正系数k值，结果如图3所示。
（2）当流体雷诺系数Re=10000，流体流动状态在紊流的情况下，同样通过改变b/d值，计算出流速修正系数k值，结果如图4所示。
（3）在超声波传播波束的直径与管径的比率b/d=0.5的情况下，通过改变流体雷诺系数Re的大小来计算流速修正系数k值。以上就是新疆流量计厂家介绍的新疆涡轮流量计相关知识，希望大家能记住，在日后的使用中有所帮助，如果还有其它问题欢迎来电咨询，我们将期待您的光临！