CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   V2f turbulence model (https://www.cfd-online.com/Forums/openfoam-solving/58440-v2f-turbulence-model.html)

qtian December 26, 2007 15:50

Hello, all, Do we have v2-f
 
Hello, all,

Do we have v2-f turbulence model implemented in OpenFoam? I checked turbulence models in src/turbulencemodel folder and there is no such model. However, I searched message board and some folks discussed the implementation of v2f model in 2005. Am I missing anything? My version is 1.4.

Thanks.

Quinn

msha December 27, 2007 07:08

Not really Foam doesn't hav
 
Not really

Foam doesn't have this model

michele January 7, 2008 06:01

Quinn, here you will find an
 
Quinn,
here you will find an implementation of v2f (Lien/Kalitzin a.k.a. code-friendly version) I made a couple of years ago.
http://www.cfd-online.com/OpenFOAM_D...hment_icon.gif DurbinV2F.tar.gz
It is not fully tested, so use at your own risk...
As a tutorial test (BC and configuration example), herebelow you find a configuration for the asymmetric diffuser, in which the v2f model is known to predict quite well separation/reattachment. No grid is included.
http://www.cfd-online.com/OpenFOAM_D...hment_icon.gif diffuser_V2F_noMesh.tar.gz
Herebelow a plot of velocity/streamlines
http://www.cfd-online.com/OpenFOAM_D...ges/1/6321.gif
A comparison with experimental friction coefficient is given below
http://www.cfd-online.com/OpenFOAM_D...ges/1/6322.gif
From the Cf plot you can see the separation/reattachment points on the bottom wall...
Separation occurs at x/H=5 for v2f and at about 7 experimentally.
Reattachement occurs at x/H=29 for v2f and at about 29 experimentally.

Regards,
Michele.

qtian January 7, 2008 11:17

hi, Rattin, It is great. Th
 
hi, Rattin,

It is great. Thanks for sharing this model with me. I am doing three dimensional separated flow simulation now and want to try all different kind of turbulence models to see how good is the prediction.

qtian January 23, 2008 17:40

Dear Rattin, I am reading
 
Dear Rattin,

I am reading the v2f model you sent me. would you mind to answer a few questions for me about this v2f model?

Is this v2f model from Davidson's v2f model? Can you please give me a reference? What is the coefficient yStarLim for? I did not see this parameter in the durbin's version?

Thanks

QT

michele January 29, 2008 05:43

Quinn, sorry for the late r
 
Quinn,
sorry for the late reply, I was abroad last week.
Here attached you will find a small reference for the model implemented:
http://www.cfd-online.com/OpenFOAM_D...s/mime_pdf.gif v2f_description.pdf
The models includes the Davidson modifications and a limiter switch for time (T) and length (L) scales (Kolmogorov scales).
In this way the limiter is always active if you set yStarLim to a large number (say 10e+10), is never active if you set yStarLim to zero and is active partially if you set to, say, 30. The last case is consistent with the validity of the theoretical Kolmogorov microscales in the near-wall region.
Hope this helps.
P.S. v2f requires to mesh the near-wall region up to y<sup>+</sup><1, like any low-Re turbulence model.

Regards,
Michele.

qtian January 29, 2008 10:22

Michele, Thanks you so much
 
Michele,

Thanks you so much for your help and suggestions. I have completed my simulation. The results from v2f are awesome. This is the best results after I tried a number of turbulence models.

michele January 30, 2008 06:22

Quinn, I'm glad that the info
 
Quinn,
I'm glad that the info were useful, and I hope you will make available to the forum the experience/results once completed the study.

Michele

sek February 20, 2008 14:21

Hi Michele, Can you please
 
Hi Michele,

Can you please tell me how you computed and plotted in OpenFOAM or any package the skin-friction coefficient along the wall for your diffuser case?

sek February 20, 2008 14:53

Hi Michele, Can you please
 
Hi Michele,

Can you please tell me how you computed and plotted in OpenFOAM or any package the skin-friction coefficient along the wall for your diffuser case?

braennstroem April 14, 2008 10:27

Hi Rattin, as you wrote, yo
 
Hi Rattin,

as you wrote, you implemented the Davidson modifications. I assume you mean Lars Davidson and you 'just' implemnted the first modification. Do you have an idea, how to implenent the second modification, i.e. using different turbulent viscosities for different compononts?

Greetings!
Fabian

anger April 30, 2008 08:17

Hi Michele, in the turbulen
 
Hi Michele,

in the turbulenceProperties file of your V2F- example, there is also an entry for the constants of a k- Omega model. Can the implementation of this model be found somewhere in the forum or on sourceforge?

Best regards,
-Thomas

michele April 30, 2008 09:14

Thomas, the k-w model availab
 
Thomas,
the k-w model available in the OpenFOAM release is the SST version.
You may also want to try the SAS modification (and several other models), available at
http://openfoamwiki.net/index.php/Turbulence_Modeling
for download.

The k-w properties you found inside the dictionary refer to an implementation I made a couple of years ago (at the time the k-w model was not yet implemented in the OpenFOAM distribution).
If you are interested, I may pack the k-w models I implemented (little testing has been made).
The particularities of the models are:
- Wilcox k-w with and without the transitional modification
- consistent treatment of the wall region ("Model-consistent universal wall-function" - Knopp)
- rough surfaces (whose equivalent sand-roughness may be specified on a patch-by-patch basis)

Regards,
Michele.

anger April 30, 2008 10:53

Michele, thanks for the hin
 
Michele,

thanks for the hint to the wiki.
Do you have any references for the transitional modification of the k-w model you mentioned?

Best regards,
-Thomas

michele April 30, 2008 11:52

Thomas, the transitional modi
 
Thomas,
the transitional modification is needed when you want to analyse laminar-to-turbulent boundary layer transition (turbulence models tipically predict transition at Reynolds numbers at least one order of magnitude too low...).
The transitional modifications are discussed in the Wilcox's book "Turbulence modeling for CFD".

Regards
Michele.

anger May 21, 2008 04:24

Hello Michele, I downloaded
 
Hello Michele,

I downloaded the k-w models and compiled them (after changing volTensorField to volSymmTensorField for R) but trying to run the k-w model gave an error:

Selecting incompressible transport model Newtonian
Selecting turbulence model kOmega


--> FOAM FATAL ERROR : LHS and RHS of + have different dimensions
dimensions : [0 0 -1 0 0 0 0] + [0 2 -3 0 0 0 0]
#0 Foam::error::printStack(Foam:: Ostream&) in "/home/anger/OpenFOAM/OpenFOAM-1.4.1-dev/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/anger/OpenFOAM/OpenFOAM-1.4.1-dev/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 Foam::operator+(Foam::dimensionSet const&, Foam::dimensionSet const&) in "/home/anger/OpenFOAM/OpenFOAM-1.4.1-dev/lib/linux64GccDPOpt/libOpenFOAM.so"
#3 Foam::tmp<foam::geometricfield<double,> > Foam::operator+<foam::fvpatchfield,>(Foam::Geometr icField<double,> const&, Foam::dimensioned<double> const&) in "/home/anger/OpenFOAM/OpenFOAM-1.4.1-dev/lib/linux64GccDPOpt/libincompressibleTu rbulenceModels.so"
#4 Foam::turbulenceModels::kOmega::kOmega(Foam::Geome tricField<foam::vector<double> , Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::transportModel&) in "/home/anger/OpenFOAM/OpenFOAM-1.4.1-dev/lib/linux64GccDPOpt/libincompressibleTu rbulenceModels.so"
#5 Foam::turbulenceModel::adddictionaryConstructorToT able<foam::turbulencemodels::k omega>::New(Foam::GeometricField<foam::vector<doub le>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::transportModel&) in "/home/anger/OpenFOAM/OpenFOAM-1.4.1-dev/lib/linux64GccDPOpt/libincompressibleTu rbulenceModels.so"
#6 Foam::turbulenceModel::New(Foam::GeometricField<fo am::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::transportModel&) in "/home/anger/OpenFOAM/OpenFOAM-1.4.1-dev/lib/linux64GccDPOpt/libincompressibleTu rbulenceModels.so"
#7 main in "/home/anger/OpenFOAM/OpenFOAM-1.4.1-dev/applications/bin/linux64GccDPOpt/simple Foam"
#8 __libc_start_main in "/lib/libc.so.6"
#9 Foam::regIOobject::readIfModified() in "/home/anger/OpenFOAM/OpenFOAM-1.4.1-dev/applications/bin/linux64GccDPOpt/simple Foam"


From function operator+(const dimensionSet& ds1, const dimensionSet& ds2)
in file dimensionSet/dimensionSet.C at line 385.

FOAM aborting (FOAM_ABORT set)

It seems that there gets something messed up with w and epsilon, but I have no idea where this could come from. Did you experience similar behaviour in one of your implementations and could you solve it?

Best regards,
-Thomas

michele May 21, 2008 06:34

Thomas, no problems on my imp
 
Thomas,
no problems on my implementations.

I suspect there may be an error in the source code.
Enstropy has frequency units.
I never tried this SAS implementation.
In the laminarOmega.H code seems to lie the error:

dimensionedScalar
(
"omega", sqr(U_.dimensions())/dimTime, 0.0
)

At a first glance you may change the above with:
dimensionedScalar
(
"omega", dimless/dimTime, 0.0
)

Now I have no possibility to try it... so may you try this for me and let me know.

Regards,
Michele.

anger June 10, 2008 02:32

Hello Michele, your above m
 
Hello Michele,

your above mentioned fix did not do the trick. I assume that you did the turbulence models implementation for an older version than 1.4.1-dev. I decided to take the SST- model provided with OF and modify it using your implementation as template to program the standard k-w model which seemed to work.
I was now trying to implement the k-w model with transitional behaviour but ran into problems with some of the terms. You mentioned that you have implemeted this. Would you be so kind to share this implementation so that I can again use it as template?

Best regards,
-Thomas

braennstroem August 10, 2008 07:09

Hi Michele, I am just tryin
 
Hi Michele,

I am just trying to implement the zeta-f model, which is based on v2f by zeta=v2/k. For this I read in the thesis of Popovac http://repository.tudelft.nl/file/354721/370010, that the boundary condition for f should be similar implemented as for epsilon. I wonder, if there is a reason, tha you did not use it!? Was it more unstable?


Regards!
Fabian

braennstroem August 10, 2008 07:28

... one more question: for wh
 
... one more question:
for what, does one need pos(yStarLim-yStar_) for the calculation of T and L?

Fabian

braennstroem October 22, 2008 10:35

Hi, I am trying to convert
 
Hi,

I am trying to convert the provide v2f model to the new 1.5 version;
until now with not to much success:

http://www.cfd-online.com/OpenFOAM_D...hment_icon.gif v2f.tar.gz

lnInclude/RASModel.H: In static member function 'static Foam::autoPtr<foam::incompressible::rasmodel> Foam::incompressible::RASModel::adddictionaryConst ructorToTable<rasmodeltype>::New(const Foam::volVectorField&, const Foam::surfaceScalarField&, Foam::transportModel&) [with RASModelType = Foam::incompressible::RASModels::DurbinV2F]':
lnInclude/RASModel.H:140: instantiated from 'Foam::incompressible::RASModel::adddictionaryCons tructorToTable<rasmodeltype>::adddictionaryConstru ctorToTable(const Foam::word&) [with RASModelType = Foam::incompressible::RASModels::DurbinV2F]'
DurbinV2F/DurbinV2F.C:44: instantiated from here
lnInclude/RASModel.H:129: error: cannot allocate an object of abstract type 'Foam::incompressible::RASModels::DurbinV2F'
DurbinV2F/DurbinV2F.H:58: note: because the following virtual functions are pure within 'Foam::incompressible::RASModels::DurbinV2F':
lnInclude/RASModel.H:269: note: virtual Foam::tmp<foam::geometricfield<foam::symmtensor<do uble>, Foam::fvPatchField, Foam::volMesh> > Foam::incompressible::RASModel::devReff() const
lnInclude/RASModel.H:272: note: virtual Foam::tmp<foam::fvmatrix<foam::vector<double> > > Foam::incompressible::RASModel::divDevReff(Foam::v olVectorField&) const
make: *** [Make/linux64GccDPOpt/DurbinV2F.o] Error 1

Does anyone know, what my mistake is?

Fabian

niklas October 22, 2008 10:57

could be because you havent de
 
could be because you havent defined the virtuals
devReff and divDevReff

braennstroem October 22, 2008 11:14

Hi Niklas, thanks, yes I fo
 
Hi Niklas,

thanks, yes I forgot those virtuals... stupid, it's actually in the message.

http://www.cfd-online.com/OpenFOAM_D...hment_icon.gif v2f.tar.gz

Fabian

ivan_cozza May 29, 2009 10:22

Hi Fabian,
did you succeed to code the v2f model in 1.5 version?
Have any test case to learn how to set up it?
Thank you...

idrama March 11, 2010 13:55

Has anybody a v2f version for OpenFOAM-1.6

sandy August 1, 2010 11:30

Quote:

Originally Posted by idrama (Post 249585)
Has anybody a v2f version for OpenFOAM-1.6

Hi idrama, did you get the v2f version for OF-1.6? Could you send a copy to me?

In addition, according to its "Notes" in the website: http://www.cfd-online.com/Wiki/V2-f_models , can this model be used to solve the Phase-change problem?


sandy August 3, 2010 04:18

I got the version ...

idrama August 5, 2010 02:44

Did you get the source? If so, could you send it to me or send me the link?

Cheers

bb_ August 23, 2010 08:09

Hello,

could anybody please tell me how to install the V2F model for OpenFoam-1.5-dev? I downloaded the sources using the above link but now I don't actually know how to proceed...

Thanks!

jitendra September 14, 2010 03:47

v2f model
 
can u tell me what is the exact meaning of v2f turbulence model.
What are its application, presently wahat is the international status of this, and what are the importance of this.Reply me soon .

thanking you

qtian316 April 12, 2011 15:15

Hello, all,

Does anyone has v2-f turbulence model implemented in OpenFoam 1.7.1? I had source code for 1.4 and having trouble to compile it. Can anyone share it with me? Thanks a lot.

Quinn

Alhasan January 6, 2012 15:06

hey,sandy
can you please send the link from where you got the source for v2f model, i wanna get it in my v2.0 :)
thank you

