考虑叶轮流动滑移的离心泵盘腔流动模型

为研究离心泵盘腔内外特性，将离心泵盘腔三维流动简化为动静腔一维流动，基于动量守恒方程和径向连续方程，构建了动静腔流动微分方程；积分该方程并应用到盘腔中，引入密封件阻力系数，考虑叶轮流动滑移，修正泵势扬程公式，建立了离心泵叶轮-盘腔-密封环-平衡孔回路的自封闭流动模型。提出了内外迭代相结合的求解方法，实现了离心泵盘腔流动快速计算。通过与盘腔压力测量结果对比表明：使用Stodola叶轮流动滑移系数的盘腔流动模型比使用Wiesner叶轮流动滑移系数的盘腔流动模型的计算精度提高了约3.5%。盘腔流动模型结果表明：泵工况从小流量调节至大流量，容积效率逐渐提高，盖板推力略微逐渐增大。该方法及结论为离心泵盘腔内外特性计算提供参考。     

跨声速风洞斜孔壁非线性流动试验

为分析跨声速风洞斜孔壁近壁区域的流动特性，评估气流偏角-压力系数的非线性关联，在0.6 m跨超声速风洞中开展了基于七孔探针的流动特性测量试验。通过气流偏角和压力系数分布分析了斜孔壁流动的差阻特性，以及马赫数、模型升力对斜孔壁流动的影响，最后基于试验结果发展了计算斜孔壁特性参数的微分法，并与经验方法结果进行对比。结果表明，斜孔壁流动呈现出明显的差阻性和非线性，在负压差范围内，近壁流动仍以出流为主；高亚声速时，空风洞模型区孔壁流动特性趋于实壁；安装模型后，随着升力的增大，升力面对应的孔壁区域流动向入流发展，孔壁流动特性趋于开口边界。     

劳动力流动动力问题分析

通过对劳动力要素流动问题的分析，试图找到劳动力流动原因以及流动意义的理论解释，从而将这一理论应用于实践中，为国际劳动力流动及区域劳动力流动的政策实践提出建议。     

运载火箭注气式蓄压器联通孔流量特性数值仿真研究

运载火箭注气式蓄压器存在主流垂直流动对联通孔内水平流动的影响，如果采用一维经典流体力学方法计算联通孔流量系数会产生很大误差。针对这一问题，采用计算流体力学方法，在不同主管流速的条件下，对单个联通孔内的流体流动进行三维稳态数值仿真研究。研究表明，联通孔流动同时具有T型三通管流动和孔板流动的特点。受主流垂直流动的影响，联通孔内流动呈现出明显的不对称性和非线性特征，流速比小于-2.54或大于1.29时的流量系数相对变化超过15%。根据仿真结果，给出了考虑主管流速的联通孔流入和流出流量系数拟合式，以及流量系数相对变化量和流速比的关系式。     

地效飞行器大迎角近地面飞行分离流动数值研究

 本文利用Fluent软件求解定常可压缩流动的质量加权平均N-S方程和标准 湍流模型，数值模拟地效飞行器大迎角近地面飞行的全机流场，研究地面效应引起的分离流动。研究结果给出了大迎角近地面飞行机翼上翼面分离流动随飞行高度的变化规律，分析了飞行高度降低引发流动分离的原因，阐述了流动分离对气动性能的影响规律。     

民机废水系统气-液-固三相流流动特性仿真计算研究

民机废水系统是保障乘员舒适性和个人卫生的重要系统，而废水系统管路气液固三相流动特性分析是其关键技术难点之一。通过对民用飞机废水系统进行研究，采用CFD-DEM模型对马桶底部和废水输送管路中的气液固三相流动过程进行仿真计算，对废水系统管路的气液固三相的流动形态和流动特性进行分析，得到了液相和固相的流动时间。结果表明：废水管路中固相和液相沿运动方向呈贴壁螺旋流动，并且气相速度大的区域湍动能越大，因此对于气液固三相流动而言，气相的运动占主导；截面平均速度气相最大，液相次之，固相最小；系统静压沿管路轴向不断下降，总压损失随时间不断减小。通过对民机废水系统流动性能的分析，为废水系统设计提供设计输入，具有指导意义。     

节流通道振动对碳氢燃料不稳定流动的影响研究

通过实验发现振动的节流通道会引起碳氢燃料的不稳定流动现象。为了探究节流通道振动引发碳氢燃料不稳定流动的机理，采用数值求解三维Reynold-Averaged Navier-Stokes（RANS）方程和k-ω湍流模型的方法研究了节流通道振动对流场动态特性的影响，重点关注节流通道振动对扩张结构附近流动分离的动态特性的影响。数值模拟节流结构两端压差与实验数据对比，验证了所采用的数值方法和湍流模型的有效性。结果表明，当燃油流动方向和节流通道振动方向一致时，会发生不稳定流动现象。当燃油流动方向和节流通道振动方向垂直时，不会发生不稳定流动现象。通过流场的动态分析，发现振动方向与流动方向一致时扩张结构附近非定常涡会发生轴向的周期性运动，从而导致了不稳定流动的现象。振动速度越大，不稳定流动程度越大。     

