MFC/examples/2D_IGR_forward_facing_step/case.py at 1b28f7622ce2618005481dac66b76cc8462d4161 · 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
import json
import math

h = 0.2

# Radius as a percentage of height (h)
rc = 0.2

gam_a = 1.4
p0 = 1
rho0 = 1.4
c0 = math.sqrt(gam_a * p0 / rho0)
v0 = 3 * c0
mu = rho0 * v0 * h / 2e5

# Configuring case dictionary
print(
    json.dumps(
        {
            # Logistics
            "run_time_info": "T",
            "x_domain%beg": 0,
            "x_domain%end": 15 * h,
            "y_domain%beg": 0,
            "y_domain%end": 5 * h,
            "cyl_coord": "F",
            "m": 1499,
            "n": 499,
            "p": 0,
            "cfl_adap_dt": "T",
            "cfl_target": 0.6,
            "n_start": 0,
            "t_save": 0.04,
            "t_stop": 4,
            # Simulation Algorithm Parameters
            "num_patches": 1,
            "model_eqns": 2,
            "alt_soundspeed": "F",
            "num_fluids": 1,
            "mpp_lim": "F",
            "mixture_err": "F",
            "time_stepper": 3,
            "igr": "T",
            "igr_pres_lim": "T",
            "igr_order": 3,
            "igr_iter_solver": 1,
            "num_igr_iters": 5,
            "num_igr_warm_start_iters": 50,
            "bc_x%beg": -3,
            "bc_x%end": -3,
            "bc_y%beg": -2,
            "bc_y%end": -2,
            "ib": "T",
            "num_ibs": 3,
            # Formatted Database Files Structure Parameters
            "format": 1,
            "precision": 2,
            "prim_vars_wrt": "T",
            "parallel_io": "T",
            # Patch 1 Background
            "patch_icpp(1)%geometry": 3,
            "patch_icpp(1)%x_centroid": 7.5 * h,
            "patch_icpp(1)%y_centroid": 2.5 * h,
            "patch_icpp(1)%length_x": 15 * h,
            "patch_icpp(1)%length_y": 5 * h,
            "patch_icpp(1)%vel(1)": v0,
            "patch_icpp(1)%vel(2)": 0.0,
            "patch_icpp(1)%pres": p0,
            "patch_icpp(1)%alpha_rho(1)": rho0,
            "patch_icpp(1)%alpha(1)": 1.0,
            # Patch: Slip rectangle with rounded corner
            "patch_ib(1)%geometry": 3,
            "patch_ib(1)%x_centroid": (9 + rc / 2) * h,
            "patch_ib(1)%y_centroid": 0.5 * h,
            "patch_ib(1)%length_x": (12 - rc) * h,
            "patch_ib(1)%length_y": h,
            "patch_ib(1)%slip": "T",
            "patch_ib(2)%geometry": 3,
            "patch_ib(2)%x_centroid": (3 + rc / 2) * h,
            "patch_ib(2)%y_centroid": (0.5 - rc / 2) * h,
            "patch_ib(2)%length_x": rc * h,
            "patch_ib(2)%length_y": (1 - rc) * h,
            "patch_ib(2)%slip": "T",
            "patch_ib(3)%geometry": 2,
            "patch_ib(3)%x_centroid": (3 + rc) * h,
            "patch_ib(3)%y_centroid": (1 - rc) * h,
            "patch_ib(3)%radius": rc * h,
            "patch_ib(3)%slip": "T",
            # Fluids Physical Parameters
            "fluid_pp(1)%gamma": 1.0 / (gam_a - 1.0),
            "fluid_pp(1)%pi_inf": 0.0,
            "viscous": "T",
            "fluid_pp(1)%Re(1)": 1 / mu,
        },
        indent=4,
    )
)