def boundary_conditions(self, hx, hy):
where = (hx == hy)
self.set_node(
where,
NTEquilibriumVelocity(
multifield((0.01 * (hx - self.gy / 2)**2, 0.0), where)))
self.set_node((hx > 10) & (hy < 5), NTFullBBWall)
def boundary_conditions(self, hx, hy):
where = (hx == hy)
self.set_node(
where,
NTEquilibriumVelocity(
multifield((0.01 * (hx - self.gy / 2)**2, 0.0), where)))
where = ((hx == 5) & (hy == 7))
self.set_node(where,
NTEquilibriumDensity(DynamicValue(0.1 * S.gx)))
# Interpolated time series.
data = np.linspace(0, 50, 10)
where = ((hx == 5) & (hy == 8))
self.set_node(
where,
NTEquilibriumDensity(
DynamicValue(0.1 * S.gx *
LinearlyInterpolatedTimeSeries(data, 40))))
# Same underlying data, but different time step.
where = ((hx == 5) & (hy == 9))
self.set_node(
where,
NTEquilibriumDensity(
DynamicValue(0.1 * S.gx *
LinearlyInterpolatedTimeSeries(data, 30))))
self.set_node((hx > 10) & (hy < 5), NTFullBBWall)
def boundary_conditions(self, hx, hy, hz):
where = np.logical_and((hx == hy), (hy == hz))
self.set_node(
where,
NTEquilibriumVelocity(
multifield((0.01 * (hy - self.gy / 2)**2, 0.03 *
(hz - self.gz / 2)**2, 0.0), where)))
self.set_node((hx > 10) & (hy < 5) & (hz < 7), NTFullBBWall)
def boundary_conditions(self, hx, hy, hz):
where = np.logical_and((hx == hy), (hy == hz))
self.set_node(
where,
NTEquilibriumVelocity(
multifield((0.01 * (hy - self.gy / 2) ** 2, 0.03 * (hz - self.gz / 2) ** 2, 0.0), where)
),
)
self.set_node((hx > 10) & (hy < 5) & (hz < 7), NTFullBBWall)
def boundary_conditions(self, hx, hy):
where = (hx == hy)
self.set_node(where, NTEquilibriumVelocity(
multifield((0.01 * (hx - self.gy / 2)**2, 0.0), where)))
where = ((hx == 5) & (hy == 7))
self.set_node(where, NTEquilibriumDensity(
DynamicValue(0.1 * S.gx)))
# Interpolated time series.
data = np.linspace(0, 50, 10)
where = ((hx == 5) & (hy == 8))
self.set_node(where, NTEquilibriumDensity(
DynamicValue(0.1 * S.gx * LinearlyInterpolatedTimeSeries(data, 40))))
# Same underlying data, but different time step.
where = ((hx == 5) & (hy == 9))
self.set_node(where, NTEquilibriumDensity(
DynamicValue(0.1 * S.gx * LinearlyInterpolatedTimeSeries(data, 30))))
self.set_node((hx > 10) & (hy < 5), NTFullBBWall)
def boundary_conditions(self, hx, hy):
where = (hx == hy)
self.set_node(where,
NTEquilibriumVelocity(
multifield((0.01 * (hx - self.gy / 2)**2, 0.0), where)))
def boundary_conditions(self, hx, hy):
where = hx == hy
self.set_node(where, NTEquilibriumVelocity(multifield((0.01 * (hx - self.gy / 2) ** 2, 0.0), where)))
self.set_node((hx > 10) & (hy < 5), NTFullBBWall)
