turbFoam
 

Hello!

I am using turbFoam to do a transient analysis. First I made a steady state analysis with simpleFoam but the results weren't good. So i decided to choose turbFoam , i changed the fvSchemes and fvSolution file but when i run , it stops :
Exec : turbFoam -parallel
Date : Jun 24 2010
Time : 10:13:20
Host : Foam1
PID : 1883
Case : /media/OpenFoam/Travaux/p-habitacle/foamProMesh-pisoFoam/pisoFoam
nProcs : 30
Slaves :
29
(
Foam1.21884
Foam1.21885
....
Foam8.20214
)

Pstream initialized with:
floatTransfer : 0
nProcsSimpleSum : 0
commsType : nonBlocking

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
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 RAS turbulence model realizableKE

Starting time loop

FunctionObject:fieldProcess:fieldProcess

FunctionObject:volumeReport:volStats
Logging volume statistics to file: "/media/OpenFoam/Travaux/p-habitacle/foamProMesh-pisoFoam/pisoFoam/processor0/../log/volStats_volumeStatistics.0"

Time = 1

Courant Number mean: 4.93168701749 max: 62294.1691783
[2]
[2]
[2]
request for surfaceScalarField UBlendingFactor from objectRegistry region0 failed
available objects of type surfaceScalarField are

3
(
weightingFactors
differenceFactors_
phi
)
#0 Foam::error::printStack(Foam::Ostream&) in "/media/OpenFoam/FOAMpro/FOAMpro-1.5-2.2/FOAM-1.5-2.2/lib/linux64GccDPOpt/libFOAM.so"
#1 Foam::error::abort() in "/media/OpenFoam/FOAMpro/FOAMpro-1.5-2.2/FOAM-1.5-2.2/lib/linux64GccDPOpt/libFOAM.so"
#2 Foam::Ostream& Foam::operator<< <Foam::error>(Foam::Ostream&, Foam::errorManip<Foam::error>) in "/media/OpenFoam/FOAMpro/FOAMpro-1.5-2.2/FOAM-1.5-2.2/applications/bin/linux64GccDPOpt/turbFoam"
#3 Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const& Foam::objectRegistry::lookupObject<Foam::Geometric Field<double, Foam::fvsPatchField, Foam::surfaceMesh> const>(Foam::word const&) const in "/media/OpenFoam/FOAMpro/FOAMpro-1.5-2.2/FOAM-1.5-2.2/lib/linux64GccDPOpt/libfiniteVolume.so"
#4 Foam::localBlended<Foam::Vector<double> >::weights(Foam::GeometricField<Foam::Vector<doubl e>, Foam::fvPatchField, Foam::volMesh> const&) const in "/media/OpenFoam/FOAMpro/FOAMpro-1.5-2.2/FOAM-1.5-2.2/lib/linux64GccDPOpt/libfiniteVolume.so"
#5 Foam::fv::gaussConvectionScheme<Foam::Vector<doubl e> >::fvmDiv(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const in "/media/OpenFoam/FOAMpro/FOAMpro-1.5-2.2/FOAM-1.5-2.2/lib/linux64GccDPOpt/libfiniteVolume.so"
#6 Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > Foam::fvm::div<Foam::Vector<double> >(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&, Foam::word const&) in "/media/OpenFoam/FOAMpro/FOAMpro-1.5-2.2/FOAM-1.5-2.2/applications/bin/linux64GccDPOpt/turbFoam"
#7 Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > Foam::fvm::div<Foam::Vector<double> >(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) in "/media/OpenFoam/FOAMpro/FOAMpro-1.5-2.2/FOAM-1.5-2.2/applications/bin/linux64GccDPOpt/turbFoam"
#8 main in "/media/OpenFoam/FOAMpro/FOAMpro-1.5-2.2/FOAM-1.5-2.2/applications/bin/linux64GccDPOpt/turbFoam"
#9 __libc_start_main in "/lib64/libc.so.6"
#10 __gxx_personality_v0 at /usr/src/packages/BUILD/glibc-2.3/csu/../sysdeps/x86_64/elf/start.S:116
[2]
[2]
[2] From function objectRegistry::lookupObject<Type>(const word&) const
[2] in file /users/echidna01/jacques/IconReleases//FOAMpro-1.5-2.2/FOAM-1.5-2.2/src/FOAM/lnInclude/objectRegistryTemplates.C at line 142.
[2]
FOAM parallel run aborting
[2]



This happens even when I lower the time step to a small value. Am I perhaps doing something wrong when going to transient? I simply set the start time, in controlDict, as the number of the latest iteration from the steady state case in whose folder all variables are saved. Will my boundary conditions be preserved now? Is this the right way to do it or am i missing something here....or somewhere else...

Any ideas appreciated!!

Are you using the localBlended scheme somewhere in your fvSchemes file?

If you are using that scheme, you have to implement a surfaceScalarField called UBlendingFactor in turbFoam, so that the localBlended scheme does know, which of the two provided schemes to use at which point of your computational domain. AFAIR, the values of UBlendingFactor should be somewhere between 1 and 0.

I hope this helps,
Jens