dynamicSmagorinskyDynamics.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 #ifndef DYNAMIC_SMAGORINSKY_DYNAMICS_H
00026 #define DYNAMIC_SMAGORINSKY_DYNAMICS_H
00027 
00028 #include "core/globalDefs.h"
00029 #include "core/dynamics.h"
00030 
00031 namespace plb {
00032 
00034 template<typename T, template<typename U> class Descriptor>
00035 class ExternalOmegaDynamics : public CompositeDynamics<T,Descriptor> {
00036 public:
00037     ExternalOmegaDynamics(Dynamics<T,Descriptor>* baseDynamics_)
00038         : CompositeDynamics<T,Descriptor>(baseDynamics_)
00039     { }
00040     virtual void prepareCollision(Cell<T,Descriptor>& cell) {
00041         // Copy relaxation parameter from external scalar.
00042         this->setOmega(*cell.getExternal(Descriptor<T>::ExternalField::omegaBeginsAt));
00043     }
00044     ExternalOmegaDynamics<T,Descriptor>* clone() const {
00045         return new ExternalOmegaDynamics(*this);
00046     }
00047 };
00048 
00049 } // namespace plb
00050 
00051 #endif  // DYNAMIC_SMAGORINSKY_DYNAMICS_H