# Shallow water + VOF

 Register Blogs Members List Search Today's Posts Mark Forums Read

 May 21, 2011, 02:29 Shallow water + VOF #1 New Member   Mitsuru Nishikawa Join Date: May 2011 Posts: 3 Rep Power: 7 I am new to this forum and beginner as a developper of CFD. Please help me in the development of openfoam solver. I am intended to merge the codes of shallow water equation and of VOF model in order to treat shallow water equations with free surface which is not height (included in the original equation) but width. (That is, shallowWaterFoam+interFoam) First, I think it is so simple, but in fact it is difficult for the difference in treatment of PISO loop or so. Have you ever treated such a case? Or, if you know the below questions, please teach me references? (1) detailed solving steps in shallowWaterFoam or interFoam (particularly, in momentum(+h)-height coupling in shallowWaterFoam) (2) In such a treatment, the equations of alpha (da/dt+v.div(a)=0) should be modified? Sincerely, MIT

 May 23, 2011, 02:04 #2 Senior Member   Bernhard Join Date: Sep 2009 Location: Delft Posts: 790 Rep Power: 14 I don't understand how this combined solver should look like? What physics do you want to describe that VOF cannot? I suppose VOF + shallowWaterFoam = VOF? Furthermore, I am missing the point why your free surface is width, and not height, what is the difference? For reference on interFoam, you can look here: http://foamcfd.org/resources/theses.html for the PhD theses of Ubbink and Rusche. Or the PhD thesis of Berberovic, Google "Investigation of Free-surface Flow Associated with Drop Impact: Numerical Simulations and Theoretical Modeling"

 May 23, 2011, 03:30 #3 New Member   Mitsuru Nishikawa Join Date: May 2011 Posts: 3 Rep Power: 7 Thank you for your reply. At first, I apologize for the lack of my explanations. And because I do not have much experience in CFD, there should be a lot of my mistakes. It will give me much pleasure if you point out them. The physics I want to describe is the flow of thin liquid films, paticularly how it spreads spatially (I defined Height as thickness of films, and Width of film width. Height is included in the schemes of SWE, I want to treat the free surface of the edge of film by VOF). As you say, the only VOF (interFoam) is enough to treat a free-surface flows even in the case of the thin liquid film, so physics does not change. However, I want to couple shallow water equations and VOF due to the two reasons, (1) the problems of high aspect ratio in height/width due to thin films (2) the reduction of computational cost (3D to 2D) Yesterday, I merged SWE and VOF, but I cannot be sure it works correctly. Therefore, I am looking of its theoretical backbones (or mistakes) or a previous researches. Sincerely, MIT

 May 26, 2011, 09:37 #4 New Member   Robin Koldeweij Join Date: Nov 2009 Posts: 19 Rep Power: 8 Could you may be explain more carefully what you're trying to do physically? I.e. start with the navier stokes equations, which assumptions do you want to make and what do you want to solve. Just merging two solvers won't bring you the joy of a physically correct solution.. May be you can just use the shallow water equations for a thin liquid film along a wall, with an extra source term for gravity effects. You'll have to think about relative viscosities and densities however. It might be I don't understand correctly, but it seems you have an idea of what to solve, but haven't looked into the physics/assumptions of both solvers.

 May 26, 2011, 20:27 #5 New Member   Mitsuru Nishikawa Join Date: May 2011 Posts: 3 Rep Power: 7 Dear RBJ, Thank you for your suggestion. And I apologize my insufficient explanations. Of course, I haven't simply merged the two codes. In swallowWatterFoam(SWF), viscosity term etc are not included. I omited the details of interactions I included because it is cumbersom to start from Navier-Stokes and it is not essential for my questions. But I write down the interations I want to inculde. (Momentum equations:time derivative of momentum,inertia,viscosity,gravity,surface tension, pressure(pressure is displaced with the derived equation of height)) and the morphology of free surface. Please image 3D plate, in which vertical length is much shorter than lateral length. Of course, it has 6 surfaces.Tｗo of them are perpendicular to vertical directions, and the four are to lateral. I mentioned the vertical ones are included in SWE scheme as height(h), and I want to treat the remainders by VOF. SWF is able to treat the free-surface of faces perpendicular to height directions, isn't it? So, ocean wave in Picked-up area of sea which is surrounded by sea can be treated. But is it able to treat movement of surf in sand beach only by SWF? (If it is correct, is the area of h!=0 regarded as the area where water is? I can't understand. If so, free surface(what I say as a faces of width) is able to be treated as the boundary between(h=0 or h!=0) and surface tensions is calculated from the curvature of this boundary, my purpose is accomplished). Last edited by MIT; May 26, 2011 at 21:03.

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post RicardoGoncalves Main CFD Forum 0 July 27, 2010 04:38 Pathway FLUENT 1 July 21, 2007 07:33 Tysor FLUENT 0 March 17, 2006 04:10 Paul Main CFD Forum 10 August 30, 2004 11:56 fredius,magige,tanzanian,(e.a) Main CFD Forum 0 January 27, 2002 08:10

All times are GMT -4. The time now is 18:39.