黏性流场中圆柱壳结构的振动能量流输入特性

以周向分布线力作用下,浸没在黏性流场中的无限长弹性薄圆柱壳为研究对象,分析流体黏性对黏性流场-圆柱壳耦合系统的振动能量流输入的影响规律,考虑周向分布线力作用下,薄圆柱壳壳体与由有黏、可压缩流体组成的流场之间的相互作用.首先采用Flügge薄壳理论分析壳体结构振动,然后把线性化的连续性方程、线性处理后的纳维-斯托克斯方程和小振幅波动下的状态方程结合起来得到了黏性流场中的声波波动方程,进而根据声场与圆柱壳外表面的运动协调条件推导得到耦合系统的声振耦合方程,提出了相应的数值计算方法,并与理想流体中的特性作了对比,研究了流体黏性对黏性流场-圆柱壳耦合系统的振动能量流输入的影响规律,为水下结构的减震降噪提供理论依据.     

具有运动边界不可压缩黏性流动的CBS有限元解法

基于特征分裂(CBS)算法,提出了一种适用于具有运动边界的不可压缩黏性流动问题的有限元解法,给出了不可压缩黏性流动问题的控制方程,并推导了其无量纲形式.将传统CBS算法与动网格技术相结合,提出了一种动边界流动问题的求解流程,并对所采用的基于弹簧近似的动网格方法进行了说明.采用该方法对静水和均匀来流中的简谐振动圆柱绕流问题进行了求解,计算所得圆柱受力、圆柱表面压力和涡量分布以及流场速度分布等结果均与已有数值或实验结果相吻合,初步验证了所提方法的可行性.     

基函数法在黏性可压缩流动中的研究与应用

研究基函数法在二维和三维轴对称情形下可压缩黏性流动中的应用, 并利用此法数值地计算了二维圆柱的超音速黏性绕流和三维球头的轴对称超音速黏性绕流两个算例。取一阶三角函数为基函数, 构造出导数的中心格式和迎风格式。对于可压缩N-S 方程中的对流项, 采用通量分裂法及中心格式与迎风格式相结合的技术;对于可压缩N-S方程中的黏性项, 则采用中心格式进行处理。由此, 构造出了数值求解黏性可压缩流动的一阶三角函数类型的基函数格式。两个算例的数值计算结果表明, 基函数法不仅在处理无黏可压缩流动时, 是一种高精度、高分辨率的新型计算方法, 而且在处理黏性可压缩流动问题时也一样行之有效。     

弱可压缩SPH方法压力初始赋值方法比较

为研究不同初始流场压力赋值方法对于弱可压缩SPH方法数值结果稳定性的影响,采用静水压力赋值方法以及弱可压缩假设赋值方法,结合不同物理模型及不同数值影响因素对弱可压缩SPH方法的初始流场压力稳定性问题进行对比分析.对粒子数的比较研究表明,相同粒子数下,弱可压缩法在流场稳定速率及流场稳定性上均优于静水压力法.对人工声速的研究表明,弱可压缩法的自身稳定性及对人工声速的不敏感性显著优于静水压力法.对人工黏性的比较表明,静水压力法较弱可压缩法对人工黏性系数的敏感性较高,易受人工黏性变化而产生压力波动.弱可压缩假设赋值方法具有良好的鲁棒性,有利于提高初始流场收敛速度,提升流场数值稳定性,对于弱可压缩SPH方法流场初始赋值方法选择具有一定意义.     

双拉格朗日模型模拟气固两相双圆柱绕流

