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.

本Matlab代码包含了氧化铬（chromia）在Na₂O-SiO₂熔体中溶解过程的方程求解器。
上述方程采用有限体积法进行离散求解，并通过搭载自适应时间步长的显式欧拉法完成时间积分。代码实现了简单的一维计算域，空间网格与求解器的计算精度均可在Input.m文件中调整。
若要执行计算，请在Inputs.m文件中修改计算条件，确保将Execute.m文件中指定的全部所需文件放置于同一文件夹下，随后运行Execute.m文件即可。
模拟结果存储于"res"矩阵（res(x,u,t)）中：
res(:,1,:) 为液/氧化物界面的距离（单位：米）
res(:,2,:) 为三价铬的原子分数
res(:,3,:) 为六价铬的原子分数
res(:,4,:) 为氧化铬厚度（单位：米）
res(:,5,:) 为氧气分压（单位：标准大气压）
res(:,6,:) 为Si Q4的原子分数
res(:,7,:) 为Si Q3的原子分数
res(:,8,:) 为Si Q2的原子分数
res(:,9,:) 为钠的原子分数
res(:,10,:) 为熔体中的速度场（单位：米/秒）
res(:,11,:) 为熔体中的电场（单位：伏/米）
最终计算结果将保存至results_tf.txt文件中。