def run(dx,
Tf,
generator=pyLBM.generator.CythonGenerator,
sorder=None,
withPlot=True):
"""
Parameters
----------
dx: double
spatial step
Tf: double
final time
generator: pyLBM generator
sorder: list
storage order
withPlot: boolean
if True plot the solution otherwise just compute the solution
"""
la = 1.
rho0 = 1.
Re = 2000
nu = 5. / Re
s1 = 1.6
s2 = 1.2
s4 = 1.6
s9 = 1. / (3 * nu + .5)
s11 = s9
s14 = 1.2
r = X**2 + Y**2 + Z**2
dico = {
'box': {
'x': [0., 2.],
'y': [0., 1.],
'z': [0., 1.],
'label': [1, 2, -1, -1, -1, -1]
},
'elements': [pyLBM.Sphere((.3, .5, .5), 0.125, 0)],
'space_step':
dx,
'scheme_velocity':
la,
'schemes': [{
'velocities':
list(range(7)) + list(range(19, 27)),
'conserved_moments': [mass, qx, qy, qz],
'polynomials': [
1, r - 2, .5 * (15 * r**2 - 55 * r + 32), X,
.5 * (5 * r - 13) * X, Y, .5 * (5 * r - 13) * Y, Z,
.5 * (5 * r - 13) * Z, 3 * X**2 - r, Y**2 - Z**2, X * Y, Y * Z,
Z * X, X * Y * Z
],
'relaxation_parameters':
[0, s1, s2, 0, s4, 0, s4, 0, s4, s9, s9, s11, s11, s11, s14],
'equilibrium': [
mass, -mass + qx**2 + qy**2 + qz**2, -mass, qx, -7. / 3 * qx,
qy, -7. / 3 * qy, qz, -7. / 3 * qz,
1. / 3 * (2 * qx**2 - (qy**2 + qz**2)), qy**2 - qz**2, qx * qy,
qy * qz, qz * qx, 0
],
'init': {
mass: 1.,
qx: 0.,
qy: 0.,
qz: 0.
},
}],
'boundary_conditions': {
0: {
'method': {
0: pyLBM.bc.Bouzidi_bounce_back
}
},
1: {
'method': {
0: pyLBM.bc.Bouzidi_bounce_back
},
'value': bc_up
},
2: {
'method': {
0: pyLBM.bc.Neumann_x
}
},
},
'parameters': {
LA: la
},
'generator':
generator,
}
sol = pyLBM.Simulation(dico, sorder=sorder)
x, y, z = sol.domain.x, sol.domain.y, sol.domain.z
im = 0
compt = 0
while sol.t < Tf:
sol.one_time_step()
compt += 1
if compt == 128 and withPlot:
if mpi.COMM_WORLD.Get_rank() == 0:
sol.time_info()
im += 1
save(x, y, z, sol.m, im)
compt = 0
return sol
from __future__ import print_function
from __future__ import division
# Authors:
# Loic Gouarin <*****@*****.**>
# Benjamin Graille <*****@*****.**>
#
# License: BSD 3 clause
"""
Example of a 3D geometry: the cube [0,1] x [0,1] x [0,1] with a spherical hole
"""
import pyLBM
dico = {
'box': {
'x': [0, 1],
'y': [0, 1],
'z': [0, 1],
'label': 0
},
'elements': [pyLBM.Sphere((.5, .5, .5), .25, label=1)],
}
geom = pyLBM.Geometry(dico)
print(geom)
geom.visualize(viewlabel=True)
# Authors:
# Loic Gouarin <*****@*****.**>
# Benjamin Graille <*****@*****.**>
#
# License: BSD 3 clause
"""
Example of the cube in 3D with a spherical hole
"""
from six.moves import range
import pyLBM
dico = {
'box': {
'x': [0, 2],
'y': [0, 2],
'z': [0, 2],
'label': list(range(1, 7))
},
'elements': [pyLBM.Sphere((1, 1, 1), 0.5, label=0)],
'space_step': 0.25,
'schemes': [{
'velocities': list(range(19))
}]
}
dom = pyLBM.Domain(dico)
print(dom)
dom.visualize(view_distance=[1, 3],
view_in=False,
view_out=False,
view_bound=True)