Data for: A new modeling 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.