甲烷-煤尘爆炸波与障碍物相互作用的数值研究

为了探讨煤矿瓦斯和煤尘爆炸的物理机制 ,基于对甲烷和煤尘爆炸传播的理论分析 ,采用数值模拟方法研究障碍物对甲烷和煤尘爆炸传播的影响 ;为了更系统地考虑颗粒相各种输运特性、相间滑移和耦合 ,采用双流体模型建立了数学模型 ,该模型的出发点是把颗粒群和气体都作为连续介质 ,两者相互渗透充满整个空间形成没有间隙的“流体”———拟流体 ,在欧拉坐标系下考虑气 -固两相流动 ,气相和颗粒相的计算网格统一 ,易于处理。数值方法采用LU分解和迎风TVD格式分别处理对流项的隐式和显式部分 ,扩散项采用中心差分 ;同时研究了不同形状的障碍物对流场的影响 ,计算结果与实验吻合较好     

Thermal Effects on the Airflow in a Street Canyon – Nantes'99 Experimental Results and Model Simulations

The thermal effects on the airflow within a street canyon, which are produced by the variation of direct solar heating of the street sides and ground, are examined in this article. The investigation is based on the experimental results of the Nantes'99 campaign and numerical simulations performed with the Computational Fluid Dynamics (CFD) code CHENSI using the standard k- model. The Nantes'99 experimental campaign was performed in a North-to-South oriented central street canyon of Nantes, France. It was observed that a thin thermal layer develops locally within a few centimetres from the heated wall. It is anticipated that, the convective flow close to the windward wall, which was visualised during the experiment, carries air masses from the street level upwards, where normally cleaner air is transported. Consequently, thermal effects may be important for the air quality in the street.Based on the temperature and wind flow measurements, the flow and temperature fields were simulated first in two dimensions with the CFD code CHENSI. It was found that CHENSI overestimates the thermal effects on the canyon airflow showing the main re-circulation simulated in the isothermal case to change into two counter-rotating vortices after the inclusion of the heating of the windward wall. A reason for this overestimation is possibly the temperature wall function implemented for such thin thermal boundary layers in conjunction with the limitations in grid resolution.     

高压淹没非自由射流压力场的数值解法

对有限空间中的淹没非自由射流流压力场的数值解法进行了较为详细的研究和分析 ,并建立了相应的流场力学模型和压力解的第二类边界条件 ,给出了流场的矩形网格几何和边界标记。文中指出 ,精确的压力解只能从泊松形式的压力方程得到 ,且与边界相邻的点和内点应分别采用不同的方程进行求解     

Multiphase transfer processes in waste rock piles producing acid mine drainage: 2. Applications of numerical simulation

Acid mine drainage (AMD) results from the oxidation of sulfides, mainly pyrite, present in mine wastes, either mill tailings or waste rock. This is the second of two papers describing the coupled physical processes taking place in waste rock piles undergoing AMD production. Since the oxidation of pyrite involves the consumption of oxygen and the production of heat, the oxidation process initiates coupled processes of gas transfer by diffusion and convection as well as heat transfer. These processes influence the supply of oxygen that is required to sustain the oxidation process. This second paper describes a numerical simulator used to represent the interaction of these coupled transfer processes. Numerical simulations are applied to two large sites with extensive characterization programs and widely different properties and behavior that were described in the first paper. The South Dump of the Doyon mine in Canada is permeable and has a high pyrite oxidation rate, thus making temperature-driven air convection the main oxygen supply mechanism. The Nordhalde of the Ronnenberg mining district in Germany contains lower permeability material which is less reactive, thus leading to a more balanced contribution of gaseous diffusion and convection as oxygen supply mechanisms. Overall, simulations allow a coherent representation of the conditions monitored within the waste rock piles and the confirmation of their physical properties. Conceptual simulations are also carried out to illustrate the potential effect of border membranes and layered co-mingling as mitigation methods used to control AMD production in either active or future waste rock piles.     

One-dimensional simulation of solute transfer in saturated–unsaturated porous media using the discontinuous finite elements method

