NavierStokes3DFlux_GPU.cu File Reference

Functions to compute the hyperbolic flux for 3D Navier-Stokes equations. More...

#include <basic_gpu.h>
#include <arrayfunctions_gpu.h>
#include <mathfunctions.h>
#include <matmult_native.h>
#include <physicalmodels/navierstokes3d.h>
#include <hypar.h>

Go to the source code of this file.

Functions
__global__ void	gpuNavierStokes3DFlux_kernel (int npoints_grid, int dir, double gamma, const double *__restrict__ u, double *__restrict__ f)
int	gpuNavierStokes3DFlux (double *__restrict__ f, double *__restrict__ u, int dir, void *__restrict__ s, double t)

Detailed Description

Functions to compute the hyperbolic flux for 3D Navier-Stokes equations.

Author

Youngdae Kim

Definition in file NavierStokes3DFlux_GPU.cu.

Function Documentation

__global__ void gpuNavierStokes3DFlux_kernel	(	int	npoints_grid,
		int	dir,
		double	gamma,
		const double *__restrict__	u,
		double *__restrict__	f
	)

Kernel for gpuNavierStokes3DFlux()

Definition at line 69 of file NavierStokes3DFlux_GPU.cu.

{
  int p = blockDim.x * blockIdx.x + threadIdx.x;

  if (p < npoints_grid) {
    double rho, vx, vy, vz, e, P;
    _NavierStokes3DGetFlowVar_((u+p),npoints_grid,rho,vx,vy,vz,e,P,gamma);
    _NavierStokes3DSetFlux_((f+p),npoints_grid,rho,vx,vy,vz,e,P,dir);
  }
  return;
}

int gpuNavierStokes3DFlux	(	double *__restrict__	f,
		double *__restrict__	u,
		int	dir,
		void *__restrict__	s,
		double	t
	)

Compute the hyperbolic flux function for the 3D Navier-Stokes equations:

\begin{eqnarray} dir = x, & {\bf f}\left({\bf u}\right) = \left[\begin{array}{c} \rho u \\ \rho u^2 + p \\ \rho u v \\ \rho u w \\ (e+p)u \end{array}\right], \\ dir = y, & {\bf f}\left({\bf u}\right) = \left[\begin{array}{c} \rho v \\ \rho u v \\ \rho v^2 + p \\ \rho v w \\ (e+p)v \end{array}\right], \\ dir = z, & {\bf f}\left({\bf u}\right) = \left[\begin{array}{c} \rho w \\ \rho u w \\ \rho v w \\ \rho w^2 + p \\ (e+p)w \end{array}\right] \end{eqnarray}

Note: the flux function needs to be computed at the ghost points as well.

Parameters

f	Array to hold the computed flux vector (same layout as u)
u	Array with the solution vector
dir	Spatial dimension (x, y, or z) for which to compute the flux
s	Solver object of type HyPar
t	Current simulation time

Definition at line 96 of file NavierStokes3DFlux_GPU.cu.

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

  int nblocks = (solver->npoints_local_wghosts-1)/GPU_THREADS_PER_BLOCK + 1;

  gpuNavierStokes3DFlux_kernel<<<nblocks, GPU_THREADS_PER_BLOCK>>>(
    solver->npoints_local_wghosts, dir, param->gamma, u, f
  );
  cudaDeviceSynchronize();

  return 0;
}