Using Vorticity as an Indicator for the Generation of Optimal Coarse Grid Distribution

An improved vorticity-based gridding technique is presented and applied to create optimal non-uniform Cartesian coarse grid for numerical simulation of two-phase flow. The optimal coarse grid distribution (OCGD) is obtained in a manner to capture variations in both permeability and fluid velocity of the fine grid using a single physical quantity called “vorticity”. Only single-phase flow simulation on the fine grid is required to extract the vorticity. Based on the fine-scale vorticity information, several coarse grid models are generated for a given fine grid model. Then the vorticity map preservation error is used to predict how well each coarse grid model reproduces the fine-scale simulation results. The coarse grid model which best preserves the fine-scale vorticity, i.e. has the minimum vorticity map preservation error is recognized as an OCGD. The performance of vorticity-based optimal coarse grid is evaluated for two highly heterogeneous 2D formations. It is also shown that two-phase flow parameters such as mobility ratio have only minor impact on the performance of the predicted OCGD.     

三维非均匀不稳定渗流方程的二维不均匀粗化算法

在无源汇条件下,根据流过某一个横截面的流体流量等于流过这一横截面内所有精细网格的流体流量之和这一特点提出了粗化网格等效渗透率的计算方法。在粗化区内,利用直接解法求解二维渗流方程,再用这些解合成粗化网格的三维合成解,并由合成解计算粗化网格的等效渗透率。根据精度的要求采用了不均匀网格粗化,在流体流速大的区域采用精细网格。利用所得等效渗透率计算了粗化网格的某三维非均匀不稳定渗流场的压降解,结果表明三维非均匀不稳定渗流方程的二维不均匀粗化解非常逼近采用精细网格的解,但计算的速度比采用精细网格提高了80倍。     

The effect of inflow conditions on the transition to turbulence in large eddy simulations of spatially developing mixing layers

Large Eddy Simulations (LES) of spatially developing turbulent mixing layers have been performed for flows of uniform density and Reynolds numbers of up to 50,000 based on the visual thickness of the layer and the velocity difference across it. On a fine LES grid, a validation simulation performed with a hyperbolic tangent inflow profile produces flow statistics that compare extremely well with reference Direct Numerical Simulation (DNS) data. An inflow profile derived from laminar Blasius profiles produces a flow that is significantly different to the reference DNS, particularly with respect to the initial development of the flow. When compared with experimental data, however, it is the boundary layer-type inflow simulation produces the better prediction of the flow statistics, including the mean transition location. It is found that the boundary layer inflow condition is more unstable than the hyperbolic tangent inlet profile. A suitably designed coarse LES grid produces good predictions of the mean transition location with boundary layer inflow conditions at a low computational cost. The results suggest that hyperbolic tangent functions may produce unreliable DNS data when used as the initial condition for studies of the transition in the mixing layer flow.     

Generation of Voronoi Grid Based on Vorticity for Coarse-Scale Modeling of Flow in Heterogeneous Formations

We present a novel unstructured coarse grid generation technique based on vorticity for upscaling two-phase flow in permeable media. In the technique, the fineness of the gridblocks throughout the domain is determined by vorticity distribution such that where the larger is the vorticity at a region, the finer are the gridblocks at that region. Vorticity is obtained from single-phase flow on original fine grid, and is utilized to generate a background grid which stores spacing parameter, and is used to steer generation of triangular and finally Voronoi grids. This technique is applied to two channelized and heterogeneous models and two-phase flow simulations are performed on the generated coarse grids and, the results are compared with the ones of fine scale grid and uniformly gridded coarse models. The results show a close match of unstructured coarse grid flow results with those of fine grid, and substantial accuracy compared to uniformly gridded coarse grid model.     

A Novel Framework for Integration of Random-Walk Particle-Tracking Simulation in Subsurface Multi-phase Flow Modeling

