Skip to content

Input File

Jump to bottom
weiqishen edited this page Apr 10, 2021 · 17 revisions

Input parameter definitions

The input file controls all the operations performed by HiFiLES. All available parameters are given below:

The general format of the parameter in the file is:

  1. For scalar input:

<parameter name> <parameter value>

  1. For vector input:

<parameter name> <number of values> <value 1> <value 2>....<value n>

Sample input file

Basic solver parameters

  • equation

    • 0 | Solve for Euler/Navier-Stokes equation
    • 1 | Solve for scalar advection/diffusion equation

  • viscous

    • 0 | inviscid(Euler)
    • 1 | viscous(Naiver-Stokes)

  • riemann_solve_type

    • 0 | Rusanov flux
    • 1 | Lax-Fredrich(advection/diffusion equation only)
    • 2 | RoeM flux
    • 3 | HLLC flux

  • vis_riemann_solve_type

    • 0 | LDG flux [default]

  • RANS

    • 0 | no RANS turbulence model [default]
    • 1 | S-A model (2D only)

  • ic_form

    • 0 | Isentropic vortex (Euler/N-S)
    • 1 | Uniform flow (Euler/N-S) [default]
    • 2 | Single sine wave (advection/diffusion equation)
    • 3 | Group sine wave (group) (advection/diffusion equation)
    • 4 | Spherical distribution (advection/diffusion equation)
    • 5 | uniform flow (advection/diffusion equation)
    • 6 | Up to 4th order polynomials for u, v, w (Euler/N-S)
    • 7 | Taylor-Green vortex (Euler/N-S)
    • 9 | Stationary shock (Euler/N-S)
    • 10| Shock tube (Euler/N-S)

  • x_shock_ic: x location of the shock for stationary shock and shock tube

  • test_case

    • 0 | do not run test case [default]
    • 1 | Isentropic vortex (Euler/N-S)
    • 2 | Single sine wave (advection/diffusion equation)
    • 3 | Group sine wave (group) (advection/diffusion equation)
    • 4 | Spherical wave (advection/diffusion equation)
    • 5 | Couette flow (Euler/N-S)

  • order: order of polynomial representation

  • dt_type

    • 0 | User defined time step (dt is needed)
    • 1 | Global time stepping (CFL is needed)
    • 2 | Local time stepping (CFL is needed)

  • n_steps: number of time steps to run

  • adv_type

    • 0 | Forward Euler
    • 1 | RK24
    • 2 | RK34
    • 3 | RK45
    • 4 | RK414

  • ldg_tau: tau for ldg flux [default 0]

  • ldg_beta: penalty factor for ldg flux [default 0.5]

Restart parameters

  • restart_flag

    • 0 | Start a new simulation [default]
    • 1 | Read ASCII restart file
    • 2 | Read HDF5 restart file

  • restart_iter: Iteration number to restart from

  • n_restart_files: number of restart file to read(only needed by ASCII restart file)

Mesh options

  • mesh_file: name of mesh file

  • dx_cyclic: Displacement (in x-direction) of periodic domain [default INFINITY]

  • dy_cyclic: Displacement (in y-direction) of periodic domain [default INFINITY]

  • dz_cyclic: Displacement (in z-direction) of periodic domain [default INFINITY]

LES parameters

  • LES

    • 0 | off [default]
    • 1 | on

  • filter_type(only needed when using SGS_model 234)

    • 0 | Vasilyev high-order commuting filter
    • 1 | Discrete Gaussian
    • 2 | Modal Vandermonde

  • filter_ratio: specify the ratio of filter width over element size. Enter any positive real.

  • SGS_model

    • 0 | Smagorinsky
    • 1 | Wall Adapting Local Eddy Viscosity (WALE)
    • 2 | WALE + Similarity Mixed (WSM)
    • 3 | Spectral Vanishing Viscosity (SVV): simply filters the solution
    • 4 | Bardina Similarity

  • wall_model Use with boundary condition option: use_wm

    • 0 | No wall model
    • 1 | Werner-Wengle model
    • 2 | Compressible log law for adiabatic wall

