def __init__(self, geo_class, options=[], args=None, defaults=None):
lbm.LBMSim.__init__(self, geo_class, options, args, defaults)
self._set_grid(self.options.grid)
# If the model has not been explicitly specified by the user, try to automatically
# select a working model.
if 'model' not in self.options.specified and (defaults is None or 'model' not in defaults.keys()):
self._set_model(self.options.model, 'mrt', 'bgk')
else:
self._set_model(self.options.model)
self.num_tracers = self.options.tracers
self.incompressible = self.options.incompressible
self.equilibrium, self.equilibrium_vars = sym.bgk_equilibrium(self.grid)
def __init__(self, geo_class, options=[], args=None, defaults=None):
super(ShanChenBinary, self).__init__(geo_class, options, args, defaults)
self.equilibrium, self.equilibrium_vars = sym.bgk_equilibrium(self.grid)
eq2, _ = sym.bgk_equilibrium(self.grid, self.S.phi, self.S.phi)
self.equilibrium.append(eq2[0])
self.add_force_coupling(0, 1, 'SCG')
def __init__(self, config):
super(LBFluidSim, self).__init__(config)
self.equilibrium, self.equilibrium_vars = sym.bgk_equilibrium(self.grid)
def __init__(self, config):
super(LBBinaryFluidShanChen, self).__init__(config)
self.equilibrium, self.equilibrium_vars = sym.bgk_equilibrium(self.grid)
eq2, _ = sym.bgk_equilibrium(self.grid, self.S.phi, self.S.phi)
self.equilibrium.append(eq2[0])
self.add_force_coupling(0, 1, 'SCG')