[dleduc]

Hi everyone,

I’m quite a beginner on OpenFoam and I would need some help about a solver modification I’m trying to make.

The aim

I would like to build a model that respresents the behavior of activated sludge inside a secondary clarifier (basicly I need to make a model about the sedimentation of a non-newtonian fluid).
For this purpose I chose the solver “driftFluxFoam” (which might be my first mistake…)
My problem is that in “driftFluxFoam” the non-newtonian behaviours follow laws from this type :

mu = A x 10^(b.alpha) [1]

Whereas, from the literature I found, the activated sludge viscosity is calculated from :

mu = K1/γ + K2 [2]

Where, γ is the strain rate.
I found that in the viscosityModels, the viscosity was calculated from this strain rate (through the strainRate function). I then tried to adapt the driftFluxFoam solver to solve [2].

The code
In the driftFluxFoam solver, I added, inside the mixtureViscosityModels :

• A “Casson” model – based on the Casson files from viscosityModel (see attached file);

• The strainRateFunction - based on the strainRateFunction files from viscosityMode.

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

#include "strainRateFunction.H"
#include "addToRunTimeSelectionTable.H"
#include "surfaceFields.H"

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

namespace Foam
{
namespace mixtureViscosityModels
{
    defineTypeNameAndDebug(strainRateFunction, 0);

    addToRunTimeSelectionTable
    (
        mixtureViscosityModel,
        strainRateFunction,
        dictionary
    );
}
}


// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //

Foam::mixtureViscosityModels::strainRateFunction::strainRateFunction
(
    const word& name,
    const dictionary& viscosityProperties,
    const volVectorField& U,
    const surfaceScalarField& phi
)
:
    mixtureViscosityModel(name, viscosityProperties, U, phi),
    strainRateFunctionCoeffs_
    (
        viscosityProperties.optionalSubDict(typeName + "Coeffs")
    ),
    strainRateFunction_
    (
        Function1<scalar>::New("function", strainRateFunctionCoeffs_)
    ),

    mu_
    (
        IOobject
        (
            "mu",
            U_.time().timeName(),
            U_.db()
        ),
        U_.mesh(),
        dimensionedScalar("mu", dimensionSet(1, -1, -1, 0, 0), 0),
        calculatedFvPatchScalarField::typeName
    )

{
    correct();
}


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

Foam::tmp<Foam::volScalarField>
Foam::mixtureViscosityModels::strainRateFunction::mu() const
{
    return mu_;
}


Foam::tmp<Foam::scalarField>
Foam::mixtureViscosityModels::strainRateFunction::mu(const label patchi) const
{
    return mu_.boundaryField()[patchi];
}


void Foam::mixtureViscosityModels::strainRateFunction::correct()
{
    tmp<volScalarField> tsigma = strainRate();
    const volScalarField& sigma = tsigma();

   mu_.primitiveFieldRef() = strainRateFunction_->value(sigma());

    volScalarField::Boundary& muBf = mu_.boundaryFieldRef();
    const volScalarField::Boundary& sigmaBf = sigma.boundaryField();

    forAll(muBf, patchi)
    {
        muBf[patchi] = strainRateFunction_->value(sigmaBf[patchi]);
    }
}


bool Foam::mixtureViscosityModels::strainRateFunction::read
(
    const dictionary& viscosityProperties
)
{
    mixtureViscosityModel::read(viscosityProperties);

    strainRateFunctionCoeffs_ = viscosityProperties.optionalSubDict
    (
        typeName + "Coeffs"
    );

    strainRateFunction_.clear();
    strainRateFunction_ = Function1<scalar>::New
    (
        "function",
        strainRateFunctionCoeffs_
    );

    return true;
}

Class
    Foam::viscosityModels::strainRateFunction


SourceFiles
    strainRateFunction.C

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

#ifndef strainRateFunction_H
#define strainRateFunction_H

