喷雾仿真中的网格依赖性

基于离散油滴模型,比较了原始喷雾模型与修正喷雾模型的网格依赖性.计算结果表明:原始喷雾模型由于采用O'Rourke碰撞模型以及最近节点法耦合算法,使得预测结果出现明显的网格依赖性.相反,修正喷雾模型采用独立的拉格朗日碰撞单元降低碰撞过程的网格依赖性,联合使用拉格朗日插值算法以及气相射流模型降低气液耦合过程对于网格的依赖性.因此修正喷雾模型能获得较好的低网格依赖性从而保证粗大网格下的计算精度.      

圆管大应变分析的自适应单元网格划分技术

以圆薄壁钢管在冲击作用下的屈曲过程为例,描述了钢管在冲击载荷作用下,局部发生的大变形情况。在ANSYS／Ls—dyna计算中采取了不同的网格划分形式．对结果进行了变形和应力的分析比较。证明在大变形情况下,由于网格划分形式的变化,对计算精度和计算结果都有较大的影响。该方法在分析薄壁防撞杆的受冲击过程中,对局部受损情况的分析,以及防撞杆的安全性分析都有帮助。     

面向混合网格高精度阻力预测的熵修正方法

Roe通量差分分裂格式具有耗散小、分辨率高等特点,在亚跨超声速流场模拟中得到广泛应用,但空间离散时需要对Roe平均矩阵的特征值进行熵修正,使用熵修正会增大耗散,影响阻力的求解精度。本文针对非结构混合网格中Roe格式熵修正的特点,通过改进传统的Roe格式Harten-Yee熵修正方法,提出了一种可提高非结构混合网格黏性计算精度的Harten-Yee熵修正改进方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析,改进后的熵修正方法残差收敛特性与原始Harten-Yee熵修正一致,计算精度提高,计算结果和无熵修正时基本一致,说明改进后的熵修正方法既保留了使用熵修正带来的程序鲁棒性等优点,同时把熵修正对阻力预测精度的影响降到了最低,验证了方法的有效性。采用改进后方法对第3届AIAA阻力会议的计算模型DLR-F6翼身组合体进行了详细的模拟,分析了网格收敛性和雷诺数的影响。结果表明：改进的Harten-Yee熵修正,更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。     

显式代数雷诺应力模型在超声速/高超声速流场计算中的应用

为了提高显式代数雷诺应力模型在高超声速流场中的计算精度,以修正后的CG k-ε模型作为湍流时间尺度决定方程,将Wallin和Johansson显式代数雷诺应力模型(WJ2000模型)应用于超声速/高超声速流场数值计算,并与线性湍流模型的计算结果进行了对比分析。结果表明WJ2000模型对高超声速压缩拐角附近壁面压力和热流分布的预测要优于线性湍流模型,同时网格依赖性明显低于后者;对CG k-ε湍流时间尺度决定方程所做修正不影响WJ2000模型对超声速流场的预测结果,但使其对高超声速流场的预测精度明显提高。     

高雷诺数槽道湍流的壁面模化大涡模拟研究

选取基准壁湍流的槽道流动,研究了多种模型的壁面模化大涡模拟.模型包括经典的大涡模拟、Spalart-Allmaras、分离涡模拟和一种动态混合模型.基于摩擦速度的雷诺数范围为395~12000,采用3组粗糙网格,流向和展向维数分别同取37,49和65,法向维数保证y+(1)~1.主要研究平均速度、雷诺切应力分布、详细分析了各模型的特性差异并展示了相应的湍流结构.研究表明:在高雷诺数粗糙网格下,大涡模拟失去求解精度,分离涡模拟出现对数律不匹配,动态混合模型的计算接近直接数值模拟,其对数率区可解应力约占雷诺切应力的93%,边界层外层可解应力约占99%.这说明合适的混合模型可以在经济成本下保证计算精度,具有解决实际问题的潜力.      

基于伴随方程的网格自适应及误差修正

基于流场方程的离散伴随优化理论和三维非结构网格,建立了网格自适应技术和目标函数误差修正方法。详细研究了用流动伴随变量进行目标函数的误差估计和修正技术,构造了适用于格心格式有限体积法的流场变量插值技术和网格单元剖分判据,初步实现了网格物面投影和空间单元优化,发展了适用于有限体积法的整套网格自适应方法。对NACA0012翼型和ONERA-M6机翼绕流进行了自适应数值模拟,并对升、阻力等目标函数进行了误差修正。数值结果表明,本文自适应方法能正确地捕捉到影响目标函数计算精度的敏感区域,网格自适应和误差修正两项技术显著提高了升、阻力等目标函数的计算精度。     

