[dealexander]

Has anyone written any code for the 2D shallow water equations? I have a dam-failure problem I would like to simulate using a finite volume technique. I've never used OpenFOAM, but it looks like if might be fun. Any advice

regards

David E. Alexander M.Sc. P.Eng
3-Sigma Consulting Inc.
david@3sigma.ca
www.3sigma.ca

[liuzhw]

I have just finished the code for the 2D shallow water equations. but it is tested only by a simple uniform channel flow. To following the rule of dimensionSet, I have the whole water depth as non-Uint. It is a little unsatisfied.
If you need, I will share it with you.

[taufik]

I hope this thread is still on. I am also in the initial phase of developing solver for shallow water. Would appreciate if you could share the code to me as well.

[liuzhw]

Hello!
I have developed a depth-averaged shallow water equation model based on Foam.
It is a two dimensional problem and the top and bottom suface is spesified as
emptyFvPatch.

however, the velocity along the z direction, U.z(), will sometime increase to a large magnitude in the computation. it is not what I expected.

now, I have to limte U.z() to zero after iteratation each time step.

any advices for this problem are grateful!

best regard!

[liuzhw]

Hello!
I have developed a depth-averaged shallow water equation model based on Foam.
It is a two dimensional problem and the top and bottom suface is spesified as
emptyFvPatch.

however, the velocity along the z direction, U.z(), will sometime increase to a large magnitude in the computation. it is not what I expected.

now, I have to limte U.z() to zero after iteratation each time step.

any advices for this problem are grateful!

best regard!

[liu]

For shallow water equation, you need to solve h (water depth), u and v (velocity in x and y direction). Do you solve them together or seperately? Or is it coupled?

For the current version of Foam, I don't think it is easy to solve shallow water equations.

Another thing is you need to implement Godunov schemes and flux routines.

Good luck!

[liuzhw]

Yes, I solve h with the water surface equation modified from continuity equation. the governing equations are solved with the PISO iterating method just like the turbForm solvers in Foam.

Now it runs well except the problem of the increasing magnitude of U.z(). although the problem does not affect the result, it is not nice.

Thanks

[liuzhw]

Next, I plan to develop a three dimensinal model for shallow water equations using sigma coordinate along the verticle direction. since the mesh of each layer is same, there are two mesh system to be chosen. one is 2D mesh and the other is a true 3D mesh.

for 2D mesh, the code will be simple but with a complicated data structure of velocity, the data structure like List<volvectorfield> will be used. more ever, it is not easy to expand a real 3D model

for 3D mesh, it will like a solver of most of foam. but the solution to water surface needs 2D mesh, which is a problem to be overcome.

any advices for this problem are grateful!

best regard!

[hahuhe]

Originally Posted by liuzhw

Yes, I solve h with the water surface equation modified from continuity equation. the governing equations are solved with the PISO iterating method just like the turbForm solvers in Foam.

Now it runs well except the problem of the increasing magnitude of U.z(). although the problem does not affect the result, it is not nice.

Thanks

it happens a lot to shallow water equation. especially when there are two waves from opposite direction. such as the incoming wave colliding with the reflected wave. the model is unable to handle the turbulence from that process. instead of adding and solving the turbulence term, i found that simple numerical filter such as Hansen (1962) really help. you can find the explanation in Kowalik's book (Numerical modeling of Ocean Dynamics vol.5)

[zhaoshiyu]

Hello. I am quite interested in your shallow water foam, I want to simulate a channel water flow from the top of slope, but it is quite difficult to use the original one in openFoam. I will appreciate it so much if you can share me with your shallow water foam code and case. Thanks.

[zhaoshiyu]

Hello. I am quite interested in your shallow water foam, I want to simulate a channel water flow from the top of slope, but it is quite difficult to use the original one in openFoam. I will appreciate it so much if you can share me with your shallow water foam code and case. Thanks.

[zhaoshiyu]

Originally Posted by liuzhw

Next, I plan to develop a three dimensinal model for shallow water equations using sigma coordinate along the verticle direction. since the mesh of each layer is same, there are two mesh system to be chosen. one is 2D mesh and the other is a true 3D mesh.

for 2D mesh, the code will be simple but with a complicated data structure of velocity, the data structure like List<volvectorfield> will be used. more ever, it is not easy to expand a real 3D model

for 3D mesh, it will like a solver of most of foam. but the solution to water surface needs 2D mesh, which is a problem to be overcome.

any advices for this problem are grateful!

best regard!

Hello. I am quite interested in your shallow water foam, I want to simulate a channel water flow from the top of slope, but it is quite difficult to use the original one in openFoam. I will appreciate it so much if you can share me with your shallow water foam code and case. Thanks.