Class StokesSecondOrderWaveBcCoef is an implementation of the strategy class SAMRAI::solv::RobinBcCoefStrategy that provides Dirichlet velocity boundary condition based upon Stokes' second order theory of water waves at the inlet of the wave tank. This class is meant to be used with INSVCStaggeredHierarchyIntegrator.

#include <ibamr/StokesSecondOrderWaveBcCoef.h>

Inheritance diagram for IBAMR::StokesSecondOrderWaveBcCoef:

[legend]

Public Member Functions
	StokesSecondOrderWaveBcCoef (std::string object_name, const int comp_idx, SAMRAI::tbox::Pointer< SAMRAI::tbox::Database > input_db, SAMRAI::tbox::Pointer< SAMRAI::geom::CartesianGridGeometry< NDIM > > grid_geom)
	Constructor. More...

	~StokesSecondOrderWaveBcCoef ()
	Destructor. More...

Functions to set boundary condition coefficients a and g.
virtual void	setBcCoefs (tbox::Pointer< pdat::ArrayData< DIM, double > > &acoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &bcoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &gcoef_data, const tbox::Pointer< hier::Variable< DIM > > &variable, const hier::Patch< DIM > &patch, const hier::BoundaryBox< DIM > &bdry_box, const double fill_time=0.0) const=0
	User-supplied function to fill arrays of Robin boundary condition coefficients at a patch boundary. More...

virtual hier::IntVector< DIM >	numberOfExtensionsFillable () const=0

Functions to set boundary condition coefficients a and g.
virtual void	setBcCoefs (tbox::Pointer< pdat::ArrayData< DIM, double > > &acoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &bcoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &gcoef_data, const tbox::Pointer< hier::Variable< DIM > > &variable, const hier::Patch< DIM > &patch, const hier::BoundaryBox< DIM > &bdry_box, const double fill_time=0.0) const=0
	User-supplied function to fill arrays of Robin boundary condition coefficients at a patch boundary. More...

virtual hier::IntVector< DIM >	numberOfExtensionsFillable () const=0

Functions to set boundary condition coefficients a and g.
virtual void	setBcCoefs (tbox::Pointer< pdat::ArrayData< DIM, double > > &acoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &bcoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &gcoef_data, const tbox::Pointer< hier::Variable< DIM > > &variable, const hier::Patch< DIM > &patch, const hier::BoundaryBox< DIM > &bdry_box, const double fill_time=0.0) const=0
	User-supplied function to fill arrays of Robin boundary condition coefficients at a patch boundary. More...

virtual hier::IntVector< DIM >	numberOfExtensionsFillable () const=0

Functions to set boundary condition coefficients a and g.
virtual void	setBcCoefs (tbox::Pointer< pdat::ArrayData< DIM, double > > &acoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &bcoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &gcoef_data, const tbox::Pointer< hier::Variable< DIM > > &variable, const hier::Patch< DIM > &patch, const hier::BoundaryBox< DIM > &bdry_box, const double fill_time=0.0) const=0
	User-supplied function to fill arrays of Robin boundary condition coefficients at a patch boundary. More...

virtual hier::IntVector< DIM >	numberOfExtensionsFillable () const=0

Functions to set boundary condition coefficients a and g.
virtual void	setBcCoefs (tbox::Pointer< pdat::ArrayData< DIM, double > > &acoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &bcoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &gcoef_data, const tbox::Pointer< hier::Variable< DIM > > &variable, const hier::Patch< DIM > &patch, const hier::BoundaryBox< DIM > &bdry_box, const double fill_time=0.0) const=0
	User-supplied function to fill arrays of Robin boundary condition coefficients at a patch boundary. More...

virtual hier::IntVector< DIM >	numberOfExtensionsFillable () const=0

Functions to set boundary condition coefficients a and g.
virtual void	setBcCoefs (tbox::Pointer< pdat::ArrayData< DIM, double > > &acoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &bcoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &gcoef_data, const tbox::Pointer< hier::Variable< DIM > > &variable, const hier::Patch< DIM > &patch, const hier::BoundaryBox< DIM > &bdry_box, const double fill_time=0.0) const=0
	User-supplied function to fill arrays of Robin boundary condition coefficients at a patch boundary. More...

virtual hier::IntVector< DIM >	numberOfExtensionsFillable () const=0

