Nektar::MovingFrameWall Class Reference

#include <MovingFrameWall.h>

Inheritance diagram for Nektar::MovingFrameWall:

Public Member Functions
	~MovingFrameWall () override
Public Member Functions inherited from Nektar::StaticWall
	~StaticWall () override
Public Member Functions inherited from Nektar::IncBaseCondition
virtual	~IncBaseCondition ()
void	Update (const Array< OneD, const Array< OneD, NekDouble > > &fields, const Array< OneD, const Array< OneD, NekDouble > > &Adv, std::map< std::string, NekDouble > &params)
void	Initialise (const LibUtilities::SessionReaderSharedPtr &pSession)

Static Public Member Functions
static IncBaseConditionSharedPtr	create (const LibUtilities::SessionReaderSharedPtr pSession, Array< OneD, MultiRegions::ExpListSharedPtr > pFields, Array< OneD, SpatialDomains::BoundaryConditionShPtr > cond, Array< OneD, MultiRegions::ExpListSharedPtr > exp, int nbnd, int spacedim, int bnddim)
Static Public Member Functions inherited from Nektar::StaticWall
static IncBaseConditionSharedPtr	create (const LibUtilities::SessionReaderSharedPtr pSession, Array< OneD, MultiRegions::ExpListSharedPtr > pFields, Array< OneD, SpatialDomains::BoundaryConditionShPtr > cond, Array< OneD, MultiRegions::ExpListSharedPtr > exp, int nbnd, int spacedim, int bnddim)

Static Public Attributes
static std::string	className
Static Public Attributes inherited from Nektar::StaticWall
static std::string	className

Protected Member Functions
void	v_Initialise (const LibUtilities::SessionReaderSharedPtr &pSession) override
void	v_Update (const Array< OneD, const Array< OneD, NekDouble > > &fields, const Array< OneD, const Array< OneD, NekDouble > > &Adv, std::map< std::string, NekDouble > &params) override
	v_Update set correct BCs (in wavespace) More...
	MovingFrameWall (const LibUtilities::SessionReaderSharedPtr pSession, Array< OneD, MultiRegions::ExpListSharedPtr > pFields, Array< OneD, SpatialDomains::BoundaryConditionShPtr > cond, Array< OneD, MultiRegions::ExpListSharedPtr > exp, int nbnd, int spacedim, int bnddim)
Protected Member Functions inherited from Nektar::StaticWall
void	v_Initialise (const LibUtilities::SessionReaderSharedPtr &pSession) override
void	v_Update (const Array< OneD, const Array< OneD, NekDouble > > &fields, const Array< OneD, const Array< OneD, NekDouble > > &Adv, std::map< std::string, NekDouble > &params) override
	StaticWall (const LibUtilities::SessionReaderSharedPtr pSession, Array< OneD, MultiRegions::ExpListSharedPtr > pFields, Array< OneD, SpatialDomains::BoundaryConditionShPtr > cond, Array< OneD, MultiRegions::ExpListSharedPtr > exp, int nbnd, int spacedim, int bnddim)
Protected Member Functions inherited from Nektar::IncBaseCondition
	IncBaseCondition (const LibUtilities::SessionReaderSharedPtr pSession, Array< OneD, MultiRegions::ExpListSharedPtr > pFields, Array< OneD, SpatialDomains::BoundaryConditionShPtr > cond, Array< OneD, MultiRegions::ExpListSharedPtr > exp, int nbnd, int spacedim, int bnddim)
virtual void	v_Initialise (const LibUtilities::SessionReaderSharedPtr &pSession)
virtual void	v_Update (const Array< OneD, const Array< OneD, NekDouble > > &fields, const Array< OneD, const Array< OneD, NekDouble > > &Adv, std::map< std::string, NekDouble > &params)
void	ExtrapolateArray (const int numCalls, Array< OneD, Array< OneD, Array< OneD, NekDouble > > > &array)
void	RollOver (Array< OneD, Array< OneD, Array< OneD, NekDouble > > > &input)
void	AddVisPressureBCs (const Array< OneD, const Array< OneD, NekDouble > > &fields, Array< OneD, Array< OneD, NekDouble > > &N, std::map< std::string, NekDouble > &params)
void	AddRigidBodyAcc (Array< OneD, Array< OneD, NekDouble > > &N, std::map< std::string, NekDouble > &params, int npts0)
void	RigidBodyVelocity (Array< OneD, Array< OneD, NekDouble > > &velocities, std::map< std::string, NekDouble > &params, int npts0)
void	InitialiseCoords (std::map< std::string, NekDouble > &params)
void	SetNumPointsOnPlane0 (int &npointsPlane0)

