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

How to simulate the gravitydriven RayleighTaylor instability

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 18, 2005, 09:08
Default What is the main difficulty in
  #1
Member
 
Luckyluke
Join Date: Mar 2009
Posts: 51
Rep Power: 8
luckyluke is on a distinguished road
What is the main difficulty included in this problem? I was told that this problem is hard.
Can I solve it using InterFoam?
luckyluke is offline   Reply With Quote

Old   January 6, 2009, 10:53
Default I have set the gravity-driven
  #2
New Member
 
Thibault Pringuey
Join Date: Mar 2009
Posts: 17
Rep Power: 8
thibault_pringuey is on a distinguished road
I have set the gravity-driven RT instability using the funkySetFields utility to give the initial sine shape to the interface and running the case with interFoam.

Unfortunately, my solution does not compare satisfactorily with previous work (Lopez et al. 2005 jcp, Popinet: http://gfs.sourceforge.net/examples/.../rt.html#htoc5) and show additional "ligaments"along the interface.

Modifying the thickness of the interface, the grid refinement and/or the surface tension did not allow to match the references given above.

Has anyone obtained a satisfactory solution on this problem?
thibault_pringuey is offline   Reply With Quote

Old   January 12, 2009, 19:37
Default Thibault, How does this loo
  #3
egp
Senior Member
 
egp's Avatar
 
Eric Paterson
Join Date: Mar 2009
Location: Blacksburg, VA
Posts: 197
Blog Entries: 1
Rep Power: 9
egp is on a distinguished road
Thibault,

How does this look?




This solution was on a fine grid of 256x1024. In general, I don't like this problem as a validation test case since there is neither an analytical solution or experimental data. In addition, the variation in published results is pretty large.

However, it is good from the point of view that it is a good tutorial on the use of funkySetFields, blockMesh, and interFoam!
egp is offline   Reply With Quote

Old   January 12, 2009, 19:43
Default Hmmm, the image-posting gremli
  #4
egp
Senior Member
 
egp's Avatar
 
Eric Paterson
Join Date: Mar 2009
Location: Blacksburg, VA
Posts: 197
Blog Entries: 1
Rep Power: 9
egp is on a distinguished road
Hmmm, the image-posting gremlin....

Try these links:

gamma field

pd field

I did do a grid study (256x1024, 128x512, 64x256, 32x128), and the coarser meshes were unable to support thin strands (as would be expected).

Also, surface tension had a pretty big impact on the solution. I used a value close to that for helium & air.

I'll post the cleaned up case later.
egp is offline   Reply With Quote

Old   January 13, 2009, 08:16
Default The case can be downloaded fro
  #5
egp
Senior Member
 
egp's Avatar
 
Eric Paterson
Join Date: Mar 2009
Location: Blacksburg, VA
Posts: 197
Blog Entries: 1
Rep Power: 9
egp is on a distinguished road
The case can be downloaded from:

fine grid case for the Rayleigh-Taylor problem

Don't forget to run blockMesh. There is a funkySetFieldsDict, but you shouldn't have to run it (0/gamma has the sin-wave I.C.).

Have fun...
egp is offline   Reply With Quote

Old   January 13, 2009, 10:28
Default Thank you very much Eric! Your
  #6
New Member
 
Thibault Pringuey
Join Date: Mar 2009
Posts: 17
Rep Power: 8
thibault_pringuey is on a distinguished road
Thank you very much Eric! Your solution looks great.

I will study your case files to identify the issue in mine.

BTW, could please let me know how you included links to images in your message?

Many thanks.
thibault_pringuey is offline   Reply With Quote

Old   January 13, 2009, 12:51
Default Thibault, See the "Formatti
  #7
egp
Senior Member
 
egp's Avatar
 
Eric Paterson
Join Date: Mar 2009
Location: Blacksburg, VA
Posts: 197
Blog Entries: 1
Rep Power: 9
egp is on a distinguished road
Thibault,

See the "Formatting" link on the left hand side of the Message Board browser window? If you follow this link, you can learn how to do all sorts of formatting, like changing the color, <u>underlining text</u>, and making text annoyingly blink!. There are also the commands,

\link
\image

However, the \image command did not work for me the last time I tried. Maybe I was using Safari instead of Firefox.

Eric
egp is offline   Reply With Quote

Old   January 16, 2009, 06:49
Default Hello Eric, The differences
  #8
New Member
 
Thibault Pringuey
Join Date: Mar 2009
Posts: 17
Rep Power: 8
thibault_pringuey is on a distinguished road
Hello Eric,

The differences with my simulation concerned the solvers and schemes employed. Thank you for making me realize this.

I have run the simulation on a 128x512 and a 256x1024 mesh with your settings. I have noticed that there are some "wiggles" along the interface in the region of the neck of the mushroom-shaped gamma distribution. These "wiggles" are amplified on the 128x512 mesh.

Would you know what is causing these "wiggles"? Are they due to spurious currents?

Many thanks.


Thibault
thibault_pringuey is offline   Reply With Quote

Old   January 19, 2009, 13:27
Default Hi Thibault and Eric, I jus
  #9
Senior Member
 
Michael Jaworski
Join Date: Mar 2009
Location: Champaign, IL, USA
Posts: 126
Rep Power: 8
mike_jaworski is on a distinguished road
Hi Thibault and Eric,

I just thought I'd drop a note about validation for this instability. There are analytical treatments of the growth rate of the instability when it is still a "small" disturbance and the initial equilibrium can be linearized. The growth rate depends on the wave number of the disturbance (the classic text on this stuff is probably Chandrasekhar's "Hydrodynamic and Hydromagnetic stability"). This could provide validation of the solver for small-scale disturbances before the problem goes non-linear.

