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

Error or bug in post-processing mapped turbulent kinetic energy

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 6, 2017, 09:48
Default Error or bug in post-processing mapped turbulent kinetic energy
  #1
New Member
 
DimitriF
Join Date: Dec 2016
Location: London
Posts: 19
Rep Power: 9
alekhine is on a distinguished road
Hey guys,

I have a weird error I have never faced before. To get the initial conditions for my main simulation, I simulate a flat plate. Afterwards, I map this results to the main simulation using timeVaryingMappedFixedValue.

I test 2 RANS models, namely Spalart-Allmaras and kOmegaSST. Both simulations converge and the results are as expected (friction and pressure coefficient). Calculating the Reynolds stess using the implemented post-processing utility works fine for the SA model but I got the following error for kOmegaSST:

Code:
--> FOAM Warning : 
--> FOAM FATAL IO ERROR: 
Cannot find patchField entry for inlet

file: /home/dima/OpenFOAM/dima-4.0/run/turbProjectUlg/SST/v2-1/6000/turbulenceProperties:R.boundaryField from line 16030 to line 16039.

    From function void Foam::GeometricField<Type, PatchField, GeoMesh>::Boundary::readField(const Foam::DimensionedField<TypeR, GeoMesh>&, const Foam::dictionary&) [with Type = Foam::SymmTensor<double>; PatchField = Foam::fvPatchField; GeoMesh = Foam::volMesh]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/OpenFOAM/lnInclude/GeometricBoundaryField.C at line 191.
I don't know why but if I change the inlet condition of the kinetic energy (.0/k) to a random fixed value, it is possible to calculae the Reynolds stress. However, the result is obvisously wrong.

What do you think?

Dimitri
alekhine is offline   Reply With Quote

Old   January 6, 2017, 09:57
Default
  #2
New Member
 
DimitriF
Join Date: Dec 2016
Location: London
Posts: 19
Rep Power: 9
alekhine is on a distinguished road
Sorry, a short update.

The above error is displayed when I execute the postProcess command for the second time. This means there is already a Reynolds stress file but which is 'unfinished'.

When executed for the first time, this error comes up

Code:
--> FOAM FATAL ERROR: 
object of type N4Foam9Function1INS_10SymmTensorIdEEEE is not allocated

    From function T* Foam::autoPtr<T>::operator->() [with T = Foam::Function1<Foam::SymmTensor<double> >]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/OpenFOAM/lnInclude/autoPtrI.H at line 176.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  Foam::timeVaryingMappedFixedValueFvPatchField<Foam::SymmTensor<double> >::write(Foam::Ostream&) const at ??:?
#3  Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>::Boundary::writeEntry(Foam::word const&, Foam::Ostream&) const at ??:?
#4  Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>::writeData(Foam::Ostream&) const at ??:?
#5  Foam::regIOobject::writeObject(Foam::IOstream::streamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const at ??:?
#6  Foam::functionObjects::regionFunctionObject::writeObject(Foam::word const&) at ??:?
#7  Foam::functionObjects::turbulenceFields::write() at ??:?
#8  Foam::functionObjects::timeControl::write() at ??:?
#9  Foam::functionObjectList::execute() at ??:?
#10  ? at ??:?
#11  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#12  ? at ??:?
./GetReStress: line 2: 10618 Aborted                 (core dumped) simpleFoam -postProcess -func 'turbulenceFields(R)' -latestTime
I checked the sourced code and libraries but couldn't find anything useful.
alekhine is offline   Reply With Quote

Old   January 6, 2017, 18:08
Default
  #3
New Member
 
DimitriF
Join Date: Dec 2016
Location: London
Posts: 19
Rep Power: 9
alekhine is on a distinguished road
No one has a clue why this happens?
From function T* Foam::autoPtr<T>:perator->() [with T = Foam::Function1<Foam::SymmTensor<double> >] what is the tensor T here? I doubt it's the temperature since I use simpleFoam.

I have to submit a report in a few day and the Reynolds stress of the kOmegaSST model is the last result that is missing

Dimitri
alekhine is offline   Reply With Quote

Old   January 9, 2017, 09:38
Default
  #4
New Member
 
DimitriF
Join Date: Dec 2016
Location: London
Posts: 19
Rep Power: 9
alekhine is on a distinguished road
I don't understand how I searched before, and I search for a long time to find something related, but finally I found a solution to my problem.

lifeinhand suggests a way to calculate the Reynold stress when having issues with timeVaryingMappedFixedValue:
Error using R utility
alekhine is offline   Reply With Quote

Reply

Tags
bug, komegasst, mapped boundary condition, reynolds stress

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Turbulent kinetic energy in k-epsilon model khunyeu FLUENT 0 October 7, 2013 01:31
Contours of turbulent kinetic energy in LES using ANSYS FLUENT 13 khosrow FLUENT 5 July 29, 2013 05:28
CFX-POST : gradients for turbulent kinetic energy syler3321 CFX 1 June 16, 2011 05:47
sgs turbulent kinetic energy blons Main CFD Forum 0 February 25, 2010 12:16
Turbulent Kinetic Energy Olga FLUENT 2 October 11, 2002 15:05


All times are GMT -4. The time now is 06:34.