首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 62 毫秒
1.
 传统的双时间方法在非定常计算中长时间的过渡迭代推进求解导致其计算效率相对较低,针对周期性非定常流动问题的流动特征,发展了一种基于离散傅里叶变换的高效时间谱方法,用于求解振荡翼型和机翼的非定常黏性绕流。在时空耦合的雷诺平均Navier-Stokes(RANS)方程的求解中,对流项的离散应用了Roe的通量差分格式,物理时间项的离散方法为时间谱方法,伪时间推进采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式。考虑到湍流的时空耦合效应,时空耦合的Spalart-Allmaras一方程湍流模型的物理时间项同样采用时间谱方法进行离散。为了进一步提高计算效率,当地时间步长和多重网格技术等加速收敛的措施均被采用。算例对俯仰振荡NACA0012翼型和Lann机翼的周期性非定常流场进行了数值计算。结果表明:对于周期性非定常流场的数值模拟,相比于传统的双时间方法,用时间谱方法近似物理时间项,不仅能够提高流场的计算精度,而且更能够大幅度提高计算效率。  相似文献   

2.
高阶精度格式WCNS-E-5在亚跨声速流动中的应用研究   总被引:3,自引:0,他引:3  
采用高阶精度非线性紧致加权格式WCNS-E-5和Baldwin-Lomax模型,求解雷诺平均Navier-Stokes方程,开展了典型翼型与机翼的湍流流动数值模拟研究.对方程中粘性项采用的四阶精度差分近似以及网格导数求解与边界格式的四阶精度,保证了高精度算法的实现.计算结果表明:本文算法能够准确地模拟这些翼型与机翼的亚跨声速流场,得到与实验测量十分吻合的壁面压力分布,计算结果对网格的依赖性小.  相似文献   

3.
使用多重网格的方法对二维翼型和三维机翼的跨声速颤振进行了数值求解.流场控制方程为N-S方程,湍流模型采用SA模型网格.计算网格采用结构网格,空间离散采用中心格式,使用双时间法进行时间推进.流场与结构之间的数据通过径向基函数(RBF)插值方法进行交换.通过耦合求解流场控制方程和结构运动方程得到Isogai(case A)机翼模型和Agard445.6机翼的跨声速颤振边界.通过与使用单重网格时得到的结果相比较得出,多重网格方法能够节省计算时间,提高颤振分析计算效率.  相似文献   

4.
用强耦合RANS方法模拟旋翼悬停流场   总被引:2,自引:2,他引:0       下载免费PDF全文
在旋转坐标系下,将Spalart-Allmaras(S-A)一方程湍流模型和Reynolds-averaged NavierStokes(RANS)方程耦合成一个新的RANS方程,并发展了基于多块重叠网格的强耦合RANS求解方法,用于直升机旋翼悬停流场的数值模拟.为了提高计算效率,针对多重网格方法在多块重叠网格上实施的困难,提出了一种基于重叠网格的多重网格实施方法.通过对Caradonna-Tung(C-T)和ONERA 7A旋翼悬停算例验证了发展的强耦合RANS方法和基于重叠网格的多重网格实施方法的有效性.研究结果表明:发展的基于重叠网格的多重网格方法有较高的计算效率,3层网格的加速比约为7.7;强耦合RANS法的计算精度明显高于传统的松耦合RANS方法,特别是在与阻力相关性能参数的预测中,强耦合RANS方法的预测结果更加精确.  相似文献   

5.
湍流模型离散精度对数值模拟影响的计算分析   总被引:1,自引:0,他引:1  
王运涛  孙岩  王光学  张玉伦  李松 《航空学报》2015,36(5):1453-1459
基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS),通过改变物面法向第一层网格间距,开展了剪切应力输运(SST)两方程模型不同离散精度的数值分析。主要目的是为高阶精度格式在复杂外形上的应用提供技术支撑。计算模型包含了低速NLR 7301两段翼型和高速RAE2822翼型,研究内容主要包括湍流模型的二阶精度离散和五阶精度离散两种方式对收敛历程、边界层湍流黏性系数分布、边界层速度分布、压力系数分布以及气动特性的影响。在与试验数据对比的基础上,计算结果表明:对于不同的第一层物面法向网格间距,湍流模型离散精度对低速绕流计算结果有比较明显的影响,对于高速小迎角附着流动计算结果影响不明显;相对于湍流模型二阶精度离散,湍流模型高阶精度离散网格敏感性较弱,具有更高的数值模拟精度,但收敛性略差。  相似文献   

