颗粒融合的相场模拟

matlab环境下的简单香相场模型，简单入门级的代码，其中包括5个文件，主程序，初始化条件，求解方程两个文件，输出结果

%== get intial wall time:
time0=clock（）;
format long;

%-- Simulation cell parameters:

Nx = 1000;
Ny = 1000;
NxNy= Nx*Ny;

dx = 0.5;
dy = 0.5;

%--- Time integration parameters:

nstep = 20000;
nprint = 100;
dtime= 1.0e-4;

%--- Material specific Parameters:
npart = 2;
%
coefm = 5.0;
coefk = 2.0;
coefl = 10.0;
%
dvol = 0.040;
dvap = 0.002;
dsur = 16.0;
dgrb = 1.6;

%----
%prepare microstructure:
%----

iflag = 2;

[npartetascon] = micro_sint_pre（NxNynpartiflag）;

%--- Get Laplacian templet:

[grad] =laplacian（NxNydxdy）;

%---
%- Evolution
%---

eta =zeros（NxNy1）;

for istep =1:nstep

%-- evolve concentration:

iflag =1;

%-- derivative of free energy:

[dfdcondfdeta]=free_energ_sint_v2（NxNyconetaetasnpartiflag）;

%

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2644  2019-03-07 10:43  fd.m

     文件        805  2019-02-10 17:33  free_energ_sint_v2.m

     文件        394  2019-02-10 17:47  laplacian.m

     文件       1639  2019-02-10 17:32  micro_sint_pre.m

     文件       1002  2019-02-27 15:45  write_vtk_grid_values.m

----------- ---------  ---------- -----  ----

                 6484                    5