From 20dcad4e57cd6c6b071c0cc66bffe2cb73beeccb Mon Sep 17 00:00:00 2001
From: henry <Henry Weller h.weller@opencfd.co.uk>
Date: Fri, 29 May 2009 11:03:20 +0100
Subject: [PATCH] Added the reciprocal of the turbulent Schmidt number.
---
.../solvers/multiphase/twoLiquidMixingFoam/alphaEqn.H | 6 +++++-
.../solvers/multiphase/twoLiquidMixingFoam/createFields.H | 3 +++
2 files changed, 8 insertions(+), 1 deletion(-)
diff --git a/applications/solvers/multiphase/twoLiquidMixingFoam/alphaEqn.H b/applications/solvers/multiphase/twoLiquidMixingFoam/alphaEqn.H
index c27a3e752e2..d152baba0bb 100644
--- a/applications/solvers/multiphase/twoLiquidMixingFoam/alphaEqn.H
+++ b/applications/solvers/multiphase/twoLiquidMixingFoam/alphaEqn.H
@@ -3,7 +3,11 @@
(
fvm::ddt(alpha1)
+ fvm::div(phi, alpha1)
- - fvm::laplacian(turbulence->nut() + Dab, alpha1)
+ - fvm::laplacian
+ (
+ Dab + alphatab*turbulence->nut(), alpha1,
+ "laplacian(Dab,alpha1)"
+ )
);
alpha1Eqn.solve();
diff --git a/applications/solvers/multiphase/twoLiquidMixingFoam/createFields.H b/applications/solvers/multiphase/twoLiquidMixingFoam/createFields.H
index a38135dc19a..3b5e064372b 100644
--- a/applications/solvers/multiphase/twoLiquidMixingFoam/createFields.H
+++ b/applications/solvers/multiphase/twoLiquidMixingFoam/createFields.H
@@ -50,6 +50,9 @@
dimensionedScalar Dab(twoPhaseProperties.lookup("Dab"));
+ // Read the reciprocal of the turbulent Schmidt number
+ dimensionedScalar alphatab(twoPhaseProperties.lookup("alphatab"));
+
// Need to store rho for ddt(rho, U)
volScalarField rho("rho", alpha1*rho1 + (scalar(1) - alpha1)*rho2);
rho.oldTime();
--
GitLab