NavierStokes.hpp

Go to the documentation of this file.

// Copyright (C) 2010-2011 von Karman Institute for Fluid Dynamics, Belgium
//
// This software is distributed under the terms of the
// GNU Lesser General Public License version 3 (LGPLv3).
// See doc/lgpl.txt and doc/gpl.txt for the license text.

#ifndef cf3_UFEM_NavierStokes_hpp
#define cf3_UFEM_NavierStokes_hpp

#include <boost/accumulators/accumulators.hpp>
#include <boost/accumulators/statistics/stats.hpp>
#include <boost/accumulators/statistics/max.hpp>

#define BOOST_PROTO_MAX_ARITY 11
#ifdef BOOST_MPL_LIMIT_METAFUNCTION_ARITY
  #undef BOOST_MPL_LIMIT_METAFUNCTION_ARITY
#endif
#define BOOST_MPL_LIMIT_METAFUNCTION_ARITY 11

#include <boost/scoped_ptr.hpp>

#include "LibUFEM.hpp"
#include "LSSActionUnsteady.hpp"
#include "NavierStokesPhysics.hpp"
#include "SUPG.hpp"

namespace cf3 {
  
  namespace solver { class ActionDirector; }
  
namespace UFEM {

class UFEM_API NavierStokes : public LSSActionUnsteady
{
public: // functions

  NavierStokes ( const std::string& name );
  
  virtual ~NavierStokes();

  static std::string type_name () { return "NavierStokes"; }

private:
  void trigger_assembly();

  virtual void on_initial_conditions_set(InitialConditions& initial_conditions);

  template<typename GenericElementsT, typename SpecializedElementsT>
  void set_assembly_expression(const std::string& action_name);
  
  void set_triag_assembly(const bool use_specialization);
  void set_quad_assembly();
  void set_tetra_assembly(const bool use_specialization);
  void set_hexa_assembly();
  void set_prism_assembly();
  

  FieldVariable<0, VectorField> u;
  FieldVariable<1, ScalarField> p;
  FieldVariable<2, VectorField> u_adv;
  FieldVariable<3, VectorField> u1;
  FieldVariable<4, VectorField> u2;
  FieldVariable<5, VectorField> u3;
  FieldVariable<6, ScalarField> nu_eff;

  FieldVariable<7, ScalarField> T;
  
  // Body force
  FieldVariable<8, VectorField> g;


  PhysicsConstant rho;
  PhysicsConstant nu;

  PhysicsConstant Tref;

  Real tau_ps, tau_su, tau_bulk;
  
  Handle<solver::ActionDirector> m_assembly;
  Handle<solver::ActionDirector> m_update;
  Handle<common::Action> m_initial_conditions;

  ComputeTau compute_tau;
};

} // UFEM
} // cf3


#endif // cf3_UFEM_NavierStokes_hpp