IBManager.h

Go to the documentation of this file.

// This file os part of FVM
// Copyright (c) 2012 FVM Authors
// See LICENSE file for terms.

#ifndef _IBMANAGER_H_
#define _IBMANAGER_H_

#include "Mesh.h"
#include "GeomFields.h"
#include <set>

class KSearchTree;
class AABB;

struct
NearestCell
{
  NearestCell():
    mesh(0),
    cell(-1),
    distanceSquared(0)
  {}

  const Mesh* mesh;
  int cell;
  double distanceSquared;
  vector<int> neighbors;
};
struct
NearestIBFace
{
  NearestIBFace():
    mesh(0),
    IBFace(-1),
    distanceSquared(0)
  {}

  const Mesh* mesh;
  int IBFace;
  double distanceSquared;
  vector<int> neighbors;
};
struct   
doubleIntStruct{
   double VALUE;
   int TAG;
};



class IBManager
{
public:
  typedef Vector<double,3> Vec3D;
  typedef Array<Vec3D> Vec3DArray;
  typedef Array<int> IntArray;

  IBManager(GeomFields& geomFields,
            Mesh& solidBoundaryMesh,
            const MeshList& fluidMeshes);

  void update();

  int fluidNeighborsPerIBFace;
  int fluidNeighborsPerSolidFace;
  int solidNeighborsPerIBFace;
  int IBNeighborsPerSolidFace;
  Mesh& _solidBoundaryMesh;
  
private:

  void markIntersections(Mesh& fluidMesh, AABB& sMeshesAABB);
  int markFluid(Mesh& fluidMesh);
  int markSolid(Mesh& fluidMesh);
  void markIBTypePlus(Mesh& fluidMesh);
  void createIBFaces(Mesh& fluidMesh);
  void createIBInterpolationStencil(Mesh& mesh,
                                    KSearchTree& fluidCellsTree,
                                    KSearchTree& solidFacesTree);
  void
  findNearestCellForSolidFaces(Mesh& mesh,
                               KSearchTree& fluidCellsTree,
                               vector<NearestCell>& nearest);
  //void
  //indNearestIBFaceForSolidFaces(Mesh& mesh,
  //                            KSearchTree& IBFacesTree,
  //                            vector<NearestIBFace>& nearestIB);
  
  void createSolidInterpolationStencil(Mesh& mesh,
                                       KSearchTree& IBFacesTree,
                                       vector<NearestCell>& nearest);
                                       // vector<NearestIBFace>& nearestIB);
  void CRConnectivityPrintFile(const CRConnectivity& conn, const string& name, const int procID) const;

  GeomFields& _geomFields;
  const MeshList _fluidMeshes;
};
#endif