#include <dgcartesianbcacousticperturbstraightductexit.h>

Inheritance diagram for DgBcAcousticPerturbStraightDuctExit< nDim, SysEqn >:

Collaboration diagram for DgBcAcousticPerturbStraightDuctExit< nDim, SysEqn >:

Public Types
using	Base = DgBoundaryCondition< nDim, SysEqn >

Public Types inherited from DgBoundaryCondition< nDim, SysEqn >
using	SysEqn = SysEqn

using	SolverType = DgCartesianSolver< nDim, SysEqn >

using	ElementCollector = maia::dg::collector::ElementCollector< nDim, SysEqn >

using	HElementCollector = maia::dg::collector::HElementCollector< nDim, SysEqn >

using	SurfaceCollector = maia::dg::collector::SurfaceCollector< nDim, SysEqn >

Public Member Functions
	DgBcAcousticPerturbStraightDuctExit (SolverType &solver_, MInt bcId)

MString	name () const override
	Returns name of boundary condition. More...

void	init () override

void	apply (const MFloat time) override
	Apply method to apply boundary condition. More...

void	applyAtSurface (const MInt surfaceId, const MFloat time)

Public Member Functions inherited from DgBoundaryCondition< nDim, SysEqn >
virtual	~DgBoundaryCondition ()
	Destructor must be virtual. More...

void	init (const MInt begin_, const MInt end_)
	Init method to initialize boundary condition for range of surfaces. More...

virtual void	apply (const MFloat time)=0
	Apply method to apply boundary condition. More...

virtual MString	name () const=0
	Returns name of boundary condition. More...

	DgBoundaryCondition (SolverType &solver_, MInt bcId)
	Constructor saves arguments to member variables. More...

MInt	id () const
	Return boundary condition if of this boundary condition. More...

MInt	begin () const
	Return index of first surface. More...

MInt	end () const
	Return index of one-past-last surface. More...

MInt	count () const
	Return number of boundary surfaces. More...

virtual MInt	noRestartVars () const

virtual MInt	getLocalNoNodes () const
	Return local number of nodes. More...

virtual MString	restartVarName (const MInt NotUsed(id)) const
	Return name of restart variable. More...

virtual void	setRestartVariable (const MInt NotUsed(id), const MFloat *const NotUsed(data))
	Copy restart variable data from pointer to boundary condition class. More...

virtual void	getRestartVariable (const MInt NotUsed(id), MFloat *const NotUsed(data)) const
	Copy restart variable data from boundary condition class to pointer. More...

virtual MInt	noBcElements () const

virtual MBool	hasBcElement (const MInt NotUsed(elementId)) const

virtual MInt	noCellDataDlb () const

virtual MInt	cellDataTypeDlb (const MInt NotUsed(dataId)) const

virtual MInt	cellDataSizeDlb (const MInt NotUsed(dataId), const MInt NotUsed(cellId)) const

virtual void	getCellDataDlb (const MInt NotUsed(dataId), MFloat *const NotUsed(data)) const

virtual void	getCellDataDlb (const MInt NotUsed(dataId), MInt *const NotUsed(data)) const

virtual void	setCellDataDlb (const MInt NotUsed(dataId), const MFloat *const NotUsed(data))

virtual void	setCellDataDlb (const MInt NotUsed(dataId), const MInt *const NotUsed(data))

Private Attributes
DgBcAcousticPerturbSolidWall< nDim, SysEqn >	m_wallBc

MFloat	m_ductPosition [2] {}

Additional Inherited Members
Protected Member Functions inherited from DgBoundaryCondition< nDim, SysEqn >
SolverType &	solver ()
	Return reference to solver. More...

SysEqn &	sysEqn ()
	Return reference to SysEqn object. More...

MFloat *	flux (const MInt i)
	Return pointer to surface flux. More...

ElementCollector &	elements ()
	Return reference to elements. More...

MInt	getElementByCellId (const MInt cellId) const
	Return element id corresponding to given cell id. More...

HElementCollector &	helements ()
	Return reference to h-elements. More...

SurfaceCollector &	surfaces ()
	Return reference to surfaces. More...

MBool	isMpiSurface (const MInt id_) const
	Return true if surface is a MPI surface. More...

MBool	needHElementForCell (const MInt cellId)
	Return if h-element is needed for given cell. More...

MInt	getHElementId (const MInt elementId)
	Return h-element id for an element. More...

const DgInterpolation &	interpolation (const MInt polyDeg, const MInt noNodes1D) const
	Return interpolation. More...

MInt	integrationMethod () const
	Return integration method. More...

MInt	timeIntegrationScheme () const
	Return time integration scheme. More...

MInt	maxPolyDeg () const
	Return maximum polynomial degree. More...

MInt	maxNoNodes1D () const
	Return maximum number of nodes. More...

MFloat	dt () const
	Return current time step size. More...

MBool	isRestart () const
	Return if a restart is performed. More...

void	subTimeStepRk (const MFloat dt_, const MInt stage, const MInt totalSize, const MFloat const rhs, MFloat const variables, MFloat *const timeIntStorage)
	Access to time integration method. More...

void	calcVolumeIntegral (const MInt noElements, ElementCollector &elem, F &fluxFct)

void	resetBuffer (const MInt totalSize, MFloat *const buffer)

void	applyJacobian (const MInt noElements, ElementCollector &elem)

void	calcSourceTerms (const MFloat t, const MInt noElements, ElementCollector &elem, F &sourceFct)

