[himanshu28]

Hi,

I am trying to add Energy source in the turbulent energy equation given below,

{
alphat = turbulence->nut()/Prt;
alphat.correctBoundaryConditions();

volScalarField alphaEff
("alphaEff",
turbulence->nu()/Pr + turbulence->nut()/Prt
);

    fvScalarMatrix TEqn
        (
            fvm::ddt(T)
            + fvm::div(phi, T)
            - fvm::laplacian(alphaEff, T)
            - q/Cp1
        );

   	TEqn.relax();
	
        TEqn.solve();
}

so if you check the dimension of the equation it is coming (kelvin/sec) , my source q is specified in terms of (j/sec) hence i divided it with the heat capacity Cp1(j/k) of my source and specified correct dimensions in the transport property as well >my solver is running properly if i remove this source term but when i am trying to run with this source term its giving me the following error.

Different dimensions for =
     dimensions : [0 2 -1 0 0 0 0] = [0 4 -3 0 0 0 0]


    From function dimensionSet::operator=(const dimensionSet&) const
    in file dimensionSet/dimensionSet.C at line 165.

i am not able to track this bug in the solver if somebody can help me on this i will be thankful

Regards
Himanshu

[Antimony]

Hi,

Are you sure Cp1 is not the specific heat capacity? If it is specific heat capacity then the units would be J/(kg.K).

Also, isn't the unit of q W/m2?

Regards,

Antimony

[himanshu28]

Originally Posted by Antimony

Hi,

Are you sure Cp1 is not the specific heat capacity? If it is specific heat capacity then the units would be J/(kg.K).

Also, isn't the unit of q W/m2?

Regards,

Antimony

i referred the following article http://www.idurun.com/?p=556 to fit turbulent energy in my case. if you see the right hand side of the equation the dimensions you will find as (K/sec) and to justify the dimension i used q in (j/sec) and Cp1 in terms of heat capacity not specific heat capacity. so q u can assume here is in Watt only.