CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Which solver to use for a dynamic incompressible multiphase case _ Dead water model

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

Like Tree1Likes
  • 1 Post By alexeym

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 15, 2017, 15:33
Default Which solver to use for a dynamic incompressible multiphase case _ Dead water model
  #1
New Member
 
Marin Ducoux
Join Date: Oct 2017
Location: France
Posts: 4
Rep Power: 8
Marin Ducoux is on a distinguished road
Hi there !

As a french student, I am attempting to modelize the so-called Dead-Water Phenomenon (for more details : see https://hal.archives-ouvertes.fr/fil...ater_final.pdf).

Picture from Matthieu Mercier's experiment


Being novice with Openfoam, I've been spending some days running basic tutorials and here is my problem :

The case I'd like to set up consists of a 2D water tank, with 2 layers of water (salt one and fresh one). I managed to modelize it using interFoam.


A hull of a boat is on the top layer (fresh water). This hull has to travel from one side to the other, being towed by a constant force. It sounds a bit like the movingCone tutorial using pimpleDyMFoam. (And with some specific conditions, it should be possible to observe my phenomenon, which consist in an interfacial wave creating a pressure drop under the hull, nearly stopping the boat.)

So which solver allows to do both at the same time ? I used to have a look at interDyMFoam solver, but I didn't find any relevant case ..

So what do you think about that ? Did you ever seen a case close to that one ?

Thanks for your wise advice

Marin D.
Marin Ducoux is offline   Reply With Quote

Old   November 2, 2017, 05:54
Unhappy
  #2
New Member
 
Marin Ducoux
Join Date: Oct 2017
Location: France
Posts: 4
Rep Power: 8
Marin Ducoux is on a distinguished road
Any idea ? :/
Marin Ducoux is offline   Reply With Quote

Old   November 3, 2017, 05:45
Default
  #3
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,938
Rep Power: 38
alexeym has a spectacular aura aboutalexeym has a spectacular aura about
Send a message via Skype™ to alexeym
Hi,

Several questions:

1. Can your water layers mix? interFoam is for two immiscible fluids.

2. What is a shape of your boat? Question here is whether you can use layer addition and removal during your boat motion or it will be only mesh deformation.

3. If you need to impose drag force on your boat, maybe it would be easier to go with available IBM implementations?
alexeym is offline   Reply With Quote

Old   November 3, 2017, 07:50
Question
  #4
New Member
 
Marin Ducoux
Join Date: Oct 2017
Location: France
Posts: 4
Rep Power: 8
Marin Ducoux is on a distinguished road
Hi,

1. The two fluides are considered as immiscible.

2. not sure to understand the layer addition; the shape of my boat can be approximated as a rectangle (it's a 2D case). I'd like to set up a moving mesh, to allow having a tight mesh following the boat in order to keep the same accuracy for the wave following the boat.

3.I don't know how the IBM works on Openfoam, i'll have a look at it, thanks !

But I wondered (and it may mean I'm quite far from really understanding how solvers work) : this Dead-water phenomenon has lead scientists to establish a new equation discribing the phenomenon :

It ain't exaclty the classic equation, so will I be able to simulate the phenomenon or should I try to modify a solver in order to solve it ?

Thanks for your answer
Marin Ducoux is offline   Reply With Quote

Old   November 6, 2017, 10:31
Default
  #5
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,938
Rep Power: 38
alexeym has a spectacular aura aboutalexeym has a spectacular aura about
Send a message via Skype™ to alexeym
Hi,

1. OK. interDyMFoam is suitable for the setup.

2. Layer addition and removal is linked with mesh motion. Imagine you move your rectangle. The mesh is deformed. And as you continue rectangle motion, in front it mesh will be denser and denser, behind it mesh will become sparser and sparser. To more-or-less keep mesh density constant layer addition and removal could be added; if cell is inflated beyond certain limit, it is split into two cell with a layer of internal faces, if cell is squeezed to much, it is merged to neighbour.

Yet, if you are OK with just mesh deformation, you can try to set up case similar to tutorials/incompressible/pimpleDyMFoam/movingCone.

3. Here you can take a look. Though, it seems they implemented it for single-phase flows.

It would be nice to know what all those greek letters mean. Especially what is epsilon and why nu is constant.
Marin Ducoux likes this.
alexeym is offline   Reply With Quote

Old   November 6, 2017, 14:36
Default
  #6
New Member
 
Marin Ducoux
Join Date: Oct 2017
Location: France
Posts: 4
Rep Power: 8
Marin Ducoux is on a distinguished road
Hi, thanks for your reply !

Quote:
Originally Posted by alexeym View Post
2. Layer addition and removal is linked with mesh motion. Imagine you move your rectangle. The mesh is deformed. And as you continue rectangle motion, in front it mesh will be denser and denser, behind it mesh will become sparser and sparser. To more-or-less keep mesh density constant layer addition and removal could be added; if cell is inflated beyond certain limit, it is split into two cell with a layer of internal faces, if cell is squeezed to much, it is merged to neighbour.

Yet, if you are OK with just mesh deformation, you can try to set up case similar to tutorials/incompressible/pimpleDyMFoam/movingCone.
Ok great I see what you mean. But for now I'll try to set up the case with mesh deformation, and I will have a look at layer addition and removal as a possible improvement once I control the solver

Quote:
3. Here you can take a look. Though, it seems they implemented it for single-phase flows.
Thanks for the link, interesting method. Though, I cannot find for now (still searching) proper source code fitting openfoam5.0 implemented for multiphase cases.. Moreover it seems that the method is not really appropriate for great Reynolds number (in my case : about 1e4). Looking more deeply for solutions.


Quote:
It would be nice to know what all those greek letters mean. Especially what is epsilon and why nu is constant.
Sorry, my bad :/.
-> Eta(x,t) is the spread pf the deformation of the surface
-> Xi(x,t) is the spread og the deformation of the interface (fresh/salt water)
-> Alpha is the (max deformation of the surface)/(max deformation of the interface)
-> Epsilon is the (max deformation of the interface)/(thickness of the top layer
-> Lambda is the characteristic length of the wave
-> Nu is the (thickness of the top layer)/(Lambda)
(alright, the notation is maybe not the most appropriate x) )

In short, I guess those parameters ain't that usefull for the comprehension of the modelization of the problem.

I try to run the case while studying in parallel the physics/applied maths of the model, so that I'm not able yet to clearly make simplifying hypotheses :/

So, unless I'm mistaken, InterDyMFoam should match with my problem ? Even with the last "Forcing" term of the equation which is quite new ? :/

(Moreover, apologies for my English, I hope you can understand me)
Best regards
Marin Ducoux is offline   Reply With Quote

Reply

Tags
dead water phenomenon, interdymfoam, interfacial wave, movingcone, openfoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
CFD analaysis of Pelton turbine amodpanthee CFX 31 April 19, 2018 18:02
fluent divergence for no reason sufjanst FLUENT 2 March 23, 2016 16:08
Is Playstation 3 cluster suitable for CFD work hsieh OpenFOAM 9 August 16, 2015 14:53
Overflow Error in Multiphase Modelling with Two Continuous Fluids ashtonJ CFX 6 August 11, 2014 14:32
Free surface boudary conditions with SOLA-VOF Fan Main CFD Forum 10 September 9, 2006 12:24


All times are GMT -4. The time now is 02:43.