mirror of
https://develop.openfoam.com/Development/openfoam.git
synced 2025-11-28 03:28:01 +00:00
ENH: upadtes to semi-implicit cloud source calcs
This commit is contained in:
@ -202,7 +202,7 @@ const Foam::vector Foam::KinematicParcel<ParcelType>::calcVelocity
|
|||||||
const vector abp = (Feff.Sp()*Uc_ + (Feff.Su() + Su))/mass;
|
const vector abp = (Feff.Sp()*Uc_ + (Feff.Su() + Su))/mass;
|
||||||
const scalar bp = Feff.Sp()/mass;
|
const scalar bp = Feff.Sp()/mass;
|
||||||
|
|
||||||
Spu = Feff.Sp();
|
Spu = Feff.Sp()*dt/td.cloud().solution().deltaT();
|
||||||
|
|
||||||
IntegrationScheme<vector>::integrationResult Ures =
|
IntegrationScheme<vector>::integrationResult Ures =
|
||||||
td.cloud().UIntegrator().integrate(U, dt, abp, bp);
|
td.cloud().UIntegrator().integrate(U, dt, abp, bp);
|
||||||
|
|||||||
@ -339,7 +339,7 @@ Foam::scalar Foam::ThermoParcel<ParcelType>::calcHeatTransfer
|
|||||||
|
|
||||||
dhsTrans += dt*htc*As*(0.5*(T + Tnew) - Tc_);
|
dhsTrans += dt*htc*As*(0.5*(T + Tnew) - Tc_);
|
||||||
|
|
||||||
Cuh = bp;
|
Cuh = bp*dt/td.cloud().solution().deltaT();
|
||||||
|
|
||||||
return Tnew;
|
return Tnew;
|
||||||
}
|
}
|
||||||
|
|||||||
Reference in New Issue
Block a user