From 4ae925db45316bff99e54670c38df920cb77f053 Mon Sep 17 00:00:00 2001
From: andy <a.heather@opencfd.co.uk>
Date: Mon, 10 Jan 2011 15:49:33 +0000
Subject: [PATCH] ENH: Updated particle surface thermo props
---
.../ReactingMultiphaseParcel.C | 19 ++++++------
.../Templates/ReactingParcel/ReactingParcel.C | 31 ++++++++++---------
.../Templates/ReactingParcel/ReactingParcel.H | 6 ++--
.../Templates/ThermoParcel/ThermoParcel.C | 12 +++----
.../Templates/ThermoParcel/ThermoParcel.H | 4 +--
5 files changed, 37 insertions(+), 35 deletions(-)
diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.C b/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.C
index 76a80c0493d..98e357950d8 100644
--- a/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.C
@@ -228,12 +228,12 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::calc
// Calc surface values
// ~~~~~~~~~~~~~~~~~~~
- scalar Ts, rhos, mus, Pr, kappa;
+ scalar Ts, rhos, mus, Prs, kappas;
ThermoParcel<ParcelType>::
- calcSurfaceValues(td, cellI, T0, Ts, rhos, mus, Pr, kappa);
+ calcSurfaceValues(td, cellI, T0, Ts, rhos, mus, Prs, kappas);
// Reynolds number
- scalar Re = this->Re(U0, d0, rhos, mus);
+ scalar Res = this->Re(U0, d0, rhos, mus);
// Sources
@@ -273,7 +273,7 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::calc
td,
dt,
cellI,
- Re,
+ Res,
Ts,
mus/rhos,
d0,
@@ -316,7 +316,8 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::calc
);
// Correct surface values due to emitted species
- this->correctSurfaceValues(td, cellI, Ts, Cs, rhos, mus, Pr, kappa);
+ this->correctSurfaceValues(td, cellI, Ts, Cs, rhos, mus, Prs, kappas);
+ Res = this->Re(U0, d0, rhos, mus);
// Surface reactions
@@ -374,9 +375,9 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::calc
td,
dt,
cellI,
- Re,
- Pr,
- kappa,
+ Res,
+ Prs,
+ kappas,
d0,
rho0,
T0,
@@ -399,7 +400,7 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::calc
td,
dt,
cellI,
- Re,
+ Res,
mus,
d0,
U0,
diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
index 74a931c16ac..2c05144d9d3 100644
--- a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
@@ -130,8 +130,8 @@ void Foam::ReactingParcel<ParcelType>::correctSurfaceValues
const scalarField& Cs,
scalar& rhos,
scalar& mus,
- scalar& Pr,
- scalar& kappa
+ scalar& Prs,
+ scalar& kappas
)
{
// No correction if total concentration of emitted species is small
@@ -177,7 +177,7 @@ void Foam::ReactingParcel<ParcelType>::correctSurfaceValues
rhos = 0;
mus = 0;
- kappa = 0;
+ kappas = 0;
scalar Cps = 0;
scalar sumYiSqrtW = 0;
scalar sumYiCbrtW = 0;
@@ -190,7 +190,7 @@ void Foam::ReactingParcel<ParcelType>::correctSurfaceValues
rhos += Xs[i]*W;
mus += Ys[i]*sqrtW*thermo.carrier().mu(i, T);
- kappa += Ys[i]*cbrtW*thermo.carrier().kappa(i, T);
+ kappas += Ys[i]*cbrtW*thermo.carrier().kappa(i, T);
Cps += Xs[i]*thermo.carrier().Cp(i, T);
sumYiSqrtW += Ys[i]*sqrtW;
@@ -199,8 +199,8 @@ void Foam::ReactingParcel<ParcelType>::correctSurfaceValues
rhos *= pc_/(specie::RR*T);
mus /= sumYiSqrtW;
- kappa /= sumYiCbrtW;
- Pr = Cps*mus/kappa;
+ kappas /= sumYiCbrtW;
+ Prs = Cps*mus/kappas;
}
@@ -252,11 +252,11 @@ void Foam::ReactingParcel<ParcelType>::calc
// Calc surface values
// ~~~~~~~~~~~~~~~~~~~
- scalar Ts, rhos, mus, Pr, kappa;
- this->calcSurfaceValues(td, cellI, T0, Ts, rhos, mus, Pr, kappa);
+ scalar Ts, rhos, mus, Prs, kappas;
+ this->calcSurfaceValues(td, cellI, T0, Ts, rhos, mus, Prs, kappas);
// Reynolds number
- scalar Re = this->Re(U0, d0, rhos, mus);
+ scalar Res = this->Re(U0, d0, rhos, mus);
// Sources
@@ -296,7 +296,7 @@ void Foam::ReactingParcel<ParcelType>::calc
td,
dt,
cellI,
- Re,
+ Res,
Ts,
mus/rhos,
d0,
@@ -313,7 +313,8 @@ void Foam::ReactingParcel<ParcelType>::calc
);
// Correct surface values due to emitted species
- correctSurfaceValues(td, cellI, Ts, Cs, rhos, mus, Pr, kappa);
+ correctSurfaceValues(td, cellI, Ts, Cs, rhos, mus, Prs, kappas);
+ Res = this->Re(U0, d0, rhos, mus);
// Update particle component mass and mass fractions
scalar mass1 = updateMassFraction(mass0, dMassPC, Y_);
@@ -330,9 +331,9 @@ void Foam::ReactingParcel<ParcelType>::calc
td,
dt,
cellI,
- Re,
- Pr,
- kappa,
+ Res,
+ Prs,
+ kappas,
d0,
rho0,
T0,
@@ -355,7 +356,7 @@ void Foam::ReactingParcel<ParcelType>::calc
td,
dt,
cellI,
- Re,
+ Res,
mus,
d0,
U0,
diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.H b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.H
index e24049dc8c6..d1eb1c58622 100644
--- a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.H
+++ b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.H
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 1991-2010 OpenCFD Ltd.
+ \\ / A nd | Copyright (C) 1991-2011 OpenCFD Ltd.
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -358,8 +358,8 @@ public:
const scalarField& Cs,
scalar& rhos,
scalar& mus,
- scalar& Pr,
- scalar& kappa
+ scalar& Prs,
+ scalar& kappas
);
//- Update parcel properties over the time interval
diff --git a/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C b/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C
index ba98abf666b..5769d77853f 100644
--- a/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 1991-2010 OpenCFD Ltd.
+ \\ / A nd | Copyright (C) 1991-2011 OpenCFD Ltd.
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -108,7 +108,7 @@ void Foam::ThermoParcel<ParcelType>::calcSurfaceValues
scalar& rhos,
scalar& mus,
scalar& Pr,
- scalar& kappa
+ scalar& kappas
) const
{
// Surface temperature using two thirds rule
@@ -144,7 +144,7 @@ void Foam::ThermoParcel<ParcelType>::calcSurfaceValues
mus = td.muInterp().interpolate(this->position(), tetIs)/TRatio;
Pr = td.cloud().constProps().Pr();
- kappa = Cpc_*mus/Pr;
+ kappas = Cpc_*mus/Pr;
}
@@ -170,8 +170,8 @@ void Foam::ThermoParcel<ParcelType>::calc
// Calc surface values
// ~~~~~~~~~~~~~~~~~~~
- scalar Ts, rhos, mus, Pr, kappa;
- calcSurfaceValues(td, cellI, T0, Ts, rhos, mus, Pr, kappa);
+ scalar Ts, rhos, mus, Pr, kappas;
+ calcSurfaceValues(td, cellI, T0, Ts, rhos, mus, Pr, kappas);
// Reynolds number
scalar Re = this->Re(U0, d0, rhos, mus);
@@ -209,7 +209,7 @@ void Foam::ThermoParcel<ParcelType>::calc
cellI,
Re,
Pr,
- kappa,
+ kappas,
d0,
rho0,
T0,
diff --git a/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.H b/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.H
index 0267e7451a1..8a954cb082c 100644
--- a/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.H
+++ b/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.H
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 1991-2010 OpenCFD Ltd.
+ \\ / A nd | Copyright (C) 1991-2011 OpenCFD Ltd.
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -381,7 +381,7 @@ public:
scalar& rhos,
scalar& mus,
scalar& Pr,
- scalar& kappa
+ scalar& kappas
) const;
//- Update parcel properties over the time interval
--
GitLab