{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,10,11]],"date-time":"2025-10-11T02:01:52Z","timestamp":1760148112513,"version":"build-2065373602"},"reference-count":26,"publisher":"MDPI AG","issue":"4","license":[{"start":{"date-parts":[[2023,3,31]],"date-time":"2023-03-31T00:00:00Z","timestamp":1680220800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"name":"NSF-CMMI Award","award":["#2043325"],"award-info":[{"award-number":["#2043325"]}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["Computation"],"abstract":"<jats:p>Although the full form of the Rayleigh\u2013Plesset (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.<\/jats:p>","DOI":"10.3390\/computation11040072","type":"journal-article","created":{"date-parts":[[2023,3,31]],"date-time":"2023-03-31T10:19:33Z","timestamp":1680257973000},"page":"72","update-policy":"https:\/\/doi.org\/10.3390\/mdpi_crossmark_policy","source":"Crossref","is-referenced-by-count":0,"title":["A Robust Bubble Growth Solution Scheme for Implementation in CFD Analysis of Multiphase Flows"],"prefix":"10.3390","volume":"11","author":[{"given":"Hao","family":"Pang","sequence":"first","affiliation":[{"name":"Department of Mechanical and Aerospace Engineering, North Carolina State University, Raleigh, NC 27695, USA"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-6670-229X","authenticated-orcid":false,"given":"Gracious","family":"Ngaile","sequence":"additional","affiliation":[{"name":"Department of Mechanical and Aerospace Engineering, North Carolina State University, Raleigh, NC 27695, USA"}]}],"member":"1968","published-online":{"date-parts":[[2023,3,31]]},"reference":[{"key":"ref_1","doi-asserted-by":"crossref","first-page":"23020","DOI":"10.1016\/j.ijhydene.2020.05.224","article-title":"Temperature effect in cavitation risk assessments of polymers for hydrogen systems","volume":"45","author":"Melnichuk","year":"2020","journal-title":"Int. J. Hydrog. Energy"},{"key":"ref_2","doi-asserted-by":"crossref","first-page":"124426","DOI":"10.1016\/j.energy.2022.124426","article-title":"Dynamic mode decomposition to classify cavitating flow regimes induced by thermodynamic effects","volume":"254","author":"Ge","year":"2022","journal-title":"Energy"},{"key":"ref_3","doi-asserted-by":"crossref","first-page":"105988","DOI":"10.1016\/j.ultsonch.2022.105988","article-title":"A correlation between cavitation bubble temperature, sonoluminescence and interfacial chemistry\u2014A minireview","volume":"85","author":"Anandan","year":"2022","journal-title":"Ultrason. Sonochemistry"},{"key":"ref_4","doi-asserted-by":"crossref","first-page":"120970","DOI":"10.1016\/j.ijheatmasstransfer.2021.120970","article-title":"Cavitation dynamics and thermodynamic effects at elevated temperatures in a small Venturi channel","volume":"170","author":"Ge","year":"2021","journal-title":"Int. J. Heat Mass Transf."},{"key":"ref_5","doi-asserted-by":"crossref","first-page":"278","DOI":"10.1016\/j.renene.2021.04.038","article-title":"Experimental investigation of temperature effect on hydrodynamic characteristics of natural cavitation in rotational supercavitating evaporator for desalination","volume":"174","author":"Zheng","year":"2021","journal-title":"Renew. Energy"},{"key":"ref_6","doi-asserted-by":"crossref","first-page":"130470","DOI":"10.1016\/j.jclepro.2022.130470","article-title":"Intensity and regimes changing of hydrodynamic cavitation considering temperature effects","volume":"338","author":"Ge","year":"2022","journal-title":"J. Clean. Prod."},{"key":"ref_7","doi-asserted-by":"crossref","first-page":"106035","DOI":"10.1016\/j.ultsonch.2022.106035","article-title":"Combined suppression effects on hydrodynamic cavitation performance in Venturi-type reactor for process intensification","volume":"86","author":"Ge","year":"2022","journal-title":"Ultrason. Sonochemistry"},{"key":"ref_8","doi-asserted-by":"crossref","first-page":"933","DOI":"10.1016\/j.camwa.2017.10.021","article-title":"The space\u2013time CESE scheme for shallow water equations incorporating variable bottom topography and horizontal temperature gradients","volume":"75","author":"Saleem","year":"2018","journal-title":"Comput. Math. Appl."},{"key":"ref_9","doi-asserted-by":"crossref","first-page":"379","DOI":"10.1016\/j.rinp.2017.12.044","article-title":"Application of discontinuous Galerkin method for solving a compressible five-equation two-phase flow model","volume":"8","author":"Saleem","year":"2018","journal-title":"Results Phys."},{"key":"ref_10","doi-asserted-by":"crossref","first-page":"617","DOI":"10.1115\/1.1486223","article-title":"Mathematical basis and validation of the full cavitation model","volume":"124","author":"Singhal","year":"2002","journal-title":"J. Fluids Eng."},{"key":"ref_11","unstructured":"Zwart, P.J., Gerber, A.G., and Belamri, T. (June, January 30). A two-phase flow model for predicting cavitation dynamics. Proceedings of the Fifth International Conference on Multiphase Flow, Yokohama, Japan."},{"key":"ref_12","first-page":"561","article-title":"Development of a new cavitation model based on bubble dynamics","volume":"81","author":"Sauer","year":"2001","journal-title":"J. Appl. Math. Mech."},{"key":"ref_13","doi-asserted-by":"crossref","first-page":"339","DOI":"10.1016\/j.euromechflu.2003.10.005","article-title":"Numerical simulation of 3-D incompressible, multi-phase flows over cavitating projectiles","volume":"23","author":"Owis","year":"2004","journal-title":"Eur. J. Mech. B\/Fluids"},{"key":"ref_14","doi-asserted-by":"crossref","first-page":"15","DOI":"10.1155\/S1023621X04000028","article-title":"Numerical and experimental investigations of the cavitating behavior of an inducer","volume":"10","author":"Bakir","year":"2004","journal-title":"Int. J. Rotating Mach."},{"key":"ref_15","doi-asserted-by":"crossref","first-page":"02018","DOI":"10.1051\/epjconf\/20146702018","article-title":"Application of a pressure based CFD code with mass transfer model based on the Rayleigh equation for the numerical simulation of the cavitating flow around a hydrofoil with circular leading edge","volume":"Volume 67","author":"Deimel","year":"2014","journal-title":"EPJ Web of Conferences"},{"key":"ref_16","first-page":"4712","article-title":"Bubble dynamics model for predicting the growth and collapse of cavitation bubbles in diesel injector","volume":"24","author":"Bicer","year":"2014","journal-title":"At. Sprays"},{"key":"ref_17","doi-asserted-by":"crossref","first-page":"116202","DOI":"10.1143\/JJAP.49.116202","article-title":"Modification of Rayleigh\u2013Plesset theory for reproducing dynamics of cavitation bubbles in liquid-phase laser ablation","volume":"49","author":"Soliman","year":"2010","journal-title":"Jpn. J. Appl. Phys."},{"key":"ref_18","doi-asserted-by":"crossref","first-page":"155","DOI":"10.1016\/j.ijmultiphaseflow.2016.03.022","article-title":"Modeling of hydrodynamic cavitating flows considering the bubble-bubble interaction","volume":"84","author":"Ye","year":"2016","journal-title":"Int. J. Multiph. Flow"},{"key":"ref_19","doi-asserted-by":"crossref","first-page":"92","DOI":"10.1016\/j.icheatmasstransfer.2017.04.021","article-title":"Application of the semi-analytical cavitation model to flows in a centrifugal pump","volume":"86","author":"Ye","year":"2017","journal-title":"Int. Commun. Heat Mass Transf."},{"key":"ref_20","doi-asserted-by":"crossref","first-page":"482","DOI":"10.1016\/j.ultsonch.2014.05.011","article-title":"Modeling cavitation in a rapidly changing pressure field\u2013application to a small ultrasonic horn","volume":"22","author":"Mettin","year":"2015","journal-title":"Ultrason. Sonochem."},{"key":"ref_21","doi-asserted-by":"crossref","first-page":"012037","DOI":"10.1088\/1757-899X\/492\/1\/012037","article-title":"Numerical prediction of the gas content effect on the cavitation characteristics of the pump using the simplified Rayleigh-Plesset equation","volume":"Volume 492","author":"Lomakin","year":"2019","journal-title":"IOP Conference Series: Materials Science and Engineering"},{"key":"ref_22","doi-asserted-by":"crossref","first-page":"1865","DOI":"10.1002\/fld.2452","article-title":"A numerical scheme for Euler\u2013Lagrange simulation of bubbly flows in complex systems","volume":"67","author":"Shams","year":"2011","journal-title":"Int. J. Numer. Methods Fluids"},{"key":"ref_23","doi-asserted-by":"crossref","first-page":"780","DOI":"10.1002\/nme.2032","article-title":"Numerical analysis of Rayleigh\u2013Plesset equation for cavitating water jets","volume":"72","author":"Alehossein","year":"2007","journal-title":"Int. J. Numer. Methods Eng."},{"key":"ref_24","doi-asserted-by":"crossref","first-page":"022009","DOI":"10.1063\/1.4942237","article-title":"Evolution of spherical cavitation bubbles: Parametric and closed-form solutions","volume":"28","author":"Mancas","year":"2016","journal-title":"Phys. Fluids"},{"key":"ref_25","doi-asserted-by":"crossref","unstructured":"Franc, J.P., and Michel, J.M. (2004). Fundamentals of Cavitation, Springer.","DOI":"10.1007\/1-4020-2233-6"},{"key":"ref_26","doi-asserted-by":"crossref","unstructured":"Hilgenfeldt, S., Brenner, M.P., Grossmann, S., and Lohse, D. (1998). Analysis of Rayleigh\u2013Plesset Dynamics for Sonoluminescing Bubbles, Cambridge University Press.","DOI":"10.1017\/S0022112098001207"}],"container-title":["Computation"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.mdpi.com\/2079-3197\/11\/4\/72\/pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,10,10]],"date-time":"2025-10-10T19:07:49Z","timestamp":1760123269000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.mdpi.com\/2079-3197\/11\/4\/72"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2023,3,31]]},"references-count":26,"journal-issue":{"issue":"4","published-online":{"date-parts":[[2023,4]]}},"alternative-id":["computation11040072"],"URL":"https:\/\/doi.org\/10.3390\/computation11040072","relation":{},"ISSN":["2079-3197"],"issn-type":[{"type":"electronic","value":"2079-3197"}],"subject":[],"published":{"date-parts":[[2023,3,31]]}}}