/* 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 SHAN_CHEN_FORCED_POST_PROCESSOR_2D_H
#define SHAN_CHEN_FORCED_POST_PROCESSOR_2D_H
#include "core/spatiallyExtendedObject2D.h"
#include "core/postProcessing.h"
#include "core/blockLattice2D.h"
namespace olb {
/**
* Multiphysics class for coupling between different lattices.
*/
// =========================================================================//
// ===========Shan Chen coupling without wall interaction===================//
// =========================================================================//
template
class ShanChenForcedPostProcessor2D : public LocalPostProcessor2D {
public:
ShanChenForcedPostProcessor2D(int x0_, int x1_, int y0_, int y1_, T G_,
std::vector rho0_, AnalyticalF1D& iP_,
std::vector partners_);
ShanChenForcedPostProcessor2D(T G_,
std::vector rho0_, AnalyticalF1D& iP_,
std::vector partners_);
int extent() const override
{
return 1;
}
int extent(int whichDirection) const override
{
return 1;
}
void process(BlockLattice2D& blockLattice) override;
void processSubDomain(BlockLattice2D& blockLattice,
int x0_, int x1_, int y0_, int y1_) override;
private:
int x0, x1, y0, y1;
T G;
std::vector rho0;
AnalyticalF1D& interactionPotential;
std::vector partners;
};
template
class ShanChenForcedGenerator2D : public LatticeCouplingGenerator2D {
public:
ShanChenForcedGenerator2D(int x0_, int x1_, int y0_, int y1_, T G_, std::vector rho0_, AnalyticalF1D& iP_);
ShanChenForcedGenerator2D(T G_, std::vector rho0_, AnalyticalF1D& iP_);
PostProcessor2D* generate(std::vector partners) const override;
LatticeCouplingGenerator2D* clone() const override;
private:
T G;
std::vector rho0;
AnalyticalF1D& interactionPotential;
};
}
#endif