GRINS-0.8.0
averaged_fan_base.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 
26 #ifndef GRINS_AVERAGED_FAN_BASE_H
27 #define GRINS_AVERAGED_FAN_BASE_H
28 
29 // GRINS
30 #include "grins_config.h"
32 
33 // libMesh
34 #include "libmesh/getpot.h"
35 #include "libmesh/parsed_function.h"
36 #include "libmesh/tensor_value.h"
37 
38 // C++
39 #include <string>
40 
41 namespace GRINS
42 {
43 
44  template<class Viscosity>
46  {
47  public:
48 
49  AveragedFanBase( const std::string& physics_name, const GetPot& input );
50 
52 
53  bool compute_force ( const libMesh::Point& point,
54  const libMesh::Real time,
55  const libMesh::NumberVectorValue& U,
56  libMesh::NumberVectorValue& F,
57  libMesh::NumberTensorValue *dFdU = NULL);
58 
59  protected:
60 
61  // Velocity of the moving fan blades as a function of x,y,z
62  // Iff there are no fan blades moving past a location, the
63  // velocity there is specified as zero.
65 
66  // "Up" direction of fan airflow, as a function of x,y,z
67  // For most fans this will be a constant, the axis of rotation.
69 
70  // Coefficients of lift and drag as a function of angle "t" in
71  // radians. Should be well defined on [-pi, pi].
72  //
73  // No, "t" is not time in these functions.
74  // Yes, I'm abusing FunctionBase.
77 
78  // The chord length of the fan wing cross-section. For fan blades
79  // with constant cross-section this will be a constant.
81 
82  // The area swept out by the local fan wing cross-section. For
83  // cylindrical areas swept out by N fan blades, this is just
84  // pi*r^2*h/N
86 
87  // The angle-of-attack between the fan wing chord line and the fan
88  // velocity vector, in radians. For fan blades with no "twist"
89  // this will be a constant.
91 
92  private:
93 
95  void read_input_options( const GetPot& input );
96 
98  };
99 
100 } // end namespace block
101 
102 #endif // GRINS_AVERAGED_FAN_BASE_H
libMesh::ParsedFunction< libMesh::Number > drag_function
Physics class for Incompressible Navier-Stokes.
bool compute_force(const libMesh::Point &point, const libMesh::Real time, const libMesh::NumberVectorValue &U, libMesh::NumberVectorValue &F, libMesh::NumberTensorValue *dFdU=NULL)
libMesh::ParsedFunction< libMesh::Number > local_vertical_function
libMesh::ParsedFunction< libMesh::Number > chord_function
libMesh::ParsedFunction< libMesh::Number > base_velocity_function
libMesh::ParsedFunction< libMesh::Number > area_swept_function
GRINS namespace.
libMesh::ParsedFunction< libMesh::Number > lift_function
libMesh::ParsedFunction< libMesh::Number > aoa_function
void read_input_options(const GetPot &input)
Read options from GetPot input file.

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