Base class for passing Dirichlet EB boundary conditions into EBHelmholtzOp. More...

#include <CD_EBHelmholtzDirichletEBBC.H>

Inheritance diagram for EBHelmholtzDirichletEBBC:

Collaboration diagram for EBHelmholtzDirichletEBBC:

Public Member Functions
	EBHelmholtzDirichletEBBC ()
	Default constructor.

virtual	~EBHelmholtzDirichletEBBC ()
	Destructor.

virtual void	setOrder (const int a_order)
	Set BC order. More...

virtual void	setWeight (const int a_weight)
	Set stencil weights. More...

virtual void	setValue (const int a_value)
	Set the value on the EB. More...

virtual void	setValue (const std::function< Real(const RealVect &a_pos)> &a_value)
	Set the value on the EB. More...

virtual void	setDomainDropOrder (const int a_domainSize)
	Drop BC order if domain size is equal or below this.

virtual void	applyEBFlux (VoFIterator &a_vofit, EBCellFAB &a_Lphi, const EBCellFAB &a_phi, const BaseIVFAB< Real > &a_Bcoef, const DataIndex &a_dit, const Real &a_beta, const bool &a_homogeneousPhysBC) const override
	Apply the EB flux. This is the version that is called by EBHelmholtzOp. More...

Public Member Functions inherited from EBHelmholtzEBBC
	EBHelmholtzEBBC ()
	Default constructor.

virtual	~EBHelmholtzEBBC ()
	Destructor.

	EBHelmholtzEBBC (const EBHelmholtzEBBC &a_other)=delete
	Disallowed – don't see why you would need it.

EBHelmholtzEBBC &	operator= (const EBHelmholtzEBBC &a_other)=delete
	Disallowed - don't see why you would need it.

void	define (const Location::Cell a_dataLocation, const EBLevelGrid &a_eblg, const RealVect &a_probLo, const Real &a_dx, const int a_ghostCF)
	Define function that is called by EBHelmholtzOp. More...

const LayoutData< BaseIVFAB< VoFStencil > > &	getGradPhiStencils () const
	Get the stencil for computing the finite-volume approximation to kappa*Div(F). More...

Protected Member Functions
virtual void	define () override
	User define function.

virtual bool	getLeastSquaresStencil (std::pair< Real, VoFStencil > &a_stencil, const VolIndex &a_vof, const VofUtils::Neighborhood a_neighborhood, const DataIndex &a_dit, const int a_order) const
	Get stencil with specified order. More...

virtual bool	getJohansenStencil (std::pair< Real, VoFStencil > &a_stencil, const VolIndex &a_vof, const DataIndex &a_dit, const int a_order) const
	Get johansen stencil with specified order. Not used but exists for backward compatibility testing and debugging. More...

virtual bool	getChomboLsqStencil (std::pair< Real, VoFStencil > &a_stencil, const VolIndex &a_vof, const DataIndex &a_dit) const
	Get Chombo's least squares gradient stencil. Not used but kept for backward compatibility testing and debugging. More...

Protected Member Functions inherited from EBHelmholtzEBBC
bool	isStencilValidCF (const VoFStencil &a_stencil, const DataIndex &a_dit) const
	Check if stencil is valid. More...

RealVect	getBoundaryPosition (const VolIndex &a_vof, const DataIndex &a_dit) const
	Returns physical position at the boundary.

Real	applyStencil (const VoFStencil &a_stencil, const EBCellFAB &a_phi) const
	Apply stencil to data holder and return result.

Protected Attributes
bool	m_useConstant
	Use constant for BC.

bool	m_useFunction
	Use function for BC value.

int	m_order
	Stencil approximation order,.

int	m_weight
	Distance weighting factor for stencils.

int	m_domainDropOrder
	Special flag for dropping stencil order when domains become coarser than this.

Real	m_constantValue
	Value on EB.

std::function< Real(const RealVect &a_pos)>	m_functionValue
	Value on the EB.

LayoutData< BaseIVFAB< Real > >	m_boundaryWeights
	Boundary weights.

Protected Attributes inherited from EBHelmholtzEBBC
Location::Cell	m_dataLocation
	Data centering.

int	m_ghostCF
	Number of ghost cells that were filled across CF interface.

Real	m_dx
	Grid resolution.

RealVect	m_probLo
	Lower-left corner of computational domain.

EBLevelGrid	m_eblg
	Level grid.

LayoutData< BaseIVFAB< VoFStencil > >	m_gradPhiStencils
	Stencils for computing the flux on a single level. This is a single-level object.

Additional Inherited Members
Static Protected Attributes inherited from EBHelmholtzEBBC
constexpr static int	m_comp = 0
	Component that everything is defined for. Always have m_comp = 0.

