Skip to content

Re-open Issue: Gauss CoBlended scheme fails for div(phi,alpha) in interFoam

Summary

Note: I re-open this issue, as it is still relevant and not fixed in the latest version. The issue was already described in #1236 (closed) for version v1812.

When using the scheme Gauss CoBlended for the term div(phi,alpha) in the solver interFoam, the simulation breaks up with the following error message: fvSchemes:

div(phi,alpha) Gauss CoBlended 1. Minmod 2. upwind;

logfile output:

\--> FOAM FATAL ERROR: Operator + is undefined for unoriented and oriented types From function Foam::orientedType Foam::operator+(const Foam::orientedType&, const Foam::orientedType&) in file orientedType/orientedType.C at line 458.

Steps to reproduce

Use tutorial /multiphase/RAS/interFoam/damBreak. Change in fvSchemes:

div(phi,alpha) Gauss CoBlended 1. Minmod 2. upwind;

Run the case with Allrun.

Example case

The damBreak tutorial (tutorials/multiphase/RAS/damBreak)

What is the current bug behaviour?

Run fails in first alpha-Subcycle. See logfile output below.

What is the expected correct behavior?

Normal interFoam-running would be expected.

Relevant logs and/or images

Courant Number mean: 0 max: 0 Interface Courant Number mean: 0 max: 0 deltaT = 0.00119048 Time = 0.00119048 PIMPLE: iteration 1 smoothSolver: Solving for alpha.water, Initial residual = 0, Final residual = 0, No Iterations 0 Phase-1 volume fraction = 0.130194 Min(alpha.water) = 0 Max(alpha.water) = 1 --> FOAM FATAL ERROR: Operator + is undefined for unoriented and oriented types From function Foam::orientedType Foam::operator+(const Foam::orientedType&, const Foam::orientedType&) in file orientedType/orientedType.C at line 458.

Environment information

  • OpenFOAM version : v1912
  • Operating system : centos
  • Hardware info :
  • Compiler : gcc

Possible fixes