二氧化碳在孔隙咸水层中多相流动数值模拟研究

There has been a growing concern about the use of fossil fuels and its adverse effects on the atmospheric greenhouse and ecological environment. A reduction in the release rate of CO2 into the atmosphere poses a major challenge to the land ecology of China. The most promising way of achieving CO2 reduction is to dispose of CO2 in deep saline aquifers. Deep aquifers have a large potential for CO2 sequestration in geological medium in terms of volume and duration. Through the numerical simulation of multiphase flow in a porous media, the transformation and motion of CO2 in saline aquifers has been implemented under various temperature and hydrostatic pressure conditions, which plays an important role to the assessment of the reliability and safety of CO2 geological storage. As expected, the calculated results can provide meaningful and scientific information for management purposes. The key problem to the numerical simulation of multiphase flow in a porous media is to accurately capture the mass interface and to deal with the geological heterogeneity. In this study, the updated CE/SE (Space and time conservation element and solution element) method has been proposed, and the Hybrid Particle Level Set method (HPLS) has extended for multiphase flows in porous medium, which can accurately trace the transformation of the mass interface. The benchmark problems have been applied to evaluate and validate the proposed method. In this study, the reliability of CO2 storage in saline aquifers in Daqingzi oil field in Sunlong basin has been discussed. The simulation code developed in this study takes into account the state for CO2 covering the triple point temperature and pressure to the supercritical region. The geological heterogeneity has been implemented, using the well known geostatistical model (GSLIB) on the base of the hard data. The 2D and 3D model have been set up to simulate the CO2 multiphase flow in the porous saline aquifer, applying the CE/SE method and the HPLS method .The main contents and results are summarized as followings. (1) The 2D CE/SE method with first and second –order accuracy has been extended to simulate the multiphase flow in porous medium, which takes into account the contribution of source and sink in the momentum equation. The 3D CE/SE method with the first accuracy has been deduced. The accuracy and efficiency of the proposed CE/SE method have been investigated, using the benchmark problems. (2) The hybrid particle level set method has been made appropriate and extended for capturing the mass interface of multiphase flows in porous media, and the numerical method for level set function calculated has been formulated. (3) The closed equations for multiphase flow in porous medium has been developed, adept to both the Darcy flow and non-Darcy flow, getting over the limitation of Reynolds number to the calculation. It is found that Darcy number has a decisive influence on pressure as well as velocity given the Darcy number. (4) The new Euler scheme for numerical simulations of multiphase flows in porous medium has been proposed, which is efficient and can accurately capture the mass interface. The artificial compressibility method has been used to couple the velocities and pressure. It is found that the Darcy number has determinant effects on the numerical convergence and stability. In terms of the different Darcy numbers, the coefficient of artificial compressibility and the time step have been obtained. (5) The time scale of the critical instability for critical CO2 in the saline aquifer has been found, which is comparable with that of completely CO2 dissolved saline aquifer. (6) The concept model for CO2 multiphase flows in the saline aquifer has been configured, based on the temperature, pressure, porosity as well as permeability of the field site .Numerical simulation of CO2 hydrodynamic trapping in saline aquifers has been performed, applying the proposed CE/SE method. The state for CO2 has been employed to take into account realistic reservoir conditions for CO2 geological sequestration. The geological heterogeneity has been sufficiently treated , using the geostatistical model. (7) It is found that the Rayleigh-Taylor instability phenomenon, which is associated with the penetration of saline fluid into CO2 fluid in the direction of gravity, has been observed in CO2 multiphase flows in the saline aquifer. Development of a mushroom-type spike is a strong indication of the formation of Kelvin-Helmholtz instability due to the developed short wavelength perturbations present along the interface and parallel to the bulk flow. Additional key findings: the geological heterogeneity can distort the flow convection. The ascending of CO2 can induce the persistent flow cycling effects. The results show that boundary conditions of the field site have determinant effects on the transformation and motion of CO2 in saline aquifers. It is confirmed that the proposed method and numerical model has the reliability to simulate the process of the hydrodynamic trapping, which is the controlling mechanism for the initial period of CO2 storage at time scale of 100 years.

两层流体热毛细对流空间实验研究=Space Experiments of Thermo Capillary Convection in Two-liquid Layers

