CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

nutSpalartAllmarasWallFunction

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   April 21, 2010, 06:56
Default nutSpalartAllmarasWallFunction
  #1
Member
 
Axel Söhngen
Join Date: Jan 2010
Location: Germany, Trier
Posts: 31
Rep Power: 7
chelvistero is on a distinguished road
Hello Foamers,
I'd like to simulate a near wall treatment with simpleFoam. In other threads I have read that there is the possibiltity to simulate KOmegaSST with near wall treatment using nutSpalartAllmarasWallFunction in the nut file. In the k -file I use kqRWallFunction and in the omega - file I use omegaWallFunction. But using this I always get a failure message. Does anybody know how to fix this?

Thank you in advance

Failure message:

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kOmegaSST
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/axel/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/axel/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/axel/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 Foam:perator/(Foam::tmp<Foam::Field<double> > const&, Foam::UList<double> const&) in "/home/axel/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libOpenFOAM.so"
#5 Foam::incompressible::RASModels::nutSpalartAllmara sWallFunctionFvPatchScalarField::calcNut() const in "/home/axel/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libincompressibleRASModels.so"
#6 Foam::incompressible::RASModels::nutWallFunctionFv PatchScalarField::updateCoeffs() in "/home/axel/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libincompressibleRASModels.so"
#7 Foam::fvPatchField<double>::evaluate(Foam::Pstream ::commsTypes) in "/home/axel/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linuxGccDPOpt/simpleFoam"
#8 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::evaluate() in "/home/axel/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linuxGccDPOpt/simpleFoam"
#9 Foam::incompressible::RASModels::kOmegaSST::kOmega SST(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/home/axel/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libincompressibleRASModels.so"
#10 Foam::incompressible::RASModel::adddictionaryConst ructorToTable<Foam::incompressible::RASModels::kOm egaSST>::New(Foam::GeometricField<Foam::Vector<dou ble>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/home/axel/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libincompressibleRASModels.so"
#11 Foam::incompressible::RASModel::New(Foam::Geometri cField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/home/axel/OpenFOAM/OpenFOAM-1.6.x/lib/linuxGccDPOpt/libincompressibleRASModels.so"
#12 main in "/home/axel/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linuxGccDPOpt/simpleFoam"
#13 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#14 _start at /build/buildd/eglibc-2.10.1/csu/../sysdeps/i386/elf/start.S:122
Floating point exception
chelvistero is offline   Reply With Quote

Old   April 22, 2010, 04:24
Default
  #2
Member
 
Axel Söhngen
Join Date: Jan 2010
Location: Germany, Trier
Posts: 31
Rep Power: 7
chelvistero is on a distinguished road
Is there a simple possibility to simulate a turbolent flow without any wall functions? And what bc's are necessary to do this? I would be very happy about help.

Thank you in advance
chelvistero is offline   Reply With Quote

Old   April 22, 2010, 04:47
Default
  #3
Senior Member
 
Florian Krause
Join Date: Mar 2009
Location: Munich
Posts: 103
Rep Power: 8
florian_krause is on a distinguished road
Hi Axel,

you can simulate a turbulent flow without any wall functions, by DNS, LES or RANS with a low-Re turbulence model. From your previous post, I assume you want to perform RANS calculations, right?! OpenFOAM offers different low-Re RANS turbulence models. I only have experience with the LaunderSharmaKE model, for which I set:

epsilon at the wall:
zeroGradient

k at the wall:
type fixedValue;
value uniform 1e-10;

your first grid point normal to the wall should have y+ < 1 in order to resolve the wall!

good luck!
Florian
florian_krause is offline   Reply With Quote

Old   April 22, 2010, 04:54
Default
  #4
Member
 
Axel Söhngen
Join Date: Jan 2010
Location: Germany, Trier
Posts: 31
Rep Power: 7
chelvistero is on a distinguished road
Thank you for your fast reply. I will try it out immediatly!!!
chelvistero is offline   Reply With Quote

Old   April 22, 2010, 16:53
Default
  #5
Senior Member
 
BastiL
Join Date: Mar 2009
Posts: 471
Rep Power: 11
bastil is on a distinguished road
I ran into the same problems with nutSpalartAllmarasWallFUnction. I still hold the opinion it should work but I don't know how.
bastil is offline   Reply With Quote

Old   April 23, 2010, 03:55
Default
  #6
Member
 
Axel Söhngen
Join Date: Jan 2010
Location: Germany, Trier
Posts: 31
Rep Power: 7
chelvistero is on a distinguished road
Hello Bastil,
when you are using OpenFoam 1.6.x check the updates. I just updated and now it worked. But there is no yPlus output. Do you know if it is advisable to check the yPlus value first with nutWallFunction and after that to turn to nutSpalartAllmarasWallFunction?

Now I have try to simulate with nutSpalartAllmarasWallFunction but the results are not so good i hoped for. For k I used fixedValue 1e-10 and for omega i used 1.zeroGradient
and 2. omegaWallFunction uniform 0 . Has anybody experience with this kind of lowRe Simulation? And are the boundaries correct?
I would be very happy about answers and suggestions!!

Last edited by chelvistero; April 23, 2010 at 14:02.
chelvistero is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On



All times are GMT -4. The time now is 09:49.