Monitoring, plotting parameters

  • p_res: number of nodes to plot on each edge of the element [default 2]

  • write_type

    • 0 | paraview
    • 1 | tecplot
    • 2 | CGNS

  • plot_freq: Frequency (in iterations) at which to output the solution [default INFINITY]

  • data_file_name: Prefix to use for all output files ["Mesh" by default]

  • restart_dump_freq: Frequency at which to output a restart file [default INFINITY]

  • monitor_res_freq: Frequency at which the residual is computed & output to the terminal [default 100]

  • calc_force

    • 0 | Don not calculate surface force
    • 1 | calculate surface force

  • area_ref: reference area to calculate surface force coefficients

  • monitor_cp_freq Frequency at which the Cp on wall boundaries is output to a file [default INFINITY]

  • res_norm_type

    • 0 | infinity norm
    • 1 | L1 norm
    • 2 | L2 norm [default]

  • error_norm_type

    • 0 | infinity norm
    • 1 | L1 norm
    • 2 | L2 norm [default]

  • diagnostic_fields

    • available fields: u v w energy mach pressure vorticity q_criterion(3D only) scaled_q_criterion(3D only) sensor

  • average_fields

    • available fields: rho_average u_average v_average w_average e_everage

  • integral_quantities

    • available fields: kineticenergy enstropy pressuredilatation straincolonproduct devstraincolonproduct

Element parameters

==== Tris ====

  • upts_type_tri

    • 0 | internal points (good quadrature points)
    • 1 | alpha optimized

  • fpts_type_tri

    • 0 | internal points (good quadrature points)
    • 1 | alpha optimized

  • vcjh_scheme_tri

    • 0 | use value of 'c' specified by user
    • 1 | 'DG'-like scheme (c = 0)
    • 2 | 'SD'-like scheme
    • 3 | Huyhn scheme
    • 4 | c+ scheme

  • c_tri: c for vcjh scheme

  • sparse_tri

    • 0 | Dense BLAS [default]
    • 1 | Sparse BLAS

=== Quads ====

  • upts_type_quad

    • 0 | Gauss
    • 1 | Gauss Lobatto

  • vcjh_scheme_quad

    • 0 | use value of 'c' specified by user
    • 1 | 'DG'-like scheme (c = 0)
    • 2 | 'SD'-like scheme
    • 3 | Huyhn scheme
    • 4 | c+ scheme

  • eta_quad: c for vcjh scheme

  • sparse_quad

    • 0 | Dense BLAS [default]
    • 1 | Sparse BLAS

==== Hexs ====

  • upts_type_hexa

    • 0 | Gauss
    • 1 | Gauss Lobatto

  • vcjh_scheme_hexa

    • 0 | use value of 'c' specified by user
    • 1 | 'DG'-like scheme (c = 0)
    • 2 | 'SD'-like scheme
    • 3 | Huyhn scheme
    • 4 | c+ scheme

  • eta_hexa: c for vcjh scheme

  • sparse_hexa

    • 0 | Dense BLAS [default]
    • 1 | Sparse BLAS

==== Tets ====

  • upts_type_tet

    • 0 | internal points (good quadrature points)
    • 1 | alpha optimized

  • vcjh_scheme_tet

    • 0 | use value of 'c' specified by user
    • 1 | 'DG'-like scheme (c = 0)
    • 2 | 'SD'-like scheme
    • 3 | Huyhn scheme
    • 4 | c+ scheme

  • eta_tet: c for vcjh scheme

  • sparse_tet

    • 0 | Dense BLAS [default]
    • 1 | Sparse BLAS

==== Prisms ====

  • upts_type_pri_tri

    • 0 | internal points (good quadrature points)
    • 1 | alpha optimized

  • upts_type_pri_1d

    • 0 | Gauss
    • 1 | Gauss Lobatto

  • vcjh_scheme_pri_1d

    • 0 | use value of 'c' specified by user
    • 1 | 'DG'-like scheme (c = 0)
    • 2 | 'SD'-like scheme
    • 3 | Huyhn scheme
    • 4 | c+ scheme

  • eta_pri: c for vcjh scheme

  • sparse_pri

    • 0 | Dense BLAS [default]
    • 1 | Sparse BLAS

