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

Wall boundary conditions for Rotational compressible flow problems. More...

#include <WallRotationalBC.h>

Inheritance diagram for Nektar::WallRotationalBC:
[legend]

Static Public Member Functions

static CFSBndCondSharedPtr create (const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble > > &pTraceNormals, const Array< OneD, Array< OneD, NekDouble > > &pGridVelocity, const int pSpaceDim, const int bcRegion, const int cnt)
 Creates an instance of this class. More...
 

Static Public Attributes

static std::string classNameRotational
 Name of the class. More...
 

Protected Member Functions

void v_Apply (Array< OneD, Array< OneD, NekDouble > > &Fwd, Array< OneD, Array< OneD, NekDouble > > &physarray, const NekDouble &time) override
 
- Protected Member Functions inherited from Nektar::CFSBndCond
 CFSBndCond (const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble > > &pTraceNormals, const Array< OneD, Array< OneD, NekDouble > > &pGridVelocity, const int pSpaceDim, const int bcRegion, const int cnt)
 Constructor. More...
 
virtual ~CFSBndCond ()=default
 
virtual void v_Apply (Array< OneD, Array< OneD, NekDouble > > &Fwd, Array< OneD, Array< OneD, NekDouble > > &physarray, const NekDouble &time)=0
 
virtual void v_ApplyBwdWeight ()
 

Private Member Functions

 WallRotationalBC (const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble > > &pTraceNormals, const Array< OneD, Array< OneD, NekDouble > > &pGridVelocity, const int pSpaceDim, const int bcRegion, const int cnt)
 
 ~WallRotationalBC (void) override=default
 

Friends

class MemoryManager< WallRotationalBC >
 

Additional Inherited Members

- Public Member Functions inherited from Nektar::CFSBndCond
void Apply (Array< OneD, Array< OneD, NekDouble > > &Fwd, Array< OneD, Array< OneD, NekDouble > > &physarray, const NekDouble &time=0)
 Apply the boundary condition. More...
 
void ApplyBwdWeight ()
 Apply the Weight of boundary condition. More...
 
- Protected Attributes inherited from Nektar::CFSBndCond
LibUtilities::SessionReaderSharedPtr m_session
 Session reader. More...
 
Array< OneD, MultiRegions::ExpListSharedPtrm_fields
 Array of fields. More...
 
Array< OneD, Array< OneD, NekDouble > > m_traceNormals
 Trace normals. More...
 
Array< OneD, Array< OneD, NekDouble > > m_gridVelocityTrace
 Grid Velocity. More...
 
int m_spacedim
 Space dimension. More...
 
VariableConverterSharedPtr m_varConv
 Auxiliary object to convert variables. More...
 
NekDouble m_diffusionAveWeight
 Weight for average calculation of diffusion term. More...
 
NekDouble m_gamma
 Parameters of the flow. More...
 
NekDouble m_rhoInf
 
NekDouble m_pInf
 
NekDouble m_pOut
 
Array< OneD, NekDoublem_velInf
 
NekDouble m_angVel
 
int m_bcRegion
 Id of the boundary region. More...
 
int m_offset
 Offset. More...
 

Detailed Description

Wall boundary conditions for Rotational compressible flow problems.

Definition at line 46 of file WallRotationalBC.h.

Constructor & Destructor Documentation

◆ WallRotationalBC()

Nektar::WallRotationalBC::WallRotationalBC ( const LibUtilities::SessionReaderSharedPtr pSession,
const Array< OneD, MultiRegions::ExpListSharedPtr > &  pFields,
const Array< OneD, Array< OneD, NekDouble > > &  pTraceNormals,
const Array< OneD, Array< OneD, NekDouble > > &  pGridVelocity,
const int  pSpaceDim,
const int  bcRegion,
const int  cnt 
)
private

Definition at line 45 of file WallRotationalBC.cpp.

