Nektar::MRFWallPressDecomp Class Reference

#include <MRFWallPressDecomp.h>

Inheritance diagram for Nektar::MRFWallPressDecomp:

Public Member Functions
	~MRFWallPressDecomp () override
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)
	MRFWallPressDecomp (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	AddExtrapCentVisPressureBCs (const Array< OneD, const Array< OneD, NekDouble > > &fields, Array< OneD, Array< OneD, NekDouble > > &N, std::map< std::string, NekDouble > &params, int npts0)
void	AddCentripetalAcc (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
bool	m_hasVels
bool	m_hasPressure
int	m_pressure
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< MRFWallPressDecomp >

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 49 of file MRFWallPressDecomp.h.

Constructor & Destructor Documentation

~MRFWallPressDecomp()

Nektar::MRFWallPressDecomp::~MRFWallPressDecomp ( )

override

Definition at line 74 of file MRFWallPressDecomp.cpp.

{
}

MRFWallPressDecomp()

Nektar::MRFWallPressDecomp::MRFWallPressDecomp ( 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 47 of file MRFWallPressDecomp.cpp.

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

AddCentripetalAcc()

void Nektar::MRFWallPressDecomp::AddCentripetalAcc ( Array< OneD, Array< OneD, NekDouble > > &  N, std::map< std::string, NekDouble > &  params, int  npts0  )

protected

Definition at line 169 of file MRFWallPressDecomp.cpp.

{
    if (npts0 == 0 || params.find("Omega_z") == params.end())
    {
        return;
    }
    // add centripetal acceleration
    NekDouble Omega = params["Omega_z"];
    NekDouble Wz2   = Omega * Omega;
    Vmath::Svtvp(npts0, Wz2, m_coords[0], 1, N[0], 1, N[0], 1);
    Vmath::Svtvp(npts0, Wz2, m_coords[1], 1, N[1], 1, N[1], 1);
    if (params.find("U") != params.end())
    {
        NekDouble mOmegaU0 = -Omega * params["U"];
        Vmath::Sadd(npts0, mOmegaU0, N[1], 1, N[1], 1);
    }
    if (params.find("V") != params.end())
    {
        NekDouble OmegaV0 = Omega * params["V"];
        Vmath::Sadd(npts0, OmegaV0, N[0], 1, N[0], 1);
    }
}

References Nektar::IncBaseCondition::m_coords, Omega, Vmath::Sadd(), and Vmath::Svtvp().

Referenced by AddExtrapCentVisPressureBCs().

AddExtrapCentVisPressureBCs()

void Nektar::MRFWallPressDecomp::AddExtrapCentVisPressureBCs ( 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 194 of file MRFWallPressDecomp.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);
    AddCentripetalAcc(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 AddCentripetalAcc(), 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::MRFWallPressDecomp::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 54 of file MRFWallPressDecomp.h.

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

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

v_Initialise()

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

overrideprotectedvirtual

Reimplemented from Nektar::IncBaseCondition.

Definition at line 78 of file MRFWallPressDecomp.cpp.

{
    IncBaseCondition::v_Initialise(pSession);
    if (!(pSession->GetSolverInfo("SolverType") == "PressDecompVCSFSI"))
    {
        ASSERTL0(false,
                 "The boundary condition MRFWallPressDecomp is only "
                 "supported for use in the SolverType 'PressDecompVCSFSI'");
    }
    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::MRFWallPressDecomp::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 108 of file MRFWallPressDecomp.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, do if not define velocity
    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 and centripetal acceleration
        AddExtrapCentVisPressureBCs(fields, rhs, params, nptsPlane0);
        m_BndExp[m_pressure]->NormVectorIProductWRTBase(
            rhs, m_BndExp[m_pressure]->UpdateCoeffs());
    }
    // velocity, do if define 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 AddExtrapCentVisPressureBCs(), 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< MRFWallPressDecomp >

friend class MemoryManager< MRFWallPressDecomp >

friend

Definition at line 1 of file MRFWallPressDecomp.h.

Member Data Documentation

className

std::string Nektar::MRFWallPressDecomp::className

static

Initial value:

=
    GetIncBCFactory().RegisterCreatorFunction(
        "MRFWallPressDecomp", MRFWallPressDecomp::create,
        "Wall boundary condition of moving reference frame with pressure "
        "decomposition")

Definition at line 68 of file MRFWallPressDecomp.h.

m_hasPressure

bool Nektar::MRFWallPressDecomp::m_hasPressure

protected

Definition at line 90 of file MRFWallPressDecomp.h.

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

m_hasVels

bool Nektar::MRFWallPressDecomp::m_hasVels

protected

Definition at line 89 of file MRFWallPressDecomp.h.

Referenced by MRFWallPressDecomp(), and v_Update().

m_pressure

int Nektar::MRFWallPressDecomp::m_pressure

protected

Definition at line 91 of file MRFWallPressDecomp.h.

Referenced by MRFWallPressDecomp(), and v_Update().