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

借助于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.

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

本文对二相流中的颗粒运动和碰撞情况进行了研究。颗粒运动和碰撞是二相流中的现象，长期以来，人们把各种情况下的颗粒无规则运动，一概称为湍流脉动，事实上，除湍流脉动外，颗粒还有因颗粒－颗粒碰撞和颗粒－壁面碰撞引起的、二相流所特有的（在单相流中不存在的）无规则运动。研究颗粒的无规则运动有助于研究二相流方程的本构关系，为更深一步的二相流研究打下一定的基础。在数值解法上，采用直接模拟蒙特卡洛方法（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新算法具有节点少、精度高,效率高,编程方便等优点,能成功模拟液液两相流动,为两相流动数值模拟提供了一种新的途径.

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.

A Simple Model For Predicting The Void Fraction Of Gas/Non-Newtonian Fluid Intermittent Flows In Upward Inclined Pipes

In this work, a simple correlation, which incorporates the mixture velocity, drift velocity, and the correction factor of Farooqi and Richardson, was proposed to predict the void fraction of gas/non-Newtonian intermittent flow in upward inclined pipes. The correlation was based on 352 data points covering a wide range of flow rates for different CMC solutions at diverse angles. A good agreement was obtained between the predicted and experimental results. These results substantiated the general validity of the model presented for gas/non-Newtonian two-phase intermittent flows.

节流器内液-固两相流固体颗粒冲蚀数值模拟

建立了考虑颗粒碰撞的颗粒冲蚀计算模型,该数学模型包括:在Eulerian坐标系下求解连续相流场;在Lagrangian坐标系下运用离散颗粒硬球模型求解颗粒碰撞;应用半实验关联式求解颗粒冲蚀速率。对水力加砂压裂施工中节流器内液-固两相流的固体颗粒运动和冲蚀特性进行了数值模拟。计算结果表明,固体颗粒密集于节流器入口到出口的一段狭长区域内,冲蚀速率随流体速度呈指数性变化。颗粒直径越大,冲蚀速率也越大。节流器内冲蚀最严重的位置发生在距离节流器出口上边缘10mm以内的局部区域。

Nonlinear optics and wavelength translation via cavity-optomechanics

The field of cavity-optomechanics explores the interaction of light with sound in an ever increasing array of devices. This interaction allows the mechanical system to be both sensed and controlled by the optical system, opening up a wide variety of experiments including the cooling of the mechanical resonator to its quantum mechanical ground state and the squeezing of the optical field upon interaction with the mechanical resonator, to name two.

In this work we explore two very different systems with different types of optomechanical coupling. The first system consists of two microdisk optical resonators stacked on top of each other and separated by a very small slot. The interaction of the disks causes their optical resonance frequencies to be extremely sensitive to the gap between the disks. By careful control of the gap between the disks, the optomechanical coupling can be made to be quadratic to first order which is uncommon in optomechanical systems. With this quadratic coupling the light field is now sensitive to the energy of the mechanical resonator and can directly control the potential energy trapping the mechanical motion. This ability to directly control the spring constant without modifying the energy of the mechanical system, unlike in linear optomechanical coupling, is explored.

Next, the bulk of this thesis deals with a high mechanical frequency optomechanical crystal which is used to coherently convert photons between different frequencies. This is accomplished via the engineered linear optomechanical coupling in these devices. Both classical and quantum systems utilize the interaction of light and matter across a wide range of energies. These systems are often not naturally compatible with one another and require a means of converting photons of dissimilar wavelengths to combine and exploit their different strengths. Here we theoretically propose and experimentally demonstrate coherent wavelength conversion of optical photons using photon-phonon translation in a cavity-optomechanical system. For an engineered silicon optomechanical crystal nanocavity supporting a 4 GHz localized phonon mode, optical signals in a 1.5 MHz bandwidth are coherently converted over a 11.2 THz frequency span between one cavity mode at wavelength 1460 nm and a second cavity mode at 1545 nm with a 93% internal (2% external) peak efficiency. The thermal and quantum limiting noise involved in the conversion process is also analyzed and, in terms of an equivalent photon number signal level, are found to correspond to an internal noise level of only 6 and 4 times 10x^-3 quanta, respectively.

We begin by developing the requisite theoretical background to describe the system. A significant amount of time is then spent describing the fabrication of these silicon nanobeams, with an emphasis on understanding the specifics and motivation. The experimental demonstration of wavelength conversion is then described and analyzed. It is determined that the method of getting photons into the cavity and collected from the cavity is a fundamental limiting factor in the overall efficiency. Finally, a new coupling scheme is designed, fabricated, and tested that provides a means of coupling greater than 90% of photons into and out of the cavity, addressing one of the largest obstacles with the initial wavelength conversion experiment.

Ultrafast dynamics in ZnO thin films irradiated by femtosecond lasers

The damage morphologies, threshold fluences in ZnO films were studied with femtosecond laser pulses. Time-resolved reflectivity and transmissivity have been measured by the pump-probe technique at different pump fluences and wavelengths. The results indicate that two-phase transition is the dominant damage mechanism, which is similar to that in narrow band gap semiconductors. The estimated energy loss rate of conduction electrons is 1.5 eV/ps. (c) 2005 Elsevier Ltd. All rights reserved.