[maalan]

Hi, Alberto!!

First of all thanks for the code! It is useful for me as I work in aerodynamics. My question is: up to the moment I used to running both k-epsilon and k-omega models, and its respectives wall functions for k, epsilon and omega. In addition, I knew how to set up the k and epsilon values in my 0 folder by using the turbulentIntensityKineticEnergyInlet and turbulentMixingLengthDissipationRateInlet boundary conditions. How should I proceed with the v2f model?? the tutorial that comes with the code just set the values for the 4 unknowns...
hope you could help!

Best regards!
Antonio

[maalan]

Could anyone send to me the v2f code for OF-1.6??

Thanks in advance!

[lakeat]

Hi Antonio,

I'll reply you here. I am not using 1.6-ext right now. But you can post your error message though. Compilation is not that difficult.

But I want to re-state, that I have never had this v2f model work for me (alberto's version and my own version and many others')
So those who have made the code really work, would you please post your results on a simple channel flow?? Because I really think that there is somthing wrong in the code.

[maalan]

Ok!! Thank you for your prompt reply, lakeat!

The point is I have just the OF-1.6 and I am trying to compile the romif codes for such OpenFOAM release. I have set my Make/files and Make/options files as:

Make/files:
DurbinV2F.C
LIB = $(FOAM_USER_LIBBIN)/libmyIncompressibleRASModels

Make/options:
EXE_INC = \
-I$(LIB_SRC)/turbulenceModels \
-I$(LIB_SRC)/transportModels \
-I$(LIB_SRC)/finiteVolume/lnInclude \
-I$(LIB_SRC)/meshTools/lnInclude \
-I$(LIB_SRC)/thermophysicalModels/specie/lnInclude \
-I$(LIB_SRC)/thermophysicalModels/reactionThermo/lnInclude \
-I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
-I$(LIB_SRC)/turbulenceModels/incompressible/RAS/lnInclude

LIB_LIBS = \
-lincompressibleTurbulenceModel \
-lfiniteVolume \
-lmeshTools \
-lcompressibleRASModels \
-lreactionThermophysicalModels \
-lspecie \
-lbasicThermophysicalModels \
-lgenericPatchFields

After doing wmake, the error which I get is (at the end of the compilation):

Make/linux64GccDPOpt/DurbinV2F.o: In function `Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::magSqr<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<Foam ::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> > const&)':
DurbinV2F.C.text._ZN4Foam6magSqrINS_10SymmTensor IdEENS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14G eometricFieldIdT0_T1_EEEERKNS5_INS6_IT_S7_S8_EEEE[Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::magSqr<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<Foam ::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> > const&)]+0x144): undefined reference to `Foam::calculatedFvPatchField<double>::typeName'
Make/linux64GccDPOpt/DurbinV2F.o.data.rel.ro._ZTVN4Foam14incompressib le9RASModels9DurbinV2FE[vtable for Foam::incompressible::RASModels:urbinV2F]+0x70): undefined reference to `Foam::incompressible::RASModel:rintCoeffs()'
Make/linux64GccDPOpt/DurbinV2F.o.data.rel.ro._ZTVN4Foam14incompressib le9RASModels9DurbinV2FE[vtable for Foam::incompressible::RASModels:urbinV2F]+0x78): undefined reference to `Foam::incompressible::RASModel::yPlus(int, double) const'
Make/linux64GccDPOpt/DurbinV2F.o.data.rel.ro._ZTIN4Foam14incompressib le9RASModels9DurbinV2FE[typeinfo for Foam::incompressible::RASModels:urbinV2F]+0x10): undefined reference to `typeinfo for Foam::incompressible::RASModel'
Make/linux64GccDPOpt/DurbinV2F.o.data.rel.ro._ZTIN4Foam11wallFvPatchE[typeinfo for Foam::wallFvPatch]+0x10): undefined reference to `typeinfo for Foam::fvPatch'
collect2: ld returned 1 exit status
make: *** [OpenFOAM.out] Error 1

Thank you for your help!
Best,
Antonio

[1/153]

Okay, these cannot be fixed in one reply, Would you mind send you v2f code to me, I'll compile here and send it back to you.

DANIEL_DOT_WEI_AT_ND_DOT_EDU

[1/153]

It works here, not sure why you cant make it. I suggest you under
/turbulenceModels/incompressible/RAS folder,
do a "wclean", and "wmake libso" again

[maalan]

Ok, I have done it: wclean and then wmake libso in the RAS folder, and the the same inside the V2F folder, and I got the next error:

[antonio@gauss0 V2F]$ wmake libso
wmakeLnInclude: linking include files to ./lnInclude
Making dependency list for source file DurbinV2F.C
SOURCE=DurbinV2F.C ; g++ -m64 -Dlinux64 -DWM_DP -Wall -Wno-strict-aliasing -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3 -DNoRepository -ftemplate-depth-40 -I/home/antonio/OpenFOAM/OpenFOAM-1.6/src/turbulenceModels -I/home/antonio/OpenFOAM/OpenFOAM-1.6/src/transportModels -I/home/antonio/OpenFOAM/OpenFOAM-1.6/src/finiteVolume/lnInclude -I/home/antonio/OpenFOAM/OpenFOAM-1.6/src/meshTools/lnInclude -I/home/antonio/OpenFOAM/OpenFOAM-1.6/src/thermophysicalModels/specie/lnInclude -I/home/antonio/OpenFOAM/OpenFOAM-1.6/src/thermophysicalModels/reactionThermo/lnInclude -I/home/antonio/OpenFOAM/OpenFOAM-1.6/src/thermophysicalModels/basic/lnInclude -I/home/antonio/OpenFOAM/OpenFOAM-1.6/src/turbulenceModels/incompressible/RAS/lnInclude -IlnInclude -I. -I/home/antonio/OpenFOAM/OpenFOAM-1.6/src/OpenFOAM/lnInclude -I/home/antonio/OpenFOAM/OpenFOAM-1.6/src/OSspecific/POSIX/lnInclude -fPIC -c $SOURCE -o Make/linux64GccDPOpt/DurbinV2F.o
/usr/bin/ld: cannot find -lgenericPatchFields
collect2: ld returned 1 exit status
make: *** [/home/antonio/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libmyIncompressibleRASModels.so] Error 1

