NavierStokes3DFunctions.c File Reference

Miscellaneous functions for the 3D Navier Stokes equations. More...

#include <stdio.h>
#include <math.h>
#include <basic.h>
#include <arrayfunctions.h>
#include <physicalmodels/navierstokes3d.h>
#include <hypar.h>

Go to the source code of this file.

Functions
int	NavierStokes3DRoeAverage (double *uavg, double *uL, double *uR, void *p)
int	NavierStokes3DComputePressure (double *P, const double *const u, void *s)
int	NavierStokes3DComputeTemperature (double *T, const double *const u, void *s)

Detailed Description

Miscellaneous functions for the 3D Navier Stokes equations.

Author

Debojyoti Ghosh

Definition in file NavierStokes3DFunctions.c.

Function Documentation

NavierStokes3DRoeAverage()

int NavierStokes3DRoeAverage	(	double *	uavg,
		double *	uL,
		double *	uR,
		void *	p
	)

Compute the Roe-averaged state for the 3D Navier Stokes equations. This function just calls the macro _NavierStokes3DRoeAverage_ and is not used by any functions within the 3D Navier Stokes module. However, it's necessary to define it and provide it to the the solver object (HyPar) so that it can then send it to interpolation functions for a characteristic-based reconstruction.

Parameters

uavg	The computed Roe-averaged state
uL	Left state (conserved variables)
uR	Right state (conserved variables)
p	Object of type NavierStokes3D with physics-related variables

Definition at line 18 of file NavierStokes3DFunctions.c.

{
  NavierStokes3D *param  = (NavierStokes3D*) p;
  _NavierStokes3DRoeAverage_(uavg,_NavierStokes3D_stride_,uL,uR,param->gamma);
  return(0);
}

NavierStokes3DComputePressure()

int NavierStokes3DComputePressure	(	double *	P,
		const double *const	u,
		void *	s
	)

Compute the pressure from the conserved solution on a grid

Parameters

P	Array to hold the computed pressure (same layout as u)
u	Array with the solution vector
s	Solver object of type HyPar

Definition at line 31 of file NavierStokes3DFunctions.c.

{
  HyPar          *solver = (HyPar*)   s;
  NavierStokes3D *param  = (NavierStokes3D*) solver->physics;
  int            i;

  int *dim    = solver->dim_local;
  int ghosts  = solver->ghosts;
  int ndims   = solver->ndims;
  int index[ndims], bounds[ndims], offset[ndims];

  /* set bounds for array index to include ghost points */
  _ArrayCopy1D_(dim,bounds,ndims);
  for (i=0; i<ndims; i++) bounds[i] += 2*ghosts;

  /* set offset such that index is compatible with ghost point arrangement */
  _ArraySetValue_(offset,ndims,-ghosts);

  int done = 0; _ArraySetValue_(index,ndims,0);
  while (!done) {
    int idx; _ArrayIndex1DWO_(ndims,dim,index,offset,ghosts,idx);
    double rho, vx, vy, vz, e, pressure;
    _NavierStokes3DGetFlowVar_((u+_MODEL_NVARS_*idx),_NavierStokes3D_stride_,rho,vx,vy,vz,e,pressure,param->gamma);
    P[idx] = pressure;
    _ArrayIncrementIndex_(ndims,bounds,index,done);
  }
  
  return(0);
}

NavierStokes3DComputeTemperature()

int NavierStokes3DComputeTemperature	(	double *	T,
		const double *const	u,
		void *	s
	)

Compute the temperature from the conserved solution on a grid

Parameters

T	Array to hold the computed pressure (same layout as u)
u	Array with the solution vector
s	Solver object of type HyPar

Definition at line 65 of file NavierStokes3DFunctions.c.

{
  HyPar          *solver = (HyPar*)   s;
  NavierStokes3D *param  = (NavierStokes3D*) solver->physics;
  int            i;

  int *dim    = solver->dim_local;
  int ghosts  = solver->ghosts;
  int ndims   = solver->ndims;
  int index[ndims], bounds[ndims], offset[ndims];

  /* set bounds for array index to include ghost points */
  _ArrayCopy1D_(dim,bounds,ndims);
  for (i=0; i<ndims; i++) bounds[i] += 2*ghosts;

  /* set offset such that index is compatible with ghost point arrangement */
  _ArraySetValue_(offset,ndims,-ghosts);

  int done = 0; _ArraySetValue_(index,ndims,0);
  while (!done) {
    int idx; _ArrayIndex1DWO_(ndims,dim,index,offset,ghosts,idx);
    double rho, vx, vy, vz, e, pressure;
    _NavierStokes3DGetFlowVar_((u+_MODEL_NVARS_*idx),_NavierStokes3D_stride_,rho,vx,vy,vz,e,pressure,param->gamma);
    T[idx] = pressure/rho;
    _ArrayIncrementIndex_(ndims,bounds,index,done);
  }
  
  return(0);
}