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

ux_uy_Lishchuk_correct.m
if(1) % surface tension as body force LISHCHUK - CARE - HALLIDAY 
    
% K_Lish
          
        relax_time=Sigma;
       % [~,ijrel]= intersect(ija,ijagr); %ijagr=ija(Flag_Grad); 
       %  K_Lish(~ijagr)=0; CGx(~ijagr)=0; CGy(~ijagr)=0;
       % maxK=max(K_Lish(ija));minK=min(K_Lish(ija));
       % K_Lish(ija)=max(K_Lish(ija),0);
       
        ux(:) = ux(:) + relax_time.*K_Lish(ija).*CGx(ija);
        uy(:) = uy(:) + relax_time.*K_Lish(ija).*CGy(ija);
        
      % imshow(K_Lish,[]);
    
    
end

Contact us