用激波管研究超音速气固两相流

在激波管里进行了可压缩性气固两相流的实验研究。测量了激波通过颗粒群时的压力的衰减过程。用纹影仪拍摄了激波与颗粒群相互干涉的照片。试验了颗粒群的不同构造对压力衰减的影响。指出了激波反射、聚焦等非线性气动因素是可压缩性气固两相流的关键问题。

具有蒸发液滴两相介质近壁流动的数学分析

在双连续介质理论框架下，采用匹配渐进展开方法导出并求解了具有蒸发液滴的汽雾流中层流边界层方程，给出了控制汽雾流的相似判据。对于沿曲面的流动，边界层方程的形式取决于是否存在液滴的惯性沉积。给出了热钝体驻点附近蒸汽-液滴边界层的数值计算结果。它们表明：由于蒸发，在边界层内近壁处形成了一个无液滴区域；在该区上边界处，液滴半径趋于零而液滴数密度急剧增高。液滴蒸发及聚集的联合效应造成了表面热流的显著增加，甚至在自由来流中液滴质量浓度很低时此效应依然存在。

气-液两相流界面迁移现象的数值模拟研究

借助于Level Set函数,建立了气-液两相的统一控制方程组,并在交错网格中进行离散.用两种格式,即Superbee-TVD格式和5阶WENO格式求解Level Set函数的输运方程,用SIMPLER算法的思想对主流场控制方程的求解方法进行改进.数值实验结果表明,在求解Level Set的控制方程时,5阶WENO方法比Superbee-TVD格式的结果更准确;用改进的数值算法可成功实现对密度比大于1 000/1的气-液两相流界面迁移问题的数值模拟.对几种典型大密度比气-液两相流问题的计算结果与实际问题的物理规律完全一致,验证了该方法的有效性和可靠性.

Design of a Self-Adaptive Down-Hole Layered Steam Injector and Its Application

Layered steam injection, widely used in Liaohe Oilfield at Present, is an effective recovery technique to heavy oil reserves. Which makes the steam front-peak push forward uniformly, the amount of steam injection be assigned rationally, and the effect of injection steam be obtained as expected. To maintain a fixed ratio of layered steam injection and solve the problem of nonadjustable hole diameter with the change of layer pressure in the existing injectors, a new method is proposed in this paper to design layered steam injectors based on the dynamic balance theory. According to gas-liquid two-phase flow theory and beat transfer theory, the energy equation and the heat conduction equation in boreholes are developed. By analyzing the energy equilibrium of water-steam passing through the injector hole, we find an expression to describe the relation between the cross-sectional area of injector hole and the layer pressure. With this expression, we provide a new set of calculation methods and write the corresponding computer program to design and calculate the main parameters of a steam injector. The actual measurement data show that the theoretically calculated results are accurate, the software runs reliably, and they provide the design of self-adjustable layered steam injectors with the theoretical foundation.

Undercooling and Rapid Solidification of Nb-Si Eutectic Alloys Studied by Long Drop Tube

Niobium-silicide alloys have great potential for high temperature turbine applications. The two-phase Nb/Nb5Si3 in situ composites exhibit a good balance in mechanical properties. Using the 52 in drop tube, the effect of undercooling and rapid solidification on the solidification process and micro-structural characterization of Nb-Si eutectic alloy was studied. The microstructures of the Nb-Si composites were investigated by optics microscope (OM), X-ray diffraction (XRD) and scanning electron microscope (SEM) equipped with X-ray energy dispersive spectrometry (EDS). Up to 480 K, deep undercooling of the Nb-Si eutectic samples was successfully obtained, which corresponds to 25% of the liquidus temperature. Contrasting to the conventional microstructure usually found in the Nb-Si eutectic alloy, the microstructure of the undercooled sample is divided into the fine and coarse regions. The most commonly observed microstructure is Nb+Nb5Si3, and the Nb3Si phase is not be found. The change of coarseness of microstructure is due to different cooling rates during and after recalescence. The large undercooling is sufficient to completely bypass the high temperature phase field.

Design of a Self-Adaptive Down-Hole Layered Steam Injector and Its Application

