research/OpenFOAM-2.4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

ENH: surface film - reverted laminar model to be equivalent to old surface shear model

Browse source
This commit is contained in:
andy 2014-03-28 16:20:40 +00:00
parent e9fdb553da
commit ca71a67b28
1 changed files with 3 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

7
src/regionModels/surfaceFilmModels/submodels/kinematic/filmTurbulenceModel/laminar/laminar.C
View file

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -132,19 +132,18 @@ tmp<fvVectorMatrix> laminar::Su(volVectorField& U) const
// local references to film fields // local references to film fields
const volScalarField& mu = film.mu(); const volScalarField& mu = film.mu();
const volVectorField& Uw = film.Uw(); const volVectorField& Uw = film.Uw();
const volVectorField& Us = film.Us();
const volScalarField& delta = film.delta(); const volScalarField& delta = film.delta();
const volVectorField& Up = film.UPrimary(); const volVectorField& Up = film.UPrimary();
const volScalarField& rhop = film.rhoPrimary(); const volScalarField& rhop = film.rhoPrimary();
// employ simple coeff-based model // 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()))); volScalarField Cw("Cw", mu/(0.3333*(delta + film.deltaSmall())));
Cw.min(5000.0); Cw.min(5000.0);
return 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 - fvm::Sp(Cw, U) + Cw*Uw // wall contribution
); );
} }