Search results

The purpose of this paper is to develop a general numerical solution for the wetting fluid spread into porous media that can be used in solving of droplet spread into soils, printing applications, fuel cells, composite processing.

A discrete capillary network model based on micro‐force balance is numerically implemented and the flow for an arbitrary capillary number can be solved. At the fluid interface, the boundary condition that accounts for the capillary pressure jump is used.

The wetting fluid spread into porous medium starts as a single‐phase flow, and after some particular number of the porous medium characteristic length scales, the multi‐phase flow pattern occurs. Hence, in the principal flow direction, the phase content (saturation) decreases, and in the lower limit for the capillary number sufficiently small, the saturation should become constant. This qualitative saturation behavior is observed irrespective of the flow dimensionality, whereas the quantitative results vary for different flow systems.

The numerical solution has to be expanded to solve the spread of the fluid in the porous medium after there is no free fluid left at the porous medium surface.

It is shown that the multi‐phase flow can develop even on a small domain due to the porous medium heterogeneity. Neglecting the medium heterogeneity and flow type can lead to a large error as shown for the droplet spread time in the porous medium.

This is believe to be the only paper relating to solving the droplet spread into porous medium as a multi‐phase flow problem.

Lattice Boltzmann method (LBM) is employed to explore friction factor of single-phase fluid flow through porous media and the effects of local porous structure including geometry of grains in porous media and specific surface of porous media on two-phase flow dynamic behavior, phase distribution and relative permeability. The paper aims to discuss this issue.

The 3D single-phase LBM model and the 2D multi-component multi-phase Shan-Chen LBM model (S-C model) are developed for fluid flow through porous media. For the solid site, the bounce back scheme is used with non-slip boundary condition.

The predicted friction factor for single-phase fluid flow agrees well with experimental data and the well-known correlation. Compared with porous media with square grains, the two-phase fluids in porous media with circle grains are more connected and continuous, and consequently the relative permeability is higher. As for the factor of specific porous media surface, the relative permeability of wetting fluids varies a little in two systems with different specific surface areas. In addition, the relative permeability of non-wetting fluid decreases with the increasing of specific surface of porous media due to the large flow resistance.

Fluid-fluid interaction and fluid-solid interaction in the SC LBM model are presented, and schemes to obtain immiscible two-phase flow and different contact angles are discussed. Two-off mechanisms acting on the wetting fluids is proposed to illustrate the relative permeability of wetting fluids varies a little in two systems with different specific surface.

The purpose of this paper is to study the multi-phase flow behaviors in solid fluidization exploitation of natural gas hydrate (NGH) and its effect on the engineering safety.

In this paper, a multi-phase flow model considering the endothermic decomposition of hydrate is established and finite difference method is used to solve the mathematical model. The model is validated by reproducing the field test data of a well in Shenhu Sea area. Besides, optimization of design parameters is presented to ensure engineering safety during the solid fluidization exploitation of NGH in South China Sea.

To ensure the engineering safety during solid fluidization exploitation of marine NGH, taking the test well as an example, a drilling flow rate range of 40–50 L/s, drilling fluid density range of 1.2–1.23 g/cm3 and rate of penetration (ROP) range of 10–20 m/h should be recommended. Besides, pre-cooled drilling fluid is also helpful for inhibiting hydrate decomposition.

Systematic research on the effect of multiphase flow behaviors on the engineering safety is scare, especially for the solid fluidization exploitation of NGH in South China Sea. With the growing demand for energy, it is of great significance to ensure the engineering safety before the large-scale extraction of commercial gas from hydrate deposits. The result of this study can provide profound theoretical bases and valuable technical guidance for the commercial solid fluidization exploitation of NGH in South China Sea.

This paper aims to present a novel numerical technique for solving the incompressible multiphase mixture model.

The multiphase mixture model contains a set of momentum and continuity equations for the mixture phase, a second phase continuity equation and the algebraic equation for the relative velocity. For solving continuity equation for the second phase and advection term of momentum, an improved approach fine grid advection-multiphase mixture flow (FGA-MMF) is developed. In the FGA-MMF method, the continuity equation for the second phase is solved with higher-order schemes in a two times finer grid. To solve the advection term of the momentum equation, the advection fluxes of the volume fraction in the continuity equation for the second phase are used.

This approach has been used in various tests to simulate unsteady flow problems. Comparison between numerical results and experimental data demonstrates a satisfactory performance. Numerical examples show that this approach increases the accuracy and stability of the solution and decreases non-monotonic results.

The solver for the multi-phase mixture model can only be adopted to solve the incompressible fluid flow.

The paper developed an innovative solution (FGA-MMF) to find multi-phase flow field value in the multi-phase mixture model. Advantages of the FGA-MMF technique are the ability to accurately determine the phases interpenetrating, decreasing the numerical diffusion of the interface and preventing instability and non-monotonicity in solution of large density variation problems.

