Search results

This paper investigates the multi‐phase behaviour of droplets injected into a nozzle at two separate wall locations. The physical features of the droplets (rate of mass, density and radius) at each injector location are identical. This system can be described by a two‐phase Eulerian—Eulerian approach that yields classical systems of equations: three for the gaseous phase and three for the dispersed droplet phase. An underlying assumption in the two phase model is that no interaction occurs between droplets. The numerical solution of the model (using the MacCormack scheme) indicates however that the opposite jets do interact to form one jet. This inconsistency is overcome in the current paper by associating the droplets from a given injection location with a separate phase and subsequently solving equations describing a multiphase system (here, three‐phase system). Comparison of numerical predications between the two‐phase and the multiphase model shows significantly different results. In particular the multiphase model shows no jet interaction.

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 a finite‐element analysis of the initiation of a slope failure in a small‐scale laboratory test due to pore pressure variation. To this aim, a fully coupled multiphase model for saturated/partially saturated solid porous materials based on porous media mechanics is used.

The slope is described as a three‐phase deforming porous continuum where heat, water and gas flow are taken into account. The gas phase is modelled as an ideal gas composed of dry air and water vapour. Phase changes of water, heat transfer through conduction and convection and latent heat transfer are considered. The independent variables are: solid displacements, capillary pressure, gas pressure and temperature. The effective stress state is limited by Drucker‐Prager yield surface for the sake of simplicity. Small strains and quasi‐static loading conditions are assumed.

The paper shows that the multiphase modelling is able to capture the main experimental observations such as the local failure zone at the onset of slope failure and the outflow appeared in that zone. It also allows understanding of the triggering mechanisms of the failure zone.

This work can be considered as a step towards a further development of a suitable numerical model for the simulation of non‐isothermal geo‐environmental engineering problems.

The multiphysics approach looks promising for the analysis of the onset of landslides, provided that the constitutive models for the multiphase porous media in saturated/unsaturated conditions and the related mechanical and hydraulic properties are described with sufficient accuracy.

Elasto‐plastic thermo‐hydro‐mechanical modelling of the initiation of slope failure subjected to variation in pore pressure boundary condition.

This paper aims to present an efficient IMPES algorithm based on a global model order reduction method, proper orthogonal decomposition (POD), to achieve the fast solution and prediction of two-phase flows in porous media.

The key point of the proposed algorithm is to establish an accurate POD reduced-order model (ROM) for two-phase porous flows. To this end, two projection methods including projecting the original governing equations (Method I) and projecting the discrete form of original governing equations (Method II) are respectively applied to construct the POD-ROM, and their distinctions are compared and analyzed in detail. It is found the POD-ROM established by Method I is inapplicable to multiphase porous flows due to its failed introduction of fluid saturation and permeability that locate on the edge of grid cell, which would lead to unphysical results.

By using Method II, an efficient IMPES algorithm that can substantially speed up the simulation of two-phase porous flows is developed based on the POD-ROM. The computational efficiency and numerical accuracy of the proposed algorithm are validated through three numerical examples, and simulation results illustrate that the proposed algorithm displays satisfactory computational speed-up (one to two orders of magnitude) without sacrificing numerical accuracy obviously when comparing to the standard IMPES algorithm that without any acceleration technique. In addition, the determination of POD modes number, the relative errors of wetting phase pressure and saturation, and the influence of POD modes number on the overall performances of the proposed algorithm, are investigated.

1. Two projection methods are applied to establish the POD-ROM for two-phase porous flows and their distinctions are analyzed. The reason why POD-ROM is difficult to be applied to multiphase porous flows is clarified firstly in this study. 2. A highly efficient IMPES algorithm based on the POD-ROM is proposed to accelerate the simulation of two-phase porous flows. 3. Satisfactory computational speed-up (one to two orders of magnitude) and prediction accuracy of the proposed algorithm are observed under different conditions.

The purpose of this paper is to the study the multiphase bubbles flow motion in a vertical channel with an electroconducting liquid without and under the influence of a magnetic field.

For numerical calculations, the lattice Boltzmann method (LBM) is used, which is based on the kinetic theory for solving fluid mechanics and other physical problems. The phase-field lattice Boltzmann model is developed to simulate the behaviour of multiphase bubble–bubble interaction while rising in the fluid with high density ratios.

The behaviour of the rising bubble flow in a rectangular column of two phases is investigated with the two-dimensional LBM.

The multiphase flow in electroconducting liquids with high ratio of density is studied using the LBM.

The purpose of this paper is to develop a coupled lattice Boltzmann model for the simulation of the freezing process in unsaturated porous media.

In the developed model, the porous structure with complexity and disorder was generated by using a stochastic growth method, and then the Shan-Chen multiphase model and enthalpy-based phase change model were coupled by introducing a freezing interface force to describe the variation of phase interface. The pore size of porous media in freezing process was considered as an influential factor to phase transition temperature, and the variation of the interfacial force formed with phase change on the interface was described.

The larger porosity (0.2 and 0.8) will enlarge the unfrozen area from 42 mm to 70 mm, and the rest space of porous medium was occupied by the solid particles. The larger specific surface area (0.168 and 0.315) has a more fluctuated volume fraction distribution.

The concept of interfacial force was first introduced in the solid–liquid phase transition to describe the freezing process of frozen soil, enabling the formulation of a distribution equation based on enthalpy to depict the changes in the water film. The increased interfacial force serves to diminish ice formation and effectively absorb air during the freezing process. A greater surface area enhances the ability to counteract liquid migration.

This paper seeks to discuss a mechanistic modeling concept for local phenomena governing two‐ and multi‐phase flows and heat transfer.

