carreauDynamics.h

00001 /* This file is part of the Palabos library.
00002  *
00003  * Copyright (C) 2011 FlowKit Sarl
00004  * Avenue de Chailly 23
00005  * 1012 Lausanne, Switzerland
00006  * E-mail contact: contact@flowkit.com
00007  *
00008  * The most recent release of Palabos can be downloaded at 
00009  * <http://www.palabos.org/>
00010  *
00011  * The library Palabos is free software: you can redistribute it and/or
00012  * modify it under the terms of the GNU Affero General Public License as
00013  * published by the Free Software Foundation, either version 3 of the
00014  * License, or (at your option) any later version.
00015  *
00016  * The library is distributed in the hope that it will be useful,
00017  * but WITHOUT ANY WARRANTY; without even the implied warranty of
00018  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
00019  * GNU Affero General Public License for more details.
00020  *
00021  * You should have received a copy of the GNU Affero General Public License
00022  * along with this program.  If not, see <http://www.gnu.org/licenses/>.
00023 */
00024 
00025 /* Main author: Orestis Malaspinas
00026  */
00027 
00028 #ifndef CARREAU_DYNAMICS_H
00029 #define CARREAU_DYNAMICS_H
00030 
00031 #include "core/globalDefs.h"
00032 #include "core/dynamics.h"
00033 #include "complexDynamics/variableOmegaDynamics.h"
00034 
00035 namespace plb {
00036 
00037         
00039 
00043 template<typename T, template<typename U> class Descriptor, int N>
00044 class CarreauDynamics : public OmegaFromPiDynamics<T,Descriptor> {
00045 public:
00046     CarreauDynamics(Dynamics<T,Descriptor>* baseDynamics_, bool automaticPrepareCollision=true);
00047     virtual T getOmegaFromPiAndRhoBar(Array<T,SymmetricTensor<T,Descriptor>::n> const& PiNeq, T rhoBar) const;
00048     CarreauDynamics<T,Descriptor,N>* clone() const;
00049     virtual int getId() const;
00050 private:
00051     static int id;
00052 };
00053 
00055 
00058 template<typename T, template<typename U> class Descriptor, int N>
00059 class BGKCarreauDynamics : public IsoThermalBulkDynamics<T,Descriptor> {
00060 public:
00061 /* *************** Construction / Destruction ************************ */
00062     BGKCarreauDynamics(T omega=(T)1);
00063 
00065     virtual BGKCarreauDynamics<T,Descriptor,N>* clone() const;
00066 
00068     virtual int getId() const;
00069 
00070 /* *************** Collision and Equilibrium ************************* */
00071 
00073     virtual void collide(Cell<T,Descriptor>& cell,
00074                          BlockStatistics& statistics_);
00075 
00077     virtual T computeEquilibrium(plint iPop, T rhoBar, Array<T,Descriptor<T>::d> const& j,
00078                                  T jSqr, T thetaBar=T()) const;
00079 private:
00080     static int id;
00081 };
00082 
00084 
00087 template<typename T, template<typename U> class Descriptor, int N>
00088 class RegularizedBGKCarreauDynamics : public IsoThermalBulkDynamics<T,Descriptor> {
00089 public:
00090 /* *************** Construction / Destruction ************************ */
00091     RegularizedBGKCarreauDynamics(T omega=(T)1);
00092 
00094     virtual RegularizedBGKCarreauDynamics<T,Descriptor,N>* clone() const;
00095 
00097     virtual int getId() const;
00098 
00099 /* *************** Collision and Equilibrium ************************* */
00100 
00102     virtual void collide(Cell<T,Descriptor>& cell,
00103                          BlockStatistics& statistics_);
00104 
00106     virtual T computeEquilibrium(plint iPop, T rhoBar, Array<T,Descriptor<T>::d> const& j,
00107                                  T jSqr, T thetaBar=T()) const;
00108 private:
00109     static int id;
00110 };
00111 
00112 } // namespace plb
00113 
00114 #endif  // VARIABLE_OMEGA_DYNAMICS_H