boundaryconditions.h

Go to the documentation of this file.

#ifndef _BC_H_
#define _BC_H_

#include <basic.h>

#define _PERIODIC_                      "periodic"

#define _EXTRAPOLATE_                   "extrapolate"

#define _DIRICHLET_                     "dirichlet"

#define _REFLECT_                       "reflect"

#define _SPONGE_                        "sponge"

/* some BC types unique to the euler/navier-stokes systems */
#define _NOSLIP_WALL_                   "noslip-wall"

#define _THERMAL_NOSLIP_WALL_           "thermal-noslip-wall"

#define _SLIP_WALL_                     "slip-wall"

#define _THERMAL_SLIP_WALL_             "thermal-slip-wall"

#define _SUBSONIC_INFLOW_               "subsonic-inflow"

#define _SUBSONIC_OUTFLOW_              "subsonic-outflow"

#define _SUBSONIC_AMBIVALENT_           "subsonic-ambivalent"

#define _SUPERSONIC_INFLOW_             "supersonic-inflow"

#define _SUPERSONIC_OUTFLOW_            "supersonic-outflow"

#define _TURBULENT_SUPERSONIC_INFLOW_   "turbulent-supersonic-inflow"

/* some BC types unique to the NUMA system */
#define _NO_FLUX_BC_                    "numa-nfbc"

/* some BC types unique to the shallow water system */
#define _SW_SLIP_WALL_                  "shallow-water-slip-wall"

#define _SW_NOSLIP_WALL_                "shallow-water-noslip-wall"

typedef struct domain_boundaries {

  char    bctype [_MAX_STRING_SIZE_];
  int     dim;
  int     face;
  double  *xmin;
  double  *xmax;

  int on_this_proc;   
  int *is, *ie;       
  int (*BCFunctionU) (void*,void*,int,int,int*,int,double*,double);

  double *DirichletValue;   
  double *SpongeValue;      
  int    *UnsteadyDirichletSize; 
  double *UnsteadyDirichletData; 
  char    UnsteadyDirichletFilename[_MAX_STRING_SIZE_];

  /* variables specific to Navier-Stokes/Euler equations BCs */
  double gamma,                                   
         FlowDensity,                             
         *FlowVelocity,                           
         FlowPressure;                            
  /* variables specific to the thermal slip-wall boundary condition */
  int    *UnsteadyTemperatureSize; 
  double *UnsteadyTimeLevels;      
  double *UnsteadyTemperatureData; 
  char  UnsteadyTemperatureFilename[_MAX_STRING_SIZE_];

#if defined(HAVE_CUDA)
  int gpu_npoints_bounds;
  int gpu_npoints_local_wghosts;
  int *gpu_bounds, *gpu_is, *gpu_ie;
  double *gpu_FlowVelocity;
#endif

} DomainBoundary;

/* Functions */
int BCInitialize(void*, int); 
int BCCleanup   (void*, int); 
/* Boundary condition implementations for the solution vector U */
int BCPeriodicU                     (void*,void*,int,int,int*,int,double*,double);
int BCExtrapolateU                  (void*,void*,int,int,int*,int,double*,double);
int BCDirichletU                    (void*,void*,int,int,int*,int,double*,double);
int BCReflectU                      (void*,void*,int,int,int*,int,double*,double);
int BCNoslipWallU                   (void*,void*,int,int,int*,int,double*,double);
int BCThermalNoslipWallU            (void*,void*,int,int,int*,int,double*,double);
int BCSlipWallU                     (void*,void*,int,int,int*,int,double*,double);
int BCThermalSlipWallU              (void*,void*,int,int,int*,int,double*,double);
int BCSubsonicInflowU               (void*,void*,int,int,int*,int,double*,double);
int BCSubsonicOutflowU              (void*,void*,int,int,int*,int,double*,double);
int BCSubsonicAmbivalentU           (void*,void*,int,int,int*,int,double*,double);
int BCSupersonicInflowU             (void*,void*,int,int,int*,int,double*,double);
int BCSupersonicOutflowU            (void*,void*,int,int,int*,int,double*,double);
int BCTurbulentSupersonicInflowU    (void*,void*,int,int,int*,int,double*,double);
int BCNoFluxU                       (void*,void*,int,int,int*,int,double*,double);
int BCSWSlipWallU                   (void*,void*,int,int,int*,int,double*,double);

int BCSpongeSource        (void*,int,int,int,int*,double*,double*,double*);
int BCSpongeUDummy        (void*,void*,int,int,int*,int,double*,double);

int BCReadTemperatureData(void*,void*,int,int,int*);
int BCReadTurbulentInflowData(void*,void*,int,int,int*);

#if defined(HAVE_CUDA)

# if __cplusplus
extern "C" {
# endif

int gpuBCPeriodicU (void*,void*,int,int,int*,int,double*,double);
int gpuBCSlipWallU (void*,void*,int,int,int*,int,double*,double);

# if __cplusplus
}
# endif

#endif

#endif