具有不同渗透的涨缩管道内微极性流动的同伦分析解

分析了壁面具有不同渗透的涨缩管道内微极性流体的流动．对于壁面的胀缩,考虑常系数和时间函数的膨胀率两种情况．对于第1种情况,应用同伦分析方法得到该问题的速度和微旋转角度的表达式．并且画图分析了各个不同参数,特别是膨胀系数和不同的渗透率对流体的动力特征的影响．可以得到第1个重要的结论:壁面的膨胀率和不同的渗透对流体的动力特征有重要的影响．根据Xu的模型,考虑了第2种也是更具有一般性的情况,假设壁面的膨胀率随时间的变化而变化．在这样的假设下,控制方程被转化成非线性偏微分方程,并且同样也可以应用HAM方法进行求解．应用代数和指数的模型来描述膨胀率从初始状态到最终状态的演变过程．然而,结果表明包含有时间的解很快地趋向于稳态的解．这样可以得到第2个重要的结论,时间在壁面的膨胀收缩中扮演着次要的角色,可以忽略不计．     

射频辉光放电CH4等离子体一维流体动力学模拟

完整建立一个关于射频辉光放电CH₄等离子体的流体动力学模型.模型包括基于迁移-扩散近似的粒子平衡方程、电子能量平衡方程,共包含了20种粒子(环境气体粒子,激发态粒子,离子和电子)和49类化学反应(电子-中性环境粒子、离子-中性环境粒子、激发态粒子-激发态粒子(中性环境粒子)).结果表明,在强电场区域有较高的电子反应率系数;等离子体中除源气体CH₄外,H₂,C₂H₆,C₃H₈,C₂H₄和C₂H₂也有较高的密度含量;激发态粒子中,CH₃含量最高,密度约为10¹⁹m^-3;在较低放电压力时(如18 Pa),CH₅⁺在离子成分中密度含量最高,当放电压力较高时(如67Pa),C₂H₅⁺在离子成分中占主导地位;除C₂H₅⁺外,其它各离子和激发态粒子在极板上的粒子流量随功率的增大逐渐升高.     

对流扩散方程QUICK格式的数值摄动高精度重构格式

利用高智提出的数值摄动算法, 把对流扩散方程的常用QUICK格式(黏性和对流项分别用二阶中心和QUICK格式离散)进行了高精度重构, 包括利用离散单元内所有结点的全域重构和分别利用离散单元内上下游结点的上下游重构, 得到两类新的更高阶精度的数值摄动重构格式, 称为高的QUICK格式(G-QUICK格式). G-QUICK格式与QUICK格式相比简单性相当, 但精度更高; 全域重构G-QUICK格式和QUICK格式均为条件稳定, 上下游重构得到一些绝对稳定的G-QUICK格式. 解析分析和数值算例均证实了G-QUICK格式的优良性能, 上下游重构的G-QUICK格式为在对流扩散方程的QUICK格式中避免使用人工黏性提供了新途径.      

Melting heat transfer effects on stagnation point flow of micropolar fluid saturated in porous medium with internal heat generation （absorption）

The effect of melting heat transfer on the two dimensional boundary layer flow of a micropolar fluid near a stagnation point embedded in a porous medium in the presence of internal heat generation/absorption is investigated. The governing non-linear partial differential equations describing the problem are reduced to a system of non-linear ordinary differential equations using similarity transformations solved numerically using the Chebyshev spectral method. Numerical results for velocity, angular velocity and temperature profiles are shown graphically and discussed for different values of the inverse Darcy number, the heat generation/absorption parameter, and the melting parameter. The effects of the pertinent parameters on the local skin-friction coefficient, the wall couple stress, and the local Nusselt number are tabulated and discussed. The results show that the inverse Darcy number has the effect of enhancing both velocity and temperature and suppressing angular velocity. It is also found that the local skin-friction coefficient decreases, while the local Nusselt number increases as the melting parameter increases.     

Dual solutions in boundary layer flow of a moving fluid over a moving permeable surface in presence of prescribed surface temperature and thermal radiation

An analysis of the heat transfer for a boundary layer forced convective flow past a moving permeable flat surface parallel to a moving fluid is presented. Prescribed surface temperature at the boundary is considered. A thermal radiation term in the energy equation is considered. The similarity solutions for the problem are obtained and the reduced ordinary differential equations are solved numerically. To support the validity of the numerical results, a comparison is made with the available results for some particular cases of this study. Dual solutions exist when the surface and the fluid move in the opposite directions.     