Maybe this one is easier to understand!

[1/153]

I dont know where did you get this Make/options file

Make/options:
EXE_INC = \
-I$(LIB_SRC)/turbulenceModels \
-I$(LIB_SRC)/transportModels \
-I$(LIB_SRC)/finiteVolume/lnInclude \
-I$(LIB_SRC)/meshTools/lnInclude \
-I$(LIB_SRC)/thermophysicalModels/specie/lnInclude \
-I$(LIB_SRC)/thermophysicalModels/reactionThermo/lnInclude \
-I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
-I$(LIB_SRC)/turbulenceModels/incompressible/RAS/lnInclude

LIB_LIBS = \
-lincompressibleTurbulenceModel \
-lfiniteVolume \
-lmeshTools \
-lcompressibleRASModels \
-lreactionThermophysicalModels \
-lspecie \
-lbasicThermophysicalModels \
-lgenericPatchFields

So, link to genericPatchFields is for Official OF version, so delete it and use the ext version's.

I might be wrong, because the ext version on my computer is not the latest one, just do a double check.

[maalan]

Ok! Thank you very much, daniel! I was trying to compile it in the cluster... Now, I will try the same with the OF-1.6 release in my PC later. Otherwise I will compile the OF-1.6-ext and I will try again. I will tell you in a short time.

[sandy]

Hi, I have no new version about V2-f for OF-1.6, and just explicitly give the coefficent values of DurbinV2's version in order to avoid a lookup file, you know. Attachment is it. Who need it?

[s.m]

hi all
is v2f a low re number turbulence model?
thank you

[Alhasan]

Hey,

From my knowledge and from the Durbin's papers i have read, it is for only high reynolds number since they are interested in heat transfer at stagnation point, the benifit from durbins model is it is very good at stagnation point for high reynolds number !!!
never come across a paper with low re and flow close to the wall for aerofoil or impinging-jet (the only stuff i know lol)
k-kl- is a good low Re i from my knowledge

and still i am a student of life so senior member please correct me if i am wrong
and btw the above model attached is a good v2f model but for some reason i am not fully able to validate all experiments with it, im cent percent sure it needs some corrections

Best Regards,
Hasan K.J

[makayasa]

Hi all,
after running paraFoam, in the terminal to get a message like this

From function dlOpen(const fileName&, const bool)
in file POSIX.C at line 1179
dlopen error : /home/yasa/OpenFOAM/yasa-2.3.0/platforms/linuxGccDPOpt/lib/libmyIncompressibleRASModels.so: undefined symbol: _ZTIN4Foam14incompressible8RASModelE

FOAM Warning :
From function dlLibraryTable:pen(const fileName&, const bool)
in file db/dynamicLibrary/dlLibraryTable/dlLibraryTable.C at line 99
could not load "libmyIncompressibleRASModels.so"

How to fix it? Please help, this is my final assignment

[hamide]

hello friends
i'm new in openfoam
i wanna model turbulence flow with v2f model. in the tutorial there isn't any example for this model. i copied one of k-epsilon tutorials and changed it's B.C to what i want. but i faced to an error!!!!
should i add any thing to fvSolution file? is it necessary to change other things?

the error is shown below:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.2.0-5be49240882f
Exec : pisoFoam
Date : Jan 13 2016
Time : 18:31:56
Host : "hamide-HP-ProBook-4530s"
PID : 4096
Case : /home/hamide/run/tutorials/incompressible/pisoFoam/ras/aaaaV2f
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model v2f
bounding v2, min: 0 max: 0.0625 average: 0.0625
v2fCoeffs
{
Cmu 0.22;
CmuKEps 0.09;
C1 1.4;
C2 0.3;
CL 0.23;
Ceta 70;
Ceps2 1.9;
sigmaK 1;
sigmaEps 1.3;
}


Starting time loop

Time = 5e-05

Courant Number mean: 0.0003125 max: 0.529571
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 7.05865e-06, No Iterations 7
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 9.57628e-06, No Iterations 7
DICPCG: Solving for p, Initial residual = 1, Final residual = 0.098608, No Iterations 539
time step continuity errors : sum local = 6.163e-05, global = -8.17162e-07, cumulative = -8.17162e-07
DICPCG: Solving for p, Initial residual = 0.0235598, Final residual = 9.65185e-07, No Iterations 582
time step continuity errors : sum local = 1.48876e-06, global = -1.1141e-10, cumulative = -8.17274e-07
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#6 Foam::incompressible::RASModels::v2f::correct() in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#7
in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/pisoFoam"
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9
in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/pisoFoam"
Floating point exception (core dumped)



would you please help me?

[dneves]

Hi everyone

I am a little puzzled because the https://www.cfd-online.com/Wiki/V2-f_models

sais that v2-f turbulence model should not be applied to solve eulerian multifase problems.

Why is that?...?

I am testing a multifase problem (sea waves breaking) and it works fine...I don't understant why should I not use it.

Can anyone enlight me?

Best regards