简介:壳核结构的微胶囊在医学药学材料食品农业等领域具有广泛的应用前景,其制备方法一直是相关领域关注的焦点.同轴流动聚焦(co-flowfocusing)是一种新型制备技术,利用复合射流的破碎制备微胶囊具有包裹率高过程量化可控参数域广产率高等诸多优势.在实验中,复合射流的破碎受到多个过程参数的影响,并涉及了多层界面的耦合效应.利用简化的物理模型,在时间和时空域中分析了三相水-油-水复合射流不稳定性的发展和演化.在黏性流体线性稳定性理论中,同轴射流和驱动液体的基本速度型分别基于管流和误差函数构造,并通过数值方法求解满足相应边界条件下的线化小扰动控制方程.结果表明:增加内外层界面的界面张力均有利于射流的破碎;流体的黏性对同轴射流的稳定性均有着促进作用;越大的黏性越小的内界面张力对应着越大的射流破碎波长;内外界面的耦合作用以及复合液滴的包裹情况均与内外射流的半径比息息相关;绝对-对流不稳定性转换的临界Weber数随Reynolds数内层界面张力的增大而增大,随内层和驱动流体的黏性增大而减小.这些结果将有助于提高液体驱动下同轴流动聚焦技术的过程控制,为实际应用提供理论指导.
简介:为了更加精确地模拟流动/运动耦合问题,建立了耦合动态混合网格生成、非定常流场计算和六自由度运动方程求解的一体化计算方法,并在统一框架内同时实现了松耦合与紧耦合方法.通过圆柱涡致自激振荡(vortexinducedvibration,VIV)的模拟,对不同时间精度的松耦合和紧耦合算法的优劣及适用范围进行了评估和分析;通过引入附加质量的概念,对耦合算法的稳定性进行了理论分析.研究表明:在流体的密度与物体的密度接近时,松耦合方法是不稳定的,必须采用紧耦合方法.最后利用耦合算法对二维鱼体的自主游动和钝锥三自由度自由飞过程进行了数值模拟,证实了理论分析的结论.
简介:流体力学界面不稳定性及其后期的界面混合现象,是一种十分复杂的多尺度非线性物理问题,在惯性约束聚变、天体物理以及水中爆炸等领域有着广泛的应用前景,对该问题的研究不仅具有很高的学术价值,而且对促进相关领域的发展具有重要意义.中国工程物理研究院流体物理研究所基于Euler有限体积方法,发展了适用于可压缩多介质黏性流动具有多亚格子尺度模型的大涡模拟程序MVFT,并评估分析了不同亚格子尺度模型对界面不稳定性及界面混合的模拟能力;提出了流场非均匀性对R-M不稳定性影响的问题,并在激波驱动轻重气体双模扰动R-M界面不稳定性实验中成功应用并解读了新的实验现象和规律,在此基础上进而开展了反射激波作用下两种初始非均匀流场界面不稳定性引起的界面混合数值模拟研究,探讨了流场非均匀性对激波反射后强非线性阶段界面不稳定性发展、演化规律的影响,近期还对非均匀流场R-M不稳定性的演化规律、初始流场非均匀性和初始扰动效应及其影响的物理机制进行了分析和研究.
简介:文章考察了相邻双侧边盖驱动方腔流动(即上壁面向右运动和左侧壁面向下运动)的三维线性整体稳定性.首先,采用Taylor—Hood有限元方法并经由Newton迭代过程计算得到双侧边盖驱动方腔流动的二维稳态基本流.其次,Taylor—Hood有限元在ChebyshevGauss配置点上进行离散,同时Gauss配置点也可以用于线性稳定性方程的高阶有限差分格式离散.然后,离散得到的矩阵形式的广义特征值问题可以结合shift-and—invert算法采用隐式重启Amoldi方法计算.最后,通过对线性稳定性方程特征值的计算,发现了一个最不稳定的驻定模态和两对对称行波模态.最不稳定的三维驻定模态的临界Reynolds数为Ree=261.5,远远小于二维不稳定的临界Revnolds数Ree2d=1061.7.通过画出这3类三维不稳定模态的流向扰动速度和扰动涡量的空间等值面图像,可以发现不稳定扰动位于稳态基本流的两个主涡区域,因此可以认为主涡区域是三维扰动失稳的主要能量来源地.
简介:实验研究复杂波形结构引起平面界面变形和反射激波冲击下的R-M不稳定性的问题.在竖直激波管中生成稳定的N2/SF6平面界面,激波在圆柱绕射后,冲击平面界面,由此研究复杂激波引起的界面变形.平面激波在圆柱绕射后的流场,演化成具有初始入射波、三波点、弯曲反射波、Mach波和Mach反射产生的滑移线等复杂结构.研究复杂结构激波对界面的作用,对认识界面扰动的生成具有较大帮助.绕柱激波冲击后,平面界面仅在两对滑移线内部发生变形.绕柱激波冲击界面后,两对滑移线将界面分成"内界面"和"外界面",界面变形形态同滑移线和界面相交位置相关.反射激波二次冲击下,界面扰动的增长与Jacobs-Sheeley涡量模型较吻合.
简介:为了更好地理解不同空间坐标系下流体界面对Rayleigh—Taylor(RT)不稳定性弱非线性阶段谐波的影响,文章采用3阶小扰动展开法,解析研究了球坐标空间经典RT不稳定性弱非线性阶段谐波的演化规律,并和柱坐标空间以及直角坐标空间相应结果进行了对比研究.当球坐标系和直角坐标系中RT不稳定性界面扰动波长相同,球坐标系中初始扰动半径为无穷大时(即球坐标下RT不稳定性初始扰动半径相对于扰动波长为无穷大时),球坐标下RT不稳定性前4次谐波的结果和直角坐标系下的相应结果相同.研究表明:由初始界面曲率引起的Bell-Plesset(BP)效应和空间效应(直角坐标空间、柱坐标空间和球坐标空间)对谐波发展有较大的影响.即在不同正交曲线坐标系下,不同曲率的流体界面效应对RT不稳定性谐波发展有较大的影响.对于柱坐标空间和球坐标空间,2阶对0次谐波的反馈加强了界面向内收缩.研究还表明:界面效应增加了2次谐波的负反馈,然而,对于基模和3次谐波却有不同的影响.
简介:采用低耗散WENO(weightedessentialnon-oscillatory)格式及锐界面方法模拟可压缩Kelvin-Helmholtz不稳定性问题.由于物质界面被描述成一种接触间断,该方法可精确求解切向速度间断.基于优化模板对原始光滑指标进行正规化后,得到一种低耗散WENO格式.修正后的方法显著降低了普通流动区域的过衰减问题,保持了良好的激波捕捉性能,并可获得与混合格式相当的求解精度.不同于以往求解单一流体或易混界面时,通过初始设定有限宽度的剪切层或快速数值耗散以抑制高波数模态,该方法允许高波数扰动的发展.计算结果表明,高波数扰动展现出与以往理想Kelvin-Helmholtz不稳定性问题数值模拟或线化理论结果不同的特征,但与有限厚度的剪切层结果相符.
简介:开式凹腔作为超燃冲压发动机中增加掺混和稳焰的装置,其流动稳定性的研究对深入理解凹腔增加掺混和稳焰机理以及凹腔的设计有着重要的学术意义和工程应用价值.基于大涡模拟方法对超燃冲压发动机开式凹腔流动进行数值模拟,分别米用动力学模态分解(dynamicmodedecomposition,DMD)和本征正交分解方法(properorthogonaldecomposition,POD)对自激振荡流动进行稳定性分析.DMD方法可准确提取凹腔的振荡频率,与Rossitei'模型以及压力脉动FFT分析得到的频率吻合较好,且DMD中对应Roster前3阶频率的模态在流动中的主导作用顺序也与FFT分析结果一致,自激振荡中RossiterH模态占据主导作用,同时DMD方法对Rossiter3阶以上模态频率的预测能力明显强于FFT分析方法.在对低频的提取方面,DMD方法比Rossiter模型更具有优势.与前6阶Rossiter模态对应DMD模态均缓慢收敛,主要表现为剪切层中的分离涡结构和中部及下游区域中的涡结构.前3阶不稳定模态中的分离涡结构主要集中在中部剪切层以及后缘附近区域.POD方法中较少的模态包含流场绝大部分的能量.但是,通过POD方法提取的模态频率在分辨率上效果不佳,提取到最低频率为Rossiter3阶模态对应的频率,且模态中均存在次频,次频与主频之间的耦合导致模态的形态相差较大.另外,与DMD方法相比POD方法无法判断所提取的模态的稳定性.
简介:文章首先采用单相浮阻力模型对不同加速度下Rayleigh-Taylor不稳定性诱发的物质渗透边界的演化过程进行了计算,揭示了该混合在常加速度和变加速度情况下不同的发展规律,并通过与实验结果的比较分析,验证了该模型的适用性.在此基础上,发展了多相浮阻力模型,采用该模型对常加速度情况下含尘气体中的RayleighTaylor不稳定性诱导混合进行了研究,发现混合区宽度随着颗粒数密度和颗粒尺寸的增大而减小,揭示了气体中所含杂质抑制混合发展的规律.
简介:传统的应用稳定性理论对横流不稳定性转捩现象的预测很难与现代CFD并行化计算结合,为了解决这个问题,文章基于SA-γ-—Re_(θt)转捩模型,使用不可压三维边界层相似性解实现横流位移损失厚度Reynolds数在流场中的当地化求解,结合亚音速试验数据-C1准则构建横流不稳定性转捩判据,从而实现了横流不稳定性转捩预测方法的当地化并行求解.首先采用SA-γ-—Re_(θt)转捩模型对NLF(1)-0416翼型进行了流向转捩预测,证实了该模型的正确性.然后应用所建立的横流转捩模型对45°前缘后掠角的NLF(2)-0415无限展长机翼和DLR-F5机翼,以及标准6∶1椭球标模进行了横流不稳定转捩数值模拟,计算结果显示转捩位置均与试验数据吻合较好,证明了文章所建立的方法在不可压边界层转捩预测具有较高的预测精度.