'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)