[lny200912]

Hi, FOAMERS,


I have met one error when I run my programming. In my programming, I redefined the way to calculate rho and phi in the following codes:

rho1 = theta_Psi * p0;
phi1 = fvc::interpolate(rho1) * fvc::interpolate(U) & mesh.Sf();


Then I used the rho1 and phi1 in my equation file:


fvScalarMatrix theaEqn
(
fvm::ddt(rho1, theta1)
+ mvConvection->fvmDiv(phi1, theta1)
- fvm::laplacian(turbulence->muEff(), theta1)
==
theta_1RR

);


And I modified discrete format in the sytem/fvScheme file:


ddtSchemes
{
// default Euler;
default backward;
}

gradSchemes
{
default Gauss linear;
}

divSchemes
{
default none;

div(phi,U) Gauss limitedLinearV 1;
div(phi1,U) Gauss limitedLinearV 1;

div(phi,Yi_h) Gauss limitedLinear 1;
div(phi1,Yi_h) Gauss limitedLinear 1;

div(phi,K) Gauss limitedLinear 1;

div(phid,p) Gauss limitedLinear 1;
div(phi,epsilon) Gauss limitedLinear 1;
// div(phi1,epsilon) Gauss limitedLinear 1;

div(phi,k) Gauss limitedLinear 1;
// div(phi1,k) Gauss limitedLinear 1;



However when I run my programming in the counterFlowFlame2D case, there is one error like this:


--> FOAM FATAL ERROR:
Operator + is undefined for unoriented and oriented types

From function Foam:rientedType Foam:perator+(const Foam:rientedType&, const Foam:rientedType&)
in file orientedType/orientedType.C at line 461.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam:perator+(Foam:rientedType const&, Foam:rientedType const&) at ??:?
#3 void Foam::add<double, double, Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<Foam::type OfSum<double, double>::type, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) at ??:?



Could someone can help me to explain the meaning of this error and how to solve it? Thank you so much.


Ning

[usv001]

Dear all,

I understand that this is an old thread but I faced a similar issue recently. While I don't understand the real meaning of this error, I found out that it occurred when trying to add the face flux (defined below) to some other quantity.

surfaceScalarField phi(Uf & mesh.Sf());

I found out that setting the flux to a nonOriented type, as shown below, solved the problem.

phi.setOriented(false);

Hope this helps future users. Maybe someone could explain what this error means and the correct way to handle it...

Best,
USV