#include <WeakPressureExtrapolate.h>

Inheritance diagram for Nektar::WeakPressureExtrapolate:

Static Public Member Functions
static ExtrapolateSharedPtr	create (const LibUtilities::SessionReaderSharedPtr &pSession, Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, MultiRegions::ExpListSharedPtr &pPressure, const Array< OneD, int > &pVel, const SolverUtils::AdvectionSharedPtr &advObject)
	Creates an instance of this class. More...

Static Public Member Functions inherited from Nektar::StandardExtrapolate
static ExtrapolateSharedPtr	create (const LibUtilities::SessionReaderSharedPtr &pSession, Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, MultiRegions::ExpListSharedPtr &pPressure, const Array< OneD, int > &pVel, const SolverUtils::AdvectionSharedPtr &advObject)
	Creates an instance of this class. More...

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

Static Public Attributes inherited from Nektar::StandardExtrapolate
static std::string	className
	Name of class. More...

Protected Member Functions
	WeakPressureExtrapolate (const LibUtilities::SessionReaderSharedPtr pSession, Array< OneD, MultiRegions::ExpListSharedPtr > pFields, MultiRegions::ExpListSharedPtr pPressure, const Array< OneD, int > pVel, const SolverUtils::AdvectionSharedPtr advObject)

	~WeakPressureExtrapolate () override=default

void	v_EvaluatePressureBCs (const Array< OneD, const Array< OneD, NekDouble > > &fields, const Array< OneD, const Array< OneD, NekDouble > > &N, NekDouble kinvis) override

void	v_MountHOPBCs (int HBCdata, NekDouble kinvis, Array< OneD, NekDouble > &Q, Array< OneD, const NekDouble > &Advection) override

void	v_AddNormVelOnOBC (const int nbcoeffs, const int nreg, Array< OneD, Array< OneD, NekDouble > > &u) override

Protected Member Functions inherited from Nektar::StandardExtrapolate
	StandardExtrapolate (const LibUtilities::SessionReaderSharedPtr pSession, Array< OneD, MultiRegions::ExpListSharedPtr > pFields, MultiRegions::ExpListSharedPtr pPressure, const Array< OneD, int > pVel, const SolverUtils::AdvectionSharedPtr advObject)

	~StandardExtrapolate () override=default

void	v_EvaluatePressureBCs (const Array< OneD, const Array< OneD, NekDouble > > &fields, const Array< OneD, const Array< OneD, NekDouble > > &N, NekDouble kinvis) override

void	v_SubSteppingTimeIntegration (const LibUtilities::TimeIntegrationSchemeSharedPtr &IntegrationScheme) override

void	v_SubStepSaveFields (int nstep) override

void	v_SubStepSetPressureBCs (const Array< OneD, const Array< OneD, NekDouble > > &inarray, NekDouble Aii_DT, NekDouble kinvis) override

void	v_SubStepAdvance (int nstep, NekDouble time) override

void	v_MountHOPBCs (int HBCdata, NekDouble kinvis, Array< OneD, NekDouble > &Q, Array< OneD, const NekDouble > &Advection) override

void	v_AccelerationBDF (Array< OneD, Array< OneD, NekDouble > > &array) override

Protected Member Functions inherited from Nektar::Extrapolate
virtual void	v_EvaluatePressureBCs (const Array< OneD, const Array< OneD, NekDouble > > &inarray, const Array< OneD, const Array< OneD, NekDouble > > &N, NekDouble kinvis)=0

virtual void	v_SubSteppingTimeIntegration (const LibUtilities::TimeIntegrationSchemeSharedPtr &IntegrationScheme)=0

virtual void	v_SubStepSaveFields (int nstep)=0

virtual void	v_SubStepSetPressureBCs (const Array< OneD, const Array< OneD, NekDouble > > &inarray, NekDouble Aii_DT, NekDouble kinvis)=0

virtual void	v_SubStepAdvance (int nstep, NekDouble time)=0

virtual void	v_MountHOPBCs (int HBCdata, NekDouble kinvis, Array< OneD, NekDouble > &Q, Array< OneD, const NekDouble > &Advection)=0

virtual std::string	v_GetSubStepName (void)

virtual void	v_AccelerationBDF (Array< OneD, Array< OneD, NekDouble > > &array)

virtual void	v_CalcNeumannPressureBCs (const Array< OneD, const Array< OneD, NekDouble > > &fields, const Array< OneD, const Array< OneD, NekDouble > > &N, NekDouble kinvis)

virtual void	v_CorrectPressureBCs (const Array< OneD, NekDouble > &pressure)

