矢性流函数涡量法充型过程数值模拟

建立了用矢性流函数－涡量法和ｋ－ε双方程湍流模型模拟铸造充型过程的数学模型,并模拟了铸造充型过程。对模拟结果的验证表明,计算结果与实测结果基本一致,证明矢性流函数－涡量法和ｋ－ε双方程湍流模型可用于模拟铸造充型过程三维非稳态湍流流动和传热过程。     

平壁附近钝体绕流的数值模拟

采用有限差分法,模拟粘性流绕平壁面附近的圆柱的流动。流动方程采用涡量-流函数形式。涡量输运方程采用ADI方法求解,流函数方程采用SOR方法求解。计算网格采用椭圆型贴体曲线网格技术数值生成。考察了缝隙比（e/D)对圆柱尾流的影响,得到了与实验基本符合的结果。     

非定常Stokes方程混合有限元方法

本文给出二维非定常Stokes方程的流函数-涡度表达式, 采用混合有限元方法分别讨论流函数方程和涡度方程,得到流函数、涡度及流速场的最优阶L2误差估计。     

圆筒燃烧室旋流流场的研究

本文介绍圆筒燃烧室内双股同心旋转射流物理模型实验与数值模拟的研究结果。数值模拟按著者最近提出的湍流涡量输运方程,采用ψ-ω方法及K-ε湍流双方程模型,Rodi提出的壁函数,以及Spalding教授等提出的计算湍流回流的方法。本文工作表明,用涡量—流函数法及K—ε双方程模型计算强旋流及回流是可行的。     

二维非定常有涡流动的数值模拟方法

为有效地模拟二维有分离现象的粘性流动,发展了求解涡量流函数方程的数值方法。对涡量输运方程的时间项运用四阶Runge-Kutta法离散,从而将方程分解为4个计算步分别求解。在每一计算步对方程进行Steger-Warming近似因式分解,从而使涡量的计算在空间的两个方向上分别进行。对流项采用Chakravaythy-Oscher总变差减少(TVD)格式离散。流函数的Poisson方程运用Tschebyscheff SLOR方法交替方向迭代求解。运用该方法对突然起动圆柱和高雷诺数时弯曲薄翼的非定常有涡流动进行了数值模拟,并将计算结果与其它方法的计算结果和试验结果进行了对比,结果表明本方法具有精度高、数值稳定性好和计算效率高的优点。     

流体力学有限元分析中的边界条件处理

阐述了流体力学有限元分析中应用流函数 -涡量法时典型边界条件的处理方法 ,并给出计算实例     

DC—EAES—LF熔池电弧电渣加热过程的三维数值模拟

建立了非规范正交曲线坐标系下ＤＣ－ＥＡＥＳ－ＬＦ炉电弧电渣加热过程三维数学模型，并进行了联立耦合求解：速度场计算使用了矢性流函数方程的涡量输运方程；湍流模型采用引入Ｒｅｃｈａｒｄｓｏｎ浮力修正贡的ｋ－ε双方程模型；温度场用非定常能量输运方程描述。现场测量的温度场与预报结果基本一致。     

一种紧致差分格式及其在新型纺纱研究上的应用

本文导出了一种与涡传输方程的ADI格式相协调的流函数方程的四阶精度五点紧致差分格式,并导出了与此相容的速度差分方程及三阶精度的壁涡公式,以及在非网格边界点上的壁涡公式与相邻内点的涡传输差分方程和流函数差分方程。本文成功地用此方法对某新型纺纱器罩壳内的流动作了数值计算。其结果与流场显示法的实验结果(照片)基本符合,对新型纺纱的理论研究以及罩壳的合理改型提供了可靠的依据。     

二维不可压N-S方程二次四边形单元有限元解

粘性不可压流体问题是众多工程中重要的力学问题.数值求解Navier-Stokes方程会遇到两大困难:非线性和不可压性.针对二维不可压Navier-Stokes方程的特点,建立了以流函数为求解变量的四阶微分控制方程,有效地避免了处理涡量边界的难题.采用8节点二次四边形单元,单元基函数为2次非线性高阶函数,建立了求解二维不可压N-S方程的有限元方程,并自主开发了二次四边形单元有限元程序.数值实验结果验证了该方法的精确性和可靠性.因此,该方法在计算流体力学中有较好的应用前景.     

圆柱绕流的边界层及其电磁控制

利用电介质溶液中圆柱体侧表面附近分布的电磁产生电磁力,作用于流体边界层,达到控制圆柱绕流的目的。在层流模式下(Re<150),采用极坐标系中NavierStokes方程的涡量-流函数形式,对置于弱电介质溶液中包覆电磁场激活板的圆柱绕流的尾迹和边界层的速度分布进行了数值研究。涡量传输方程采用ADI(隐式交替)格式,流函数方程采用FFT(快速傅里叶变换)格式,计算结果具有二阶空间精度和一阶时间精度,作了相应的实验,对实验结果进行了分析和讨论。     

