zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

Added limitWith scheme. Currently this evaluated the limiter twice if the scheme is corrected.

Browse Source
This commit is contained in:
henry
2008-09-13 11:11:19 +01:00
parent 165bac9745
commit 78352100f2
2 changed files with 198 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

37
src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/limitWith/limitWith.C Normal file
View File

@ -0,0 +1,37 @@
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 1991-2008 OpenCFD Ltd.
\\/ M anipulation |
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
OpenFOAM is free software; you can redistribute it and/or modify it
under the terms of the GNU General Public License as published by the
Free Software Foundation; either version 2 of the License, or (at your
option) any later version.
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
for more details.
You should have received a copy of the GNU General Public License
along with OpenFOAM; if not, write to the Free Software Foundation,
Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
\*---------------------------------------------------------------------------*/
#include "fvMesh.H"
#include "limitWith.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
namespace Foam
{
makeSurfaceInterpolationScheme(limitWith)
}
// ************************************************************************* //

161
src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/limitWith/limitWith.H Normal file
View File

@ -0,0 +1,161 @@
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 1991-2008 OpenCFD Ltd.
\\/ M anipulation |
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
OpenFOAM is free software; you can redistribute it and/or modify it
under the terms of the GNU General Public License as published by the
Free Software Foundation; either version 2 of the License, or (at your
option) any later version.
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
for more details.
You should have received a copy of the GNU General Public License
along with OpenFOAM; if not, write to the Free Software Foundation,
Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
Class
Foam::limitWith
Description
limitWith differencing scheme limits the specified scheme with the
specified limiter.
SourceFiles
limitWith.C
\*---------------------------------------------------------------------------*/
#ifndef limitWith_H
#define limitWith_H
#include "surfaceInterpolationScheme.H"
#include "limitedSurfaceInterpolationScheme.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
namespace Foam
{
/*---------------------------------------------------------------------------*\
Class limitWith Declaration
\*---------------------------------------------------------------------------*/
template<class Type>
class limitWith
:
public surfaceInterpolationScheme<Type>
{
// Private Member Functions
//- Interpolation scheme
tmp<surfaceInterpolationScheme<Type> > tInterp_;
//- Limiter
tmp<limitedSurfaceInterpolationScheme<Type> > tLimiter_;
//- Disallow default bitwise copy construct
limitWith(const limitWith&);
//- Disallow default bitwise assignment
void operator=(const limitWith&);
public:
//- Runtime type information
TypeName("limitWith");
// Constructors
//- Construct from mesh and Istream.
// The name of the flux field is read from the Istream and looked-up
// from the mesh objectRegistry
limitWith
(
const fvMesh& mesh,
Istream& is
)
:
surfaceInterpolationScheme<Type>(mesh),
tInterp_
(
surfaceInterpolationScheme<Type>::New(mesh, is)
),
tLimiter_
(
limitedSurfaceInterpolationScheme<Type>::New(mesh, is)
)
{}
//- Construct from mesh, faceFlux and Istream
limitWith
(
const fvMesh& mesh,
const surfaceScalarField& faceFlux,
Istream& is
)
:
surfaceInterpolationScheme<Type>(mesh),
tInterp_
(
surfaceInterpolationScheme<Type>::New(mesh, faceFlux, is)
),
tLimiter_
(
limitedSurfaceInterpolationScheme<Type>::New(mesh, faceFlux, is)
)
{}
// Member Functions
//- Return the interpolation weighting factors
virtual tmp<surfaceScalarField> weights
(
const GeometricField<Type, fvPatchField, volMesh>& vf
) const
{
return tLimiter_().weights
(
vf,
tInterp_().weights(vf),
tLimiter_().limiter(vf)
);
}
//- Return true if this scheme uses an explicit correction
virtual bool corrected() const
{
return tInterp_().corrected();
}
//- Return the explicit correction to the face-interpolate
// for the given field
virtual tmp<GeometricField<Type, fvsPatchField, surfaceMesh> >
correction(const GeometricField<Type, fvPatchField, volMesh>& vf) const
{
return tLimiter_().limiter(vf)*tInterp_().correction(vf);
}
};
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
} // End namespace Foam
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
#endif
// ************************************************************************* //