壁面凸体层流绕流流场的数值计算

采用贴体坐标变换方法将直角坐标系中不规则的流动区域变换成规则的计算区域。对求解原始变量的ＳＩＭＰＬＥＲ方法进行了改进，应用有限分析法给出了由正常网格求解原始变量的迭代方法，然后计算了台阶绕流和任意凸体绕流的流速、压力场。计算结果表明该方法简单有效，压力、速度收剑效果好，计算精度较高，可以用于计算各种流动区域的流场。     

直角网格壁函数的二维高雷诺数机翼绕流计算

针对传统直角网格方法无法准确计算高雷诺数下近壁面区域的流场及剪切应力问题，本文基于自主开发的二维直角网格CFD求解器，采用添加壁面函数的方式模拟了高雷诺数下的机翼绕流。数值模型基于直角网格有限差分方法对N-S方程进行离散，整体求解器的精度在空间和时间上分别达到二阶。通过对低雷诺数下圆柱绕流及机翼绕流的模拟计算，验证了此求解器的精度和可靠性。在高雷诺数下利用壁面函数修正近壁面的速度剖面，使近壁面区域的速度分布更符合实际流动。研究结果表明：本文方法能够有效求解二维高雷诺数机翼绕流问题。同时利用商业计算流体力学软件Star-CCM+采用贴体网格对同样问题进行模拟，数值模拟结果表明本文方法的网格需求要大于贴体网格。     

基于微分求积法的三维非恒定、不可压N-S方程的数值计算模型

采用SIMPLE算法与用高阶多项式来逼近的微分求积法（PDQ）结合，建立了在非交错网格系统上求解三维、非恒定、不可压、原始变量形式Navier-Stokes方程的计算方法．在本算法中，流动控制方程中的空间导数用PDQ法进行离散，流动控制方程以SIMPLE算法的思想进行求解；并提出了在非交错网格系统中保证在边界上满足连续性方程以及补充压力校正方程的边界条件的方法．     

用贴体坐标方法计算后向台阶内圆柱绕流问题

本文采用涡量一流函数方法，应用贴体坐标技术来研究后向台阶粘性流动，利用作者提出的多连通区域分块计算方法研究了流动区域内含一圆柱体的粘性流动问题。在圆柱体表面上用粘性条件确定流函数值。得到后向台阶内绕流问题的数值解。此法可推广到求解有限区域内多体统流问题。     

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

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

椭球体有攻角绕流对称面上三维层流边界层的数值计算

本文提出了求解回转体绕流初值问题的一种简便的新方法,并利用这种方法对椭球体的有攻角绕流进行了计算分析,讨论了来流攻角和椭球体几何形状对流动的影响,同时还对流动分离和周向回流进行了分析,计算结果与实验结果吻合良好.     

激波绕尖劈障碍物运动粘性流场的数值模拟研究

本文以数值方法研究了运动激波绕射的非定常粘性流的发展过程．计算中采用二阶ＴＶＤ差分格式求解二维非定常薄层Ｎａｖｉｏｒ－Ｓｔｏｋｅｓ方程，得到了清晰的流场结构图谱，分析了产生局部粘性分离的现象．研究结果表明，本文的方法可以很好地描述激波绕射的复杂流动，对流场间断有较高的捕获能力，同时表明在激波障碍物干扰的绕射场研究中局部粘性效应作用是不可忽视的重要因素     

原始变量法计算封闭腔内自然对流

在文献[1，2]的基础上，利用原始变量法求解封闭腔内自然对流换热问题．对求解域采用非等距交错网格剖分，利用泰勒级数于网格点展开取二阶精度进行(u,v,θ)方程各项的离散，采用SIMPLE方法对压力及压力修正进行求解．不同瑞利数(Ra)条件下的数值试验显示，该数值方法物理概念清晰，计算稳定,可通过调节网格的疏密等方法，达到改善求解的收敛及提高求解的精度等目的．     

绕流柱形对层板流阻特性影响的数值模拟研究

利用有限体积法对三维可压缩的N—S方程进行离散，对内部绕流柱形分别为圆形、方形和菱形的多孔层板冷却结构进行了数值模拟。网格划分采用非结构化网格，湍流模型为Realizableκ-ε双方程模型，近壁处湍流采用壁面函数法处理，采用SIMPLE算法求解速度与压力的耦合。计算获得了3种冷却结构内部各气流参数的三维分布和流动阻力特性。结果表明层板内部的流场结构十分复杂，射流冲击后在扰流柱前反卷形成驻涡。3种绕流结构中方形结构流动阻力最小，菱形最大。将计算结果与实验结果进行了比较，两符合良好。     

黏性非定常圆柱绕流的升阻力研究

对于单圆柱绕流问题，以前的分块耦合法一般将整个流场区域划分为八个子区域，同时每一子块的内部拟边界与相邻块重叠一层网格，且节点数一致，以便实现数据传递．这种分块方法紧贴圆柱面上也会有拟边界点，所以不利于以后的动网格计算，使用新的分块耦合求解方法来数值模拟单圆柱以及不同间距下的串列双圆柱绕流情况，克服了过去分块方法造成的边界附近的数值奇性问题，并将该方法扩展用于串列双圆柱绕流的数值模拟，研究分析了改变双圆柱中心间距对上下游圆柱的升阻力系数和脉动频率所产生的影响，为进一步研究涡致振动提供了依据，使用壁面涡量积分法计算升阻力系列化脉动，精确计算出了脉动周期和幅值，计算结果与实验数据符合较好。     

起伏地表条件下有限差分地震波数值模拟 ——基于广义正交曲线坐标系

