[alekhine]

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:

--> 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]

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

--> 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]

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]

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