飞行器操纵面嗡鸣的非结构网格并行计算方法

为了数值模拟飞行器操纵面的嗡鸣现象,在集群计算机MPI并行计算环境下建立基于三维非定常欧拉方程耦合结构运动方程的嗡鸣计算方法.气动流场求解采用基于非结构网格的中心有限体积法进行空间离散,时间推进采用双时间方法,结构运动方程采用Adams预估校正方法求解.针对翼面与操纵面缝隙间存在的网格运动问题,在非结构网格系统上采用Delaunay图映射方法实现网格的运动变形.最后,使用飞行器操纵面标准嗡鸣计算模型对计算方法进行验证,结果表明:所建立的并行计算方法正确,程序具有很好的计算效率,能够对飞行器操纵面嗡鸣进行高效的数值分析.     

新的多重网格算法在压气机流场计算中的应用

根据在流场数值计算中使用多重网格方法可以加速收敛的原理,在叶栅流场控制方程使用有限体积差分法进行离散的前提下,提出了一种新的4重网格方法来加速求解压气机叶栅可压缩流动,该方法可以与显式时间推进法有效配合使用。通过对NASA67超音速轴流压气机叶栅流场的数值计算表明,该方法明显加快了流场计算的收敛速度,并且计算结果与试验结果吻合较好。     

薛定谔方程向后欧拉全离散两网格混合有限元方法

针对二维依赖于时间的线性薛定谔方程,在空间方向采用混合有限元方法,时间方向利用向后欧拉方法,得到一种全离散混合有限元格式.为了将薛定谔方程耦合的实部和虚部解耦,提出了一种全离散混合有限元的两网格算法,将方程在细网格上的求解问题,简化为在一个相对更粗的网格上求解原问题以及在细网格上求解两个泊松方程,从而减小计算工作量,节省计算时间.数值实验结果验证了两网格混合有限元方法的高效性.     

一类导热型微分方程的数值求解方法及其比较

针对一维、稳态且含内热源的导热型微分方程,分别采用区域离散方法A和B划分网格,且每种网格划分方法中又分别采用Taylor展开法和控制容积平衡法进行数值求解。计算结果表明,虽然Taylor展开法中右端点离散方程采用二阶截差展开时数值求解精度较高,但离散方程推导过程复杂,尤其在计算节点间距不等的情况下更是如此;控制容积平衡法物理意义明晰,数值求解精度高,无论对于区域离散方法A或B,其离散方程推导过程均简单明了。     

二阶时域波动方程的无网格方法求解

将径向基函数配点型无网格方法引入二阶时域波动方程的求解中,方程的空间导数采用径向基函数逼近,时间导数采用Crank-Nicolson方法离散,对应的边界条件直接施加在离散的边界数据点上.采用该方法对二维非规则求解域内的波传播问题进行了数值计算,并与有限元计算结果进行了对比分析.结果表明:基于径向基函数配点的无网格方法不但形式简单、易于实施,而且能够有效解决复杂求解域高维的波动问题.     

180°方形弯管中湍流流动的大涡模拟

把大涡模拟应用于复杂几何形状的流域是一个重要研究课题。对１８０°方形弯管中的湍流流动作大涡模拟,以网格节点为中心的交错网格作为离散网格,用有限体积方法和中心差分进行空间离散,用Ａｄａｍｓ－Ｂａｓｈｆｏｒｔｈ格式进行时间离散,用简单的Ｓｍａｇｏｒｉｎｓｋｙ模型来作为次网格湍流模型。用多重网格方法求解压力方程以提高计算效率。数值模拟结果与实验结果取得了很好的一致。     

预处理方法在含运动边界问题中的应用

通过把预处理方法引入ALE(arbitrary Lagrangian Eulerian)控制方程的求解中,使得对含有动边界的从低速到高速的广泛流动问题可采用统一的数值算法.采用MUSCL格式在非结构网格上进行空间离散;采用隐式双时间推进,以提高计算效率并保持足够的时间精度.通过几个算例验证了该算法的计算能力.     

带端板地效应翼性能的数值研究

通过求解不可压缩流体ＲＡＮＳ方程,数值模拟带端板三维地效应翼的性能及周围流场。数值方法引进了Ｃｈｏｒｉｎ的人工可压缩性概念,应用近似因式分解技术同时求解速度和压力场,动量方程对流项用二阶迎风差分格式离散,其余空间导数项均采用四阶精度的中心差分格式离散,时间离散采用欧拉隐式格式,计算在非交错网格上进行,为了避免压力场的振荡,在连续方程中隐式地加入了压力的四阶数值耗散项,湍流计算采用了Ｂａｌｄｗｉｎ－     

基于非结构化网格有限体积的LBM

提出了一种基于非结构化网格有限体积的LBM.采用Cell-vertex有限体积法离散控制方程.该方法在时间上采用伪、实二时间步,其中伪时间步采用向前差分,实时间步采用二阶向后差分方法;空间上采用edge-based通量计算方法,采用高阶TVD格式计算控制体边界通量.离散后的控制方程采用隐式迭代,控制变量采用五层二阶Runge-Kutta方法求解.二维同心圆环内圆柱间Couette流与顶盖驱动方腔流的数值结果显示该方法为一种有效求解不规则边界流体动力特性的实用工具.     

基于二维非结构网格的全隐式,迎风格式欧拉方程研究

对二维非结构网格进行了网格重构,使用格点格式有限体积法离散流场方程,采用Roe矢通量分裂方法计算通量,时间推进采用Lu-SGS方法.通过对NACA0012翼型和RAE2822翼型的流场模拟,验证了本算法的可行性和高效性.     

多重网格法求解原始变量形式的Navier-Stokes方程

