Problem in implementation of Shallow water equation
I am trying to implement shallow water equations as described in Finite-Volume Model for Shallow_water Flooding of Arbitrary Topography(Bradford and sanders). I got matlab code from the site of Sanders( http://sanders.eng.uci.edu/matlabcodes.html). For 2D it didn't take into account for bed elevation and wet/dry boundary. I extended the implementation to take account for bed elevation and wet/dry boundary. But the system gets unstable after sometime(Courant no. increases more than 1).
For Wet/Dry Boundary,
instead of F = 0.5(FL + FR - R*A*dW), equation for equating fluxes at the interface
i m doing , F = 0.5*(FL + FR), for extrapolation, is it correct?
and when height goes below 0.0, i m setting F = 0;
Is there anything else i need to do? Can anyone guide me how to apply conditions for wet/dry boundary? or provide me link to the code in which case of wet/dry boundary is implemented.
Thanks in advance
did you manage to solve the problem with the wetdrying front? Currrently, I'm working on that problem, too. So I would be happy, if you could share your progress!
|All times are GMT -4. The time now is 16:38.|