def main():
model = Problem().model(
model="2d",
dx=1.0e-1,
dt=1.0e-3,
t0=0.0,
u0=0.1,
order=5,
num_nodes=400,
num_iters=5000,
pumping=GaussianRingPumping2D(power=20.0, radius=10.0, variation=3.14),
original_params={
"R": 0.0242057488654,
"gamma": 0.0242057488654,
"g": 0.00162178517398,
"tilde_g": 0.0169440242057,
"gamma_R": 0.242057488654,
},
dimless_params={},
)
# Obtain steady state solution
solution = model.solve()
solution.report()
solution.visualize()
solution.show()
def solve_problem(desc):
desc = loads(desc.encode('utf8'))
model = Problem().model(model='2d',
dx=2.0e-1,
dt=2.0e-3,
t0=0.0,
u0=0.1,
order=5,
num_nodes=200,
num_iters=2000,
pumping=GaussianPumping2D(power=desc['power'],
variation=5.0),
original_params={
'R': 0.0242057488654,
'gamma': 0.0242057488654,
'g': 0.00162178517398,
'tilde_g': 0.0169440242057,
'gamma_R': 0.242057488654,
},
dimless_params={})
solution = model.solve()
solution.report()
solution.store()
return True
def solve_problem(desc):
desc = loads(desc.encode('utf8'))
model = Problem().model(
model = '2d',
dx = 2.0e-1,
dt = 2.0e-3,
t0 = 0.0,
u0 = 0.1,
order = 5,
num_nodes = 200,
num_iters = 2000,
pumping = GaussianPumping2D(power=desc['power'], variation=5.0),
original_params = {
'R': 0.0242057488654,
'gamma': 0.0242057488654,
'g': 0.00162178517398,
'tilde_g': 0.0169440242057,
'gamma_R': 0.242057488654,
},
dimless_params = {
})
solution = model.solve()
solution.report()
solution.store()
return True
def main():
model = Problem().model(model='1d',
dx=1.0e-1,
dt=1.0e-3,
t0=0.0,
u0=0.1,
order=5,
num_nodes=400,
num_iters=10000,
pumping=GaussianRingPumping1D(power=20.0,
radius=10.0,
variation=3.14),
original_params={
'R': 0.0242057488654,
'gamma': 0.0242057488654,
'g': 0.00162178517398,
'tilde_g': 0.0169440242057,
'gamma_R': 0.242057488654,
},
dimless_params={})
# Obtain steady state solution
solution = model.solve()
solution.report()
solution.visualize()
solution.show()