virtual void	v_AddNormVelOnOBC (const int nbcoeffs, const int nreg, Array< OneD, Array< OneD, NekDouble > > &u)

void	CalcOutflowBCs (const Array< OneD, const Array< OneD, NekDouble > > &fields, NekDouble kinvis)

void	RollOver (Array< OneD, Array< OneD, NekDouble > > &input)

Static Protected Attributes
static std::string	solverTypeLookupId

Static Protected Attributes inherited from Nektar::StandardExtrapolate
static NekDouble	DuDt_Coeffs [3][4]

Static Protected Attributes inherited from Nektar::Extrapolate
static NekDouble	StifflyStable_Betaq_Coeffs [3][3]

static NekDouble	StifflyStable_Alpha_Coeffs [3][3]

static NekDouble	StifflyStable_Gamma0_Coeffs [3] = {1.0, 1.5, 11.0 / 6.0}

Friends
class	MemoryManager< WeakPressureExtrapolate >

Additional Inherited Members
Public Member Functions inherited from Nektar::Extrapolate
	Extrapolate (const LibUtilities::SessionReaderSharedPtr pSession, Array< OneD, MultiRegions::ExpListSharedPtr > pFields, MultiRegions::ExpListSharedPtr pPressure, const Array< OneD, int > pVel, const SolverUtils::AdvectionSharedPtr advObject)

virtual	~Extrapolate ()=default

void	SubSteppingTimeIntegration (const LibUtilities::TimeIntegrationSchemeSharedPtr &IntegrationScheme)

void	SubStepSaveFields (const int nstep)

void	SubStepSetPressureBCs (const Array< OneD, const Array< OneD, NekDouble > > &inarray, const NekDouble Aii_DT, NekDouble kinvis)

void	SubStepAdvance (const int nstep, NekDouble time)

void	MountHOPBCs (int HBCdata, NekDouble kinvis, Array< OneD, NekDouble > &Q, Array< OneD, const NekDouble > &Advection)

void	EvaluatePressureBCs (const Array< OneD, const Array< OneD, NekDouble > > &fields, const Array< OneD, const Array< OneD, NekDouble > > &N, NekDouble kinvis)

void	AddNormVelOnOBC (const int nbcoeffs, const int nreg, Array< OneD, Array< OneD, NekDouble > > &u)

void	CorrectPressureBCs (const Array< OneD, NekDouble > &pressure)

void	CalcNeumannPressureBCs (const Array< OneD, const Array< OneD, NekDouble > > &fields, const Array< OneD, const Array< OneD, NekDouble > > &N, NekDouble kinvis)

std::string	GetSubStepName (void)

void	SetForcing (const std::vector< SolverUtils::ForcingSharedPtr > &forcing)

void	GenerateHOPBCMap (const LibUtilities::SessionReaderSharedPtr &pSsession)

void	UpdateRobinPrimCoeff (void)

void	AddDuDt (void)

void	AddVelBC (void)

void	ExtrapolatePressureHBCs (void)

void	CopyPressureHBCsToPbndExp (void)

Array< OneD, NekDouble >	GetMaxStdVelocity (const Array< OneD, Array< OneD, NekDouble > > inarray)

void	IProductNormVelocityOnHBC (const Array< OneD, const Array< OneD, NekDouble > > &Vel, Array< OneD, NekDouble > &IprodVn)

void	IProductNormVelocityBCOnHBC (Array< OneD, NekDouble > &IprodVn)

void	ExtrapolateArray (Array< OneD, Array< OneD, NekDouble > > &array)

void	EvaluateBDFArray (Array< OneD, Array< OneD, NekDouble > > &array)

void	ExtrapolateArray (Array< OneD, Array< OneD, NekDouble > > &oldarrays, Array< OneD, NekDouble > &newarray, Array< OneD, NekDouble > &outarray)

void	AddPressureToOutflowBCs (NekDouble kinvis)

void	GenerateBndElmtExpansion (void)

Protected Attributes inherited from Nektar::Extrapolate
LibUtilities::SessionReaderSharedPtr	m_session

LibUtilities::CommSharedPtr	m_comm

Array< OneD, HBCType >	m_hbcType
	Array of type of high order BCs for splitting shemes. More...

Array< OneD, MultiRegions::ExpListSharedPtr >	m_fields
	Velocity fields. More...

MultiRegions::ExpListSharedPtr	m_pressure
	Pointer to field holding pressure field. More...

Array< OneD, int >	m_velocity
	int which identifies which components of m_fields contains the velocity (u,v,w); More...

SolverUtils::AdvectionSharedPtr	m_advObject