Gas Parameters

  • gamma: ratio of specific heat [default 1.4]
  • prandtl: prandtl number [default 0.72]
  • S_gas: "S" from Sutherland's Law for viscosity [default 120.]
  • T_gas: temperature of the gas for Sutherland's Law in Kelvin [default 291.15]
  • R_gas: Ideal gas constant [default 286.9]
  • mu_gas: Dynamic viscosity (kg/m-s) [default 1.827E-05]
  • fix_vis
    • 0 | fixed viscosity
    • 1 | Sutherland's Law
  • prandtl_t: turbulent Prandtl number for turbulent models [default 0.9]

Boundary conditions

The boundary condition are in the following form

  • bc_[name of boundary]_type: type of the boundary in the mesh file with name [name of boundary]
  • bc_[name of boundary]_[param]: paramter needed by the boundary, can be found in source code
  • bc_[name of no-slip wall boundary]_use_wm: turn on wall model for this boundary

Supported boundary condition types are:

  • sub_in_simp: Subsonic inflow simple (free pressure)
  • sub_out_simp: Subsonic outflow simple (fixed pressure)
  • sub_in_char: Subsonic inflow characteristic
  • sub_out_char: Subsonic outflow characteristic
  • sup_in: Supersonic inflow
  • sup_out: Supersonic outflow
  • slip_wall: Slip wall
  • cyclic: Cyclic
  • isotherm_wall: Isothermal, no-slip wall
  • adiabat_wall: Adiabatic, no-slip wall
  • char: Characteristic
  • slip_wall_dual: Dual consistent BC
  • ad_wall: Advection, Advection-Diffusion Boundary Conditions

Parameters required by different types of boundary conditions can be found in the source file: src/input.cpp, function: void input::read_boundary_param(void)

Reference values

  • Mach_free_stream: reference Mach number to define refence velocity $U_{ref}$ [default 1.0]
  • L_free_stream: reference length $L_{ref}$ [default 1.0]
  • T_free_stream: reference temperature $T_{ref}$ [default 300.0]
  • rho_free_stream: reference densith $\rho_{ref}$ [default 1.17723946]

NOTE: These values are used to non-dimensionalize the equation for viscous simualtion only before the simualtion start.

Initial conditions

  1. For viscous flow, HiFiLES need:
  • Mach_c_ic: initial Mach number
  • nx_c_ic: x component of velocity direction vector [default 1]
  • ny_c_ic：y component of velocity direction vector [default 0]
  • nz_c_ic：z component of velocity direction vector [default 0]
  • T_c_ic: initial temperature
  • rho_c_ic: initial density
  1. For inviscid flow, HiFiLES need:
  • u_c_ic: initial x velocity
  • v_c_ic: initial y velocity
  • w_c_ic： iniital z velocity
  • p_c_ic: initial pressure
  • rho_c_ic: initial density

Solution patch

  • patch [default 0]
    • 0: do not use solution patch
    • 1: use solution patch
  • patch_type [default 0]
    • 0: vortex
    • 1: uniform flow patch

Parameter for composite vortex patch

  • Mv: Max Mach number of the vortex [default 0.5]
  • ra: inner radius for the vortex [default 0.075]
  • rb: outer radius for the vortex [default 0.175]
  • xc: x position of the vortex centroid [default 0.25]
  • yc: y postiton of the vortex centroid [default 0.5]

Parameters for uniform patch

  • patch_x: uniform patch coordinate x > patch_x with initial conditions

Shock capturing / dealiasing

  • over_int: over-integration de-aliasing mehtod [default 0]

    • 0: off
    • 1: on

  • over_int_order: order of over integration (recommend 1.5~2*Order)

  • shock_cap: shock capturing method [default 0]

    • 0: off
    • 1: exponential filter

  • shock_det: shock detector

    • 0: Persson's method
    • 1: concentration method

  • s0:threshold to apply artifitial dissipation

  • expf_fac: the strength of the expoenetial filter [default 36.0]

  • expf_order: the order of the exponential filter [default 4.0]

Advection-Diffusion Parameters

  • wave_speed_x: wave speed in x direction
  • wave_speed_y: wave speed in y direction [default 0]
  • wave_speed_z: wave speed in z direction [default 0]
  • diff_coeff: diffusion coefficient [default 0]
  • lambda: parameter for Lax-Fredrich flux

Uncategorized / Other

  • body_forcing [default 0]
  • perturb_ic [default 0]

Clone this wiki locally