[ldxfx]

Dear OFers,
I am working on the subject of periodic flows in a pipe. I set the inlet boundary type "mappedPatch" and outlet "patch". It works when the gridding is excessive, but if I reduce number of gridding, errors appear as follows.
I don't know where is the mistake.Oh,really hope you can bail me out!
Thanks!

Codes:
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 laminar
No finite volume options present


SIMPLE: convergence criteria
field p tolerance 0.01
field U tolerance 0.001


Starting time loop

Time = 0.001

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigSegv::sigHandler(int) at ??:?
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::mapDistribute::mapDistribute(Foam::List<int> const&, Foam::List<int> const&) at ??:?
#4 Foam::mappedPatchBase::calcMapping() const at ??:?
#5 void Foam::mappedPatchBase::distribute<Foam::Vector<dou ble> >(Foam::List<Foam::Vector<double> >&) const at ??:?
#6 Foam::mappedPatchFieldBase<Foam::Vector<double> >::mappedField() const at ??:?
#7 Foam::mappedFixedValueFvPatchField<Foam::Vector<do uble> >::updateCoeffs() at ??:?
#8 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::updateCoef fs() at ??:?
#9 Foam::fvMatrix<Foam::Vector<double> >::fvMatrix(Foam::GeometricField<Foam::Vector<doub le>, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) at ??:?
#10
at ??:?
#11
at ??:?
#12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13
at ??:?

[coccone]

Hi,

I guess I have a similar error output, I am working with chtMultiRegionFoam solver on coupled solid elements. I am not experienced in OF so I will not be too helpful, but hopefully sharing info could enlighten you or someone else.

What do you mean by gridding? Mesh fineness?

Anyway this is the log output:

[1] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[1] #1  Foam::sigSegv::sigHandler(int) at ??:?
[1] #2   in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3  Foam::mapDistribute::mapDistribute(Foam::List<int> const&, Foam::List<int> const&) at ??:?
[1] #4  Foam::mappedPatchBase::calcMapping() const at ??:?
[1] #5  void Foam::mappedPatchBase::distribute<double>(Foam::List<double>&) const at ??:?
[1] #6  Foam::compressible::turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::updateCoeffs() at ??:?
[1] #7  Foam::mixedFvPatchField<double>::evaluate(Foam::UPstream::commsTypes) at ??:?
[1] #8  Foam::mixedEnergyFvPatchScalarField::updateCoeffs() at ??:?
[1] #9  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::updateCoeffs() at ??:?
[1] #10  Foam::fvMatrix<double>::fvMatrix(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) at ??:?
[1] #11  
[1]  at ??:?
[1] #12  
[1]  at ??:?
[1] #13  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #14  
[1]  at ??:?

A few things I noticed on my case:

- I have this error output if I run it in parallel only.

-On a single processor it runs without crashing but the heat-flux calculated with snGrad() on one side of a coupled patch is not consistent with the opposite neither in magnitude nor in sign. In the same case there are other coupled patches which eventually tend to converge though...

-This happens only if I include one particular region on its coupled boundary, otherwise no issues with single or parallel run.

- I am not sure whether it could be mesh-dependent since I have tried several settings and geometries with no success. I must admit from checkMesh my mesh is quite bad though (8500 underDeterminedCells and 40 lowWeightFaces over 1M cells, only 300 underDeterminedCells come from that last region I want to add)

Have fun!


PS: maybe posting in the solving section could have been more appropriate than in the pre-processing one

[wanghongjie]

Originally Posted by ldxfx

Dear OFers,
I am working on the subject of periodic flows in a pipe. I set the inlet boundary type "mappedPatch" and outlet "patch". It works when the gridding is excessive, but if I reduce number of gridding, errors appear as follows.
I don't know where is the mistake.Oh,really hope you can bail me out!
Thanks!

Codes:
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 laminar
No finite volume options present


SIMPLE: convergence criteria
field p tolerance 0.01
field U tolerance 0.001


Starting time loop

Time = 0.001

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigSegv::sigHandler(int) at ??:?
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::mapDistribute::mapDistribute(Foam::List<int> const&, Foam::List<int> const&) at ??:?
#4 Foam::mappedPatchBase::calcMapping() const at ??:?
#5 void Foam::mappedPatchBase::distribute<Foam::Vector<dou ble> >(Foam::List<Foam::Vector<double> >&) const at ??:?
#6 Foam::mappedPatchFieldBase<Foam::Vector<double> >::mappedField() const at ??:?
#7 Foam::mappedFixedValueFvPatchField<Foam::Vector<do uble> >::updateCoeffs() at ??:?
#8 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::updateCoef fs() at ??:?
#9 Foam::fvMatrix<Foam::Vector<double> >::fvMatrix(Foam::GeometricField<Foam::Vector<doub le>, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) at ??:?
#10
at ??:?
#11
at ??:?
#12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13
at ??:?

Did you have a solution? I have a similar error output, maybe you can help me. Thanks!