lakeat January 18, 2012 19:45

Hi,

Two questions:

1. Has anyone written the zeta-f code?
2. How to set the wall b.c condition for v2 and f?

Alhasan February 15, 2012 12:26

they are downloading as a unk. file... how do you even open them...

lakeat February 15, 2012 12:34

You can have whatever extension you want, whatever name you want, it doesn't change any contents in the file.

In a word, just use any tar software to untar it as usual, that's it.

Alhasan February 15, 2012 12:36

I'm having troubles even opening it... :(... on top of that i need to compile a v2f model for open foam v.2.0 ...... :(

Alhasan February 15, 2012 12:38

hey sorry i managed to open it... .. now any good ideas for compiling it in open foam v2.0

lakeat February 15, 2012 12:40

Quote:

Originally Posted by Alhasan (Post 344625)
hey sorry i managed to open it... .. now any good ideas for compiling it in open foam v2.0

As usual.
Read the manual before ask, please

Alhasan February 15, 2012 13:52

parallels@ubuntu:~/Desktop/DurbinV2F$ wmake libso
linux64GccDPOpt/options:10: *** missing separator. Stop.
wmake error: file 'Make/linux64GccDPOpt/objectFiles' could not be created in /home/parallels/Desktop/DurbinV2F

Alhasan February 15, 2012 13:53

any idea... :)


All times are GMT -4. The time now is 01:58.