From 432f72a88a0810d68fbf9eee74719dcaa721cd31 Mon Sep 17 00:00:00 2001 From: andy Date: Wed, 2 Mar 2011 15:41:24 +0000 Subject: [PATCH] ENH: upadtes to semi-implicit cloud source calcs --- .../parcels/Templates/KinematicParcel/KinematicParcel.C | 2 +- .../intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C index 9a83e69908..0a3e05d200 100644 --- a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C +++ b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C @@ -202,7 +202,7 @@ const Foam::vector Foam::KinematicParcel::calcVelocity const vector abp = (Feff.Sp()*Uc_ + (Feff.Su() + Su))/mass; const scalar bp = Feff.Sp()/mass; - Spu = Feff.Sp(); + Spu = Feff.Sp()*dt/td.cloud().solution().deltaT(); IntegrationScheme::integrationResult Ures = td.cloud().UIntegrator().integrate(U, dt, abp, bp); diff --git a/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C b/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C index 064f945ad1..9df054c2d0 100644 --- a/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C +++ b/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C @@ -339,7 +339,7 @@ Foam::scalar Foam::ThermoParcel::calcHeatTransfer dhsTrans += dt*htc*As*(0.5*(T + Tnew) - Tc_); - Cuh = bp; + Cuh = bp*dt/td.cloud().solution().deltaT(); return Tnew; }