||April 14, 2010 07:05
problem when simulating overflow at a sharp-crested weir
I am trying to reproduce the overflow at a sharp-crested vertical weir but experiencing problem. The computational domain is sketched in the attached picture. When the overflow jet hit the right bottom, the air region (identified as D in the pic) enclosed by weir, water jet and ground will be disconnected from the atmosphere in 2D case. At the mean time, the jet will split into two part basically, one of them (the left-going part) will climb the weir and collapse finally. The problem is, the total area of the closed air region D will not be likely to change, but as the water increases, the area occupied by air will have to decrease, but this is not quite possible neither for interFoam nor physical facts. Therefore, I changed the boundary condition of the edge I to pressureInletOutlet for U and Totalpressure = 0 for P. And I applied non-slip wall condition on all other wall surfaces. It seems doesn't solve the problem and the solver will always crash soon after the moment the snapshot shows. Please help to provide some suggestions on addressing this problem, many thanks!