Coarse-scale models are generally preferred in the numerical simulation of multi-phase flow due to computational constraints. However, capturing the effects of fine-scale heterogeneity on flow and isolating the impacts of numerical (artificial) dispersion, which increases with scale, are not trivial. In this paper, a particle-tracking method is devised and integrated in a scale-up workflow to estimate the conditional probability distributions of multi-phase flow functions, which can be considered as inputs in coarse-scale simulations with existing commercial packages. First, a novel particle-tracking method is developed to solve the saturation transport equation. The transport calculation is coupled with a velocity update, following the implicit pressure, explicit saturation framework, to solve the governing equations of two-phase immiscible flow. Each phase particle is advanced in a deterministic convection step according to the phase velocity, as well as in a stochastic dispersion step based on the random Brownian motion. A kernel-based formulation is proposed for computation of fluid saturation in accordance with the phase particle distribution. A novel aspect is that this method employs the kernel approach to construct saturation from phase particle distribution, which is an important improvement to the conventional box method that necessitates a large number of particles per grid cell for consistent saturation interpolation. The model is validated against various analytical solutions. Finally, the validated model is integrated in a statistical scale-up procedure to calibrate effective, or “pseudo,” multi-phase flow functions (e.g., relative permeability functions). The proposed scale-up framework does not impose any length scale requirement regarding the distribution of sub-grid heterogeneities.     

Three-dimensional analysis of the flow in a curved hydraulic turbine draft tube

The three-dimensional turbulent flow in a curved hydraulic turbine draft tube is studied numerically. The analysis is based on the steady Reynolds-averaged Navier–Stokes equations closed with the κ-ε model. The governing equations are discretized by a conservative finite volume formulation on a non-orthogonal body-fitted co-ordinate system. Two grid systems, one with 34 × 16 × 12 nodes and another with 50 × 30 × 22 nodes, have been used and the results from them are compared. In terms of computing effort, the number of iterations needed to yield the same degree of convergence is found to be proportional to the square root of the total number of nodes employed, which is consistent with an earlier study made for two-dimensional flows using the same algorithm. Calculations have been performed over a wide range of inlet swirl, using both the hybrid and second-order upwind schemes on coarse and fine grids. The addition of inlet swirl is found to eliminate the stalling characteristics in the downstream region and modify the behaviour of the flow markedly in the elbow region, thereby affecting the overall pressure recovery noticeably. The recovery factor increases up to a swirl ratio of about 0˙75, and then drops off. Although the general trends obtained with both finite difference operators are in agreement, the quantitative values as well as some of the fine flow structures can differ. Many of the detailed features observed on the fine grid system are smeared out on the coarse grid system, pointing out the necessity of both a good finite difference operator and a good grid distribution for an accurate result.     

基于离散缝洞网络模型的缝洞型油藏混合多尺度有限元数值模拟

针对缝洞型油藏具有多尺度特征,基于离散缝洞网络模型,建立了Darcy/Stokes-Brinkman多尺度耦合数学模型,采用多尺度混合有限元方法,对缝洞型油藏的流体流动问题进行了研究。阐述了多尺度混合有限元方法的基本原理,推导得到了Darcy/Stokes-Brinkman多尺度模型的多尺度混合有限元计算格式。数值计算结果表明,在大尺度模型上进行计算能够捕捉到小尺度上的流动特征;相对于传统有限元,多尺度混合有限元能够捕捉小尺度上的非均质特征而具有更高的计算精度,在保证计算精度的同时能够减少计算量。     

Large-eddy Simulation of Triangular-stabilized Lean Premixed Turbulent Flames: Quality and Error Assessment

In this numerical study, an algebraic flame surface wrinkling (AFSW) reaction submodel based on the progress variable approach is implemented in the large-eddy simulation (LES) context and validated against the triangular stabilized bluff body flame configuration measurements i.e. in VOLVO test rig. The quantitative predictability of the AFSW model is analyzed in comparison with another well validated turbulent flame speed closure (TFC) combustion model in order to help assess the behaviour of the present model and to further help improve the understanding of the flow and flame dynamics. Characterization of non-reacting (or cold) and reacting flows are performed using various subgrid scale models for consistent grid size variation with 300,000 (coarse), 1.2 million (intermediate) and 2.4 million (fine) grid cells. For non-reacting flows at inlet velocity of 17?m/s and inlet temperature 288?K, coarse grid leads to over prediction of turbulence quantities due to low dissipation at the early stage of flow development behind the bluff body that convects downstream eventually polluting the resulting solution. The simulated results with the intermediate (and fine) grid for mean flow and turbulence quantities, and the vortex shedding frequency (f_s) closely match experimental data. For combusting flows for lean propane/air mixtures at 35?m/s and 600?K, the vortex shedding frequency increase threefold compared with cold scenario. The predicted results of mean, rms velocities and reaction progress variable are generally in good agreement with experimental data. For the coarse grid the combustion predictions show a shorter recirculation region due to higher turbulent burning rate. Finally, both cold and reacting LES data are analyzed for uncertainty in the solution using two quality assessment techniques: two-grid estimator by Celik, and model and grid variation by Klein. For both approaches, the resolved turbulent kinetic energy is used to estimate the grid quality and error assessment. The quality assessment reveals that the cold flows are well resolved even on the intermediate mesh, while for the reacting flows even the fine mesh is locally not sufficient in the flamelet region. The Klein approach estimates that depending on the recirculation region in cold scenario both numerical and model errors rise near the bluff-body region, while in combusting flows these errors are significant behind the stabilizing point due to preheating of unburned mixture and reaction heat release. The total error mainly depends on the numerical error and the influence of model error is low for this configuration.     

