高超声速前台阶绕流的计算分析

气动力和气动热的精确计算是高超声速流动分析的关键和难点，一直以来没有得到解决。本 文选择前台阶绕流这一经典问题，分析来流马赫数为10 的氩气流动以及台阶壁面的气动性能。通 过求解Navier-Stokes 方程和直接Monte-Carlo 粒子模拟，得到了高超声速前台阶绕流的精确流 动，并详细分析了计算网格对Navier-Stokes 方程求解结果的影响，同时解释了高超声速流动中热 流和摩阻计算困难的原因。通过壁面边界层的准一维假设，首次提出了壁面热流和剪切力的层流壁 面函数。研究发现，壁面函数很好地描述了高超声速层流边界层内的温度分布和速度分布，为高超 声速流动的气动力和气动热预估提供了有力的工具。

Simulation of gas flow over a micro cylinder

Gas flow over a micro cylinder is simulated using both a compressible Navier-Stokes solver and a hybrid continuum /particle approach. The micro cylinder flow has low Reynolds number because of the small length scale and the low speed, which also indicates that the rarefied gas effect exists in the flow. A cylinder having a diameter of 20 microns is simulated under several flow conditions where the Reynolds number ranges from 2 to 50 and the Mach number varies from 0.1 to 0.8. It is found that the low Reynolds number flow can be compressible even when the Mach number is less than 0.3, and the drag coefficient of the cylinder increases when the Reynolds number decreases. The compressible effect will increase the pressure drag coefficient although the friction coefficient remains nearly unchanged. The rarefied gas effect will reduce both the friction and pressure drag coefficients, and the vortex in the flow may be shrunk or even disappear.

Numerical simulation of condensation of bubbles under microgravity conditions by moving mesh method in the double-staggered grid

A numerical 2D method for simulation of two-phase flows including phase change under microgravity conditions is presented in this paper, with a level set method being coupled with the moving mesh method in the double-staggered grid systems. When the grid lines bend very much in a curvilinear grid, great errors may be generated by using the collocated grid or the staggered grid. So the double-staggered grid was adopted in this paper. The level set method is used to track the liquid-vapor interface. The numerical analysis is fulfilled by solving the Navier-Stokes equations using the SIMPLER method, and the surface tension force is modeled by a continuum surface force approximation. A comparison of the numerical results obtained with different numerical strategies shows that the double-staggered grid moving-mesh method presented in this paper is more accurate than that used previously in the collocated grid system. Based on the method presented in this paper, the condensation of a single bubble in the cold water under different level of gravity is simulated. The results show that the condensation process under the normal gravity condition is different from the condensation process under microgravity conditions. The whole condensation time is much longer under the normal gravity than under the microgravity conditions.

Smoothed Particle Hydrodynamics (SPH): an Overview and Recent Developments

Smoothed particle hydrodynamics (SPH) is a meshfree particle method based on Lagrangian formulation, and has been widely applied to different areas in engineering and science. This paper presents an overview on the SPH method and its recent developments, including (1) the need for meshfree particle methods, and advantages of SPH, (2) approximation schemes of the conventional SPH method and numerical techniques for deriving SPH formulations for partial differential equations such as the Navier-Stokes (N-S) equations, (3) the role of the smoothing kernel functions and a general approach to construct smoothing kernel functions, (4) kernel and particle consistency for the SPH method, and approaches for restoring particle consistency, (5) several important numerical aspects, and (6) some recent applications of SPH. The paper ends with some concluding remarks.

Numerical analysis of LEC growth of GaAs with an axial magnetic field

A set of numerical analyses for momentum and heat transfer For a 3 in. (0.075 m) diameter Liquid Encapsulant Czochralski (LEC) growth of single-crystal GaAs with or without all axial magnetic field was carried Out using the finite-element method. The analyses assume a pseudosteady axisymmetric state with laminar floats. Convective and conductive heat transfers. radiative heat transfer between diffuse surfaces and the Navier-Stokes equations for both melt and encapsulant and electric current stream function equations Cor melt and crystal Lire considered together and solved simultaneously. The effect of the thickness of encapsulant. the imposed magnetic field strength as well as the rotation rate of crystal and crucible on the flow and heat transfer were investigated. (C) 2002 Published by Elsevier Science Ltd.

Modelling of submarine landslides and generated water waves

The slide of unstable sedimentary bodies and their hydraulic effects are studied by numerical means. A two-dimensional fluid mechanics model based on Navier-Stokes equations has been developed considering the sediments and water as a mixture. Viscoplastic and diffusion laws for the sediments have been introduced into the model. The numerical model is validated with an analytical solution for a Bingham flow. Laboratory experiments consisting in the slide of gravel mass have been carried out. The results of these experiments have shown the importance of the sediment rheology and the diffusion. The model parameters are adjusted by trial and error to match the observed “sandflow”.

Two-dimensional numerical simulations of shock waves in micro convergent-divergent nozzles