Layered steam injection, widely used in Liaohe Oilfield at present, is an effective recovery technique to heavy oil reserves. Which makes the steam front-peak push forward uniformly, the amount of steam injection be assigned rationally, and the effect of injection steam be obtained as expected. To maintain a fixed ratio of layered steam injection and solve the problem of nonadjustable hole diameter with the change of layer pressure in the existing injectors, a new method is proposed in this paper to design layered steam injectors based on the dynamic balance theory According to gas-liquid two-phase flow theory and heat transfer theory, the energy equation and the heat conduction equation in boreholes are developed. By analyzing the energy equilibrium of water-steam passing through the injector hole, we find an expression to describe the relation between the cross-sectional area of injector hole and the layer pressure. With this expression, we provide a new set of calculation methods and write the corresponding computer program to design and calculate the main parameters of a steam injector. The actual measurement data show that the theoretically calculated results are accurate, the software runs reliably, and they provide the design of self-adjustable layered steam injectors with the theoretical foundation.

Chemical equilibrium analysis and CFD simulation of coal combustion and NOx formation

A full two-fluid model of reacting gas-particle flows with an algebraic unified second-order moment (AUSM) turbulence-chemistry model is used to simulate Beijing coal combustion and NOx formation. The sub-models are the k-epsilon-kp two-phase turbulence model, the EBU-Arrhenius volatile and CO combustion model, the six-flux radiation model, coal devolatilization model and char combustion model. The blocking effect on NOx formation is discussed. In addition, the chemical equilibrium analysis is used to predict NOx concentration at different temperature. Results of CID simulation and chemical equilibrium analysis show that, optimizing air dynamic parameters can delay the NOx formation and decrease NOx emission, but it is effective only in a restricted range. In order to decrease NOx emission near to zero, the re-burning or other chemical methods must be used.

Numerical Sensitivity Analysis and Partial Similarity of Porous Media Flows

A numerical optimisation approach to identify dominant dimensionless variables in porous media flows by sensitivity analysis is proposed. We have validated the approach at first by examining a simple oil reservoir theoretically and numerically as well. A more complex water-flooding reservoir is examined based on sensitivity analysis of oil recovery to the similarity parameters, thus demonstrating the feasibility of the proposed approach to identify dominant similarity parameters for water-oil two-phase flows.

Pool Boiling in Microgravity: Recent Results and Perspectives for the Project DEPA-SJ10

Two research projects on pool boiling in microgravity have been conducted aboard the Chinese recoverable satellites. Ground-based experiments have also been performed both in normal gravity and in short-term microgravity in the Drop Tower Beijing. Steady boiling of R113 on thin platinum wires was studied with a temperature-controlled heating method, while quasi-steady boiling of FC-72 on a plane plate was investigated with an exponentially increasing heating voltage. In the first case, slight enhancement of heat transfer is observed in microgravity, while diminution is evident for high heat flux in the second one. Lateral motions of bubbles on the heaters are observed before their departure in microgravity. The surface oscillation of the merged bubbles due to lateral coalescence between adjacent bubbles drives it to detach from the heaters. The Marangoni effect on the bubble behavior is also discussed. The perspectives for a new project DEPA-SJ10, which has been planned to be flown aboard the Chinese recoverable satellite SJ-10 in the future, are also presented.

二相流中颗粒脉动和阵性泥石流立面二维运动特性研究

本文对二相流中的颗粒运动和碰撞情况进行了研究。颗粒运动和碰撞是二相流中的现象，长期以来，人们把各种情况下的颗粒无规则运动，一概称为湍流脉动，事实上，除湍流脉动外，颗粒还有因颗粒－颗粒碰撞和颗粒－壁面碰撞引起的、二相流所特有的（在单相流中不存在的）无规则运动。研究颗粒的无规则运动有助于研究二相流方程的本构关系，为更深一步的二相流研究打下一定的基础。在数值解法上，采用直接模拟蒙特卡洛方法（DSMC），用一定数目的模拟颗粒代表流场中的大量实际颗粒，每个模拟颗粒具有一确定的速度。颗粒的运动和碰撞解耦，先计算出颗粒一个时间步长内运动的轨迹和最终速度，下一步就是判断在这段进间内那些颗粒发生了碰撞，然后用碰撞理论同碰撞后的速度，而碰撞过程中的颗粒的坐标不变。我们计算采用的颗粒模型是光滑的硬球模型，不考虑旋转效应，只计二体碰撞，模拟颗粒的直径一定。判断两颗粒是否发生碰撞主要取决于碰撞概率，而不是运动轨迹。计算结果表明随着颗粒的驰豫时间的增加，颗粒的聚团现象变得剧烈。颗粒的脉动强度跟颗粒的体积浓度成一反比例函数关系，颗粒的各向脉动强度跟总脉动强度成正比关系。随着驰豫时间的增加，颗粒的总脉动强度增加，而主流方向的脉动强度在总脉动强度中占的比例减少。本文的研究表明，DSMC方法模拟颗粒运动和碰撞是可行的，得出的规律是可信的。

