From a86b1ce4a37e543d5cfb64ceccf0b581c761ab2b Mon Sep 17 00:00:00 2001
From: henry <Henry Weller h.weller@opencfd.co.uk>
Date: Thu, 10 Jul 2008 15:53:10 +0100
Subject: [PATCH] Corrected for SP.
---
.../solvers/multiphase/compressibleLesInterFoam/alphaEqns.H | 2 +-
.../multiphase/compressibleLesInterFoam/createFields.H | 2 +-
.../solvers/multiphase/compressibleLesInterFoam/pEqn.H | 5 ++++-
3 files changed, 6 insertions(+), 3 deletions(-)
diff --git a/applications/solvers/multiphase/compressibleLesInterFoam/alphaEqns.H b/applications/solvers/multiphase/compressibleLesInterFoam/alphaEqns.H
index 12eea31127e..819cd0f538b 100644
--- a/applications/solvers/multiphase/compressibleLesInterFoam/alphaEqns.H
+++ b/applications/solvers/multiphase/compressibleLesInterFoam/alphaEqns.H
@@ -28,7 +28,7 @@
),
// Divergence term is handled explicitly to be
// consistent with the explicit transport solution
- divU*min(alpha1, 1.0)
+ divU*min(alpha1, scalar(1))
);
forAll(dgdt, celli)
diff --git a/applications/solvers/multiphase/compressibleLesInterFoam/createFields.H b/applications/solvers/multiphase/compressibleLesInterFoam/createFields.H
index c6dd6c0f48a..6fa5049de47 100644
--- a/applications/solvers/multiphase/compressibleLesInterFoam/createFields.H
+++ b/applications/solvers/multiphase/compressibleLesInterFoam/createFields.H
@@ -140,7 +140,7 @@
);
volScalarField dgdt =
- pos(alpha2)*fvc::div(phi)/max(alpha2, 0.0001);
+ pos(alpha2)*fvc::div(phi)/max(alpha2, scalar(0.0001));
// Construct interface from alpha1 distribution
interfaceProperties interface(alpha1, U, twoPhaseProperties);
diff --git a/applications/solvers/multiphase/compressibleLesInterFoam/pEqn.H b/applications/solvers/multiphase/compressibleLesInterFoam/pEqn.H
index b794d4df2d5..ebf24498ade 100644
--- a/applications/solvers/multiphase/compressibleLesInterFoam/pEqn.H
+++ b/applications/solvers/multiphase/compressibleLesInterFoam/pEqn.H
@@ -40,7 +40,10 @@
solve
(
- (max(alpha1, 0.0)*(psi1/rho1) + max(alpha2, 0.0)*(psi2/rho2))
+ (
+ max(alpha1, scalar(0))*(psi1/rho1)
+ + max(alpha2, scalar(0))*(psi2/rho2)
+ )
*pdEqnComp()
+ pdEqnIncomp
);
--
GitLab