A Robust Bubble Growth Solution Scheme for Implementation in CFD Analysis of Multiphase Flows
Although the full form of the Rayleigh–Plesset (RP) equation more accurately depicts the bubble behavior in a cavitating flow than its reduced form, it finds much less application than the latter in the computational fluid dynamic (CFD) simulation due to its high stiffness. The traditional variable...
| Main Authors: | , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
MDPI AG
2023-03-01
|
| Series: | Computation |
| Subjects: | |
| Online Access: | https://www.mdpi.com/2079-3197/11/4/72 |
| _version_ | 1797605868184797184 |
|---|---|
| author | Hao Pang Gracious Ngaile |
| author_facet | Hao Pang Gracious Ngaile |
| author_sort | Hao Pang |
| collection | DOAJ |
| description | Although the full form of the Rayleigh–Plesset (RP) equation more accurately depicts the bubble behavior in a cavitating flow than its reduced form, it finds much less application than the latter in the computational fluid dynamic (CFD) simulation due to its high stiffness. The traditional variable time-step scheme for the full form RP equation is difficult to be integrated with the CFD program since it requires a tiny time step at the singularity point for convergence and this step size may be incompatible with time marching of conservation equations. This paper presents two stable and efficient numerical solution schemes based on the finite difference method and Euler method so that the full-form RP equation can be better accepted by the CFD program. By employing a truncation bubble radius to approximate the minimum bubble size in the collapse stage, the proposed schemes solve for the bubble radius and wall velocity in an explicit way. The proposed solution schemes are more robust for a wide range of ambient pressure profiles than the traditional schemes and avoid excessive refinement on the time step at the singularity point. Since the proposed solution scheme can calculate the effects of the second-order term, liquid viscosity, and surface tension on the bubble evolution, it provides a more accurate estimation of the wall velocity for the vaporization or condensation rate, which is widely used in the cavitation model in the CFD simulation. The legitimacy of the solution schemes is manifested by the agreement between the results from these schemes and established ones from the literature. The proposed solution schemes are more robust in face of a wide range of ambient pressure profiles. |
| first_indexed | 2024-03-11T05:07:10Z |
| format | Article |
| id | doaj.art-91bf514cb612494684b52cfd7732f5fc |
| institution | Directory Open Access Journal |
| issn | 2079-3197 |
| language | English |
| last_indexed | 2024-03-11T05:07:10Z |
| publishDate | 2023-03-01 |
| publisher | MDPI AG |
| record_format | Article |
| series | Computation |
| spelling | doaj.art-91bf514cb612494684b52cfd7732f5fc2023-11-17T18:49:05ZengMDPI AGComputation2079-31972023-03-011147210.3390/computation11040072A Robust Bubble Growth Solution Scheme for Implementation in CFD Analysis of Multiphase FlowsHao Pang0Gracious Ngaile1Department of Mechanical and Aerospace Engineering, North Carolina State University, Raleigh, NC 27695, USADepartment of Mechanical and Aerospace Engineering, North Carolina State University, Raleigh, NC 27695, USAAlthough the full form of the Rayleigh–Plesset (RP) equation more accurately depicts the bubble behavior in a cavitating flow than its reduced form, it finds much less application than the latter in the computational fluid dynamic (CFD) simulation due to its high stiffness. The traditional variable time-step scheme for the full form RP equation is difficult to be integrated with the CFD program since it requires a tiny time step at the singularity point for convergence and this step size may be incompatible with time marching of conservation equations. This paper presents two stable and efficient numerical solution schemes based on the finite difference method and Euler method so that the full-form RP equation can be better accepted by the CFD program. By employing a truncation bubble radius to approximate the minimum bubble size in the collapse stage, the proposed schemes solve for the bubble radius and wall velocity in an explicit way. The proposed solution schemes are more robust for a wide range of ambient pressure profiles than the traditional schemes and avoid excessive refinement on the time step at the singularity point. Since the proposed solution scheme can calculate the effects of the second-order term, liquid viscosity, and surface tension on the bubble evolution, it provides a more accurate estimation of the wall velocity for the vaporization or condensation rate, which is widely used in the cavitation model in the CFD simulation. The legitimacy of the solution schemes is manifested by the agreement between the results from these schemes and established ones from the literature. The proposed solution schemes are more robust in face of a wide range of ambient pressure profiles.https://www.mdpi.com/2079-3197/11/4/72Rayleigh–Plesset equationfinite differenceeuler methodcavitation |
| spellingShingle | Hao Pang Gracious Ngaile A Robust Bubble Growth Solution Scheme for Implementation in CFD Analysis of Multiphase Flows Computation Rayleigh–Plesset equation finite difference euler method cavitation |
| title | A Robust Bubble Growth Solution Scheme for Implementation in CFD Analysis of Multiphase Flows |
| title_full | A Robust Bubble Growth Solution Scheme for Implementation in CFD Analysis of Multiphase Flows |
| title_fullStr | A Robust Bubble Growth Solution Scheme for Implementation in CFD Analysis of Multiphase Flows |
| title_full_unstemmed | A Robust Bubble Growth Solution Scheme for Implementation in CFD Analysis of Multiphase Flows |
| title_short | A Robust Bubble Growth Solution Scheme for Implementation in CFD Analysis of Multiphase Flows |
| title_sort | robust bubble growth solution scheme for implementation in cfd analysis of multiphase flows |
| topic | Rayleigh–Plesset equation finite difference euler method cavitation |
| url | https://www.mdpi.com/2079-3197/11/4/72 |
| work_keys_str_mv | AT haopang arobustbubblegrowthsolutionschemeforimplementationincfdanalysisofmultiphaseflows AT graciousngaile arobustbubblegrowthsolutionschemeforimplementationincfdanalysisofmultiphaseflows AT haopang robustbubblegrowthsolutionschemeforimplementationincfdanalysisofmultiphaseflows AT graciousngaile robustbubblegrowthsolutionschemeforimplementationincfdanalysisofmultiphaseflows |