For the non-linear growth etc, I'm not sure about the best approach to validating this except by comparison to literature and experiments.

Best,
Mike J.
mike_jaworski is offline   Reply With Quote

Old   January 19, 2009, 14:50
Default @Thibault: I also saw the int
  #10
egp
Senior Member
 
egp's Avatar
 
Eric Paterson
Join Date: Mar 2009
Location: Blacksburg, VA
Posts: 197
Blog Entries: 1
Rep Power: 9
egp is on a distinguished road
@Thibault: I also saw the interface wiggles in the neck region. It is in the high-velocity region of the flow, and is sensitive to the value of surface tension. Without more study, I don't think I can ascribe it to either numerical instability or something physical such as capillary waves.

@Mike: Good suggestion on testing the linear regime of the instability. However, I'll leave that to someone else! Maybe a young energetic student can chew on this as part of the verification section of their thesis.
egp is offline   Reply With Quote

Old   March 3, 2009, 07:49
Default Hello. Talking about "wiggl
  #11
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
Hello.

Talking about "wiggles".
Do you mean something like this?



This is a calculation using interDyMFoam.

And what are these strange non-orthogonal cells in the region some distance away from the interface? Or is this due to the visualization within ParaView?
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   March 3, 2009, 07:56
Default Well hope it will work this ti
  #12
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
Well hope it will work this time.
Like this:


__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   March 5, 2009, 06:34
Default This behaviour was due to the
  #13
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
This behaviour was due to the usage of CrankNicholson 1 for ddtSchemes.

Using the Euler Scheme solved the problem.
Have a look:


__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
KH instability kumar Main CFD Forum 0 December 9, 2006 19:30
1D shock instability with Roe FDS? ma Main CFD Forum 4 January 22, 2005 00:10
crow instability A.Jaafar Main CFD Forum 3 October 30, 2000 07:19
Nonlinear instability Yap Wen Jiun Main CFD Forum 3 September 12, 2000 23:16
Instability of convection Zeng Main CFD Forum 0 May 26, 2000 05:39


All times are GMT -4. The time now is 12:29.