共查询到20条相似文献,搜索用时 156 毫秒
1.
首先,简单介绍了基于粘接元的无重叠区域分裂方法.这种方法利用变分原理,非常适合有限元近似.然后,着重讨论了这种区域分裂方法在求解不可压Navier-Stokes方程中的应用,具体包括等价变分公式的建立、通过算子分裂的时间离散、区域分裂情形下广义Stokes问题的共轭梯度迭代求解方法、空间的有限元离散.最后,以数值实验结果验证了这种区域分裂方法应用于不可压Navier-Stokes方程求解时的可靠性. 相似文献
2.
3.
隐式格式求解拟压缩性非定常不可压Navier-Stokes方程 总被引:1,自引:0,他引:1
采用Rogers发展的双时间步拟压缩方法,数值求解不可压非定常问题.数值通量分别采用三阶精度Roe格式和二阶精度Harten-Yee的TVD格式离散.为了加快收敛,提高求解效率,试验了几种隐式格式(ADI-LU,LGS,LU-SGS).针对经典的低雷诺数(Re=200)圆柱绕流问题,比较了不同隐式方法的计算结果和求解效率,以及两种数值离散格式计算结果的异同.最后采用Roe格式数值求解了两种典型的低速非定常流动问题:绕转动圆柱(ω=1)低雷诺数流动;NACA0015翼型等速拉起数值模拟. 相似文献
4.
多孔介质中高Péclet数和大黏性比下混溶流体的流动和扩散广泛存在于二氧化碳驱油、化工生产等工业过程中.用数值方法对该问题进行研究时,关键在于如何正确描述高Péclet数和大黏性比下多孔介质内流体的行为.为此,提出了一种基于多松弛模型和格子动理模型的耦合格子Boltzmann模型.通过Chapman-Enskog分析,证明该模型能有效求解不可压Navier-Stokes方程和对流扩散方程.数值结果表明,该模型不仅具有二阶精度和良好的稳健性,而且对于高Péclet数和大黏性比的问题具有良好的数值稳定性,为模拟此类问题提供了有效工具. 相似文献
5.
提出了一种改进的基于相场理论的两相流格子Boltzmann模型.通过引入一种新的更加简化的外力项分布函数,使得此模型克服了前人工作中界面力尺度与理论分析不一致的问题,并且通过Chapman-Enskog多尺度分析表明,所提出的模型能够准确恢复到追踪界面的Cahn-Hilliard方程和不可压的Navier-Stokes方程,并且宏观速度的计算更为简化.利用所提模型对几个经典两相流问题,包括静态液滴测试、液滴合并问题、亚稳态分解以及瑞利-泰勒不稳定性进行了数值模拟,发现本模型可以获得量级为10-9极小的虚假速度,并且这些算例获取的数值解与解析解或已有的文献结果相吻合,从而验证了模型的准确性和可行性.最后,利用所发展的两相流格子Boltzmann模型研究了随机扰动的瑞利-泰勒不稳定性问题,并着重分析了雷诺数对流体相界面的影响.发现对于高雷诺数情形,在演化前期,流体界面出现一排“蘑菇”形状,而在演化后期,流体界面呈现十分复杂的混沌拓扑结构.不同于高雷诺数情形,低雷诺数时流体界面变得相对光滑,在演化后期未观察到混沌拓扑结构. 相似文献
6.
谱元方法求解正方形封闭空腔内的自然对流换热 总被引:6,自引:0,他引:6
提出谱元方法计算正方形截面封闭空腔内的自然对流问题,具体求解了原始变量速度和压力的不可压Navier-Stokes方程和温度方程,所有的求为量均采用Chebyshev谱逼近,Navier-Stokes方程和温度方程的时间离散采用时间分裂法,非线性项用4阶Runge-Kutta法,扩散项用Crank-Nicolson半隐方法,获得了与文献发表的基准解较一致的计算结果。 相似文献
7.
基于两重网格离散方法,提出三种求解大雷诺数定常Navier-Stokes方程的两水平亚格子模型稳定化有限元算法.其基本思想是首先在一粗网格上求解带有亚格子模型稳定项的Navier-Stokes方程,然后在细网格上分别用三种不同的校正格式求解一个亚格子模型稳定化的线性问题,以校正粗网格解.通过适当的稳定化参数和粗细网格尺寸的选取,这些算法能取得最优渐近收敛阶的有限元解.最后,用数值模拟验证三种算法的有效性. 相似文献
8.
VOF法模拟剪切流动下液滴的变形和断裂运动 总被引:1,自引:0,他引:1
本文对剪切作用下悬浮液滴在另一种不相融的液体中的变形和断裂过程进行了数值模拟.采用VOF(Volume ofFluid)法中的三维PLIC(Piecewise Linear Interface Calculation)算法实现界面的重构和输运,交错网格下投影法离散表面张力为源项的不可压缩Navier-Stokes方程.... 相似文献
9.
10.
求解Navier-Stokes方程组的组合紧致迎风格式 总被引:1,自引:0,他引:1
给出一种新的至少有四阶精度的组合紧致迎风(CCU)格式,该格式有较高的逼近解率,利用该组合迎风格式,提出一种新的适合于在交错网格系统下求解Navier-Stokes方程组的高精度紧致差分投影算法.用组合紧致迎风格式离散对流项,粘性项、压力梯度项以及压力Poisson方程均采用四阶对称型紧致差分格式逼近,算法的整体精度不低于四阶.通过对Taylor涡列、对流占优扩散问题和双周期双剪切层流动问题的计算表明,该算法适合于对复杂流体流动问题的数值模拟. 相似文献
11.
提出二维可压缩流体力学问题的拉格朗日有限点方法,将求解区域离散为适当的点集.在每个时间步,每个离散点与其周围适当的五个邻点组成一个基本计算单元.在每个计算单元上,利用有限点方法中的典型微分算子的五点近似公式直接离散流体力学方程中的微分算子,并在每个方程中加上一个人为拉普拉斯粘性项,达到稳定格式的目的.给出时间步长的自动选取算法.数值算例结果验证了算法的有效性,初步展示了其计算大变形流体问题的良好发展潜力. 相似文献
12.
本文继Zheng-huan Teng[1]中介绍的解Navier-stokes方程的椭圆涡团法,研究了一种新的变形涡团法,用以模拟不可压粘性流体绕圆柱的不定常流动。圆柱在静止流体中突然起动并做匀速直线运动。对整个流动区域构造完全Navier-Stokes方程的解并不容易,近十年出现很多数值研究,本文对算法有所推进。把本文的方法称作变形涡团法是因为圆柱边界附近的流体中用椭圆涡团,远离边界时用圆形涡团。计算圆柱绕流比平板绕流在满足附着条件上更为困难,本文分析了怎样在圆柱边界上给出适当的附着条件的数值方法。在算例中雷诺数分别取200、550、3000,得到了不定常边界层分离,二次涡等复杂的物理现象,这些数值结果与近年实验结果[2]是一致的。 相似文献
13.
14.
本文通过破开算子方法,把二维输运问题的控制方程破开为对流问题和扩散问题。在任意四边形网格的离散下,用特征线法解对流问题,并采用伽辽金加权余量法,从而有效地减少插值所引起的数值阻尼,提高计算精度。用有限单元法和迭代计算格式解扩散问题。由于采用了辛普生积分公式,在每个时间步长都不需要求逆矩律,节省了计算时间。算例表明,本文数值模拟结果与精确的理论解吻合较好。 相似文献
15.
16.
本文提出了一种数值求解对流扩散方程的分步杂交方法。在不规则的三角形网格上,采用迎风离散格式或改型特征线方法处理对流算子;采用集中质量的有限元方法处理扩散算子。详细分析了这种算法的稳定性同题,在数学上严格证明了在满足①Δt≤min((2d)/v,(d2)/(3K)),其中d是三角形网格中最短垂线的长度,V和K分别为流场中的最大速度和扩散系数。②所有三角形的内角θ≤π/2的条件下,整个计算格式是L∞稳定的,从而保证了在海洋环境和水质的数值模拟中海水的盐度、污染物的浓度和核电站冷却水系统中的超温不会出现负值。应用非线性的对流扩散方程对此方法的精度和收敛性进行了检验。通过数值解与精确解的比较,表明本方法的数值耗散很小,用改型特征线方法处理对流算子较迎风离散格式有更高的精度;两种处理对流算子的方法都没有伪振荡现象发生。本方法由于具有算法简单、L∞稳定、计算网格灵活等优点,可推广使用于实际的海洋环境(潮波、海流、海洋污染)、港口和海湾的数值模拟以及不可压粘性流和对流传热同题的数值计算。 相似文献
17.
18.
用DES数值模拟分离绕流中的旋涡运动 总被引:1,自引:0,他引:1
脱体涡模拟(Detached-Eddy Simulation,DES)是近年来出现的一种结合雷诺平均方法和大涡数值模拟两者优点的湍流模拟方法.采用基于Spalart-Allmaras方程模型的DES方法,数值求解Navier-Stokes方程,模拟绕流发生分离后的旋涡运动.其中空间区域离散采用有限体积法,方程空间项和时间项的数值离散分别采用Jameson中心格式和双时间步长推进方法.通过模拟圆柱绕流以及翼型失速绕流,观察到了与物理现象一致的旋涡结构,得到与实验数据相吻合的计算结果. 相似文献
19.
一、引言 边界元计算方法是七十年代迅速发展起来的一种数值计算方法,其主要优点是:将求解区域微分方程的问题转化成求解边界积分方程的问题,因而一般都把物理问题降了一维求解,使该方法计算效率和求解精度都较高.但它用于时关问题和非线性问题时,积分方程中还含有物理量的区域积分项,该方法的优点几乎全部消失.另外,在边界积分方程离散后,代数方程的系数矩阵为满阵.如果边界单元划分很多,其效率不如具 相似文献
20.
流体系统在重力作用下,流体压力可与重力相平衡,达到静力平衡状态.常规的有限差分和有限体积方法无法在离散尺度上保持这种平衡态,常常产生虚假速度等非物理现象.通过将重力源项的体积分改写为网格界面值的加权平均,实现离散条件下压力与重力的精细平衡,得到一种Navier-Stokes (NS)方程的精细平衡气体动理学格式.该格式可以在机器精度上精确地保持等温条件下的静力平衡态以及捕捉平衡态附近的小扰动传播.同时,该格式还能求解自然界中更常见的非等温平衡态.此类平衡态除了要求静力平衡还要求热流平衡.利用提出的格式求解NS方程,可以得到静止(机器精度下的零速度)的非等温平衡态,并且密度和温度分布都具有二阶精度.通过多个算例验证格式的有效性,表明该格式可更好地模拟重力场下的温度、密度的分层流动. 相似文献