std::vector< SolverUtils::ForcingSharedPtr >	m_forcing

Array< OneD, Array< OneD, NekDouble > >	m_previousVelFields

int	m_curl_dim
	Curl-curl dimensionality. More...

int	m_bnd_dim
	bounday dimensionality More...

Array< OneD, const SpatialDomains::BoundaryConditionShPtr >	m_PBndConds
	pressure boundary conditions container More...

Array< OneD, MultiRegions::ExpListSharedPtr >	m_PBndExp
	pressure boundary conditions expansion container More...

Array< OneD, MultiRegions::ExpListSharedPtr >	m_bndElmtExps
	Boundary expansions on each domain boundary. More...

int	m_pressureCalls
	number of times the high-order pressure BCs have been called More...

int	m_numHBCDof

int	m_HBCnumber

int	m_intSteps
	Maximum points used in pressure BC evaluation. More...

NekDouble	m_timestep

Array< OneD, Array< OneD, NekDouble > >	m_pressureHBCs
	Storage for current and previous levels of high order pressure boundary conditions. More...

Array< OneD, Array< OneD, NekDouble > >	m_iprodnormvel
	Storage for current and previous levels of the inner product of normal velocity. More...

Array< OneD, Array< OneD, NekDouble > >	m_traceNormals

HighOrderOutflowSharedPtr	m_houtflow

Detailed Description

Definition at line 58 of file WeakPressureExtrapolate.h.

Constructor & Destructor Documentation

◆ WeakPressureExtrapolate()

Nektar::WeakPressureExtrapolate::WeakPressureExtrapolate	(	const LibUtilities::SessionReaderSharedPtr	pSession,
		Array< OneD, MultiRegions::ExpListSharedPtr >	pFields,
		MultiRegions::ExpListSharedPtr	pPressure,
		const Array< OneD, int >	pVel,
		const SolverUtils::AdvectionSharedPtr	advObject
	)

protected

Definition at line 49 of file WeakPressureExtrapolate.cpp.

    : StandardExtrapolate(pSession, pFields, pPressure, pVel, advObject)
{
}

◆ ~WeakPressureExtrapolate()

Nektar::WeakPressureExtrapolate::~WeakPressureExtrapolate ( )

overrideprotecteddefault

Member Function Documentation

◆ create()

static ExtrapolateSharedPtr Nektar::WeakPressureExtrapolate::create	(	const LibUtilities::SessionReaderSharedPtr &	pSession,
		Array< OneD, MultiRegions::ExpListSharedPtr > &	pFields,
		MultiRegions::ExpListSharedPtr &	pPressure,
		const Array< OneD, int > &	pVel,
		const SolverUtils::AdvectionSharedPtr &	advObject
	)

inlinestatic

Creates an instance of this class.

Definition at line 64 of file WeakPressureExtrapolate.h.

    {
        ExtrapolateSharedPtr p =
            MemoryManager<WeakPressureExtrapolate>::AllocateSharedPtr(
                pSession, pFields, pPressure, pVel, advObject);
        return p;
    }

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

◆ v_AddNormVelOnOBC()

void Nektar::WeakPressureExtrapolate::v_AddNormVelOnOBC	(	const int	nbcoeffs,
		const int	nreg,
		Array< OneD, Array< OneD, NekDouble > > &	u
	)

overrideprotectedvirtual

Reimplemented from Nektar::Extrapolate.

Definition at line 95 of file WeakPressureExtrapolate.cpp.

{
    if (!m_houtflow.get()) // no outflow on partition so just return
    {
        return;
    }
 
    int nbcoeffs = m_PBndExp[nreg]->GetNcoeffs();
    // int nqb      = m_PBndExp[nreg]->GetTotPoints();
 
    Array<OneD, NekDouble> IProdVnTmp(nbcoeffs);
 
#if 0
        Array<OneD, Array<OneD, NekDouble> > ubnd(m_curl_dim);
 
        for(int i = 0; i < m_curl_dim; ++i)
        {
            EvaluateBDFArray(m_houtflow->m_outflowVelBnd[noutflow][i]);
 
            ubnd[i] = m_houtflow->m_outflowVelBnd[noutflow][i][m_intSteps-1];
 
            // point input u to the first part of the array for later uee.
            // u[i] = m_houtflow->m_outflowVelBnd[noutflow][i][0];
            u[i] = ubnd[i];
        }
 
        m_PBndExp[nreg]->NormVectorIProductWRTBase(ubnd,IProdVnTmp);
#endif
    m_PBndExp[nreg]->NormVectorIProductWRTBase(u, IProdVnTmp);
 
    Vmath::Svtvp(nbcoeffs, -1.0 / m_timestep, IProdVnTmp, 1,
                 m_PBndExp[nreg]->UpdateCoeffs(), 1,
                 m_PBndExp[nreg]->UpdateCoeffs(), 1);
}