void	calcSurfaceIntegral (const MInt begin_, const MInt end_, ElementCollector &elem, SurfaceCollector &surf, HElementCollector &helem, const MInt noHElements)

void	calcRegularSurfaceFlux (const MInt begin_, const MInt end_, SurfaceCollector &surf, F &riemannFct)

Detailed Description

template<MInt nDim, class SysEqn>
class DgBcAcousticPerturbStraightDuctExit< nDim, SysEqn >

Definition at line 20 of file dgcartesianbcacousticperturbstraightductexit.h.

Member Typedef Documentation

◆ Base

template<MInt nDim, class SysEqn >

using DgBcAcousticPerturbStraightDuctExit< nDim, SysEqn >::Base = DgBoundaryCondition<nDim, SysEqn>

Definition at line 23 of file dgcartesianbcacousticperturbstraightductexit.h.

Constructor & Destructor Documentation

◆ DgBcAcousticPerturbStraightDuctExit()

template<MInt nDim, class SysEqn >

DgBcAcousticPerturbStraightDuctExit< nDim, SysEqn >::DgBcAcousticPerturbStraightDuctExit	(	SolverType &	solver_,
		MInt	bcId
	)

inline

Definition at line 33 of file dgcartesianbcacousticperturbstraightductexit.h.

33: Base(solver_, bcId), m_wallBc(solver_, bcId) {}

DgBcAcousticPerturbStraightDuctExit::Base

DgBoundaryCondition< nDim, SysEqn > Base

Definition: dgcartesianbcacousticperturbstraightductexit.h:23

DgBcAcousticPerturbStraightDuctExit::m_wallBc

DgBcAcousticPerturbSolidWall< nDim, SysEqn > m_wallBc

Definition: dgcartesianbcacousticperturbstraightductexit.h:46

Member Function Documentation

◆ apply()

template<MInt nDim, class SysEqn >

void DgBcAcousticPerturbStraightDuctExit< nDim, SysEqn >::apply ( const MFloat time )

inlineoverridevirtual

Implements DgBoundaryCondition< nDim, SysEqn >.

Definition at line 38 of file dgcartesianbcacousticperturbstraightductexit.h.

                                         {
    maia::dg::bc::loop(&DgBcAcousticPerturbStraightDuctExit::applyAtSurface, this, begin(), end(), time);
  }

◆ applyAtSurface()

template<MInt nDim, class SysEqn >

void DgBcAcousticPerturbStraightDuctExit< nDim, SysEqn >::applyAtSurface	(	const MInt	surfaceId,
		const MFloat	time
	)

Definition at line 81 of file dgcartesianbcacousticperturbstraightductexit.h.

                                                                                                              {
  const MFloat* const coordinates = &surfaces().coords(surfaceId, 0);
  const MInt sOrientation = surfaces().orientation(surfaceId);
 
  const MFloat pos = coordinates[(sOrientation + 1) % 2];
  const MFloat lb = m_ductPosition[0];
  const MFloat ub = m_ductPosition[1];
 
  // Check position of surface center
  if(lb > pos || pos > ub) {
    // Apply wall boundary condition
    m_wallBc.applyAtSurface(surfaceId, time);
  } else {
    // Calculate flux with the inner state of a surface and boundary state
    const MInt noNodes1D = surfaces().noNodes1D(surfaceId);
 
    const MFloat* surfaceStateL = &surfaces().variables(surfaceId, 0);
    const MFloat* surfaceStateR = &surfaces().variables(surfaceId, 1);
    const MFloat* outerState = (surfaces().nghbrElementIds(surfaceId, 0) == -1) ? surfaceStateL : surfaceStateR;
    MFloatTensor oS(const_cast<MFloat*>(outerState), noNodes1D, SysEqn::noVars());
 
    // Check the orientation of the Duct-exit and set the outer state
    if(sOrientation == 1) {
      for(MInt n = 0; n < noNodes1D; n++) {
        // for y-orientation x-velocity is zero
        oS(n, 0) = 0.0;
        oS(n, 1) = sin(2 * PI * time);
        oS(n, 2) = sin(2 * PI * time);
      }
    } else {
      for(MInt n = 0.0; n < noNodes1D; n++) {
        // for x orientation y-velocity is zero
        oS(n, 0) = sin(2 * PI * time);
        oS(n, 1) = 0;
        oS(n, 2) = sin(2 * PI * time);
      }
    }
 
    // Call regular Riemann solver
    const MFloat* nodeVarsL = &surfaces().nodeVars(surfaceId, 0);
    const MFloat* nodeVarsR = &surfaces().nodeVars(surfaceId, 1);
    const MFloat* innerNodeVars = (surfaces().nghbrElementIds(surfaceId, 0) == -1) ? nodeVarsR : nodeVarsL;
    sysEqn().calcRiemann(innerNodeVars, innerNodeVars, surfaceStateL, surfaceStateR, noNodes1D, sOrientation,
                         flux(surfaceId));
  }
}

◆ init()

template<MInt nDim, class SysEqn >

void DgBcAcousticPerturbStraightDuctExit< nDim, SysEqn >::init

overridevirtual

Init method that may be implemented by the derived classes. This is called by the init(MInt, MInt) method after the begin/end indices are saved and accessible by begin()/end(). If a derived class does not need to initialize anything, it just does not provide an implementation.

Reimplemented from DgBoundaryCondition< nDim, SysEqn >.

Definition at line 62 of file dgcartesianbcacousticperturbstraightductexit.h.