Implementing the He/Lee model was a real crash-test for Palabos. The model requires couplings between many scalar-fields, vector-fields, and lattices, and uses next-to-nearest-neighbor interaction stencils. Test result positive: the model works fine, runs fast, and scales impressively. This adds a new multi-phase model to Palabos which handles large densit and viscosity ratios between the two phases.

Use the data processors SetExternalScalarFunctionalXD or SetExternalVectorFunctionalXD, respectively their wrapper functions setExternalScalar and setExternalVector, to initialize external scalar- or vector-fields of a lattice. Example: an external force term may be initialized to a constant value through a command like