6.
DG/FV混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度DG方法计算量、存储量小等优点,自提出以来已成功应用于一维、二维标量方程和Euler/N-S方程的求解。综述了DG/FV混合方法的研究进展,重点介绍了DG/FV混合方法的空间重构算法、针对RANS方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析,并给出了系列粘性流算例的计算结果,包括用于验证混合方法数值精度的库埃特流,以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明DG/FV混合方法达到了设计的精度阶,且相比同阶DG方法计算量减少约40%,而隐式方法能大幅提高定常流的收敛历程,较显式Runge-Kutta的收敛速度提高1~2个量级。  相似文献   

7.
绕翼型低雷诺数流动的数值分析研究   总被引:5,自引:0,他引:5  
分离泡的产生是绕翼型低雷诺数流动的一个重要特征,分离泡通常都是非定常的,会对整个流场产生极大的影响,由于分离附面层的不稳定性,很快诱发转捩,并产生湍流再附。文中通过求解雷诺平均N-S方程,数值模拟了绕Eppler387翼型的低雷诺数非定常流动,并对两方程SSTk-ω湍流模型,代数B-L模型和层流的计算结果作了比较。N-S方程和两方程湍流模型的控制方程非耦合求解,时间推进均采用近似因子方法(AF);空间离散无粘项采用ROE格式,粘性项用中心差分方法计算;计算非定常流场时,采用伪时间子迭代(τt-s)方法保证二阶时间精度,湍流模型计算时都是在固定点转捩。最后分析了转捩对低雷诺数流动的影响、分离泡的存在引起的流动不稳定性和周期性的脱出涡,数值计算给出的时间平均升力系数、阻力系数和压力分布与实验结果比较吻合。  相似文献   

8.
在Reynolds-Averaged Navier-Stokes(RANS)方程计算中耦合Michel经验判据和Chen-Thy-son转捩模型对风力机翼型S809进行计算。由RANS方程求得翼型表面压力分布作为层流边界层方程求解的输入参数,然后使用Michel经验判据分析判断层流边界层的解得到转捩点的位置,使用Chen-Thyson转捩模型得到从层流到湍流的转捩过渡区,这样随着流场的迭代求解,求解器自动更新判断出转捩点位置。对S809翼型进行全湍流和耦合转捩判断的RANS方程计算,可以看到考虑转捩判断后得到的升阻力系数与实验值吻合较好,验证了方法的可行性。  相似文献   

9.
两种湍流模型在跨声速绕流计算的应用研究   总被引:2,自引:0,他引:2  
湍流模型在对复杂流场的数值计算中起着非常重要的作用.本文采用SA一方程湍流模型和SST k-ω二方程湍流模型,通过解耦求解雷诺平均N-S方程和湍流模型方程实现对亚跨超声速湍流流场的数值模拟.对NACA0012翼型和ONERA-M6机翼跨声速绕流流场进行了计算,对压力分布和激波位置与实验结果进行了细致的比较,并分析了不同离散格式、不同网格疏密及壁面函数对计算结果的影响,在计算过程中这两个模型体现出了较好的简捷性和健壮性.  相似文献   

10.
NACA4412翼型低速绕流数值计算中湍流模型对比   总被引:1,自引:0,他引:1  
使用Spalart-Allmaras(S-A)、SSTk-ω、Gao-Yong 3种湍流模型对NACA4412翼型低速绕流进行了数值计算,研究了尾迹流动松弛效应。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。计算中对翼型尾缘流松弛效应进行了分析,比较了翼型表面压力系数、速度剖面、雷诺应力等的分布,3种湍流模型总体上能够较好地模拟NACA4412翼型低速绕流。SSTk-ω模型对流动细节及升力系数计算最好,Gao-Yong模型对翼型平均速度剖面及雷诺剪切应力分布计算最准确。  相似文献   

11.
谱方法的基本问题与有限谱法   总被引:6,自引:2,他引:6  
本文概述了谱方法的特点,讨论了谱方法的四个基本问题及解决方法,由此得出谱方法局域化的重要性,重点介绍了作者提出的有限谱法两个发展阶段,特别是最近的基于截断傅立叶积分的有限谱法,列举了有限谱法的一系列计算格式,给出了计算高超声速流,各向同性湍流,地震波的传播,分区算法,复杂形状物体,激波管问题,KdV方程,湍流和激波相干扰等问题的计算结果。  相似文献   

