Solid (slip) wall boundary condition. More...

#include <dgcartesianbcacousticperturbsolidwall.h>

Inheritance diagram for DgBcAcousticPerturbSolidWall< nDim, SysEqn, slipWall >:

Collaboration diagram for DgBcAcousticPerturbSolidWall< nDim, SysEqn, slipWall >:

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
	DgBcAcousticPerturbSolidWall (SolverType &solver_, MInt bcId)

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

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

void	applyAtSurface (const MInt surfaceId, const MFloat NotUsed(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))

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, MBool slipWall = false>
class DgBcAcousticPerturbSolidWall< nDim, SysEqn, slipWall >

Author: Lev Liberson, Hsun-Jen Cheng

Date: 2014

Definition at line 17 of file dgcartesianbcacousticperturbsolidwall.h.

Member Typedef Documentation

◆ Base

template<MInt nDim, class SysEqn , MBool slipWall = false>

using DgBcAcousticPerturbSolidWall< nDim, SysEqn, slipWall >::Base = DgBoundaryCondition<nDim, SysEqn>

Definition at line 20 of file dgcartesianbcacousticperturbsolidwall.h.

Constructor & Destructor Documentation

◆ DgBcAcousticPerturbSolidWall()

template<MInt nDim, class SysEqn , MBool slipWall = false>

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

inline

Definition at line 29 of file dgcartesianbcacousticperturbsolidwall.h.

29: Base(solver_, bcId) {}

DgBcAcousticPerturbSolidWall::Base

DgBoundaryCondition< nDim, SysEqn > Base

Definition: dgcartesianbcacousticperturbsolidwall.h:20

Member Function Documentation

◆ apply()

template<MInt nDim, class SysEqn , MBool slipWall = false>

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

inlineoverridevirtual

Implements DgBoundaryCondition< nDim, SysEqn >.

Definition at line 34 of file dgcartesianbcacousticperturbsolidwall.h.

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

◆ applyAtSurface()

template<MInt nDim, class SysEqn , MBool slipWall>

void DgBcAcousticPerturbSolidWall< nDim, SysEqn, slipWall >::applyAtSurface	(	const MInt	surfaceId,
		const MFloat	NotUsedtime
	)

Definition at line 43 of file dgcartesianbcacousticperturbsolidwall.h.

                                                                                                      {
  // Storing the surface-information and define the inner state
  MFloat* stateL = &surfaces().variables(surfaceId, 0);
  MFloat* stateR = &surfaces().variables(surfaceId, 1);
  MFloat* innerState = (surfaces().internalSideId(surfaceId) == 0) ? stateL : stateR;
  MFloat* nodeVarsL = &surfaces().nodeVars(surfaceId, 0);
  MFloat* nodeVarsR = &surfaces().nodeVars(surfaceId, 1);
  MFloat* innerNodeVars = (surfaces().nghbrElementIds(surfaceId, 0) == -1) ? nodeVarsR : nodeVarsL;
 
  // Collect data for flux-computation
  const MInt noNodes1D = surfaces().noNodes1D(surfaceId);
  const MInt noNodes1D3 = (nDim == 3) ? noNodes1D : 1;
  const MInt dirId = surfaces().orientation(surfaceId);
 
  MFloatTensor variables(innerState, noNodes1D, noNodes1D3, SysEqn::noVars());
  MFloatTensor nodeVars(innerNodeVars, noNodes1D, noNodes1D3, SysEqn::noNodeVars());
  MFloatTensor f(flux(surfaceId), noNodes1D, noNodes1D3, SysEqn::noVars());
 
  // Computing the flux at all nodes considering zero-velocities at the wall
  //
  // Flux-components for boundary-surfaces (no-slip wall version):
  //       u'  v'  w'  p'
  // Fx = (p', 0,  0,  0)
  // Fy = (0,  p', 0,  0)
  // Fz = (0,  0,  p', 0)
  //
  // Reference: M. Bauer, Airframe noise prediction using a discontinuous
  //            Galerkin method, PhD thesis, p. 18, 2011
  //
  // Additionally, if slipWall == true, additional velocity terms are added to
  // account for non-zero mean flow in the wall-parallel direction.
  //
  // Flux-components for boundary-surfaces (slip-wall version):
  //       u'                v'                w'                p'
  // Fx = (p' + v'v0 + w'w0, 0,                0,                0)
  // Fy = (0,                p' + u'u0 + w'w0, 0,                0)
  // Fz = (0,                0,                p' + u'u0 + v'v0, 0)
 
  // Set the flux to zero
  f.set(0.0);
 
  // Set the x- or y- or z-component of the velocity to p', depending on the
  // surface orientation
  if(slipWall) {
    // If it is a slip wall, add pressure and additional velocity terms
    for(MInt i = 0; i < noNodes1D; i++) {
      for(MInt j = 0; j < noNodes1D3; j++) {
        // This following line is just so that std::inner_product can be used
        variables(i, j, dirId) = 0.0;
 
        // Set flux normal to wall (see above for equation)
        f(i, j, dirId) =
            std::inner_product(&variables(i, j, 0), &variables(i, j, nDim), &nodeVars(i, j, 0), variables(i, j, nDim));
      }
    }
  } else {
    // If it is a no-slip wall, just add the pressure term and implicity assume
    // that all velocity compontents are zero at the wall
    for(MInt i = 0; i < noNodes1D; i++) {
      for(MInt j = 0; j < noNodes1D3; j++) {
        // Set flux normal to wall (see above for equation)
        f(i, j, dirId) = variables(i, j, nDim);
      }
    }
  }
}

◆ name()

template<MInt nDim, class SysEqn , MBool slipWall = false>

MString DgBcAcousticPerturbSolidWall< nDim, SysEqn, slipWall >::name ( ) const

inlineoverridevirtual

Implements DgBoundaryCondition< nDim, SysEqn >.

Definition at line 30 of file dgcartesianbcacousticperturbsolidwall.h.

                                {
    return slipWall ? "solid wall (mean flow slip wall)" : "solid wall (mean flow no-slip wall)";
  }

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

/home/gitlab-runner/scratch/builds/oxpnswJ6/1/aia/m-AIA/m-AIA/src/DG/dgcartesianbcacousticperturbsolidwall.h

Public Types

Public Member Functions

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ Base

Constructor & Destructor Documentation

◆ DgBcAcousticPerturbSolidWall()

Member Function Documentation

◆ apply()

◆ applyAtSurface()

◆ name()