From 4dda7cbea2eddcaebb0a18cdb3c6f5f507948bc1 Mon Sep 17 00:00:00 2001 From: danielque Date: Tue, 12 Oct 2021 10:37:36 +0200 Subject: [PATCH] fix heat transfer for general heat capacities same as for cfdemSolverRhoPimple in commit 0fa1f023d62cd43898b81e04b6c6bd3529e8ce03 --- .../solvers/cfdemSolverRhoPimpleChem/EEqn.H | 16 +++++++++------- 1 file changed, 9 insertions(+), 7 deletions(-) diff --git a/applications/solvers/cfdemSolverRhoPimpleChem/EEqn.H b/applications/solvers/cfdemSolverRhoPimpleChem/EEqn.H index 9a5f10b8..efcf6fa1 100644 --- a/applications/solvers/cfdemSolverRhoPimpleChem/EEqn.H +++ b/applications/solvers/cfdemSolverRhoPimpleChem/EEqn.H @@ -8,9 +8,11 @@ particleCloud.energyCoefficients(QCoeff); Cpv = he.name() == "e" ? thermo.Cv() : thermo.Cp(); -// correct source for the thermodynamic reference temperature -// dimensionedScalar Tref("Tref", dimTemperature, T[0]-he[0]/(Cpv[0]+SMALL)); -// Qsource += QCoeff*Tref; +// For implict T terms in the energy/enthalpy transport equation, use +// (he_n+1 - he_n) / (T_n+1 - T_n) = Cpv to eliminate T_n+1 with he_n+1. +// This formula is valid for ideal gases with e=e(T) and h=h(T). For +// incompressible fluids, e=e(T) holds, too, but enthalpy would need correction +// terms accounting for pressure variations. fvScalarMatrix EEqn ( @@ -26,13 +28,13 @@ fvScalarMatrix EEqn ) : -dpdt ) - // net heat transfer from particles to fluid - Qsource + - QCoeff*T - fvm::Sp(QCoeff/Cpv, he) - // thermal conduction of the fluid with effective conductivity + + QCoeff/Cpv*he + - fvc::laplacian(voidfraction*thCond,T) - fvm::laplacian(voidfraction*thCond/Cpv,he) - // + particle-fluid energy transfer due to work - // + fluid energy dissipation due to shearing + + fvc::laplacian(voidfraction*thCond/Cpv,he) == // + combustion->Sh() fvOptions(rho, he)