fluid {
  species = air_perfect
}

reference conditions {
  mach        = 1.44
  Re          = 2.3e6
  temperature = 292.81 # K
  alpha       = 0.
  beta        = 0.
  length      = 1.47702127659574E+00 # m  (diameter)
  area        = 1.71341823346392E+00 # m2 (frontal area)
}

time {
  steady      = false
  finish ctu  = 15
  initial cfl = 0.05
  cfl         = 2.
}

unst {

  explicit {
    method = ssprk33
  }

  # This will restart the solution from the specified restart file
  #restart {
  #  file = restart/Unst_restart.000005000.hdf5
  #}

  grid {
    file = voro_genesis_capsule.cgns
  }

  diagnostics {
    target size filter   = 0.007385
    max percent filtered = 0.01
  }

  convective terms {
    scheme = KEEPBLEND
    order = 2
    limiter = minmod
    flux function = ausmpw+
    add shock dissipation = true
    upwind blending sensors = [ducrosmod]
    upwind sensor enstrophy factor = 1
  }

  navierstokes {
    viscous = true
  }

  turbulence {
    sgs model = vreman
    sgs const = 0.08
  }

  wall function {
    type = blendedlog
  }

  output {

    # Screen and .iter file output interval
    monitor {
      interval  = 10
      variables = [nt, time, ctu, dt, maxcfl, maxmach, wallt, mups, cd_capsule, mups, elapsed]
    }

    volume {
      write initial = false
      format        = cgns
      interval      = 5000
      variables     = [q, density, mach, upwind]
    }

    # Outputs the solution at boudary condition regions only
    surface {
      interval  = 1000
      variables = [q, mach, tauw, pcoef, cf];
    }

    # Output restart files at specified interval
    restart {
      interval = 5000
    }

    # Output integrated forces and moments
    # multiple loads groups can be setup
    loads {
      moment center = [0., 0., 0.]
      interval      = 10
      # Provide non-dimensional force and moments
      normalize = true
      # Use gauge pressure (p-p_ref) in pressure forces
      gauge = true
      group_1 {
        name = capsule
        # Which BC zones to include
        include = [1-5]
      }
    }

    isosurfaces_1 {
      contour    = cutplane
      normal     = [0., 0., 1]
      point      = [0., 0., 0.]
      colorfield = fullq
      interval   = 500
    }

    isosurfaces_2 {
      contour    = cutplane
      normal     = [0., 0., 1]
      point      = [0., 0., 0.]
      colorfield = densitygradmag
      interval   = 500
    }

  }

  BC_1  {type = wallfunction}
  BC_2  {type = wallfunction}
  BC_3  {type = wallfunction}
  BC_4  {type = wallfunction}
  BC_5  {type = wallfunction}
  BC_6 {type  = farfield}
  BC_7 {type  = farfield}
  BC_8 {type  = farfield}
  BC_9 {type  = farfield}
  BC_10 {type = farfield}
  BC_11 {type = farfield}

}