10 likes | 104 Views
% Face 1--values of velocity vxf=(5*vxL(1)+5*vxL(2)+2*vxL(3))/12; vyf=(5*vyL(1)+5*vyL(2)+2*vyL(3))/12;. qf=vxf*dely-vyf*delx; %Face 1 volume flux. Face 1--Advection Coefficients ap(i)= ap(i) +max(qf,0); a(i,j)=a(i,j)+max(-qf,0);.
E N D
% Face 1--values of velocity vxf=(5*vxL(1)+5*vxL(2)+2*vxL(3))/12; vyf=(5*vyL(1)+5*vyL(2)+2*vyL(3))/12; qf=vxf*dely-vyf*delx; %Face 1 volume flux Face 1--Advection Coefficients ap(i)= ap(i) +max(qf,0); a(i,j)=a(i,j)+max(-qf,0);