[ZeinabM]

Hello, I hope you're doing well. So, now there is a task required from me to edit rhoSimpleFoam solver to accommodate the SRF model. I copied the solver to my application file and performed some edits. Also, I have done some edits in the SRFModel.H and SRFModel.C. I compile the new solver and everything is good, when I get to use it with a problem, it gives me an error of inconsistent units. I've successfully known where the inconsistency happens but I don't why.

--> FOAM FATAL ERROR: (openfoam-2012)

    [Urel[1 -2 -2 0 0 0 0] ] + [(Fcoriolis+Fcentrifugal)[0 1 -2 0 0 0 0] ]

    From void Foam::checkMethod(const Foam::fvMatrix<Type>&, const Foam::DimensionedField<Type, Foam::volMesh>&, const char*) [with Type = Foam::Vector<double>]
    in file /opt/OpenFOAM/OpenFOAM-v2012/src/finiteVolume/lnInclude/fvMatrix.C at line 1351.

FOAM aborting

I put some flags to know where the error is, and it turned out that it happens when solving the Urel equation.

  tmp<fvVectorMatrix> tUrelEqn
   (
        fvm::div(phi, Urel)
        + SRF -> SuCompressible()
        + turbulence->divDevRhoReff(Urel)
     ==
        fvOptions(rho, Urel)
    );

Below I am sharing SRFModel.C, the lines in bold are the edited lines:

\*---------------------------------------------------------------------------*/

#include "SRFModel.H"
#include "SRFVelocityFvPatchVectorField.H"

// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //

namespace Foam
{
    namespace SRF
    {
        defineTypeNameAndDebug(SRFModel, 0);
        defineRunTimeSelectionTable(SRFModel, dictionary);
    }
}


// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
Foam::SRF::SRFModel::SRFModel
(
    const word& type,
    const volVectorField& Urel,
    const volScalarField& rho 
)
:
    IOdictionary
    (
        IOobject
        (
            "SRFProperties",
            Urel.time().constant(),
            Urel.db(),
            IOobject::MUST_READ_IF_MODIFIED,
            IOobject::NO_WRITE
        )
    ),
    Urel_(Urel),
    rho_ (rho),
    mesh_(Urel_.mesh()),
    origin_("origin", dimLength, get<vector>("origin")),
    axis_(normalised(get<vector>("axis"))),
    SRFModelCoeffs_(optionalSubDict(type + "Coeffs")),
    omega_(dimensionedVector("omega", dimless/dimTime, Zero))
{}


// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //

Foam::SRF::SRFModel::~SRFModel()
{}


// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

bool Foam::SRF::SRFModel::read()
{
    if (regIOobject::read())
    {
        // Re-read origin
        readEntry("origin", origin_);

        // Re-read axis
        readEntry("axis", axis_);
        axis_.normalise();

        // Re-read sub-model coeffs
        SRFModelCoeffs_ = optionalSubDict(type() + "Coeffs");

        return true;
    }

    return false;
}


const Foam::dimensionedVector& Foam::SRF::SRFModel::origin() const
{
    return origin_;
}


const Foam::vector& Foam::SRF::SRFModel::axis() const
{
    return axis_;
}


const Foam::dimensionedVector& Foam::SRF::SRFModel::omega() const
{
    return omega_;
}


Foam::tmp<Foam::DimensionedField<Foam::vector, Foam::volMesh>>
Foam::SRF::SRFModel::FcoriolisCompressible() const
{
    return tmp<volVectorField::Internal>
    (
        new volVectorField::Internal
        (
            IOobject
            (
                "FcoriolisComp",
                mesh_.time().timeName(),
                mesh_,
                IOobject::NO_READ,
                IOobject::NO_WRITE
            ),
            rho_* 2.0*omega_ ^ Urel_
        )
    );
}

Foam::tmp<Foam::DimensionedField<Foam::vector, Foam::volMesh>>
Foam::SRF::SRFModel::FcentrifugalCompressible() const
{
    return tmp<volVectorField::Internal>
    (
        new volVectorField::Internal
        (
            IOobject
            (
                "FcentrifugalComp",
                mesh_.time().timeName(),
                mesh_,
                IOobject::NO_READ,
                IOobject::NO_WRITE
            ),
            rho_ *omega_ ^ (omega_ ^ (mesh_.C() - origin_))
        )
    );
}

Foam::tmp<Foam::DimensionedField<Foam::vector, Foam::volMesh>>
Foam::SRF::SRFModel::SuCompressible() const
{
    return FcoriolisCompressible() + FcentrifugalCompressible();
}

The problem might be in the process of incorporating the rho in the SRFModel files, but I don't know exactly where it is.

[Ununul]

The problem is coming from dimensions inconsistency. Are you sure that your Urel variable should have dimension kg m^-2 s^-2?

Greetings