GRINS-0.8.0
low_mach_navier_stokes_stab_helper.h
Go to the documentation of this file.
1 //-----------------------------------------------------------------------bl-
2 //--------------------------------------------------------------------------
3 //
4 // GRINS - General Reacting Incompressible Navier-Stokes
5 //
6 // Copyright (C) 2014-2017 Paul T. Bauman, Roy H. Stogner
7 // Copyright (C) 2010-2013 The PECOS Development Team
8 //
9 // This library is free software; you can redistribute it and/or
10 // modify it under the terms of the Version 2.1 GNU Lesser General
11 // Public License as published by the Free Software Foundation.
12 //
13 // This library is distributed in the hope that it will be useful,
14 // but WITHOUT ANY WARRANTY; without even the implied warranty of
15 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
16 // Lesser General Public License for more details.
17 //
18 // You should have received a copy of the GNU Lesser General Public
19 // License along with this library; if not, write to the Free Software
20 // Foundation, Inc. 51 Franklin Street, Fifth Floor,
21 // Boston, MA 02110-1301 USA
22 //
23 //-----------------------------------------------------------------------el-
24 
25 #ifndef LOW_MACH_NAVIER_STOKES_STAB_HELPER_H
26 #define LOW_MACH_NAVIER_STOKES_STAB_HELPER_H
27 
28 // GRINS
30 
31 namespace GRINS
32 {
34  {
35  public:
36 
38  ( const std::string& helper_name,
39  const GetPot& input);
40 
42 
43  libMesh::Real compute_tau_energy( AssemblyContext& c,
44  unsigned int qp,
45  libMesh::RealGradient& g,
46  libMesh::RealTensor& G,
47  libMesh::Real rho,
48  libMesh::Gradient U,
49  libMesh::Real k,
50  libMesh::Real cp,
51  bool is_steady ) const;
52 
53  }; // class LowMachNavierStokesStabilizationHelper
54 
55  /* ------------- Inline Functions ---------------*/
56  inline
58  unsigned int qp,
59  libMesh::RealGradient& g,
60  libMesh::RealTensor& G,
61  libMesh::Real rho,
62  libMesh::Gradient U,
63  libMesh::Real k,
64  libMesh::Real cp,
65  bool is_steady ) const
66  {
67  /*
68  libMesh::Real tau = (rho*U)*(G*(rho*U)) + this->_C*k*k*G.contract(G);
69 
70  if(!is_steady)
71  tau += (2.0*rho/c.get_deltat_value())*(2.0*rho/c.get_deltat_value());
72 
73  return this->_tau_factor/std::sqrt(tau);
74  */
75  return this->compute_tau( c, qp, k*k, g, G, rho*cp, U, is_steady );
76  }
77 
78 }
79 #endif // LOW_MACH_NAVIER_STOKES_STAB_HELPER_H
GRINS namespace.
LowMachNavierStokesStabilizationHelper(const std::string &helper_name, const GetPot &input)
libMesh::Real compute_tau(AssemblyContext &c, unsigned int qp, libMesh::Real mat_prop_sq, libMesh::RealGradient &g, libMesh::RealTensor &G, libMesh::Real rho, libMesh::Gradient U, bool is_steady) const
libMesh::Real compute_tau_energy(AssemblyContext &c, unsigned int qp, libMesh::RealGradient &g, libMesh::RealTensor &G, libMesh::Real rho, libMesh::Gradient U, libMesh::Real k, libMesh::Real cp, bool is_steady) const

Generated on Tue Dec 19 2017 12:47:28 for GRINS-0.8.0 by  doxygen 1.8.9.1