An overview is given of selected issues concerning the formulation of multidimensional models of two‐phase flow and heat transfer. A complete computational multiphase fluid dynamics (CMFD) model of two‐phase flow is presented, including local constitutive models applicable to two‐phase flows in heated channels. Results are shown of model testing and validation.

It has been demonstrated that the overall model is capable of capturing various local flow and heat transfer phenomena in general, and the onset of temperature excursion (CHF) in low quality forced‐convection boiling, in particular.

Whereas the multiphase model formulation is applicable to a large class of problems, geometries and operating conditions, the closure laws and results are focused on forced‐convection boiling in heated channels.

The proposed approach can be used to predict multidimensional velocity field and phase distribution in two‐phase flow devices and components used in thermal power plants, nuclear power plants and chemical processing plants.

A complete mechanistic multidimensional model of forced‐convection boiling in heated channels is given. The potential of a CMFD approach is demonstrated to perform virtual experiments that can be used in system design and optimization, and in safety analysis.

In plethora of petroleum, chemical and heat transfer applications, T-junction is often used to partially separate gas from other fluids, to reduce work burden on other separating equipment. The abundance of liquid carryovers from the T-junction side arm is the cause of production downtime in terms of frequent tripping of downstream equipment train. Literature review revealed that regular and reduced T-junctions either have high peak liquid carryovers (PLCs) or the liquid appears early in the side arm [liquid carryover threshold (LCT)]. The purpose of this study is to harvest the useful features of regular and reduced T-junction and analyze diverging T-junction having upstream and downstream pipes.

Volume of fluid as a multiphase model, available in ANSYS Fluent, was used to simulate air–water slug flow in five diverging T-junctions for eight distinct velocity ratios. PLCs and LCT were chosen as key performance indices.

The results indicated that T (0.5–1) and (0.8–1) performed better as low liquid carryovers and high LCT were achieved having separation efficiencies of 96% and 94.5%, respectively. These two diverging T-junctions had significantly lower PLCs and high LCT when compared to other three T-junctions. Results showed that the sudden reduction in the side arm diameter results in high liquid carryovers and lower LCT. Low water and air superficial velocities tend to have low PLC and high LCT.

This study involved working fluids air and water but applies to other types of fluids as well.

The novel T-junction design introduced in this study has significantly higher LCT and lower PLC. This is an indication of higher phase separation performance as compared to other types of T-junctions. Because of lower liquid take-offs, there will be less frequent downstream equipment tripping resulting in lower maintenance costs. Empirical correlations presented in this study can predict fraction of gas and liquid in the side arm without having to repeat the experiment.

Maintenance costs and production downtime can be significantly reduced with the implication of diverging T-junction design.

The presented study revealed that the diameter ratio has a significant impact on PLC and LCT. It can be concluded that novel T-junction designs, T2 and T3, achieved high phase separation; therefore, it is favorable to use in the industry. Furthermore, a few limitations in terms of diameter ratio are also discussed in detail.

The purpose of this paper is to test the feasibility of lattice Boltzmann method (LBM) for numerical simulation of nucleate boiling and transition boiling. In addition, the processes of nucleate and transition boiling on vertical wall are simulated. The heat transfer mechanism is discussed based on the evolution of temperature field.

In this paper, nucleate boiling and transition boiling are numerically investigated by LBM. A lattice Boltzmann (LB) multiphase model combining with a LB thermal model is used to predict the phase-change process.

Numerical results are in good agreement with existing experimental results. Numerical results confirm the feasibility of the hybrid LBM for direct simulations of nucleate and transition boiling. The data exhibit correct parametric dependencies of bubble departure diameter compared with experimental correlation and relevant references.

All the simulations are performed in two-dimensions in this paper. In the future work, the boiling process will be simulated in three-dimensional.

This study demonstrated a potential model that can be applied to the investigation of phase change heat transfer, which is one of the effective techniques for enhance the heat transfer in engineering. The numerical results can be considered as a basic work or a reference for generalizing LB method in the practical application about nucleate boiling and transition boiling.

The hybrid LBM is first used for simulation of nucleate and transition boiling on vertical surface. Heat transfer mechanism during boiling is discussed based on the numerical results.

In this work, a new two-phase version of the finite element-based Artificial Compressibility (AC) Characteristic-Based Split (CBS) algorithm is developed and applied for the first time to heat and mass transfer phenomena in porous media with associated phase change. The purpose of this study is to provide an alternative for the theoretical analysis and numerical simulation of multiphase transport phenomena in porous media. Traditionally, the more complex Separate Flow Model was used in which the vapour and liquid phases were considered as distinct fluids and mathematically described by the conservation laws for each phase separately, resulting in a large number of governing equations.

Even though the adopted mathematical model presents analogies with the conventional multicomponent mixture flow model, it is characterized by a considerable reduction in the number of the differential equations for the primary variables. The fixed-grid numerical formulation can be applied to the resolution of general problems that may simultaneously include a superheated vapour region, a two-phase zone and a sub-cooled liquid region in a single physical domain with irregular and moving phase interfaces in between. The local thermal non-equilibrium model is introduced to consider the heat exchange between fluid and solid within the porous matrix.

The numerical model is verified considering the transport phenomena in a homogenous and isotropic porous medium in which water is injected from one side and heated from the other side, where it leaves the computational domain in a superheated vapour state. Dominant forces are represented by capillary interactions and two-phase heat conduction. The obtained results have been compared with the numerical data available in the scientific literature.

The present algorithm provides a powerful routine tool for the numerical modelling of complex two-phase transport processes in porous media.

For the first time, the stabilized AC-CBS scheme is applied to the resolution of compressible viscous flow transport in porous materials with associated phase change. A properly stabilized matrix inversion-free procedure employs an adaptive local time step that allows acceleration of the solution process even in the presence of large source terms and low diffusion coefficients values (near the phase change point).