首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 78 毫秒
1.
超声速平板边界层斜波失稳转捩过程研究   总被引:6,自引:0,他引:6  
马汉东  潘宏禄  王强 《力学学报》2007,39(2):153-157
以5阶迎风和6阶对称紧致格式混合差分求解三维可压缩滤波Navier-Stokes方程,对Mach 数为4.5, Reynolds数为10000的空间发展平板边界层湍流进行了大涡模拟. 时间推进采用 紧致存储3阶Runge-Kutta方法,亚格子尺度模型为修正Smagorinsky涡黏性模型. 通过在 入口边界叠加一对线性最不稳定第一模态斜波扰动,数值模拟得到了平板层流边界层失稳转 捩直至湍流的演化过程. 对流场转捩过程中瞬时量及统计平均量的分析表明,数值模拟结果 与理论吻合,得到的Y型剪切层、交替\Lambda涡结构以及转捩后期的发卡涡结构的发展 变化与相关文献结果一致,湍流流谱定性合理.  相似文献   

2.
激励小尺度模式在湍流圆管射流中的应用   总被引:1,自引:0,他引:1  
严红  苏铭德 《力学学报》2000,32(5):513-522
采用非涡黏性的激励小尺度(Stimulated Small Scale)模式对空间发展的轴对称湍流圆管射流进行了大涡模拟。以雷诺数为10000的流动为例,考证了激励小尺度模式在自由剪切流模拟中的可行性,描述了湍流强度、雷诺应力和湍流耗散量的变化,同时与标准的Smagorinsky涡黏性模式的计算结果进行了比较。数值结果显示,激励小尺度模式能够更为合理地描述湍流的耗散特性和能量传输特性,从而较为准确地展示出空间发展射流中由于流动不稳定而出现的旋涡产生、发展、破碎及合并等过程。  相似文献   

3.
王涛  李平  柏劲松  汪兵  陶钢 《爆炸与冲击》2013,33(5):487-493
采用拉伸涡亚格子尺度应力模型对湍流输运中的亚格子作用项进行模式化处理,发展了适用于可压多介质黏性流动和湍流的大涡模拟方法和代码MVFT(multi-viscous flow and turbulence)。利用MVFT代码对低密度流体界面不稳定性及其诱发的湍流混合问题进行了数值模拟。详细分析了扰动界面的发展,流场中冲击波的传播、相互作用、湍流混合区边界的演化规律,以及流场瞬时密度和湍动能的分布和发展。数值模拟获得的界面演化图像和流场中波系结构与实验结果吻合较好。三维和二维模拟结果的比较显示,两者得到的扰动界面位置、波系及湍流混合区边界基本一致,只是后期的界面构型有所不同,这也正说明湍流具有强三维效应。  相似文献   

4.
均匀各向同性湍流的频率波数能量谱   总被引:1,自引:0,他引:1  
张曙光  雷磊  何国威 《力学学报》2003,35(3):317-320
对均匀各向同性湍流的频率波数能量谱进行了初步研究,得出均匀各向同性湍流中能量随时间、空间尺度分布的解析表达式,并运用直接数值模拟的方法对理论预测进行验证.数值结果很好地支持了理论预测:横扫速度是决定频率波数能量谱的特征量.  相似文献   

5.
饶瑞中 《力学学报》2002,34(5):682-687
建立了一种数值算法,从光波在湍流大气中传播的闪烁时间频谱的高频特征反演大气光学湍流谱.此算法具有两点优越性:(1)反演出的湍流谱直接表征大气折射率的起伏而非大气温度场的起伏;(2)可以获得高空间波数、大部分位于耗散区的湍流谱,而用其它方法不易获得此区域的湍流谱.对Hill湍流谱以及一种复杂的湍流谱进行了大量数值模拟来检验此种算法,结果表明:如果知道大气风速与湍流内尺度就可以相当有效地获得湍流谱.  相似文献   

6.
大涡模拟的壁模型及其应用   总被引:3,自引:0,他引:3  
大涡模拟是研究湍流的非定常特性的重要方法. 但解析壁面层的大涡模拟所需的计算量与直接数值模拟相当,是大涡模拟在高雷诺数壁湍流数值模拟中所面临的主要困难. 解析壁面层所需的网格尺度与壁面黏性长度同量级,是引起壁湍流大涡模拟计算量增加的主要原因. 壁模型通过模化近壁流动避免了完全解析壁面层,可以显著地降低壁湍流大涡模拟的计算量,是克服上述困难的有效方法. 本文介绍了大涡模拟壁模型的主要类型;详细讨论了常用的壁面应力模型,特别是平衡层模型和双层模型的构建思路和特点;基于近壁流动的特征讨论了应力边界条件的必要性和适用性;指出了壁面应力模型的局限性以及考虑非平衡效应修正的各种方法;讨论了壁面应力模型的研究历史、最新进展和发展趋势,给出了常用的壁面应力模型的分支与发展关系图;并基于Werner-Wengle模型实现了周期山状流的大涡模拟.   相似文献   

7.
风口区空气流动的有限元模拟与实验研究   总被引:1,自引:1,他引:0       下载免费PDF全文
通风和空调出风口的射流特性主要取决于出风口截面之前的一段管道内的流动状况,研究送风口的射流偏转特性有重要实用价值.应用湍流大涡模拟技术结合Taylor-Galerkin有限单元法对工程中常见工况下的风口区管道内的三维流场进行了数值模拟和实验研究,数值计算结果与实验结果符合较好.表明湍流大涡模拟方法适合于边界形状复杂,存在各向异性的大尺度涡的内流情况,能可靠地预报风口区空气流动的射流偏转特性.  相似文献   

8.
大涡模拟是研究湍流的非定常特性的重要方法.但解析壁面层的大涡模拟所需的计算量与直接数值模拟相当,是大涡模拟在高雷诺数壁湍流数值模拟中所面临的主要困难.解析壁面层所需的网格尺度与壁面黏性长度同量级,是引起壁湍流大涡模拟计算量增加的主要原因.壁模型通过模化近壁流动避免了完全解析壁面层,可以显著地降低壁湍流大涡模拟的计算量,是克服上述困难的有效方法.本文介绍了大涡模拟壁模型的主要类型;详细讨论了常用的壁面应力模型,特别是平衡层模型和双层模型的构建思路和特点;基于近壁流动的特征讨论了应力边界条件的必要性和适用性;指出了壁面应力模型的局限性以及考虑非平衡效应修正的各种方法;讨论了壁面应力模型的研究历史、最新进展和发展趋势,给出了常用的壁面应力模型的分支与发展关系图;并基于Werner-Wengle模型实现了周期山状流的大涡模拟.  相似文献   

9.
研究二维矩形管道中底部加热的不可压缩Poiseuille-Benard流的谱元法数值计算问题.讨论各种不同的出口边界条件的处理及其对谱元法数值模拟的影响.通过干扰区、干扰幅度和计算时间的比较,确定比较理想的出口边界条件.  相似文献   

10.
采用大涡模拟方法和Smagorinsky亚格子模型,求解三维Navier-Stokes方程,研究了波流边界层中的湍流特性.将大涡模拟结果与相应的直接数值模拟结果和实验数据进行比较,吻合较好.获得了不同波雷诺数,不同波流比情况下的大涡模拟数据库,并由此分析了波流边界层中各种湍流统计量,如速度廓线、剪应力、湍流强度等的变化规律.  相似文献   

11.
采用谱单元方法推导带阻尼梁的传递函数矩阵. 采用一维连续梁的解析解作为动力形状函 数. 与有限元法相比,谱单元方法的自由度和单元数目小且计算精度高. 传递函数表示为梁 的几何和物理参数的超越隐函数,不同于用模态参数表示的传统传递函数. 提出了采用遗传 算法的结构物理参数识别方法. 以变截面悬臂梁为数值算例,显示该方法的有效性.  相似文献   

12.
关于Orr—Sommerfeld方程的Chebyshev谱方法的讨论   总被引:3,自引:0,他引:3  
本文讨论了Orr-Sommereld方程的各种Chebyshev谱离散方法,数值证明了Chebyshev配置法离散Orr-Sommerfeld方程没有伪谱,并以此构造了适于任意平面平行速度剖面情形,对时间和时空稳定性模式一致有效的无伪谱的离离散方法,其中,对时空稳定性问题本文给出了一种新的迭代法可以快速有效地求出复频率的鞍点,对平面Poiseuille流,Blasius边界层流和Gauss模型尾迹  相似文献   

13.
计算响应谱矩的复模态方法   总被引:2,自引:0,他引:2  
方同  张天舒 《力学学报》1989,21(1):65-74
  相似文献   

14.
基于豪斯道夫导数扩散模型的空间谱熵推导 描述反常扩散过程时空复杂程度的空间累积谱熵,并考察 个体谱熵、谱熵、累积谱熵随时 空豪斯道夫导数、扩散系数和扩散时间的变化情况. 计算结果表明,谱熵与累积谱熵随时间豪斯道夫导数α或空间豪斯道夫导数β的减小而增大,且具有拖尾特征. 此外,随着扩散时间t或扩散系数Dα,β的减小,正常扩散对应的个体谱熵衰减的速率比反常扩散快,且对应的谱密度更窄. 因此,豪斯道夫导数扩散模型的谱熵和累积谱熵均能够反映复杂介质的非均质特征和内部扩散过程的不确定性.  相似文献   

