smoothIndicatorF2D.h

Go to the documentation of this file.

/*  This file is part of the OpenLB library
 *
 *  Copyright (C) 2016 Thomas Henn, Cyril Masquelier, Jan Marquardt, Mathias J. Krause
 *  E-mail contact: info@openlb.net
 *  The most recent release of OpenLB can be downloaded at
 *  <http://www.openlb.net/>
 *
 *  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 SMOOTH_INDICATOR_F_2D_H
#define SMOOTH_INDICATOR_F_2D_H
 
#include <vector>
 
#include "smoothIndicatorBaseF2D.h"
#include "io/xmlReader.h"
#include "utilities/functorPtr.h"
 
#include "core/blockData.h"
#include "core/unitConverter.h"
#include "smoothIndicatorBaseF3D.h"
#include "functors/analytical/indicator/indicatorBaseF2D.h"
#include "functors/analytical/indicator/indicatorBaseF3D.h"
#include "functors/analytical/indicator/indicatorF2D.h"
 
namespace olb {
 
 
 
template <typename T, typename S, bool PARTICLE=false>
class SmoothIndicatorCuboid2D final : public SmoothIndicatorF2D<T,S,PARTICLE> {
private:
  IndicatorCuboid2D<S> _ind;
public:
  SmoothIndicatorCuboid2D(IndicatorCuboid2D<S>& ind, S epsilon, S theta=0);
  SmoothIndicatorCuboid2D(Vector<S,2>center, S xLength, S yLength, S epsilon, S theta=0);
  IndicatorCuboid2D<S>& getIndicator();
  Vector<S,2> surfaceNormal(const Vector<S,2>& pos, const S meshSize) override;
  const S signedDistance( const Vector<S,2> input ) override;
  bool distance(S& distance, const Vector<S,2>& origin, const Vector<S,2>& direction, S precision, S pitch) override;
  S getArea( ) override;
  Vector<S,2> calcMofiAndMass(const S density) override;
};
 
template <typename T, typename S, bool PARTICLE=false>
class SmoothIndicatorCircle2D final : public SmoothIndicatorF2D<T,S,PARTICLE> {
private:
  IndicatorCircle2D<S> _ind;
public:
  SmoothIndicatorCircle2D(IndicatorCircle2D<S>& indPtr, S epsilon);
  SmoothIndicatorCircle2D(Vector<S,2>center, S radius, S epsilon);
  IndicatorCircle2D<S>& getIndicator();
  Vector<S,2> surfaceNormal(const Vector<S,2>& pos, const S meshSize) override;
  const S signedDistance( const Vector<S,2> input ) override;
  bool distance(S& distance, const Vector<S,2>& origin, const Vector<S,2>& direction, S precision, S pitch) override;
  S getArea( ) override;
  Vector<S,2> calcMofiAndMass(const S density) override;
};
 
template <typename T, typename S, bool PARTICLE=false>
class SmoothIndicatorTriangle2D final : public SmoothIndicatorF2D<T,S,PARTICLE> {
private:
  IndicatorEquiTriangle2D<S> _ind;
public:
  SmoothIndicatorTriangle2D(IndicatorEquiTriangle2D<S>& indPtr, S epsilon, S theta=0);
  SmoothIndicatorTriangle2D(Vector<S,2>center, S radius, S epsilon, S theta=0);
  IndicatorEquiTriangle2D<S>& getIndicator();
  Vector<S,2> surfaceNormal(const Vector<S,2>& pos, const S meshSize) override;
  const S signedDistance( const Vector<S,2> input ) override;
  bool distance(S& distance, const Vector<S,2>& origin, const Vector<S,2>& direction, S precision, S pitch) override;
  S getArea( ) override;
  Vector<S,2> calcMofiAndMass(const S density) override;
};
 
/*
template <typename T, typename S, bool PARTICLE=false>
class SmoothIndicatorTriangle2D final : public SmoothIndicatorF2D<T,S,PARTICLE> {
private:
  Vector<S, 2> _PointA, _PointB, _PointC;
  Vector<S, 2> _ab, _bc, _ca;
  S _ab_d, _bc_d, _ca_d;
public:
  SmoothIndicatorTriangle2D(Vector<S,2> center, S radius, S epsilon, S theta=0, S density=0, Vector<S,2> vel = Vector<S,2> (0.,0.), S omega = 0);
  bool operator() (T output[], const S input[]) override;
};
*/
 