用于涡流系统双分布函数的格子Boltzmann模型

提出一个格子Boltzmann模型应用于不可压缩流的涡流系统. 用两个分布函数分别定义涡量和流函数, 得到用两个格子Boltzmann方程建立的模型. 以数值为例, 得到一个比较好的方腔流动模式.     

涡方法数值模拟高雷诺数圆柱绕流

圆柱绕流问题是二维物体粘性绕流的基本问题.二维N-S方程的各种差分求解格式往往难以用来实际计算高雷诺数流动问题。80年代中期以来,Stansby等发展了70年代初Chorin提出的随机点涡法,用网格涡方法计算点涡的对流速度,大大减小了计算量,笔者采用Stansby等发展了的随机点涡法,对雷诺数2000的圆柱绕流进行了数值模拟。     

沿岸回流区污染物对流扩散的数值模拟

为研究回流区污染物堆积的机理及规律,采用高阶三点紧致差分格式对涡量-流函数方程和污染物对流扩散方程进行高精度数值求解,得到了污染物堆积与回流区岸线角度、排放点位置之间的定量关系.该方法有效抑制了通常的数值方法中对流项差分离散造成数值不稳定的现象,且程序编制简洁,运行效率高。数值计算表明该法有效、可靠.研究结果对沿岸排污点布置、减少海域污染具有一定的现实指导意义,并具有重要的理论价值.     

用有限分析法对连续弯道岸边排污的数值模拟

基于边界拟合坐标系,运用涡量-流函数方程和标准k-ε双方程湍流模型,通过有限分析方法对连续弯道的水流和污染物的扩散输移进行了数值模拟,计算了流场及岸边污染物排放的浓度场,并与实验室结果和有限体积法的计算结果进行了比较.     

圆柱绕流的格子Boltzmann模拟

针对流体力学中模拟圆柱绕流的边界层内部流动问题,采用格子Boltzmann方法,用两个分布函数分别定义涡量和流函数,得到用两个格子Boltzmann方程建立的模型。以数值为例,圆柱绕流的数值模拟结果符合经典的理论结果。与直接模拟Navier-Stokes方程相比,该方法计算模型简单,分布函数简单,易于计算。     

离散涡法模拟建筑物绕流流场

应用离散涡数值方法对高雷诺数下建筑物绕流流场进行了数值模拟，得到了不同相对位置、相对距离的建筑物在不同来流方向的分离流场和涡量分布，总结出受建筑物影响的流函数变化规律，取得了较为合理和可靠的计算结果，说明离散涡方法是研究建筑物的绕流问题的有效手段     

数值模拟圆柱统流旋涡生成、分离及演化

本文应用有限差分方法求解涡量－流函数形式的Ｎ－Ｓ方程，数值模拟了圆柱瞬时起动后流动分离，旋涡生成、脱落以及随时间推进涡街产生和长期非定常演化过程，包括对称涡的生成，二次涡形成的α结构及卡门涡街等．本文方法效率较高，数值模拟结果与实验结果吻合较好，可推广到模拟其他复杂钝体在中等或高雷诺数下的非定常绕流问题．     

有限变形理论的若干进展及其在流体力学中的相关应用

概要性地叙述了作者新近提出的"当前物理构型对应之曲线坐标系显含时间的有限变形理论"、"几何形态为曲面的连续介质的有限变形理论",前者针对介质几何形态为Euclid流形(体积形态),后者针对Riemann流形(曲面形态).类比于一般有限变形理论,上述理论均包括物理及参数构型构造,变形梯度定义及其基本性质,变形刻画,输运定理以及守恒律方程.基于上述理论提出对应曲线坐标系显含时间的流函数涡量解法,固定曲面上二维不可压缩流动的流函数涡量解法以及海面油污扩散控制方程,并给出了相关数值研究结果.     

二维N-S方程的有限分析数值解

本文导出了以涡量-流函数表示的二维恒定Navier-Stokes(N-S)方程的有限分析解,并对原有的线性化N-S方程的局部解析解进行了改进,所得的解在四个角点上满足相容性条件.对放大比为2.2:1, 2.0:1, 1.8:1, 6:1的平面突然放大问题,采用边长为0.2的正方形网格,求得了有限分析法的数值解,并研究了对称区域的不对称流问题.根据计算结果,给出了几组拟合经验公式,即平面对称放大后的旋涡长度、中心线流速的衰减变化及过渡段长度公式等.     

瞬态二维不可压缩流动N-S方程数值解

用有限元方法对粘性不可压缩流体做二维瞬态流的Ｎ－Ｓ方程的变形公式－流函数－涡量式进行求解,壁面上的涡量边界条件用时间迭代法加以确定。并分别对雷诺数Ｒｅ＝４００,８００,１０００时方腔驱动流上到了收敛的效果,所得流线图和速度矢量图均符合物理规律。说明本文所用方法正确,对分析其它流动问题有一定的参考价值。