MFC/src/common/m_constants.fpp at a2bf4ac82c2d80e47d7a641f4de5242609cbf278 · MFlowCode/MFC · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
!>
!! @file
!! @brief Contains constant values used throughout the code(s).

!> @brief Compile-time constant parameters: default values, tolerances, and physical constants
module m_constants

    use m_precision_select

    character, parameter :: dflt_char = ' ' !< Default string value

    real(wp), parameter :: dflt_real = -1.e6_wp                !< Default real value
    real(wp), parameter :: sgm_eps = 1.e-16_wp               !< Segmentation tolerance
    real(wp), parameter :: Chem_Tolerance = 1.e-16_wp               !< Speed of Sound Tolerance in Chemistry
    real(wp), parameter :: small_alf = 1.e-11_wp                !< Small alf tolerance
    real(wp), parameter :: pi = 3.141592653589793_wp !< Pi
    real(wp), parameter :: verysmall = 1.e-12_wp              !< Very small number

    integer, parameter :: num_stcls_min = 5                       !< Minimum # of stencils
    integer, parameter :: path_len = 400                          !< Maximum path length
    integer, parameter :: name_len = 50                           !< Maximum name length
    integer, parameter :: dflt_int = -100                         !< Default integer value
    integer, parameter :: fourier_rings = 5                       !< Fourier filter ring limit
    integer, parameter :: num_fluids_max = 10                     !< Maximum number of fluids in the simulation
    integer, parameter :: num_probes_max = 10                     !< Maximum number of flow probes in the simulation
    integer, parameter :: num_patches_max = 1000
    integer, parameter :: num_bc_patches_max = 10
    integer, parameter :: pathlen_max = 400
    integer, parameter :: nnode = 4    !< Number of QBMM nodes
    integer, parameter :: dflt_num_igr_iters = 2 !< number of iterations for IGR elliptic solve
    integer, parameter :: dflt_num_igr_warm_start_iters = 50 !< default number of iterations for IGR elliptic solve
    real(wp), parameter :: dflt_alf_factor = 10._wp !< scaling factor for IGR alpha
    integer, parameter :: gp_layers = 3 !< Number of ghost point layers for IBM
    real(wp), parameter :: capillary_cutoff = 1.e-6 !< color function gradient magnitude at which to apply the surface tension fluxes
    real(wp), parameter :: acoustic_spatial_support_width = 2.5_wp !< Spatial support width of acoustic source, used in s_source_spatial
    real(wp), parameter :: dflt_vcfl_dt = 100._wp !< value of vcfl_dt when viscosity is off for computing adaptive timestep size
    real(wp), parameter :: broadband_spectral_level_constant = 20._wp !< The constant to scale the spectral level at the lower frequency bound
    real(wp), parameter :: broadband_spectral_level_growth_rate = 10._wp !< The spectral level constant to correct the magnitude at each frequency to ensure the source is overall broadband

    ! Reconstruction Types
    integer, parameter :: WENO_TYPE = 1 !< Using WENO for reconstruction type
    integer, parameter :: MUSCL_TYPE = 2 !< Using MUSCL for reconstruction type

    ! Interface Compression
    real(wp), parameter :: dflt_ic_eps = 1e-4_wp !< Ensure compression is only applied to surface cells in THINC
    real(wp), parameter :: dflt_ic_beta = 1.6_wp !< Sharpness parameter's default value used in THINC
    integer, parameter :: moncon_cutoff = 1e-8_wp !< Monotonicity constraint's limiter to prevent extremas in THINC

    ! Chemistry
    real(wp), parameter :: dflt_T_guess = 1200._wp ! Default guess for temperature (when a previous value is not available)

    ! IBM+STL interpolation constants
    integer, parameter :: num_ray = 20 !< Default number of rays traced per cell
    real(wp), parameter :: ray_tracing_threshold = 0.9_wp !< Threshold above which the cell is marked as the model patch
    real(wp), parameter :: threshold_vector_zero = 1.e-10_wp !< Threshold to treat the component of a vector to be zero
    real(wp), parameter :: threshold_edge_zero = 1.e-10_wp !< Threshold to treat two edges to be overlapped
    real(wp), parameter :: initial_distance_buffer = 1.e12_wp !< Initialized levelset distance for the shortest path pair algorithm

    ! Lagrange bubbles constants
    integer, parameter :: mapCells = 3         !< Number of cells around the bubble/particle where the smoothening function will have effect
    real(wp), parameter :: R_uni = 8314._wp    !< Universal gas constant - J/kmol/K
    integer, parameter :: lag_io_vars = 21     !< Number of variables per particle for MPI_IO

    ! Strang Splitting constants
    real(wp), parameter :: dflt_adap_dt_tol = 1.e-4_wp !< Default tolerance for adaptive step size
    integer, parameter :: dflt_adap_dt_max_iters = 100 !< Default max iteration for adaptive step size

    ! Constants of the algorithm described by Heirer, E. Hairer, S. P.Nørsett, G. Wanner, Solving Ordinary Differential Equations I, Chapter II.4
    ! to choose the initial time step size for the adaptive time stepping routine
    real(wp), parameter :: threshold_first_guess = 1.e-5_wp
    real(wp), parameter :: threshold_second_guess = 1.e-15_wp
    real(wp), parameter :: scale_first_guess = 1.e-3_wp
    real(wp), parameter :: scale_guess = 1.e-2_wp
    real(wp), parameter :: small_guess = 1.e-6_wp

    ! Relativity
    integer, parameter :: relativity_cons_to_prim_max_iter = 100

    ! Pseudo-random number generator
    integer, parameter :: modulus = 2**30 - 1
    integer, parameter :: multiplier = 1664525
    integer, parameter :: increment = 1013904223
    integer, parameter :: amplifier = 3**13
    real(wp), parameter :: decimal_trim = 1.e5_wp

    ! System constants
    integer, parameter :: CASE_FILE_ERROR_CODE = 22

    ! Boundary condition enumeration
    ! Abbreviations
    !  CHAR - Characteristic
    !  NR   -  Non-reflecting
    !  SUB  - subsonic
    !  SUP  - supersonic
    !  FF   - Force-free
    !  CP   - Constant pressure
    integer, parameter :: BC_PERIODIC = -1
    integer, parameter :: BC_REFLECTIVE = -2
    integer, parameter :: BC_GHOST_EXTRAP = -3
    integer, parameter :: BC_RIEMANN_EXTRAP = -4
    integer, parameter :: BC_CHAR_SLIP_WALL = -5
    integer, parameter :: BC_CHAR_NR_SUB_BUFFER = -6
    integer, parameter :: BC_CHAR_NR_SUB_INFLOW = -7
    integer, parameter :: BC_CHAR_NR_SUB_OUTFLOW = -8
    integer, parameter :: BC_CHAR_FF_SUB_OUTFLOW = -9
    integer, parameter :: BC_CHAR_CP_SUB_OUTFLOW = -10
    integer, parameter :: BC_CHAR_SUP_INFLOW = -11
    integer, parameter :: BC_CHAR_SUP_OUTFLOW = -12
    integer, parameter :: BC_NULL = -13
    integer, parameter :: BC_AXIS = -14
    integer, parameter :: BC_SLIP_WALL = -15
    integer, parameter :: BC_NO_SLIP_WALL = -16
    integer, parameter :: BC_DIRICHLET = -17

end module m_constants