[sita]

Hi all,


I just ran a case using simpleFoam, and now I'd like to compute the wall shear stress in my geometry. I tried this using

simpleFoam -postProcess -func wallShearStress

which seems to work fine for time 0, but when it starts on the next time step, I get a fatal error saying

request for volVectorField wallShearStress from objectRegistry region0 failed
available objects of type volVectorField are
1(U)

It looks like OpenFOAM is telling me it can't find the wallShearStress field, but that's exactly what it's supposed to be computing, so no wonder it isn't there yet, right? Or not? Does anyone know what's going on here, and how I should avoid this error?


Thanks in advance!
Sita


P.S. I'm using OpenFOAM 6 on Ubuntu 18.04


Edit:

-other functions seem to work fine, I tried postProcess -func "components(U)" and that didn't give any problems
- when I try using the function wallShearStress under Foam Extend 4.0 I get the following error:

--> FOAM FATAL IO ERROR: 

    Cannot find 'value' entry on patch symmetry of field U in file "/home/administrator/foam/administrator-4.0/run/cases/of6/FJ26C_turbulent/0/U"
    which is required to set the values of the generic patch field.
    (Actual type symmetry)

    Please add the 'value' entry to the write function of the user-defined boundary-condition
    or link the boundary-condition into libfoamUtil.so

file: /home/administrator/foam/administrator-4.0/run/cases/of6/FJ26C_turbulent/0/U::boundaryField::symmetry from line 4514636 to line 4514636.

    From function genericFvPatchField<Type>::genericFvPatchField(const fvPatch&, const Field<Type>&, const dictionary&)
    in file fields/fvPatchFields/basic/generic/genericFvPatchField.C at line 71.

FOAM exiting

So apparently something's wrong with my symmetry boundary patch. I already noticed before that I seem to have two of these in the same location. Does anyone know how to fix that? I used snappyHexMesh to create my mesh


Edit 2: I fixed the mesh now, but still get the same fatal error in OF6 upon running

simpleFoam -postProcess -func wallShearStress

--> FOAM FATAL ERROR: 

    request for volVectorField wallShearStress from objectRegistry region0 failed
    available objects of type volVectorField are
1(U)

    From function const Type& Foam::objectRegistry::lookupObject(const Foam::word&) const [with Type = Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 193.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const& Foam::objectRegistry::lookupObject<Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> >(Foam::word const&) const in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
#3  Foam::functionObjects::wallShearStress::execute() at ??:?
#4  Foam::functionObjects::timeControl::execute() at ??:?
#5  Foam::functionObjectList::execute() at ??:?
#6  ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
#7  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8  ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/simpleFoam"
Aborted (core dumped)

[sita]

Alright, apparently this is some sort of bug: when I use the -time option and process the time folders one by one, it runs without trouble. So at least I have the results I wanted now

[schuyler]

Did you submit a bug report for this? I have had the exact same problem.

[wyldckat]

Originally Posted by schuyler

Did you submit a bug report for this? I have had the exact same problem.

Quick question: Which OpenFOAM version/fork are you using?

[sita]

@Wyldcat: as I mentioned in my first post in this thread, I'm using OpenFOAM 6 on Ubuntu 18.04

Originally Posted by schuyler

Did you submit a bug report for this? I have had the exact same problem.

No, I didn't. Good idea, I will do that.

[wyldckat]

@Sita: My apologies, I was in a hurry and only asked schuyler, because I wanted to look into both situations at the same time...

Regarding OpenFOAM 6, the bug is known and has already been fixed in the git version: https://github.com/OpenFOAM/OpenFOAM-6/commit/0f14683c - but the cumulative fixes haven't been release as Deb packages yet: https://bugs.openfoam.org/view.php?id=3055
Therefore, you will have to rebuild OpenFOAM 6 from source code, in order to fix this issue.

Or as an alternative, you can use the workaround detailed on this post: request for volScalarField k from objectRegistry region0 failed+(DPMFoam) - post #25

[sita]

Thanks, that's good to know!