#include <StandardExtrapolate.h>

Inheritance diagram for Nektar::StandardExtrapolate:

Collaboration diagram for Nektar::StandardExtrapolate:

Public Member Functions
	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 ()

void	SubSteppingTimeIntegration (const int intMethod, const LibUtilities::TimeIntegrationWrapperSharedPtr &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 LibUtilities::TimeIntegrationSolutionSharedPtr &integrationSoln, 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	CalcExplicitDuDt (const Array< OneD, const Array< OneD, NekDouble > > &fields)

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)

LibUtilities::TimeIntegrationMethod	GetSubStepIntegrationMethod (void)

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

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

virtual void	v_SubSteppingTimeIntegration (int intMethod, const LibUtilities::TimeIntegrationWrapperSharedPtr &IntegrationScheme)

virtual void	v_SubStepSaveFields (int nstep)

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

virtual void	v_SubStepAdvance (const LibUtilities::TimeIntegrationSolutionSharedPtr &integrationSoln, int nstep, NekDouble time)

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

Protected Member Functions inherited from Nektar::Extrapolate
virtual LibUtilities::TimeIntegrationMethod	v_GetSubStepIntegrationMethod (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)

void	CurlCurl (Array< OneD, Array< OneD, const NekDouble > > &Vel, Array< OneD, Array< OneD, NekDouble > > &Q, const int j)

Additional Inherited Members
Protected Attributes inherited from Nektar::Extrapolate
LibUtilities::SessionReaderSharedPtr	m_session

LibUtilities::CommSharedPtr	m_comm

Array< OneD, MultiRegions::ExpListSharedPtr >	m_fields

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

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_pressureBCsMaxPts
	Maximum points used in pressure BC evaluation. More...

int	m_pressureBCsElmtMaxPts
	Maximum points used in Element adjacent to pressure BC evaluation. More...

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

NekDouble	m_timestep

bool	m_SingleMode
	Flag to determine if single homogeneous mode is used. More...

bool	m_HalfMode
	Flag to determine if half homogeneous mode is used. More...

bool	m_MultipleModes
	Flag to determine if use multiple homogenenous modes are used. More...

NekDouble	m_LhomZ
	physical length in Z direction (if homogeneous) More...

int	m_npointsX
	number of points in X direction (if homogeneous) More...

int	m_npointsY
	number of points in Y direction (if homogeneous) More...

int	m_npointsZ
	number of points in Z direction (if homogeneous) More...

Array< OneD, int >	m_pressureBCtoElmtID
	Id of element to which pressure boundary condition belongs. More...

Array< OneD, int >	m_pressureBCtoTraceID
	Id of edge (2D) or face (3D) to which pressure boundary condition belongs. More...

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_acceleration
	Storage for current and previous levels of the acceleration term. More...

Array< OneD, HBCInfo >	m_HBCdata
	data structure to old all the information regarding High order pressure BCs More...

Array< OneD, NekDouble >	m_wavenumber
	wave number 2 pi k /Lz More...

Array< OneD, NekDouble >	m_negWavenumberSq
	minus Square of wavenumber More...

Array< OneD, Array< OneD, Array< OneD, NekDouble > > >	m_outflowVel
	Storage for current and previous velocity fields at the otuflow for high order outflow BCs. More...

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

Array< OneD, Array< OneD, Array< OneD, NekDouble > > >	m_PhyoutfVel
	Storage for current and previous velocity fields in physical space at the otuflow for high order outflow BCs. More...

Array< OneD, NekDouble >	m_nonlinearterm_phys
	(if homogeneous) More...

Array< OneD, NekDouble >	m_nonlinearterm_coeffs
	(if homogeneous) More...

Array< OneD, unsigned int >	m_expsize_per_plane
	(if homogeneous) More...

Array< OneD, NekDouble >	m_PBndCoeffs
	(if homogeneous) More...

Array< OneD, Array< OneD, NekDouble > >	m_UBndCoeffs
	(if homogeneous) More...

int	m_totexps_per_plane
	(if homogeneous) More...

Static Protected Attributes inherited from Nektar::Extrapolate
static NekDouble	StifflyStable_Betaq_Coeffs [3][3]
	total number of expansion for each plane (if homogeneous) More...

static NekDouble	StifflyStable_Alpha_Coeffs [3][3]

static NekDouble	StifflyStable_Gamma0_Coeffs [3]

Detailed Description

