Wall boundary conditions for Rotational compressible flow problems. More...

#include <WallRotationalBC.h>

Inheritance diagram for Nektar::WallRotationalBC:

Static Public Member Functions
static CFSBndCondSharedPtr	create (const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble > > &pTraceNormals, const Array< OneD, Array< OneD, NekDouble > > &pGridVelocity, const int pSpaceDim, const int bcRegion, const int cnt)
	Creates an instance of this class. More...

Static Public Attributes
static std::string	classNameRotational
	Name of the class. More...

Protected Member Functions
void	v_Apply (Array< OneD, Array< OneD, NekDouble > > &Fwd, Array< OneD, Array< OneD, NekDouble > > &physarray, const NekDouble &time) override

Protected Member Functions inherited from Nektar::CFSBndCond
	CFSBndCond (const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble > > &pTraceNormals, const Array< OneD, Array< OneD, NekDouble > > &pGridVelocity, const int pSpaceDim, const int bcRegion, const int cnt)
	Constructor. More...

virtual	~CFSBndCond ()=default

virtual void	v_Apply (Array< OneD, Array< OneD, NekDouble > > &Fwd, Array< OneD, Array< OneD, NekDouble > > &physarray, const NekDouble &time)=0

virtual void	v_ApplyBwdWeight ()

Private Member Functions
	WallRotationalBC (const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble > > &pTraceNormals, const Array< OneD, Array< OneD, NekDouble > > &pGridVelocity, const int pSpaceDim, const int bcRegion, const int cnt)

	~WallRotationalBC (void) override=default

Friends
class	MemoryManager< WallRotationalBC >

Additional Inherited Members
Public Member Functions inherited from Nektar::CFSBndCond
void	Apply (Array< OneD, Array< OneD, NekDouble > > &Fwd, Array< OneD, Array< OneD, NekDouble > > &physarray, const NekDouble &time=0)
	Apply the boundary condition. More...

void	ApplyBwdWeight ()
	Apply the Weight of boundary condition. More...

Protected Attributes inherited from Nektar::CFSBndCond
LibUtilities::SessionReaderSharedPtr	m_session
	Session reader. More...

Array< OneD, MultiRegions::ExpListSharedPtr >	m_fields
	Array of fields. More...

Array< OneD, Array< OneD, NekDouble > >	m_traceNormals
	Trace normals. More...

Array< OneD, Array< OneD, NekDouble > >	m_gridVelocityTrace
	Grid Velocity. More...

int	m_spacedim
	Space dimension. More...

VariableConverterSharedPtr	m_varConv
	Auxiliary object to convert variables. More...

NekDouble	m_diffusionAveWeight
	Weight for average calculation of diffusion term. More...

NekDouble	m_gamma
	Parameters of the flow. More...

NekDouble	m_rhoInf

NekDouble	m_pInf

NekDouble	m_pOut

Array< OneD, NekDouble >	m_velInf

NekDouble	m_angVel

int	m_bcRegion
	Id of the boundary region. More...

int	m_offset
	Offset. More...

Detailed Description

Wall boundary conditions for Rotational compressible flow problems.

Definition at line 46 of file WallRotationalBC.h.

Constructor & Destructor Documentation

◆ WallRotationalBC()

Nektar::WallRotationalBC::WallRotationalBC	(	const LibUtilities::SessionReaderSharedPtr &	pSession,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	pFields,
		const Array< OneD, Array< OneD, NekDouble > > &	pTraceNormals,
		const Array< OneD, Array< OneD, NekDouble > > &	pGridVelocity,
		const int	pSpaceDim,
		const int	bcRegion,
		const int	cnt
	)

private

Definition at line 45 of file WallRotationalBC.cpp.

    : CFSBndCond(pSession, pFields, pTraceNormals, pGridVelocity, pSpaceDim,
                 bcRegion, cnt)
{
    m_diffusionAveWeight = 0.5;
 
    // Set up rotational boundary edge velocities
    const Array<OneD, const int> &traceBndMap = m_fields[0]->GetTraceBndMap();
 
    int eMax = m_fields[0]->GetBndCondExpansions()[m_bcRegion]->GetExpSize();
    for (int e = 0; e < eMax; ++e)
    {
        int nBCEdgePts = m_fields[0]
                             ->GetBndCondExpansions()[m_bcRegion]
                             ->GetExp(e)
                             ->GetTotPoints();
        int id2 =
            m_fields[0]->GetTrace()->GetPhys_Offset(traceBndMap[m_offset + e]);
 
        Array<OneD, NekDouble> x(nBCEdgePts, 0.0);
        Array<OneD, NekDouble> y(nBCEdgePts, 0.0);
        m_fields[0]->GetBndCondExpansions()[m_bcRegion]->GetExp(e)->GetCoords(
            x, y);
        for (int j = 0; j < nBCEdgePts; ++j)
        {
            m_gridVelocityTrace[0][id2 + j] = m_angVel * -y[j];
            m_gridVelocityTrace[1][id2 + j] = m_angVel * x[j];
        }
    }
}

References Nektar::CFSBndCond::m_angVel, Nektar::CFSBndCond::m_bcRegion, Nektar::CFSBndCond::m_diffusionAveWeight, Nektar::CFSBndCond::m_fields, Nektar::CFSBndCond::m_gridVelocityTrace, and Nektar::CFSBndCond::m_offset.