15.
透射边界条件在波动谱元模拟中的实现:二维波动   总被引:1,自引:0,他引:1  
邢浩洁  李鸿晶 《力学学报》2017,49(4):894-906
将邢浩洁和李鸿晶提出的多次透射公式(multi-transmitting formula,MTF)的谱元格式应用于均匀介质中线弹性SH波动问题的谱元模拟.假定紧邻人工边界的一层谱单元为具有直线边界的四边形单元,以保证每个人工边界节点都唯一对应一条指向内域的离散网格线.人工边界节点在某时刻的位移由该离散网格线上的节点在前若干时刻的位移确定,按照MTF谱元格式进行计算.通过平面波以一定角度传播的外源问题算例和点源脉冲自由扩散的内源问题算例,验证了方法的可行性以及对实际复杂波动问题的适用性.通过不同类型初值问题算例,在时域内分析了插值多项式阶次、人工波速和透射阶次三个参数对反射误差的影响.结果表明:插值多项式阶次较高的格式会表现出更好的精度,但总体上对反射误差的影响较小;人工波速对反射误差具有显著影响,当人工波速小于介质物理波速时反射误差较大,而当人工波速等于或稍大于介质物理波速时反射误差处于较低水平;透射阶次对反射误差具有决定性影响,表现在不失稳的情形下提高透射阶次能够迅速降低反射误差,但内源问题从三阶MTF开始出现飘移失稳,外源问题从二阶MTF开始出现轻微的飘移失稳.  相似文献   

16.
We extend a large-eddy simulation (LES) methodology, based on using the spectral vanishing viscosity (SVV) method to stabilize spectral collocation approximations, from the Cartesian to the cylindrical geometry. The capabilities of the SVV–LES approach are illustrated for two very different physical problems: (1) the influence of thermal stratification on the wake of a cylinder, and (2) the instabilities that develop in transitional and fully turbulent rotor–stator flows.  相似文献   

17.
将比例边界坐标插值方法引入谱元法, 构成比例边界谱单元, 对无穷域Euler方程进行数值模拟.阐述了比例边界谱单元的基本使用方法以及基于比例边界谱元的Runge-Kutta间断Galerkin方法求解Euler方程的过程;计算了无穷域圆柱和NACA0012翼型绕流问题, 并与已有结果进行了比较, 显示了计算结果的正确性.用基于比例边界谱元的间断Galerkin方法求解无穷域Euler方程时, 最多只需将求解域划分为2个子域, 避免了一般谱方法将求解域划分为9个或者27个子域的麻烦. 比例边界谱单元为无穷域Euler方程的直接求解提供了一个可供参考的方法.  相似文献   

18.
透射边界条件在波动谱元模拟中的实现:一维波动   总被引:1,自引:0,他引:1  
邢浩洁  李鸿晶 《力学学报》2017,49(2):367-379
多次透射公式(multi-transmitting formula,MTF)是一种具有普适性的局部人工边界条件,但其在近场波动数值模拟中一般与有限元法结合.由于波动谱元模拟的数值格式与有限元格式有极大的不同,传统的MTF在谱元离散格式中无法直接实现.为了使物理概念清楚、精度可控的多次透射人工边界条件能够适应波动谱元模拟的需求,首先指出多次透射边界与谱元离散格式结合的基本问题,并分析了空间内插和时间内插两种方案的可行性.然后从空间内插角度出发,提出基于拉格朗日多项式插值模式的MTF谱元格式,并采用一种简单内插方法实现高阶MTF.最后通过一维波动数值试验检验这些MTF谱元格式的精度,并讨论其数值稳定性.结果表明:对于一、二阶MTF,几种格式的精度相当;对于三、四阶MTF,基于谱单元位移模式插值的格式精度最高.相反,随着插值多项式阶次的升高,不同MTF格式的稳定临界值逐步降低,但是所有格式均在人工波速大大超过物理波速时才可能发生失稳.  相似文献   

19.
This paper presents a global Galerkin spectral method for solving the incompressible Navier–Stokes equations in three‐dimensional bounded domains. The method is based on helical‐wave decomposition (HWD), which uses the vector eigenfunctions of the curl operator as orthogonal basis functions. We shall first review the general theory of HWD in an arbitrary simply connected domain, along with some new developments. We then employ the HWD to construct a Galerkin spectral method. The current method innovates the existing HWD‐based spectral method by (a) adding a series of auxiliary fields to the HWD of the velocity field to fulfill the no‐slip boundary condition and to settle the convergence problem of the HWD of the curl fields, and (b) providing a pseudo‐spectral method that utilizes a fast spherical harmonic transform algorithm and Gaussian quadrature to calculate the nonlinear term in the Navier–Stokes equations. The auxiliary fields are uniquely determined by solving the Stokes and Stokes‐like equations under adequate boundary conditions. The implementation of the method under the spherical geometry is presented in detail. Several numerical examples are provided to validate the proposed method. The method can be easily extended to other domains once the helical‐wave bases, which depend only on the geometry of the domains, are available. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

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

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