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 and LES (http://www.cfd-online.com/Forums/openfoam-solving/63667-twophaseeulerfoam-les.html)

raagh77 April 15, 2009 14:43

twoPhaseEulerFoam and LES
 
Hello all,

anyone tried implementing LES in twoPhaseEulerFoam solver?


awaiting for replies.. :)

Regards
Raghavendra

alberto April 16, 2009 02:13

Hi,

yes, I worked on that during my PhD. I essentially used Smagorynsky type models for the gas phase to simulated gas-particle flows. The work was done with OF 1.4.0, and does not compile with more recent versions of OpenFOAM. I did not port it to OpenFOAM 1.5.x for lack of time (I'm working on a different topic right now).

What I can suggest to you is to reuse the single phase classes for LES, and adapt them at the multiphase case, which is essentially what I did.

Best regards,

juho April 16, 2009 03:45

Hi,

I actually started modifying the single phase incompressible LES-models for use with twoPhaseEulerFoam yesterday evening. Based on Albertos' example. Didn't get that far with it yet though.

I'd be interested in including some sort of model to account for modulation of turbulence due to particles in them. (gas-solid flows)

Regards,
Juho

raagh77 April 16, 2009 04:32

Hi Alberto and Juho,

Thanks for your replies..Its good to know about LES implementation in twoPhaseEulerFoam.

Are there any references about LES in mutiphase flows which I can go through before working on it?

Is it possible to share the LES implementation for twoPhaseEulerFoam (in OF-1.4) and I can continue from there to implement the same in OF-1.5 ..?


Regards
Raghavendra

alberto April 17, 2009 02:28

Hi Raghavendr,

I have no problem to share the code, but I don't really know how useful it can be for your purposes, because starting from the single-phase LES classes in OpenFOAM and adapting them to the multiphase case is not a long work (1-2 days should be enough if you know OF a bit). You mainly need to rename the classes, make some changes to the submodels according to the literature (for gas-particle flows there are many paper, tell me if you need references) and make some minor change to the twoPhaseEulerFoam to replace the turbulent viscosity with the SGS viscosity.

An Juho said, it would be necessary to account for the effect of particles on the sub-grid phenomena description (See for example the very interesting work of Prof. Simonin on this purpose), and also include the effect of sub-grid phenomena on the dispersed phase (drag correction). This was not done in my code.


Best,

raagh77 April 18, 2009 16:45

Hi Alberto,

Yes, it seems to be quite simple and just to change the way nutb is calculated to implement LES..
but
I am very much new to OF and I haven't tried modifying the source code. Till now (past 2months) I was just using the available solvers :)

If its not a problem is it OK to share the LES code you implemented for twoPhaseEulerFoam?..
are there any references available with you regarding implementation of LES in twoPhaseEulerFoam for liquid gas flows?..
(Similar to Phd thesis of Henrik Rusche what you mentioned earlier for K-epsilon model and twoPhaseEulerFoam)

alberto April 19, 2009 16:16

Hi, I don't see any contact information in your profile where I can send the code.

Best,

raagh77 April 19, 2009 16:22

ohh!!

u can send it to me through

raagh77@gmail.com

raagh77 April 19, 2009 16:25

in the old forum my mail id was visible in myy profile..
don't know why its not showing in the new forum!

alberto April 20, 2009 01:14

Check you email. :-)

About your address, I have no idea. Also mine should be visible, and it should be possible to download a vCard, but the forum is so bloated I didn't find it anywhere! :-P

raagh77 April 20, 2009 03:08

Hi Alberto,

Yes, I received your mail..
Thanks :)

I made changes in my profile (which I made in the old forum also) so now its fine :)
thanks for notifyiong me and also for the code


Regards
Raghavendra

gwierink April 21, 2009 02:23

Hi Alberto and Raghavendra,

I am playing around with twoPhaseEulerFoam as well and I am interested in using LES with that solver. Could one of you please send me the code as well? It would help me a lot!

Many thanks in advance!

Regards,

Gijsbert

gijsbert.wierink@gmail.com

