Nektar++
Go to the documentation of this file.
1 ///////////////////////////////////////////////////////////////////////////////
2 //
3 // File UnsteadyInviscidBurger.h
4 //
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 // License for the specific language governing rights and limitations under
14 // Permission is hereby granted, free of charge, to any person obtaining a
15 // copy of this software and associated documentation files (the "Software"),
16 // to deal in the Software without restriction, including without limitation
17 // the rights to use, copy, modify, merge, publish, distribute, sublicense,
18 // and/or sell copies of the Software, and to permit persons to whom the
19 // Software is furnished to do so, subject to the following conditions:
20 //
21 // The above copyright notice and this permission notice shall be included
22 // in all copies or substantial portions of the Software.
23 //
24 // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
25 // OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
26 // FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
27 // THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
28 // LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
29 // FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
30 // DEALINGS IN THE SOFTWARE.
31 //
32 // Description: Unsteady inviscid Burger solve routines
33 //
34 ///////////////////////////////////////////////////////////////////////////////
35
38
42
43 using namespace Nektar::SolverUtils;
44
45 namespace Nektar
46 {
48  {
49  public:
51
52  /// Creates an instance of this class
53  static EquationSystemSharedPtr create(
54  const LibUtilities::SessionReaderSharedPtr& pSession) {
56  p->InitObject();
57  return p;
58  }
59  /// Name of class
60  static std::string className;
61
62  /// Destructor
64
65  protected:
68  Array<OneD, NekDouble> m_traceVn;
69
70  /// Session reader
72
73  /// Evaluate the flux at each solution point
74  void GetFluxVector(
75  const Array<OneD, Array<OneD, NekDouble> > &physfield,
76  Array<OneD, Array<OneD, Array<OneD, NekDouble> > > &flux);
77
78  /// Compute the RHS
79  void DoOdeRhs(
80  const Array<OneD, const Array<OneD, NekDouble> >&inarray,
81  Array<OneD, Array<OneD, NekDouble> >&outarray,
82  const NekDouble time);
83
84  /// Compute the projection
85  void DoOdeProjection(
86  const Array<OneD, const Array<OneD, NekDouble> > &inarray,
87  Array<OneD, Array<OneD, NekDouble> > &outarray,
88  const NekDouble time);
89
90  /// Get the normal velocity
91  Array<OneD, NekDouble> &GetNormalVelocity();
92
93  /// Initialise the object
94  virtual void v_InitObject();
95
96  private:
97  };
98 }
99
100 #endif