复合材料固化过程树脂流动模型

本文以Ｄａｒｃｙ定律为基础，以复合材料固化过程中的树脂流动为条件，推导和建立了树流动模型。该模型形式简单，易于应用于实际。推导结果表明：树脂流动仅同三个因有关，即树脂粘度，加压压力和时间。特别是树脂流动质量与时间有简单的玉方根关系。另外，树脂粘度是影响流动的唯一物质参数，也是影响流动的最重要的因素。模型可用来讨算脂流动速度和质量。同粘度关系式联系起来，用以选择加压时间和加压大小。有关实验证实，模型     

轮背凹槽结构对半开式向心涡轮性能影响数值模拟

为了厘清半开式向心涡轮转子轮背间隙内的流动机理及轮背间隙泄漏流动对涡轮性能的影响，以某半开式向心涡轮为 研究对象，采用数值模拟的方法对轮背间隙泄漏流动进行分析，深入剖析了封严气、泄漏流及主流的运动轨迹，并明确了轮背泄漏损失分布。针对轮背间隙内泄漏流动的特点，对轮背结构进行了改进设计，并研究了轮背凹槽的形状、长度、宽度和深度对泄漏损失的影响。结果表明：半开式向心涡轮轮背存在横向泄漏流动，并与主流和封严气进行掺混，造成流动损失，恶化了涡轮性能，采用轮背凹槽结构后，可有效抑制轮背的泄漏流，降低流动损失。凹槽的形状、长度、宽度和深度对泄漏流的影响程度不同，椭圆形凹槽可使涡轮效率提高0.2%，矩形凹槽可使涡轮效率最大提高0.23%。     

低速风洞设计中的导流片损失问题

低速回流式风洞中拐角部位的流动损失占总损失的比重很大，该位置导流片绕流的流动状况对于流动损失有明显的影响。与常规的工程估算方法不同，风洞实验与数值计算相结合，并通过正交优化能够对不同工作条件下导流片流动及其性能指标的变化规律进行研究，对特定的导流片找到了最佳的工作条件。     

球柱三维分离及涡流的数值模拟

对三维可压非定常Ｎ－Ｓ方程进行有限体积中心离散,并用Ｒｕｎｇｅ－Ｋｕｔａ显式格式进行数值求解。数值模拟了球柱分离及涡流的形态结构,给出了具有激波、边界层、大分离区、涡流及其间相互干扰的复杂流型     

MLP高阶重构格式应用

在有限体积框架下,利用MLP(multi-dimensional limiting process)系列重构格式结合HLLHLLC(Harten-Lax-van Leer with contact)近似黎曼求解器,同时引入激波探测函数进一步降低MLP在光滑流动区域的数值耗散,数值模拟了超声速前台阶流动、结冰翼型的非定常流动、高超声速双楔流动和DLR F6-WB跨声速流场,研究了MLP系列格式在可压缩复杂流场中的表现.结果表明:在多维空间中,MLP格式能够在如强斜激波与网格线不重合等复杂流场数值模拟中保持严格的流场单调性;具有和传统MUSCL(monotone upstream-centered schemes for conservation laws)格式类似的计算效率,可以实现5阶,甚至更高阶重构;数值耗散更低,捕获更准确的激波位置,对航空工程数值模拟具有重要意义.     

运载火箭复杂管路随机性压力跳变机理及抑制

针对3种典型的管路随机性压力跳变现象，通过仿真分析，提出管路中的直角多通结构使管内流动具有非线性的流动多解特性，在特定扰动条件下，会产生多通内的大涡结构及沿分支管的螺旋流动，引起额外的流阻，并导致管路局部压力发生跳变，从机理上复现了复杂管路内的随机性压力跳变现象；在此基础上，提出了增加导流结构/隔板结构抑制多通内大涡流动的改进方案，并针对真实管路结构，开展了管路压力跳变现象的地面复现及抑制改进试验，试验结果表明:3种直角多通结构均存在流动多解特性，增加导流结构可有效消除管路中的随机性压力跳变现象。      

NND格式在多维理想磁流体方程组中的应用

采用修正的四步Runge-Kutta方法求解三维一般曲线坐标系下的理想磁流体方程组,为克服数值振荡,加特征型NND格式进行后处理.特征型NND格式推广到求解三维磁流体(MHD)问题需要知道雅可比通量的左右特征矩阵,在具体计算时需要克服矩阵的奇性.本文用三维程序采用推广的特征NND格式计算了一维MHD激波管和二维(MHD)喷管流动,计算结果表明,特征NND格式保持了TVD格式高精度的优点,又具有计算简单的特点,在包括强弱间断等复杂波系的定常和非定常MHD流场数值模拟中是成功的.     

DCS5在计算气动声学中的应用研究