51 : CFSBndCond(pSession, pFields, pTraceNormals, pGridVelocity, pSpaceDim,
52 bcRegion, cnt)
53{
55
56 // Set up rotational boundary edge velocities
57 const Array<OneD, const int> &traceBndMap = m_fields[0]->GetTraceBndMap();
58
59 int eMax = m_fields[0]->GetBndCondExpansions()[m_bcRegion]->GetExpSize();
60 for (int e = 0; e < eMax; ++e)
61 {
62 int nBCEdgePts = m_fields[0]
63 ->GetBndCondExpansions()[m_bcRegion]
64 ->GetExp(e)
65 ->GetTotPoints();
66 int id2 =
67 m_fields[0]->GetTrace()->GetPhys_Offset(traceBndMap[m_offset + e]);
68
69 Array<OneD, NekDouble> x(nBCEdgePts, 0.0);
70 Array<OneD, NekDouble> y(nBCEdgePts, 0.0);
71 m_fields[0]->GetBndCondExpansions()[m_bcRegion]->GetExp(e)->GetCoords(
72 x, y);
73 for (int j = 0; j < nBCEdgePts; ++j)
74 {
75 m_gridVelocityTrace[0][id2 + j] = m_angVel * -y[j];
76 m_gridVelocityTrace[1][id2 + j] = m_angVel * x[j];
77 }
78 }
79}
NekDouble m_diffusionAveWeight
Weight for average calculation of diffusion term.
Definition: CFSBndCond.h:102
Array< OneD, Array< OneD, NekDouble > > m_gridVelocityTrace
Grid Velocity.
Definition: CFSBndCond.h:96
int m_bcRegion
Id of the boundary region.
Definition: CFSBndCond.h:113
CFSBndCond(const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble > > &pTraceNormals, const Array< OneD, Array< OneD, NekDouble > > &pGridVelocity, const int pSpaceDim, const int bcRegion, const int cnt)
Constructor.
Definition: CFSBndCond.cpp:46
int m_offset
Offset.
Definition: CFSBndCond.h:115
Array< OneD, MultiRegions::ExpListSharedPtr > m_fields
Array of fields.
Definition: CFSBndCond.h:92
NekDouble m_angVel
Definition: CFSBndCond.h:110

References Nektar::CFSBndCond::m_angVel, Nektar::CFSBndCond::m_bcRegion, Nektar::CFSBndCond::m_diffusionAveWeight, Nektar::CFSBndCond::m_fields, Nektar::CFSBndCond::m_gridVelocityTrace, and Nektar::CFSBndCond::m_offset.

◆ ~WallRotationalBC()

Nektar::WallRotationalBC::~WallRotationalBC ( void  )
overrideprivatedefault

Member Function Documentation

◆ create()

static CFSBndCondSharedPtr Nektar::WallRotationalBC::create ( const LibUtilities::SessionReaderSharedPtr pSession,
const Array< OneD, MultiRegions::ExpListSharedPtr > &  pFields,
const Array< OneD, Array< OneD, NekDouble > > &  pTraceNormals,
const Array< OneD, Array< OneD, NekDouble > > &  pGridVelocity,
const int  pSpaceDim,
const int  bcRegion,
const int  cnt 
)
inlinestatic

Creates an instance of this class.

Definition at line 52 of file WallRotationalBC.h.

58 {
61 pSession, pFields, pTraceNormals, pGridVelocity, pSpaceDim,
62 bcRegion, cnt);
63 return p;
64 }
static std::shared_ptr< DataType > AllocateSharedPtr(const Args &...args)
Allocate a shared pointer from the memory pool.
std::shared_ptr< CFSBndCond > CFSBndCondSharedPtr
A shared pointer to a boundary condition object.
Definition: CFSBndCond.h:52

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

◆ v_Apply()

void Nektar::WallRotationalBC::v_Apply ( Array< OneD, Array< OneD, NekDouble > > &  Fwd,
Array< OneD, Array< OneD, NekDouble > > &  physarray,
const NekDouble time 
)
overrideprotectedvirtual

Implements Nektar::CFSBndCond.

Definition at line 81 of file WallRotationalBC.cpp.

