PLUTO Test Problems
4.4-patch2
|
Viscous compressible flow past a cylinder. More...
Functions | |
void | InitDomain (Data *d, Grid *grid) |
Set initial and boundary conditions for a flow past a cylinder in 2D cylindrical polar coordinates . The cylinder has radius 1 and the domain is initially filled with constant-density and pressure gas with value
.
The velocity field is initialized using the potential flow solution for an inviscid incompressible flow around a cylinder,
where U
, the far-field velocity, is given by the Mach number. The boundary conditions in phi
are periodic while the outer radial boundary is set to inflow for negative values of x while outflow for positive values. A no-slip boundary condition is used at the fluid-solid interface.
The flow past the cylinder, no matter how small the viscosity, will acquire vorticity in a thin boundary layer adjacent to the cylinder. Boundary layer separation may occur leading to the formation of a trailing wake behind the cylinder.
The input parameters are:
g_inputParam
[MACH]: sets the upstream Mach numberg_inputParam
[NU_VISC]: sets the viscosity