From 13fb959f0d108f4a821caaa9713c345d7fbc4f30 Mon Sep 17 00:00:00 2001
From: andy <a.heather@opencfd.co.uk>
Date: Wed, 17 Jun 2009 13:08:33 +0100
Subject: [PATCH] lagrangian solver updates
---
applications/solvers/combustion/coalChemistryFoam/hEqn.H | 2 +-
applications/solvers/combustion/dieselEngineFoam/YEqn.H | 5 ++---
applications/solvers/lagrangian/reactingParcelFoam/YEqn.H | 3 ++-
3 files changed, 5 insertions(+), 5 deletions(-)
diff --git a/applications/solvers/combustion/coalChemistryFoam/hEqn.H b/applications/solvers/combustion/coalChemistryFoam/hEqn.H
index 434177e136e..2eaf6315cfb 100644
--- a/applications/solvers/combustion/coalChemistryFoam/hEqn.H
+++ b/applications/solvers/combustion/coalChemistryFoam/hEqn.H
@@ -2,7 +2,7 @@
fvScalarMatrix hEqn
(
fvm::ddt(rho, h)
- + fvm::div(phi, h)
+ + mvConvection->fvmDiv(phi, h)
- fvm::laplacian(turbulence->alphaEff(), h)
==
DpDt
diff --git a/applications/solvers/combustion/dieselEngineFoam/YEqn.H b/applications/solvers/combustion/dieselEngineFoam/YEqn.H
index 2f74884bbde..9d46c0d6111 100644
--- a/applications/solvers/combustion/dieselEngineFoam/YEqn.H
+++ b/applications/solvers/combustion/dieselEngineFoam/YEqn.H
@@ -14,7 +14,7 @@ tmp<fv::convectionScheme<scalar> > mvConvection
label inertIndex = -1;
volScalarField Yt = 0.0*Y[0];
- for(label i=0; i<Y.size(); i++)
+ for (label i=0; i<Y.size(); i++)
{
if (Y[i].name() != inertSpecie)
{
@@ -39,8 +39,7 @@ tmp<fv::convectionScheme<scalar> > mvConvection
inertIndex = i;
}
}
-
+
Y[inertIndex] = scalar(1) - Yt;
Y[inertIndex].max(0.0);
-
}
diff --git a/applications/solvers/lagrangian/reactingParcelFoam/YEqn.H b/applications/solvers/lagrangian/reactingParcelFoam/YEqn.H
index c7ade57fbea..5cf44cb09fe 100644
--- a/applications/solvers/lagrangian/reactingParcelFoam/YEqn.H
+++ b/applications/solvers/lagrangian/reactingParcelFoam/YEqn.H
@@ -26,7 +26,8 @@ tmp<fv::convectionScheme<scalar> > mvConvection
- fvm::laplacian(turbulence->muEff(), Yi)
==
parcels.Srho(i)
- + kappa*chemistry.RR(i)().dimensionedInternalField()
+ + kappa*chemistry.RR(i)().dimensionedInternalField(),
+ mesh.solver("Yi")
);
Yi.max(0.0);
--
GitLab