FuncDef.h

Go to the documentation of this file.

 
#define  pi 3.14159265359
extern char maindir[125];
 
extern double porosity;
 
extern double density;
 
extern unsigned long int timesteps;
 
extern double thickness;
 
extern double saturation;
 
extern double timeunit;
 
extern char controlfile[120];
 
extern double totalFluxIn;
 
/* parts_fracture: initial number of particles per fracture */
/* Global variables: nnodes - number of nodes     */
/* ncells - number of triangles */
/* nfract - number of fractures */
/* max_neighb - maximum number of edges in Voronoi polygon, node connections */
 
extern unsigned int pflotran;
 
extern unsigned int fehm;
 
extern unsigned int nnodes;
 
extern unsigned int ncells;
 
extern unsigned int nfract;
 
extern unsigned int max_neighb;
 
extern unsigned int npart;
 
extern unsigned int np;
 
extern unsigned int nzone_in;
 
extern unsigned int *nodezonein;
 
extern unsigned int *nodezoneout;
 
extern unsigned int flag_w;
 
 
 
struct material {
 
    float theta;
    
    float nvect_xy[2];
    
    float nvect_z[2];
    
    unsigned int firstnode;
    
    unsigned int lastnode;
    
    unsigned int firstcell;
    
    unsigned int numbcells;
    
    double rot2mat[3][3];
    
    double rot3mat[3][3];
};
 
 
struct vertex {
 
    unsigned int typeN;
    
    unsigned int fracture[2];
    
    double coord[3];
    
    double coord_xy[6];
    
    double pvolume;  // volume of Voronoi polygon
    
    double pressure; //pressure
    
    unsigned int numneighb;  //number of neighboring nodes
    
    double velocity[4][2]; //Darcy's velocity
    
    double timestep[4]; // time step
    
    unsigned int* indnodes; //array of neighboring node's indices
    
    unsigned int** cells; //array of neighboring cells's indices
    
    unsigned int** fracts;//array of fracture's numbers
    
    unsigned int* type; //array of node's type
    
    double* flux;//array of fluxes
    
    double* area;//area of Voronoi facea
    
    double residtime; //residential time of particles
    
    double aperture; //aperture of the cell
};
 
 
/* element structure contains the information about trianglar cells in DFN mesh */
struct element {
    unsigned int fracture;
    
    unsigned int node_ind[3];// three nodes indices
    
    unsigned int veloc_ind[3]; // index of velocity of each node
    
    
};
 
 
struct contam {
 
    double velocity[2];
    
    double position[2];
    
    double prev_pos[2];
    
    unsigned int fracture;
    
    unsigned int cell;
    
    unsigned int intcell;
    
    double weight[3];
    
    double time;
    
    double fl_weight;
    
    double pressure; //fluid pressure
    
    double t_adv_diff; //in case of TDRW - total travel time
    
    double t_diff; // in case of TDRW - accumulative diffusion time
};
 
extern    struct material *fracture;
 
extern    struct vertex *node;
 
extern    struct contam *particle;
 
extern    struct element *cell;
 
 
void ReadInit();
void ReadDataFiles();
void AdjacentCells(int  ln, int i, int  j, int  k);
void Convertto2d();
void Convertto3d();
void DarcyVelocity();
struct matr   MatrixProducts (double normxarea[][2],  int number);
struct lb  DefineBoundaryAngle(int i, unsigned int edge_1, unsigned int edge_2, int f1, int coorf);
void VelocityInteriorNode (double normxarea[][2], int i, int number, unsigned int indj[max_neighb], int vi);
void VelocityExteriorNode (double normxarea[][2], int i, int number, unsigned int indj[max_neighb], struct lb lbound, int vi) ;
void CheckNewCell();
void ParticleTrack();
struct intcoef  CalculateWeights(int nn1, int nn2, int nn3);
void SearchNeighborCells(int nn1, int nn2, int nn3);
int InsideCell (unsigned int numc);
void NeighborCells (int k);
void PredictorStep();
void CorrectorStep();
void DefineTimeStep();
int CheckDistance();
void AcrossIntersection (int prevcell, int int1, int int2, int mixing_rule);
void ChangeFracture(int cell_win);
struct posit3d CalculatePosition3D();
int InitCell ();
int InitPos();
void Moving2Center(int nnp, int cellnumber);
int Moving2NextCell(int stuck, int k);
void BoundaryCells();
int InitParticles_np (int k_current, int firstn, int lastn, int parts_fracture, int first_ind, int last_ind);
int InitParticles_eq (int k_current, int firstn, int lastn, double parts_dist, int first_ind, int last_ind);
int CornerVelocity(int i, int m1, int m2, int m3, int s1, int s2, int s3);
void ReadBoundaryNodes();
FILE *OpenFile(char filen[120], char fileopt[2]);
void ReadFEHMfile(int nedges);
void ReadPFLOTRANfile(int nedges);
void WritingInit();
void Velocity3D();
double CalculateCurrentDT();
int Yindex(int nodenum, int np);
int Xindex(int nodenum, int np);
int CompleteMixingRandomSampling(double products[4], double speedsq[4], int indj, int int1, int indk);
int StreamlineRandomSampling(double products[4], double speedsq[4], int indj, int int1, int indk, int neighborcellind[4], int neighborfracind[4], int prevfrac, int prevcell);
void OutputVelocities();
int XindexC(int nodenum, int ii);
int YindexC(int nodenum, int ii);
double DefineAngle(double u1, double u2, double v1, double v2);
void HalfPolygonVelocity(int i, int k, int fractn, int indc, unsigned int fractj[max_neighb]);
struct posit3d CalculateVelocity3D();
void BoundaryLine(int n1, int n2, int n3);
void CheckGrid();
int BVelocityDirection(int b1, int b2);
double InOutFlowCell(int indcell, int int1, double nposx, double nposy);
void Moving2NextCellBound(int prevcell);
struct inpfile Control_File(char fileobject[], int ctr);
struct inpfile Control_Data(char fileobject[], int ctr);
void ParticleOutput (int currentt, int frac_p);
struct inpfile Control_Param(char fileobject[], int ctr);
void FlowInWeight(int numbpart);
int InitParticles_ones (int k_current, double inter_p[][4], int fracture, int parts_fracture, int ii, double thirdcoor, int zonenumb_in, int first_ind, int last_ind);
void Coordinations2D ();
void ReadAperture();
void InitInMatrix();
double TimeFromMatrix(double pdist);
void FinalPosition();
struct lagrangian CalculateLagrangian(double xcurrent, double ycurrent, double zcurrent, double xprev, double yprev, double zprev);
void OutputMarPlumDisp (int currentnum, char path[125]);
int String_Compare(char string1[], char string2[]);
struct inpfile Control_File_Optional(char fileobject[], int ctr);
double TimeDomainRW (double time_advect);
int InitParticles_flux (int k_current, int firstn, int lastn, double weight_p);
int InitInWell(int nodepart);