[Utkan]

I am a newbie in OpenFOAM and I desperately need any type of help or guidance in my situation. I want to simulate the forced cooling of a radiator with air by a fan. I am working with the heatExchanger case under chtMultiRegionFoam. However, the amount of fin/capillary details of this radiator makes it impossible to mesh with a commercial meshing tool. Hence, I'm going for the porosity approach. Water is considered a porous fluid of which the porous parameters are given in the fvOptions. I will include fvOptions of both regions below. The momentum source term due to the porosity blockage is I think already being added correctly, but from what I understood, I have to find a way to add the energy source term to get rid of the warning. But I don't know how to do that. I would really appreciate any help overcoming this error. Thank you in advance.


This is the error that I get, when I try to solve the case

Solving for fluid region air
--> FOAM Warning : 
    From function virtual void Foam::fv::option::checkApplied() const
    in file cfdTools/general/fvOptions/fvOption.C at line 125
    Source airToporous defined for field h but never used
--> FOAM Warning : 
    From function virtual void Foam::fv::option::checkApplied() const
    in file cfdTools/general/fvOptions/fvOption.C at line 125
    Source airToporous defined for field h but never used
DILUPBiCGStab:  Solving for Ux, Initial residual = 0.474622, Final residual = 0.00170632, No Iterations 1
DILUPBiCGStab:  Solving for Uy, Initial residual = 0.403169, Final residual = 0.00159073, No Iterations 1
DILUPBiCGStab:  Solving for Uz, Initial residual = 0.573588, Final residual = 0.00181622, No Iterations 1
--> FOAM Warning : 
    From function virtual void Foam::fv::option::checkApplied() const
    in file cfdTools/general/fvOptions/fvOption.C at line 125
    Source airToporous defined for field h but never used
--> FOAM Warning : 
    From function virtual void Foam::fv::option::checkApplied() const
    in file cfdTools/general/fvOptions/fvOption.C at line 125
    Source airToporous defined for field h but never used
DILUPBiCGStab:  Solving for e, Initial residual = 0.00898525, Final residual = 5.40341e-05, No Iterations 1
Min/max T:292.92 325.612
--> FOAM Warning : 
    From function virtual void Foam::fv::option::checkApplied() const
    in file cfdTools/general/fvOptions/fvOption.C at line 125
    Source airToporous defined for field h but never used
GAMG:  Solving for p_rgh, Initial residual = 0.441588, Final residual = 0.00311991, No Iterations 12
time step continuity errors : sum local = 0.000507967, global = 3.8117e-05, cumulative = 3.8117e-05
--> FOAM Warning : 
    From function virtual void Foam::fv::option::checkApplied() const
    in file cfdTools/general/fvOptions/fvOption.C at line 125
    Source airToporous defined for field h but never used
--> FOAM Warning : 
    From function virtual void Foam::fv::option::checkApplied() const
    in file cfdTools/general/fvOptions/fvOption.C at line 125
    Source airToporous defined for field h but never used
DILUPBiCGStab:  Solving for epsilon, Initial residual = 0.0690558, Final residual = 0.000399087, No Iterations 1
--> FOAM Warning : 
    From function virtual void Foam::fv::option::checkApplied() const
    in file cfdTools/general/fvOptions/fvOption.C at line 125
    Source airToporous defined for field h but never used
--> FOAM Warning : 
    From function virtual void Foam::fv::option::checkApplied() const
    in file cfdTools/general/fvOptions/fvOption.C at line 125
    Source airToporous defined for field h but never used

the fvOptions for both regions is as follows;

For air region;

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

airToporous
{
    type            constantHeatTransfer;

    interpolationMethod cellVolumeWeight;
    nbrRegionName   porous;
    master          false;

    nbrModel        porousToair;
    fields          (h);
    semiImplicit    no;
}

porosityBlockage
{
    type            interRegionExplicitPorositySource;

    interRegionExplicitPorositySourceCoeffs
    {
        interpolationMethod cellVolumeWeight;
        nbrRegionName   porous;

        type            DarcyForchheimer;

        d   (-1000 -1000 1e4);
        f   (0 0 0);

        coordinateSystem
        {
            type    cartesian;
            origin  (0 0 0);
            coordinateRotation
            {
                type    axesRotation;
                e1      (0 1 0);
                e2      (0 0 1);
            }
        }
    }
}


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

For the porous region;

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

porousToair
{
    type            constantHeatTransfer;

    interpolationMethod cellVolumeWeight;
    nbrRegionName   air;
    master          true;

    nbrModel        airToporous;
    fields          (e);
    semiImplicit    no;
}


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

[Raphael_Santos]

Hello Utkan,

have you find a solution to that issue?

Cheers,

[Pappelau]

he defined energy as h in fvoptions but is using in the theromprops the formulation with epsilon this visible in the solver log: solving for epsilon .... this cant work out