euler2d.h File Reference

#include <basic.h>

Go to the source code of this file.

Data Structures
struct	Euler2D

Macros
#define	_EULER_2D_   "euler2d"
#define	_MODEL_NDIMS_   2
#define	_MODEL_NVARS_   4
#define	_ROE_   "roe"
#define	_RF_   "rf-char"
#define	_LLF_   "llf-char"
#define	_SWFS_   "steger-warming"
#define	_XDIR_   0
#define	_YDIR_   1
#define	_Euler2DGetFlowVar_(u, rho, vx, vy, e, P, p)
#define	_Euler2DSetFlux_(f, rho, vx, vy, e, P, p, dir)
#define	_Euler2DRoeAverage_(uavg, uL, uR, p)
#define	_Euler2DEigenvalues_(u, D, p, dir)
#define	_Euler2DLeftEigenvectors_(u, L, p, dir)
#define	_Euler2DRightEigenvectors_(u, R, p, dir)

Functions
int	Euler2DInitialize (void *, void *)
int	Euler2DCleanup (void *)

---

Data Structure Documentation

Euler2D

struct Euler2D

Definition at line 244 of file euler2d.h.

Data Fields
double	gamma
char	upw_choice[_MAX_STRING_SIZE_]

Macro Definition Documentation

_EULER_2D_

#define _EULER_2D_   "euler2d"

Definition at line 25 of file euler2d.h.

_MODEL_NDIMS_

#define _MODEL_NDIMS_   2

Definition at line 30 of file euler2d.h.

_MODEL_NVARS_

#define _MODEL_NVARS_   4

Definition at line 31 of file euler2d.h.

_ROE_

#define _ROE_   "roe"

Definition at line 34 of file euler2d.h.

_RF_

#define _RF_   "rf-char"

Definition at line 35 of file euler2d.h.

_LLF_

#define _LLF_   "llf-char"

Definition at line 36 of file euler2d.h.

_SWFS_

#define _SWFS_   "steger-warming"

Definition at line 37 of file euler2d.h.

_XDIR_

#define _XDIR_   0

Definition at line 40 of file euler2d.h.

_YDIR_

#define _YDIR_   1

Definition at line 41 of file euler2d.h.

_Euler2DGetFlowVar_

#define _Euler2DGetFlowVar_	(		u,
			rho,
			vx,
			vy,
			e,
			P,
			p
	)

Value:

{ \
    double  gamma = p->gamma, vsq; \
    rho = u[0]; \
    vx  = u[1] / rho; \
    vy  = u[2] / rho; \
    e   = u[3]; \
    vsq  = (vx*vx) + (vy*vy); \
    P   = (e - 0.5*rho*vsq) * (gamma-1.0); \
  }

Definition at line 44 of file euler2d.h.

_Euler2DSetFlux_

#define _Euler2DSetFlux_	(		f,
			rho,
			vx,
			vy,
			e,
			P,
			p,
			dir
	)

Value:

{ \
    if (dir == _XDIR_) { \
      f[0] = rho * vx; \
      f[1] = rho * vx * vx + P; \
      f[2] = rho * vx * vy; \
      f[3] = (e + P) * vx; \
    } else if (dir == _YDIR_) { \
      f[0] = rho * vy; \
      f[1] = rho * vy * vx; \
      f[2] = rho * vy * vy + P; \
      f[3] = (e + P) * vy; \
    } \
  }

Definition at line 55 of file euler2d.h.

_Euler2DRoeAverage_

#define _Euler2DRoeAverage_	(		uavg,
			uL,
			uR,
			p
	)

Definition at line 70 of file euler2d.h.

_Euler2DEigenvalues_

#define _Euler2DEigenvalues_	(		u,
			D,
			p,
			dir
	)

Value:

{ \
    double  gamma = p->gamma; \
    double  rho,vx,vy,e,P,c,vn,vsq; \
    rho = u[0]; \
    vx  = u[1] / rho; \
    vy  = u[2] / rho; \
    e   = u[3]; \
    vsq  = (vx*vx) + (vy*vy); \
    P   = (e - 0.5*rho*vsq) * (gamma-1.0); \
    c    = sqrt(gamma*P/rho); \
    if      (dir == _XDIR_) vn = vx; \
    else if (dir == _YDIR_) vn = vy; \
    else               vn = 0.0; \
    if (dir == _XDIR_) {\
      D[0*_MODEL_NVARS_+0] = vn-c;   D[0*_MODEL_NVARS_+1] = 0;    D[0*_MODEL_NVARS_+2] = 0;      D[0*_MODEL_NVARS_+3] = 0; \
      D[1*_MODEL_NVARS_+0] = 0;      D[1*_MODEL_NVARS_+1] = vn+c; D[1*_MODEL_NVARS_+2] = 0;      D[1*_MODEL_NVARS_+3] = 0; \
      D[2*_MODEL_NVARS_+0] = 0;      D[2*_MODEL_NVARS_+1] = 0;    D[2*_MODEL_NVARS_+2] = vn;     D[2*_MODEL_NVARS_+3] = 0; \
      D[3*_MODEL_NVARS_+0] = 0;      D[3*_MODEL_NVARS_+1] = 0;    D[3*_MODEL_NVARS_+2] = 0;      D[3*_MODEL_NVARS_+3] = vn; \
    } else if (dir == _YDIR_) { \
      D[0*_MODEL_NVARS_+0] = vn-c;   D[0*_MODEL_NVARS_+1] = 0;    D[0*_MODEL_NVARS_+2] = 0;      D[0*_MODEL_NVARS_+3] = 0; \
      D[1*_MODEL_NVARS_+0] = 0;      D[1*_MODEL_NVARS_+1] = vn;   D[1*_MODEL_NVARS_+2] = 0;      D[1*_MODEL_NVARS_+3] = 0; \
      D[2*_MODEL_NVARS_+0] = 0;      D[2*_MODEL_NVARS_+1] = 0;    D[2*_MODEL_NVARS_+2] = vn+c;   D[2*_MODEL_NVARS_+3] = 0; \
      D[3*_MODEL_NVARS_+0] = 0;      D[3*_MODEL_NVARS_+1] = 0;    D[3*_MODEL_NVARS_+2] = 0;      D[3*_MODEL_NVARS_+3] = vn; \
    }\
  }