Protected Attributes
int	m_pressure
bool	m_hasVels
bool	m_hasPressure
Protected Attributes inherited from Nektar::IncBaseCondition
int	m_spacedim
int	m_bnddim
	bounday dimensionality More...
int	m_nbnd
int	m_numCalls
int	m_intSteps
std::map< int, SpatialDomains::BoundaryConditionShPtr >	m_BndConds
std::map< int, MultiRegions::ExpListSharedPtr >	m_BndExp
int	m_npoints
Array< OneD, Array< OneD, NekDouble > >	m_coords
MultiRegions::ExpListSharedPtr	m_bndElmtExps
MultiRegions::ExpListSharedPtr	m_field
Array< OneD, Array< OneD, Array< OneD, NekDouble > > >	m_viscous
int	m_pressure
std::string	classname

Friends
class	MemoryManager< MovingFrameWall >

Additional Inherited Members
Static Protected Attributes inherited from Nektar::IncBaseCondition
static NekDouble	StifflyStable_Betaq_Coeffs [3][3]
static NekDouble	StifflyStable_Alpha_Coeffs [3][3]
static NekDouble	StifflyStable_Gamma0_Coeffs [3]

Detailed Description

Definition at line 48 of file MovingFrameWall.h.

Constructor & Destructor Documentation

~MovingFrameWall()

Nektar::MovingFrameWall::~MovingFrameWall ( )

override

Definition at line 72 of file MovingFrameWall.cpp.

{
}

MovingFrameWall()

Nektar::MovingFrameWall::MovingFrameWall ( const LibUtilities::SessionReaderSharedPtr  pSession, Array< OneD, MultiRegions::ExpListSharedPtr >  pFields, Array< OneD, SpatialDomains::BoundaryConditionShPtr >  cond, Array< OneD, MultiRegions::ExpListSharedPtr >  exp, int  nbnd, int  spacedim, int  bnddim  )

protected

Definition at line 45 of file MovingFrameWall.cpp.

    : StaticWall(pSession, pFields, cond, exp, nbnd, spacedim, bnddim)
{
    classname = "MovingFrameWall";
    m_hasVels = false;
    for (size_t i = 0; i < m_bnddim; ++i)
    {
        if (cond[i]->GetUserDefined() == classname)
        {
            m_hasVels   = true;
            m_BndExp[i] = exp[i];
        }
    }
    m_pressure    = cond.size() - 1;
    m_hasPressure = cond[m_pressure]->GetUserDefined() == classname;
    if (m_hasPressure)
    {
        m_BndExp[m_pressure] = exp[m_pressure];
    }
}

References Nektar::IncBaseCondition::classname, Nektar::IncBaseCondition::m_bnddim, Nektar::IncBaseCondition::m_BndExp, m_hasPressure, m_hasVels, and m_pressure.

Member Function Documentation

create()

static IncBaseConditionSharedPtr Nektar::MovingFrameWall::create ( const LibUtilities::SessionReaderSharedPtr  pSession, Array< OneD, MultiRegions::ExpListSharedPtr >  pFields, Array< OneD, SpatialDomains::BoundaryConditionShPtr >  cond, Array< OneD, MultiRegions::ExpListSharedPtr >  exp, int  nbnd, int  spacedim, int  bnddim  )

inlinestatic

Definition at line 53 of file MovingFrameWall.h.

    {
        IncBaseConditionSharedPtr p =
            MemoryManager<MovingFrameWall>::AllocateSharedPtr(
                pSession, pFields, cond, exp, nbnd, spacedim, bnddim);
        p->Initialise(pSession);
        return p;
    }

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

v_Initialise()

void Nektar::MovingFrameWall::v_Initialise ( const LibUtilities::SessionReaderSharedPtr &  pSession )

overrideprotectedvirtual

Reimplemented from Nektar::IncBaseCondition.

Definition at line 76 of file MovingFrameWall.cpp.

{
    IncBaseCondition::v_Initialise(pSession);
    m_field->GetBndElmtExpansion(m_nbnd, m_bndElmtExps, false);
    if (m_hasPressure)
    {
        m_viscous =
            Array<OneD, Array<OneD, Array<OneD, NekDouble>>>(m_intSteps);
        for (int n = 0; n < m_intSteps; ++n)
        {
            m_viscous[n] = Array<OneD, Array<OneD, NekDouble>>(m_bnddim);
            for (int i = 0; i < m_bnddim; ++i)
            {
                m_viscous[n][i] = Array<OneD, NekDouble>(m_npoints, 0.0);
            }
        }
    }
}