12.
连续体结构拓扑优化方法评述   总被引:8,自引:2,他引:6  
连续体结构拓扑优化在优化中能产生新的构型,对实现自动化智能结构设计具有重要意义。目前,连续体结构拓扑优化方法主要有:均匀化方法、变厚度法、变密度法、渐进结构优化方法、水平集法、独立连续映射方法。本文首先系统回顾了以上方法的发展历程,介绍了它们的研究现状。其次,通过对比以上拓扑优化方法对若干典型算例的优化结果,表明以上方法都有较好的减重效果。最后,对以上方法进行了总结,列出了它们的优缺点和发展方向。  相似文献   

13.
柔性决策   总被引:3,自引:0,他引:3  
石永恒 《航空学报》1994,15(4):468-471
 提出了柔性决策的思想,并对决策科学的发展进行了回顾。为解决柔性决策问题,提出了极大-极小方法、人机交互方法和集包含方法。  相似文献   

14.
一种改进型洞映射法   总被引:1,自引:0,他引:1  
针对传统洞映射法存在挖洞曲面必须封闭的局限性和消耗大量内存的缺点,提出了改进型洞映射法。主要做了两点改进:1)采用“反向标记法”判别洞映射笛卡尔网格单元属性,使得挖洞曲面不再局限于封闭的情况;2)运用“矢量射线法”进一步判断洞映射笛卡尔网格边缘单元内的计算网格点属性,使得笛卡尔网格分辨率低时仍能精确挖洞。研究表明:相比于传统洞映射法,改进型洞映射法保持了很高的计算效率和精准挖洞的能力,还极大地降低了内存消耗,并且具有较强的可靠性。  相似文献   

15.
基于水平集方法的连续体结构拓扑优化方法研究   总被引:1,自引:0,他引:1  
考虑到现有各种连续体结构拓扑优化方法的优缺点,提出了一种新的连续体结构拓扑优化新策略。该方法将无网格法和水平集方法结合起来,充分发挥各自的优点,使得该拓扑优化新策略具有良好的收敛性和计算精度。  相似文献   

16.
本文提出了选择复连通域边界(多环路)计算J积分的一种新方法,成功地解决了单连通域内有铆接点时计算J积分的问题。文中还提出了一种“虚设铆钉法”,采用这种方法,可利用同一网格有限元模型计算不同裂纹长度下裂纹尖端的应力强度因子,这大大简化了模型准备工作并减少了计算结构总刚度矩阵的工作量。本文采用刚度导数法、J积分法及能量释放率法,计算了含裂纹铆接加筋板的应力强度因子,通过计算结果对比,对各种方法作了粗浅的评价。  相似文献   

17.
18.
海面(或地面)环境中目标电磁散射特性,对实际战场环境下军用目标的雷达探测、雷达特性形成以及雷达目标识别等技术具有基础支撑作用。因此,快速准确地分析目标和粗糙面复合电磁散射特性具有重要理论意义和应用价值。文章采用基于矩量法(MOM)和物理光学法(PO)的混合方法MOM-PO分析了目标和粗糙面复合电磁散射问题。与传统的矩量法相比,该方法在保证精确度的前提下大大加快了计算速度。  相似文献   

19.
关于电磁场的自然单元法   总被引:1,自引:0,他引:1  
有限元法和无网格法在电磁场数值计算中已经得到了广泛的应用,然而有限元法存在前处理网格剖分问题,无网格法存在计算时间长、边界条件和不连续面难处理等问题。针对以上问题,本文提出了利用自然单元法求解电磁场的方法,该方法根据场域中的离散点的信息构造自然邻点插值函数,最后算出场值。通过算例证明了该方法在电磁场计算中的可行性。  相似文献   

20.
典型翼型参数化方法的翼型外形控制能力评估   总被引:3,自引:0,他引:3  
现代飞行器的精细化设计需求要求翼型参数化方法具有精准的翼型外形控制能力。依据标准的风洞模型公差要求,以五个典型翼型作为测试翼型,分别对Hicks-Henne参数化方法、样条参数化方法、PARSEC参数化方法和CST翼型参数化方法的翼型外形控制能力进行评估,依据每个翼型参数化方法的评估结果并结合其固有特点进行适应性分析。结果表明:四种翼型参数化方法对翼型外形的控制能力存在较大差别,具有不同的适用范围。评估典型翼型参数化方法的所得结果可指导设计人员依据具体的设计需求选择适当的翼型参数化方法。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号