preInlet.h

Go to the documentation of this file.

/*
This file is part of the HemoCell library
 
HemoCell is developed and maintained by the Computational Science Lab
in the University of Amsterdam. Any questions or remarks regarding this library
can be sent to: info@hemocell.eu
 
When using the HemoCell library in scientific work please cite the
corresponding paper: https://doi.org/10.3389/fphys.2017.00563
 
The HemoCell library is free software: you can redistribute it and/or
modify it under the terms of the GNU Affero General Public License as
published by the Free Software Foundation, either version 3 of the
License, or (at your option) any later version.
 
The library 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 Affero General Public License for more details.
 
You should have received a copy of the GNU Affero General Public License
along with this program.  If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef PREINLET_H
#define PREINLET_H
 
class PreInlet;
#include "config/constant_defaults.h"
#include "core/hemoCellFunctional.h"
 
#include "core/geometry3D.h"
#include "atomicBlock/atomicBlock3D.h"
#include "config.h"
#include "hemocell.h"
 
#ifndef HEMOCELL_H
namespace hemo {
  class HemoCell;
}
#endif
 
enum Direction : int {
  Xpos, Xneg, Ypos, Yneg, Zpos, Zneg
};
 
 
#define DSET_SLICE 1000
 
namespace hemo {
 
 
inline plint cellsInBoundingBox(plb::Box3D const & box) {
  return abs((box.x1 - box.x0)*(box.y1-box.y0)*(box.z1-box.z0));
}
 
 
class PreInlet {
public:
    class CreateDrivingForceFunctional: public HemoCellFunctional {
    void processGenericBlocks(plb::Box3D, std::vector<plb::AtomicBlock3D*>);
    CreateDrivingForceFunctional * clone() const;
    plint & count;
    public:
      CreateDrivingForceFunctional(plint & count_) :
                                count(count_) {count = 0;}
  };
  class CreatePreInletBoundingBox: public HemoCellFunctional {
    void processGenericBlocks(plb::Box3D, std::vector<plb::AtomicBlock3D*>);
    CreatePreInletBoundingBox * clone() const;
    plb::Box3D & boundingBox;
    bool & foundPreInlet;
    public:
      CreatePreInletBoundingBox(plb::Box3D & b_, bool & fp_) :
                                boundingBox(b_), foundPreInlet(fp_) {}
  };
 
  class FillFlagMatrix: public HemoCellFunctional {
    void processGenericBlocks(plb::Box3D, std::vector<plb::AtomicBlock3D*>);
    FillFlagMatrix * clone() const;
  };
 
  PreInlet(hemo::HemoCell * hemocell_, plb::MultiScalarField3D<int> * flagMatrix_);
  PreInlet(hemo::HemoCell * hemocell_, plb::MultiBlockManagement3D & management);
  inline plint getNumberOfNodes() { return cellsInBoundingBox(location);}
  void createBoundary();
  bool readNormalizedVelocities();
  void setDrivingForce();
  void setDrivingForceTimeDependent(double t);
  void calculateDrivingForce();
  double interpolate(vector<double> &xData, vector<double> &yData, double x, bool extrapolate);
  double average(vector<double> values);
  void applyPreInletVelocityBoundary();
  void applyPreInletParticleBoundary();
  void applyPreInlet() { applyPreInletVelocityBoundary();
                         applyPreInletParticleBoundary(); };
  void initializePreInletParticleBoundary();
  void initializePreInletVelocityBoundary();
  void initializePreInlet() { initializePreInletVelocityBoundary(); initializePreInletParticleBoundary(); };
 
  void autoPreinletFromBoundary(Direction);
  void preInletFromSlice(Direction direction_, Box3D boundary);
 
  Direction direction = Direction::Zneg;
  plb::Box3D location;
  plb::Box3D fluidInlet;
  int nProcs = 0;
  bool initialized = false;
  double drivingForce = 0.0;
  double average_vel = 0.0;
  double pulseEndTime = 1.0;
  double pFrequency = 1.0;
  std::vector<double> normalizedVelocityTimes;
  std::vector<double> normalizedVelocityValues;
  std::map<plint,plint> BlockToMpi;
  std::map<int,plint> particleSendMpi;
  std::map<int,bool> particleReceiveMpi;
  std::vector<plint> communicating_blocks;
  std::map<int, Box3D> domain_at_rank;
  std::vector<int> my_send_blocks;
  std::vector<int> my_recv_blocks;
  int sendingBlocks = 0;
  bool partOfpreInlet = false;
  int inflow_length = 0;
  int preinlet_length = 0;
  bool communications_mapped = false;
  std::vector<int> particle_receivers;
  std::vector<int> particle_senders;
  HemoCell * hemocell;
  MultiScalarField3D<int> *flagMatrix = 0;
};
 
}
#endif /* PREINLET_H */