Commit eecbe110 authored by Henry's avatar Henry
Browse files

CentredFitSnGrad: Change deltaCoeffs -> nonOrthDeltaCoeffs for very non-orthogonal meshes

parent da535a44
......@@ -25,9 +25,7 @@ License
#include "CentredFitSnGradData.H"
#include "surfaceFields.H"
#include "volFields.H"
#include "SVD.H"
#include "syncTools.H"
#include "extendedCentredCellToFaceStencil.H"
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
......@@ -226,7 +224,7 @@ void Foam::CentredFitSnGradData<Polynomial>::calcFit()
// find the fit coefficients for every face in the mesh
const surfaceScalarField& w = mesh.surfaceInterpolation::weights();
const surfaceScalarField& dC = mesh.deltaCoeffs();
const surfaceScalarField& dC = mesh.nonOrthDeltaCoeffs();
for (label facei = 0; facei < mesh.nInternalFaces(); facei++)
{
......
......@@ -25,7 +25,7 @@ Class
Foam::CentredFitSnGradData
Description
Data for the quadratic fit correction interpolation scheme
Data for centred fit snGrad schemes
SourceFiles
CentredFitSnGradData.C
......@@ -64,6 +64,7 @@ class CentredFitSnGradData
// values of the stencil to obtain the gradient for each direction
List<scalarList> coeffs_;
public:
TypeName("CentredFitSnGradData");
......
......@@ -100,7 +100,7 @@ public:
const GeometricField<Type, fvPatchField, volMesh>&
) const
{
return this->mesh().deltaCoeffs();
return this->mesh().nonOrthDeltaCoeffs();
}
//- Return true if this scheme uses an explicit correction
......@@ -124,18 +124,18 @@ public:
);
const CentredFitSnGradData<Polynomial>& cfd =
CentredFitSnGradData<Polynomial>::New
(
mesh,
stencil,
linearLimitFactor_,
centralWeight_
);
const List<scalarList>& f = cfd.coeffs();
CentredFitSnGradData<Polynomial>::New
(
mesh,
stencil,
linearLimitFactor_,
centralWeight_
);
tmp<GeometricField<Type, fvsPatchField, surfaceMesh> > sft
= stencil.weightedSum(vf, f);
(
stencil.weightedSum(vf, cfd.coeffs())
);
sft().dimensions() /= dimLength;
......@@ -143,6 +143,7 @@ public:
}
};
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
} // End namespace fv
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment