Phase Field Method for Simulation of Multiphase Flow

Author(s):  
Harish Ganapathy ◽  
Ebrahim Al-Hajri ◽  
Michael M. Ohadi

The present paper reports a comprehensive study on the numerical simulation of Taylor flow in microchannels by the phase field method. Additionally, a comparative study was also performed against an alternative volume of fluid model based on which the phase field method was found to be more advantageous in key aspects such as the absence of unphysical interfacial pressure oscillations and the ability to account for variations in the surface tension force and thus predict several bubble lengths under constant flow conditions while observing the physics of homogeneous two-phase flow. Different bubble formation mechanisms were simulated and compared against experimental findings in literature. The simulation of a thin liquid film at the channel wall was found to be a limitation of most works pertaining to Taylor flow, including the present. This was ascribed to be more likely due to limited dimensional and spatial resolution as well as inaccurate contact angle dynamics rather than limitations of the modeling approach itself. The effect of wall adhesion was studied with respect to the flow and pressure field in the channel. A validation of the model was achieved through a favorable comparison of the numerically predicted gas void fraction and bubble lengths with existing models and correlations. On the whole, the phase field method was concluded to have improved predictive accuracy with respect to certain aspects as compared to conventional multiphase flow models.

Author(s):  
Harish Ganapathy ◽  
Ebrahim Al-Hajri ◽  
Michael M. Ohadi

The present paper reports a comprehensive study on the numerical simulation of Taylor flow in microchannels by the phase field method. Additionally, a comparative study was also performed against an alternative volume of fluid model based on which the phase field method was found to be more advantageous in key aspects such as the absence of unphysical interfacial pressure oscillations and the ability to account for variations in the surface tension force and thus predict several bubble lengths under constant flow conditions while observing the physics of homogeneous two-phase flow. Different bubble formation mechanisms were simulated and compared against experimental findings in literature. The simulation of a thin liquid film at the channel wall was found to be limitation of most works pertaining to Taylor flow, including the present. This was ascribed to be more likely due to limited dimensional and spatial resolution as well as inaccurate contact angle dynamics rather than limitations of the modeling approach itself. The effect of wall adhesion was studied with respect to the flow and pressure field in the channel. A validation of the model was achieved through a favorable comparison of the numerically predicted gas void fraction and bubble lengths with existing models and correlations. On the whole, the phase field method was concluded to have improved predictive accuracy with respect to certain aspects as compared to conventional multiphase flow models.


Author(s):  
Akinori Tamura ◽  
Kenichi Katono

Abstract Two-phase flows including a phase change such as liquid-vapor flows play an important role in many industrial applications. A deeper understanding of the phase change phenomena is required to improve performance and safety of nuclear power plants. For this purpose, we developed a phase change simulation method based on the phase field method (PFM). Low computational efficiency of the conventional PFM based on the Cahn-Hilliard equation is an obstacle in practical simulations. To resolve this problem, we presented a new PFM based on the conservative Allen-Cahn equation including a phase change model. The wettability also needs to be considered in the phase change simulation. When we apply the conventional wetting boundary condition to the conservative Allen-Cahn equation, there is a problem that the mass of each phase is not conserved on the boundary. To resolve this issue, we developed the mass correction method which enables mass conservation in the wetting boundary. The proposed PFM was validated in benchmark problems. The results agreed well with the theoretical solution and other simulation results, and we confirmed that this PFM is applicable to the two-phase flow simulation including the phase change. We also investigated the computational efficiency of the PFM. In a comparison with the conventional PFM, we found that our proposed PFM was more than 100 times faster. Since computational efficiency is an important factor in practical simulations, the proposed PFM will be preferable in many industrial simulations.


Author(s):  
Naoki Takada

For interface-tracking simulation of two-phase flows in various micro-fluidics devices, the applicability of two versions of Navier-Stokes phase-field method (NS-PFM) was examined, combining NS equations for a continuous fluid with a diffuse-interface model based on the van der Waals-Cahn-Hilliard free-energy theory. Through the numerical simulations, the following major findings were obtained: (1) The first version of NS-PFM gives good predictions of interfacial shapes and motions in an incompressible, isothermal two-phase fluid with high density ratio on solid surface with heterogeneous wettability. (2) The second version successfully captures liquid-vapor motions with heat and mass transfer across interfaces in phase change of a non-ideal fluid around the critical point.


2018 ◽  
Vol 108 ◽  
pp. 169-180 ◽  
Author(s):  
Nazia Talat ◽  
Boštjan Mavrič ◽  
Grega Belšak ◽  
Vanja Hatić ◽  
Saša Bajt ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document