Nektar::IsentropicVortexBC Class Reference

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

#include <IsentropicVortexBC.h>

Inheritance diagram for Nektar::IsentropicVortexBC:

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

Private Member Functions
	IsentropicVortexBC (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)
void	EvaluateIsentropicVortex (const Array< OneD, NekDouble > &x, const Array< OneD, NekDouble > &y, const Array< OneD, NekDouble > &z, Array< OneD, Array< OneD, NekDouble > > &u, NekDouble time, const int o=0)
	~IsentropicVortexBC (void) override

Friends
class	MemoryManager< IsentropicVortexBC >

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...
void	ApplyBwdWeight ()
	Apply the Weight of 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_diffusionAveWeight
	Weight for average calculation of diffusion term. More...
NekDouble	m_gamma
	Parameters of the flow. More...
NekDouble	m_rhoInf
NekDouble	m_pInf
NekDouble	m_pOut
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 46 of file IsentropicVortexBC.h.

Constructor & Destructor Documentation

IsentropicVortexBC()

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

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

~IsentropicVortexBC()

Nektar::IsentropicVortexBC::~IsentropicVortexBC ( void  )

inlineoverrideprivate

Definition at line 85 of file IsentropicVortexBC.h.

{};

Member Function Documentation

create()

static CFSBndCondSharedPtr Nektar::IsentropicVortexBC::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 52 of file IsentropicVortexBC.h.

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

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

EvaluateIsentropicVortex()

void Nektar::IsentropicVortexBC::EvaluateIsentropicVortex ( const Array< OneD, NekDouble > &  x, const Array< OneD, NekDouble > &  y, const Array< OneD, NekDouble > &  z, Array< OneD, Array< OneD, NekDouble > > &  u, NekDouble  time, const int  o = 0  )

private

Definition at line 99 of file IsentropicVortexBC.cpp.

{
    int nq = x.size();
 
    // Flow parameters
    const NekDouble x0    = 5.0;
    const NekDouble y0    = 0.0;
    const NekDouble beta  = 5.0;
    const NekDouble u0    = 1.0;
    const NekDouble v0    = 0.5;
    const NekDouble gamma = m_gamma;
    NekDouble r, xbar, ybar, tmp;
    NekDouble fac = 1.0 / (16.0 * gamma * M_PI * M_PI);
 
    // In 3D zero rhow field.
    if (m_spacedim == 3)
    {
        Vmath::Zero(nq, &u[3][o], 1);
    }
 
    // Fill storage
    for (int i = 0; i < nq; ++i)
    {
        xbar = x[i] - u0 * time - x0;
        ybar = y[i] - v0 * time - y0;
        r    = sqrt(xbar * xbar + ybar * ybar);
        tmp  = beta * exp(1 - r * r);
        u[0][i + o] =
            pow(1.0 - (gamma - 1.0) * tmp * tmp * fac, 1.0 / (gamma - 1.0));
        u[1][i + o] = u[0][i + o] * (u0 - tmp * ybar / (2 * M_PI));
        u[2][i + o] = u[0][i + o] * (v0 + tmp * xbar / (2 * M_PI));
        u[m_spacedim + 1][i + o] =
            pow(u[0][i + o], gamma) / (gamma - 1.0) +
            0.5 * (u[1][i + o] * u[1][i + o] + u[2][i + o] * u[2][i + o]) /
                u[0][i + o];
    }
}

References Nektar::LibUtilities::beta, Nektar::CFSBndCond::m_gamma, Nektar::CFSBndCond::m_spacedim, tinysimd::sqrt(), and Vmath::Zero().

Referenced by v_Apply().

v_Apply()

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

overrideprotectedvirtual

Implements Nektar::CFSBndCond.

Definition at line 56 of file IsentropicVortexBC.cpp.

{
    int nvariables = physarray.size();
 
    const Array<OneD, const int> &bndTraceMap = m_fields[0]->GetTraceBndMap();
    // loop over Boundary Regions
    int npoints, id1, id2, e_max;
 
    e_max = m_fields[0]->GetBndCondExpansions()[m_bcRegion]->GetExpSize();
 
    for (int e = 0; e < e_max; ++e)
    {
        npoints = 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(bndTraceMap[m_offset + e]);
 
        Array<OneD, NekDouble> x(npoints, 0.0);
        Array<OneD, NekDouble> y(npoints, 0.0);
        Array<OneD, NekDouble> z(npoints, 0.0);
 
        m_fields[0]->GetBndCondExpansions()[m_bcRegion]->GetExp(e)->GetCoords(
            x, y, z);
 
        EvaluateIsentropicVortex(x, y, z, Fwd, time, id2);
 
        for (int i = 0; i < nvariables; ++i)
        {
            Vmath::Vcopy(npoints, &Fwd[i][id2], 1,
                         &(m_fields[i]
                               ->GetBndCondExpansions()[m_bcRegion]
                               ->UpdatePhys())[id1],
                         1);
        }
    }
}

References EvaluateIsentropicVortex(), Nektar::CFSBndCond::m_bcRegion, Nektar::CFSBndCond::m_fields, Nektar::CFSBndCond::m_offset, Vmath::Vcopy(), and Nektar::UnitTests::z().

Friends And Related Function Documentation

MemoryManager< IsentropicVortexBC >

friend class MemoryManager< IsentropicVortexBC >

friend

Definition at line 1 of file IsentropicVortexBC.h.

Member Data Documentation

className

std::string Nektar::IsentropicVortexBC::className

static

Initial value:

=
    GetCFSBndCondFactory().RegisterCreatorFunction(
        "IsentropicVortex", IsentropicVortexBC::create,
        "Isentropic vortex boundary condition.")

Name of the class.

Definition at line 65 of file IsentropicVortexBC.h.