Compound Wall Treatment for RANS Computation of Complex Turbulent Flows and Heat Transfer

We present a generalised treatment of the wall boundary conditions for RANS computation of turbulent flows and heat transfer. The method blends the integration up to the wall (ItW) with the generalised wall functions (GWF) that include non-equilibrium effects. Wall boundary condition can thus be defined irrespective of whether the wall-nearest grid point lies within the viscous sublayer, in the buffer zone, or in the fully turbulent region. The computations with fine and coarse meshes of a steady and pulsating flow in a plane channel, in flow behind a backward-facing step and in a round impinging jet using the proposed compound wall treatment (CWT) are all in satisfactory agreement with the available experiments and DNS data. The method is recommended for computations of industrial flows in complex domains where it is difficult to generate a computational grid that will satisfy a priori either the ItW or WF prerequisites.     

Accurate Representation of Near-well Effects in Coarse-Scale Models of Primary Oil Production

Near-well effects can have a strong impact on many subsurface flow processes. In oil production, because dissolved gas is released from the oil phase when the pressure falls below the bubble point, the detailed pressure field in the immediate vicinity of a production well strongly impacts gas (and thus oil) production. This effect is complicated by the interplay of fine-scale heterogeneity and two-phase flow physics, and can be difficult to capture in coarse-grid simulations. In this article, we develop and apply a new upscaling (coarse-graining) procedure to capture such near-well subgrid effects in coarse-scale flow simulation models. The method entails the use of preprocessing computations over near-well domains [referred to as local well models (LWM)] for the determination of upscaled single-phase and two-phase near-well parameters. These parameters are computed by minimizing the mismatch between fine and coarse-scale flows over the LWM. Minimization is accomplished using a gradient-based optimization procedure, with gradients calculated through solution of adjoint equations. The boundary conditions applied on the LWM can impact the upscaled parameters, but these boundary conditions depend on the global flow and are not, therefore, known a priori. In order to circumvent this difficulty, an adaptive local–global procedure is applied. This entails performing a global coarse-scale simulation with initial estimates for well-block parameters. The resulting pressure and saturation fields are then used to define local boundary conditions for the near-well computations. The overall procedure is applied to several example problems and is shown to provide results in close agreement with reference fine-scale computations. Significant improvement in accuracy over existing near-well upscaling treatments is demonstrated, particularly for a heavy oil case with oil viscosity of ~10⁴ cp.     

Effect of heat transfer augmentation on two-phase flow instabilities in a vertical boiling channel

The effect of different heater surface configurations on two-phase flow instabilities has been investigated in a single channel, forced convection, open loop, up-flow system. Freon-11 is used as the test fluid, and six different heater tubes with various inside surface configurations have been tested at five different heat inputs. In addition to temperature and pressure recordings, high speed motion pictures of the two-phase flow were taken for some of the experiments to study the two-phase flow behavior at different operating points. Experimental results are shown on system pressure drop versus mass flow rate curves, and stability boundaries are also indicated on these curves. Comparison of different heater tubes is made by the use of the stability boundary maps and the plots of inlet throttling necessary to stabilize the system versus mass flow rate. Tubes with internal springs were found to be more stable than the other tubes.     

Dual Mesh Method for Upscaling in Waterflood Simulation

