Contributions of 3D Printed Fracture Networks to Development of Flow and Transport Models

Transport in Porous Media - Conventional experiments using natural rock samples have trouble in observing rock structures and controlling fracture properties. Taking advantage of 3D printing...     

Large-Eddy and Direct Numerical Simulations of the Bachalo-Johnson Flow with Shock-Induced Separation

The Bachalo-Johnson experiment on an axisymmetric bump has been a primary validation case for turbulence models in shock-boundary-layer interactions since the 1980’s. In the present work, Wall-Modelled Large-Eddy Simulations (WMLES) of this flow were conducted using Improved Delayed Detached-Eddy Simulation (IDDES) as the sub-grid-scale (SGS) and wall model, with a synthetic turbulence generator, expecting close enough agreement with experiment. However, the WMLES results are disappointing, even in terms of the shock position, even though the results from two grids with 4.7 × 10⁸ and 1.6 × 10⁹ cells respectively agree well with each other. This failure of grid refinement to warn of an inaccurate simulation is of great concern, and the reasons for it are explored. We then conducted a Direct Numerical Simulation (DNS) embedded in the LES over a reduced domain, with 8 × 10⁹ grid cells. The DNS has a far more accurate shock position and overall pressure distribution. The skin friction in the favourable pressure gradient is also much higher than in the LES; thus, wide differences appear upstream of the shock wave, most probably caused by the rapid acceleration which leads to atypical shear-stress profiles. Other SGS models were tried, and performed worse than IDDES. The DNS essentially fulfils the initial expectations although in a reduced domain and provides data for turbulence-modelling studies, for instance by extracting an effective eddy viscosity from it. The most noticeable remaining disagreement with experiment is over the Reynolds shear stress.     

OH Concentration Measurements by Resonant Holographic Interferometry and Comparison with Direct Numerical Simulations

We apply a novel laser diagnostic technique — Resonant Holographic Interferometry (RHI) to measure the concentration of hydroxyl radical (∼2000 ppm) in a co-flow diffusion flame of diluted hydrogen and air stabilized on a Wolfhard-Parkerburner. This methodology is based upon the dispersion of light of frequency close to an electronic transition of a target molecule. The two-color setup utilized in RHI provides a two-dimensional distribution of the target species concentration and quantitative information can be obtained from the interferogram without requiring any calibration. To provide independent flame data for comparison, a two-dimensional numerical simulation was performed taking into account the effects of detailed chemical kinetics and transport phenomena. In spite of a number of simplifying assumptions made in the simulation, computational and experimental results are in good agreement with respect to the magnitude and width of the region where OH is found. We do observe a difference of approximately 1 mm in the flame position due to the simplifying assumptions made in the simulation. The comparison between the experimental and numerical results clearly demonstrated the potential of RHI in flame diagnostics. This revised version was published online in July 2006 with corrections to the Cover Date.     

Bulk Flow Regimes and Fractional Flow in 2D Porous Media by Numerical Simulations

We investigate a two-dimensional network simulator that models the dynamics of two-phase immiscible bulk flow where film flow can be neglected. We present a method for simulating the detailed dynamical process where the two phases are allowed to break up into bubbles, and bubbles are allowed to merge together. The notions of drainage and imbibition are not adequate to describe this process since there is no clear front between the fluids. In fact, the simulator is constructed so that one can study the behaviour of the system far from inlets and outlets, where the two fluids have been mixed together so much that all initial fronts have broken up. The simulator gives the fractional flow as a function of the saturation of each of the fluids. For the case of two fluids with equal viscosity, we classify flow regimes that are parametrized by the capillary number.     

Direct Numerical Simulations of Turbulent Flow Over Various Riblet Shapes in Minimal-Span Channels

Flow, Turbulence and Combustion - Riblets reduce skin-friction drag until their viscous-scaled size becomes large enough for turbulence to approach the wall, leading to the breakdown of...     

Fluid Flow Across Capillary Heterogeneities,Experiments and Simulations

