[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

[MMRC]

Hi JSSS, did you solve this?
I'm taking a look and LAD and Cd are defined as scalarfied to apply as B.C, so it looks like they don't have to be declared inside the fvOption,

Did you find your way to use this tool?

[JSSS]

Originally Posted by MMRC

Hi JSSS, did you solve this?
I'm taking a look and LAD and Cd are defined as scalarfied to apply as B.C, so it looks like they don't have to be declared inside the fvOption,

Did you find your way to use this tool?

Hello MMRC,

Yes I did. I was using OpenFOAM.org when I posted this which apparently doesn't have these fvOptions available (as seen in the list provided in the error raised, they are not there). The fix to this was switching from the OpenFOAM.org implementation to the OpenFOAM.com one.

Thank you for posting, I had completely forgotten to update my situation here for people who may potentially come across the same roadblock