mesh.h

用于快速的估计任意形状的物体的运动情况受力分子

/*************************************************************************************************/
/*      变量的声明                                                                               */
/*-----------------------------------------------------------------------------------------------*/
/*      struct mesh_data			  网格数据                                                   */
/*************************************************************************************************/


#ifndef MESH_H
#define MESH_H

#include <math.h>
#include <stdio.h>
#include <malloc.h>

#define PI 3.1415926535


//---------------------------------------------------------------------
//      任意外形物体网格数据

struct mesh
{
    int NBLOCK;                             // Block的数目
    int *IMAX, *JMAX, *KMAX;                // 每个Block在3个方向的节点数
    int gridnum;                            // 节点数目
    int facenum;                            // 面元数目
    double *gridx, *gridy, *gridz;          // 每个节点的(x,y,z)坐标
    int *node1, *node2, *node3, *node4;     // 每个面元四个顶点的编号
    double xmax, xmin;	                      // x坐标的最大值与最小值
    double ymax, ymin;                      // y坐标的最大值与最小值
    double zmax, zmin;                      // z坐标的最大值与最小值
};

int trans_calc_force(double, double, double, double, double, double, double, double, double, mesh, double *, double *);
int env_calc_atmos(double, double *, double *, double *, double *, double *);
double env_interp(double *,double *,int,double);
double erf(double);

#endif