A comparative study of numerical modelling and laboratory experiments of two-phase immiscible displacements in a 33 cm × 10 × 3 cm thick cross-bedded reservoir model is reported. Dynamic two-dimensional fluid saturation development was obtained from experiments by use of a nuclear tracer imaging technique and compared to numerical predictions using a full-field black oil simulator.The laboratory cross-bedded reservoir model was a sandpack consisting of two strongly waterwet sands of different grain sizes, packed in sequential layers. The inlet and outlet sand consisted of low permeable, high capillary, sand while the central crosslayer with a dip angle of 30° was a high permeable, low capillary, sand. Results on moderate contrasts in permeability and capillary heterogeneities in the cross-bedded reservoir model at different mobility ratios and capillary number floods temporarily showed a bypass of oil, resulting in a prolonged two-phase production. The final remaining oil saturations, however, were as for isolated samples. Hence, permanently trapped oil was not observed.Simulations of waterfloods, using a commercial software package, displayed correct water breakthrough at low flow rate and unity viscosity ratio, but failed in predicting local saturation development in detail, probably due to numerical diffusion.The simulator was used to test several cases of heterogeneity contrasts, and influence from different relative permeability curves. Further, by altering the capillary pressure at the outlet, the end effects were proven important.     

Large Eddy Simulations and Experiments of Flow Around Finite-Height Cylinders

The paper presents experimental and numerical results for the flow around a surface-mounted circular cylinder at the two height-to-diameter ratios of 2.5 and 5. The Reynolds number based on approach flow velocity and cylinder diameter is 43,000 and 22,000 for these two cases and the boundary layer of the approach flow has a thickness of about 10% of the cylinder height. The experiments comprise both flow visualizations with dye and laser Doppler velocimeter measurements of all mean velocity and fluctuation components. The numerical study is performed by an elaborate large eddy simulation on a staggered Cartesian grid using the immersed boundary method. The instantaneous flow behaviour including the shedding is analysed with information based on animations. For the long cylinder alternating shedding is found to occur over most of the height while for the shorter cylinder the shedding is observed mainly near the ground where it is also mostly alternating but intermittently also symmetrical. The mean-flow behaviour is analysed with the aid of streamlines and contour plots of mean-velocity and fluctuation components in various planes and a detailed comparison of LES and LDV results is provided, showing generally good agreement. The LES with very fine resolution near the free end allow a detailed study of the complex flow in this region with owl-face topology on the end wall previously observed in experiments. Behind the cylinder, the longitudinal recirculation region, the downstream development of tip vortices and the emergence of trailing vortices further downstream are analysed. The sum of the results, together with those from previous studies that were reviewed extensively, provides a comprehensive picture of the very complex flow behaviour.     

Validation of the Slow Compressional Wave in Porous Media: Comparison of Experiments and Numerical Simulations

We perform numerical simulation of ultrasonic experiments on poroelastic samples, in which Biot's slow compressional wave had been observed. The simulation is performed using OASES modeling code, which allows to compute elastic wave fields in layered poroelastic media. Modeled were the experiments of Plona (1980), Rasolofosaon (1988), and our own measurements. In all the three situations, a good agreement between experiment and simulations has been observed. This further confirms the fact that Biot's theory of poroelasticity, on which the simulations were based, adequately describes the behavior of the porous materials under investigations at ultrasonic frequencies.     

Simulation of Fluid Flow on Fractures and Implications for Reactive Transport Simulations

The diverse numerical simulation techniques employed to predict fluid flow properties of fractures yield differing results which limits their applicability for reactive transport simulations. Basically the fluid flow simulation techniques can be divided in two groups: (i) techniques that yield average fluid flow characteristics and (ii) techniques that produce space-resolved properties. These differences may have substantial impact on the reactive transport simulations but may also depend on the fracture characteristics. For this reason, a sensitivity analysis of the geometrical properties of fractures on the fluid flow properties is conducted and evaluated with respect to their impact on reactive transport modeling. Although employing space-resolved simulation techniques, the results of the tests show average values for permeability and fluid velocity that are comparable to previous studies that used other simulation techniques. Observed fluid flow channeling appears to be related to fracture surfaces matching and anisotropy. However, average flow velocities at potential sites for reactive transport differ up to a factor of five from the average ones for the entire fracture. Furthermore, extreme values at reactive transport sites may differ even more and the flow may be directed against the applied pressure gradient. For studies concerned with simulation of reactive transport, these deviations are crucial and should be explicitly considered in the calculations. Hence space-resolved fluid flow simulations should be employed for the simulation of reactive transport.     

FAE爆炸波对地面目标作用的三维数值研究

利用二阶迎风TVD格式 ,对多组分全N S方程进行全耦合数值求解 ,研究了近地空中环氧丙烷 (C3H6 O)形成的圆柱形云雾爆炸波和地面作用问题 ,得到了对目标作用前、后冲击波的传播和波系结构的演化以及不同时刻爆源组分的分布。对目标表面计算单元的压力进行面积加权平均 ,得到目标各表面平均超压随时间变化的曲线。压力值的动画显示表明了目标表面所受载荷的大小和分布随时间变化的规律。修改初、边值条件 ,本文的方法和开发的软件可对任意形状云雾爆炸波与目标作用问题进行综合数值分析。     

