Nektar::WallViscousBC Class Reference

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

#include <WallViscousBC.h>

Inheritance diagram for Nektar::WallViscousBC:

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. More...

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

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. 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 ()

Protected Attributes
Array< OneD, Array< OneD, NekDouble > >	m_bndPhys
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...
Array< OneD, Array< OneD, NekDouble > >	m_gridVelocityTrace
	Grid Velocity. 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
NekDouble	m_angVel
int	m_bcRegion
	Id of the boundary region. More...
int	m_offset
	Offset. More...

Private Member Functions
	WallViscousBC (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)
	~WallViscousBC (void) override

Friends
class	MemoryManager< WallViscousBC >

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

Detailed Description

Wall boundary conditions for viscous compressible flow problems.

Definition at line 46 of file WallViscousBC.h.

Constructor & Destructor Documentation

WallViscousBC()

Nektar::WallViscousBC::WallViscousBC ( 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 52 of file WallViscousBC.cpp.

    : CFSBndCond(pSession, pFields, pTraceNormals, pGridVelocity, pSpaceDim,
                 bcRegion, cnt)
{
    m_diffusionAveWeight = 0.5;
 
    m_bndPhys = Array<OneD, Array<OneD, NekDouble>>(m_fields.size());
}

References m_bndPhys, Nektar::CFSBndCond::m_diffusionAveWeight, and Nektar::CFSBndCond::m_fields.

~WallViscousBC()

Nektar::WallViscousBC::~WallViscousBC ( void  )

inlineoverrideprivate

Definition at line 85 of file WallViscousBC.h.

{};

Member Function Documentation

create()

static CFSBndCondSharedPtr Nektar::WallViscousBC::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 WallViscousBC.h.

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

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

v_Apply()

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

overrideprotectedvirtual

Implements Nektar::CFSBndCond.

Definition at line 66 of file WallViscousBC.cpp.

{
    int i;
    int nVariables = physarray.size();
 
    // Find the fields whose WallViscous/Adiabatic-BC is time-dependent
    // Update variables on the boundaries of these fields
    // Get the updated variables on the WallViscous/Adiabatic boundary
    //
    // Maybe the EvaluateBoundaryConditions() should be put upstream to
    // CompressibleFlowSystem::NumCalRiemFluxJac(), So that the BCs will not
    // be repeatedly updated when there are more than one time-dependent BC.
    std::string varName;
    for (i = 0; i < nVariables; ++i)
    {
        if (m_fields[i]->GetBndConditions()[m_bcRegion]->IsTimeDependent())
        {
            varName = m_session->GetVariable(i);
            m_fields[i]->EvaluateBoundaryConditions(time, varName);
 
            m_bndPhys[i] =
                m_fields[i]->GetBndCondExpansions()[m_bcRegion]->UpdatePhys();
        }
    }
 
    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.
        if (nVariables == m_spacedim + 3)
        {
            Vmath::Zero(nBCEdgePts, &Fwd[nVariables - 1][id2], 1);
        }
 
        // V = - Vin
        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];
            }
        }
 
        // Superimpose the perturbation
        for (i = 0; i < nVariables; ++i)
        {
            if (m_fields[i]->GetBndConditions()[m_bcRegion]->IsTimeDependent())
            {
                Vmath::Vadd(nBCEdgePts, &m_bndPhys[i][id1], 1, &Fwd[i][id2], 1,
                            &Fwd[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);
        }
    }
}

References Nektar::CFSBndCond::m_bcRegion, m_bndPhys, Nektar::CFSBndCond::m_fields, Nektar::CFSBndCond::m_gridVelocityTrace, Nektar::CFSBndCond::m_offset, Nektar::CFSBndCond::m_session, Nektar::CFSBndCond::m_spacedim, Vmath::Vadd(), Vmath::Vcopy(), and Vmath::Zero().

Friends And Related Function Documentation

MemoryManager< WallViscousBC >

friend class MemoryManager< WallViscousBC >

friend

Definition at line 1 of file WallViscousBC.h.

Member Data Documentation

classNameAdiabatic

std::string Nektar::WallViscousBC::classNameAdiabatic

static

Initial value:

=
    GetCFSBndCondFactory().RegisterCreatorFunction(
        "WallAdiabatic", WallViscousBC::create,
        "Adiabatic wall boundary condition.")

Definition at line 67 of file WallViscousBC.h.

classNameViscous

std::string Nektar::WallViscousBC::classNameViscous

static

Initial value:

=
    GetCFSBndCondFactory().RegisterCreatorFunction(
        "WallViscous", WallViscousBC::create,
        "No-slip (viscous) wall boundary condition.")

Name of the class.

Definition at line 66 of file WallViscousBC.h.

m_bndPhys

Array<OneD, Array<OneD, NekDouble> > Nektar::WallViscousBC::m_bndPhys

protected

Definition at line 72 of file WallViscousBC.h.

Referenced by v_Apply(), and WallViscousBC().