Data for: A new modeling of the dissolution of chromia in Na2O-SiO2 liquids

This matlab code contains the solver for the equations of the dissolution of chromia in Na2O-SiO2 liquids.
The equations are solved using the finite volume method and integrated over time using an Euler explicit method with an adaptative time step. A simple 1D domaine is implemented. The spatial mesh and the accuracy of the solver can be modified in the Input.m file.
To performed a calculation, modify the conditions into the Inputs.m file, be sure to get all of the required files specified in the Execute.m file in the same folder and execute the file Execute.m file.
The results of the simulation are given in the "res" matrix (res(x,u,t)).
res(:,1,:) = distance from the liquid/oxide interface (m)
res(:,2,:) = atomic fraction of CrIII
res(:,3,:) = atomic fraction of CrVI
res(:,4,:) = chromia thickness (m)
res(:,5,:) = partial pressure of O2 (atm)
res(:,6,:) = atomic fraction of Si Q4
res(:,7,:) = atomic fraction of Si Q3
res(:,8,:) = atomic fraction of Si Q2
res(:,9,:) = atomic fraction of Na
res(:,10,:) = velocity field in the liquid (m/s)
res(:,11,:) = electric field in the liquid (V/m)
The final result is save into the results_tf.txt file.