Numerical simulations of the latest caldera-forming eruption of Okmok volcano, Alaska

Raw outputs of the numerical simulations of an explosive volcanic eruption.

The file "source.zip" includes all the source files modified compared to the defaultMFIX 2016-1 files (www.mfix.netl.doe.gov, v.2016). It also includes an example of the input file mfix.dat (run 3). Theseare raw file intended for enabling result reproductibility. They contain unused optional variables and are not intented (and commented) to serve as tutorials. Thesubfolder "postmfix" includes the source file for option 6 of the post processingprogram. Some options are user-defined, and others are hard coded.

The file "AllRuns.zip" inculde all the Okmok runs. All runs were generated with MFIX 2016-1 in TFM mode. They are organized in folders corresponding to the namingconvention of Table 1. Some runs were generated in several sequences (e.g., 0-300 s,then 300-700 s) that were stitched into a single output bundle. Other sequences werekept separated.

The files OK.RES can be openend with Paraview >5.6 and selecting the MFIXReader.

The variables are:

EP_g = Gas volume fractionGas Velocity = Gas velocity vector (m/s)P_g = Pressure (Pa)P_star = Solid pressure (Pa)ROP_s_m = Solid density times particle volume fraction of solid phase m (kg/m3)RRates_1 = Gas density (kg/m3)RRates_2 = Gas viscosity (Pa s)Solids_Velocity_x = Velocity vector of solid phase x (m/s)T_g = Gas temperature (K)T_s_x = Temperature of solid phase x (K)Theta_m_x = Granular temperature of solid phase x (m2/s2)U_g = Horizontal component of gas velocity (m/s)V_g = Vertical component of gas velocity (m/s)W_g = Z component of gas velocity (0 because axisymmetric domain)U_s_x = Horizontal component of velocity of solid phase x (m/s)V_s_x = Vertical component of velocity of solid phase x (m/s)W_s_x = Z component of velocity of solid phase x (0 because axisymmetric domain)X_g_x, X_s_x = unusedk_turb_g (not always recorded) = Turbulent kinetic energy (J/kg)e_turb_g = duplicate of k_turb_g

Refer to the MFIX documentation for more details on the format and on the model itself.