Nektar++
Static Public Member Functions | Static Public Attributes | Protected Member Functions | Private Member Functions | Private Attributes | Friends | List of all members
Nektar::SolverUtils::ForcingBody Class Reference

#include <ForcingBody.h>

Inheritance diagram for Nektar::SolverUtils::ForcingBody:
[legend]

Static Public Member Functions

static SOLVER_UTILS_EXPORT ForcingSharedPtr create (const LibUtilities::SessionReaderSharedPtr &pSession, const std::weak_ptr< EquationSystem > &pEquation, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const unsigned int &pNumForcingFields, const TiXmlElement *pForce)
 Creates an instance of this class. More...
 
- Static Public Member Functions inherited from Nektar::SolverUtils::Forcing
static SOLVER_UTILS_EXPORT std::vector< ForcingSharedPtrLoad (const LibUtilities::SessionReaderSharedPtr &pSession, const std::weak_ptr< EquationSystem > &pEquation, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const unsigned int &pNumForcingFields=0)
 

Static Public Attributes

static std::string classNameBody
 Name of the class. More...
 
static std::string classNameField
 

Protected Member Functions

virtual SOLVER_UTILS_EXPORT void v_InitObject (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const unsigned int &pNumForcingFields, const TiXmlElement *pForce)
 
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)
 
virtual SOLVER_UTILS_EXPORT void v_ApplyCoeff (const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray, const NekDouble &time)
 
- Protected Member Functions inherited from Nektar::SolverUtils::Forcing
SOLVER_UTILS_EXPORT Forcing (const LibUtilities::SessionReaderSharedPtr &pSession, const std::weak_ptr< EquationSystem > &pEquation)
 Constructor. More...
 
virtual SOLVER_UTILS_EXPORT void v_PreApply (const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray, const NekDouble &time)
 
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)
 

Private Member Functions

 ForcingBody (const LibUtilities::SessionReaderSharedPtr &pSession, const std::weak_ptr< EquationSystem > &pEquation)
 
virtual ~ForcingBody (void)
 
void Update (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble >> &inarray, const NekDouble &time)
 

Private Attributes

std::string m_funcName
 
bool m_hasTimeFcnScaling
 
LibUtilities::EquationSharedPtr m_timeFcnEqn
 
bool m_transform
 

Friends

class MemoryManager< ForcingBody >
 

Additional Inherited Members

- Public Member Functions inherited from Nektar::SolverUtils::Forcing
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 void PreApply (const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray, const NekDouble &time)
 Change the advection velocity before applying the forcing. For example, subtracting the frame velocity from the advection velocity in the MovingRefercenceFrame. More...
 