以湍流剪切层中大尺度旋涡结构与颗粒群的相互作用及运动规律为研究对象,提出了一种基于离散涡方法与颗粒碰撞模型相耦合的双拉格朗日气固两相流模型.用次循环调整颗粒运动计算的时间步长,将两相计算耦合在一起,从而解决了离散涡方法求解单相流场与颗粒运动方程的时间相容性问题.首先用颗粒在兰金涡中的运动验证了所提的算法,随后将该方法应用于高雷诺数湍流旋涡结构中的颗粒水平输运过程研究,模拟了圆柱尾迹中颗粒的沉降过程及其被旋涡卷吸、扬起和团聚的过程,说明在有大尺度旋涡存在的流场中颗粒更容易被流体向下游榆运.当来流速度与重力方向一致时,颗粒分布特性受颗粒Stokes数(St)影响.当St小于1时,颗粒响应流场变化的时间短,随St增加颗粒对流场变化的响应变得迟钝,重力所起的作用更加明显;当St远远大于1时,颗粒的分布几乎不受旋涡的影响,而与自由沉降的分布十分接近.     

列车进入隧道时产生涡流的数值模拟

给出了列车穿越隧道时在隧道入口形成涡流场的数值模拟过程．控制方程为三维黏性、可压缩、等熵、非定常流的N-S（Navier—Stokes）方程,空间离散采用了中心有限体积法格式,时间采用预处理二阶精度多步后差分格式进行离散,对列车与隧道之间的相对移动采用移动网格技术处理．计算结果与国外的试验结果基本一致．研究结果表明,在整个涡流的形成过程中,可以将这个过程分为连续的4个阶段,即涡的形成、发展、传播及破坏阶段．研究结果还表明,涡受到列车头部形状的影响要比受到列车速度的影响要大;喷射流速度与列车速度有关系．     

黏性流中的链环涡管演化与螺旋度分析

本文通过直接数值模拟研究了链环涡管在不可压缩黏性流中的演化过程.在初始时刻的链环涡管由两个变形的涡环组合而成,其螺旋度为依赖于涡轴参数方程的解析表达式,进而可利用该初始流场进行涡管演化研究与螺旋度分析.发现当初始链环涡管的涡量方向具有相同手性时,链环涡管和环面纽结涡管的演化具有类似的涡动力学过程;而当它们具有相反手性时,涡管间的强涡量梯度会使两涡环在短时间内产生剧烈的涡重联,从而导致涡环由快速的尺度级串过程达到类湍流状态.     

带哨嘴喷嘴射流流场的离散涡模拟

采用离散涡方法模拟了两种不同形状哨嘴的喷嘴出口流场旋涡结构和压力分布。模拟结果表明 ,射流在喷嘴出口和哨嘴出口形成的旋涡结构在流场及其他涡元的作用下 ,向下游运动 ,在流场中某个地方形成低压区 ,这有助于诱发空化现象。对两种不同喷嘴的出口流场进行比较发现 ,哨嘴形状对形成的射流流场的旋涡结构有很大影响。因此 ,研究和应用空化射流必须考虑喷嘴出口哨嘴形状的影响     

结构网格下求解非规则区域流场的一种差分法

对于粘性不可压缩流场的计算，提出一种标识边界的差分方法，对不规则区域的流场边界做标记，应用SIMPLE算法在直角坐标系的结构网格下离散求解N-s方程；对圆柱绕流，不规则腔体流等问题进行数值模拟，取得了较好的效果。     

着火物在火焰旋涡流场中运动规律研究

在森林火灾中,各种旋涡大量存在,着火物在旋转流场作用下运动导致火灾蔓延．本文采用二相流单颗粒模型,研究了着火物在火焰旋涡流场中的运动,给出了环境气流分别以Ｂｕｒｇｅｒｓ涡和Ｓｕｌｉｖａｎ涡作旋转运动时,着火物的运动规律．结果表明,着火物将在火焰旋涡流场作用下,被加速带到高空,并飞离流场,形成飞火传播．     

二维浅水流动的涡运动学特性及动力学特征

建立了二维浅水流动的旋涡诱导速度公式，在已知旋涡场和胀量场时，由该表达式可确定流动速度场，对比分析了二维浅水流动中的涡运动与可压缩流动中涡运动的异同点。     

用LBM方法模拟壁面驱动粘性不可压半圆形空腔流

