一种模拟桥墩周期受力的数值方法技术

本发明专利技术是一种模拟桥墩周期受力的数值方法。根据不可压缩流的特点，在动量方程中通过加入两种不同形式的力，以提高一类以旋涡运动为主的流场的数值模拟精度。这两种形式的力分别是涡量在变化梯度方向的螺旋力和涡量在变化梯度方向的粘性耗散力。该方法使计算网格内的涡量在变化梯度方向的螺旋力的积分计算转化为计算网格边界上的上的力的通量计算；同时动量方程的源项保留涡量在变化梯度方向的粘性耗散力。

【技术实现步骤摘要】
一种模拟桥墩周期受力的数值方法
本专利技术涉及计算流体力学应用领域中的一种数值方法，具体是一种模拟桥墩周期受力的数值方法。
技术介绍
当河流经过大型桥墩物时，因为桥墩的迎风面和背风面的压力差的存在，在桥墩的后方会产生水流的漩涡。漩涡随着流水向下方继续流动，形成一系列脱落涡。这一系列脱落涡的流动也被称为卡门涡街，是一个充满了旋涡的旋流场，河流水流流场可以被视为是一个不可压缩流动。实际上，脱落涡会对桥墩周期地产生诱导作用力，严重的会使建筑物晃动，如果涡脱落的周期和桥墩建筑的自激频率相接近，建筑物的晃动会被放大，甚至发生毁坏。对于桥墩周期受力有多种研究手段。其中，计算机的数值模拟技术有着重要地位，是计算流体力学在该领域的一个扩展应用。计算流体力学综合了流体力学、应用数学、计算机科学，是一门应用性极强的学科。流体力学问题的数值模拟以其低成本、直观性强的优势，在流体流动的机理探索、工业产品设计等各个相关领域占据重要地位。桥墩后方脱落涡的数值模拟面临的最大的问题即是如何提高数值模拟的精度、降低误差，忠实地表现桥墩后方脱落涡流动的特性。影响流体力学问题的数值模拟的精度的重要因素之一是：当使用数值方法求解流体控制方程，即欧拉(Euler)方程或者纳维尔-斯托克斯(Navier-Stokes)方程时，会产生数值耗散(numericaldiffusion)，造成数值解的误差。例如数值方法中对控制方程的对流项的空间离散方法(如中心差分、迎风差分)、时间离散方法(如显示时间积分、隐式时间积分)、湍流模型(如双方程模型、大涡模拟)的使用，以及计算网格正交性都会产生不同程度的数值耗散。此外，数值模拟中还经常要使用一种人工数值耗散(artificialdiffusion)技术，其目的是通过适当降低计算精度而获得稳定的数值解。数值耗散对流场的数值模拟结果的最明显的影响体现在对流动变量的不连续界面(discontinuity)的捕捉。流场中一种强不连续的界面的捕捉，例如激波(shock)的捕捉，即是依靠加入适量的人工耗散项，以避免数值解在流动变量梯度变化较大的地方出现数值振荡现象。数值耗散可以理解为是流场中的一种能量损失，这种能量损失在某种程度上使得数值模拟结果不能忠实的体现流体的流动特性，降低了计算精度。先进的数值方法应该是在保证获得稳定性的数值解的前提下，将数值耗散减至最小。激波是强间断界面，对其捕捉必须加入一定的人工数值耗散。因为激波前后存在熵增，即能量的损失，所以，通过加入适当的人工数值耗散捕捉激波具有合理的物理意义。但是，诸如桥墩后方脱落涡一类的流动，存在着流场中另一类流动不连续现象，即接触不连续(contactdiscontinuity)。旋涡的产生自然和其周围的流体产生一个不连续面。这个接触不连续面相对激波而言是弱不连续，跨过不连续界面，压力和法向速度是连续的。数值模拟中对于这种接触不连续的捕捉更加困难，因为数值方法中的数值耗散即使很小也会使弱不连续界面变得模糊，降低数值解对流场的预测精度，这也是诸如桥墩后方脱落涡一类的旋流场的数值模拟技术成为CFD领域的重大挑战的原因。为了提高桥墩后方脱落涡流场的数值模拟的精度，一种方法是加密计算网格，在更加细小的空间尺度内求解流体控制方程。加密计算网格首先会使计算量加大，增加计算成本。此外，数值计算的误差随着计算网格的增加会不断积累，在一定程度上造成相反的效果。另一种方法是在流场中采用物理模型来增加流场中描述旋流流动的变量—涡量(vorticity)的强度。例如在流场中加入点涡模型，可以人为地增加涡量；或者在流场局部直接求解涡量方程，以减小涡量的输运过程中的耗散。但是，这些方法在应用上仍受到一定限制。点涡模型是在预先明确旋涡发生位置的前提下才能使用，仅适合一些简单的流动现象。除了二维不可压缩正压流场，涡量方程比与欲求解的Euler、Navier-Stokes方程更为复杂。二十世纪初期，美国科学家JohnSteinhoff提出了一种提高不可压缩旋流场的求解精度的数值方法，涡量限制法(VorticityConfinement)。该方法的原理是依靠在流场中加入限定的涡量以抵消数值耗散来模拟旋流场的流动状态。具体表现形式是在流体控制方程的动量方程中的源项位置，加入一个涡量形式的体积力项，从数值耗散中将涡量减去，克服数值耗散造成的旋涡场的接触不连续界面的模糊，从而更精确地捕捉旋涡结构，实现提高旋涡场的计算精度的目的。原始的涡量限制法是公知的，这里不再叙述。尽管该方法在捕捉不可压缩流场中的接触不连续的方面已经取得了明显的改进效果，但存在以下缺陷：1.对加入的涡量调整完全靠常系数；2.加入的涡量的空间离散精度是被限定的，无法进一步提高；3.源项对动量方程的数值解的收敛的稳定性影响是不确定的。为了更加精确地模拟桥墩后方脱落涡流动—不可压缩无粘流的旋涡运动，需要进一步改进对于流场中接触不连续的捕捉机制。一种途径是根据不可压缩流的特点，改进在流场中通过加入涡量来抵消数值耗散的内在工作机理，通过保持涡量的精度，更精确地用模拟流场中的旋涡运动，形成一种新的，针对桥墩后方的涡脱落长生的周期受力的数值模拟技术。该技术可以使涡量的空间离散具有高阶精度的格式，同时该可以利用源项增进收敛进程的稳定性。
技术实现思路
本专利技术涉及计算流体力学的应用领域中一种模拟桥墩周期受力(不可压缩无粘旋流流动)的数值方法，具体是根据不可压缩流的特点，在流场中通过加入两种不同形式的涡量力来抵消数值耗散的数值方法。该方法可以使加入的涡量的空间离散具有高阶精度的格式，同时还可以利用源项增进数值解的收敛进程的稳定性。通过保持涡量的精度，更精确地用模拟流场中的旋涡运动，是一种新的不可压缩流的旋涡运动的数值模拟技术。首先写出不可压缩、无粘流流动的控制方程，包括连续方程和动量方程，分别为式中，V是速度矢量，包含直角坐标系下的三维i，j，k分量u，v，w；ρ、p、t分别是密度、压力和时间。式中算子表示为符号·代表内积计算。动量方程(2)的右边是体积力项形式的源项按照涡量限制法中的定义式中，代表涡量，在直角坐标系下有按照涡量的定义，式中符号×代表叉乘运算，代表涡量梯度变化的方向，即其中，φ是涡量的模；是涡量的模的梯度的模，即动量方程(2)中的源项需要公式(3)中的乘以一个比例因子ε，则其物理意义是指向涡量中心的力。原始的涡量限制法中给出ε为常数，为0.01-0.05，其作用是用来调整涡量限制的大小。该方法形式简单，不需要加密计算网格，在不可压缩旋流场的数值模拟中得到广泛应用。对该方法在随后的改进主要集中在参数ε的表达。如公式(2)表示的，涡量限制法中在不可压缩流动的动量方程的等号右边加入了一个体积力项，它代表涡量在其变化梯度相反的方向的受力，如公式(7)所示。实际上，如果将公式(4)带入公式示(7)并运用算子的运算规律可得在上式中再次运用算子的运算规律可得式中是拉普拉斯算子，定义为又根据公式(1)不可压缩流的连续方程，将公式(9)代入公式(8)，可得对于不可压缩流，公式(10)与公式(7)完全等价。现将公式(10)中的重新写成如下形式，即其中，因为是螺旋度的定以，所以被称作涡量在变化梯度方向的螺旋力；因为拉普拉斯算子的耗散特性，被称作在涡量在变化本文档来自技高网...

【技术保护点】
一种模拟桥墩周期受力的数值方法，其特征在于在浅水流动的动量方程中加入两种不同形式的力，分别是涡量在变化梯度方向的螺旋力

【技术特征摘要】
1.一种模拟桥墩周期受力的数值方法，其特征在于在浅水流动的动量方程中加入两种不同形式的力，分别是涡量在变化梯度方向的螺旋力和涡量在变化梯度方向的粘性耗散力其形式分别为

【专利技术属性】
技术研发人员：路明，
申请(专利权)人：成都金景盛风科技有限公司，
类型：发明
国别省市：四川,51