
[Sponsors] 
November 7, 2005, 19:44 
Has anyone written any code fo

#1 
New Member
david e. alexander
Join Date: Mar 2009
Location: North Vancouver, BC, Canada
Posts: 1
Rep Power: 0 
Has anyone written any code for the 2D shallow water equations? I have a damfailure 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 3Sigma Consulting Inc. david@3sigma.ca www.3sigma.ca 

November 8, 2005, 21:01 
I have just finished the code

#2 
New Member
Liu Zhaowei
Join Date: Mar 2009
Posts: 9
Rep Power: 9 
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 nonUint. It is a little unsatisfied.
If you need, I will share it with you. 

April 18, 2006, 03:35 
I hope this thread is still on

#3 
New Member
Achmad Taufik
Join Date: Mar 2009
Posts: 2
Rep Power: 0 
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.


January 3, 2007, 02:16 
Hello!
I have developed a dep

#4 
New Member
Liu Zhaowei
Join Date: Mar 2009
Posts: 9
Rep Power: 9 
Hello!
I have developed a depthaveraged 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! 

January 3, 2007, 02:20 
Hello!
I have developed a dep

#5 
New Member
Liu Zhaowei
Join Date: Mar 2009
Posts: 9
Rep Power: 9 
Hello!
I have developed a depthaveraged 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! 

January 3, 2007, 16:33 
For shallow water equation, yo

#6 
Senior Member
Xiaofeng Liu
Join Date: Mar 2009
Location: State College, PA, USA
Posts: 118
Rep Power: 9 
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!
__________________
Xiaofeng Liu, Ph.D., P.E., Assistant Professor Department of Civil and Environmental Engineering Penn State University 223B Sackett Building University Park, PA 16802 Web: http://water.engr.psu.edu/liu/ 

January 3, 2007, 20:51 
Yes, I solve h with the water

#7 
New Member
Liu Zhaowei
Join Date: Mar 2009
Posts: 9
Rep Power: 9 
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 

January 3, 2007, 21:33 
Next, I plan to develop a thre

#8 
New Member
Liu Zhaowei
Join Date: Mar 2009
Posts: 9
Rep Power: 9 
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! 

April 8, 2009, 01:55 

#9  
New Member
Adityawan, Bagus
Join Date: Apr 2009
Posts: 2
Rep Power: 0 
Quote:


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
shallow water problem  seeking help  FLUENT  2  December 5, 2011 10:05 
shallow water equations  Nalan ANtar  Main CFD Forum  2  February 28, 2011 06:03 
Looking for stable integration scheme for 2D Shallow Water Equations  maddhi  OpenFOAM Running, Solving & CFD  1  August 23, 2008 06:24 
shallow water kepsilon equation  Ali  Main CFD Forum  0  January 5, 2006 16:56 
shallow water VS deep water  Paul  Main CFD Forum  10  August 30, 2004 11:56 