資源簡介
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
評論
共有 條評論