Code covered by the BSD License  

Highlights from
immiscible LB

5.0

5.0 | 2 ratings Rate this file 38 Downloads (last 30 days) File Size: 356 KB File ID: #24824

immiscible LB

by

 

23 Jul 2009 (Updated )

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

| Watch this File

File Information
Description

The MatLab code uses lattice Boltzmann (LB) method to simulate two phase flows for immiscible fluids (blue & red fluids) in 2D according to the D2Q9 scheme.
The model is also referred to as color model or chromodynamic model or ILB (immiscible LB).
This MatLab implementation of ILB saves an AVI file in '' dir_avi='C:\' '' that records frames generated by the time evolution iterations.

The code implements in D2Q9 the model originally developed by Gunstensen [Gunstensen AK, Rothman D. Lattice Boltzmann model of immiscible fluids. Phys Rev A 1991;43(8):4320?4327], which is based on the original Rothman?Keller model [Rothman DH, Keller JM. Immiscible cellular automaton fluids. J Stat Phys 1988;52:1119?1127].
see also : Gunstensen, A., Rothman, D.H., Zaleski, S., and Zanetti, G., Physical Review A 43, 4320-4327, April 15, 1991.
and
page 190-191 of the reference book by S. Succi, The Lattice Boltzmann Equation, Oxford Univ. Press, 2001.
Only periodic boundary conditions are implemented in this version.
Also the method suggested by Dr Sergey Lishchuk for including surface tension is implemented, see :
S. V. Lishchuk, C. M. Care, I. Halliday. Lattice Boltzmann algorithm for surface tension with greatly reduced micro-currents. Phys. Rev. E, 2003,
In order to run this code you need MatLab 2009 or later with the image processing toolbox (IPT).
The files Run* are those to execute (eg. Run_Demo_ILBGK_D2Q9_Spinodal_Decomposition OR Run_ILB_Flow_inPM_Periodic_sigle_bubble).
If you stop the execution you can restart with >Restart_Time_Evolution

Required Products Image Processing Toolbox
MATLAB release MATLAB 7.9 (R2009b)
Other requirements IMAGE PROCESSING TOOLBOX tested on Windows XP 64
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (2)
03 Jul 2014 Isaac  
14 Aug 2012 Victoria  
Updates
28 Jul 2010

faster than previous submission

22 Sep 2010

new examples in Run* files

28 Sep 2010

corrects minor bug surface tension treatment

Contact us