The steady two-dimensional Navier-Stokes equations with the slip wall boundary conditions were used to simulate the supersonic flow in micro convergent-divergent nozzles. It is observed that shock waves can take place inside or outside of the micronozzles under the earth environment. For the over-expanded flows, there is a boundary layer separation point, downstream of which a wave interface separates the viscous boundary layer with back air flow and the inviscid core flow. The oblique shock wave is followed by the bow shock and shock diamond. The viscous boundary layer thickness relative to the whole nozzle width on the exit plane is increased but attains the maximum value around of 0.5 and oscillates against this value with the continuous increasing of the nozzle upstream pressures. The viscous effect either changes the normal shock waves outside of the nozzle for the inviscid flow to the oblique shock waves inside the nozzle, or transfers the expansion jet flow without shock waves for the inviscid flow to the oblique shock waves outside of the nozzle. 

多相聚合物体系相分离动力学的理论和模拟

本文提出了将自洽场理论与多相流格子Boltzmann方法相结合的模型，并从该模型出发推导出了描述动量守恒的Navier-Stokes方程和扩散方程，验证了模型理论上的正确性。应用此模型，对聚合物的相分离过程进行了模拟。 首先证实了本模型最终能够得到正确的热力学平衡结果。对于高分子共混物和嵌段共聚物相分离的动力学过程。在分相各个阶段，对高分子共混物和嵌段共聚物都进行了验证。 其次，应用所提出的格子Boltzmann模型，分别对二元聚合物共混物和二嵌段共聚物的相分离后期相区增长过程进行了研究。 最后，通过模型的进一步扩展，实现了对具有复杂结构的嵌段共聚物和复杂共混物的模拟。

单晶硅化学气相沉积反应器流场初探

对单晶硅化学气相沉积(CVD)反应器在沉积过程中的流场进行了初步分析。通过数值求解三维层流Navier-Stokes方程，研究了反应器内浮力效应所引起的流场对称性破坏。结果表明，由于存在浮力效应，轴对称几何体中也会发生非轴对称流场分布，从而影响单晶硅的均匀生长。

Statistical analysis of microscale gas flows and stock price changes

First, recent studies on the information preservation (IP) method, a particle approach for low-speed micro-scale gas flows, are reviewed. The IP method was validated for benchmark issues such as Couette, Poiseuille and Rayleigh flows, compared well with measured data for typical internal flows through micro-channels and external flows past micro flat plates, and combined with the Navier-Stokes equations to be a hybrid scheme for subsonic, rarefied gas flows. Second, the focus is moved to the microscopic characteristic of China stock market, particularly the price correlation between stock deals. A very interesting phenomenon was found that showed a reverse transition behaviour between two neighbouring price changes. This behaviour significantly differs from the transition rules for atomic and molecular energy levels, and it is very helpful to understand the essential difference between stock markets and nature.

零攻角NADTPS重构高超声速流场气动力/热研究

长时间、远距离高超声速飞行要求飞行器设计兼顾减阻和防热特点,论文提出无烧蚀自适应减阻防热新概念(NADTPS)并展开相关研究。NADTPS基于支杆-钝头体结构重构流场,利用侧向射流重点防护支杆头部。本文通过数值求解二维轴对称Navier-Stokes方程,对零攻角下NADTPS重构流场进行模拟,研究流场基本特征、气动力/热特性,分析新概念的减阻和防热机理。

超声速平板绕流的全流域阻力特性研究

超声速流动中摩擦阻力的准确预测一直是气动研究中的难点。本文以零攻角的平板绕流为对象，结合传统的CFD方法求解Navier-Stokes方程和直接模拟Monte-Carlo（DSMC）方法模拟粒子的微观运动，研究平板在超声速来流下的全流域阻力特性。研究发现，在超声速流中会产生激波-边界层相互干扰作用并随着来流马赫数的增加而加强，同时平板前缘存在非平衡效应显著的区域；努森数的增大使激波和边界层均增厚，粘性作用增强，非平衡区的范围增大，非平衡强度增强，而激波强度减弱；最后通过拟合连续流条件下的摩阻系数和修正自由分子流结果，结合桥函数，获得了全流域范围内均有效的平板阻力系数表达式。

支杆-钝头体带攻角飞行流场和“军刺”挡板作用机理研究

本文通过三维Navier-Stokes方程的数值求解，研究了支杆-钝头体结构在10°攻角M∞=6.0飞行条件下的流场结构和特点，指出其气动力特性恶化的原因，提出采用“军刺”挡板改善流场和气动力特性，并通过对比两种不同挡板作用下的流场和气动力特性变化分析其作用机理，发现“军刺”挡板结构分割流场抑制三维效应形成的周向流动，迎风面形成稳定的回流区和剪切层结构，将迎风面锥激波推离轴线，降低钝头体肩部流动结构相互作用强度，并在一定程度上缓解背风面流动干扰，明显改善支杆-钝头体带攻角飞行时的气动力特性。

Heat Transfer Problems Induced By Multi-shocks Interaction

Numerical simulations of the multi-shock interactions observable around hypersonic vehicles were carried out by solving Navier-Stokes equations with the AUSMPW scheme and the new type of the IV interaction created by two incident shock waves was investigated in detail. Numerical results show that the intersection point of the second incident shock with the bow shock plays important role on the flow pattern, peak pressures and heat fluxes. In the case of two incident shocks interacting with the bow shock at the same position, the much higher peak pressure and more severe heat transfer rate are induced than the classical IV interaction. The phenomenon is referred to as the multi-shock interaction and higher requirements will be imposed on thermal protection systems.