基于格子Boltzmann方法(LBM)数值模拟壁面驱动的粘性不可压半圆形空腔流. 采用具有二阶精度的曲线边界处理方法, 得到了不同雷诺数下的流线图、 涡线图及速度分量沿半圆形中心线的分布. 在小雷诺数的条件下, 流动状态仅由一个涡组成; 随着雷诺数的增加, 出现一个二级涡, 涡的大小与雷诺数有关. 数值结果表明, 格子Boltzmann方法简单有效, 适合处理该问题.     

超临界雷诺数下圆柱分离旋涡运动的数值研究

本文在文献[1]的基础上,用离散涡模型与湍流边界层理论相结合的方法,研究了超临界雷诺数下圆柱突然起动后的分离旋涡运动。文中考虑了二次涡的影响和涡旋的粘性扩散。计算结果与实验结果相符。     

用离散涡方法模拟喷嘴出口射流流场的数值模型

根据流体力学原理,用离散涡方法建立了在复平面上模拟喷嘴出口流场的一套完整模型。在此基础上,用ＦＯＲＴＲＡＮ语言编写了计算机程序,模拟了喷嘴出口流场的速度分布,利用涡元矢量图对计算结果进行了分析,得出了射流在喷嘴出口分离流动及流场中旋涡形成及发展的规律。射流在有锐缘的地方会发生流动分离,形成的剪切边界层逐渐卷起,形成旋涡结构。在喷嘴出口与哨嘴出口处形成的旋涡结构,在势流流场和其它涡元的作用下逐渐向下游运动,并向两边扩散。由于粘性扩散及计算误差的影响,在远喷距处会出现涡旋结构变形     

圆柱绕流远场涡对的数值模拟

基于离散涡方法求得非定常、不稳定流场,数值模拟了三种不同时刻高雷诺数下圆柱绕流结构的发展,从流谱图、等涡量线图和涡谱图可以清晰地看出从近场的初生卡门涡街,过渡到远场的二次涡街的过程,计算结果发现：远场离散涡有形成二个涡的涡对及三个涡的涡对的趋势,计算结果说明了流体运动中涡对结构的本质：由于来流是均匀的,没有加入任何拔动,当流体流过钝体时产生具有剧烈分离的不稳定流动,因此在远场形成的二次涡对及卡门涡     

摇荡物体周围粘性流的数值模拟

提出了一种新的求解非定常不可压Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的数值方法，在每个时间步长上采用人工可压性技术，通过对动量方程取散度，从而获得离散化的Ｐｏｉｓｓｏｎ方程，当关于虚拟时间的迭代过程达到收敛时，便得到了无散速率场，采用该方法模拟了一个纵摇椭圆柱体周围的非定常不可压性流动。     

单圆柱微通道内流动特性实验研究

采用Micro-PIV实验系统和压差测试系统,研究了含有单个微圆柱的通道内去离子水在10相似文献   

飞行器头部的涡系结构研究

通过三阶精度离散格式求解N-S方程,模拟了飞行器头部涡系流动过程,展示了不同层次的涡系结构.结果表明,较小攻角时,表面不同的鞍点与结点拓扑组合,在飞行器头部空间演化出不同的涡系.稳定螺旋点在空间演化为羊角涡,促使对称主涡向非对称发展.较大攻角时,头部空间演化出U型马蹄涡,通过牵制主涡涡核,抑制主涡向非对称发展.中等攻角时,飞行器头部既形成羊角涡,又形成马蹄涡,二者对主涡的作用是相反的.     

Generalization of Hopf Functional Equation

IntroductionThe techniques ofcharacteristic functionals of fluidmechanical fields are widely used in solving thegeneral problem of turbulence. Thesecharacteristic functionals uniquely define theprobability distributions P(dω) in the phase spaceof turbulent flow,which would give a completesolution of the problem of turbulence.Thetechniques of using a characteristic functional todefine the random field of any fluid dynamicquantity were first pointed out by Kolmogorov(1 93 5 ) ,since then a num…