/* This file is part of the OpenLB library * * Copyright (C) 2008 Orestis Malaspinas, Andrea Parmigiani, Jonas Latt * E-mail contact: info@openlb.net * The most recent release of OpenLB can be downloaded at * * * This program is free software; you can redistribute it and/or * modify it under the terms of the GNU General Public License * as published by the Free Software Foundation; either version 2 * of the License, or (at your option) any later version. * * This program is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * GNU General Public License for more details. * * You should have received a copy of the GNU General Public * License along with this program; if not, write to the Free * Software Foundation, Inc., 51 Franklin Street, Fifth Floor, * Boston, MA 02110-1301, USA. */ #ifndef NAVIER_STOKES_ADVECTION_DIFFUSION_COUPLING_POST_PROCESSOR_2D_H #define NAVIER_STOKES_ADVECTION_DIFFUSION_COUPLING_POST_PROCESSOR_2D_H #include "core/spatiallyExtendedObject2D.h" #include "core/postProcessing.h" #include "core/blockLattice2D.h" #include namespace olb { /** * Class for the coupling between a Navier-Stokes (NS) lattice and an * Advection-Diffusion (AD) lattice. */ // !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!// // This coupling must be necessarily be put on the Navier-Stokes lattice!! // !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!// //====================================================================== // ========Regularized NSDiffusion Coupling 2D ====================// //====================================================================== template class NavierStokesAdvectionDiffusionCouplingPostProcessor2D : public LocalPostProcessor2D { public: NavierStokesAdvectionDiffusionCouplingPostProcessor2D(int x0_, int x1_, int y0_, int y1_, T gravity_, T T0_, T deltaTemp_, std::vector dir_, std::vector partners_); int extent() const override { return 0; } int extent(int whichDirection) const override { return 0; } void process(BlockLattice2D& blockLattice) override; void processSubDomain(BlockLattice2D& blockLattice, int x0_, int x1_, int y0_, int y1_) override; private: typedef DESCRIPTOR L; int x0, x1, y0, y1; T gravity, T0, deltaTemp; std::vector dir; BlockLattice2D> *tPartner; T forcePrefactor[L::d]; std::vector partners; }; template class NavierStokesAdvectionDiffusionCouplingGenerator2D : public LatticeCouplingGenerator2D { public: NavierStokesAdvectionDiffusionCouplingGenerator2D(int x0_, int x1_, int y0_, int y1_, T gravity_, T T0_, T deltaTemp_, std::vector dir_); PostProcessor2D* generate(std::vector partners) const override; LatticeCouplingGenerator2D* clone() const override; private: T gravity, T0, deltaTemp; std::vector dir; }; template class SmagorinskyBoussinesqCouplingPostProcessor2D : public LocalPostProcessor2D { public: SmagorinskyBoussinesqCouplingPostProcessor2D(int x0_, int x1_, int y0_, int y1_, T gravity_, T T0_, T deltaTemp_, std::vector dir_, T PrTurb_, std::vector partners_); int extent() const override { return 0; } int extent(int whichDirection) const override { return 0; } void process(BlockLattice2D& blockLattice) override; void processSubDomain(BlockLattice2D& blockLattice, int x0_, int x1_, int y0_, int y1_) override; private: typedef DESCRIPTOR L; int x0, x1, y0, y1; T gravity, T0, deltaTemp; std::vector dir; T PrTurb; BlockLattice2D> *tPartner; T forcePrefactor[L::d]; T tauTurbADPrefactor; std::vector partners; enum { velOffset = descriptors::D2Q5::template index(), forceOffset = descriptors::D2Q9::template index(), tauADoffset = descriptors::D2Q5::template index(), tauNSoffset = descriptors::D2Q9::template index() }; }; template class SmagorinskyBoussinesqCouplingGenerator2D : public LatticeCouplingGenerator2D { public: SmagorinskyBoussinesqCouplingGenerator2D(int x0_, int x1_, int y0_, int y1_, T gravity_, T T0_, T deltaTemp_, std::vector dir_, T PrTurb_); PostProcessor2D* generate(std::vector partners) const override; LatticeCouplingGenerator2D* clone() const override; private: T gravity, T0, deltaTemp; std::vector dir; T PrTurb; }; template class MixedScaleBoussinesqCouplingPostProcessor2D : public LocalPostProcessor2D { public: MixedScaleBoussinesqCouplingPostProcessor2D(int x0_, int x1_, int y0_, int y1_, T gravity_, T T0_, T deltaTemp_, std::vector dir_, T PrTurb_, std::vector partners_); int extent() const override { return 0; } int extent(int whichDirection) const override { return 0; } void process(BlockLattice2D& blockLattice) override; void processSubDomain(BlockLattice2D& blockLattice, int x0_, int x1_, int y0_, int y1_) override; private: typedef DESCRIPTOR L; int x0, x1, y0, y1; T gravity, T0, deltaTemp, PrTurb; std::vector dir; BlockLattice2D> *tPartner; T forcePrefactor[L::d]; T tauTurbADPrefactor; std::vector partners; }; template class MixedScaleBoussinesqCouplingGenerator2D : public LatticeCouplingGenerator2D { public: MixedScaleBoussinesqCouplingGenerator2D(int x0_, int x1_, int y0_, int y1_, T gravity_, T T0_, T deltaTemp_, std::vector dir_, T PrTurb_); PostProcessor2D* generate(std::vector partners) const override; LatticeCouplingGenerator2D* clone() const override; private: T gravity, T0, deltaTemp, PrTurb; std::vector dir; }; } #endif