#include <WeakPressureExtrapolate.h>

Inheritance diagram for Nektar::WeakPressureExtrapolate:

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

virtual	~WeakPressureExtrapolate ()

Public 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)

virtual	~StandardExtrapolate ()

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

void	GenerateHOPBCMap (const LibUtilities::SessionReaderSharedPtr &pSsession)

void	UpdateRobinPrimCoeff (void)

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	SetForcing (const std::vector< SolverUtils::ForcingSharedPtr > &forcing)

void	AddDuDt (void)

void	AddVelBC (void)

void	ExtrapolatePressureHBCs (void)

void	CopyPressureHBCsToPbndExp (void)

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

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

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

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

std::string	GetSubStepName (void)

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	AddNormVelOnOBC (const int nbcoeffs, const int nreg, Array< OneD, Array< OneD, NekDouble >> &u)

void	AddPressureToOutflowBCs (NekDouble kinvis)

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
virtual void	v_EvaluatePressureBCs (const Array< OneD, const Array< OneD, NekDouble >> &fields, const Array< OneD, const Array< OneD, NekDouble >> &N, NekDouble kinvis) override

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

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

Protected Member Functions inherited from Nektar::StandardExtrapolate
virtual void	v_SubSteppingTimeIntegration (const LibUtilities::TimeIntegrationSchemeSharedPtr &IntegrationScheme) override

virtual void	v_SubStepSaveFields (int nstep) override

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

virtual void	v_SubStepAdvance (int nstep, NekDouble time) override

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

Protected Member Functions inherited from Nektar::Extrapolate
virtual std::string	v_GetSubStepName (void)

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

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)

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

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

Additional Inherited Members
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...

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

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}

Detailed Description

Definition at line 57 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
	)

Definition at line 47 of file WeakPressureExtrapolate.cpp.

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

◆ ~WeakPressureExtrapolate()

Nektar::WeakPressureExtrapolate::~WeakPressureExtrapolate ( )

virtual

Definition at line 56 of file WeakPressureExtrapolate.cpp.

57 {

58 }

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 61 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 97 of file WeakPressureExtrapolate.cpp.

 {
     boost::ignore_unused(noutflow);
  
     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 69 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 138 of file WeakPressureExtrapolate.cpp.

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

References Vmath::Smul().

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 74 of file WeakPressureExtrapolate.h.

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ WeakPressureExtrapolate()

◆ ~WeakPressureExtrapolate()

Member Function Documentation

◆ create()

◆ v_AddNormVelOnOBC()

◆ v_EvaluatePressureBCs()

◆ v_MountHOPBCs()

Member Data Documentation

◆ className