Nektar::MRFWall Class Reference

#include <MRFWall.h>

Inheritance diagram for Nektar::MRFWall:

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

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)
	MRFWall (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)
void	AddExtrapAcceVisPressureBCs (const Array< OneD, const Array< OneD, NekDouble > > &fields, Array< OneD, Array< OneD, NekDouble > > &N, std::map< std::string, NekDouble > &params, int npts0)
Protected Member Functions inherited from Nektar::StaticWall
	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)
	~StaticWall () override=default
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
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)
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	AddExtrapVisPressureBCs (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
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_extrapArray
int	m_pressure
std::string	classname

Friends
class	MemoryManager< MRFWall >

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 MRFWall.h.

Constructor & Destructor Documentation

~MRFWall()

Nektar::MRFWall::~MRFWall ( )

overridedefault

MRFWall()

Nektar::MRFWall::MRFWall ( 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 MRFWall.cpp.

    : StaticWall(pSession, pFields, cond, exp, nbnd, spacedim, bnddim)
{
    classname = "MRFWall";
    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

AddExtrapAcceVisPressureBCs()

void Nektar::MRFWall::AddExtrapAcceVisPressureBCs ( const Array< OneD, const Array< OneD, NekDouble > > &  fields, Array< OneD, Array< OneD, NekDouble > > &  N, std::map< std::string, NekDouble > &  params, int  npts0  )

protected

Definition at line 160 of file MRFWall.cpp.

{
    for (int i = 0; i < m_bnddim; ++i)
    {
        Vmath::Zero(m_npoints, m_extrapArray[m_intSteps - 1][i], 1);
    }
    AddVisPressureBCs(fields, m_extrapArray[m_intSteps - 1], params);
    AddRigidBodyAcc(m_extrapArray[m_intSteps - 1], params, npts0);
    ExtrapolateArray(m_numCalls, m_extrapArray);
    for (int i = 0; i < m_bnddim; i++)
    {
        Vmath::Vadd(m_npoints, m_extrapArray[m_intSteps - 1][i], 1, N[i], 1,
                    N[i], 1);
    }
}

References Nektar::IncBaseCondition::AddRigidBodyAcc(), Nektar::IncBaseCondition::AddVisPressureBCs(), Nektar::IncBaseCondition::ExtrapolateArray(), Nektar::IncBaseCondition::m_bnddim, Nektar::IncBaseCondition::m_extrapArray, Nektar::IncBaseCondition::m_intSteps, Nektar::IncBaseCondition::m_npoints, Nektar::IncBaseCondition::m_numCalls, Vmath::Vadd(), and Vmath::Zero().

Referenced by v_Update().

create()

static IncBaseConditionSharedPtr Nektar::MRFWall::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 MRFWall.h.

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

References Nektar::MemoryManager< DataType >::AllocateSharedPtr().

v_Initialise()

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

overrideprotectedvirtual

Reimplemented from Nektar::IncBaseCondition.

Definition at line 72 of file MRFWall.cpp.

{
    IncBaseCondition::v_Initialise(pSession);
    if (!(pSession->GetSolverInfo("SolverType") == "VCSFSI"))
    {
        ASSERTL0(false, "The boundary condition MRFWall is only "
                        "supported for use in the SolverType 'VCSFSI'");
    }
    m_field->GetBndElmtExpansion(m_nbnd, m_bndElmtExps, false);
    if (m_hasPressure)
    {
        m_extrapArray =
            Array<OneD, Array<OneD, Array<OneD, NekDouble>>>(m_intSteps);
        for (int n = 0; n < m_intSteps; ++n)
        {
            m_extrapArray[n] = Array<OneD, Array<OneD, NekDouble>>(m_bnddim);
            for (int i = 0; i < m_bnddim; ++i)
            {
                m_extrapArray[n][i] = Array<OneD, NekDouble>(m_npoints, 0.0);
            }
        }
    }
}

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

v_Update()

void Nektar::MRFWall::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 MRFWall.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 (params.find("pressure") != params.end() && 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
        AddExtrapAcceVisPressureBCs(fields, rhs, params, nptsPlane0);
        m_BndExp[m_pressure]->NormVectorIProductWRTBase(
            rhs, m_BndExp[m_pressure]->UpdateCoeffs());
    }
    // velocity
    if (params.find("velocity") != params.end() && 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 AddExtrapAcceVisPressureBCs(), 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 Symbol Documentation

MemoryManager< MRFWall >

friend class MemoryManager< MRFWall >

friend

Definition at line 1 of file MRFWall.h.

Member Data Documentation

className

std::string Nektar::MRFWall::className

static

Initial value:

= GetIncBCFactory().RegisterCreatorFunction(
    "MRFWall", MRFWall::create,
    "Wall boundary condition of moving reference frame")

Definition at line 66 of file MRFWall.h.

m_hasPressure

bool Nektar::MRFWall::m_hasPressure

protected

Definition at line 87 of file MRFWall.h.

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

m_hasVels

bool Nektar::MRFWall::m_hasVels

protected

Definition at line 86 of file MRFWall.h.

Referenced by MRFWall(), and v_Update().

m_pressure

int Nektar::MRFWall::m_pressure

protected

Definition at line 76 of file MRFWall.h.

Referenced by MRFWall(), and v_Update().