immersedboundaries.h

Go to the documentation of this file.

#ifndef _IB_H_
#define _IB_H_

#define _IB_NDIMS_ 3

#define _IB_NNODES_ 8

#define _IB_XY_ "2d (xy)"

#define _IB_XZ_ "2d (xz)"

#define _IB_YZ_ "2d (yz)"

#define _IB_3D_ "3d"

#include <basic.h>

typedef struct _facet_3d_{
    double x1, 
         x2, 
         x3, 
           y1, 
         y2, 
         y3, 
           z1, 
         z2, 
         z3, 
           nx, 
         ny, 
         nz; 
} Facet3D;

typedef struct _facet_map_{
  Facet3D   *facet; 
  int       index;  
  int       interp_nodes    [_IB_NNODES_];
  double    interp_coeffs   [_IB_NNODES_];
  int       interp_nodes_ns [_IB_NNODES_];
  double    interp_coeffs_ns[_IB_NNODES_];
  double    xc, 
            yc, 
            zc, 
            xns,
            yns,
            zns;
  double    dx, 
            dy, 
            dz; 
} FacetMap;

typedef struct _body_3d_{
  int     nfacets;    
  Facet3D *surface;   
  /* coordinates of bounding box */
  double xmin, 
         xmax, 
         ymin, 
         ymax, 
         zmin, 
         zmax; 
} Body3D;

typedef struct _boundary_node_{
  int     i,      
          j,      
          k,      
          p;      
  double  x,      
          y,      
          z;      
  Facet3D *face;  
  int     interp_nodes[_IB_NNODES_];   
  double  interp_coeffs[_IB_NNODES_],  
          interp_node_distance,        
          surface_distance;            
} IBNode;

typedef struct immersed_boundary{
  Body3D    *body;      
  IBNode    *boundary;  
  FacetMap  *fmap;      
  double  tolerance; 
  double  delta;     
  int     itr_max;   
  int     n_boundary_nodes; 
  int     nfacets_local;    
  char    mode[_MAX_STRING_SIZE_]; 
} ImmersedBoundary;


int IBReadBodySTL (Body3D**,char*,void*,int*);
int IBWriteBodySTL(Body3D*,char*,void*,int,int*);

int IBCleanup           (void*);
int IBComputeBoundingBox(Body3D*);
int IBCreateFacetMapping(void*,void*,double*,int*,int);
int IBIdentifyBody      (void*,int*,int*,int,void*,double*,double*);
int IBIdentifyBoundary  (void*,void*,int*,int,double*);
int IBIdentifyMode      (double*,int*,void*);
int IBNearestFacetNormal(void*,void*,double*,double,int*,int);
int IBInterpCoeffs      (void*,void*,double*,int*,int,double*);

int IBAssembleGlobalFacetData(void*,void*,const double* const, double** const,int);

int IBComputeNormalGradient(void*,void*,const double* const, int, double** const);
int IBComputeFacetVar(void*,void*,const double* const, int, double** const);

#endif