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

k epsilon boundary conditions for pisofoam ras solver

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 18, 2009, 06:11
Question k epsilon boundary conditions for pisofoam ras solver
  #1
New Member
 
Join Date: Dec 2009
Posts: 5
Rep Power: 7
chamoun is on a distinguished road
Hi,

I am struggling with understanding the proper boundary conditions for inlet, outlet and fixed walls for the epsilon, k, nut, nutilda and R files located in the 0 folder. Any feedback would be appreciated. I'm not sure when to use epsilonWallFunction, kqRWallFunction, and nutWallFunction. This is the error message that I'm getting:







----------------------------------------

Invalid wall function specification
Patch type for patch inlet must be wall
Current patch type is patch

#0 Foam::error:rintStack(Foam::Ostream&) in "/home/george/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/george/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 Foam::fvPatchField<double>::adddictionaryConstruct orToTable<Foam::incompressible::RASModels::kqRWall FunctionFvPatchField<double> >::New(Foam::fvPatch const&, Foam:imensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/home/george/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#3 Foam::fvPatchField<double>::New(Foam::fvPatch const&, Foam:imensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/home/george/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/pisoFoam"
#4 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::GeometricB oundaryField(Foam::fvBoundaryMesh const&, Foam:imensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/home/george/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/pisoFoam"
#5 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readField(Foam::Istream&) in "/home/george/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/pisoFoam"
#6 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&) in "/home/george/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/pisoFoam"
#7 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::incompressible::autoCreateWallFunctionField< double, Foam::incompressible::RASModels::kqRWallFunctionFv PatchField<double> >(Foam::word const&, Foam::fvMesh const&) in "/home/george/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#8 Foam::incompressible::autoCreateK(Foam::word const&, Foam::fvMesh const&) in "/home/george/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#9 Foam::incompressible::RASModels::kEpsilon::kEpsilo n(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/home/george/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#10 Foam::incompressible::RASModel::adddictionaryConst ructorToTable<Foam::incompressible::RASModels::kEp silon>::New(Foam::GeometricField<Foam::Vector<doub le>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/home/george/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/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/george/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#12 Foam::incompressible::turbulenceModel::addturbulen ceModelConstructorToTable<Foam::incompressible::RA SModel>::NewturbulenceModel(Foam::GeometricField<F oam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/home/george/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#13 Foam::incompressible::turbulenceModel::New(Foam::G eometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in "/home/george/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleTurbulenceModel.so"
#14 main in "/home/george/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/pisoFoam"
#15 __libc_start_main in "/lib64/libc.so.6"
#16 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116


From function kqRWallFunctionFvPatchField::checkType()
in file derivedFvPatchFields/wallFunctions/kqRWallFunctions/kqRWallFunction/kqRWallFunctionFvPatchField.C at line 48.

FOAM aborting
chamoun is offline   Reply With Quote

Old   December 18, 2009, 07:57
Default
  #2
Member
 
Julien Schaguene
Join Date: Apr 2009
Location: France
Posts: 55
Rep Power: 8
Schag is on a distinguished road
I'm no expert in turbulence, but I can quickly give you what I have...

each time a function is called 'variable'WallFunction, it can only be applied on a boundary which type is WALL (and no patch), what is mentionned in constant/polyMesh/boundary.

Unfortunatly, I cannot really help further.

Just my 2 cents, regards
Schag is offline   Reply With Quote

Old   December 18, 2009, 10:40
Default
  #3
New Member
 
Join Date: Apr 2009
Posts: 26
Rep Power: 8
askjak is on a distinguished road
All patches with a wall function need to specified as "wall" and not "patch" (see the very first error message). But you should not supply a wallfunction on the inlet.

For example my 0/epsilon looks like this
Quote:
inlet
{
type fixedValue;
value uniform <some value that you specify>;
}
Best regards, Ask
askjak is offline   Reply With Quote

Old   April 13, 2010, 14:05
Default
  #4
Senior Member
 
Nilesh Rane
Join Date: Apr 2010
Posts: 122
Rep Power: 7
nileshjrane is on a distinguished road
Quote:
Originally Posted by chamoun View Post
Hi,

I am struggling with understanding the proper boundary conditions for inlet, outlet and fixed walls for the epsilon, k, nut, nutilda and R files located in the 0 folder. Any feedback would be appreciated. I'm not sure when to use epsilonWallFunction, kqRWallFunction, and nutWallFunction. This is the error message that I'm getting:


The value of k, eps and nu goes to zero at wall. and the inlet valu depends upon the kind of case one is solving. Basically its specifying the incoming turbulence. There are other ways to specify the valu of k and eps as well, like turbulent intencity and length scale. I am not very much sure about the inlet values but wall conditions should be zero.
__________________
Imagination is more important than knowledge..
nileshjrane is offline   Reply With Quote

Old   May 10, 2011, 05:30
Default R Boundary Condition
  #5
New Member
 
Rogelio Chovet
Join Date: Feb 2011
Posts: 4
Rep Power: 6
chichovet is on a distinguished road
I have already specified k and epsilon conditions by solving the respectives equations. But I'm struggling with the R condition what value should I put for an turbulent developed flow inlet??
Thanks in advance
chichovet is offline   Reply With Quote

Reply

Tags
boundary conditions, k-epsilon, nutilda, pisofoam, ras

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Wind turbine simulation Saturn CFX 34 October 16, 2014 05:27
TwoPhaseEulerFoam and Boundary conditions raagh77 OpenFOAM Running, Solving & CFD 98 September 23, 2013 06:31
Boundary Conditions Thomas P. Abraham Main CFD Forum 20 July 7, 2013 05:05
No results for solid domain Gary Holland CFX 10 March 13, 2009 04:30
Boundary Conditions and solver for low Mach laminar air jet fbisetti OpenFOAM Running, Solving & CFD 1 September 13, 2006 12:49


All times are GMT -4. The time now is 19:08.