网格生成方程中非线性系数对网格正交和分布的影响研究

通过数学分析，对椭圆型方程非线性系数对生成网格的正交性和空间分布的影响进行了研究，理论分析和数值结果均表明，系数β具有控制网格正交性的作用，系数α和γ具有调节网格空间分布的作用，且其调节的强增和作用与其比值有直接关系，基于此提出了一种混合型边界条件的给定方法，较好地实现了网格正交和疏密分布的控制。     

基于局部无网格的混合算法研究

研究了一种求解任意外形绕流的局部无网格混合算法.该算法计算区域整体采用了直角网格,只在物面附近嵌入局部无网格处理.整体直角网格的规则性与正交性,不仅使得网格生成快,而且计算格式实施简单,避免了与非正交性相关的计算项,从而引入误差相对较小.局部无网格处理只要求物面附近布点离散,具有灵活性,使得发展的算法适合处理任意外形.先以理想不可压流动问题为例,给出了其控制方程Laplace方程的具体求解实施过程,并给出了直槽道内柱体绕流的计算结果;接着把方法推广用于求解可压缩流动的Euler方程,对喷管等典型流动问题进行了数值模拟,并与传统有限体积法结果进行了比较,在捕捉的波系结构、表面压力系数分布等方面具有良好的一致性,表明发展的新方法是可行的.     

H型叶栅通道流计算网格的改进

提出一种壁面正交、尾迹区网格节点分布逐渐均匀化叶栅通道网格生成方法。为了研究壁面处网格正交性对数值计算结果的影响，推导出壁面附近流动参数一阶导数离散误差与网格线交角的数学关系。采用N-S方程分别对平板紊流附面层流和叶栅通道流进行计算分析。结果表明：壁面处网格正交性对附面层内流动计算预测精度影响较大，而对势流区流动影响很小；尾迹区切向网格节点逐渐均匀化分布能更准确计算尾迹区流动。     

复杂空气动力外形网格生成技术

通过对用于网格生成的Ｌａｐｌａｃｅ方程进行抛物化差分离散,发展了一种快速抛物型方程贴体及与边界正交网格生成方法,另外,采用源项修正技术求解用于网格生成的Ｐｏｉｓｓｏｎ方程,发展了一种可实现对网格线与边界间距及正交性进行双重控制的椭圆型方程方法,运用这种两种基本方法以及这两种方法的混合,分别生成了真实歼击机全机外形,带十字叉翼导弹以及双发,下单翼运输机外形的计算网格。并已成功地秀于跨声速得杂流场的娄     

结冰计算中翼面网格重构研究

对结冰数值计算中的网格生成进行了深入研究,针对结冰后几何边界的复杂性及边界依赖于时间步长变化对网格生成的影响展开了分析.网格生成是结冰计算的基础,文中基于求解椭圆形偏微分方程(PDEs)生成了未结冰时翼面周围的背景网格,采用Hilgenstock源项构造方法,合理控制了结冰边界周围的网格间距和正交性.模型表面结冰后,为...     

一种基于求解椭圆型方程的结构动网格生成方法

基于椭圆型网格生成法,实现了一种简单高效的贴体结构动网格生成方法,可用于具有移动边界问题的非定常流动数值模拟。该方法提出,在网格变形过程中,Poisson方程需要的控制网格间距和正交性的源项可以通过提取已知的静态网格源项直接得到,并在整个动网格生成过程中保持不变。因此,在椭圆型网格生成中需要通过外迭代确定源项的过程可以得到省略,而且该方法不需要人工指定参数。这使得方法具有高效和易于嵌入到已有程序中的特点。数值模拟结果证明,采用这种方法获得的网格能够较好地保持静态网格原有的正交性和光滑性,在相同迭代步数约束下,网格求解效率低于传统弹簧模拟法,但鲁棒性优于弹簧模拟法。     

旋转坐标系N-S 方程无网格/网格混合算法

为了模拟涡桨发动机等绕固定轴旋转的桨叶流场,发展了一种求解旋转体非定常黏性绕流的无网格/网格混合算法。算法基于求解旋转坐标系纳维尔-斯托克斯方程展开,避免了旋转角速度特征物理量的插值运算;计算域采用整体网格和物面附近局部无网格离散,通过引入无网格和网格对偶点,实现了混合算法绕流信息的跨区交换,并给出了一种对偶点调整选取的方法;基于无网格点云空间导数逼近方法离散控制方程,并结合隐式LU-SGS算法,给出了混合算法求解非定常问题的双时间步隐式推进格式。结果表明:所提算法通过2维振荡翼型非定常绕流、旋转圆柱黏性绕流和模拟发动机桨叶旋转运动的3维悬停旋翼绕流进行了考核,所得升力系数等重要气动数据的变化趋势与文献试验或计算值一致,典型截面处所示的激波强度和位置亦与文献值吻合,展现出算法在模拟旋转部件绕流问题方面具有广阔的工程应用前景。     

强制偏流喷管的三维网络生成

