From 94d3e79a8617f88dc0219cfdeedfa3147833719d Mon Sep 17 00:00:00 2001 From: Adrian Kummerlaender Date: Mon, 24 Jun 2019 14:43:36 +0200 Subject: Initialize at openlb-1-3 --- ...okesAdvectionDiffusionCouplingPostProcessor3D.h | 244 +++++++++++++++++++++ 1 file changed, 244 insertions(+) create mode 100644 src/dynamics/navierStokesAdvectionDiffusionCouplingPostProcessor3D.h (limited to 'src/dynamics/navierStokesAdvectionDiffusionCouplingPostProcessor3D.h') diff --git a/src/dynamics/navierStokesAdvectionDiffusionCouplingPostProcessor3D.h b/src/dynamics/navierStokesAdvectionDiffusionCouplingPostProcessor3D.h new file mode 100644 index 0000000..e69a87a --- /dev/null +++ b/src/dynamics/navierStokesAdvectionDiffusionCouplingPostProcessor3D.h @@ -0,0 +1,244 @@ +/* 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_INTO_ADVECTION_DIFFUSION_COUPLING_POST_PROCESSOR_3D_H +#define NAVIER_STOKES_INTO_ADVECTION_DIFFUSION_COUPLING_POST_PROCESSOR_3D_H + +#include + +#include "core/spatiallyExtendedObject3D.h" +#include "core/postProcessing.h" +#include "core/blockLattice3D.h" +#include "advectionDiffusionForces.hh" +#include "advectionDiffusionForces.h" + +namespace olb { + + + +/** +* Multiphysics class for coupling between different lattices. +*/ +//====================================================================== +// ========Regularized NSDiffusion Coupling 3D ====================// +//====================================================================== +template +class NavierStokesAdvectionDiffusionCouplingPostProcessor3D : + public LocalPostProcessor3D { +public: + NavierStokesAdvectionDiffusionCouplingPostProcessor3D( + int x0_, int x1_, int y0_, int y1_, int z0_, int z1_, + T gravity_, T T0_, T deltaTemp_, std::vector dir_, + std::vector partners_); + int extent() const override + { + return 1; + } + int extent(int whichDirection) const override + { + return 1; + } + void process(BlockLattice3D& blockLattice) override; + void processSubDomain(BlockLattice3D& blockLattice, + int x0_, int x1_, int y0_, int y1_, int z0_, int z1_) override; +private: + typedef DESCRIPTOR L; + int x0, x1, y0, y1, z0, z1; + T gravity, T0, deltaTemp; + std::vector dir; + T forcePrefactor[L::d]; + + std::vector partners; + BlockLattice3D> *tPartner; +}; + +template +class NavierStokesAdvectionDiffusionCouplingGenerator3D : + public LatticeCouplingGenerator3D { +public: + NavierStokesAdvectionDiffusionCouplingGenerator3D( + int x0_, int x1_, int y0_, int y1_, int z0_, int z1_, + T gravity_, T T0_, T deltaTemp_, std::vector dir_); + PostProcessor3D* generate( + std::vector partners) const override; + LatticeCouplingGenerator3D* clone() const override; + +private: + T gravity, T0, deltaTemp; + std::vector dir; +}; + + +//====================================================================== +// ============= SmagorinskyBoussinesqCoupling 3D ===================// +//====================================================================== +template +class SmagorinskyBoussinesqCouplingPostProcessor3D : public LocalPostProcessor3D { +public: + SmagorinskyBoussinesqCouplingPostProcessor3D(int x0_, int x1_, int y0_, int y1_, int z0_, int z1_, + 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(BlockLattice3D& blockLattice) override; + void processSubDomain(BlockLattice3D& blockLattice, + int x0_, int x1_, int y0_, int y1_, int z0_, int z1_) override; +private: + typedef DESCRIPTOR L; + int x0, x1, y0, y1, z0, z1; + T gravity, T0, deltaTemp; + std::vector dir; + T PrTurb; + BlockLattice3D> *tPartner; + T forcePrefactor[L::d]; + T tauTurbADPrefactor; + + std::vector partners; +}; + +template +class SmagorinskyBoussinesqCouplingGenerator3D : public LatticeCouplingGenerator3D { +public: + SmagorinskyBoussinesqCouplingGenerator3D( + int x0_, int x1_, int y0_, int y1_, int z0_, int z1_, + T gravity_, T T0_, T deltaTemp_, std::vector dir_, T PrTurb_); + PostProcessor3D* generate(std::vector partners) const override; + LatticeCouplingGenerator3D* clone() const override; + +private: + T gravity, T0, deltaTemp; + std::vector dir; + T PrTurb; +}; + +//================================================================================================== +// ========Coupling 3D of Navier-Stokes on Advection-Diffusion====================// +//================================================================================================== +template> +class AdvectionDiffusionParticleCouplingPostProcessor3D : + public LocalPostProcessor3D { +public: + AdvectionDiffusionParticleCouplingPostProcessor3D( + int x0_, int x1_, int y0_, int y1_, int z0_, int z1_, int iC_, int offset_, + std::vector partners_, + std::vector > > forces_); + int extent() const override + { + return 1; + } + int extent(int whichDirection) const override + { + return 1; + } + void process(BlockLattice3D& blockLattice) override; + void processSubDomain(BlockLattice3D& blockLattice, + int x0_, int x1_, int y0_, int y1_, int z0_, int z1_) override; +private: + int x0, x1, y0, y1, z0, z1, iC; + T dragCoeff; + int offset; + T *vel, *vel_new, *velXp, *velXn, *velYp, *velYn, *velZp, *velZn; + bool par = true; + Cell *adCell; + Cell *nsCell; + +protected: + std::vector > > forces; +}; + +template> +class AdvectionDiffusionParticleCouplingGenerator3D : + public LatticeCouplingGenerator3D { +public: + AdvectionDiffusionParticleCouplingGenerator3D(int offset_); + PostProcessor3D* generate( + std::vector partners) const override; + LatticeCouplingGenerator3D* clone() const override; + void addForce(AdvectionDiffusionForce3D &force); + +private: + int offset; + +protected: + std::vector > > ADforces; +}; + + +//====================================================================== +// ======== Porous Regularized NSDiffusion Coupling 3D =================// +//====================================================================== +template +class PorousNavierStokesAdvectionDiffusionCouplingPostProcessor3D : + public LocalPostProcessor3D { +public: + PorousNavierStokesAdvectionDiffusionCouplingPostProcessor3D( + int x0_, int x1_, int y0_, int y1_, int z0_, int z1_, + T gravity_, T T0_, T deltaTemp_, std::vector dir_, + std::vector partners_); + int extent() const override + { + return 1; + } + int extent(int whichDirection) const override + { + return 1; + } + void process(BlockLattice3D& blockLattice) override; + void processSubDomain(BlockLattice3D& blockLattice, + int x0_, int x1_, int y0_, int y1_, int z0_, int z1_) override; +private: + int x0, x1, y0, y1, z0, z1; + T gravity, T0, deltaTemp; + std::vector dir; + + std::vector partners; +}; + +template +class PorousNavierStokesAdvectionDiffusionCouplingGenerator3D : + public LatticeCouplingGenerator3D { +public: + PorousNavierStokesAdvectionDiffusionCouplingGenerator3D( + int x0_, int x1_, int y0_, int y1_, int z0_, int z1_, + T gravity_, T T0_, T deltaTemp_, std::vector dir_); + PostProcessor3D* generate( + std::vector partners) const override; + LatticeCouplingGenerator3D* clone() const override; + +private: + T gravity, T0, deltaTemp; + std::vector dir; +}; + + +} + +#endif -- cgit v1.2.3