no slip BC
0001 function bc=rbconvbc_noslip(p,u,lam) 0002 % no slip BC 0003 0004 stifffac=10^6; % stiff spring factor 0005 qd = stifffac*[[1;0;0] [0;1;0] [0;0;1]]; % Dirichlet 0006 qn = zeros(p.neq); % zero flux: 0007 gv = zeros(p.neq,1); % values 0008 % mixed: zero flux for psi and T, and psi=0 0009 qm = qn; qm(2,1)=stifffac; 0010 % ordering: bottom, right, top, left 0011 bc = gnbc(p.neq,qd,gv,qm,gv,qd,gv,qm,gv);