SOLVER_UTILS_EXPORT void ApplyCoeff (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 ()
 
- Protected Attributes inherited from Nektar::SolverUtils::Forcing
LibUtilities::SessionReaderSharedPtr m_session
 Session reader. More...
 
const std::weak_ptr< EquationSystemm_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::SessionFunctionSharedPtrm_sessionFunctions
 Map of known SessionFunctions. More...
 

Detailed Description

Definition at line 50 of file ForcingBody.h.

Constructor & Destructor Documentation

◆ ForcingBody()

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

Definition at line 52 of file ForcingBody.cpp.

54  : Forcing(pSession, pEquation), m_hasTimeFcnScaling(false)
55 {
56 }
SOLVER_UTILS_EXPORT Forcing(const LibUtilities::SessionReaderSharedPtr &pSession, const std::weak_ptr< EquationSystem > &pEquation)
Constructor.
Definition: Forcing.cpp:50

◆ ~ForcingBody()

virtual Nektar::SolverUtils::ForcingBody::~ForcingBody ( void  )
inlineprivatevirtual

Definition at line 96 of file ForcingBody.h.

96 {};

Member Function Documentation

◆ create()

static SOLVER_UTILS_EXPORT ForcingSharedPtr Nektar::SolverUtils::ForcingBody::create ( const LibUtilities::SessionReaderSharedPtr pSession,
const std::weak_ptr< EquationSystem > &  pEquation,
const Array< OneD, MultiRegions::ExpListSharedPtr > &  pFields,
const unsigned int &  pNumForcingFields,
const TiXmlElement *  pForce 
)
inlinestatic

Creates an instance of this class.

Definition at line 56 of file ForcingBody.h.

61  {
64  p->InitObject(pFields, pNumForcingFields, pForce);
65  return p;
66  }
static std::shared_ptr< DataType > AllocateSharedPtr(const Args &...args)
Allocate a shared pointer from the memory pool.
SOLVER_UTILS_EXPORT typedef std::shared_ptr< Forcing > ForcingSharedPtr
A shared pointer to an EquationSystem object.
Definition: Forcing.h:52

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

◆ Update()

void Nektar::SolverUtils::ForcingBody::Update ( const Array< OneD, MultiRegions::ExpListSharedPtr > &  pFields,
const Array< OneD, Array< OneD, NekDouble >> &  inarray,
const NekDouble time 
)
private

Definition at line 123 of file ForcingBody.cpp.

126 {
127  for (int i = 0; i < m_NumVariable; ++i)
128  {
130  m_session->GetFunctionType(m_funcName, m_session->GetVariable(i)))
131  {
133  m_session->GetFunction(m_funcName, m_session->GetVariable(i));
134  if (!boost::iequals(eqn->GetVlist(), "x y z t"))
135  {
136  // Coupled forcing
137  int nq = pFields[0]->GetNpoints();
138  Array<OneD, NekDouble> xc(nq), yc(nq), zc(nq), t(nq, time);
139  std::string varstr = "x y z";
140  std::vector<Array<OneD, const NekDouble>> fielddata = {xc, yc,
141  zc, t};
142 
143  for (int j = 0; j < m_NumVariable; ++j)
144  {
145  varstr += " " + m_session->GetVariable(j);
146  fielddata.push_back(inarray[j]);
147  }
148 
149  // Evaluate function
150  m_session->GetFunction(m_funcName, m_session->GetVariable(i))
151  ->Evaluate(fielddata, m_Forcing[i]);
152  continue;
153  }
154  }
155  std::string s_FieldStr = m_session->GetVariable(i);
156  ASSERTL0(m_session->DefinesFunction(m_funcName, s_FieldStr),
157  "Variable '" + s_FieldStr + "' not defined.");
158  GetFunction(pFields, m_session, m_funcName, true)
159  ->Evaluate(s_FieldStr, m_Forcing[i], time);
160  }
161 
162  // If singleMode or halfMode, transform the forcing term to be in
163  // physical space in the plane, but Fourier space in the homogeneous
164  // direction
165  if (m_transform)
166  {
167  for (int i = 0; i < m_NumVariable; ++i)
168  {
169  pFields[0]->HomogeneousFwdTrans(m_Forcing[i], m_Forcing[i]);
170  }
171  }
172 }
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:215
int m_NumVariable
Number of variables.
Definition: Forcing.h:122
Array< OneD, Array< OneD, NekDouble > > m_Forcing
Evaluated forcing function.
Definition: Forcing.h:120
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.
Definition: Forcing.cpp:190
LibUtilities::SessionReaderSharedPtr m_session
Session reader.
Definition: Forcing.h:116
std::shared_ptr< Equation > EquationSharedPtr
Definition: Equation.h:130

References ASSERTL0, Nektar::LibUtilities::eFunctionTypeExpression, Nektar::SolverUtils::Forcing::GetFunction(), Nektar::SolverUtils::Forcing::m_Forcing, m_funcName, Nektar::SolverUtils::Forcing::m_NumVariable, Nektar::SolverUtils::Forcing::m_session, and m_transform.

Referenced by v_Apply(), v_ApplyCoeff(), and v_InitObject().

◆ v_Apply()

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

Implements Nektar::SolverUtils::Forcing.

Definition at line 174 of file ForcingBody.cpp.

