Nektar++
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
WallViscousBC.cpp
Go to the documentation of this file.
1///////////////////////////////////////////////////////////////////////////////
2//
3// File: WallViscousBC.cpp
4//
5// For more information, please see: http://www.nektar.info
6//
7// The MIT License
8//
9// Copyright (c) 2006 Division of Applied Mathematics, Brown University (USA),
10// Department of Aeronautics, Imperial College London (UK), and Scientific
11// Computing and Imaging Institute, University of Utah (USA).
12//
13// Permission is hereby granted, free of charge, to any person obtaining a
14// copy of this software and associated documentation files (the "Software"),
15// to deal in the Software without restriction, including without limitation
16// the rights to use, copy, modify, merge, publish, distribute, sublicense,
17// and/or sell copies of the Software, and to permit persons to whom the
18// Software is furnished to do so, subject to the following conditions:
19//
20// The above copyright notice and this permission notice shall be included
21// in all copies or substantial portions of the Software.
22//
23// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
24// OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
25// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
26// THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
27// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
28// FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
29// DEALINGS IN THE SOFTWARE.
30//
31// Description: No-slip wall boundary condition
32//
33///////////////////////////////////////////////////////////////////////////////
34
35#include <boost/core/ignore_unused.hpp>
36
37#include "WallViscousBC.h"
38
39using namespace std;
40
41namespace Nektar
42{
43
46 "WallViscous", WallViscousBC::create,
47 "No-slip (viscous) wall boundary condition.");
48
51 "WallAdiabatic", WallViscousBC::create,
52 "Adiabatic wall boundary condition.");
53
57 const Array<OneD, Array<OneD, NekDouble>> &pTraceNormals,
58 const int pSpaceDim, const int bcRegion, const int cnt)
59 : CFSBndCond(pSession, pFields, pTraceNormals, pSpaceDim, bcRegion, cnt)
60{
62
64}
65
68 const NekDouble &time)
69{
70 boost::ignore_unused(time);
71
72 int i;
73 int nVariables = physarray.size();
74
75 // Find the fields whose WallViscous/Adiabatic-BC is time-dependent
76 // Update variables on the boundaries of these fields
77 // Get the updated variables on the WallViscous/Adiabatic boundary
78 //
79 // Maybe the EvaluateBoundaryConditions() should be put upstream to
80 // CompressibleFlowSystem::NumCalRiemFluxJac(), So that the BCs will not
81 // be repeatedly updated when there are more than one time-dependent BC.
82 std::string varName;
83 for (i = 0; i < nVariables; ++i)
84 {
85 if (m_fields[i]->GetBndConditions()[m_bcRegion]->IsTimeDependent())
86 {
87 varName = m_session->GetVariable(i);
88 m_fields[i]->EvaluateBoundaryConditions(time, varName);
89
90 m_bndPhys[i] =
91 m_fields[i]->GetBndCondExpansions()[m_bcRegion]->UpdatePhys();
92 }
93 }
94
95 const Array<OneD, const int> &traceBndMap = m_fields[0]->GetTraceBndMap();
96
97 // Take into account that for PDE based shock capturing, eps = 0 at the
98 // wall. Adjust the physical values of the trace to take user defined
99 // boundaries into account
100 int e, id1, id2, nBCEdgePts, eMax;
101
102 eMax = m_fields[0]->GetBndCondExpansions()[m_bcRegion]->GetExpSize();
103
104 for (e = 0; e < eMax; ++e)
105 {
106 nBCEdgePts = m_fields[0]
107 ->GetBndCondExpansions()[m_bcRegion]
108 ->GetExp(e)
109 ->GetTotPoints();
110 id1 =
111 m_fields[0]->GetBndCondExpansions()[m_bcRegion]->GetPhys_Offset(e);
112 id2 =
113 m_fields[0]->GetTrace()->GetPhys_Offset(traceBndMap[m_offset + e]);
114
115 // Boundary condition for epsilon term.
116 if (nVariables == m_spacedim + 3)
117 {
118 Vmath::Zero(nBCEdgePts, &Fwd[nVariables - 1][id2], 1);
119 }
120
121 // V = - Vin
122 for (i = 0; i < m_spacedim; i++)
123 {
124 Vmath::Neg(nBCEdgePts, &Fwd[i + 1][id2], 1);
125 }
126
127 // Superimpose the perturbation
128 for (i = 0; i < nVariables; ++i)
129 {
130 if (m_fields[i]->GetBndConditions()[m_bcRegion]->IsTimeDependent())
131 {
132 Vmath::Vadd(nBCEdgePts, &m_bndPhys[i][id1], 1, &Fwd[i][id2], 1,
133 &Fwd[i][id2], 1);
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}
148
149} // namespace Nektar
Encapsulates the user-defined boundary conditions for compressible flow solver.
Definition: CFSBndCond.h:70
NekDouble m_diffusionAveWeight
Weight for average calculation of diffusion term.
Definition: CFSBndCond.h:99
LibUtilities::SessionReaderSharedPtr m_session
Session reader.
Definition: CFSBndCond.h:89
int m_spacedim
Space dimension.
Definition: CFSBndCond.h:95
int m_bcRegion
Id of the boundary region.
Definition: CFSBndCond.h:109
int m_offset
Offset.
Definition: CFSBndCond.h:111
Array< OneD, MultiRegions::ExpListSharedPtr > m_fields
Array of fields.
Definition: CFSBndCond.h:91
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
Definition: NekFactory.hpp:198
static std::string classNameAdiabatic
Definition: WallViscousBC.h:65
WallViscousBC(const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble > > &pTraceNormals, const int pSpaceDim, const int bcRegion, const int cnt)
Array< OneD, Array< OneD, NekDouble > > m_bndPhys
Definition: WallViscousBC.h:70
virtual void v_Apply(Array< OneD, Array< OneD, NekDouble > > &Fwd, Array< OneD, Array< OneD, NekDouble > > &physarray, const NekDouble &time) override
static CFSBndCondSharedPtr create(const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble > > &pTraceNormals, const int pSpaceDim, const int bcRegion, const int cnt)
Creates an instance of this class.
Definition: WallViscousBC.h:52
static std::string classNameViscous
Name of the class.
Definition: WallViscousBC.h:64
std::shared_ptr< SessionReader > SessionReaderSharedPtr
The above copyright notice and this permission notice shall be included.
Definition: CoupledSolver.h:2
CFSBndCondFactory & GetCFSBndCondFactory()
Declaration of the boundary condition factory singleton.
Definition: CFSBndCond.cpp:41
double NekDouble
void Neg(int n, T *x, const int incx)
Negate x = -x.
Definition: Vmath.cpp:513
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:354
void Zero(int n, T *x, const int incx)
Zero vector.
Definition: Vmath.cpp:487
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1191