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

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

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.     

Increasing accuracy of dispersal kernels in grid-based population models

Dispersal kernels in grid-based population models specify the proportion, distance and direction of movements within the model landscape. Spatial errors in dispersal kernels can have large compounding effects on model accuracy. Circular Gaussian and Laplacian dispersal kernels at a range of spatial resolutions were investigated, and methods for minimizing errors caused by the discretizing process were explored. Kernels of progressively smaller sizes relative to the landscape grid size were calculated using cell-integration and cell-center methods. These kernels were convolved repeatedly, and the final distribution was compared with a reference analytical solution. For large Gaussian kernels (σ > 10 cells), the total kernel error was <10⁻¹¹ compared to analytical results. Using an invasion model that tracked the time a population took to reach a defined goal, the discrete model results were comparable to the analytical reference. With Gaussian kernels that had σ ≤ 0.12 using the cell integration method, or σ ≤ 0.22 using the cell center method, the kernel error was greater than 10%, which resulted in invasion times that were orders of magnitude different than theoretical results. A goal-seeking routine was developed to adjust the kernels to minimize overall error. With this, corrections for small kernels were found that decreased overall kernel error to <10⁻¹¹ and invasion time error to <5%.     

Quantifying the effect of strong ignition sources on particle preconditioning and distribution in the 20-L chamber

Computational fluid dynamics is used to investigate the preconditioning aspect of overdriving in dust explosion testing. The results show that preconditioning alters both the particle temperature and distribution prior to flame propagation in the 20-L chamber. A parametric study gives the fluid pressure and temperature, and particle temperature and concentration at an assumed flame kernel development time (10 ms) for varying ignitor size and particle diameter. For the 10 kJ ignitor with 50% efficiency, polyethylene particles under 50 μm reach 400 K and may melt prior to flame propagation. Gases from the ignitor detonation displace the dust from the center of the chamber and may increase local particle concentration up to two times the nominal value being tested. These effects have important implications for explosive testing of dusts in the 20-L chamber and comparing to larger 1-m³ testing, where these effects may be negligible.     

底层设有商业服务网点的高层居民住宅楼 消防安全问题研究

高层建筑已经成为城市建筑中最主要的建筑形式,底层设有商业服务网点的高层居民住宅楼由于兼具商业和居住的功能而受到开发商和用户的青睐。由于对底层商业服务网点的使用需求不同进行的改建、扩建带来的较为突出的消防安全问题,给城市的消防安全管理增加了不确定因素。采用数值模拟的方法构建了底层商业服务网点改扩建后的火灾场景,研究了此类场所改扩建后的火灾危险性,对增设自动喷水灭火系统的消防设计方法可行性进行了论证。模拟结果显示,仅在改扩建网点增设自动喷水灭火系统可以达到阻止火灾蔓延并扑灭火灾的安全要求。