PLUTO Test Problems  4.4-patch2
Functions
init.c File Reference

Viscous compressible flow past a cylinder. More...

Functions

void InitDomain (Data *d, Grid *grid)
 

Detailed Description

Set initial and boundary conditions for a flow past a cylinder in 2D cylindrical polar coordinates $(r,\phi)$. The cylinder has radius 1 and the domain is initially filled with constant-density and pressure gas with value $\rho = 1,\, p = 1/\Gamma$.

The velocity field is initialized using the potential flow solution for an inviscid incompressible flow around a cylinder,

\[ V_r = U\left(1 - \frac{1}{r^2}\right)\cos\theta \,,\qquad V_\phi = -U\left(1 + \frac{1}{r^2}\right)\sin\theta \]

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:

hd_flow_past_cylinder.01.jpg
Density map for configuration #01
hd_flow_past_cylinder.04.jpg
Entropy distribution for configuration #04 using 3 levels of refinement.
Author
A. Mignone (migno.nosp@m.ne@t.nosp@m.o.inf.nosp@m.n.it)
Date
Feb 25, 2019

Function Documentation

◆ InitDomain()

void InitDomain ( Data *  d,
Grid *  grid 
)

Assign initial condition by looping over the computational domain. Called after the usual Init() function to assign initial conditions on primitive variables. Value assigned here will overwrite those prescribed during Init().