Nektar::WallRotationalBC Class Reference

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.

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

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.
virtual	~CFSBndCond ()=default
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.
void	ApplyBwdWeight ()
	Apply the Weight of boundary condition.
Protected Attributes inherited from Nektar::CFSBndCond
LibUtilities::SessionReaderSharedPtr	m_session
	Session reader.
Array< OneD, MultiRegions::ExpListSharedPtr >	m_fields
	Array of fields.
Array< OneD, Array< OneD, NekDouble > >	m_traceNormals
	Trace normals.
Array< OneD, Array< OneD, NekDouble > >	m_gridVelocityTrace
	Grid Velocity.
int	m_spacedim
	Space dimension.
VariableConverterSharedPtr	m_varConv
	Auxiliary object to convert variables.
NekDouble	m_diffusionAveWeight
	Weight for average calculation of diffusion term.
NekDouble	m_gamma
	Parameters of the flow.
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.
int	m_offset
	Offset.

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);
        Array<OneD, NekDouble> z(nBCEdgePts, 0.0);
        m_fields[0]->GetBndCondExpansions()[m_bcRegion]->GetExp(e)->GetCoords(
            x, y, z);
        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().

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 82 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 Symbol 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.