The purpose of this paper is to present detailed algorithms for simulation of individual and group control of production wells in hydrocarbon reservoirs which are implemented in a finite volume-based reservoir simulator.

The algorithm for individual control is described for the multi-lateral multi-connection ones based on the multi-segment model considering cross-flow. Moreover, a general group control algorithm is proposed which can be coupled with any well model that can handle a constraint and returns the flow rates. The performance of oil production process based on the group control criteria is investigated and compared for various cases.

The proposed algorithm for group control of production wells is a non-optimization iterative scheme converging within a few number of iterations. The numerical results of many computer runs indicate that the nominal power of the production wells, in general, is the best group control criterion for the proposed algorithm. The production well group control with a proper criterion can generally improve the oil recovery process at negligible computational costs when compared with individual control of production wells.

Although the group control algorithm is implemented for both production and injection wells in the developed simulator, the numerical algorithm is here described only for production wells to provide more details.

The proposed algorithm can be coupled with any well model providing the fluid flow rates and can be efficiently used for group control of production wells. In addition, the calculated flow rates of the production wells based on the group control algorithm can be used as candidate solutions for the optimizer in the simulation-optimization models. It may reduce the total number of iterations and consequently the computational cost of the simulation-optimization models for the well control problem.

A complete and detailed description of ingredients of an efficient well group control algorithm for the hydrocarbon reservoir is presented. Five group control criteria are extracted from the physical, geometrical and operating conditions of the wells/reservoir. These are the target rate, weighted potential, ultimate rate and introduced nominal power of the production wells. The performance of the group control of production wells with different group control criteria is compared in three different oil production scenarios from a black-oil and highly heterogeneous reservoir.

Natural gas hydrate (NGH) has been regarded as one of the most important resources due to NGH's large amounts of reserve. However, NGH development still faces many technical challenges, such as low production rate and reservoir instability resulting from NGH decomposition. Therefore, developing a fully coupled THMC model for simulating the hydrate decomposition and studying its mechanical behavior is very important and necessary. The purpose of this article is to develop and solve a multi-phase, strong nonlinearity and large-scale fully coupled thermal-hydro-mechanical–chemical (THMC) model for simulating the multi-physics processes involving solid-liquid-gas flow, heat transfer, NGH phase change and rock deformation during NGH decomposition.

In this paper, a multi-phase, strong nonlinearity and large-scale fully coupled THMC model is developed for simulating the multi-physics processes involving solid-liquid-gas flow, heat transfer, NGH phase change and rock deformation during NGH dissociation. The fully coupled THMC model is solved by using a fully implicit finite element method, in which the gas pressure, water pressure, temperature and displacement are taken as basic unknown variables. The proposed model is validated against with the experimental data, showing high accuracy and reliability.

A multi-phase, strong nonlinearity and large-scale fully coupled THMC model is developed for simulating the multi-physics processes involving solid-liquid-gas flow, heat transfer, NGH phase change and rock deformation during NGH decomposition. The proposed model is validated against with the experimental data, showing high accuracy and reliability.

Some assumptions are made to make the model tractable, including (1) the composition gas of hydrate is pure methane; (2) the gas-liquid multi-phase flow in the pore obeys Darcy's law; (3) hydrate occurs on the surface of soil particles, both of them form the composite consolidation material; (4) the small-strain assumption is applied to composite solid materials, which are treated as skeletons and cannot be moved; (5) momentum change caused by phase change is not considered.

NGH has been regarded as one of the most important resources due to its large amounts of reserve. However, NGH development still faces many technical challenges, such as low production rate and reservoir instability resulting from NGH decomposition. Most of the existing studies decouple the process with solid deformation and seepage behavior, but the accuracy of the numerical results will be sacrificed to certain extent. Therefore, it is very important and necessary to develop a fully coupled THMC model for simulating the hydrate decomposition and studying its mechanical behavior.

NGH, widely distributed in shallow seabed or permanent frozen region, has the characteristics of high energy density and high combustion efficiency (Yan et al., 2020). A total of around 7.5 × 1,018 m³ has been proved to exist around the world and 1 m³ of NGH can release about 160–180 m³ of natural gas (Kvenvolden and Lorenson) under normal conditions. Safely and sustainably extracting NGH commercially can effectively relieve global energy pressure and contribute to achieving carbon reduction goals.

The novelty of the present work lies in mainly two aspects. First, a fully coupled THMC model is developed for studying the multi-physics processes involving solid-liquid-gas flow, heat transfer, NGH phase change and solid deformation during NGH dissociation. Second, the numerical solution is obtained by using a fully implicit finite element method (FEM) and is validated against experimental data.