template <typename T, typename W>
class AnalyticalFfromBlockF2D;
template <typename T, typename BaseType>
class BlockDataF2D;
 
//implements a custom shaped smooth particle //TODO: Check for consistency
//ALSO: adap .hh
template <typename T, typename S, bool PARTICLE=false>
class SmoothIndicatorCustom2D : public SmoothIndicatorF2D<T, S, PARTICLE> {
private:
  std::shared_ptr<IndicatorF2D<T>> _indPtr;
  const T _latticeSpacing;
  // _center is the local center
  Vector<T,2> _center;
  std::unique_ptr<Cuboid2D<T>> _cuboid;
  std::unique_ptr<BlockData<2,T,BaseType<T>>> _blockData;
  std::unique_ptr<BlockDataF2D<T,BaseType<T>>> _cacheFunctor;
  std::unique_ptr<AnalyticalFfromBlockF2D<T,T>> _interpolateCache;
  T _area;
 
  void initData(IndicatorF2D<T>& ind);
  void initBlockData(IndicatorF2D<T>& ind);
  void calcCenter();
  void calcCircumRadius();
 
public:
  SmoothIndicatorCustom2D(T latticeSpacing,
                          std::shared_ptr<IndicatorF2D<T>> indPtr, Vector<T,2> pos, T epsilon,
                          T theta=0.);
 
  Vector<T,2> getLocalCenter();
  S getArea( ) override;
  Vector<T,2> calcMofiAndMass(T rhoP) override;
  Vector<S,2> surfaceNormal(const Vector<S,2>& pos, const S meshSize) override;
  const S signedDistance( const Vector<S,2> input ) override;
  bool regardCell(int input[2]);
  bool operator()(T output[], const S input[]) override;
};
 
//Geng2019:
template <typename T, typename S, bool PARTICLE=false>
class SmoothIndicatorHTCircle2D final : public SmoothIndicatorF2D<T,S,PARTICLE> {
private:
  S _radius;
public:
  SmoothIndicatorHTCircle2D(Vector<S,2> center, S radius, S epsilon, S density=0, Vector<S,2> vel = Vector<S,2> (0.,0.), S omega = 0);
  bool operator() (T output[], const S input[]) override;
  Vector<S,2> calcMofiAndMass(const S density) override;
};
 
template <typename T, typename S, bool PARTICLE=false>
class SmoothIndicatorFactoredCircle2D final : public SmoothIndicatorF2D<T,S,PARTICLE> {
private:
  S _radius;
  S _factor;
public:
  SmoothIndicatorFactoredCircle2D(Vector<S,2> center, S radius, S epsilon, S density=0, Vector<S,2> vel = Vector<S,2> (0.,0.), S omega = 0, S factor = 1.);
  bool operator() (T output[], const S input[]) override;
  Vector<S,2> calcMofiAndMass(const S density) override;
};
 
template <typename T, typename S, bool PARTICLE=false>
class SmoothIndicatorFactoredCuboid2D final : public SmoothIndicatorF2D<T,S,PARTICLE> {
private:
  S _xLength;
  S _yLength;
  S _factor;
public:
  SmoothIndicatorFactoredCuboid2D(Vector<S,2> center, S xLength, S yLength, S epsilon, S density=0, Vector<S,2> vel = Vector<S,2> (0.,0.), S omega = 0, S factor = 1.);
  bool operator() (T output[], const S input[]) override;
  Vector<S,2> calcMofiAndMass(const S density) override;
};
 
}
 
#endif