Dual solutions in MHD stagnation-point flow of Prandtl fluid impinging on shrinking sheet

The present article investigates the dual nature of the solution of the magneto- hydrodynamic （MHD） stagnation-point flow of a Prandtl fluid model towards a shrinking surface. The self-similar nonlinear ordinary differential equations are solved numerically by the shooting： method. It is found that the dual solutions of the flow exist for cer- tain values of tile velocity ratio parameter. The special case of the first branch solutions （the classical Newtonian fluid model） is compared with the present numerical results of stretching flow. The results are found to be in good agreement. It is also shown that the boundary layer thickness for the second solution is thicker than that for the first solution.     

Flux vector splitting solutions for coupling hydraulic transient of gas-liquid-solid three-phase flow in pipelines

The gas-liquid-solid three-phase mixed flow is the most general in multiphase mixed transportation. It is significant to exactly solve the coupling hydraulic transient problems of this type of multiphase mixed flow in pipelines. Presently, the method of characteristics is widely used to solve classical hydraulic transient problems. However, when it is used to solve coupling hydraulic transient problems, excessive interpolation errors may be introduced into the results due to unavoidable multiwave interpolated calculations. To deal with the problem, a finite difference scheme based on the Steger-Warming flux vector splitting is proposed. A flux vector splitting scheme is established for the coupling hydraulic transient model of gas-liquid-solid three-phase mixed flow in the pipelines. The flux subvectors are then discretized by the Lax-Wendroff central difference scheme and the Warming-Beam upwind difference scheme with second-order precision in both time and space. Under the Rankine-Hugoniot conditions and the corresponding boundary conditions, an effective solution to those points located at the boundaries is developed, which can avoid the problem beyond the calculation region directly induced by the second-order discrete technique. Numerical and experimental verifications indicate that the proposed scheme has several desirable advantages including high calculation precision, excellent shock wave capture capability without false numerical oscillation, low sensitivity to the Courant number, and good stability.     

MULTI-OBJECTIVE OPTIMIZATION AND AERODYNAMIC PERFORMANCE ANALYSIS OF THE UPPER SURFACE FOR HYPERSONIC VEHICLES

为分析小攻角巡航条件下吸气式高超声速飞行器上壁面的变化对其气动性能和容积的影响, 以参数化后的飞行器上壁面对称面型线为设计变量, 在飞行马赫数6.5, 飞行高度27 km, 飞行攻角为4°的条件下, 采用计算流体力学为性能分析工具, Pareto多目标遗传算法为优化设计方法, 开展了二维条件下的升阻比/容积双目标优化设计. 在此基础上, 选择典型的二维优化结果, 重构生成对应的三维构型并进行数值分析, 获得了飞行器气动性能和容积间的相互关系. 结果表明在巡航条件下, 尽管二维/三维条件下飞行器的气动参数数值有较大差别, 但在这2种条件下, 飞行器的升阻比和容积间的关系均近似呈线性反比例关系. 同时, 对于三维构型而言, 在给定容积不变的条件下, 通过改变上壁面对称面型线的形状仅能使升阻比获得较小的增量(约0.36%). 相比之下, 当给定升阻比基本不变的条件下, 飞行器容积可调空间相对较大, 约为1.93%. 此外, 计算结果还表明, 在飞行器的容积基本不变情况下, 通过调节上壁面对称面型线, 可使飞行器的俯仰力矩获得5%左右的调节空间, 且其升阻比基本不变.     

Nonlinear water wave interaction with floating bodies in SPH

A weakly compressible SPH solver is presented for applications involving nonlinear interaction between water waves and floating bodies. A complete algorithm able to compute fully coupled viscous Fluid–Solid interactions is described. No slip boundary condition on the solid surface is enforced through a ghost–fluid technique and the global loads are evaluated through the momentum exchange between fluid and ghost particles. A dedicated algorithm is developed to manage the intersection between the free surface and the solid profile. An explicit synchronous algorithm is proposed for the full coupling between fluid and rigid bodies. Stability, convergence and conservation properties are tested on several freely floating test cases and a final validation of the full algorithm is performed for the interaction between a 2-D box and a wave packet.