[JSSS]

So I was trying to run a model in simpleFoam (RAS k-E model) and the hope was to have a plant canopy as a source in the domain

I am getting the following error:

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  10
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 10-c4cf895ad8fa
Exec   : simpleFoam
Date   : Jun 06 2023
Time   : 08:10:21
Host   : "amun"
PID    : 3138812
I/O    : uncollated
Case   : 
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


SIMPLE: Convergence criteria found
        p: tolerance 0.0001
        U: tolerance 0.0001
        "(k|omega|epsilon)": tolerance 0.0001

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting viscosity model constant
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
RAS
{
    model           kEpsilon;
    turbulence      on;
    printCoeffs     on;
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              0;
    sigmak          1;
    sigmaEps        1.3;
}

No MRF models present

--> FOAM Warning : Creating fvModels from "constant/fvOptions"

Selecting finite volume model type atmPlantCanopyUSource


--> FOAM FATAL IO ERROR: 
Unknown fvModel atmPlantCanopyUSource

Valid fvModels are:

22
(
accelerationSource
actuationDiskSource
buoyancyEnergy
buoyancyForce
coded
effectivenessHeatExchangerSource
explicitPorositySource
heatSource
heatTransfer
interRegionExplicitPorositySource
interRegionHeatTransfer
isotropicDamping
massSource
phaseLimitStabilisation
radialActuationDiskSource
rotorDisk
semiImplicitSource
sixDoFAccelerationSource
solidEquilibriumEnergySource
solidificationMeltingSource
verticalDamping
volumeFractionSource
)

    From function static Foam::autoPtr<Foam::fvModel> Foam::fvModel::New(const Foam::word&, const Foam::dictionary&, const Foam::fvMesh&)
    in file cfdTools/general/fvModels/fvModel.C at line 117.

FOAM exiting

I understand that I need to include the lib in the controlDict which currently is as follows

FoamFile
{
    version         2;
    format          ascii;
    class           dictionary;
    location        "system";
    object          controlDict;
}

application     simpleFoam;

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         1000;

deltaT          0.1;

writeControl    timeStep;

writeInterval   10;

purgeWrite      0;

writeFormat     ascii;

writePrecision  7;

writeCompression on;

timeFormat      general;

timePrecision   6;

runTimeModifiable yes;

libs ("atmosphericModels.so");
libs ("libatmosphericModels.so");

are these the correct library names I need to include?


For context this is my fvOptions

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2212                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      fvOptions;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //



atmPlantCanopyUSource1
{
    // Mandatory entries (unmodifiable)
    type                  atmPlantCanopyUSource;
    atmPlantCanopyUSourceCoeffs
    {
        // Mandatory (inherited) entries (unmodifiable)
        selectionMode    barrier1Blockage;
        plantCd          0.2;   // Plant canopy drag coefficient          [-]
        leafAreaDensity  2;     // Leaf area density                      [1/m]
    }
}
atmPlantCanopyUSource2
{
    // Mandatory entries (unmodifiable)
    type                  atmPlantCanopyUSource;
    atmPlantCanopyUSourceCoeffs
    {
        // Mandatory (inherited) entries (unmodifiable)
        selectionMode    barrier2Blockage;
        plantCd          0.2;   // Plant canopy drag coefficient          [-]
        leafAreaDensity  2;     // Leaf area density                      [1/m]
    }
}

atmPlantCanopyTurbSource1
{
    // Mandatory entries (unmodifiable)
    type                  atmPlantCanopyTurbSource;
    atmPlantCanopyUSourceCoeffs
    {
        // Mandatory (inherited) entries (unmodifiable)
        selectionMode    barrier1Blockage;
        plantCd          0.2;   // Plant canopy drag coefficient              [-]
        leafAreaDensity  2;     // Leaf area density                          [1/m]
    }
}
atmPlantCanopyTurbSource2
{
    // Mandatory entries (unmodifiable)
    type                  atmPlantCanopyTurbSource;
    atmPlantCanopyUSourceCoeffs
    {
        // Mandatory (inherited) entries (unmodifiable)
        selectionMode    barrier2Blockage;
        plantCd          0.2;   // Plant canopy drag coefficient              [-]
        leafAreaDensity  2;     // Leaf area density                          [1/m]
    }
}

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

This is my first post here, please let me know if there's any more info that I can provide