References Nektar::IncBaseCondition::m_bnddim, Nektar::IncBaseCondition::m_bndElmtExps, Nektar::IncBaseCondition::m_field, m_hasPressure, Nektar::IncBaseCondition::m_intSteps, Nektar::IncBaseCondition::m_nbnd, Nektar::IncBaseCondition::m_npoints, Nektar::IncBaseCondition::m_viscous, and Nektar::IncBaseCondition::v_Initialise().

v_Update()

void Nektar::MovingFrameWall::v_Update ( const Array< OneD, const Array< OneD, NekDouble > > &  fields, const Array< OneD, const Array< OneD, NekDouble > > &  Adv, std::map< std::string, NekDouble > &  params  )

overrideprotectedvirtual

v_Update set correct BCs (in wavespace)

Parameters

fields
Adv	is in wavespace for 3DH1D
params

Reimplemented from Nektar::IncBaseCondition.

Definition at line 100 of file MovingFrameWall.cpp.

{
    int nptsPlane0 = 0;
    SetNumPointsOnPlane0(nptsPlane0);
    if (params.find("Omega_x") != params.end() ||
        params.find("Omega_y") != params.end() ||
        params.find("Omega_z") != params.end())
    {
        InitialiseCoords(params);
    }
    // pressure
    if (m_hasPressure && fields.size() > 0)
    {
        ++m_numCalls;
 
        Array<OneD, Array<OneD, NekDouble>> rhs(m_bnddim);
        for (int i = 0; i < m_bnddim; ++i)
        {
            rhs[i] = Array<OneD, NekDouble>(m_npoints, 0.);
        }
        // add viscous term
        AddVisPressureBCs(fields, rhs, params);
        // add DuDt term
        AddRigidBodyAcc(rhs, params, nptsPlane0);
        m_BndExp[m_pressure]->NormVectorIProductWRTBase(
            rhs, m_BndExp[m_pressure]->UpdateCoeffs());
    }
    // velocity
    if (m_hasVels && nptsPlane0)
    {
        Array<OneD, Array<OneD, NekDouble>> velocities(m_bnddim);
        for (size_t k = 0; k < m_bnddim; ++k)
        {
            if (m_BndExp.find(k) != m_BndExp.end())
            {
                velocities[k] = Array<OneD, NekDouble>(nptsPlane0, 0.0);
            }
        }
        RigidBodyVelocity(velocities, params, nptsPlane0);
        for (int k = 0; k < m_bnddim; ++k)
        {
            if (m_BndExp.find(k) != m_BndExp.end())
            {
                if (m_BndExp[k]->GetExpType() == MultiRegions::e2DH1D)
                {
                    m_BndExp[k]->GetPlane(0)->FwdTransBndConstrained(
                        velocities[k],
                        m_BndExp[k]->GetPlane(0)->UpdateCoeffs());
                }
                else
                {
                    m_BndExp[k]->FwdTransBndConstrained(
                        velocities[k], m_BndExp[k]->UpdateCoeffs());
                }
            }
        }
    }
}

References Nektar::IncBaseCondition::AddRigidBodyAcc(), Nektar::IncBaseCondition::AddVisPressureBCs(), Nektar::MultiRegions::e2DH1D, Nektar::IncBaseCondition::InitialiseCoords(), Nektar::IncBaseCondition::m_bnddim, Nektar::IncBaseCondition::m_BndExp, m_hasPressure, m_hasVels, Nektar::IncBaseCondition::m_npoints, Nektar::IncBaseCondition::m_numCalls, m_pressure, Nektar::IncBaseCondition::RigidBodyVelocity(), and Nektar::IncBaseCondition::SetNumPointsOnPlane0().

Friends And Related Function Documentation

MemoryManager< MovingFrameWall >

friend class MemoryManager< MovingFrameWall >

friend

Definition at line 1 of file MovingFrameWall.h.

Member Data Documentation

className

std::string Nektar::MovingFrameWall::className

static

Initial value:

=
    GetIncBCFactory().RegisterCreatorFunction(
        "MovingFrameWall", MovingFrameWall::create,
        "Far field boundary condition of moving reference frame")

Definition at line 67 of file MovingFrameWall.h.

m_hasPressure

bool Nektar::MovingFrameWall::m_hasPressure

protected

Definition at line 84 of file MovingFrameWall.h.

Referenced by MovingFrameWall(), v_Initialise(), and v_Update().

m_hasVels

bool Nektar::MovingFrameWall::m_hasVels

protected

Definition at line 83 of file MovingFrameWall.h.

Referenced by MovingFrameWall(), and v_Update().

m_pressure

int Nektar::MovingFrameWall::m_pressure

protected

Definition at line 77 of file MovingFrameWall.h.

Referenced by MovingFrameWall(), and v_Update().