[svramana]

Hi Everyone,

I am new to OpenFoam. I am doing numerical simulation of fluid flow(incompressible)and Heat transfer from a porous square cylinder placed in the middle of a channel (maintained at constant Temp. > fluid Temp.)to a flowing wind.

For fluid flow(incompressible) simulation, the existing Poroussimplefoam(Darcy-Forchheimer Model) is giving fairly reasonable results on Cd and wake length.

Now I need to simulate heat transfer from porous square cylinder.
I found a thread on this forum where in, the procedure for adding TEqn. to icoFoam is available.

In my case the flow filed is only fluid and porosity is defined as source/sink ter using "topoSet.
The porous square cylinder is maintained at fixed Temp.using scalarFiXedValueConstarint in fvOptions
fvOptionsDict

s1
{
    type            scalarFixedValueConstraint;
    active          true;
    scalarFixedValueConstraintCoeffs
    {
        selectionMode   cellZone;
        cellZone        porousBlockage;
        volumeMode      uniform;
        fieldValues
        {
            T              353;
        }
    }
}

and TEqn. H included in PorousSimpleFoam solver

TEqn:fvScalarMatrix TEqn
(
 //fvm::ddt(T) // <--- steady state simulation, so ddt disabled...
//+
            fvm::div(phi, T)
         - fvm::laplacian(DT, T)
         ==
    fvOptions(T) // porouszone is maintained at const. temperature
        );
    TEqn.relax();
    fvOptions.constrain(TEqn);
    TEqn.solve();
    fvOptions.correct(T);

TEqn.H included in main solver file after governing equations. Compilation was done without any error msg.

following are the flow and heat transfer parameters considered
Re 40
Darcy no (Da ): 10e-2
porosity value: 0.993
Pr :0.71

I have few observations in this

1). There is no entry for porosity value as an input parameter.is it okay or am i missing something?

2). I want to calculate Nusselt Number ,i have general idea about calculation of Nusset number for solid cylinder, in my case the porous cylinder defined as source/sink (no physical surfaces/patches)

3). Is there available code for calculation Nusselt Number for porous blockage.?

I have searched many threads on this forum,but no luck.It will be great if someone can help me to run this simulation
Thank You,
Svramana