三维裂纹前缘参数化网格模型研究

为实现三维疲劳裂纹扩展自动模拟及获得更高精度的应力强度因子,在实体参数化建模的基础上,自编程开发了两种参数化裂尖网格模型,分别对应奇异元法和虚拟裂纹扩展法(VCCT法)计算应力强度因子,且两种参数化网格可互换以便相互验证。两种裂纹模型均精心设计,裂纹前缘网格正交,过渡均匀,疏密可调。与解析解的对比结果表明,两种参数化网格模型均有很高的计算精度,可为一般工程结构获得较高精度的应力强度因子结果提供技术支持。     

轮盘径向破裂转速计算方法分析及修正

分析了轮盘径向破裂转速理论计算方法及常用平均径向应力计算方法的含义及其特点.且依据理论计算方法的含义,提出了基于有限元计算结果的轮盘径向破裂转速计算修正方法.通过实例应用,证明了相对于理论方法计算结果,修正方法较常用平均径向应力方法计算精度更高,其计算精度由6.6%提高到0.5%以内,且修正方法较理论方法使用更方便、工程适用性更强.      

风扇叶片飞失显式动力学仿真网格与时间步长研究

为了研究大涵道比涡扇发动机风扇叶片飞失仿真中网格尺度与时间步长对于瞬态显式动力学分析结果的影响,采用有 限元法中单元计算理论分析的方法分析了网格尺度和时间步长在计算稳定性方面的理论联系。采用多因素试验设计方法组织了 针对叶片应力分布、转子轴心轨迹和支点外传振动应力3个风扇叶片飞失的典型输出参数的网格尺度和时间步长的影响研究方 案,开展了典型物理过程、模型简化和数值计算无关性研究。针对时间步长为7×10-7、5×10-7、3×10-7和2×10-7s,网格尺度为40 、30、 20和15 mm的计算结果进行分析,结果表明：对于给定物理过程和目标时长,存在着网格尺度和时间步长的门槛值,如：2×10-7 s 时间 步长和30 mm网格尺度,超过该值进一步细化网格和减小时间步长对于精度提升不明显;显式动力学的计算原理决定了目标时间 越长偏差累计越大,所需的网格精度就越高;碰摩过程对于风扇叶片飞失仿真的应力和外传振动偏差影响较大。     

带旋转修正的弹簧-TFI混合动网格方法

网格存在大变形时,传统的超限插值(TFI)动网格方法易于引起网格正交性问题.依据几何关系并考虑其插值特点,提出了一种带旋转修正的TFI动网格方法.对结构网格块进行分块,在子块角点间建立弹簧元,并基于弹簧类比法计算角点位移;子网格内部结点位移采用修正后的TFI动网格方法进行计算,从而形成具有弹簧-TFI混合特征的动网格新...     

气体动理论BGK格式的网格自适应方法

为了提高气体动理论BGK（Bhatnagr-Gross-Krook）格式在超声速流动问题计算时激波捕捉的准确性与计算效率,提出了一种适用于气体动理论BGK格式的网格自适应加密方法。该方法采用基于四边形的链表技术来描述网格的拓扑结构,在物理量重构过程中,使用了在四边形网格中表现优异的van Leer限制器,以保证粗细网格过渡处物理量重构的精度。用跨声速翼型绕流（马赫数Ma=0.85）、超声速前台阶流（Ma=3）和高超声速圆柱绕流（Ma=8.03）等多个典型算例验证了BGK自适应网格方法。计算结果表明,自适应网格BGK方法在保证数值精度的前提下,可大幅度提高计算效率。这为该方法用于高效地解决复杂问题提供了一种选择。     

A mesh optimization method using machine learning technique and variational mesh adaptation

Computational mesh is an important ingredient that affects the accuracy and efficiency of CFD numerical simulation. In light of the introduced large amount of computational costs for many adaptive mesh methods, moving mesh methods keep the number of nodes and topology of a mesh unchanged and do not increase CFD computational expense. As the state-of-the-art moving mesh method, the variational mesh adaptation approach has been introduced to CFD calculation. However, quickly estimating the flow fi...     

