Defines a forcing term to be explicitly applied. More...

#include <Forcing.h>

Inheritance diagram for Nektar::SolverUtils::Forcing:

Public Member Functions
virtual SOLVER_UTILS_EXPORT	~Forcing ()

SOLVER_UTILS_EXPORT void	InitObject (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const unsigned int &pNumForcingFields, const TiXmlElement *pForce)
	Initialise the forcing object. More...

SOLVER_UTILS_EXPORT void	Apply (const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble &time)
	Apply the forcing. More...

SOLVER_UTILS_EXPORT const Array< OneD, const Array< OneD, NekDouble > > &	GetForces ()

SOLVER_UTILS_EXPORT Array< OneD, Array< OneD, NekDouble > > &	UpdateForces ()

Static Public Member Functions
static SOLVER_UTILS_EXPORT std::vector< ForcingSharedPtr >	Load (const LibUtilities::SessionReaderSharedPtr &pSession, const std::weak_ptr< EquationSystem > &pEquation, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const unsigned int &pNumForcingFields=0)

Protected Member Functions
SOLVER_UTILS_EXPORT	Forcing (const LibUtilities::SessionReaderSharedPtr &pSession, const std::weak_ptr< EquationSystem > &pEquation)
	Constructor. More...

virtual SOLVER_UTILS_EXPORT void	v_InitObject (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const unsigned int &pNumForcingFields, const TiXmlElement *pForce)=0

virtual SOLVER_UTILS_EXPORT void	v_Apply (const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble &time)=0

SOLVER_UTILS_EXPORT SessionFunctionSharedPtr	GetFunction (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const LibUtilities::SessionReaderSharedPtr &pSession, std::string pName, bool pCache=false)
	Get a SessionFunction by name. More...

SOLVER_UTILS_EXPORT void	EvaluateTimeFunction (LibUtilities::SessionReaderSharedPtr pSession, std::string pFieldName, Array< OneD, NekDouble > &pArray, const std::string &pFunctionName, NekDouble pTime=NekDouble(0))

SOLVER_UTILS_EXPORT void	EvaluateTimeFunction (const NekDouble pTime, const LibUtilities::EquationSharedPtr &pEqn, Array< OneD, NekDouble > &pArray)

Protected Attributes
LibUtilities::SessionReaderSharedPtr	m_session
	Session reader. More...

const std::weak_ptr< EquationSystem >	m_equ
	Weak pointer to equation system using this forcing. More...

Array< OneD, Array< OneD, NekDouble > >	m_Forcing
	Evaluated forcing function. More...

int	m_NumVariable
	Number of variables. More...

std::map< std::string, SolverUtils::SessionFunctionSharedPtr >	m_sessionFunctions
	Map of known SessionFunctions. More...

Detailed Description

Defines a forcing term to be explicitly applied.

Definition at line 72 of file Forcing.h.

Constructor & Destructor Documentation

◆ ~Forcing()

virtual SOLVER_UTILS_EXPORT Nektar::SolverUtils::Forcing::~Forcing ( )

inlinevirtual

Definition at line 75 of file Forcing.h.

References Apply(), GetForces(), InitObject(), Load(), SOLVER_UTILS_EXPORT, and UpdateForces().

75 {}

◆ Forcing()

Nektar::SolverUtils::Forcing::Forcing	(	const LibUtilities::SessionReaderSharedPtr &	pSession,
		const std::weak_ptr< EquationSystem > &	pEquation
	)

protected

Constructor.

Definition at line 50 of file Forcing.cpp.

                 : m_session(pSession), m_equ(pEquation)
         {
 
         }

Member Function Documentation

◆ Apply()

void Nektar::SolverUtils::Forcing::Apply	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		Array< OneD, Array< OneD, NekDouble > > &	outarray,
		const NekDouble &	time
	)

Apply the forcing.

Parameters

fields	Expansion lists corresponding to input arrays
inarray	u^n from previous timestep
outarray	output array to append forcing to

Definition at line 72 of file Forcing.cpp.

References v_Apply().

Referenced by ~Forcing().

         {
             v_Apply(fields, inarray, outarray, time);
         }

◆ EvaluateTimeFunction() [1/2]

void Nektar::SolverUtils::Forcing::EvaluateTimeFunction	(	LibUtilities::SessionReaderSharedPtr	pSession,
		std::string	pFieldName,
		Array< OneD, NekDouble > &	pArray,
		const std::string &	pFunctionName,
		NekDouble	pTime = `NekDouble(0)`
	)

protected

Definition at line 131 of file Forcing.cpp.

References ASSERTL0.

