[Lorddrinkalot]

Hey guys,

I could use some help setting up my Volumetric heat source. I am trying to model a uniform volumetric thermal source in the fluid.

For that, I am using fvOptions. I already wrote my option file and the solver also works; however, I get a warning "Source heatSource defined for field he but never used", and I can't find out why. I have tried different fvOptions.

In both cases, I get the same error, and I can't find out why.

I am using a solver named buoyantPimpleMhdFoam (combination to use MHD and Heat). Further down, I'll post the EEqn.H.

Thanks, Christoph

first one

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvOptions;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
options
{
heatSource
{
    type            scalarSemiImplicitSource;
    active          true;
    selectionMode   all;
	volumeMode      specific;

		sources
		{
            he           (19.4487 0); 
        }
    
}
}

second one I used

    HeatSource
    {
        type scalarSemiImplicitSource;
        active          true;
        selectionMode   all; 

        scalarSemiImplicitSourceCoeffs 
        {
            volumeMode absolute; //alt. specific
			selectionMode   all;
            injectionRateSuSp 
            {
                he (20 0);
            }
        }
   }

{
    volScalarField& he = thermo.he();

    fvScalarMatrix EEqn
    (
        fvm::ddt(rho, he) + fvm::div(phi, he)
      + fvc::ddt(rho, K) + fvc::div(phi, K)
      + (
            he.name() == "e"
          ? fvc::div
            (
                fvc::absolute(phi/fvc::interpolate(rho), U),
                p,
                "div(phiv,p)"
            )
          : -dpdt
        )
      - fvm::laplacian(turbulence->alphaEff(), he)
     ==
        rho*(U&g)
      + radiation->Sh(thermo, he)
      + fvOptions(rho, he)
    );

    EEqn.relax();

    fvOptions.constrain(EEqn);

    EEqn.solve();

    fvOptions.correct(he);

    thermo.correct();
    radiation->correct();
}

The warning says:
PIMPLE: iteration 1
--> FOAM Warning :
From virtual void Foam::fv:ption::checkApplied() const
in file cfdTools/general/fvOptions/fvOption.C at line 137
Source heatSource defined for field he but never used

[Yann]

Hello Christoph,

Here:

he           (19.4487 0);

Try replacing he by h or e depending on your thermophysical setup.
If you are using enthalpy it should be h, if you are using internal energy it should be e.

In the code the variable is named he because it can be both, but the simulation actually solves for h or e depending on the case settings, and your source must be defined for the variable which is solved.

I hope this will solve your problem.
Yann