Data Documentation for the Original Research Paper (ORP) entitled: "An improved, Rayleigh-Plesset based homogeneous cavitation model accounting for microbubble behaviour and turbulent interaction".


General Information
Simulation results obtained from CFD (Computer Fluid Dynamics) calculations associated with the aforementioned ORP authored by Álvaro Pardo Vigil, Laura Suárez Fernández, José González Pérez and Adrián Pandal Blanco, all members of the FERES Research Group in the University of Oviedo. These datasets are protected under Creative Commons License: CC BY-NC-SA.

Name of data set
"Raw Data of improved RP model"

Name of data files in data set
Variable information at each time step is denoted XXXX-YYYY, where:
	XXXX: Name of the working point, characterized by its cavitation number "sigma". When the nomenclature is 	"Sigma055", this stands for sigma = 0.55.
	YYYY: Simulation time step at which the data are collected. This is easily translated to time units 	multiplying by the fixed time step of 0.5 microseconds (5e-7 [s])
Information in these files is given as a function of the coordinates of the mesh point where the results are evaluated.
Data included in ".out" files for grid convergence represents vertical coordinates in [m] versus the variable of interest in each case

Dataset language
English (although some sentences in Spanish may appear)

Date the data set was last modified
14/01/2025

Funder
At data collection date, Álvaro Pardo Vigil was supported by the Spanish Ministry of Science, Innovation and Universities within the “FPU” Program (grant reference FPU23/03450).


Methodology for data collection
The data set is composed by numerical data obtained from CFD simulations carried out in the commercial software ANSYS Fluent 2023 R2.

Data collector
Álvaro Pardo Vigil

Date of data collection
06/11/2024

Person to contact with questions
Adrián Pandal Blanco (pandaladrian@uniovi.es)

Data entry
06/11/2024

Software used to prepare data set
ANSYS Fluent, Maltab, MS Excel and Windows Notepad

Data processing that was performed
Grid Convergence: Computational meshes of increasing cell count were tested and the velocity, turbulent kinetic energy and vapour volume fraction profiles in the injector throat were exported (in .out format) to Matlab in order to assess grid convergence. Based on these results, the Grid Convergence Index (GCI) was found.
Simulation Results: ASCII-type ".out" files were exported every 50 simulation steps. Not only vapour volume fraction, but also pressure, velocity and turbulence variables were exported for further postprocessing. In particular, for the slip velocity simulations the X and Y components of liquid and bubbly phases were treated in Matlab to obtain the non-dimensional relative velocity fields of interest. Each of these files is treated as a video frame, and exported to .png format in the separate "Contours" folder.


Variables
. x: X coordinate in the mesh [m]
. x: Y coordinate in the mesh [m]
. p: Static pressure [Pa]
. pabs: Absolute pressure [Pa]
. pT: Total pressure [Pa]
. rho: Mixture density [kg/m3]
. vx: Velocity magnitude in X [m/s]
. vy: Velocity magnitude in Y [m/s]
. k: Turbulent kinetic energy in [m2/s2]
. e: Turbulent dissipation rate in [m2/s3]
. rhob: Bubbly phase (vapour) density [kg/m3]
. vof: Volume fraction of the bubbly phase [-]

File Overview
Inside the main folder "Raw Dataof improved RP model", three folders are available. In addition, a ".xlsx" spreadsheet including the boundary conditions for each of the working points analyzed is provided.
In the Grid Convergence folder, the following files are available:
	The ".out" files exported from ANSYS Fluent indicating the profile of different variables in the throat 	section.
	The "GridConvergence.m" Matlab file plots these profiles in order to facilitate visualization

In the Simulation Results folder, the results for the different working points under study are denoted:
	"Sigma055" and "Sigma065", respectively corresponding to sigma = 0.55 (hydraulic flip) and sigma = 0.65 	cavitation numbers, contain the text files with the simulation results and the matlab code that allows to 	obtain the images available in the "Contours" folder.
	"Sigma078" corresponds to a sigma = 0.78 cavitation number, and is subdivided in the following folders:
		"Standard R0" for the simulation results obtained using the proposed model.
		"Increased R0" for the simulation results obtained using the proposed model and a greater initial 		bubble radius.
		"Slip" for the simulation results obtained using the proposed model and a Manninen et al. 			slip formulation.
		"FullSourceTerm" for the simulation results obtained using the proposed model taking into account 		gas compressibility in the source term.

Finally, the UDFs folder contains the User-Defined Functions that were used to implement the improved cavitation model in ANSYS Fluent:
	"FullSourceTerm" is the full source term for the vapour transport equation including gas compressibility 	effects.
	"fullRP_DEFINITIVA" is the source term for the vapour transport equation neglecting gas compressibility.
	"TurbulentViscosity" is the modified turbulent viscosity formula for the k-e model.
	"VapourEqDensityNuevo" is the modified vapour density equation that includes both evaporated water and 	noncondensable gases in a single phase.
	"BubbleDiameterImproved" is the UDF that calculates the volumetric average bubble diameter from the void 	fraction.