Hi
Could you please help me with a changing boundary condition during simulation?
It has been mentioned that the "update_node" do this but it seems not working.
class LDCSim(LBFluidSim):
subdomain = LDCBlock # DuctSubdomain
@classmethod
def update_defaults(cls, defaults):
defaults.update({
'lat_nx': 128,
'lat_ny': 64,
'lat_nz': 64,
'grid': 'D3Q19'})
prev_f = None
every = 500
def after_step(self, runner):
if self.iteration % self.every == 0:
wallBC = NTFullBBWall
hx, hy, hz = runner._subdomain._get_mgrid()
whr = (np.square(hx - 80) + np.square(hy - 16) + np.square(hz - 32)) <= np.square(8)
runner._subdomain.update_node(whr,wallBC)