'omega_momentum': 0.7,
'omega_p_corr': 0.3,
'advection_scheme': 'upwind',
'maxiters': 200,
'n_p_corrs': 3,
}
grid = FvRectilinearGrid((60, 180), (1, 3))
if __name__ == '__main__':
from solvers.multiphase_piso import MultiphasePiso
import grid.viewers as viewers
import initial_conditions.circle
solver = MultiphasePiso(grid, **input)
f, = grid.get_cell_fields('f')
coords = grid.cell_nodes
initial_conditions.circle.init_circle(grid, 'f', (0.5, 1), 0.2, 1.)
out_format = 'Solution completion: {}%'
"""
for iter_no in xrange(input['maxiters']):
solver.solve()
print out_format.format((iter_no + 1.)/input['maxiters']*100.)
"""
solver._compute_curvature()
viewers.quiver(grid, 'df_x', 'df_y', show=True)
'bcs': bcs,
'rho': 1.,
'mu': 0.001,
'omega_momentum': 0.7,
'omega_p_corr': 0.3,
'advection_scheme': 'upwind',
'maxiters': 200,
'n_p_corrs': 3,
}
grid = FvRectilinearGrid((180, 60), (3, 1))
if __name__ == '__main__':
from solvers.piso import Piso
import grid.viewers as viewers
import numpy as np
solver = Piso(grid, **input)
out_format = 'Solution completion: {}%'
for iter_no in xrange(input['maxiters']):
solver.solve()
print out_format.format((iter_no + 1.)/input['maxiters']*100.)
vel, = grid.add_cell_fields('vel')
u, v = grid.get_cell_fields('u', 'v')
vel[:, :] = np.sqrt(u**2 + v**2)
viewers.display_fv_solution(grid, 'vel', show=True, show_grid=True)