Edit: Thanks Raghavendra, I received the file. Now for some homework to get it to compile nicely on OpenFOAM-1.5.

alberto April 22, 2009 20:27

It seems Raghavendra sent it to you without even asking me...

Just to clarify, that code is NOT released publicly, so I would like to know who is using it.

gwierink April 23, 2009 00:57

Hi Alberto,

Yes, Raghavendra has sent me the code. I understand you want to know what goes on with it and where, and it's not so nice if other people throw your stuff around. I am sure no harm was intended. As for me, I am developing a model to simulate three-phase stirred tank reactors (continuous liquid with dispersed bubbles and solids) for my PhD. The issue with models in literature is that most of them simulate a (passive) scalar field on top of just a single phase. Another big thing is that commonly the k-epsilon model is used, so that turbulence is isotropic. For bubble-particle interaction in my case the anisotropy, in my opinion, is significant. Hence my interest in a LES two-phase solver. I have understood that your PhD thesis is on CFD of two-phase flows. Would it be possible to obtain your thesis as a pdf from somewhere, or should I perhaps ask your former university about it? (I sent you a message about this already, sorry if you get it twice ...). Your thesis would be very useful and of course, if I use your solver I will refer to your work.

For now I'm having some compilation trouble with twoPhaseLESEulerFoam, probably due to a difference in versions of OpenFOAM. I'll let you know if I come up with something interesting :). And please inform me if you would just not like me to use your solver the way it is. If you are interested I can send you a more detailed mail about what I do and my progress, in any case your feedback is very welcome.

Have a great day! :)

Best regards,

Gijsbert

raagh77 April 23, 2009 02:55

@ Alberto
 
I had forwarded LEStwophaseEulerfoam code to him! (the same mail, without asking you)
am sorry I should have asked you before as you are the sole owner of the code..

:confused:

Regards
Raghavendra

alberto April 23, 2009 08:21

Hello,

no problem. I simply wanted to clarify the code is not published. The reasons are many:

a) it is not up to date
b) some parts are not tested (see interstitial fluid effects, for example)
c) there are some quite rough edges (instabilities when the granular temperature goes to alphaMax, limitations on SGS models (no account for particle influence, no SGS model for gas-liquid flow)).
d) finally, I would like to keep track of what goes on with the code to eventually put all the efforts together and avoid repeating work done by others ;-)

Just an idea: why don't we keep in touch about the progress you have? I should work soon on the basic solver again to make it more stable.

Best,
Alberto

raagh77 April 24, 2009 08:26

Yes,
thats a good idea..
(is that through Forums we can be in touch? or anyother means)

Hopefully we can see twoPhaseLESEulerFoam in the next version of OF :)

alberto April 25, 2009 01:08

Hi Raghavendra,

to keep in touch email is OK I think.

Best,

piterw82 January 28, 2010 11:48

hello,

I am new in OpenFOAM and I am starting my work on LES multiphase with EE approach. I see you guys are already working on that topic and probably have a lot of experience in that field. I was wondering Alberto if it would be possible to have a look on your code? I read it is not checked yet but I am sure it would help me to understand how I can make some changes inside the source terms in OpenFOAM. Thank you in advance.

Peter

alberto January 28, 2010 12:33

Quote:

Originally Posted by piterw82 (Post 244088)
hello,

I am new in OpenFOAM and I am starting my work on LES multiphase with EE approach. I see you guys are already working on that topic and probably have a lot of experience in that field. I was wondering Alberto if it would be possible to have a look on your code? I read it is not checked yet but I am sure it would help me to understand how I can make some changes inside the source terms in OpenFOAM. Thank you in advance.

Peter

Hi Peter,

the main problem is that the code is outdated. I have written it using OF 1.3/1.4, and I did not adapt it to the new versions of OpenFOAM(r).

However if you drop me an email I can send it to you.

P.S. Juho Peltola used my code as a starting point to write his own, and he used OF 1.5. Maybe you want to contact him too ;-)

Best,


All times are GMT -4. The time now is 05:45.