From ca71a67b2815c0ed943f5e296b2a7616921e155e Mon Sep 17 00:00:00 2001 From: andy Date: Fri, 28 Mar 2014 16:20:40 +0000 Subject: [PATCH] ENH: surface film - reverted laminar model to be equivalent to old surface shear model --- .../kinematic/filmTurbulenceModel/laminar/laminar.C | 7 +++---- 1 file changed, 3 insertions(+), 4 deletions(-) diff --git a/src/regionModels/surfaceFilmModels/submodels/kinematic/filmTurbulenceModel/laminar/laminar.C b/src/regionModels/surfaceFilmModels/submodels/kinematic/filmTurbulenceModel/laminar/laminar.C index a64205df..6b015cd8 100644 --- a/src/regionModels/surfaceFilmModels/submodels/kinematic/filmTurbulenceModel/laminar/laminar.C +++ b/src/regionModels/surfaceFilmModels/submodels/kinematic/filmTurbulenceModel/laminar/laminar.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -132,19 +132,18 @@ tmp laminar::Su(volVectorField& U) const // local references to film fields const volScalarField& mu = film.mu(); const volVectorField& Uw = film.Uw(); - const volVectorField& Us = film.Us(); const volScalarField& delta = film.delta(); const volVectorField& Up = film.UPrimary(); const volScalarField& rhop = film.rhoPrimary(); // employ simple coeff-based model - volScalarField Cs("Cs", Cf_*rhop*mag(Up - Us)); + volScalarField Cs("Cs", Cf_*rhop*mag(Up - U)); volScalarField Cw("Cw", mu/(0.3333*(delta + film.deltaSmall()))); Cw.min(5000.0); return ( - - fvm::Sp(Cs, U) + Cs*Us // surface contribution + - fvm::Sp(Cs, U) + Cs*Up // surface contribution - fvm::Sp(Cw, U) + Cw*Uw // wall contribution ); }