3D Finite different method using matlab

3 replies · 3K views · last reply
nnadiya
Messages
2
Reaction score
0
Can anyone show me how to solve the 3D diffusion equation which has been modeled into FDM by using matlab?
 
Suppose your stencil only uses nearest and next nearest neighbors on a cubic periodic domain with N lattice points in each direction. Define I am = (1:N)-1;Im(Im==0)=N; and Ip = 2:N+1;Ip(Ip==N+1)=1; Define your system and matrices and use A(lp,lp,lm) ## \sim A_{i+1,i+1,i-1} ## in your stencil.
 
Strum said:
Suppose your stencil only uses nearest and next nearest neighbors on a cubic periodic domain with N lattice points in each direction. Define I am = (1:N)-1;Im(Im==0)=N; and Ip = 2:N+1;Ip(Ip==N+1)=1; Define your system and matrices and use A(lp,lp,lm) ## \sim A_{i+1,i+1,i-1} ## in your stencil.
Can I learn further from you?
 
You are sure you don't want to try yourself first? Also you haven't mentioned which time update algorithm you are using.