Definition at line 108 of file euler2d.h.

_Euler2DLeftEigenvectors_

#define _Euler2DLeftEigenvectors_	(		u,
			L,
			p,
			dir
	)

Definition at line 135 of file euler2d.h.

_Euler2DRightEigenvectors_

#define _Euler2DRightEigenvectors_	(		u,
			R,
			p,
			dir
	)

Definition at line 189 of file euler2d.h.

Function Documentation

Euler2DInitialize()

int Euler2DInitialize	(	void *	,
		void *
	)

Definition at line 21 of file Euler2DInitialize.c.

{
  HyPar         *solver  = (HyPar*)          s;
  MPIVariables  *mpi     = (MPIVariables*)   m; 
  Euler2D       *physics = (Euler2D*) solver->physics;
  int           ferr     = 0;

  static int count = 0;

  if (solver->nvars != _MODEL_NVARS_) {
    fprintf(stderr,"Error in Euler2DInitialize(): nvars has to be %d.\n",_MODEL_NVARS_);
    return(1);
  }
  if (solver->ndims != _MODEL_NDIMS_) {
    fprintf(stderr,"Error in Euler2DInitialize(): ndims has to be %d.\n",_MODEL_NDIMS_);
    return(1);
  }

  /* default values */
  physics->gamma = 1.4; 
  strcpy(physics->upw_choice,"roe");

  /* reading physical model specific inputs - all processes */
  if (!mpi->rank) {
    FILE *in;
    if (!count) printf("Reading physical model inputs from file \"physics.inp\".\n");
    in = fopen("physics.inp","r");
    if (!in) printf("Warning: File \"physics.inp\" not found. Using default values.\n");
    else {
      char word[_MAX_STRING_SIZE_];
      ferr = fscanf(in,"%s",word); if (ferr != 1) return(1);
      if (!strcmp(word, "begin")){
          while (strcmp(word, "end")){
              ferr = fscanf(in,"%s",word); if (ferr != 1) return(1);
          if (!strcmp(word, "gamma")) { 
            ferr = fscanf(in,"%lf",&physics->gamma); if (ferr != 1) return(1);
          } else if (!strcmp(word,"upwinding")) {
            ferr = fscanf(in,"%s",physics->upw_choice); if (ferr != 1) return(1);
          } else if (strcmp(word,"end")) {
            char useless[_MAX_STRING_SIZE_];
            ferr = fscanf(in,"%s",useless); if (ferr != 1) return(ferr);
            printf("Warning: keyword %s in file \"physics.inp\" with value %s not ",word,useless);
            printf("recognized or extraneous. Ignoring.\n");
          }
        }
        } else {
          fprintf(stderr,"Error: Illegal format in file \"physics.inp\".\n");
        return(1);
        }
    }
    fclose(in);
  }

#ifndef serial
  IERR MPIBroadcast_double    (&physics->gamma,1,0,&mpi->world);                      CHECKERR(ierr);
  IERR MPIBroadcast_character (physics->upw_choice,_MAX_STRING_SIZE_,0,&mpi->world);  CHECKERR(ierr);
#endif

  if (!strcmp(solver->SplitHyperbolicFlux,"yes")) {
    if (!mpi->rank) {
      fprintf(stderr,"Error in Euler2DInitialize: This physical model does not have a splitting ");
      fprintf(stderr,"of the hyperbolic term defined.\n");
    }
    return(1);
  }

  /* initializing physical model-specific functions */
  solver->ComputeCFL  = Euler2DComputeCFL;
  solver->FFunction   = Euler2DFlux;
  if      (!strcmp(physics->upw_choice,_ROE_ )) solver->Upwind = Euler2DUpwindRoe;
  else if (!strcmp(physics->upw_choice,_RF_  )) solver->Upwind = Euler2DUpwindRF;
  else if (!strcmp(physics->upw_choice,_LLF_ )) solver->Upwind = Euler2DUpwindLLF;
  else if (!strcmp(physics->upw_choice,_SWFS_)) solver->Upwind = Euler2DUpwindSWFS;
  else {
    fprintf(stderr,"Error in Euler2DInitialize(): %s is not a valid upwinding scheme.\n",
            physics->upw_choice);
    return(1);
  }
  solver->AveragingFunction     = Euler2DRoeAverage;
  solver->GetLeftEigenvectors   = Euler2DLeftEigenvectors;
  solver->GetRightEigenvectors  = Euler2DRightEigenvectors;

  /* set the value of gamma in all the boundary objects */
  int n;
  DomainBoundary  *boundary = (DomainBoundary*) solver->boundary;
  for (n = 0; n < solver->nBoundaryZones; n++)  boundary[n].gamma = physics->gamma;

  count++;
  return(0);
}

Euler2DCleanup()

int Euler2DCleanup

(

void *

)

Definition at line 4 of file Euler2DCleanup.c.

{
  return(0);
}