Referenced by Nektar::SolverUtils::ForcingBody::v_Apply(), and Nektar::SolverUtils::ForcingAbsorption::v_Apply().

         {
             ASSERTL0(pSession->DefinesFunction(pFunctionName),
                      "Function '" + pFunctionName + "' does not exist.");
             
             LibUtilities::EquationSharedPtr ffunc =
                 pSession->GetFunction(pFunctionName, pFieldName);
             
             EvaluateTimeFunction(pTime,ffunc,pArray);
         }

◆ EvaluateTimeFunction() [2/2]

void Nektar::SolverUtils::Forcing::EvaluateTimeFunction	(	const NekDouble	pTime,
		const LibUtilities::EquationSharedPtr &	pEqn,
		Array< OneD, NekDouble > &	pArray
	)

protected

Definition at line 148 of file Forcing.cpp.

         {
             // dummy array of zero pts.
             Array<OneD, NekDouble> x0(pArray.num_elements(),0.0);
 
             pEqn->Evaluate(x0, x0, x0, pTime, pArray);
         }

◆ GetForces()

const Nektar::Array< OneD, const Array< OneD, NekDouble > > & Nektar::SolverUtils::Forcing::GetForces ( )

Definition at line 126 of file Forcing.cpp.

References m_Forcing.

Referenced by ~Forcing().

         {
             return m_Forcing;
         }

◆ GetFunction()

SessionFunctionSharedPtr Nektar::SolverUtils::Forcing::GetFunction	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	pFields,
		const LibUtilities::SessionReaderSharedPtr &	pSession,
		std::string	pName,
		bool	pCache = `false`
	)

protected

Get a SessionFunction by name.

Definition at line 159 of file Forcing.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), and m_sessionFunctions.

Referenced by Nektar::SolverUtils::ForcingAbsorption::CalcAbsorption(), Nektar::SolverUtils::ForcingBody::Update(), Nektar::ForcingMovingBody::v_Apply(), Nektar::ForcingQuasi1D::v_InitObject(), and Nektar::SolverUtils::ForcingAbsorption::v_InitObject().

         {
             if (pCache)
             {
                 if ((m_sessionFunctions.find(pName) == m_sessionFunctions.end())
                     || (m_sessionFunctions[pName]->GetSession() != pSession)
                     || (m_sessionFunctions[pName]->GetExpansion() != pFields[0])
                 )
                 {
                     m_sessionFunctions[pName] =
                         MemoryManager<SessionFunction>::AllocateSharedPtr(
                             pSession, pFields[0], pName, pCache);
                 }
 
                 return m_sessionFunctions[pName];
             }
             else
             {
                 return SessionFunctionSharedPtr(new SessionFunction(pSession, pFields[0], pName, pCache));
             }
         }

◆ InitObject()

void Nektar::SolverUtils::Forcing::InitObject	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	pFields,
		const unsigned int &	pNumForcingFields,
		const TiXmlElement *	pForce
	)

Initialise the forcing object.

Definition at line 58 of file Forcing.cpp.

References v_InitObject().

Referenced by ~Forcing().

         {
             v_InitObject(pFields, pNumForcingFields, pForce);
         }

◆ Load()

vector< ForcingSharedPtr > Nektar::SolverUtils::Forcing::Load	(	const LibUtilities::SessionReaderSharedPtr &	pSession,
		const std::weak_ptr< EquationSystem > &	pEquation,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	pFields,
		const unsigned int &	pNumForcingFields = `0`
	)

static

Definition at line 85 of file Forcing.cpp.

References Nektar::SolverUtils::GetForcingFactory().

Referenced by Nektar::UnsteadyReactionDiffusion::v_InitObject(), Nektar::Dummy::v_InitObject(), Nektar::AcousticSystem::v_InitObject(), Nektar::CompressibleFlowSystem::v_InitObject(), Nektar::UnsteadyViscousBurgers::v_InitObject(), Nektar::IncNavierStokes::v_InitObject(), and ~Forcing().

         {
             vector<ForcingSharedPtr> vForceList;
 
             if (!pSession->DefinesElement("Nektar/Forcing"))
             {
                 return vForceList;
             }
 
             TiXmlElement* vForcing = pSession->GetElement("Nektar/Forcing");
             if (vForcing)
             {
                 unsigned int vNumForcingFields = pNumForcingFields;
                 if (!pNumForcingFields)
                 {
                     vNumForcingFields = pFields.num_elements();
                 }
 
                 TiXmlElement* vForce = vForcing->FirstChildElement("FORCE");
                 while (vForce)
                 {
                     string vType = vForce->Attribute("TYPE");
 
                     vForceList.push_back(GetForcingFactory().CreateInstance(
                                         vType, pSession, pEquation, pFields,
                                         vNumForcingFields, vForce));
                     vForce = vForce->NextSiblingElement("FORCE");
                 }
             }
             return vForceList;
         }