运用泊松方程实强了制偏流喷管内的三维网络生成，方程中的源项采用了与边界形状和边界网格点分布无关的数学形式，源项内的自由参数由横向网格线与边界的交角及横向网格线在边界处的局部曲率两个约束条件来确定，实现了由边界上的网格分布直接控制内部的网格点，使生成的网格具有良好的正交性和合理性。     

多部件组合体分块网格生成技术及应用

给出多部件组合体分块网格生在的一种方案及实践,包括交界面网格生成技术、空间网格生成方法及质量控制。以一个平尾和机翼具有一定高度差的机身－机翼－平尾－立尾组合体外形为例,生成与公共交界面或边界面正交或接近正交的分块网格。用于Ｅｕｌｅｒ方程跨声速来流几个迎角的计算,结果与实验符合较好。     

双边界代数法在喷管内流场网格生成中的应用

用一种双边界代数方法生成二维网格，通过在超限插值函数中引入边界的一阶导数实现网格与边界是的正交，采用几何一致性条件保证物理区域不变形，发展一系列拉伸方程控制网格在边界或给定位置的间距。应用此方法生成了满意的喷管内流场网格，网格生成效率明显优于微分方程网格生成方法。     

一种适用于飞行器复杂外形的表面与空间网格生成方法

本文通过建立反映机身横截面外形特征的控制点概念，及在控制点间分配网格点数方式，给出了一种非常实用、交互式的表面网格生成方法，所生成的网格失真度小，物体保形好，生成适用于CFD计算需要的表面网络花费的时间少，适用性强。基于无限插值理论，引入本文提出的网格正交控制、物面法向量控制、加权平均光顺措施，给出了一种改进的空间网格生成方法，有效克服了传统无限插值网格生成方法在复杂外形网格生成方面的缺陷，有效改善了网格的法向疏密性、贴体性及周向网格的均匀性，对真实飞行器部件网格生成非常迅速。以某飞行器复杂机身表面与空间网格生成，及无人机翼身组合体空间网格生成为例，检验了本文网格生成方法的实用性和有效性。     

基于控制理论和NS方程的气动设计方法研究

研究了基于控制理论和NS方程的气动设计方法,针对给定的目标函数表达形式,应用该设计理论在计算坐标下详细推导了相应的共轭方程及边界条件具体表达形式,以及梯度方程求解表达式,通过合理的数学变换,导出了共轭方程在笛卡尔坐标系下的直观表达形式,发展了有效的共轭方程数值求解方法,通过流动控制方程数值求解、共轭方程数值求解、目标函数对设计变量的梯度求解和优化算法等方面的有效结合,研究与发展了一种新的气动设计方法,以二维机翼气动设计为例,成功进行了亚、跨音速情形下的相关设计算例研究,研究结果表明应用控制理论和NS方程的气动设计方法在设计理论、适用性以及时间花费等方面都有着很好的特色和优点,且设计结果也更为可靠.     

运动激波自适应网格算法中权函数问题的研究

变分原理的自适应网格技术被应用到运动激波问题的求解上,在解的大梯度区自动加密网格,从而非常成功地算出了激波。通过分析发现权函数选择和取值对自适应网格技术至关重要     

应用嵌入式H型网格技术计算孤立翼型绕流场

提出一种嵌入式H型网格技术,用于生成孤立翼型和叶栅的计算网格。与传统的H型网格相比,采用嵌入式H型网格时,翼型头部及尾部附近的网格畸变很容易得到控制,有益于提高计算的准确性。      

解跨音速升力翼型的有限元法

 在应用解全速位方程的最小压强积分有限元法求解绕升力翼型的跨音速流动时,将不可压流中求解绕升力翼型的耦合单位环量流动和无环量流动的解法推广到可压流中。为了确定环量,本文所用Kutta条件是:在后缘处,气流流向平行于后缘角二等分线。因有限元法对网格无正交性要求,因而可在椭圆变换前后进行剪切和延伸变换。这种网格生成法易于构成适用于复杂形状的有限元网格。通过计算并将其结果与文献中的数据比较,表明这种方法应用方便且有较快的计算速度和较高的计算精度。     

气冷涡轮级流场的数值模拟方法与试验验证

采用具有TVD(total variation diminishing)性质的三阶精度Godunov格式、自由型曲面网格生成技术以及分区网格算法,对某型涡轮级进行考虑冷气掺混的全三维Navier-Stokes(N-S)方程数值求解,并将所得的结果与试验数据进行对比分析.通过分析数值模拟结果研究了该涡轮级所具有的气动特点.结果表明:所开发的快速、高效冷气掺混网格生成方法以及任意分区的流场求解算法可以满足工程上对气冷涡轮级的总体性能的快速估算及流场结构的详细描述.      

固体火箭发动机燃烧室喷管一体化网格生成

给出固体火箭发动机燃烧室喷管一体化网格策略,用二维方法生成平面网格经旋转或堆砌形成空间网格。对网格正交性受几何域形状限制的情况进行了讨论。给出了一种燃烧室－喷管外形的网格结果。