1D second order upwind
I am working on modeling hydrodynamic flow using a 2nd order upwind scheme (2nd order Runge-Kutta for the temporal discretization). The code mainly involves solving the Saint-Venantís equation (SWE) for a one-dimensional flow with a given boundary condition. Here I have attached the code. While running the model it has a problem both in stability and accuracy. But, according to the classical boundary condition that I used, it has to converge to a uniform 1m flow depth throughout the profile.
Could anyone help me please?
|All times are GMT -4. The time now is 22:33.|