Multicut Piecewise-Linear Interface Calculation (MPLIC) corrected scheme is a surface interpolation scheme for flux calculation in advection of a bounded variable, e.g. phase fraction and for interface capturing in the volume of fluid (VoF) method. More...

Inheritance diagram for MPLIC:

[legend]

Collaboration diagram for MPLIC:

[legend]

Public Member Functions
	TypeName ("MPLIC")
	Runtime type information. More...

	MPLIC (const fvMesh &mesh, const surfaceScalarField &faceFlux, Istream &is)
	Construct from faceFlux and Istream. More...

virtual tmp< surfaceScalarField >	weights (const GeometricField< scalar, fvPatchField, volMesh > &) const
	Return the interpolation weighting factors. More...

virtual tmp< surfaceScalarField >	interpolate (const GeometricField< scalar, fvPatchField, volMesh > &vf) const
	Return the face-interpolate of the given cell field. More...

void	operator= (const MPLIC &)=delete
	Disallow default bitwise assignment. More...

Public Member Functions inherited from surfaceInterpolationScheme< scalar >
	TypeName ("surfaceInterpolationScheme")
	Runtime type information. More...

	declareRunTimeSelectionTable (tmp, surfaceInterpolationScheme, Mesh,(const fvMesh &mesh, Istream &schemeData),(mesh, schemeData))

	declareRunTimeSelectionTable (tmp, surfaceInterpolationScheme, MeshFlux,(const fvMesh &mesh, const surfaceScalarField &faceFlux, Istream &schemeData),(mesh, faceFlux, schemeData))

	surfaceInterpolationScheme (const fvMesh &mesh)
	Construct from mesh. More...

	surfaceInterpolationScheme (const surfaceInterpolationScheme &)
	Disallow default bitwise copy construction. More...

virtual	~surfaceInterpolationScheme ()
	Destructor. More...

const fvMesh &	mesh () const
	Return mesh reference. More...

tmp< GeometricField< scalar, fvsPatchField, surfaceMesh > >	interpolate (const tmp< GeometricField< scalar, fvPatchField, volMesh >> &) const
	Return the face-interpolate of the given tmp cell field. More...

virtual tmp< GeometricField< typename innerProduct< vector, scalar >::type, fvsPatchField, surfaceMesh > >	dotInterpolate (const surfaceVectorField &Sf, const GeometricField< scalar, fvPatchField, volMesh > &vf) const
	Return the face-interpolate of the given cell field. More...

tmp< GeometricField< typename innerProduct< vector, scalar >::type, fvsPatchField, surfaceMesh > >	dotInterpolate (const surfaceVectorField &Sf, const tmp< GeometricField< scalar, fvPatchField, volMesh >> &) const
	Return the face-interpolate of the given tmp cell field. More...

virtual bool	corrected () const
	Return true if this scheme uses an explicit correction. More...

virtual tmp< GeometricField< scalar, fvsPatchField, surfaceMesh > >	correction (const GeometricField< scalar, fvPatchField, volMesh > &) const
	Return the explicit correction to the face-interpolate. More...

void	operator= (const surfaceInterpolationScheme &)=delete
	Disallow default bitwise assignment. More...

Public Member Functions inherited from refCount
int	count () const
	Return the current reference count. More...

bool	unique () const
	Return true if the reference count is zero. More...

void	operator++ ()
	Increment the reference count. More...

void	operator++ (int)
	Increment the reference count. More...

void	operator-- ()
	Decrement the reference count. More...

void	operator-- (int)
	Decrement the reference count. More...

Protected Member Functions
void	setCellAlphaf (const label celli, const scalarField &phi, scalarField &alphaf, boolList &correctedFaces, const DynamicList< scalar > &cellAlphaf, const fvMesh &mesh) const
	Set alphaPhi for the faces of the given cell. More...

tmp< surfaceScalarField >	surfaceAlpha (const volScalarField &alpha, const surfaceScalarField &phi, scalarField &spicedTvff, const bool unweighted, const scalar tol, const bool isMPLIC=true) const
	Return alpha interpolation. More...