Functions to set boundary condition coefficients a and g.
virtual void	setBcCoefs (tbox::Pointer< pdat::ArrayData< DIM, double > > &acoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &bcoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &gcoef_data, const tbox::Pointer< hier::Variable< DIM > > &variable, const hier::Patch< DIM > &patch, const hier::BoundaryBox< DIM > &bdry_box, const double fill_time=0.0) const=0
	User-supplied function to fill arrays of Robin boundary condition coefficients at a patch boundary. More...

virtual hier::IntVector< DIM >	numberOfExtensionsFillable () const=0

Functions to set boundary condition coefficients a and g.
virtual void	setBcCoefs (tbox::Pointer< pdat::ArrayData< DIM, double > > &acoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &bcoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &gcoef_data, const tbox::Pointer< hier::Variable< DIM > > &variable, const hier::Patch< DIM > &patch, const hier::BoundaryBox< DIM > &bdry_box, const double fill_time=0.0) const=0
	User-supplied function to fill arrays of Robin boundary condition coefficients at a patch boundary. More...

virtual hier::IntVector< DIM >	numberOfExtensionsFillable () const=0

Functions to set boundary condition coefficients a and g.
virtual void	setBcCoefs (tbox::Pointer< pdat::ArrayData< DIM, double > > &acoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &bcoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &gcoef_data, const tbox::Pointer< hier::Variable< DIM > > &variable, const hier::Patch< DIM > &patch, const hier::BoundaryBox< DIM > &bdry_box, const double fill_time=0.0) const=0
	User-supplied function to fill arrays of Robin boundary condition coefficients at a patch boundary. More...

virtual hier::IntVector< DIM >	numberOfExtensionsFillable () const=0

Functions to set boundary condition coefficients a and g.
virtual void	setBcCoefs (tbox::Pointer< pdat::ArrayData< DIM, double > > &acoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &bcoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &gcoef_data, const tbox::Pointer< hier::Variable< DIM > > &variable, const hier::Patch< DIM > &patch, const hier::BoundaryBox< DIM > &bdry_box, const double fill_time=0.0) const=0
	User-supplied function to fill arrays of Robin boundary condition coefficients at a patch boundary. More...

virtual hier::IntVector< DIM >	numberOfExtensionsFillable () const=0

Functions to set boundary condition coefficients a and g.
virtual void	setBcCoefs (tbox::Pointer< pdat::ArrayData< DIM, double > > &acoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &bcoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &gcoef_data, const tbox::Pointer< hier::Variable< DIM > > &variable, const hier::Patch< DIM > &patch, const hier::BoundaryBox< DIM > &bdry_box, const double fill_time=0.0) const=0
	User-supplied function to fill arrays of Robin boundary condition coefficients at a patch boundary. More...

virtual hier::IntVector< DIM >	numberOfExtensionsFillable () const=0

Functions to set boundary condition coefficients a and g.
virtual void	setBcCoefs (tbox::Pointer< pdat::ArrayData< DIM, double > > &acoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &bcoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &gcoef_data, const tbox::Pointer< hier::Variable< DIM > > &variable, const hier::Patch< DIM > &patch, const hier::BoundaryBox< DIM > &bdry_box, const double fill_time=0.0) const=0
	User-supplied function to fill arrays of Robin boundary condition coefficients at a patch boundary. More...

virtual hier::IntVector< DIM >	numberOfExtensionsFillable () const=0

Functions to set boundary condition coefficients a and g.
virtual void	setBcCoefs (tbox::Pointer< pdat::ArrayData< DIM, double > > &acoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &bcoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &gcoef_data, const tbox::Pointer< hier::Variable< DIM > > &variable, const hier::Patch< DIM > &patch, const hier::BoundaryBox< DIM > &bdry_box, const double fill_time=0.0) const=0
	User-supplied function to fill arrays of Robin boundary condition coefficients at a patch boundary. More...

virtual hier::IntVector< DIM >	numberOfExtensionsFillable () const=0

Functions to set boundary condition coefficients a and g.
virtual void	setBcCoefs (tbox::Pointer< pdat::ArrayData< DIM, double > > &acoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &bcoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &gcoef_data, const tbox::Pointer< hier::Variable< DIM > > &variable, const hier::Patch< DIM > &patch, const hier::BoundaryBox< DIM > &bdry_box, const double fill_time=0.0) const=0
	User-supplied function to fill arrays of Robin boundary condition coefficients at a patch boundary. More...

virtual hier::IntVector< DIM >	numberOfExtensionsFillable () const=0