◆ ~WallRotationalBC()

Nektar::WallRotationalBC::~WallRotationalBC ( void )

overrideprivatedefault

Member Function Documentation

◆ create()

static CFSBndCondSharedPtr Nektar::WallRotationalBC::create	(	const LibUtilities::SessionReaderSharedPtr &	pSession,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	pFields,
		const Array< OneD, Array< OneD, NekDouble > > &	pTraceNormals,
		const Array< OneD, Array< OneD, NekDouble > > &	pGridVelocity,
		const int	pSpaceDim,
		const int	bcRegion,
		const int	cnt
	)

inlinestatic

Creates an instance of this class.

Definition at line 52 of file WallRotationalBC.h.

    {
        CFSBndCondSharedPtr p =
            MemoryManager<WallRotationalBC>::AllocateSharedPtr(
                pSession, pFields, pTraceNormals, pGridVelocity, pSpaceDim,
                bcRegion, cnt);
        return p;
    }

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), and CellMLToNektar.cellml_metadata::p.

◆ v_Apply()

void Nektar::WallRotationalBC::v_Apply	(	Array< OneD, Array< OneD, NekDouble > > &	Fwd,
		Array< OneD, Array< OneD, NekDouble > > &	physarray,
		const NekDouble &	time
	)

overrideprotectedvirtual

Implements Nektar::CFSBndCond.

Definition at line 81 of file WallRotationalBC.cpp.

{
    int nVariables = physarray.size();
 
    // @TODO: ALE we subtract the grid velocity ? "Set u = to ug for this one" -
    // Dave
 
    int i;
    const Array<OneD, const int> &traceBndMap = m_fields[0]->GetTraceBndMap();
 
    // Take into account that for PDE based shock capturing, eps = 0 at the
    // wall. Adjust the physical values of the trace to take user defined
    // boundaries into account
    int e, id1, id2, nBCEdgePts, eMax;
 
    eMax = m_fields[0]->GetBndCondExpansions()[m_bcRegion]->GetExpSize();
 
    for (e = 0; e < eMax; ++e)
    {
        nBCEdgePts = m_fields[0]
                         ->GetBndCondExpansions()[m_bcRegion]
                         ->GetExp(e)
                         ->GetTotPoints();
        id1 =
            m_fields[0]->GetBndCondExpansions()[m_bcRegion]->GetPhys_Offset(e);
        id2 =
            m_fields[0]->GetTrace()->GetPhys_Offset(traceBndMap[m_offset + e]);
 
        // Boundary condition for epsilon term. @TODO: Is this correct, or
        // should I do E = p/(gamma -1) + 1/2*rho(u^2 +v^2 + w^2)... or
        if (nVariables == m_spacedim + 3)
        {
            Vmath::Zero(nBCEdgePts, &Fwd[nVariables - 1][id2], 1);
        }
 
        for (i = 0; i < m_spacedim; i++)
        {
            // V = -Vin
            // Vmath::Neg(nBCEdgePts, &Fwd[i+1][id2], 1);
 
            // This now does Vg * rho + Vin
            // Vmath::Vvtvp(nBCEdgePts, &m_gridVelocityTrace[i][id2], 1,
            // &Fwd[0][id2], 1, &Fwd[i+1][id2], 1, &Fwd[i+1][id2], 1);
 
            for (int j = 0; j < nBCEdgePts; ++j)
            {
                // Fwd[i+1][id2 + j] = 2 * m_gridVelocityTrace[i][id2 + j] *
                // Fwd[0][id2 + j] - Fwd[i+1][id2 + j];
                Fwd[i + 1][id2 + j] =
                    2 * m_gridVelocityTrace[i][id2 + j] * Fwd[0][id2 + j] -
                    Fwd[i + 1][id2 + j];
            }
        }
 
        // Copy boundary adjusted values into the boundary expansion
        for (i = 0; i < nVariables; ++i)
        {
            Vmath::Vcopy(nBCEdgePts, &Fwd[i][id2], 1,
                         &(m_fields[i]
                               ->GetBndCondExpansions()[m_bcRegion]
                               ->UpdatePhys())[id1],
                         1);
        }
    }
}

References Nektar::CFSBndCond::m_bcRegion, Nektar::CFSBndCond::m_fields, Nektar::CFSBndCond::m_gridVelocityTrace, Nektar::CFSBndCond::m_offset, Nektar::CFSBndCond::m_spacedim, Vmath::Vcopy(), and Vmath::Zero().

Friends And Related Function Documentation

◆ MemoryManager< WallRotationalBC >

friend class MemoryManager< WallRotationalBC >

friend

Definition at line 1 of file WallRotationalBC.h.

Member Data Documentation

◆ classNameRotational

std::string Nektar::WallRotationalBC::classNameRotational

static

Initial value:

=
    GetCFSBndCondFactory().RegisterCreatorFunction(
        "WallRotational", WallRotationalBC::create,
        "Adiabatic rotational wall boundary condition.")

Name of the class.

Definition at line 67 of file WallRotationalBC.h.

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Private Member Functions

Friends

Additional Inherited Members