1999年,在我国实践5号卫星上完成了两层流体空间实验,实验研究两层不相混合流体的纯Marangoni对流（温度梯度与界面垂直）与热毛细对流（温度梯度方向与流体界面平行）.前者存在发生Marangoni对流的最小临界温差值△Tc,低于该值流体系统处于静止状态;后者中只要存在沿界面的温度梯度便会产生热毛细对流.空间实验采用石蜡和氟化液两层流体新体系,实现了平整的液-液交界面,并从卫星上传回上万幅数字图像.通过多幅图像叠加处理得到了定量的流速场.数值模拟计算分析了相应工况时对流流动的速度场,两者的流场结构和速度大小基本一致,实验验证了理论模型.

Boundary knot method for 2D and 3D Helmholtz and convection-diffusion problems under complicated geometry

The boundary knot method (BKM) of very recent origin is an inherently meshless, integration-free, boundary-type, radial basis function collocation technique for the numerical discretization of general partial differential equation systems. Unlike the method of fundamental solutions, the use of non-singular general solution in the BKM avoids the unnecessary requirement of constructing a controversial artificial boundary outside the physical domain. The purpose of this paper is to extend the BKM to solve 2D Helmholtz and convection-diffusion problems under rather complicated irregular geometry. The method is also first applied to 3D problems. Numerical experiments validate that the BKM can produce highly accurate solutions using a relatively small number of knots. For inhomogeneous cases, some inner knots are found necessary to guarantee accuracy and stability. The stability and convergence of the BKM are numerically illustrated and the completeness issue is also discussed.

Characteristics of Argon Laminar DC Plasma Jet at Atmospheric Pressure

利用特殊设计的等离子体发生器,选择等离子体产生的工艺参数,实现工艺过程的精确控制,在大气压环境下获得了性能稳定的氖气直流层流等离子体射流。与湍流等离子体射流长度较短、径向尺寸较大、工作噪音高等特点相比,层流等离于体射流长度可达到550mm,而且沿整个射流长度其径向尺寸维持不变,工作噪音很小。当气流量为120cm~3/s、弧电流在70-200A的范围时,射流长度随弧电流的增加而增加,热效率起初略有降低然后维持平稳。随气流量的增加,层流等离子体射流的热效率也增加,在弧电流为200A时,可以达到40%。实验中测

Instability of Two-Layer Rayleigh-Benard Convection with Interfacial Thermocapillary Effect

The linear instability analysis of the Rayleigh-Allarangoni-Benard convection in a two-layer system of silicon oil 10cS and fluorinert FC70 liquids are performed in a larger range of two-layer depth ratios H, from 0.2 to 5.0 for different total depth H less than or equal to 12 mm. Our results are different from the previous study on the Rayleigh-Benard instability and show strong effects of thermocapillary force at the interface on the time-dependent oscillations arising from the onset of instability convection.

A Continuation Method of Parameter Inversion for Non-Equilibrium Convection-Dispersion Equation

Based on the homotopy mapping, a globally convergent method of parameter inversion for non-equilibrium convection-dispersion equations (CDEs) is developed. Moreover, in order to further improve the computational efficiency of the algorithm, a properly smooth function, which is derived from the sigmoid function, is employed to update the homotopy parameter during iteration. Numerical results show the feature of global convergence and high performance of this method. In addition, even the measurement quantities are heavily contaminated by noises, and a good solution can be found.

Experimental Study on Thermocapillary Convection

In the present paper, the experimental studies on thermocapillary convection are reviewed. The author's interest is mainly focused on the onset of oscillatory thermocapillary convection, the features of oscillatory flow pattern, and the critical Marangoni number related with temperature and free surface oscillation. The coordinated measurement in a microgravity environment of a drop shaft is also addressed.

Experimental investigation on the chaotic phenomena in the wake of a natural thermal convection flow

Chaotic phenomena in the wake of thermal convection flow fields above a heating flat plate were investigated experimentally. A newly developed electron beam fluorescence technique (EBF) was used to simultaneously measure density fluctuation at 7 points in a cross section above the plate. Correlation dimensions, intermittence coefficients, Fourier spectrum have been obtained for different Grashof numbers. Spatial distribution of correlation dimensions are presented. The experimental result shows that there is a certain relationship between the density fluctuation and the Gr number. And time-spacial characteristic of chaos evolution is also given.

Numerical simulation of magnetic field design for damping thermocapillary convection in a floating half-zone