为了解决常规有限差分法在处理起伏地表自由边界条件时存在的需进行复杂坐标旋转和插值运算,以及对起伏地表进行阶梯状离散近似所产生的虚假绕射波等问题,把待求解的物理域离散成贴体正交曲线网格,将该网格映射到计算域的矩形规则网格上,在该计算域内求解一阶弹性波速度-应力方程及相应的自由边界条件,实现了广义正交曲线坐标系下的起伏地表地震波数值模拟问题.该算法简化了自由边界条件的实施,避免了因起伏地表阶梯化网格离散近似所产生的虚假波.数值模拟结果表明,该方法是可行的,可以用于起伏地表条件下地震波传播规律的研究.     

基于正交贴体网格的VTI介质地震波数值模拟

为了解决规则网格有限差分法在处理起伏地表自由边界条件时需进行复杂坐标转换和插值运算的问题,以及对起伏地表进行阶梯状离散近似所产生的虚假散射波问题,将计算流体力学中的正交贴体网格生成技术引入到起伏地表下VTI介质的网格剖分中,采用基于同位网格的DRP/opt MacCormack有限差分法模拟起伏地表下VTI介质中的地震波场。网格的正交性使实施自由边界条件时无需做复杂的坐标转换和插值运算。数值算例表明,本文提出的网格剖分法能够有效消除阶梯状网格产生的虚假散射波,从而提高起伏地表下VTI介质地震波数值模拟精度。     

射流泵内部流动的数值模拟

全面考虑了射流泵的不规则几何形状和各结构尺寸对其性能的影响 ,采用贴体坐标变换技术 ,结合混合有限分析法和k-ε紊流模型 ,给出了一套有效的射流泵内部流动计算的数值方法 .算例表明 ,其结果与实验值吻合较好 .     

采用非规则区域有限分析5点格式计算均质渗流

针对有限分析法边界适应性差的问题,采用坐标旋转技术,建立了拉普拉斯方程非规则区域的有限分析5点格式,对均质渗流问题进行了计算,给出了土坝渗流的浸润线和等水头线,并与有限元方法的计算成果进行了比较,验证了该方法的实用性.该方法可以作为有限分析方法对解决不规则区域的一个扩充.     

二维扩散方程的5点格式有限近似解法

提出求解二维扩散方程的曲线网格有限近似解法.对非规则区域定常扩散方程和规则区域非定常扩散方程的计算结果与精确解比较,表明该方法既可以用于求解规则区域矩形网格的扩散问题,又可以用于求解非规则区域曲线网格的扩散问题.最后模拟计算了典型的坝基渗流流场,计算结果与实验结果吻合较好,表明该方法具有计算简便、精度高、适应强等特点.     

Characteristic-Based Time Domain Method for Cylindrically Conformal Microstrip Patch Antennas

The characteristics of a cylindrical conformal microstrip patch antenna are analyzed by using the characteristic-based time domain （CBTD） method. A governing equation in the cylindrical coordinate system is formulated directly to facilitate the analysis of cylindrically conformal microstrip patch antennas. The algorithm has second-order accuracy both in time and space domain and has the potential to eliminate the spurious wave reflection from the numerical boundaries of the computational domain, Numerical results demonstrate the important merits and accuracy of the proposed technique in computational electromagnetics,     

非正交曲线坐标系下梯形区域的流场计算

利用了生成贴体坐标的代数法将一梯形区域从物理平面转化到计算平面,将物理区域中的边界条件向计算区域进行了一一对应转换.运用贴体坐标转换方程对其流场控制方程进行离散和求解,生成了梯形区域物理平面的贴体网格,同时应用非正交曲线坐标系对任一角度的梯形区域的流场进行模拟计算,并在同一条件下与正交坐标系所模拟的流场进行比较,说明了非正交曲线坐标方法的正确性和可行性.利用非正交曲线坐标可以更准确地模拟中包内钢水流动的变化规律,从而更好地确定最佳包内钢水浇注温度.     

非正交曲线坐标系下的流场计算

利用生成贴体坐标的代数法将一梯形区域从物理平面转化到计算平面,将物理区域中的边界条件向计算区域进行了一一对应转换.运用贴体坐标转换方程对其流场控制方程进行离散和求解,生成了梯形区域物理平面的贴体网格,同时应用非正交曲线坐标系对任一角度的梯形区域的流场进行模拟计算.并在同一条件下与正交坐标系所模拟的流场进行比较,说明了非正交曲线坐标方法的正确性和可行性.     

极坐标系下黏性流动的拟谱方法

针对拟谱方法在极坐标下遇到的坐标奇异以及压力边界条件难处理等问题,提出了一种新的基于极坐标系下的拟谱方法.该方法在一致网格上直接求解非定常原始变量形式的N-S方程,采用算子分裂法来解耦速度和压力,在径向采用Chebyshev-Radau配置点,在角方向采用标准的Fourier配置点.同时对压力的谱导数矩阵进行了变换,使边界上的压力点不出现在计算方程中.把该方法应用于圆截面的曲线管道的层流流动中,数值实验结果表明,该方法成功解决了极坐标系中的坐标奇异问题;对压力谱导数矩阵的变换,避免了对压力提非物理的边界条件,保证了解的精度,同时也消除了虚假压力模式     

非正交曲线坐标系下的梯形物体导热过程温度场的数值计算

利用生成贴体坐标的代数法将一梯形区域从物理平面转化到计算平面,将物理区域中的边界条件向计算区域进行了一一对应转换.运用贴体坐标转换方程对其温度场控制方程进行离散和求解,生成了梯形区域物理平面的贴体网格,同时应用非正交曲线坐标系对任一角度的梯形区域的温度场进行模拟计算.并在同一条件下与正交坐标系所模拟的温度场进行比较,说明了非正交曲线坐标方法的正确性和可行性.