对流项采用五阶线性耗散紧致格式(DCS5)离散,粘性项采用6阶中心紧致格式离散,时间项采用5级5阶龙格-库塔时间推进格式求解流体动力学方程组,并由此建立近声场的直接模拟方法。为了应用结构网格在复杂几何形状条件下对流场进行计算,发展了高精度的多块特征对接边界技术。完成了低速无粘振荡圆柱和粘性圆柱绕流等标准计算声学算例,获得了令人满意的结果。     

无波动,无自由参数的耗散差分格式（NND）在喷流计算中的应用

本文用文献[1]建立的无波动、无自由参数的耗散差分格式(NND格式),对无粘喷流流动进行了数值模拟,并就一典型的自由射流流动,与国外高精度、高分辨率的PLM差分格式(Piecewise-Linear Method)的计算结果和实验结果作了比较,得到了很一致的结果,且所用的计算时间少于PLM格式。     

A hybrid multidimensional Riemann solver to couple self-similar method with MULTV method for complex flows

Since proposed, the self-similarity variables based genuinely multidimensional Riemann solver is attracting more attentions due to its high resolution in multidimensional complex flows. However, it needs numerous logical operations in supersonic cases, which limit the method’s applicability in engineering problems greatly. In order to overcome this defect, a hybrid multidimensional Riemann solver, called HMTHS (Hybrid of MulTv and multidimensional HLL scheme based on Self-similar structures), is proposed. It simulates the strongly interacting zone by adopting the MHLLES (Multidimensional Harten-Lax-van Leer-Eifeldt scheme based on Self-similar structures) scheme at subsonic speeds, which is with a high resolution by considering the second moment in the similarity variables. Also, it adopts the MULTV (Multidimensional Toro and Vasquez) scheme, which is with a high resolution in capturing discontinuities, to simulate the flux at supersonic speeds. Systematic numerical experiments, including both one-dimensional cases and two-dimensional cases, are conducted. One-dimensional moving contact discontinuity case and sod shock tube case suggest that HMTHS can accurately capture one-dimensional expansion waves, shock waves, and linear contact discontinuities. Two-dimensional cases, such as the double Mach reflection case, the supersonic shock / boundary layer interaction case, the hypersonic flow over the cylinder case, and the hypersonic viscous flow over the double-ellipsoid case, indicate that the HMTHS scheme is with a high resolution in simulating multidimensional complex flows. Therefore, it is promising to be widely applied in both scholar and engineering areas.     

再人体复杂外形无粘流场数值计算研究

本文在文[1,2,3]的基础上研究了非定常Euler方程的推进迭代方法,并将无波动、无自由参数的耗散差分格式(下称NND格式)发展为隐式、迎风格式,用以计算有差曲控制翼的再入复杂飞行器的育攻角、有侧滑角的超声速无粘绕流流场。其结果是准确的。本文的方法具有以下特点: 1.可以在同一计算程序中实现亚跨超流场的计算:在超声速区域,采用推进技术;在亚跨声速区,则自动选择推进迭代求解技术。 2.NND格式不仅具有公式简单,无自由参数的特点,而且过激波时满足熵条件,不会产生非物理解,也不发生振荡。 3.隐式差分格式的无条件稳定和Gauss-Seidel迭代的快速收敛,计算省时。 本文的方法曾用于烧蚀头部凹陷外形和航天飞机简化外形的计算,能自动选择迭代和推进区域。可以预言,如计算机条件许可,本方法可以用于大攻角的无粘流场计算。     

再入弹抛壳气动特性研究

本文用有限差分方法求解再入母弹头与其所抛壳体之间的无粘干扰流场,获得了体—壳的干扰气动力,并在准定常假设的基础上近似模拟了弹体与壳体的分离过程。数值研究的结果表明,体—壳间存在的干扰会导致壳的压心位置急剧变化,危及母弹安全;抛壳速度是体—壳安全分离的关键参数。     

气体运动论BGK格式的翼型绕流数值模拟

以气体运动论BGK(Bhatnagaar-Gross-Krook)格式为基础,通过改进非均匀网格下的通量计算方法,把气体运动论BGK方法应用到二维贴体网格上模拟了可压缩的无黏和黏性流动。在重构阶段,以交界面周围4个网格单元为一组,先把每个网格单元的宏观量(密度、动量和总能量)转换在相对于交界面的坐标系下,然后使用van Leer限制器插值求出交界面的宏观量,在求解阶段再使用BGK理论求出流过交界面的通量。通过上述方法使贴体网格和均匀网格,在计算格式上一致,与其他求解二维贴体网格的气体运动论BGK方法相比,该方法相对简单,计算量减小。结合贴体网格通量计算和气体运动论BGK格式方法理论,分别对亚、跨、超声速下的NACA0012翼型绕流进行了数值模拟和对比研究,所得可压缩无黏流动的数值结果与北大西洋公约组织航天研究与发展咨询组(AGARD)的结果吻合良好。本文还模拟了低雷诺数状态下的可压缩黏性流动,所得结果与文献中的数值计算结果吻合良好,证明了该方法在数值模拟及激波捕捉方面具有鲁棒性,为在复杂流动中的应用提供了一种简单的途径。