Gauss CoBlended scheme fails for div(phi,alpha) in interFoam
Summary
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. Logfile output:
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.
What is the expected correct behavior?
Normal interFoam-running would be expected.
Relevant logs and/or images
Environment information
OpenFOAM version :v1812 Operating system :centos Compiler :gcc
Possible fixes
## Reattaching the author to the issue ticket: @Lydia ##