Our ILB (Immiscible Lattice Boltzmann) is the Rothman -Gunstensen (in D2Q9) .

We would like to introduce superficial tension values in our code.

We thought of converting the physical value of the superficial tension (sigma_p= 0.035 N/m) using converting factors

sigma_lb=sigma_p.*(Time_fact).^-2;

Time_fact=(U_lb./U_p).*(L_p/L_lb);

where U is fluid velocity and L is length, the subscripts p and lb are used to identify physical and lattice units respectively.

or (alternatively) to set a Capillary Number (Ca_N) and calculate the superficial tension in lattice units (sigma_lb) as

sigma_lb= (U_lb.*nuK_lb)./Ca_N;

where nuK_lb is our kinematic viscosity in lattice units and U_lb is the fluid velocity in lb units.

Rothman - Keller - Gunstensen etc suggest to use the relationship:

lambda=1./tau; density=1;

A=(sigma_p*lambda)./(192*density);

where , as far as we understand, the parameter A is the superficial tension in lattice units.

tau is the relaxation factor.

The two approaches lead to very different values and sigma_lb is orders of magnitude larger than A. Typical values of A are 1E-4. We calculate sigma_lb values of order 1.E2.

Can you help us to understand? why conversion factor does not work in this specific case?