pdd_plast_flow_supg/include/pdh_plast_flow_supg_bc.h

Go to the documentation of this file.

/************************************************************************
File pdh_plast_flow_bc.h - types and functions related to boundary conditions 
               for Navier-Stokes SUPG approximation
               
Contains definition of types used by bcs               

Contains declarations of routines:
  pdr_plast_flow_bc_read - read bc data from config file
  pdr_plast_flow_bc_free - free bc resources
  pdr_plast_flow_get_bc_count - get num of boundaries with conditions set in file
  pdr_plast_flow_get_stress_pins_count
  pdr_plast_flow_get_velocity_pins_count
  pdr_plast_flow_get_stress_pin - get stress pin data
  pdr_plast_flow_get_velocity_pin - get velocity pin data
  pdr_plast_flow_get_bc_type - get type of flow bc for boundary
  pdr_plast_flow_get_bc_data - get flow bc data
  pdr_plast_flow_update_timedep_bc - update timedependent boundary conditions

------------------------------
History:
    2015    - Krzysztof Banas (pobanas@cyf-kr.edu.pl) - based on ns_supg
    2016    - Aleksander Siwek (Aleksander.Siwek@agh.edu.pl)

*************************************************************************/

#ifndef PDH_PLAST_FLOW_BC_H
#define PDH_PLAST_FLOW_BC_H

#include <stdio.h>

#ifdef __cplusplus
extern "C" 
{
#endif

/**************************************/
/* TYPES                              */
/**************************************/
/* Rules:
/* - type name starts always with pdt_ */

/* types of flow bcs */
typedef enum {
  BC_PLAST_FLOW_NONE,
  BC_PLAST_FLOW_NOSLIP,
  BC_PLAST_FLOW_FREE,
  BC_PLAST_FLOW_FREE_SLIP,
  BC_PLAST_FLOW_SYMMETRY,
  BC_PLAST_FLOW_VELOCITY,
  BC_PLAST_FLOW_TRACTION, //TODO  
  BC_PLAST_FLOW_FRICTION // TODO
} pdt_plast_flow_bctype;

  typedef enum {
    BC_PLAST_FLOW_VEL_ORIENT_GLOB=0,
    BC_PLAST_FLOW_VEL_ORIENT_LOC=1
  } pdt_plast_flow_bc_vel_orient;
  
/* FLOW BCs and pins */

/*
every structure contains data related to
particular flow bc
*/


//noslip - don't need struct for noslip
//symm - don't need struct for symm

typedef struct {
  int bnum;
  double v[3];
  int vel_orient;
} pdt_plast_flow_bc_velocity;

typedef struct {
  int bnum;
  double t[3];
} pdt_plast_flow_bc_traction;

typedef struct {
  int bnum;
  double n[3];   //versor (normalized vector) normal to the (>free<)surface
  // ??? - friction type, friction coefficient, etc. - TODO
} pdt_plast_flow_bc_friction;

typedef struct {
  int bnum;
  double v[3];
  int vel_orient;
} pdt_plast_flow_bc_symmetry;

typedef struct {
  int bnum;
  double v[3];
  int vel_orient;
} pdt_plast_flow_bc_free_slip;

typedef struct {
  double p;
  double pin_node_coor[3];
} pdt_plast_flow_pin_stress;

typedef struct {
  double v[3];
  double pin_node_coor[3]; 
} pdt_plast_flow_pin_velocity;



/* utility types - not for direct access!*/

typedef struct {
  int bnum;
  pdt_plast_flow_bctype bc_plast_flow;
  int bc_plast_flow_data_idx;
} pdt_plast_flow_bc_assignments;

/* main structure containing bc data */
typedef struct {
  pdt_plast_flow_bc_assignments *bc_assignments;
  int bc_assignments_count;
  
  pdt_plast_flow_bc_velocity *bc_velocity;
  int bc_velocity_count;
  
  pdt_plast_flow_bc_traction *bc_traction;
  int bc_traction_count;

  pdt_plast_flow_bc_friction *bc_friction;
  int bc_friction_count;

  pdt_plast_flow_bc_symmetry *bc_symmetry;
  int bc_symmetry_count;
  
  pdt_plast_flow_bc_free_slip *bc_free_slip;
  int bc_free_slip_count;
  
  pdt_plast_flow_pin_stress *pin_stress;
  int pin_stress_count;

  pdt_plast_flow_pin_velocity *pin_velocity;
  int pin_velocity_count;


} pdt_plast_flow_bc;


/**************************************/
/* INTERNAL PROCEDURES                */
/**************************************/
/* Rules:
/* - name always begins with pdr_ */
/* - argument names start uppercase */

int pdr_plast_flow_bc_read(
  char *Work_dir,
  char *Filename,
  FILE *Interactive_output,
  pdt_plast_flow_bc *Bc_db);

int pdr_plast_flow_bc_free(pdt_plast_flow_bc *Bc_db);

int pdr_plast_flow_get_stress_pins_count(const pdt_plast_flow_bc *Bc_db); 
  /* in: pdt_bc structure to read from */

int pdr_plast_flow_get_velocity_pins_count(const pdt_plast_flow_bc *Bc_db); 
  /* in: pdt_bc structure to read from */

int pdr_plast_flow_get_bc_assign_count(const pdt_plast_flow_bc *Bc_db); 
  /* in: pdt_bc structure to read from */

pdt_plast_flow_pin_stress* pdr_plast_flow_get_stress_pin(
  const pdt_plast_flow_bc *Bc_db, /* in: pdt_bc structure to read from */
  int Idx  /* in: pin stress no. */
  );

pdt_plast_flow_pin_velocity* pdr_plast_flow_get_velocity_pin(
  const pdt_plast_flow_bc *Bc_db, /* in: pdt_bc structure to read from */
  int Idx /* in: pin stress no. */
  );

pdt_plast_flow_bctype pdr_plast_flow_get_bc_type(
  const pdt_plast_flow_bc *Bc_db, /* in: pdt_bc structure to read from */
  int Bnum /* in: boundary number */
  );


void* pdr_plast_flow_get_bc_data(
  const pdt_plast_flow_bc *Bc_db, /* in: pdt_bc structure to read from */
  int Bnum /* in: boundary number */
  );


int pdr_plast_flow_update_timedep_bc(
  const pdt_plast_flow_bc *Bc_db, /* in: pdt_bc structure to read from */
  double Time
  );


#ifdef __cplusplus
}
#endif

#endif

---