84{
85 int nVariables = physarray.size();
86
87 // @TODO: ALE we subtract the grid velocity ? "Set u = to ug for this one" -
88 // Dave
89
90 int i;
91 const Array<OneD, const int> &traceBndMap = m_fields[0]->GetTraceBndMap();
92
93 // Take into account that for PDE based shock capturing, eps = 0 at the
94 // wall. Adjust the physical values of the trace to take user defined
95 // boundaries into account
96 int e, id1, id2, nBCEdgePts, eMax;
97
98 eMax = m_fields[0]->GetBndCondExpansions()[m_bcRegion]->GetExpSize();
99
100 for (e = 0; e < eMax; ++e)
101 {
102 nBCEdgePts = m_fields[0]
103 ->GetBndCondExpansions()[m_bcRegion]
104 ->GetExp(e)
105 ->GetTotPoints();
106 id1 =
107 m_fields[0]->GetBndCondExpansions()[m_bcRegion]->GetPhys_Offset(e);
108 id2 =
109 m_fields[0]->GetTrace()->GetPhys_Offset(traceBndMap[m_offset + e]);
110
111 // Boundary condition for epsilon term. @TODO: Is this correct, or
112 // should I do E = p/(gamma -1) + 1/2*rho(u^2 +v^2 + w^2)... or
113 if (nVariables == m_spacedim + 3)
114 {
115 Vmath::Zero(nBCEdgePts, &Fwd[nVariables - 1][id2], 1);
116 }
117
118 for (i = 0; i < m_spacedim; i++)
119 {
120 // V = -Vin
121 // Vmath::Neg(nBCEdgePts, &Fwd[i+1][id2], 1);
122
123 // This now does Vg * rho + Vin
124 // Vmath::Vvtvp(nBCEdgePts, &m_gridVelocityTrace[i][id2], 1,
125 // &Fwd[0][id2], 1, &Fwd[i+1][id2], 1, &Fwd[i+1][id2], 1);
126
127 for (int j = 0; j < nBCEdgePts; ++j)
128 {
129 // Fwd[i+1][id2 + j] = 2 * m_gridVelocityTrace[i][id2 + j] *
130 // Fwd[0][id2 + j] - Fwd[i+1][id2 + j];
131 Fwd[i + 1][id2 + j] =
132 2 * m_gridVelocityTrace[i][id2 + j] * Fwd[0][id2 + j] -
133 Fwd[i + 1][id2 + j];
134 }
135 }
136
137 // Copy boundary adjusted values into the boundary expansion
138 for (i = 0; i < nVariables; ++i)
139 {
140 Vmath::Vcopy(nBCEdgePts, &Fwd[i][id2], 1,
141 &(m_fields[i]
142 ->GetBndCondExpansions()[m_bcRegion]
143 ->UpdatePhys())[id1],
144 1);
145 }
146 }
147}
int m_spacedim
Space dimension.
Definition: CFSBndCond.h:98
void Zero(int n, T *x, const int incx)
Zero vector.
Definition: Vmath.hpp:273
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.hpp:825

References Nektar::CFSBndCond::m_bcRegion, Nektar::CFSBndCond::m_fields, Nektar::CFSBndCond::m_gridVelocityTrace, Nektar::CFSBndCond::m_offset, Nektar::CFSBndCond::m_spacedim, Vmath::Vcopy(), and Vmath::Zero().

Friends And Related Function Documentation

◆ MemoryManager< WallRotationalBC >

friend class MemoryManager< WallRotationalBC >
friend

Definition at line 1 of file WallRotationalBC.h.

Member Data Documentation

◆ classNameRotational

std::string Nektar::WallRotationalBC::classNameRotational
static
Initial value:
=
"WallRotational", WallRotationalBC::create,
"Adiabatic rotational wall boundary condition.")
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
static CFSBndCondSharedPtr create(const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble > > &pTraceNormals, const Array< OneD, Array< OneD, NekDouble > > &pGridVelocity, const int pSpaceDim, const int bcRegion, const int cnt)
Creates an instance of this class.
CFSBndCondFactory & GetCFSBndCondFactory()
Declaration of the boundary condition factory singleton.
Definition: CFSBndCond.cpp:40

Name of the class.

Definition at line 67 of file WallRotationalBC.h.