constexpr static int	m_nComp = 1
	Number of components. Always have m_nComp = 1.

Detailed Description

Base class for passing Dirichlet EB boundary conditions into EBHelmholtzOp.

Member Function Documentation

◆ applyEBFlux()

void EBHelmholtzDirichletEBBC::applyEBFlux	(	VoFIterator &	a_vofit,
		EBCellFAB &	a_Lphi,
		const EBCellFAB &	a_phi,
		const BaseIVFAB< Real > &	a_Bcoef,
		const DataIndex &	a_dit,
		const Real &	a_beta,
		const bool &	a_homogeneousPhysBC
	)		const

overridevirtual

Apply the EB flux. This is the version that is called by EBHelmholtzOp.

Parameters

[in,out]	a_vofit	Iterator for cut-cells
[in]	a_Lphi	Operator kappa*L(phi)
[in]	a_phi	Phi, cell-centered.
[in]	a_Bcoef	EB-centered B-coefficient
[in]	a_dit	Data index
[in]	a_beta	Beta. Passed in from operator.
[in]	a_homogeneousPhysBC	Homogeneous BC or not.

Implements EBHelmholtzEBBC.

◆ getChomboLsqStencil()

bool EBHelmholtzDirichletEBBC::getChomboLsqStencil	(	std::pair< Real, VoFStencil > &	a_stencil,
		const VolIndex &	a_vof,
		const DataIndex &	a_dit
	)		const

protectedvirtual

Get Chombo's least squares gradient stencil. Not used but kept for backward compatibility testing and debugging.

Parameters

[out]	a_stencil	Weight and stencil
[out]	a_vof	Grid cell
[in]	a_dit	Data index
[in]	a_order	Attempted order for stencil. This is also the stencil radius.

Returns: If this returns true, a_stencil contains a valid stencil and weight.

◆ getJohansenStencil()

bool EBHelmholtzDirichletEBBC::getJohansenStencil	(	std::pair< Real, VoFStencil > &	a_stencil,
		const VolIndex &	a_vof,
		const DataIndex &	a_dit,
		const int	a_order
	)		const

protectedvirtual

Get johansen stencil with specified order. Not used but exists for backward compatibility testing and debugging.

Parameters

[out]	a_stencil	Weight and stencil
[out]	a_vof	Grid cell
[in]	a_dit	Data index
[in]	a_order	Attempted order for stencil. This is also the stencil radius.

Returns: If this returns true, a_stencil contains a valid stencil and weight.

◆ getLeastSquaresStencil()

bool EBHelmholtzDirichletEBBC::getLeastSquaresStencil	(	std::pair< Real, VoFStencil > &	a_stencil,
		const VolIndex &	a_vof,
		const VofUtils::Neighborhood	a_neighborhood,
		const DataIndex &	a_dit,
		const int	a_order
	)		const

protectedvirtual

Get stencil with specified order.

Parameters

[out]	a_stencil	Weight and stencil
[out]	a_vof	Grid cell
[out]	a_neighborhood	Neighboorhood specifier
[in]	a_dit	Data index
[in]	a_order	Attempted order for stencil. This is also the stencil radius.

Returns: If this returns true, a_stencil contains a valid stencil and weight.

Note: Stencils are computed using weighted least squares with squared inverse distance weighting.

◆ setOrder()

void EBHelmholtzDirichletEBBC::setOrder ( const int a_order )

virtual

Set BC order.

Parameters

[in] a_order Order

Note: The order is also the attempted radius of the stencil.

◆ setValue() [1/2]

void EBHelmholtzDirichletEBBC::setValue ( const int a_value )

virtual

Set the value on the EB.

Parameters

[in] a_value Constant value on EB

◆ setValue() [2/2]

void EBHelmholtzDirichletEBBC::setValue ( const std::function< Real(const RealVect &a_pos)> & a_value )

virtual

Set the value on the EB.

Parameters

[in] a_value Spatially dependent value on EB

◆ setWeight()

void EBHelmholtzDirichletEBBC::setWeight ( const int a_weight )

virtual

Set stencil weights.

Parameters

[in] a_weight Weighting factor for least squares stencils

Note: a_weight=0 can lead to ill-conditioned systems.

The documentation for this class was generated from the following files:

Source/Elliptic/CD_EBHelmholtzDirichletEBBC.H
Source/Elliptic/CD_EBHelmholtzDirichletEBBC.cpp

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Member Function Documentation

◆ applyEBFlux()

◆ getChomboLsqStencil()

◆ getJohansenStencil()

◆ getLeastSquaresStencil()

◆ setOrder()

◆ setValue() [1/2]

◆ setValue() [2/2]

◆ setWeight()