The purpose of this paper is to present a new numerical approach for modeling the multi‐phase flow during an alloy solidification process. In many solidification processes, advection of solid may have a dramatic effect on bulk convection field as well as on the solid front growth and hence on the macro‐segregation pattern. In the present work, a numerical model is developed to simulate directional solidification in presence of melt convection as well as solid advection in the form of sedimentation. A 2D cavity filled with hyper‐eutectic aqueous ammonium chloride solution (25 wt.% of ammonium chloride) being chilled from one of the side walls has been chosen as the model problem for the numerical simulation.

A fixed grid volume averaging technique has been used for solving mass, momentum, energy, and species equation while taking into account the solid phase advection and local re‐melting. Two different criteria have been identified for the solid particles in the mushy zone to be mobile. These two criteria are represented by a critical solid fraction, and a critical velocity. Based on these two criteria, the mushy zone has been subdivided into two different regions namely, an immobile coherent zone consisting of packed equiaxed crystals and a mobile non‐coherent zone where the solid crystals are able to move.

The numerical results are compared with corresponding experimental observations.

The solid advection velocity and source terms dealing with solid velocity have been calculated using an explicit scheme, whereas the main conservation equations are solved using an implicit scheme.

This paper aims to present the implementation of a multi-layer radiation propagation model in simulations of multi-phase flow and heat transfer, for a dissociating methane hydrate reservoir subjected to microwave heating.

To model the induced heterogeneity due to dissociation of hydrates in the reservoir, a multiple homogeneous layer approach, used in food processes modelling, is suggested. The multi-layer model is incorporated in an in-house, multi-phase, multi-component hydrate dissociation simulator based on the finite volume method. The modified simulator is validated with standard experimental results in the literature and subsequently applied to a hydrate reservoir to study the effect of water content and sand dielectric nature on radiation propagation and hydrate dissociation.

The comparison of the multi-layer model with experimental results show a maximum difference in temperature estimation to be less than 2.5 K. For reservoir scale simulations, three homogeneous layers are observed to be sufficient to model the induced heterogeneity. There is a significant contribution of dielectric properties of sediments and water content of the reservoir in microwave radiation attenuation and overall hydrate dissociation. A high saturation reservoir may not always provide high gas recovery by dissociation of hydrates in the case of microwave heating.

The multi-layer approach to model microwave radiation propagation is introduced and tested for the first time in dissociating hydrate reservoirs. The multi-layer model provides better control over reservoir heterogeneity and interface conditions compared to existing homogeneous models.

Modeling of multi-phase flows for Rayleigh-Taylor instability and natural convection in a square cavity has been investigated using an incompressible smoothed particle hydrodynamics (ISPH) technique. In this technique, incompressibility is enforced by using SPH projection method and a stabilized incompressible SPH method by relaxing the density invariance condition is applied. The paper aims to discuss these issues.

The Rayleigh-Taylor instability is introduced in two and three phases by using ISPH method. The author simulated natural convection in a square/cubic cavity using ISPH method in two and three dimensions. The solutions represented in temperature, vertical velocity and horizontal velocity have been studied with different values of Rayleigh number Ra parameter (103=Ra=105). In addition, characteristic based scheme in Finite Element Method is introduced for modeling the natural convection in a square cavity.

The results for Rayleigh-Taylor instability and natural convection flow had been compared with the previous researches.

Modeling of multi-phase flows for Rayleigh-Taylor instability and natural convection in a square cavity has been investigated using an ISPH technique. In ISPH method, incompressibility is enforced by using SPH projection method and a stabilized incompressible SPH method by relaxing the density invariance condition is introduced. The Rayleigh-Taylor instability is introduced in two and three phases by using ISPH method. The author simulated natural convection in a square/cubic cavity using ISPH method in two and three dimensions.

The purpose of this paper is the development of a new turbocharger compressor is a challenging task particularly when both wider operating range and higher efficiency are required. However, the cumbersome design effort and the inherent calculus burden can be significantly reduced by using appropriate design optimization approaches as an alternative to conventional design techniques.

This paper presents an optimization-based preliminary-design (OPD) approach based on a judicious coupling between evolutionary optimization techniques and a modified one-dimensional mean-line model. Two optimization strategies are considered. The first one is mono-objective and is solved using genetic algorithms. The second one is multi-objective and it is handled using the non-dominated sorting genetic algorithm-II. The proposed approach constitutes an automatic search process to select the geometrical parameters of the compressor, ensuring the most common requirements of the preliminary-design phase, with a minimum involvement of the designer.

The obtained numerical results demonstrate that the proposed tool can rapidly produce nearly optimal designs as an excellent basis for further refinement in the phase by using more complex analysis methods such as computational fluid dynamics and meta-modeling.

This paper outlines a new fast OBPD approach for centrifugal compressor turbochargers. The proposal adopts an inverse design method and consists of two main phases: a formulation phase and a solution phase. The complexity of the formulated problem is reduced by using a sensitivity analysis. The solution phase requires to link, in an automatic way, three processes, namely, optimization, design and analysis.