[kaaja]

Hi!

I am trying to incorporate a stationary external inlet BC for U in a channel flow problem in openFoam using timeVaryingMappedFixedValue.

I have managed to solve the problem with a parabolic inlet BC, but when I change to the external BC using timeVaryingMappedFixedValue, which looks a little bit more flat than the parabolic profile, the program will not solve. All velocity values for the exernal BC are in the same ballpark as for the parabolic velocity profile, which solved, so the physics should not be much changed.

My files look very much like the original PitzDaily from the tutorial. What I have changed is:
- Copied the folder boundaryData from pitzDailyExptInlet-tutorial
- Deleted the files k and u in the above mentioned folder, since I do not change the inlet BC's for these variables (also tried keeping these files, but that does not change anything)
- Edited the files in "U" and "points" in boundaryData to my physical problem
- In the file "U" in the normal folder "0", I changed the U-inlet to be as in the corresponding file in
pitzDailyExptInlet.

Doing the above gives the error message below. I hope someone maybe see what is going wrong here!

I tried doing the same changes in the original PitzDaily, and that worked.

The error message I get for my channel problem:

Build  : 4.1
Exec   : simpleFoam

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


SIMPLE: convergence criteria
    field p     tolerance 0.01
    field U     tolerance 0.001
    field "(k|epsilon|omega|f|v2)"     tolerance 0.001

Reading field p

Reading field U

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::pointToPointPlanarInterpolation::calcCoordinateSystem(Foam::Field<Foam::Vector<double> > const&) const at ??:?
#4  Foam::pointToPointPlanarInterpolation::pointToPointPlanarInterpolation(Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, double, bool) at ??:?
#5  Foam::timeVaryingMappedFixedValueFvPatchField<Foam::Vector<double> >::checkTable() at ??:?
#6  Foam::timeVaryingMappedFixedValueFvPatchField<Foam::Vector<double> >::updateCoeffs() at ??:?
#7  Foam::timeVaryingMappedFixedValueFvPatchField<Foam::Vector<double> >::timeVaryingMappedFixedValueFvPatchField(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#8  Foam::fvPatchField<Foam::Vector<double> >::adddictionaryConstructorToTable<Foam::timeVaryingMappedFixedValueFvPatchField<Foam::Vector<double> > >::New(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#9  Foam::fvPatchField<Foam::Vector<double> >::New(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#10  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::Boundary::readField(Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#11  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readFields(Foam::dictionary const&) at ??:?
#12  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readFields() at ??:?
#13  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&) at ??:?
#14  ? at ??:?
#15  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#16  ? at ??:?
Floating point exception (core dumped)

[kaaja]

I think I found one bug in my program.
For "points" I used the exact same points for "minz" and "maxz".
Looking at what was done in PitzDaily, there was a slight difference in the z-coordinate.
Now I have changed the z-coordinate slightly, but still it will not run.
The new error message is:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


SIMPLE: convergence criteria
    field p     tolerance 0.01
    field U     tolerance 0.001
    field "(k|epsilon|omega|f|v2)"     tolerance 0.001

Reading field p

Reading field U



--> FOAM FATAL IO ERROR: 
"ill defined primitiveEntry starting at keyword 'version' on line 10 and ending at line 98"

file: /home/.../constant/boundaryData/inlet/0/U at line 98.

    From function void Foam::primitiveEntry::readEntry(const Foam::dictionary&, Foam::Istream&)
    in file db/dictionary/primitiveEntry/primitiveEntryIO.C at line 189.

FOAM exiting

[kaaja]

Found the error.
One file included a parenthesis that was not supposed to be there.