A one-dimensional transport model for simulating water flow and solute transport in homogeneous–heterogeneous, saturated–unsaturated porous media is presented. The model is composed of a combination of accurate numerical algorithms for solving the nonlinear Richard's and advection–dispersion equations (ADE). The mixed form of Richard's equation is solved using a standard finite element method (FEM) with primary variable switching. The transport equation is solved using operator splitting, with the discontinuous finite element method (DFE) for discretization of the advective term. A slope limiting procedure for DFE avoids numerical instabilities but creates very limited numerical dispersion for high Peclet numbers. An implicit finite differences scheme (FD) is used for the dispersive term.The unsaturated flow and transport model (Wamos-T) is applied to a variety of rigorous problems including transient flow, heterogeneous medium and abrupt variations of velocity in magnitude and direction due to time-varying boundary conditions. It produces accurate and mass-conservative solutions for a very large range of grid Peclet numbers. The Wamos-T model is a good and robust alternative for the simulation of mass transport in unsaturated domain.     

Influence of non-reactive particle cloud on heterogeneous detonation propagation

An interaction of a detonation wave propagating in the cellular detonation mode with a cloud of inert particles is investigated numerically. The analysis of results allows the regimes of propagation of the heterogeneous plane Chapman–Jouguet and cellular detonations and their suppression to be identified. The influence of various parameters of the inert cloud is demonstrated. The critical length of the cloud sufficient for detonation suppression is determined. It is shown that the disperse composition and the nonuniform distribution of particles in the cloud are important parameters affecting the detonation propagation mode.     

Simulating on rollover phenomenon in LNG storage tanks and determination of the rollover threshold

Rollover is a potential risk to the safety of LNG storage tanks during the LNG storage process, so study of its prevention method is very important. In this paper, rollover phenomenon in a liquefied natural gas (LNG) storage tank is modeled physically and mathematically. Its evolution is simulated using FLUENT™ software from the breakdown of stratification to the occurrence of rollover. Results show that the evolution consists of three phases: the initial phase where rollover occurs near the side wall of the storage tank; the turbulent phase where rollover transfers to the center of the tank; and the final phase where new layers evolve. Based on these phases, rollovers in 160,000, 30,000, and 5000 m³ LNG storage tanks are simulated at varying initial density differences, and a rollover coefficient is defined to describe rollover intensity. The simulations show that the rollover coefficient initially increases within a small scope and then increases rapidly with the increment of initial density difference. This turning point is chosen to be the rollover threshold, which is regarded as the critical density difference in this study. The critical density differences obtained from the simulation results of the 160,000, 30,000, and 5000 m³ LNG storage tanks are 3, 5, and 7 kg/m³, respectively, which can be used as their rollover criteria to ensure the safety of LNG storage tanks.     

Numerical analysis of keyhole geometry and temperature profiles in plasma arc welding

In keyhole plasma arc welding (PAW), the formation of the keyhole involves complicated thermo-physical mechanism. The shape and dimensions of the keyhole directly determine the heat deposition along the thickness direction of the test plates, the penetration depth and the weld bead quality. Based on analysis of the forces acting on the weld pool, a keyhole model is established according to the force-balance condition on the keyhole wall. The shape and dimensions of quasi-steady state keyhole are numerically predicted. The keyhole shape and size are used to calibrate some distribution parameters of the combined volumetric heat source which takes into consideration the geometrical feature of welds with larger ratio of penetration depth to width and volumetric distribution characteristic of the plasma heat intensity along the plate thickness direction. The temperature fields on stainless steel plates are numerically simulated. The calculated fusion line at the transverse cross-section of PAW welds are compared with the experimental measurements. It is found that the predicted results are in agreement with the experimental data.     

Scenario simulations of CO2 injection feasibility, plume migration and storage in a saline aquifer, Scania, Sweden

Deep saline aquifers have large capacity for geological CO₂ storage, but are generally not as well characterized as petroleum reservoirs. We here aim at quantifying effects of uncertain hydraulic parameters and uncertain stratigraphy on CO₂ injectivity and migration, and provide a first feasibility study of pilot-scale CO₂ injection into a multilayered saline aquifer system in southwest Scania, Sweden. Four main scenarios are developed, corresponding to different possible interpretations of available site data. Simulation results show that, on the one hand, stratigraphic uncertainty (presence/absence of a thin mudstone/claystone layer above the target storage formation) leads to large differences in predicted CO₂ storage in the target formation at the end of the test (ranging between 11% and 98% of injected CO₂ remaining), whereas other parameter uncertainty (in formation and cap rock permeabilities) has small impact. On the other hand, the latter has large impact on predicted injectivity, on which stratigraphic uncertainty has small impact. Salt precipitation at the border of the target storage formation affects CO₂ injectivity for all considered scenarios and injection rates. At low injection rates, salt is deposited also within the formation, considerably reducing its availability for CO₂ storage.