用于气动声学计算的非均匀网格紧致差分格式

为克服传统紧致差分格式在数值求解非均匀网格问题时产生的寄生波,构造了一种新的高精度紧致差分格式.通过泰勒展开分析方法,详细给出了格式系数的通用形式;利用傅里叶分析方法,分析了数值耗散、色散误差.以3对角6阶精度紧致差分格式求解均匀扰动网格问题为例,计算表明:色散值和耗散值随扰动因子的增加而更加趋近于精确值;当扰动因子大于0.213时,格式不稳定,当扰动因子小于等于0.213时,格式渐近稳定;对一维对流波和二维波传播的模拟计算所得数值解与精确解吻合,体现了该格式在求解非均匀网格问题时的优越性.      

三维弹塑性问题的有限元分析

本文用有限元法对三维弹塑性的应力与应变进行了分析,文中将非线性求解和数值积分做了改进,使计算时间缩短;三维元素的各种信息自动生成,提高了输入数据的精度和效率。用本法计算的各种实例与实验数据进行比较,二者吻合甚佳,能满足工程上所要求的精度。     

Mesh generation: Art or science?

Mesh generation has evolved to the point where highly complicated domains can be covered by a variety of mesh types including hexahedral, tetrahedral and overset meshes. The application of these methods to computational aerodynamics has become a routine exercise and numerical predictions over complete aircraft now complement experimental results obtained from wind tunnels. This paper surveys the main developments that have taken place and traces the evolution of mesh generation over the last 35 years. This is followed by an assessement of the accuracy of Navier Stokes codes that are currently in use for predicting the drag of an aircraft at transonic cruise. The relationship between solution accuracy, mesh size and mesh type is examined in some detail and the implications for further research are discussed.     

循环加载下的热弹塑蠕变有限元分析

本文考虑在循环加载条件下, 采用混合强化模型, 导出相应的热弹塑性蠕变增量本构方程及有限元格式。研制出平面应力程序, 用牛顿迭代法求弹塑性权, 并提出精度高的屈服吸面校正方法, 程序经考核, 可供工程设计使用。      

Euler方程的自适应多重非结构网格计算

介绍一种基于四面体网格的多重网格技术,用于求解三维Euler方程定常解,具有十分明显的加速收敛效果。求解Euler方程的基本算法采用有限体积中心差分格点格式。多重网格技术采用一组几何上相互独立的粗细不同的网格。流动变量、余量和修正量在相邻粗细网格间的相互传递通过线性插值实现。采用自适应技术对网格进行加密,可以获得较优的网格,并作为一套细网格,进行多重网格计算。数值实验表明本文自适应多重网格法十分有效。     

中长寿命轮盘应力寿命及可靠性分析方法

针对轮盘榫槽结构理论应力集中系数难确定、采用传统应力寿命法估寿误差大问题, 提出一种改进应力寿命预估法.改进方法包括:提出应力修正系数概念和一种适用于轮盘榫槽结构的理论应力集中系数确定方法, 以及应力比修正方法.一旦确定了试棒的应力修正系数, 改进应力寿命法可采用光滑棒试验数据预估复杂结构寿命, 并考虑应力集中影响.进而建立了基于改进应力寿命法的应力寿命可靠性分析方法, 并对某风扇盘进行了疲劳寿命可靠性分析, 计算概率寿命与试验寿命吻合良好.      

表面裂纹的分区混合元分析

 本文对三维体表面裂纹问题进行分析,应用应力和位移的渐近解在表面裂纹区域构造一个奇异应力单元,采用分区混合有限元法计算裂纹顶端各个位置的应力强度因子。     

圆柱筒轴向载荷下单元尺寸效应数值模拟

利用显式非线性有限元方法对薄壁圆柱筒轴向载荷下的动力学响应进行了研究,模拟了圆柱筒在轴向冲击作用下的动力屈曲变形及其发展过程。分析比较了不同单元尺度下冲击模型的数值计算解与理论解,揭示了在进行有限元模拟时存在的网格尺寸效应,为研究有限元网格尺寸与有限元求解精度的内在联系提供参考;为在保证数值解满足工程实际精度要求的前提下,确定合理的单元尺度,提高有限元分析效率进行了有益的探索。