Code covered by the BSD License  

Highlights from
immiscible LB

immiscible LB

by

 

23 Jul 2009 (Updated )

Implements Immiscible Lattice Boltzmann (ILB, D2Q9) method for two phase flows

Red_and_Blue_Quiver.m
       %set(gca,'nextplot','replacechildren');
       iR_domina=(rhoR(ija)>rhoB(ija)); 
       idom=find(iR_domina); 
       idom=idom(find(jump(1:length(idom))));
       %[jawet(idom),Nr-iawet(idom)]
       %[ux(idom)./vmax,-uy(idom)./vmax]
       quiver(double(jawet(idom)),double(Nr-iawet(idom)),ux(idom)./vmax,-uy(idom)./vmax,expf,'r');
       idom=find(~iR_domina); idom=idom(find(jump(1:length(idom))));
       hold on
       quiver(double(jawet(idom)),double(Nr-iawet(idom)),ux(idom)./vmax,-uy(idom)./vmax,expf,'b');
       %axis off , title('Velocity arrows')

Contact us