Large-eddy simulation of turbulent collision of heavy particles in isotropic turbulence

The small-scale motions relevant to the collision of heavy particles represent a general challenge to the conventional large-eddy simulation (LES) of turbulent particle-laden flows. As a first step toward addressing this challenge, we examine the capability of the LES method with an eddy viscosity subgrid scale (SGS) model to predict the collision-related statistics such as the particle radial distribution function at contact, the radial relative velocity at contact, and the collision rate for a wide range of particle Stokes numbers. Data from direct numerical simulation (DNS) are used as a benchmark to evaluate the LES using both a priori and a posteriori tests. It is shown that, without the SGS motions, LES cannot accurately predict the particle-pair statistics for heavy particles with small and intermediate Stokes numbers, and a large relative error in collision rate up to 60% may arise when the particle Stokes number is near St_K=0.5. The errors from the filtering operation and the SGS model are evaluated separately using the filtered-DNS (FDNS) and LES flow fields. The errors increase with the filter width and have nonmonotonic variations with the particle Stokes numbers. It is concluded that the error due to filtering dominates the overall error in LES for most particle Stokes numbers. It is found that the overall collision rate can be reasonably predicted by both FDNS and LES for St_K>3. Our analysis suggests that, for St_K<3, a particle SGS model must include the effects of SGS motions on the turbulent collision of heavy particles. The spectral analysis of the concentration fields of the particles with different Stokes numbers further demonstrates the important effects of the small-scale motions on the preferential concentration of the particles with small Stokes numbers.

New Numerical Algorithms in SUPER CE/SE and Their Applications in Explosion Mechanics

The new numerical algorithms in SUPER/CESE and their applications in explosion mechanics are studied. The researched algorithms and models include an improved CE/SE (space-time Conservation Element and Solution Element) method, a local hybrid particle level set method, three chemical reaction models and a two-fluid model. Problems of shock wave reflection over wedges, explosive welding, cellular structure of gaseous detonations and two-phase detonations in the gas-droplet system are simulated by using the above-mentioned algorithms and models. The numerical results reveal that the adopted algorithms have many advantages such as high numerical accuracy, wide application field and good compatibility. The numerical algorithms presented in this paper may be applied to the numerical research of explosion mechanics.

摄动有限体积算法及其在两相流中应用

采用一种全新的摄动有限体积(PFV)算法和水平集(Level Set)技术对液液两相系统中液滴坠落进行数值模拟,数值结果表明,PFV新算法具有节点少、精度高,效率高,编程方便等优点,能成功模拟液液两相流动,为两相流动数值模拟提供了一种新的途径.

Space experimental studies of microgravity fluid science in China

Microgravity fluid physics is an important part of microgravity sciences, which consists of simple fluids of many new systems, gas-liquid two-phase flow and heat transfer, and complex fluid mechanics. In addition to the importance of itself in sciences and applications, microgravity fluid physics closely relates to microgravity combustion, space biotechnology and space materials science, and promotes the developments of interdisciplinary fields. Many space microgravity experiments have been per- formed on board the recoverable satellites and space ships of China and pushed the rapid development of microgravity sciences in China. In the present paper, space experimental studies and the main re- sults of the microgravity fluid science in China in the last 10 years or so are introduced briefly.

Investigation on mechanism of critical cavitating flow in liquid jet pumps under operating limits

The critical cavitating flow in liquid jet pumps under operating limits is investigated in this paper. Measurements on the axial pressure distribution along the wall of jet pumps indicate that two-phase critical flow occurs in the throat pipe under operating limits. The entrained flow rate and the distribution of the wall pressure upstream lowest pressure section does not change when the outlet pressure is lower than a critical value. A liquid-vapor mixing shockwave is also observed under operating limits. The wave front moves back and forth in low frequency around the position of the lowest pressure. With the measured axial wall pressures, the Mach number of the two-phase cavitating flow is calculated. It's found that the maximum Mach number is very close to I under operating limits. Further analysis infers a cross-section where Mach number approaches to I near the wave front. Thus, the liquid-vapor mixture velocity should reach the local sound velocity and resulting in the occurrence of operating limits.