Functions to set boundary condition coefficients a and g.
virtual void	setBcCoefs (tbox::Pointer< pdat::ArrayData< DIM, double > > &acoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &bcoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &gcoef_data, const tbox::Pointer< hier::Variable< DIM > > &variable, const hier::Patch< DIM > &patch, const hier::BoundaryBox< DIM > &bdry_box, const double fill_time=0.0) const=0
	User-supplied function to fill arrays of Robin boundary condition coefficients at a patch boundary. More...

virtual hier::IntVector< DIM >	numberOfExtensionsFillable () const=0

Functions to set boundary condition coefficients a and g.
virtual void	setBcCoefs (tbox::Pointer< pdat::ArrayData< DIM, double > > &acoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &bcoef_data, tbox::Pointer< pdat::ArrayData< DIM, double > > &gcoef_data, const tbox::Pointer< hier::Variable< DIM > > &variable, const hier::Patch< DIM > &patch, const hier::BoundaryBox< DIM > &bdry_box, const double fill_time=0.0) const=0
	User-supplied function to fill arrays of Robin boundary condition coefficients at a patch boundary. More...

virtual hier::IntVector< DIM >	numberOfExtensionsFillable () const=0

Implementation of SAMRAI::solv::RobinBcCoefStrategy interface.
std::string	d_object_name

int	d_comp_idx = IBTK::invalid_index

IBTK::muParserRobinBcCoefs	d_muparser_bcs

SAMRAI::tbox::Pointer< SAMRAI::geom::CartesianGridGeometry< NDIM > >	d_grid_geom

double	d_depth
	Wave parameters. More...

double	d_omega

double	d_wave_number

double	d_amplitude

double	d_gravity

double	d_num_interface_cells

void	setBcCoefs (SAMRAI::tbox::Pointer< SAMRAI::pdat::ArrayData< NDIM, double > > &acoef_data, SAMRAI::tbox::Pointer< SAMRAI::pdat::ArrayData< NDIM, double > > &bcoef_data, SAMRAI::tbox::Pointer< SAMRAI::pdat::ArrayData< NDIM, double > > &gcoef_data, const SAMRAI::tbox::Pointer< SAMRAI::hier::Variable< NDIM > > &variable, const SAMRAI::hier::Patch< NDIM > &patch, const SAMRAI::hier::BoundaryBox< NDIM > &bdry_box, double fill_time=0.0) const override
	Function to fill arrays of Robin boundary condition coefficients at a patch boundary. More...

SAMRAI::hier::IntVector< NDIM >	numberOfExtensionsFillable () const override

	StokesSecondOrderWaveBcCoef ()=delete
	Default constructor. More...

	StokesSecondOrderWaveBcCoef (const StokesSecondOrderWaveBcCoef &from)=delete
	Copy constructor. More...

StokesSecondOrderWaveBcCoef &	operator= (const StokesSecondOrderWaveBcCoef &that)=delete
	Assignment operator. More...

void	getFromInput (SAMRAI::tbox::Pointer< SAMRAI::tbox::Database > db)

Constructor & Destructor Documentation

◆ StokesSecondOrderWaveBcCoef() [1/3]

IBAMR::StokesSecondOrderWaveBcCoef::StokesSecondOrderWaveBcCoef	(	std::string	object_name,
		const int	comp_idx,
		SAMRAI::tbox::Pointer< SAMRAI::tbox::Database >	input_db,
		SAMRAI::tbox::Pointer< SAMRAI::geom::CartesianGridGeometry< NDIM > >	grid_geom
	)

◆ ~StokesSecondOrderWaveBcCoef()

IBAMR::StokesSecondOrderWaveBcCoef::~StokesSecondOrderWaveBcCoef ( )

◆ StokesSecondOrderWaveBcCoef() [2/3]

IBAMR::StokesSecondOrderWaveBcCoef::StokesSecondOrderWaveBcCoef ( )

privatedelete

Note: This constructor is not implemented and should not be used.

◆ StokesSecondOrderWaveBcCoef() [3/3]

IBAMR::StokesSecondOrderWaveBcCoef::StokesSecondOrderWaveBcCoef ( const StokesSecondOrderWaveBcCoef & from )

privatedelete

Note: This constructor is not implemented and should not be used.

Parameters

from	The value to copy to this object.

Member Function Documentation

◆ setBcCoefs() [1/2]

void IBAMR::StokesSecondOrderWaveBcCoef::setBcCoefs	(	SAMRAI::tbox::Pointer< SAMRAI::pdat::ArrayData< NDIM, double > > &	acoef_data,
		SAMRAI::tbox::Pointer< SAMRAI::pdat::ArrayData< NDIM, double > > &	bcoef_data,
		SAMRAI::tbox::Pointer< SAMRAI::pdat::ArrayData< NDIM, double > > &	gcoef_data,
		const SAMRAI::tbox::Pointer< SAMRAI::hier::Variable< NDIM > > &	variable,
		const SAMRAI::hier::Patch< NDIM > &	patch,
		const SAMRAI::hier::BoundaryBox< NDIM > &	bdry_box,
		double	fill_time = `0.0`
	)		const

override

See also: SAMRAI::solv::RobinBcCoefStrategy::setBcCoefs()

Parameters

acoef_data	Boundary coefficient data. The array will have been defined to include index range for corresponding to the boundary box bdry_box and appropriate for the alignment of the given variable. If this is a null pointer, then the calling function is not interested in a, and you can disregard it.
bcoef_data	Boundary coefficient data. This array is exactly like acoef_data, except that it is to be filled with the b coefficient.
gcoef_data	Boundary coefficient data. This array is exactly like acoef_data, except that it is to be filled with the g coefficient.
variable	Variable to set the coefficients for. If implemented for multiple variables, this parameter can be used to determine which variable's coefficients are being sought.
patch	Patch requiring bc coefficients.
bdry_box	Boundary box showing where on the boundary the coefficient data is needed.
fill_time	Solution time corresponding to filling, for use when coefficients are time-dependent.

◆ numberOfExtensionsFillable() [1/2]

SAMRAI::hier::IntVector<NDIM> IBAMR::StokesSecondOrderWaveBcCoef::numberOfExtensionsFillable ( ) const

override

◆ operator=()

StokesSecondOrderWaveBcCoef& IBAMR::StokesSecondOrderWaveBcCoef::operator= ( const StokesSecondOrderWaveBcCoef & that )

privatedelete

Note: This operator is not implemented and should not be used.

Parameters

that	The value to assign to this object.

Returns: A reference to this object.

◆ getFromInput()

void IBAMR::StokesSecondOrderWaveBcCoef::getFromInput ( SAMRAI::tbox::Pointer< SAMRAI::tbox::Database > db )

private

Get wave parameters from input db.

◆ setBcCoefs() [2/2]

virtual void SAMRAI::solv::RobinBcCoefStrategy< DIM >::setBcCoefs	(	tbox::Pointer< pdat::ArrayData< DIM, double > > &	acoef_data,
		tbox::Pointer< pdat::ArrayData< DIM, double > > &	bcoef_data,
		tbox::Pointer< pdat::ArrayData< DIM, double > > &	gcoef_data,
		const tbox::Pointer< hier::Variable< DIM > > &	variable,
		const hier::Patch< DIM > &	patch,
		const hier::BoundaryBox< DIM > &	bdry_box,
		const double	fill_time = `0.0`
	)		const

pure virtualinherited

This class specifies the Robin boundary condition coefficients at discrete locations on the patch boundary. Though these locations are defined by boundary box object, they do not necessarily coincide with the centers of the cells referred to by those boxes. These locations typically coincide with the nodes or face centers which do lie on the patch boundary. Accordingly, you use this function to provide the boundary coefficients at those locations by filling an array at indices corresponding to those locations.

When setting the values of the boundary condition coefficients it is useful to note that for any cell (i,j,k), the indices of the sides, edges and nodes are easily determined. The index on the lower side of the cell is the same as the index of the cell, whereas the index on the upper side of the cell has the next higher value. In 2D, the cell and its surrounding nodes and faces has the following indices:

*
*       (i,j+1)----(i,j+1)---(i+1,j+1)
*          |                     |
*          |                     |
*          |                     |
*          |                     |
*        (i,j)      (i,j)     (i+1,j)
*          |                     |
*          |                     |
*          |                     |
*          |                     |
*        (i,j)------(i,j)-----(i+1,j)
*
*

Once this is understood, translation between the index in the boundary box index space to the index of things on the boundary is simple.

The boundary condition coefficients should be placed in the pdat::ArrayData<DIM> objects, acoef_data and gcoef_data (see argument list), which are dimensioned to contain the indices of the points alligned with variable and lying on the the boundary defined by bdry_box.

