|
[Sponsors] |
March 6, 2015, 06:12 |
heat transfer equation of parcel
|
#1 |
Member
Neilson Whit
Join Date: Aug 2011
Posts: 74
Rep Power: 14 |
Dear All
Could you please write mathematical formulation of heat transfer equation for a droplet? It is interpreted as follows in src/lagrangian/intermediate/parcels/Templates/ThermoParcel.C Code:
template<class ParcelType> template<class TrackData> Foam::scalar Foam::ThermoParcel<ParcelType>::calcHeatTransfer ( TrackData& td, const scalar dt, const label cellI, const scalar Re, const scalar Pr, const scalar kappa, const scalar NCpW, const scalar Sh, scalar& dhsTrans, scalar& Sph ) { if (!td.cloud().heatTransfer().active()) { return T_; } const scalar d = this->d(); const scalar rho = this->rho(); // Calc heat transfer coefficient scalar htc = td.cloud().heatTransfer().htc(d, Re, Pr, kappa, NCpW); if (mag(htc) < ROOTVSMALL && !td.cloud().radiation()) { return max ( T_ + dt*Sh/(this->volume(d)*rho*Cp_), td.cloud().constProps().TMin() ); } htc = max(htc, ROOTVSMALL); const scalar As = this->areaS(d); scalar ap = Tc_ + Sh/As/htc; scalar bp = 6.0*(Sh/As + htc*(Tc_ - T_)); if (td.cloud().radiation()) { tetIndices tetIs = this->currentTetIndices(); const scalar Gc = td.GInterp().interpolate(this->position(), tetIs); const scalar sigma = physicoChemical::sigma.value(); const scalar epsilon = td.cloud().constProps().epsilon0(); ap = (ap + epsilon*Gc/(4.0*htc))/(1.0 + epsilon*sigma*pow3(T_)/htc); bp += 6.0*(epsilon*(Gc/4.0 - sigma*pow4(T_))); } bp /= rho*d*Cp_*(ap - T_) + ROOTVSMALL; // Integrate to find the new parcel temperature IntegrationScheme<scalar>::integrationResult Tres = td.cloud().TIntegrator().integrate(T_, dt, ap*bp, bp); scalar Tnew = min ( max ( Tres.value(), td.cloud().constProps().TMin() ), td.cloud().constProps().TMax() ); Sph = dt*htc*As; dhsTrans += Sph*(Tres.average() - Tc_); return Tnew; } |
|
June 23, 2016, 04:24 |
|
#2 |
New Member
Amir
Join Date: Feb 2012
Posts: 12
Rep Power: 14 |
Dear Neilson,
Have you found any mathematics formulation for parcel energy balance or any reference that the model is written based on? Regards, |
|
November 24, 2016, 09:49 |
|
#3 |
Member
Chris Cloney
Join Date: Jun 2016
Location: Halifax, Canada
Posts: 62
Rep Power: 9 |
One note here is that Spu, Sh, Sph, and dhstrans are all passed in as zeros to the calcHeatTransfer function. That starts to get you close to something that looks like the textbook equations, but I am still working on it.
|
|
December 9, 2016, 08:54 |
|
#4 |
Member
Chris Cloney
Join Date: Jun 2016
Location: Halifax, Canada
Posts: 62
Rep Power: 9 |
As you look more you learn more.
Sh, Sph, ext are only passed in as zero when the parcels are defined as a thermoCloud. ReactingCloud and ReactingMultiphaseCloud may pass in nonzero values. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Question about heat transfer coefficient setting for CFX | Anna Tian | CFX | 1 | June 16, 2013 06:28 |
energy equation show abnormal behavior in conjugate heat transfer | Muhammad imtiaz | FLUENT | 0 | June 13, 2013 02:00 |
radiation heat transfer | GeHa | FLUENT | 1 | September 5, 2012 14:56 |
error message | cuteapathy | CFX | 14 | March 20, 2012 06:45 |
Convective / Conductive Heat Transfer in Hypersonic flows | enigma | Main CFD Forum | 2 | November 1, 2009 22:53 |