[luckyluke]

What is the main difficulty included in this problem? I was told that this problem is hard.
Can I solve it using InterFoam?

[thibault_pringuey]

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?

[egp]

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]

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]

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...

[thibault_pringuey]

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.

[egp]

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

[thibault_pringuey]

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

[mike_jaworski]

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.

[egp]

@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.

[sega]

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?

[sega]

Well hope it will work this time.
Like this:

[sega]

This behaviour was due to the usage of CrankNicholson 1 for ddtSchemes.

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

[kk415]

Hi all


I am struggling in solving this Rayleigh-Taylor Stability problem. My results I am attaching for the reference.


I read all the the documents available in the net and tried all the boundary condition suggested but still getting same answer. I am using OpenFoam 4.1. where the VOF is solved algebraically. Is that the reason for the error? IsoAdvector is a must for such problems ?


My domain 1 * 4 divided in 64 *256. I have also tried 128 * 512 and I am getting the same result in fine mesh as well.