采用多重网格法 ,求解原始变量形式的Navier Stokes方程 .控制方程在交错网格上离散 ,并采用SIMPLE算法计算 .为加速收敛 ,结合使用适合非线性方程的全近似格式的多重网格法FAS .文章计算了雷诺数Re =10 0 ,10 0 0时不可压粘性流体绕圆柱的流动 .计算结果与实验结果和国外的计算结果基本吻合 ,说明了该方法的准确性和有效性     

Cahn-Hilliard方程的时间双层网格混合有限元方法

针对数值求解Cahn-Hilliard方程时非线性项引起的时间耗时问题,提出了时间双层网格混合有限元方法.首先,在时间粗网格上,通过非线性牛顿迭代方法求解非线性混合有限元系统,其中空间离散采用混合有限元方法,时间离散采用隐式欧拉格式;其次,基于初始迭代数值解和拉格朗日插值公式,在时间细网格上求解线性混合有限元系统;最后,分析了该方法的稳定性和误差估计,并通过数值算例进行验证.结果表明,与传统的混合有限元方法相比,该方法可以节省计算时间.     

一种无积分任意四边形非结构化网格节点间断Galerkin方法

节点间断Galerkin方法是近年来得到迅速发展的高精度数值方法,可以采用任意多边形网格对平面求解域进行离散.针对任意四边形非结构化网格,传统的节点间断Galerkin方法采用数值积分对离散方程进行计算,需要较大的计算量与存储空间.为了提高任意四边形非结构化网格上节点间断Galerkin方法的计算效率,提出了一种新的无积分格式实现方法,即将积分节点与插值节点定义为同一节点集,并利用节点基函数的插值性质,推导出每个单元内控制方程的无积分离散格式.通过在任意四边形非结构化网格中对二维对流方程进行数值求解,验证了新提出的无积分方法的准确性和计算效率.结果表明,无积分方法与传统数值积分方法计算误差和收敛精度基本相同,而其计算效率提高1倍以上.     

管式反应器中交混流场的数值模拟

对管式反应器中的交叉混合流场进行了三维数值模拟,采用分块结构化网格技术以及网络拼接技术实现了适用于复杂形体的分块结构化网格的计算程序,采用修正的ｋ－ε双方程湍流模型描述反应器的湍流流动,用类似ＳＩＭＰＬＥ法对控制方程进行求解,经冷态模拟试验检验,结果令人满意     

弯曲圆管中矿浆湍流场的三维数值模拟及分析

采用三维、定常、不可压缩,时均Ｎ－Ｓ方程和混合长湍流模型,建立了弯曲圆管中矿浆流动的控制方程和模型,使用有限元法对计算区域进行了离散,并对流动方程进行了求解,为减小用有限元法计算三维湍流场的计算量而采用了加罚有限元法。     

非结构化网格上N-S方程的离散

采用有限容积法离散控制方程,实现了二维非结构化同位网格上的SIMPLE算法,并编写了四边形网格求解不可压缩N-S方程组的程序.分别采用三角形网格与四边形网格,对计算传热学中有基准解的顶盖驱动空腔流问题进行了模拟.结果表明:所提算法既可用于三角形网格,又可用于四边形网格；三角形网格与四边形网格的模拟解均与基准解吻合,验证...     

基于高阶紧致格式的二维不可压NS方程求解

本文介绍了一种基于原始变量的用于求解二维非定常不可压Navier-Stokes方程的高阶紧致格式。这种紧致格式最初是用于计算声学（CAA）的高精度格式,相对于传统的紧致格式,使用该格式的优点在于减少计算量的同时降低了边界模板的处理难度。这种方法建立在非交错网格上,空间离散具有六阶精度。压力Poisson方程基于九基点模板的四阶紧致格式进行离散,超松弛迭代进行求解。时间推进上采用四阶Runge-Kutta方法。为验证该方法的精度和有效性,利用该格式计算了一个具有解析解的问题,以及二维非定常情况下的方腔驱动流动问题,并且和传统的紧致格式进行了计算时间的对比。     

模拟圆柱涡致振动问题的重叠网格法

使用重叠网格法求解二维不可压N—S方程,数值模拟了雷诺数在90—150之间的圆柱涡致振动。N—S方程对流项采用三阶迎风格式,扩散项使用四阶中心紧致格式离散,圆柱计算网格随圆柱一起运动,采用运动的非惯性参照系下的控制方程,振动控制方程采用四阶龙枷库塔法求解。计算成功地捕捉到了“拍”现象和“锁定”现象,并给出在不同折合来流速度下圆柱振幅的实验值和数值模拟结果的对比,计算结果与实验值基本吻合。     

导体电磁散射特性无网格方法研究

传统数值方法的展开依赖于基函数网格生成,当离散网格不满足其构成条件时,算法的计算精度和效率会降低。针对这一问题,该文提出了不受离散网格质量限制的无网格方法。该方法采用无规律、自由生成的任意网格,基于点离散基础对目标的电磁特性进行数值分析。给出了无网格法在求解导体表面电磁散射时的矩阵求解过程,包括方程离散过程、矩阵元素计算过程及求解表面电场积分方程时对奇异积分的处理过程。几种典型算例的分析结果表明,无网格方法结合该文预条件提高了近相互作用的准确性。     

多段翼型无网格算法及其布点技术

用显式无网格算法实现了多段翼型的数值模拟。与传统的网格算法不同,求解区域用“点云”离散代替通常的网格划分。基于当地点云离散结构,用二次极小曲面逼近计算空间导数。在研究该算法的基础上,给出了Euler方程无网格离散形式,运用Runge-Kutta显式时间推进格式推进求解。此外,还描述了一种区域离散布点方法,研究了点云生成的选点准则,并成功地数值模拟了复杂多段翼型的绕流。