◆ UpdateForces()

Nektar::Array< OneD, Array< OneD, NekDouble > > & Nektar::SolverUtils::Forcing::UpdateForces ( )

Definition at line 121 of file Forcing.cpp.

References m_Forcing.

Referenced by ~Forcing().

         {
             return m_Forcing;
         }

◆ v_Apply()

virtual SOLVER_UTILS_EXPORT void Nektar::SolverUtils::Forcing::v_Apply	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		Array< OneD, Array< OneD, NekDouble > > &	outarray,
		const NekDouble &	time
	)

protectedpure virtual

Implemented in Nektar::SolverUtils::ForcingAbsorption, Nektar::SolverUtils::ForcingProgrammatic, Nektar::ForcingMovingBody, Nektar::SolverUtils::ForcingBody, Nektar::SolverUtils::ForcingMovingReferenceFrame, Nektar::SolverUtils::ForcingNoise, Nektar::ForcingStabilityCoupledLNS, Nektar::ForcingAxiSymmetric, and Nektar::ForcingQuasi1D.

Referenced by Apply().

◆ v_InitObject()

virtual SOLVER_UTILS_EXPORT void Nektar::SolverUtils::Forcing::v_InitObject	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	pFields,
		const unsigned int &	pNumForcingFields,
		const TiXmlElement *	pForce
	)

protectedpure virtual

Implemented in Nektar::SolverUtils::ForcingAbsorption, Nektar::SolverUtils::ForcingProgrammatic, Nektar::ForcingMovingBody, Nektar::SolverUtils::ForcingBody, Nektar::SolverUtils::ForcingMovingReferenceFrame, Nektar::SolverUtils::ForcingNoise, Nektar::ForcingStabilityCoupledLNS, Nektar::ForcingAxiSymmetric, and Nektar::ForcingQuasi1D.

Referenced by InitObject().

Member Data Documentation

◆ m_equ

const std::weak_ptr<EquationSystem> Nektar::SolverUtils::Forcing::m_equ

protected

Weak pointer to equation system using this forcing.

Definition at line 105 of file Forcing.h.

Referenced by Nektar::ForcingMovingBody::InitialiseFilter().

◆ m_Forcing

Array<OneD, Array<OneD, NekDouble> > Nektar::SolverUtils::Forcing::m_Forcing

protected

Evaluated forcing function.

Definition at line 107 of file Forcing.h.

◆ m_NumVariable

int Nektar::SolverUtils::Forcing::m_NumVariable

protected

Number of variables.

Definition at line 109 of file Forcing.h.

◆ m_session

LibUtilities::SessionReaderSharedPtr Nektar::SolverUtils::Forcing::m_session

protected

Session reader.

Definition at line 103 of file Forcing.h.

Referenced by Nektar::SolverUtils::ForcingAbsorption::CalcAbsorption(), Nektar::ForcingMovingBody::CheckIsFromFile(), Nektar::ForcingMovingBody::EvaluateStructDynModel(), Nektar::ForcingMovingBody::InitialiseCableModel(), Nektar::ForcingMovingBody::Newmark_betaSolver(), Nektar::ForcingMovingBody::SetDynEqCoeffMatrix(), Nektar::SolverUtils::ForcingBody::Update(), Nektar::ForcingMovingBody::v_Apply(), Nektar::SolverUtils::ForcingAbsorption::v_Apply(), Nektar::ForcingQuasi1D::v_InitObject(), Nektar::ForcingAxiSymmetric::v_InitObject(), Nektar::SolverUtils::ForcingNoise::v_InitObject(), Nektar::SolverUtils::ForcingMovingReferenceFrame::v_InitObject(), Nektar::SolverUtils::ForcingBody::v_InitObject(), Nektar::ForcingMovingBody::v_InitObject(), and Nektar::SolverUtils::ForcingAbsorption::v_InitObject().

◆ m_sessionFunctions

std::map<std::string, SolverUtils::SessionFunctionSharedPtr> Nektar::SolverUtils::Forcing::m_sessionFunctions

protected

Map of known SessionFunctions.

Definition at line 111 of file Forcing.h.

Referenced by GetFunction().

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ~Forcing()

◆ Forcing()

Member Function Documentation

◆ Apply()

◆ EvaluateTimeFunction() [1/2]

◆ EvaluateTimeFunction() [2/2]

◆ GetForces()

◆ GetFunction()

◆ InitObject()

◆ Load()

◆ UpdateForces()

◆ v_Apply()

◆ v_InitObject()

Member Data Documentation

◆ m_equ

◆ m_Forcing

◆ m_NumVariable

◆ m_session

◆ m_sessionFunctions