The magnetic fields produced by electrical coils are designed for damping the the thermocapillary convection in a floating half-zone in microgravity. The fields are designed specially to reduce the flow near the free surface and then in the melt zone by adjusting the longitudinal coil positions close to the melt zone. The effects of the designed magnetic fields on reducing the flow velocity and temperature distribution non-uniformity in the melt zone are stronger than those of the case of an uniform longitudinal magnetic field obtained by numerical simulation, particularly at the melt-rod interface. It brings fundamental insights into the heat and mass transfer control at the solidification interface by the magnetic field design for crystal growth by the floating full-zone method.

Study on Flow Patterns of Oscillatory Thermocapillary Convection In Floating Half Zone by PIV

The velocity fields of oscillatory convection have been measured using the techniques of Particle Image Velocimetry (PIV) in a liquid bridge of half floating zone with small typical scales of a few millimeters for emphasizing the thermocapillary effect in comparison with the effect of buoyancy. The flow patterns of the oscillatory flow have been studied experimentally in a liquid bridge. The flow patterns in the liquid bridge are classified with mode numbers according to oscillatory flow characteristics. Results of the experiment show that the mode depends on the aspect ratio as well as the volume ratio of the liquid bridge. The experimental results are helpful for studying the structure of flow at the onset of oscillatory thermocapillary convection in a liquid bridge.

Effects of Natural Convection on the Characteristics of a Long Laminar Argon Plasma Jet Issuing Horizontally into Ambient Air

Three-dimensional modeling results show that the appearance of the long laminar plasma jet is less influenced by natural convection even as it is issuing into ambient air horizontally. However, plasma parameter distributions may deviate from axi-symmetry

Influence of the thickness ratio of air layer to liquid layer on the onset of Benard-Marangoni convection

An experimental investigation of the onset of Benard-Marangoni convection has been performed in a liquid layer of rectangular configuration. The critical temperature difference was measured via the detections of both temperature field pattern (IR-imaging) on the free surface and fluid convection (PIV) in the liquid layer. The critical temperature difference or the critical Marangoni number was given. The experiments were performed for a fixed depth of air layer and a changeable depth of the liquid layer, and then the influence of the thickness ratio of the air layer to liquid layer on the Marangoni instability was studied.

Characters of Surface Deformation and Surface Wave in Thermal Capillary Convection

In the field of fluid mechanics, free surface phenomena is one of the most important physical processes. In the present research work, the surface deformation and surface wave caused by temperature difference of sidewalls in a rectangular cavity have been investigated. The horizontal cross-section of the container is 52 mmx42 mm, and there is a silicon oil layer of height 3.5 mm in the experimental cavity. Temperature difference between the two side walls of the cavity is increased gradually, and the flow on the liquid layer will develop from stable convection to un-stable convection. An optical diagnostic system consisting of a modified Michelson interferometer and image processor has been developed for study of the surface deformation and surface wave of thermal capillary convection. The Fourier transformation method is used to interferometer fringe analysis. The quantitative results of surface deformation and surface wave have been calculated from a serial of the interference fringe patterns.The characters of surface deformation and surface wave have been obtained. They are related with temperature gradient and surface tension. Surface deformation is fluctuant with time, which shows the character of surface wave. The cycle period of the wave is 4.8 s, and the amplitudes are from 0 to 0.55 mu m. The phase of the wave near the cool side of the cavity is opposite and correlative to that near the hot side. The present experiment proves that the surface wave of thermal capillary convection exists on liquid free surface, and it is wrapped in surface deformation.

Experimental Study on the Thermal Argon Plasma Generation and Jet Length Change Characteristics at Atmospheric Pressure

The generation, jet length and flow-regime change characteristics of argon plasma issuing into ambient air have been experimentally examined. Different torch structures have been used in the tests. Laminar plasma jets can be generated within a rather wide range of working-gas flow rates, and an unsteady transitional flow state exists between the laminar and turbulent flow regimes. The high-temperature region length of the laminar plasma jet can be over an order longer than that of the turbulent plasma jet and increases with increasing argon flow rate or arc current, while the jet length of the turbulent plasma is less influenced by the generating parameters. The flow field of the plasma jet has very high radial gradients of plasma parameters, and a Reynolds number alone calculated in the ordinary manner may not adequately serve as a criterion for transition. The laminar plasma jet can have a higher velocity than that of an unsteady or turbulent jet. The long laminar plasma jet has good stiffness to withstand the impact of laterally injected cold gas and particulate matter. It could be used as a rather ideal object for fundamental studies and be applied to novel materials processing due to its attractive stable and adjustable properties.