Detailed geological models typically contain many more cells than can be accommodated by reservoir simulation due to computer time and memory constraints. However, recovery predictions performed on a coarser upscaled mesh are inevitably less accurate than those performed on the initial fine mesh. Recent studies have shown how to use both coarse and fine mesh information during waterflooding simulations. In this paper, we present an extension of the dual mesh method (Verdière and Guérillot, 1996) which simulates water flooding injection using both the coarse and the original fine mesh information. The pressure field is first calculated on the coarse mesh. This information is used to estimate the pressure field within each coarse cell and then phase saturations are updated on the fine mesh. This method avoids the most time consuming step of reservoir simulation, namely solving for the pressure field on the fine grid. A conventional finite difference IMPES scheme is used considering a two phase fluid with gravity and vertical wells. Two upscaling methodologies are used and compared for averaging the coarse grid properties: geometric average and the pressure solve method. A series of test cases show that the method provides predictions similar to those of full fine grid simulations but using less computer time.     

A semi-empirical model of two-phase flow of refrigerant-134a through short tube orifices

Measurements were conducted on Refrigerant-134a flowing through short tube orifices with length-to-diameter (L/D) ratios ranging from 5 to 20. Both two-phase and subcooled liquid flow conditions entering the short tube were examined for upstream pressures ranging from 896 to 1448 kPa and for qualities as high as 10% and subcoolings as high as 13.9°C. Data were analyzed as a function of the main operating variables and tube geometry. Semi-empirical models for both single- and two-phase flow at the inlet of the short tubes were developed to predict the mass flow of Refrigerant-134a through short tube orifices.

Choked flow conditions for Refrigerant-134a were typically established when downstream pressures were reduced below the saturation pressure corresponding to the inlet temperature. The flow rate strongly depended on the upstream pressure and upstream subcooling/quality. The mass flow also depended on cross-sectional area and short tube length. The mass flow model utilized a modified orifice equation that formulated the mass flow as a function of normalized operating variables and short tube geometry. For a two-phase flow entering the short tube, the modified orifice equation was corrected using a theoretically derived expression that related the liquid portion of the mass flow under two-phase conditions to a flow that would occur if the flow were a single-phase liquid. It was found that for sharp-edged short tubes with single- and two-phase flow, approximately 95% of the measured data and model's prediction were within ±15% of each other.     

A multiscale approach to hybrid RANS/LES wall modeling within a high-order discontinuous Galerkin scheme using function enrichment

We present a novel approach to hybrid Reynolds-averaged Navier-Stokes (RANS)/ large eddy simulation (LES) wall modeling based on function enrichment, which overcomes the common problem of the RANS-LES transition and enables coarse meshes near the boundary. While the concept of function enrichment as an efficient discretization technique for turbulent boundary layers has been proposed in an earlier article by Krank & Wall (A new approach to wall modeling in LES of incompressible flow via function enrichment. J Comput Phys. 2016;316:94-116), the contribution of this work is a rigorous derivation of a new multiscale turbulence modeling approach and a corresponding discontinuous Galerkin discretization scheme. In the near-wall area, the Navier-Stokes equations are explicitly solved for an LES and a RANS component in one single equation. This is done by providing the Galerkin method with an independent set of shape functions for each of these two methods; the standard high-order polynomial basis resolves turbulent eddies, where the mesh is sufficiently fine and the enrichment automatically computes the ensemble-averaged flow if the LES mesh is too coarse. As a result of the derivation, the RANS model is applied solely to the RANS degrees of freedom, which effectively prevents the typical issue of a log-layer mismatch in attached boundary layers. As the full Navier-Stokes equations are solved in the boundary layer, spatial refinement gradually yields wall-resolved LES with exact boundary conditions. Numerical tests show the outstanding characteristics of the wall model regarding grid independence, superiority compared to equilibrium wall models in separated flows, and achieve a speed-up by two orders of magnitude compared to wall-resolved LES.     

A local adaptive grid procedure for incompressible flows with multigridding and equidistribution concepts

An adaptive grid solution procedure is developed for incompressible flow problems in which grid refinement based on an equidistribution law is performed in high-error-estimate regions that are flagged from a preliminary coarse grid solution. Solutions on the locally refined and equidistributed meshes are obtained using boundary conditions interpolated from the preliminary coarse grid solution, and solutions on both the refined and coarse grid regions are successively improved using a multigrid approach. For this purpose, suitable correction terms for the coarse grid equations are derived for all variables in the flagged regions. This procedure with Local Adaptation, Multigridding and Equidistribution (LAME) concepts is applied to various flow problems to demonstrate the accuracy improvements obtained using this method.     

Improved subgrid scale model for dense turbulent solid-liquid two-phase flows

