next up previous contents
Next: Summary Up: Gravity in a Lattice Previous: Immiscible Binary Fluid

Galilean Invariance

  To test the Galilean invariance of the model a drop of radius tex2html_wrap_inline15559 lattice units was initialised at the centre of a 128 by 128 grid with the whole fluid moving with speed tex2html_wrap_inline12169 in the x-direction and with gravity acting in the z-direction with strength g = 0.0005. The bubble was then allowed to equilibrate for a number of different values of the relaxation time associated with the order parameter, tex2html_wrap_inline15569 . The fluid relaxation parameter tex2html_wrap_inline15571 was set to 1.1 throughout. The ratio of tex2html_wrap_inline15573 , the drop diameter in the z-direction, to tex2html_wrap_inline15575 , the diameter in the x-direction, is shown in figure 5-14.

  figure5080

When tex2html_wrap_inline14895 the radio is independent of tex2html_wrap_inline12169 as observed elsewhere [38] in the absence of gravity. For other values of tex2html_wrap_inline15569 the ratio is dependent on tex2html_wrap_inline12169 and the model is not Galilean invariant, the further tex2html_wrap_inline15569 is from tex2html_wrap_inline14897 the worse the lack of Galilean invariance is. The ratio of the two diameters appears independent of the velocity when tex2html_wrap_inline14895 and is slightly smaller than unity. The difference from unity is due to the density gradient across the drop. This does not affect the Galilean invariance because the free energy binary-fluid model is Galilean invariant [38] for all densities and the transport coefficients are independent of density. The constant value of tex2html_wrap_inline15611 shows that this model, incorporating the gravitational interaction, is, at worst, very close to Galilean invariant provided tex2html_wrap_inline14895 .



James Buick
Tue Mar 17 17:29:36 GMT 1998