[bijan darbari]

Dear all

I,m trying to incorporate Temperature field and the Energy equation in "DPMFoam" solver. So I do these main steps.

1- create "EEqn.H":

fvScalarMatrix EEqn
        (
            fvm::ddt(T)
            + fvm::div(alphaPhic, T)
            - fvm::laplacian(DT, T)
        );
EEqn.relax();

2- add "#include EEqn.H" in "DPMFoam.c",just after pimple loop

3- create temperature field and "DT" in "createFields.H"

IOdictionary transportProperties
(
    IOobject
    (
        "transportProperties",
        runTime.constant(),
        mesh,
        IOobject::MUST_READ_IF_MODIFIED,
        IOobject::NO_WRITE
    )
);
word continuousPhaseName
(
    transportProperties.lookup("continuousPhaseName")
);

dimensionedScalar DT
(
    transportProperties.lookup("DT")
);

Info<< "Reading field T\n" << endl;
volScalarField T
(
    IOobject
    (
        "T",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::AUTO_WRITE
    ),
    mesh
);

4- rename this solver to "DPMH"

this solver compiled successfully without any error.

but when I try to solve a simple test case, I get this error:

--> FOAM FATAL IO ERROR: 
wrong token type - expected word, found on line 19 the label 100

file: /home/bijan/Desktop/mycase1/constant/transportProperties.DT at line 19.

    From function operator>>(Istream&, word&)
    in file primitives/strings/word/wordIO.C at line 74.

FOAM exiting

The solver and test case attached.
DPMH.zip

mycase1.zip

Please kindly guide me.

[floquation]

I'm no expert, but given the error message I'd first try the following:

Change (constant/transportProperties)

DT 0.001;

into

DT.air 0.001;

as it is apparently looking for a "word" after "DT", I suspect ".air" must be written there.

[bijan darbari]

Thanks very much floquation.

I used this code and the case run without any problem. so calculate the temperature filed of continuous phase.

DT        DT [0 2 -1 0 0 0 0] 1.42e-7;

the modified solver:
DPMH.zip

but still my modified solver don't calculate the temperature field of particles.

in the following, I'll try to incorporate heat transfer equation to particles cloud.

[bijan darbari]

Dear all

As I mentioned earlier, I want to solve the heat transfer equation for
lagrangian particles in DPMFoam solver.

Fortunately, it can be done by openfoam standard "basicThermoCloud"
(without considering the colliding effect). Therefore, I do following
step in in my modified solver:

1- in "DPMH.C", I replace:

#ifdef MPPIC
    #include "basicThermoCloud.H"
    #define basicKinematicTypeCloud basicThermoCloud
#else
    #include "basicThermoCloud.H"
    #define basicKinematicTypeCloud basicThermoCloud
#endif

2- in "createFields", I add:

Info<< "Constructing kinematicCloud " << kinematicCloudName << endl;
basicKinematicTypeCloud basicThermoCloud
(
    kinematicCloudName,
    rhoc,
    Uc,
    muc,
    slgThermo,
   g
    
);

3- in "make/options", I add:

-I$(LIB_SRC)/thermophysicalModels/properties/liquidProperties/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/properties/liquidMixtureProperties/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/properties/solidProperties/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/properties/solidMixtureProperties/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/thermophysicalFunctions/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/reactionThermo/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/SLGThermo/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/chemistryModel/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/radiation/lnInclude \

and

 -lliquidProperties \
    -lliquidMixtureProperties \
    -lsolidProperties \
    -lsolidMixtureProperties \
    -lthermophysicalFunctions \
    -lreactionThermophysicalModels \
    -lSLGThermo \
    -lchemistryModel \
    -lradiationModels \

4- in the the beginning of "DPMH.C", I add:

#include "SLGThermo.H"

and I compile this solver. but I get this error:

In file included from DPMHP.C:68:0:
createFields.H: In function ‘int main(int, char**)’:
createFields.H:159:5: error: ‘slgThermo’ was not declared in this scope
     slgThermo,

What is the solution??

In the other words, How I can import the thermophysical properties of
continuous phase (K,Cp,...) to "basicThermoCloud" as input parameters
when I use the "constant transport properties" in the Constant Directory??

This solver have been attached.

[bijan darbari]

I'm looking forward to hear from you soon

[Scram_1]

Hey Bijan,
did you find a work-around for solving the energy equation in DPM Foam?

[tri3345]

Dear Bijan ,
could you solve the error?, I have a problem same as you, in a simulation I need the temperature of particles in DPMFoam, Can you guide me?