[dvcauwe]

Hello Foamers, I have been breaking my head over this for the past 2 days, so I would really appreciate it if I could get some feedback on this.

I'm simulating flow through a pipe using a modified version of rhoPimpleFoam with cyclic boundaries at the in- and outlet of my domain while the flow is heated with a constant temperature gradient imposed on the tube wall. The problem however is that when I impose a gradient that is twice as big, the rate at which my gas is heating is only slightly higher and takes a considerable amount of time (~1s) to reach the double value of my initial deltaT/time with half the gradient imposed.

I understand that it takes a while to develop a temperature profile but still energy should be conserved, right? Performing additional subiterations offers a slight improvement but still doesn't balance out the energy equation. I thought a good way to monitor the heat going into the volume would be this function:

const fvPatchField<scalar>& heatFluxPatchField = T.boundaryField()[patchI];

// Print patch total heat transfer rate
const tmp<volScalarField> tkappa = thermo.kappa();
const scalarField& kappa = tkappa().boundaryField()[patchI];

scalar patchHeatRate = gSum
(
	kappa * heatFluxPatchField.snGrad() *
	heatFluxPatchField.patch().magSf()
);

, which shows me the desired total heat input but unfortunately this energy is not actually absorbed by the gas.

The energy equation is standard I'd say:

    volScalarField& he = thermo.he();

    fvScalarMatrix EEqn
    (
        fvm::ddt(rho, he) + fvm::div(phi, he)
      + fvc::ddt(rho, K) + fvc::div(phi, K)
      - fvm::laplacian(turbulence->alphaEff(), he)
     ==
        fvOptions(rho, he)
    );
    EEqn.relax();
    fvOptions.constrain(EEqn);
    EEqn.solve();
    fvOptions.correct(he);
    thermo.correct();

My only remaining guess at this moment is that setting a constant gradient for temperature simply does not guarantee that amount of energy going into the fluid during this time interval but that would imply setting an energy flux for...h? Then I'm not entirely sure how to use this field again in my thermo calculations...

Anyway, thanks to anyone that might be able to shed some light on this!

Best regards,
David