[Vishsel]

Hi all,

I need to create heat flux (heat source) on the surface of the solid body.

Actually now my requirement is to run Incompressible fluid, Heat transfer, Steady state with porous media for the geometry with both solid and fluid. Till now I worked out for Incompressible fluid, Steady state with porous media only for fluid using solver porousSimpleFoam.

Now as an addition, I need to create heat source on the surface of the solid body (@wall).

For example - These are my wall surfaces where I have to create heat flux

wall-mp-1
wall-mp-2
wall-dp-1
wall-dp-2

And I have used cyclicAMI type for interfacing the fluid regions and for measure the pressure values.
So I succeed with cyclicAMI, now i have to create a heat source at wall(ie @surface).
I need to only check heat distribution from those BC's.Solver i am using is porousSimpleFoam.. running by parallel

And i have used porousSimpleFoam solver for fluid, and also have to see the temperature distribution in wall of the solid body. How could i run for both combined fluid and solid?
which Solver is suitable for my case?

For that i have to add energy equation EEqn.H with my case. I have an input as heat flux (W/m^2) and heat generation (W) value.
Q1) Is it possible to give these as an input in energy equation? once again, my case is incompressible fluid with steady state
Q2) buoyantBoussinesqPimpleFoam & buoyantBoussinesqSimpleFoam is for incompressible fluid but we can assign only in Temperature (T) term, not in heat flux or heat generation. is it correct? And buoyantSimpleFoam & buoyantPimpleFoam for compressible fluid.

Please guide me to add energy equation (EEqn.H) with my case in (fvOptions)


why I am facing this error? for solver porousSimpleFoam

Creating finite volume options from fvOptions

Selecting finite volume options model type scalarSemiImplicitSource
    Source: wall-mp-1
    - applying source for all time
    - selecting all cells
    - selected 5338536 cell(s) with volume 0.000333435

Selecting finite volume options model type scalarSemiImplicitSource
    Source: wall-mp-2
    - applying source for all time
    - selecting all cells
    - selected 5338536 cell(s) with volume 0.000333435

Selecting finite volume options model type scalarSemiImplicitSource
    Source: wall-dp-1
    - applying source for all time
    - selecting all cells
    - selected 5338536 cell(s) with volume 0.000333435

Selecting finite volume options model type scalarSemiImplicitSource
    Source: wall-dp-2
    - applying source for all time
    - selecting all cells
    - selected 5338536 cell(s) with volume 0.000333435

No MRF models present

Creating porosity model list from porosityProperties

Porosity region f-turbo:
    selecting model: DarcyForchheimer
    creating porous zone: f-turbo
Using pressure explicit porosity

Starting time loop

Time = 1

DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 0.0135476, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 0.00853288, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 0.00907314, No Iterations 1
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.00736855, No Iterations 8
time step continuity errors : sum local = 0.00220987, global = -0.000476259, cumulative = -0.000476259
DILUPBiCG:  Solving for epsilon, Initial residual = 0.0494679, Final residual = 0.000589308, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 1, Final residual = 0.0183162, No Iterations 1
ExecutionTime = 64.352 s  ClockTime = 64 s

Time = 2

--> FOAM Warning :
    From function void option::checkApplied() const
    in file fvOptions/fvOption.C at line 368
    Source wall-mp-1 defined for field h but never used

thermophysicalProperties for solid

thermoType
{
    type            heSolidThermo;
    mixture         pureMixture;
    transport       constIso;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleEnthalpy;
}

mixture
{
    specie
    {
        nMoles      1;
        molWeight   12;
    }

    transport
    {
        kappa   350;
    }

    thermodynamics
    {
        Hf      0;
        Cp      500;
    }

    equationOfState
    {
        rho     4000;
    }
}

fvOptions file

wall-mp-1
{
    type            scalarSemiImplicitSource;
    active          true;
    selectionMode   all;
 
    scalarSemiImplicitSourceCoeffs
    {
        volumeMode      absolute; 
        injectionRateSuSp //Semi-implicit source term S(x) = S_u + S_p x
        {
            h     (50 0);
        }
    }
}
wall-mp-2
{
    type            scalarSemiImplicitSource;
    active          true;
    selectionMode   all;
 
    scalarSemiImplicitSourceCoeffs
    {
        volumeMode      absolute; 
        injectionRateSuSp
        {
            h     (50 0);
        }
    }
}
wall-dp-1
{
    type            scalarSemiImplicitSource;
    active          true;
    selectionMode   all;
 
    scalarSemiImplicitSourceCoeffs
    {
        volumeMode      absolute; 
        injectionRateSuSp
        {
            h     (137.5 0);
        }
    }
}
wall-dp-2
{
    type            scalarSemiImplicitSource;
    active          true;
    selectionMode   all;
    
 
    scalarSemiImplicitSourceCoeffs
    {
        volumeMode      absolute;
        injectionRateSuSp
        {
            h     (137.5 0);
        }
    }
}

It will be very appreciable if you can help me out in this.

[Vishsel]

I still have no idea about this.
Please give me some hints if are seeing this post.

thank you very much!

[ashishvshelke]

You have added heat source in "all" regions