The dense solid-phase governing equations for two-phase flows are obtained by using the kinetic theory of gas molecules. Assuming that the solid-phase velocity distributions obey the Maxwell equations, the collision term for particles under dense two-phase flow conditions is also derived. In comparison with the governing equations of a dilute two-phase flow, the solid-particle‘s governing equations are developed for a dense turbulent solid-liquid flow by adopting some relevant terms from the dilute two-phase governing equations. Based on Cauchy-Helmholtz theorem and Smagorinsky model, a second-order dynamic sub-grid-scale (SGS) model, in which the sub-grid-scale stress is a function of both the strain-rate tensor and the rotation-rate tensor, is proposed to model the two-phase governing equations by applying dimension analyses. Applying the SIMPLEC algorithm and staggering grid system to the two-phase discretized governing equations and employing the slip boundary conditions on the walls, the velocity and pressure fields, and the volumetric concentration are calculated. The simulation results are in a fairly good agreement with experimental data in two operating cases in a conduit with a rectangular cross-section and these comparisons imply that these models are practical.     

Multiscale simulation of flow and heat transfer of nanofluid with lattice Boltzmann method

Based on the lattice Boltzmann (LB) approach, a novel hybrid method has been proposed for getting insight into the microscale characteristics of the multicomponent flow of nanofluid. In this method, the whole computational domain is divided into two regions in which different-sized meshes are involved for simulation (fine mesh and coarse mesh). The multicomponent LB method is adopted in the fine mesh region, and the single-component LB approach is applied to the coarse mesh region where the nanofluid is treated as a mixed single-component fluid. The conservation principles of mass, momentum and energy are used to derive a hybrid scheme across the different scaled regions. Numerical simulation is carried out for the Couette flow and convective heat transfer in a parallel plate channel to validate the hybrid method. The computational results indicate that by means of the present method, not only the microscopic characteristics of the nanofluid flow can be simulated, but also the computational efficiency can be remarkably improved compared with the pure multicomponent LB method.     

Fractional Flow Theory Applicable to Non-Newtonian Behavior in EOR Processes

The method of characteristics, or fractional-flow theory, is extremely useful in understanding complex Enhanced Oil Recovery (EOR) processes and in calibrating simulators. One limitation has been its restriction to Newtonian rheology except in rectilinear flow. Its inability to deal with non-Newtonian rheology in polymer and foam EOR has been a serious limitation. We extend fractional flow methods for two-phase flow to non-Newtonian fluids in one-dimensional cylindrical flow, where rheology changes with distance from injection well. The fractional flow curve is then a function of position and we analyze the characteristic equations for two applications—polymer and foam floods. For polymer flooding, we present a semi-analytical solution for the changing fractional flow curve where characteristics and shocks collide. The semi-analytical solution is shown to give good agreement with the finite-difference simulation thus helping us understand the development and resolution of shocks. We discuss two separate cases of foam injection with or without preflush. We observe that the fractional flow solutions are more accurate than finite-difference simulations on a comparable grid and hence the method can be used to calibrate simulators. For SAG (alternating-slug) foam injection, characteristics and shocks collide, making the fractional-flow solution complex. Nonetheless, one can solve exactly for changing mobility near the well, to greater accuracy than with conventional simulation. The fractional-flow method extended to non-Newtonian flow can be useful both for its insights for scale-up of laboratory experiments and to calibrate computer simulators involving non-Newtonian EOR. It can also be an input to streamline simulations.     

Application of Effective Flux Boundary Conditions to Two-Phase Upscaling in Porous Media

A new algorithm is introduced for upscaling relative permeabilities, and tested in simulations of two-dimensional reservoir displacement processes. The algorithm is similar to existing algorithms for computing upscaled relative permeabilities from subgrid simulations, but uses new boundary conditions for the pressure field. The new 'effective flux boundary conditions' were introduced in a previous paper and provide a more accurate estimate of flux through high permeability channels. The algorithm was tested in conjunction with uniform grid coarsening and upscaled absolute permeabilities for a broad range of coarsenings. The permeability fields were highly heteroge-neous and layered, and were obtained from synthetic data and from conditioned realizations of actual oil reservoirs. The algorithm was tested for a wide variety of grid aspect ratios, and for both viscous-and gravity-dominated flow. Typical fine grids were of the order of 100×100 cells; the coarsest scaled-up grids were on the order of 5×5 cells. The quality of scale up was evaluated by comparing oil cut curves for the fine and coarse grid simulations. We consistently obtained excellent agreement, even at the coarsest levels of scale up.