178 {
180  {
181  Array<OneD, NekDouble> TimeFcn(1);
182 
183  for (int i = 0; i < m_NumVariable; i++)
184  {
185  EvaluateTimeFunction(time, m_timeFcnEqn, TimeFcn);
186 
187  Vmath::Svtvp(outarray[i].size(), TimeFcn[0], m_Forcing[i], 1,
188  outarray[i], 1, outarray[i], 1);
189  }
190  }
191  else
192  {
193  Update(fields, inarray, time);
194 
195  for (int i = 0; i < m_NumVariable; i++)
196  {
197  Vmath::Vadd(outarray[i].size(), outarray[i], 1, m_Forcing[i], 1,
198  outarray[i], 1);
199  }
200  }
201 }
LibUtilities::EquationSharedPtr m_timeFcnEqn
Definition: ForcingBody.h:90
void Update(const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble >> &inarray, const NekDouble &time)
SOLVER_UTILS_EXPORT void EvaluateTimeFunction(LibUtilities::SessionReaderSharedPtr pSession, std::string pFieldName, Array< OneD, NekDouble > &pArray, const std::string &pFunctionName, NekDouble pTime=NekDouble(0))
Definition: Forcing.cpp:166
void Svtvp(int n, const T alpha, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
svtvp (scalar times vector plus vector): z = alpha*x + y
Definition: Vmath.cpp:622
void Vadd(int n, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
Add vector z = x+y.
Definition: Vmath.cpp:359

References Nektar::SolverUtils::Forcing::EvaluateTimeFunction(), Nektar::SolverUtils::Forcing::m_Forcing, m_hasTimeFcnScaling, Nektar::SolverUtils::Forcing::m_NumVariable, m_timeFcnEqn, Vmath::Svtvp(), Update(), and Vmath::Vadd().

◆ v_ApplyCoeff()

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

Reimplemented from Nektar::SolverUtils::Forcing.

Definition at line 203 of file ForcingBody.cpp.

207 {
208  int ncoeff = outarray[m_NumVariable - 1].size();
209  Array<OneD, NekDouble> tmp(ncoeff, 0.0);
210 
212  {
213  Array<OneD, NekDouble> TimeFcn(1);
214 
215  for (int i = 0; i < m_NumVariable; ++i)
216  {
217  EvaluateTimeFunction(time, m_timeFcnEqn, TimeFcn);
218 
219  fields[i]->FwdTrans(m_Forcing[i], tmp);
220 
221  Vmath::Svtvp(ncoeff, TimeFcn[0], tmp, 1, outarray[i], 1,
222  outarray[i], 1);
223  }
224  }
225  else
226  {
227  Update(fields, inarray, time);
228 
229  for (int i = 0; i < m_NumVariable; ++i)
230  {
231  fields[i]->FwdTrans(m_Forcing[i], tmp);
232 
233  Vmath::Vadd(ncoeff, outarray[i], 1, tmp, 1, outarray[i], 1);
234  }
235  }
236 }

References Nektar::SolverUtils::Forcing::EvaluateTimeFunction(), Nektar::SolverUtils::Forcing::m_Forcing, m_hasTimeFcnScaling, Nektar::SolverUtils::Forcing::m_NumVariable, m_timeFcnEqn, Vmath::Svtvp(), Update(), and Vmath::Vadd().

◆ v_InitObject()

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

Implements Nektar::SolverUtils::Forcing.

Definition at line 58 of file ForcingBody.cpp.

61 {
62  m_NumVariable = pNumForcingFields;
63 
64  const TiXmlElement *funcNameElmt = pForce->FirstChildElement("BODYFORCE");
65  if (!funcNameElmt)
66  {
67  funcNameElmt = pForce->FirstChildElement("FIELDFORCE");
68 
69  ASSERTL0(funcNameElmt,
70  "Requires BODYFORCE or FIELDFORCE tag "
71  "specifying function name which prescribes body force.");
72  }
73 
74  m_funcName = funcNameElmt->GetText();
75  ASSERTL0(m_session->DefinesFunction(m_funcName),
76  "Function '" + m_funcName + "' not defined.");
77 
78  bool singleMode, halfMode;
79  m_session->MatchSolverInfo("ModeType", "SingleMode", singleMode, false);
80  m_session->MatchSolverInfo("ModeType", "HalfMode", halfMode, false);
81  bool homogeneous = pFields[0]->GetExpType() == MultiRegions::e3DH1D ||
82  pFields[0]->GetExpType() == MultiRegions::e3DH2D;
83  m_transform = (singleMode || halfMode || homogeneous);
84 
85  // Time function is optional
86  funcNameElmt = pForce->FirstChildElement("BODYFORCETIMEFCN");
87  if (!funcNameElmt)
88  {
89  funcNameElmt = pForce->FirstChildElement("FIELDFORCETIMEFCN");
90  }
91 
92  // Load time function if specified
93  if (funcNameElmt)
94  {
95  std::string funcNameTime = funcNameElmt->GetText();
96 
97  ASSERTL0(!funcNameTime.empty(),
98  "Expression must be given in BODYFORCETIMEFCN or "
99  "FIELDFORCETIMEFCN.");
100 
101  m_session->SubstituteExpressions(funcNameTime);
103  m_session->GetInterpreter(), funcNameTime);
104 
105  m_hasTimeFcnScaling = true;
106  }
107 
108  m_Forcing = Array<OneD, Array<OneD, NekDouble>>(m_NumVariable);
109  for (int i = 0; i < m_NumVariable; ++i)
110  {
111  m_Forcing[i] = Array<OneD, NekDouble>(pFields[0]->GetTotPoints(), 0.0);
112  }
113 
114  Array<OneD, Array<OneD, NekDouble>> tmp(pFields.size());
115  for (int i = 0; i < m_NumVariable; ++i)
116  {
117  tmp[i] = pFields[i]->GetPhys();
118  }
119 
120  Update(pFields, tmp, 0.0);
121 }

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), ASSERTL0, Nektar::MultiRegions::e3DH1D, Nektar::MultiRegions::e3DH2D, Nektar::SolverUtils::Forcing::m_Forcing, m_funcName, m_hasTimeFcnScaling, Nektar::SolverUtils::Forcing::m_NumVariable, Nektar::SolverUtils::Forcing::m_session, m_timeFcnEqn, m_transform, and Update().

