Nektar::WallBC Class Reference

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

#include <WallBC.h>

Inheritance diagram for Nektar::WallBC:

Collaboration diagram for Nektar::WallBC:

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 int pSpaceDim, const int bcRegion, const int cnt)
	Creates an instance of this class. More...

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

Protected Member Functions
virtual void	v_Apply (Array< OneD, Array< OneD, NekDouble > > &Fwd, Array< OneD, Array< OneD, NekDouble > > &physarray, const NekDouble &time)
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 int pSpaceDim, const int bcRegion, const int cnt)
	Constructor. More...

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

Friends
class	MemoryManager< WallBC >

Additional Inherited Members
Public Member Functions inherited from Nektar::CFSBndCond
virtual	~CFSBndCond ()
void	Apply (Array< OneD, Array< OneD, NekDouble > > &Fwd, Array< OneD, Array< OneD, NekDouble > > &physarray, const NekDouble &time=0)
	Apply the 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...
int	m_spacedim
	Space dimension. More...
VariableConverterSharedPtr	m_varConv
	Auxiliary object to convert variables. More...
NekDouble	m_gamma
	Parameters of the flow. More...
NekDouble	m_rhoInf
NekDouble	m_pInf
Array< OneD, NekDouble >	m_velInf
int	m_bcRegion
	Id of the boundary region. More...
int	m_offset
	Offset. More...

Detailed Description

Wall boundary conditions for compressible flow problems.

Definition at line 48 of file WallBC.h.

Constructor & Destructor Documentation

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

private

Definition at line 48 of file WallBC.cpp.

    : CFSBndCond(pSession, pFields, pTraceNormals, pSpaceDim, bcRegion, cnt)
{
}

virtual Nektar::WallBC::~WallBC ( void  )

inlineprivatevirtual

Definition at line 85 of file WallBC.h.

{};

Member Function Documentation

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

inlinestatic

Creates an instance of this class.

Definition at line 55 of file WallBC.h.

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

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

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

protectedvirtual

Implements Nektar::CFSBndCond.

Definition at line 58 of file WallBC.cpp.

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

{
    int i;
    int nVariables = physarray.num_elements();

    const Array<OneD, const int> &traceBndMap
        = m_fields[0]->GetTraceBndMap();

    // 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.
        if (nVariables == m_spacedim+3)
        {
            Vmath::Zero(nBCEdgePts, &Fwd[nVariables-1][id2], 1);
        }

        // For 2D/3D, define: v* = v - 2(v.n)n
        Array<OneD, NekDouble> tmp(nBCEdgePts, 0.0);

        // Calculate (v.n)
        for (i = 0; i < m_spacedim; ++i)
        {
            Vmath::Vvtvp(nBCEdgePts,
                         &Fwd[1+i][id2], 1,
                         &m_traceNormals[i][id2], 1,
                         &tmp[0], 1,
                         &tmp[0], 1);
        }

        // Calculate 2.0(v.n)
        Vmath::Smul(nBCEdgePts, -2.0, &tmp[0], 1, &tmp[0], 1);

        // Calculate v* = v - 2.0(v.n)n
        for (i = 0; i < m_spacedim; ++i)
        {
            Vmath::Vvtvp(nBCEdgePts,
                         &tmp[0], 1,
                         &m_traceNormals[i][id2], 1,
                         &Fwd[1+i][id2], 1,
                         &Fwd[1+i][id2], 1);
        }

        // 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);
        }
    }
}

Friends And Related Function Documentation

friend class MemoryManager< WallBC >

friend

Definition at line 52 of file WallBC.h.

Member Data Documentation

std::string Nektar::WallBC::className

static

Initial value:

= GetCFSBndCondFactory().
    RegisterCreatorFunction("Wall",
                            WallBC::create,
                            "Slip wall boundary condition.")

Name of the class.

Definition at line 68 of file WallBC.h.