CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   twoPhaseEulerFoam - sudden enlargement of circular pipe validation case (http://www.cfd-online.com/Forums/openfoam-solving/110856-twophaseeulerfoam-sudden-enlargement-circular-pipe-validation-case.html)

yanxiang December 22, 2012 19:41

twoPhaseEulerFoam - sudden enlargement of circular pipe validation case
 
1 Attachment(s)
Hi all,

In Rusche's thesis, he validated the two-fluid method with the sudden enlargement of circular pipe case. Could anyone send me a copy of this original experimental paper? I can't find it through google or our library.

Bel Fdhila R, Simonin O. Eulerian prediction of a turbulent bubbly ow downstream of a sudden pipe expansion. In Proceedings of the 6th Workshop on Two Phase Flow Predictions, Sommerfeld M (ed), University of Erlangen, 1992)

I followed whatever is available in Rusche's thesis, but I got a high gas fraction zone right after the enlargement (see attachment).

Thanks,
yanxiang

alberto December 23, 2012 05:09

Are you using the latest 2.1.x code with MULES?

yanxiang December 23, 2012 11:55

Nope. This is OF211. Do you think it would be different if I I use 2.1.x code? Nevertheless, I think, OF211 shouldn't behave like that either.

alberto December 24, 2012 05:42

The MULES version of twoPhaseEulerFoam should do a much better job at ensuring that alpha is bounded, which is less robust in the old version. Maybe give it a try :-)

yanxiang December 24, 2012 12:34

2 Attachment(s)
Absolutely. I gave it a try and the result is attached. Gas still tends to go to the corner. The velocity field looks fine though. I have also attached the case tar ball for your review.

Thanks,
yanxiang

sharonyue December 24, 2012 23:59

Quote:

Originally Posted by alberto (Post 398868)
The MULES version of twoPhaseEulerFoam should do a much better job at ensuring that alpha is bounded, which is less robust in the old version. Maybe give it a try :-)

Sorry for hijacking this thread.But I really need you help,Prof.Alberto.
Much details has been depicted here
http://www.cfd-online.com/Forums/ope...ured-mesh.html

Look forward to you replay.Thanks in advance.

alberto December 26, 2012 09:36

Quote:

Originally Posted by yanxiang (Post 398916)
Absolutely. I gave it a try and the result is attached. Gas still tends to go to the corner. The velocity field looks fine though. I have also attached the case tar ball for your review.

Thanks,
yanxiang

I don't see problems in the setup of your case. The accumulation of gas in the corner does not seem too different from Rusche's results (Fig. 3.3 c, page 131).

It is possible you have less numerical diffusion than in those results (it seems so from your picture). You should run with the same schemes to have an apple to apple comparison.

Best,

alberto December 26, 2012 10:07

Quote:

Originally Posted by sharonyue (Post 398943)
Sorry for hijacking this thread.But I really need you help,Prof.Alberto.
Much details has been depicted here
http://www.cfd-online.com/Forums/ope...ured-mesh.html

Look forward to you replay.Thanks in advance.

I commented on the other thread.

yanxiang January 2, 2013 12:39

Hi Alberto,

I tried the same schemes as used by Rusche, but the results stay pretty much unchanged, with a high gas zone in the corner, and this is not there in Rusche's results, or in Oliveira's (Int. J. Numer. Meth. Fluids 2003; 43:1177–1198). I played around with the BC's, but still couldn't get rid of that. Also, sadly, I was not able to find any of the references (except for Oliveira's) used by Rusche in that particular test case section where others validated their models with it. Do you have any of those papers?

Thanks,
Yanxiang

P.S: Happy New Year!!!!

alberto January 2, 2013 12:50

How did you decide the BC's for the turbulent quantities? They aren't provided in the article...

alberto January 2, 2013 13:22

1 Attachment(s)
I ran the case you attached turning turbulence off, just to see the effect, and segregation is much lower, as expected.

Anyways, I am working on a validation case for twoPhaseEulerFoam/multiphaseEulerFoam for bubbly flows... we'll see :-)

yanxiang January 2, 2013 16:18

4 Attachment(s)
Alberto,

So after reading your question about the turbulent quantities, I made some changes to my very original test case at work, and ran the OF211 version of tPEF on it. Weirdly, it worked :-$. By that I mean, I got similar results to those in Rusche's and Oliveira's works. I will try to reproduce that with OF21x version when I get home (although a diff on the two cases didn't give me any hints why the results would differ). Anyways, I attached the screenshots and case for OF211.

Best,
yanxiang

alberto January 2, 2013 17:18

There is a significant difference in the two algorithms used to solve for alpha in twoPhaseEulerFoam before and after the introduction of MULES, which might explain the differences. I am running some validation cases exactly to see if accuracy was maintained (for now, it is, at least in the case of a simple bubble column).

sharonyue January 2, 2013 20:01

Quote:

Originally Posted by alberto (Post 399795)

Anyways, I am working on a validation case for twoPhaseEulerFoam/multiphaseEulerFoam for bubbly flows... we'll see :-)

I am eager to see validation cases for twoPhaseEulerFoam on tet mesh~:rolleyes:

alberto January 2, 2013 21:06

Quote:

Originally Posted by sharonyue (Post 399830)
I am eager to see validation cases for twoPhaseEulerFoam on tet mesh~:rolleyes:

I am not considering tetrahedral grids, because my geometry does not require them (it is a box with rectangular section!). I consider this test-case because it is the precursor for new model developments.

If you want to perform a verification and validation study on tetrahedral meshes, however, it should not be that hard. You should consider a simple geometry with well-defined results (there are many on bubble columns in the literature, pick a case where simulation results are also available, and check the model is the same), do a grid-independence study, and then do the experimental validation on the converged grid.

P.S. Have you tried to use more appropriate numerical schemes, as I suggested in the other thread, with a good quality mesh?

sharonyue January 2, 2013 22:04

5 Attachment(s)
Quote:

Originally Posted by alberto (Post 399832)
I am not considering tetrahedral grids, because my geometry does not require them (it is a box with rectangular section!). I consider this test-case because it is the precursor for new model developments.

If you want to perform a verification and validation study on tetrahedral meshes, however, it should not be that hard. You should consider a simple geometry with well-defined results (there are many on bubble columns in the literature, pick a case where simulation results are also available, and check the model is the same), do a grid-independence study, and then do the experimental validation on the converged grid.

P.S. Have you tried to use more appropriate numerical schemes, as I suggested in the other thread, with a good quality mesh?


I am sorry about that question. because it looks like there is no substantial change to the result.but I dont know if its still because of the mesh.


So I refine the mesh, now the cell is half than that cell in the other thread.so the total cell number reachs to 1,79 million,its too many, so I shorten the height, the other thing is not changed.the total cells number is about 55000.

I attached my image.I think that dont need to depict it.My primary language is not english so~.

btw,I turn to the other thread,I am sorry for that.
http://www.cfd-online.com/Forums/ope...ured-mesh.html

Thanks for you consistent attention.Alberto

yanxiang January 4, 2013 00:22

Hmmm... interesting. I tried the OF21x version on the same case again, and as you would expect, the results just look like what we had previously. So does that mean the new solver with the MULES method is not solving the alpha equations correctly?

Thanks,
yanxiang


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