#include "mixtureViscosityModel.H"
#include "volFields.H"
#include "Function1.H"
#include "plastic.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

namespace Foam
{
namespace mixtureViscosityModels
{

/*---------------------------------------------------------------------------*\
                           Class strainRateFunction Declaration
\*---------------------------------------------------------------------------*/

class strainRateFunction
:
    public mixtureViscosityModel
{
    // Private data

        //- Coefficients dictionary
        dictionary strainRateFunctionCoeffs_;

        //- Strain-rate function
        autoPtr<Function1<scalar>> strainRateFunction_;

        //- Current viscosity field
        const volScalarField mu_;



public:

    //- Runtime type information
    TypeName("strainRateFunction");


    // Constructors

        //- Construct from components
        strainRateFunction
        (
            const word& name,
            const dictionary& viscosityProperties,
            const volVectorField& U,
            const surfaceScalarField& phi
        );


    //- Destructor
    virtual ~strainRateFunction()
    {}


    // Member Functions

        //- Return the viscosity
        virtual tmp<volScalarField> mu() const;

        //- Return the viscosity for patch
        virtual tmp<scalarField> mu(const label patchi) const;

        //- Correct the laminar viscosity
        virtual void correct();

        //- Read transportProperties dictionary
        virtual bool read(const dictionary& viscosityProperties);
};


// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

} // End namespace viscosityModels
} // End namespace Foam

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

#endif

// ************************************************************************* //

You will also find the mixtureViscosityModel code which I modified a bit :

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

#include "mixtureViscosityModel.H"
#include "volFields.H"
#include "fvcGrad.H"

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

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


// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //

Foam::mixtureViscosityModel::mixtureViscosityModel
(
    const word& name,
    const dictionary& viscosityProperties,
    const volVectorField& U,
    const surfaceScalarField& phi
)
:
    name_(name),
    viscosityProperties_(viscosityProperties),
    U_(U),
    phi_(phi)
{}


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

Foam::tmp<Foam::volScalarField> Foam::mixtureViscosityModel::strainRate() const
{
    return sqrt(2.0)*mag(symm(fvc::grad(U_)));
}


bool Foam::mixtureViscosityModel::read(const dictionary& viscosityProperties)
{
    viscosityProperties_ = viscosityProperties;

    return true;
}


// ************************************************************************* //

SourceFiles
    mixtureViscosityModel.C
    mixtureViscosityModelNew.C

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

#ifndef mixtureViscosityModel_H
#define mixtureViscosityModel_H

#include "dictionary.H"
#include "volFieldsFwd.H"
#include "surfaceFieldsFwd.H"
#include "dimensionedScalar.H"
#include "runTimeSelectionTables.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

namespace Foam
{

/*---------------------------------------------------------------------------*\
                Class mixtureViscosityModel Declaration
\*---------------------------------------------------------------------------*/

class mixtureViscosityModel
{

protected:

    // Protected data

        word name_;
        dictionary viscosityProperties_;

        const volVectorField& U_;
        const surfaceScalarField& phi_;


    // Private Member Functions

        //- Disallow copy construct
        mixtureViscosityModel(const mixtureViscosityModel&);

        //- Disallow default bitwise assignment
        void operator=(const mixtureViscosityModel&);


public:

    //- Runtime type information
    TypeName("mixtureViscosityModel");


    // Declare run-time constructor selection table

        declareRunTimeSelectionTable
        (
            autoPtr,
            mixtureViscosityModel,
            dictionary,
            (
                const word& name,
                const dictionary& viscosityProperties,
                const volVectorField& U,
                const surfaceScalarField& phi
            ),
            (name, viscosityProperties, U, phi)
        );


    // Selectors

        //- Return a reference to the selected viscosity model
        static autoPtr<mixtureViscosityModel> New
        (
            const word& name,
            const dictionary& viscosityProperties,
            const volVectorField& U,
            const surfaceScalarField& phi
        );