References Nektar::Extrapolate::EvaluateBDFArray(), Nektar::Extrapolate::m_curl_dim, Nektar::Extrapolate::m_houtflow, Nektar::Extrapolate::m_intSteps, Nektar::Extrapolate::m_PBndExp, Nektar::Extrapolate::m_timestep, and Vmath::Svtvp().

◆ v_EvaluatePressureBCs()

void Nektar::WeakPressureExtrapolate::v_EvaluatePressureBCs	(	const Array< OneD, const Array< OneD, NekDouble > > &	fields,
		const Array< OneD, const Array< OneD, NekDouble > > &	N,
		NekDouble	kinvis
	)

overrideprotectedvirtual

Function to extrapolate the new pressure boundary condition. Based on the velocity field and on the advection term. Acceleration term is also computed. This routine is a general one for 2d and 3D application and it can be called directly from velocity correction scheme. Specialisation on dimensionality is redirected to the CalcNeumannPressureBCs method.

Reimplemented from Nektar::StandardExtrapolate.

Definition at line 67 of file WeakPressureExtrapolate.cpp.

{
    m_pressureCalls++;
    if (m_HBCnumber > 0)
    {
        // Calculate just viscous BCs at current level and put in
        // m_pressureHBCs[nlevels-1]
        CalcNeumannPressureBCs(fields, N, kinvis);
 
        // Extrapolate to m_pressureHBCs to n+1
        ExtrapolateArray(m_pressureHBCs);
 
        // \int_bnd q  n.u^{n} ds update current normal of field
        // add m_pressureHBCs to gamma_0/Dt * m_acceleration[0]
        AddVelBC();
 
        // Copy m_pressureHBCs to m_PbndExp
        CopyPressureHBCsToPbndExp();
    }
 
    // Evaluate High order outflow conditions if required.
    CalcOutflowBCs(fields, kinvis);
}

References Nektar::Extrapolate::AddVelBC(), Nektar::Extrapolate::CalcNeumannPressureBCs(), Nektar::Extrapolate::CalcOutflowBCs(), Nektar::Extrapolate::CopyPressureHBCsToPbndExp(), Nektar::Extrapolate::ExtrapolateArray(), Nektar::Extrapolate::m_HBCnumber, Nektar::Extrapolate::m_pressureCalls, and Nektar::Extrapolate::m_pressureHBCs.

◆ v_MountHOPBCs()

void Nektar::WeakPressureExtrapolate::v_MountHOPBCs	(	int	HBCdata,
		NekDouble	kinvis,
		Array< OneD, NekDouble > &	Q,
		Array< OneD, const NekDouble > &	Advection
	)

overrideprotectedvirtual

vritual function which only puts in the curl operator into the bcs

Reimplemented from Nektar::StandardExtrapolate.

Definition at line 135 of file WeakPressureExtrapolate.cpp.

{
    Vmath::Smul(HBCdata, -kinvis, Q, 1, Q, 1);
}

References Vmath::Smul().

Friends And Related Function Documentation

◆ MemoryManager< WeakPressureExtrapolate >

friend class MemoryManager< WeakPressureExtrapolate >

friend

Definition at line 56 of file WeakPressureExtrapolate.h.

Member Data Documentation

◆ className

std::string Nektar::WeakPressureExtrapolate::className

static

Initial value:

=
    GetExtrapolateFactory().RegisterCreatorFunction(
        "WeakPressure", WeakPressureExtrapolate::create, "WeakPressure")

Name of class.

Registers the class with the Factory.

Definition at line 77 of file WeakPressureExtrapolate.h.

◆ solverTypeLookupId

std::string Nektar::WeakPressureExtrapolate::solverTypeLookupId

staticprotected

Initial value:

=
    LibUtilities::SessionReader::RegisterEnumValue("SolverType", "WeakPressure",
                                                   eWeakPressure)

Definition at line 99 of file WeakPressureExtrapolate.h.

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Static Protected Attributes

Friends

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ WeakPressureExtrapolate()

◆ ~WeakPressureExtrapolate()

Member Function Documentation

◆ create()

◆ v_AddNormVelOnOBC()

◆ v_EvaluatePressureBCs()

◆ v_MountHOPBCs()

Friends And Related Function Documentation

◆ MemoryManager< WeakPressureExtrapolate >

Member Data Documentation

◆ className

◆ solverTypeLookupId