Friends And Related Function Documentation

◆ MemoryManager< ForcingBody >

friend class MemoryManager< ForcingBody >
friend

Definition at line 1 of file ForcingBody.h.

Member Data Documentation

◆ classNameBody

std::string Nektar::SolverUtils::ForcingBody::classNameBody
static
Initial value:
=
"Body Forcing")
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
Definition: NekFactory.hpp:198
static SOLVER_UTILS_EXPORT ForcingSharedPtr create(const LibUtilities::SessionReaderSharedPtr &pSession, const std::weak_ptr< EquationSystem > &pEquation, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const unsigned int &pNumForcingFields, const TiXmlElement *pForce)
Creates an instance of this class.
Definition: ForcingBody.h:56
ForcingFactory & GetForcingFactory()
Declaration of the forcing factory singleton.
Definition: Forcing.cpp:44

Name of the class.

Definition at line 69 of file ForcingBody.h.

◆ classNameField

std::string Nektar::SolverUtils::ForcingBody::classNameField
static
Initial value:

Definition at line 70 of file ForcingBody.h.

◆ m_funcName

std::string Nektar::SolverUtils::ForcingBody::m_funcName
private

Definition at line 88 of file ForcingBody.h.

Referenced by Update(), and v_InitObject().

◆ m_hasTimeFcnScaling

bool Nektar::SolverUtils::ForcingBody::m_hasTimeFcnScaling
private

Definition at line 89 of file ForcingBody.h.

Referenced by v_Apply(), v_ApplyCoeff(), and v_InitObject().

◆ m_timeFcnEqn

LibUtilities::EquationSharedPtr Nektar::SolverUtils::ForcingBody::m_timeFcnEqn
private

Definition at line 90 of file ForcingBody.h.

Referenced by v_Apply(), v_ApplyCoeff(), and v_InitObject().

◆ m_transform

bool Nektar::SolverUtils::ForcingBody::m_transform
private

Definition at line 91 of file ForcingBody.h.

Referenced by Update(), and v_InitObject().