    // Constructors

        //- Construct from components
        mixtureViscosityModel
        (
            const word& name,
            const dictionary& viscosityProperties,
            const volVectorField& U,
            const surfaceScalarField& phi
        );


    //- Destructor
    virtual ~mixtureViscosityModel()
    {}


    // Member Functions

        //- Return the phase transport properties dictionary
        const dictionary& viscosityProperties() const
        {
            return viscosityProperties_;
        }

        //- Return the strain rate
        tmp<volScalarField> strainRate() const;

        //- Return the mixture viscosity
        //  given the viscosity of the continuous phase
        virtual tmp<volScalarField> mu(const volScalarField& muc) const;
        //virtual tmp<volScalarField> mu(const volScalarField& muc) const = 0;

        //- Return the laminar viscosity for patch
        //virtual tmp<scalarField> mu(const label patchi) const = 0;   Cet ajout fait bugger la compilation de plastic.C...

        //- Read transportProperties dictionary
        virtual bool read(const dictionary& viscosityProperties) = 0;
};


// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

}

The problem
When I try to compile the solver, I’m having some troubles with a part of the strainRateFunction (see attached file)
From what I understood, this problems comes from the fact that my code risks to modify a variable that shouldn’t be modified (const). But I have no idea about how to fix this…

What I have tried
I tried to avoid the parts of the code at the origin of my problem (with “//”), then the solver compiles, but when I try to run it, I’m getting an error message (“symbol look up error”).
I also tried to make my solver without the strainRateFunction replacing “strainRate()” with “sqrt(2.0)*mag(symm(fvc::grad(U_)))” directly inside “Casson.C”. The solver also compiles but crashes when I try to run it.
This time the problem came from LHS and RHS which didn’t have the same dimension.
As I said I’m very new to OpenFoam, and to be perfectly honest, I don’t really understand how the programming process… That is why I expect quite some huge mistakes and apologize about this…


Finally, I don't know if it's usefull, but attached the code about incompressibleTwoPhaseInteratingmixture.


As I said, maybe my first mistake was to choose driftFluxFoam… If you know another solver which would enable be to solve my problem, it would be a pleasure to try it !
I know it's quite time consuming to read all this... So thanks in advance for your help !

[dleduc]

Problem solved...

When I declared my variable mu_ in strainRateFunction.H, I used "const volScalarFiel mu_" where I should only have used "volScalarFiel mu_" in order to enable the correct() function to modify the mu_ variable !

[dleduc]

Unfortunately, I'm running into another type of error :

symbol lookup error: /home/leduc/OpenFOAM/leduc-5.0/platforms/linux64GccDPInt32Opt/lib/libdriftFluxTransportModels.so: undefined symbol: _ZNK4Foam22mixtureViscosityModels18strainRateFunction2muERKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEE

When I check, where it comes from :

c++filt _ZNK4Foam22mixtureViscosityModels18strainRateFunction2muERKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEE
Foam::mixtureViscosityModels::strainRateFunction::mu(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const

I believe the part of my code which causes this error is in my file strainRateFunction.C :

Foam::tmp<Foam::volScalarField>
Foam::mixtureViscosityModels::strainRateFunction::mu() const
{
    return mu_;
}


Foam::tmp<Foam::scalarField>
Foam::mixtureViscosityModels::strainRateFunction::mu(const label patchi) const
{
    return mu_.boundaryField()[patchi];
}

Do you have any idea about the source of this particular problem ?
Any idea about how to fix this ?

[mAlletto]

did you try to include the library you compiled in the controlDict

[dleduc]

Hi mAlleto,

Thank you very much for your answer, I just tried this :

oamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

libs            ("home/OpenFOAM/leduc-5.0/platforms/linux64GccDPInt32Opt/lib/libdriftFluxRelativeVelocityModels.so" "home/OpenFOAM/leduc-5.0/platforms/linux64GccDPInt32Opt/lib/libdriftFluxTransportModels.so");

application     driftFluxFoam_Casson;

Unfortunatly, it didn't solve the problem.
From what I read, generally those problems come from the wmake files.

Here are the wmake files of my solver :

File

incompressibleTwoPhaseInteractingMixture/incompressibleTwoPhaseInteractingMixture.C
compressibleTurbulenceModels.C
driftFluxFoam_Casson.C

EXE = $(FOAM_USER_APPBIN)/driftFluxFoam_Casson
LIB = $(FOAM_USER_LIBBIN)/driftFluxFoam_Casson

Options

"EXE_INC = \
    -IincompressibleTwoPhaseInteractingMixture \
    -ImixtureViscosityModels/lnInclude \
    -IrelativeVelocityModels/lnInclude \
    -I$(LIB_SRC)/finiteVolume/lnInclude \
    -I$(LIB_SRC)/meshTools/lnInclude \
    -I$(LIB_SRC)/sampling/lnInclude \
    -I$(LIB_SRC)/transportModels \
    -I$(LIB_SRC)/transportModels/compressible/lnInclude \
    -I$(LIB_SRC)/transportModels/twoPhaseMixture/lnInclude \
    -I$(LIB_SRC)/TurbulenceModels/turbulenceModels/lnInclude \
    -I$(LIB_SRC)/TurbulenceModels/compressible/lnInclude

EXE_LIBS = \
    -ldriftFluxRelativeVelocityModels \
    -ldriftFluxTransportModels \
    -L$(FOAM_USER_LIBBIN) \
    -ltwoPhaseMixture \
    -lfiniteVolume \
    -lmeshTools \
    -lsampling \
    -lfvOptions \
    -lcompressibleTransportModels \
    -lturbulenceModels \
    -lcompressibleTurbulenceModels

Also, the incompressibleTwoPhaseInteractingMixture.C calls for mixtureViscosityModel.H which is in the mixtureViscosityModels folder.
The wmake files for the mixtureViscosityModels solver are :

Files

mixtureViscosityModel/mixtureViscosityModel.C
mixtureViscosityModel/mixtureViscosityModelNew.C
Casson/Casson.C
strainRateFunction/strainRateFunction.C

LIB = $(FOAM_USER_LIBBIN)/libdriftFluxTransportModels

Options

EXE_INC = \
    -I../incompressibleTwoPhaseInteractingMixture \
    -I$(LIB_SRC)/finiteVolume/lnInclude \
    -I$(LIB_SRC)/transportModels \
    -I$(LIB_SRC)/transportModels/twoPhaseMixture/lnInclude \
    -I$(LIB_SRC)/transportModels/incompressible/lnInclude \


LIB_LIBS = \
    -ltwoPhaseMixture \
    -lfiniteVolume

The error I get tells me that it comes from the compiledlibdriftFluxTransportModels.so file. More precisely from the strainRateFunction.C (see previous post).

I'm not really used to link different libraries, and I suspect an error in my options files but can't find it...

Thanks again for the time you took mAlleto !

[mAlletto]

Originally Posted by dleduc

Unfortunately, I'm running into another type of error :

symbol lookup error: /home/leduc/OpenFOAM/leduc-5.0/platforms/linux64GccDPInt32Opt/lib/libdriftFluxTransportModels.so: undefined symbol: _ZNK4Foam22mixtureViscosityModels18strainRateFunction2muERKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEE

When I check, where it comes from :

c++filt _ZNK4Foam22mixtureViscosityModels18strainRateFunction2muERKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEE
Foam::mixtureViscosityModels::strainRateFunction::mu(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const

I believe the part of my code which causes this error is in my file strainRateFunction.C :

Foam::tmp<Foam::volScalarField>
Foam::mixtureViscosityModels::strainRateFunction::mu() const
{
    return mu_;
}


Foam::tmp<Foam::scalarField>
Foam::mixtureViscosityModels::strainRateFunction::mu(const label patchi) const
{
    return mu_.boundaryField()[patchi];
}

Do you have any idea about the source of this particular problem ?
Any idea about how to fix this ?

The error you get means you declared the function but you don't defined it.

Did you define the function:

mu(const volScalarField& muc) const;

somewhere in the code a pice of code tries the execute the function mu which takes a const volScalarField as argument but it does not find it. So the error is generated from it.

Or better said it has found the declaration but not the definition. So I ask if you declared it

[dleduc]

You were perfectly right :

In the mixtureViscosityModel.H, I declared :

        //- Return the mixture viscosity
        //- given the viscosity of the continuous phase
        virtual tmp<volScalarField> mu(const volScalarField &muc) const = 0

and in the strainRateFunction.H I declared :

        //- Return the mixture viscosity
        //  given the viscosity of the continuous phase
        virtual tmp<volScalarField> mu(const volScalarField& muc) const

whereas in strainRateFunction.C mu(const volScalarField& muc) was not defined...

I solved the problem by supressing the declaration of mu(const volScalarField& muc) in strainRateFunction.H and by modifying the declaration in mixtureViscosityModel.H as follows :

mixtureViscosityModel.H

      //- Return the mixture viscosity
        //- given the viscosity of the continuous phase
        virtual tmp<volScalarField> mu(const volScalarField &muc) const;

Thank you very much mAlletto !!

Unfortunately, as soon as a problem is solved, another comes up... The error message I'm getting now is :

driftFluxFoam_Casson: symbol lookup error: /home/leduc/OpenFOAM/leduc-5.0/platforms/linux64GccDPInt32Opt/lib/libdriftFluxTransportModels.so: undefined symbol: _ZTIN4Foam21mixtureViscosityModelE

c++filt  _ZTIN4Foam21mixtureViscosityModelE 
typeinfo for Foam::mixtureViscosityModel

It seems like the new problem is similar to the one with mu() definition... So, I checked definition and declaration of mixtureViscosityModel but it is exactly the same as the initial mixtureViscosityModel in driftFluxFoam solver... I then don't understand why the typeinfo would be missing...

[mAlletto]

Originally Posted by dleduc

Hi mAlleto,

Thank you very much for your answer, I just tried this :

oamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

libs            ("home/OpenFOAM/leduc-5.0/platforms/linux64GccDPInt32Opt/lib/libdriftFluxRelativeVelocityModels.so" "home/OpenFOAM/leduc-5.0/platforms/linux64GccDPInt32Opt/lib/libdriftFluxTransportModels.so");

application     driftFluxFoam_Casson;

Unfortunatly, it didn't solve the problem.
From what I read, generally those problems come from the wmake files.

Here are the wmake files of my solver :

File

incompressibleTwoPhaseInteractingMixture/incompressibleTwoPhaseInteractingMixture.C
compressibleTurbulenceModels.C
driftFluxFoam_Casson.C

EXE = $(FOAM_USER_APPBIN)/driftFluxFoam_Casson
LIB = $(FOAM_USER_LIBBIN)/driftFluxFoam_Casson

Options

"EXE_INC = \
    -IincompressibleTwoPhaseInteractingMixture \
    -ImixtureViscosityModels/lnInclude \
    -IrelativeVelocityModels/lnInclude \
    -I$(LIB_SRC)/finiteVolume/lnInclude \
    -I$(LIB_SRC)/meshTools/lnInclude \
    -I$(LIB_SRC)/sampling/lnInclude \
    -I$(LIB_SRC)/transportModels \
    -I$(LIB_SRC)/transportModels/compressible/lnInclude \
    -I$(LIB_SRC)/transportModels/twoPhaseMixture/lnInclude \
    -I$(LIB_SRC)/TurbulenceModels/turbulenceModels/lnInclude \
    -I$(LIB_SRC)/TurbulenceModels/compressible/lnInclude

EXE_LIBS = \
    -ldriftFluxRelativeVelocityModels \
    -ldriftFluxTransportModels \
    -L$(FOAM_USER_LIBBIN) \
    -ltwoPhaseMixture \
    -lfiniteVolume \
    -lmeshTools \
    -lsampling \
    -lfvOptions \
    -lcompressibleTransportModels \
    -lturbulenceModels \
    -lcompressibleTurbulenceModels

Also, the incompressibleTwoPhaseInteractingMixture.C calls for mixtureViscosityModel.H which is in the mixtureViscosityModels folder.
The wmake files for the mixtureViscosityModels solver are :

Files

mixtureViscosityModel/mixtureViscosityModel.C
mixtureViscosityModel/mixtureViscosityModelNew.C
Casson/Casson.C
strainRateFunction/strainRateFunction.C

LIB = $(FOAM_USER_LIBBIN)/libdriftFluxTransportModels

Options

EXE_INC = \
    -I../incompressibleTwoPhaseInteractingMixture \
    -I$(LIB_SRC)/finiteVolume/lnInclude \
    -I$(LIB_SRC)/transportModels \
    -I$(LIB_SRC)/transportModels/twoPhaseMixture/lnInclude \
    -I$(LIB_SRC)/transportModels/incompressible/lnInclude \


LIB_LIBS = \
    -ltwoPhaseMixture \
    -lfiniteVolume

The error I get tells me that it comes from the compiledlibdriftFluxTransportModels.so file. More precisely from the strainRateFunction.C (see previous post).

I'm not really used to link different libraries, and I suspect an error in my options files but can't find it...

Thanks again for the time you took mAlleto !

By the way why did you put

incompressibleTwoPhaseInteractingMixture/incompressibleTwoPhaseInteractingMixture.C
compressibleTurbulenceModels.C
driftFluxFoam_Casson.C

EXE = $(FOAM_USER_APPBIN)/driftFluxFoam_Casson
LIB = $(FOAM_USER_LIBBIN)/driftFluxFoam_Casson

the lib in the Make/files ?

usually you put only the EXE

[dleduc]

The error I get tells me that it comes from the compiledlibdriftFluxTransportModels.so file. More precisely from the strainRateFunction.C (see previous post).

I'm not really used to link different libraries, and I suspect an error in my options files but can't find it...

Thanks again for the time you took mAlleto !

The LIB is what remains of what I tried to do to solve my problem. When you gave me the solution I forgot to remove it. It is done now !

For the actual error :

driftFluxFoam_Casson: symbol lookup error: /home/leduc/OpenFOAM/leduc-5.0/platforms/linux64GccDPInt32Opt/lib/libdriftFluxTransportModels.so: undefined symbol: _ZTIN4Foam21mixtureViscosityModelE

c++filt  _ZTIN4Foam21mixtureViscosityModelE 
typeinfo for Foam::mixtureViscosityModel

Do you think that it could come from my libraries paths which could be wrong ?
Or is it again a declaration problem in the code ?

Just in case... Here's a wetransfer link to download my solver an the test case I use to see if it works : https://we.tl/t-7VZHXDPdKU

[mAlletto]

what version of OF did you use to compile your code?

[dleduc]

Hi mAlleto,


I'm using OF5.0 ! Do you think I should use the last version ?

[mAlletto]

I compiled your code with OF6 and a bunch of errors popped up. Probably some classes have changed from OF5 and OF6.



Can you try to post a code which compiles with OF6? Otherwise it would take me too long to figure out what's going wrong.


best


michael

[mAlletto]

Ok i was able to compile your code with OF6 by changing alphaEqn.H with the one of OF. I get the the same error as you when executing it.



when looking through Make/options I found that you included twice a the viscosity model:




EXE_INC = \
-IincompressibleTwoPhaseInteractingMixture \
-ImixtureViscosityModels/lnInclude \
-I./relativeVelocityModels/lnInclude \
-I$(LIB_SRC)/finiteVolume/lnInclude \
-I$(LIB_SRC)/meshTools/lnInclude \
-I$(LIB_SRC)/sampling/lnInclude \
-I$(LIB_SRC)/transportModels \
-I$(LIB_SRC)/transportModels/compressible/lnInclude \
-I$(LIB_SRC)/transportModels/twoPhaseMixture/lnInclude \
-I$(LIB_SRC)/TurbulenceModels/turbulenceModels/lnInclude \
-I$(LIB_SRC)/TurbulenceModels/compressible/lnInclude

EXE_LIBS = \
-ldriftFluxTransportModels \ -> the one of OF
-ldriftFluxRelativeVelocityModels \
-ltwoPhaseMixture \
-lfiniteVolume \
-lmeshTools \
-lsampling \
-lfvOptions \
-lcompressibleTransportModels \
-lturbulenceModels \
-lcompressibleTurbulenceModels \
-L$(FOAM_USER_LIBBIN)/libdriftFluxTransportModels -> your


when removing the library of OF i get similar errors as you but already at compile time. There must be a naming conflict somewhere.

[dleduc]

Hi mAlletto,

indeed there was a naming problem... I began to suspect it when I tried to run my test case with the original driftFluxFoam and got the same error...

driftFluxFoam tried to access my user library :

/home/leduc/OpenFOAM/leduc-5.0/platforms/linux64GccDPInt32Opt/lib/libdriftFluxTransportModels.so

For some reason, even if not stored in the same place, the fact that the original and the user library had the same name was the source of some confusion...

I then changed the names of my solver and library to Cassonsed and libCassonsed.so and the error went away...

Now the error I get is quite more basic :

-> FOAM FATAL ERROR: 
Unknown mixtureViscosityModel type transportModel

Valid mixtureViscosityModels are : 

3
(
BinghamPlastic
plastic
slurry
)

This error is the very reason why I wanted to compile a new solver... I thought that by compiling a solver only with my Casson viscosity model it would be recognised by Foam::mixtureViscosityModel::New, but it's unfortunatly not the case... For some reason, my Casson model is not considered as a transportModel :

{
    const word modelType(viscosityProperties.lookup("transportModel"));

    Info<< "Selecting incompressible transport model " << modelType << endl;

    dictionaryConstructorTable::iterator cstrIter =
        dictionaryConstructorTablePtr_->find(modelType);

    if (cstrIter == dictionaryConstructorTablePtr_->end())
    {
        FatalErrorInFunction
            << "Unknown viscosityModel type "
            << modelType << nl << nl
            << "Valid viscosityModels are : " << endl
            << dictionaryConstructorTablePtr_->sortedToc()
            << exit(FatalError);
    }

I don't know how to implement my Casson viscosity model so it can be found in the dictionaryConstructorTablePtr_...

Thank you very much for what you're doing, you're helping me a lot !

[mAlletto]

Just looking at the other viscosity models I would have just copied one and renamed it.

Afterward I would have modified the mu and read function (letting the input arguments as they are).


By the way you don't need any strain rate functions. You can use the velocity U to calculate it.

By compiling the new viscosity model it should by selectable at run time.

[dleduc]

Yes, I already tried that :

I also tried to make my solver without the strainRateFunction replacing “strainRate()” with “sqrt(2.0)*mag(symm(fvc::grad(U_)))” directly inside “Casson.C”. The solver also compiles but crashes when I try to run it. 
This time the problem came from LHS and RHS which didn’t have the same dimension.

I admit that I wanted to do something "cleaner" with a new solver (it would have been some kind of practice with OpenFoam programmation)...

With my C++ level, it might indeed be more appriopriate just to modify the plastic solver just to fulfill my needs !
I'll try this again and hope I'll be able to find the solution about the LHS and RHS dimensions !

Thanks a lot michael !