Hi every one
I have a flow containing reacting chemical species can be modelled within LBM.
I have a good background in LBM for fluid flow but no idea about chemical reaction. Who can help me to accomplish this code?
To be specific, here is a simple problem:

A cubic simulation cell Lx=Ly=Lz=L

Periodic boundary conditions along x,y,z

A solid sphere of radius R=L/4 at the center (x=L/2, y=L/2, z=L/2)

A gas enters from the bottom (z=0) with constant velocity: v(x,y, z=0)=v_0 for all x,y

At the entrance (z=0) two spatially separated reacting species also enter the channel. This means:
C_A(x,y,0)=C_1 for x<L/2 for all times t
C_A(x,y,0)=0 for x>L/2 for all times t
AND
C_B(x,y,0)=0 for x<L/2 for all times t
C_B(x,y,0)=C_2 for x>L/2 for all times t
At time t=0 there is no chemical species within the simulation cell. A and B undergo chemical reaction with oneanother with a rate k_AB.
I have purposely left out some important ingredients which are necessary to make the problem statement complete.