Protected Member Functions inherited from refCount
	refCount ()
	Construct null initializing count to 0. More...

Protected Attributes
const surfaceScalarField &	phi_

Additional Inherited Members
Static Public Member Functions inherited from surfaceInterpolationScheme< scalar >
static tmp< surfaceInterpolationScheme< scalar > >	New (const fvMesh &mesh, Istream &schemeData)
	Return new tmp interpolation scheme. More...

static tmp< surfaceInterpolationScheme< scalar > >	New (const fvMesh &mesh, const surfaceScalarField &faceFlux, Istream &schemeData)
	Return new tmp interpolation scheme. More...

static tmp< GeometricField< scalar, fvsPatchField, surfaceMesh > >	interpolate (const GeometricField< scalar, fvPatchField, volMesh > &, const tmp< surfaceScalarField > &, const tmp< surfaceScalarField > &)
	Return the face-interpolate of the given cell field. More...

static tmp< GeometricField< scalar, fvsPatchField, surfaceMesh > >	interpolate (const GeometricField< scalar, fvPatchField, volMesh > &, const tmp< surfaceScalarField > &)
	Return the face-interpolate of the given cell field. More...

static tmp< GeometricField< typename innerProduct< typename SFType::value_type, scalar >::type, fvsPatchField, surfaceMesh > >	dotInterpolate (const SFType &Sf, const GeometricField< scalar, fvPatchField, volMesh > &vf, const tmp< surfaceScalarField > &tlambdas)
	Return the face-interpolate of the given cell field. More...

Detailed Description

Multicut Piecewise-Linear Interface Calculation (MPLIC) corrected scheme is a surface interpolation scheme for flux calculation in advection of a bounded variable, e.g. phase fraction and for interface capturing in the volume of fluid (VoF) method.

The interface is represented by multiple cuts which split each cell to match the volume fraction of the phase in the cell. The cut planes are oriented according to the point field of the local phase fraction. The phase fraction at each cell face - the interpolated value - is then calculated from the face area on either side of the cuts.

Three progressively more complex algorithms are used to ensure the cell volume fraction is accurately reproduced:

single cut: cuts all the cell faces regardless the order
multi cut: topological face-edge-face walk which can split cell into multiple sub-volumes
tetrahedron cut: decomposes cell into tetrahedrons which are cut

Example:

divSchemes
{
    .
    .
    div(phi,alpha)      Gauss MPLIC;
    .
    .
}

See also: Foam::MPLICU Foam::PLIC Foam::PLICU Foam::interfaceCompression

Source files

Definition at line 83 of file MPLIC.H.

Constructor & Destructor Documentation

◆ MPLIC()

MPLIC	(	const fvMesh &	mesh,
		const surfaceScalarField &	faceFlux,
		Istream &	is
	)

inline

Construct from faceFlux and Istream.

Definition at line 129 of file MPLIC.H.

References MPLIC::weights().

Here is the call graph for this function:

Member Function Documentation

◆ setCellAlphaf()

void setCellAlphaf	(	const label	celli,
		const scalarField &	phi,
		scalarField &	alphaf,
		boolList &	correctedFaces,
		const DynamicList< scalar > &	cellAlphaf,
		const fvMesh &	mesh
	)		const

protected

Set alphaPhi for the faces of the given cell.

Definition at line 47 of file MPLIC.C.

References primitiveMesh::cells(), polyMesh::faceOwner(), forAll, Foam::sign(), and MPLIC::surfaceAlpha().

Here is the call graph for this function:

◆ surfaceAlpha()

Foam::tmp< Foam::surfaceScalarField > surfaceAlpha	(	const volScalarField &	alpha,
		const surfaceScalarField &	phi,
		scalarField &	spicedTvff,
		const bool	unweighted,
		const scalar	tol,
		const bool	isMPLIC = `true`
	)		const