This function is only used with type-1 boundary boxes, such as faces in 3D. Other types of boundaries do not have a well-defined surface normal.

The parameter variable is passed through to tell the implementation of this function what variable to set the coefficients for. You may wish to ignore it if your implementation is intended for a specific variable.

Parameters

acoef_data	boundary coefficient data. The array will have been defined to include index range for corresponding to the boundary box `bdry_box` and appropriate for the alignment of the given variable. If this is a null pointer, then the calling function is not interested in a, and you can disregard it.
bcoef_data	boundary coefficient data. This array is exactly like `acoef_data`, except that it is to be filled with the b coefficient.
gcoef_data	boundary coefficient data. This array is exactly like `acoef_data`, except that it is to be filled with the g coefficient.
variable	variable to set the coefficients for. If implemented for multiple variables, this parameter can be used to determine which variable's coefficients are being sought.
patch	patch requiring bc coefficients
bdry_box	boundary box showing where on the boundary the coefficient data is needed.
fill_time	solution time corresponding to filling, for use when coefficients are time-dependent.

◆ numberOfExtensionsFillable() [2/2]

virtual hier::IntVector<DIM> SAMRAI::solv::RobinBcCoefStrategy< DIM >::numberOfExtensionsFillable ( ) const

pure virtualinherited

Member Data Documentation

◆ d_object_name

std::string IBAMR::StokesSecondOrderWaveBcCoef::d_object_name

private

Book-keeping.

◆ d_comp_idx

int IBAMR::StokesSecondOrderWaveBcCoef::d_comp_idx = IBTK::invalid_index

private

Velocity component index.

◆ d_muparser_bcs

IBTK::muParserRobinBcCoefs IBAMR::StokesSecondOrderWaveBcCoef::d_muparser_bcs

private

muparser object for filling boundary conditions other than wave inlet.

◆ d_grid_geom

SAMRAI::tbox::Pointer<SAMRAI::geom::CartesianGridGeometry<NDIM> > IBAMR::StokesSecondOrderWaveBcCoef::d_grid_geom

private

The Cartesian grid geometry object provides the extents of the computational domain.

◆ d_depth

double IBAMR::StokesSecondOrderWaveBcCoef::d_depth

private

Parameters

d_wave_number	: Wave number of dominant wave component [$2\pi/m$]
d_amplitude	: Amplitude of the dominant wave component [m]
d_depth	: Depth of water, from sea bed to still water level [m]
d_gravity	: Acceleration due to gravity [$m/s^2$]
d_omega	: Angular frequency [$2 \pi/s$]

\NOTE d_omega, d_wave_number, and d_depth obtained from the input file should satisfy the dispersion relation to get a consistent Stokes wave.

◆ d_omega

double IBAMR::StokesSecondOrderWaveBcCoef::d_omega

private

◆ d_wave_number

double IBAMR::StokesSecondOrderWaveBcCoef::d_wave_number

private

◆ d_amplitude

double IBAMR::StokesSecondOrderWaveBcCoef::d_amplitude

private

◆ d_gravity

double IBAMR::StokesSecondOrderWaveBcCoef::d_gravity

private

◆ d_num_interface_cells

double IBAMR::StokesSecondOrderWaveBcCoef::d_num_interface_cells

private

Number of interface cells.

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

include/ibamr/StokesSecondOrderWaveBcCoef.h

Public Member Functions

Implementation of SAMRAI::solv::RobinBcCoefStrategy interface.

Constructor & Destructor Documentation

◆ StokesSecondOrderWaveBcCoef() [1/3]

◆ ~StokesSecondOrderWaveBcCoef()

◆ StokesSecondOrderWaveBcCoef() [2/3]

◆ StokesSecondOrderWaveBcCoef() [3/3]

Member Function Documentation

◆ setBcCoefs() [1/2]

◆ numberOfExtensionsFillable() [1/2]

◆ operator=()

◆ getFromInput()

◆ setBcCoefs() [2/2]

◆ numberOfExtensionsFillable() [2/2]

Member Data Documentation

◆ d_object_name

◆ d_comp_idx

◆ d_muparser_bcs

◆ d_grid_geom

◆ d_depth

◆ d_omega

◆ d_wave_number

◆ d_amplitude

◆ d_gravity

◆ d_num_interface_cells