Definition at line 58 of file StandardExtrapolate.h.

Constructor & Destructor Documentation

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

Definition at line 49 of file StandardExtrapolate.cpp.

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

Nektar::StandardExtrapolate::~StandardExtrapolate ( )

virtual

Definition at line 59 of file StandardExtrapolate.cpp.

60 {

61 }

Member Function Documentation

static ExtrapolateSharedPtr Nektar::StandardExtrapolate::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 63 of file StandardExtrapolate.h.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr().

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

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

protectedvirtual

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.

Implements Nektar::Extrapolate.

Definition at line 73 of file StandardExtrapolate.cpp.

References Nektar::Extrapolate::CalcExplicitDuDt(), Nektar::Extrapolate::CalcNeumannPressureBCs(), Nektar::Extrapolate::CalcOutflowBCs(), Nektar::Extrapolate::CopyPressureHBCsToPbndExp(), Nektar::Extrapolate::ExtrapolatePressureHBCs(), Nektar::Extrapolate::m_HBCdata, Nektar::Extrapolate::m_pressureCalls, Nektar::Extrapolate::m_pressureHBCs, and Nektar::Extrapolate::RollOver().

     {
         if(m_HBCdata.num_elements()>0)
         {
             m_pressureCalls++;
 
             // Rotate HOPBCs storage
             RollOver(m_pressureHBCs);
             
             // Calculate non-linear and viscous BCs at current level
             // and put in m_pressureHBCs[0]
             CalcNeumannPressureBCs(fields,N,kinvis);
             
             // calculate (phi,du/dt) and level n which will then be
             // extrpolated.
             CalcExplicitDuDt(fields);
             
             // Extrapolate to m_pressureHBCs to n+1
             ExtrapolatePressureHBCs();
             
             // Copy m_pressureHBCs to m_PbndExp
             CopyPressureHBCsToPbndExp();            
 
             // Evaluate High order outflow conditiosn if required. 
             CalcOutflowBCs(fields, kinvis);
         }
     }

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

protectedvirtual

Implements Nektar::Extrapolate.

Definition at line 165 of file StandardExtrapolate.cpp.

References Vmath::Svtvp().

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

void Nektar::StandardExtrapolate::v_SubStepAdvance	(	const LibUtilities::TimeIntegrationSolutionSharedPtr &	integrationSoln,
		int	nstep,
		NekDouble	time
	)

protectedvirtual

Implements Nektar::Extrapolate.

Definition at line 146 of file StandardExtrapolate.cpp.

150 {

151 }

void Nektar::StandardExtrapolate::v_SubSteppingTimeIntegration	(	int	intMethod,
		const LibUtilities::TimeIntegrationWrapperSharedPtr &	IntegrationScheme
	)

protectedvirtual

Implements Nektar::Extrapolate.

Definition at line 108 of file StandardExtrapolate.cpp.

References Nektar::LibUtilities::eIMEXOrder1, Nektar::LibUtilities::eIMEXOrder2, Nektar::LibUtilities::eIMEXOrder3, and Nektar::Extrapolate::m_intSteps.

     {
         switch(intMethod)
         {
             case LibUtilities::eIMEXOrder1:
             {
                 m_intSteps = 1; 
             }
             break;
             case LibUtilities::eIMEXOrder2:
             {
                 m_intSteps = 2;
             }
             break;
             case LibUtilities::eIMEXOrder3:
             {
                 m_intSteps = 3;
             }
             break;
         }
     }

void Nektar::StandardExtrapolate::v_SubStepSaveFields ( int nstep )

protectedvirtual

Implements Nektar::Extrapolate.

Definition at line 157 of file StandardExtrapolate.cpp.

159 {

160 }

void Nektar::StandardExtrapolate::v_SubStepSetPressureBCs	(	const Array< OneD, const Array< OneD, NekDouble > > &	inarray,
		NekDouble	Aii_DT,
		NekDouble	kinvis
	)

protectedvirtual

Implements Nektar::Extrapolate.

Definition at line 135 of file StandardExtrapolate.cpp.

139 {

140 }

Member Data Documentation

std::string Nektar::StandardExtrapolate::className

static

Initial value:

= GetExtrapolateFactory().RegisterCreatorFunction(
        "Standard",
        StandardExtrapolate::create,
        "Standard")

Name of class.

Registers the class with the Factory.

Definition at line 76 of file StandardExtrapolate.h.

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation