MFC_PR/examples/2D_patch_modal_shape_exp/case.py at a53c1a84df101245ca5ff893a01e5f6c99145f21 · ChrisZYJ/MFC_PR · 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
116
117
118
119
120
121
122
123
124
125
126
127
#!/usr/bin/env python3
"""Minimal 2D acoustic case with geometry 13 in exponential form (modal_use_exp_form)."""
import math
import json

Nx, Ny = 64, 64
Lx = 8.0
dx = Lx / (Nx + 1)
p_inf = 101325.0
rho_inf = 1.0
gam = 1.4
# Patch 2: same pressure, higher sound speed (different gas) -> lower density
c_ratio = 1.5
rho2 = rho_inf / (c_ratio * c_ratio)
c = math.sqrt(gam * p_inf / rho_inf)
cfl = 0.3
mydt = cfl * dx / c
Tfinal = 10.0 / c
Nt = int(Tfinal / mydt)

# Acoustic source (similar to 2D_acoustic_support5)
t_cross = Lx / c  # domain crossing time
gauss_sigma_time = 0.03 * t_cross
acoustic_delay = 0.15 * t_cross
acoustic_loc_x = -3.5
acoustic_loc_y = 0.0
acoustic_aperture = 6.0
acoustic_foc_length = 3.5

print(
    json.dumps(
        {
            "run_time_info": "T",
            "x_domain%beg": -4.0,
            "x_domain%end": 4.0,
            "y_domain%beg": -4.0,
            "y_domain%end": 4.0,
            "m": Nx,
            "n": Ny,
            "p": 0,
            "dt": mydt,
            "t_step_start": 0,
            "t_step_stop": Nt,
            "t_step_save": max(1, int(Nt / 10)),
            "num_patches": 2,
            "model_eqns": 2,
            "alt_soundspeed": "F",
            "num_fluids": 2,
            "mpp_lim": "F",
            "mixture_err": "F",
            "time_stepper": 3,
            "weno_order": 5,
            "weno_eps": 1.0e-16,
            "mapped_weno": "T",
            "null_weights": "F",
            "mp_weno": "F",
            "riemann_solver": 2,
            "wave_speeds": 1,
            "avg_state": 2,
            "bc_x%beg": -8,
            "bc_x%end": -8,
            "bc_y%beg": -8,
            "bc_y%end": -8,
            "format": 1,
            "precision": 2,
            "prim_vars_wrt": "T",
            "parallel_io": "T",
            "omega_wrt(3)": "T",
            "fd_order": 2,
            "patch_icpp(1)%geometry": 3,
            "patch_icpp(1)%x_centroid": 0.0,
            "patch_icpp(1)%y_centroid": 0.0,
            "patch_icpp(1)%length_x": Lx,
            "patch_icpp(1)%length_y": Lx,
            "patch_icpp(1)%vel(1)": 0.0,
            "patch_icpp(1)%vel(2)": 0.0,
            "patch_icpp(1)%pres": p_inf,
            "patch_icpp(1)%alpha_rho(1)": rho_inf,
            "patch_icpp(1)%alpha_rho(2)": 0.0,
            "patch_icpp(1)%alpha(1)": 1.0,
            "patch_icpp(1)%alpha(2)": 0.0,
            "patch_icpp(2)%geometry": 13,
            "patch_icpp(2)%x_centroid": 0.0,
            "patch_icpp(2)%y_centroid": 0.0,
            "patch_icpp(2)%radius": 1.0,
            "patch_icpp(2)%modal_use_exp_form": "T",
            "patch_icpp(2)%fourier_cos(1)": -0.3,
            "patch_icpp(2)%fourier_sin(1)": 0.1,
            "patch_icpp(2)%fourier_sin(3)": -0.4,
            "patch_icpp(2)%fourier_cos(5)": 0.2,
            "patch_icpp(2)%fourier_sin(5)": 0.05,
            "patch_icpp(2)%vel(1)": 0.0,
            "patch_icpp(2)%vel(2)": 0.0,
            "patch_icpp(2)%pres": p_inf,
            "patch_icpp(2)%alpha_rho(1)": 0.0,
            "patch_icpp(2)%alpha_rho(2)": rho2,
            "patch_icpp(2)%alpha(1)": 0.0,
            "patch_icpp(2)%alpha(2)": 1.0,
            "patch_icpp(2)%alter_patch(1)": "T",
            "bc_x%grcbc_in": "F",
            "bc_x%grcbc_out": "T",
            "bc_x%grcbc_vel_out": "F",
            "bc_x%pres_out": p_inf,
            "bc_y%grcbc_in": "F",
            "bc_y%grcbc_out": "T",
            "bc_y%grcbc_vel_out": "F",
            "bc_y%pres_out": p_inf,
            "fluid_pp(1)%gamma": 1.0 / (gam - 1.0),
            "fluid_pp(1)%pi_inf": 0.0,
            "fluid_pp(2)%gamma": 1.0 / (gam - 1.0),
            "fluid_pp(2)%pi_inf": 0.0,
            # Acoustic source (similar to 2D_acoustic_support5)
            "acoustic_source": "T",
            "num_source": 1,
            "acoustic(1)%support": 5,
            "acoustic(1)%loc(1)": acoustic_loc_x,
            "acoustic(1)%loc(2)": acoustic_loc_y,
            "acoustic(1)%pulse": 2,
            "acoustic(1)%npulse": 1,
            "acoustic(1)%mag": 1.0,
            "acoustic(1)%gauss_sigma_time": gauss_sigma_time,
            "acoustic(1)%foc_length": acoustic_foc_length,
            "acoustic(1)%aperture": acoustic_aperture,
            "acoustic(1)%delay": acoustic_delay,
        }
    )
)