A Direct Numerical Study of Transport and Anisotropy in a Stably Stratified Turbulent Flow with Uniform Horizontal Shear

Direct numerical simulations of homogeneous turbulence in stably stratified shear flow have been performed to aid the understanding of turbulence and turbulent mixing in geophysical flow. Two cases are compared. In the first case, which has been studied in the past, the mean velocity has vertical shear and the mean density is vertically stably stratified. In the second case, which has not been studied systematically before, the mean velocity has horizontal shear and the mean density is again vertically stably stratified. The critical value of the gradient Richardson number, for which a constant turbulence level is obtained, is found to be an order of magnitude larger in the horizontal shear case. The turbulent transport coefficients of momentum and vertical mass transfer are also an order of magnitude larger in the horizontal shear case. The anisotropy of the turbulence intensities are found to be in the range expected of flows with mean shear with no major qualitative change in the range of Richardson numbers studied here. However, the anisotropy of the turbulent dissipation rate is strongly affected by stratification with the vertical component dominating the others. This revised version was published online in July 2006 with corrections to the Cover Date.     

CO2 Injection into Saline Carbonate Aquifer Formations II: Comparison of Numerical Simulations to Experiments

Sequestration of carbon dioxide in geological formations is an alternative way of managing extra carbon. Although there are a number of mathematical modeling studies related to this subject, experimental studies are limited and most studies focus on injection into sandstone reservoirs as opposed to carbonate ones. This study describes a fully coupled geochemical compositional equation-of-state compositional simulator (STARS) for the simulation of CO₂ storage in saline aquifers. STARS models physical phenomena including (1) thermodynamics of sub- and supercritical CO₂, and PVT properties of mixtures of CO₂ with other fluids, including (saline) water; (2) fluid mechanics of single and multiphase flow when CO₂ is injected into aquifers; (3) coupled hydrochemical effects due to interactions between CO₂, reservoir fluids, and primary mineral assemblages; and (4) coupled hydromechanical effects, such as porosity and permeability change due to the aforementioned blocking of pores by carbonate particles and increased fluid pressures from CO₂ injection. Matching computerized tomography monitored laboratory experiments showed the uses of the simulation model. In the simulations dissolution and deposition of calcite as well as adsorption of CO₂ that showed the migration of CO₂ and the dissociation of CO₂ into HCO₃ and its subsequent conversion into carbonate minerals were considered. It was observed that solubility and hydrodynamic storage of CO₂ is larger compared to mineral trapping.     

A Pore-Scale Model for Permeable Biofilm: Numerical Simulations and Laboratory Experiments

Transport in Porous Media - In this paper, we derive a pore-scale model for permeable biofilm formation in a two-dimensional pore. The pore is divided into two phases: water and biofilm. The...     

非稳态不可压缩自由表面粘性流体的三维数值计算

扩展了相对体积算法，计算了变速旋转的敞口圆筒内水的真实非稳态流动。采用了原始变量、交错非等分网格和显式迭代。计算结果与实验现象相吻合。当圆筒长时间等速旋转，其内流体与筒体一起作刚性旋转时，计算自由表面形状与流体力学理论公式的预言吻合得很好。     

Direct Numerical Simulations of Premixed Turbulent Flames with Reduced Chemistry: Validation and Flamelet Analysis

Direct numerical simulation is a very powerful tool to evaluate the validity of new models and theories for turbulent combustion. In this paper, direct numerical simulations of spherically expanding premixed turbulent flames in the thin reaction zone regime and in the broken reaction zone regime are performed. The flamelet-generated manifold method is used in order to deal with detailed reaction kinetics. The computational results are analyzed by using an extended flame stretch theory. It is investigated whether this theory is able to describe the influence of flame stretch and curvature on the local burning velocity of the flame. It is found that if the full profiles of flame stretch and curvature through the flame front are included in the theory, the local mass burning rate is well predicted. The influence of several approximations, which are used in other existing theories, is studied. When flame stretch is assumed constant through the flame front or when curvature of the flame front is neglected, the theory fails to predict the local mass burning rate. The influence of using a reduced chemistry model is investigated by comparing flamelet simulations with reduced and detailed chemistry.