|
23 | 23 |
|
24 | 24 | dx = 1.0 / (1.0 * (Nx + 1)) |
25 | 25 |
|
26 | | -Tend = 64e-06 |
| 26 | +Tend = 1e-4 |
27 | 27 | Nt = 200 |
28 | 28 | mydt = Tend / (1.0 * Nt) |
29 | 29 |
|
|
48 | 48 | "t_step_stop": int(20 * (95 * size + 5)), |
49 | 49 | "t_step_save": int(20 * (95 * size + 5)), |
50 | 50 | # Simulation Algorithm Parameters |
51 | | - "num_patches": 2, |
| 51 | + "num_patches": 1, |
52 | 52 | "model_eqns": 2, |
53 | 53 | "alt_soundspeed": "F", |
54 | | - "num_fluids": 1, |
| 54 | + "num_fluids": 2, |
55 | 55 | "mpp_lim": "F", |
56 | 56 | "mixture_err": "F", |
57 | 57 | "time_stepper": 3, |
58 | 58 | "weno_order": 5, |
59 | 59 | "weno_eps": 1.0e-16, |
60 | | - "weno_Re_flux": "F", |
61 | | - "weno_avg": "F", |
62 | | - "mapped_weno": "F", |
| 60 | + "weno_Re_flux": "T", |
| 61 | + "weno_avg": "T", |
| 62 | + "mapped_weno": "T", |
63 | 63 | "null_weights": "F", |
64 | 64 | "mp_weno": "F", |
65 | | - "riemann_solver": 1, |
| 65 | + "riemann_solver": 2, |
66 | 66 | "wave_speeds": 1, |
67 | 67 | "avg_state": 2, |
68 | 68 | "bc_x%beg": -3, |
|
71 | 71 | "bc_y%end": -3, |
72 | 72 | "bc_z%beg": -3, |
73 | 73 | "bc_z%end": -3, |
74 | | - # Turning on Hypoelasticity |
75 | | - "hypoelasticity": "T", |
76 | | - "fd_order": 4, |
| 74 | + # Turn on IBM |
| 75 | + "ib": "T", |
| 76 | + "num_ibs": 1, |
| 77 | + "viscous": "T", |
77 | 78 | # Formatted Database Files Structure Parameters |
78 | 79 | "format": 1, |
79 | 80 | "precision": 2, |
80 | 81 | "prim_vars_wrt": "T", |
81 | 82 | "parallel_io": "T", |
82 | 83 | # Patch 1 L |
83 | 84 | "patch_icpp(1)%geometry": 9, |
84 | | - "patch_icpp(1)%x_centroid": 0.25, |
| 85 | + "patch_icpp(1)%x_centroid": 0.5, |
85 | 86 | "patch_icpp(1)%y_centroid": 0.25, |
86 | 87 | "patch_icpp(1)%z_centroid": 0.25, |
87 | | - "patch_icpp(1)%length_x": 0.5, |
| 88 | + "patch_icpp(1)%length_x": 1.0, |
88 | 89 | "patch_icpp(1)%length_y": 0.5, |
89 | 90 | "patch_icpp(1)%length_z": 0.5, |
90 | | - "patch_icpp(1)%vel(1)": 0.0, |
| 91 | + "patch_icpp(1)%vel(1)": 0.1, |
91 | 92 | "patch_icpp(1)%vel(2)": 0, |
92 | 93 | "patch_icpp(1)%vel(3)": 0, |
93 | | - "patch_icpp(1)%pres": 1.0e8, |
94 | | - "patch_icpp(1)%alpha_rho(1)": 1000, |
95 | | - "patch_icpp(1)%alpha(1)": 1.0, |
96 | | - "patch_icpp(1)%tau_e(1)": 0.0, |
97 | | - # Patch 2 R |
98 | | - "patch_icpp(2)%geometry": 9, |
99 | | - "patch_icpp(2)%x_centroid": 0.75, |
100 | | - "patch_icpp(2)%y_centroid": 0.25, |
101 | | - "patch_icpp(2)%z_centroid": 0.25, |
102 | | - "patch_icpp(2)%length_x": 0.5, |
103 | | - "patch_icpp(2)%length_y": 0.5, |
104 | | - "patch_icpp(2)%length_z": 0.5, |
105 | | - "patch_icpp(2)%vel(1)": 0, |
106 | | - "patch_icpp(2)%vel(2)": 0, |
107 | | - "patch_icpp(2)%vel(3)": 0, |
108 | | - "patch_icpp(2)%pres": 1.0e05, |
109 | | - "patch_icpp(2)%alpha_rho(1)": 1000, |
110 | | - "patch_icpp(2)%alpha(1)": 1.0, |
111 | | - "patch_icpp(2)%tau_e(1)": 0.0, |
| 94 | + "patch_icpp(1)%pres": 1.0, |
| 95 | + "patch_icpp(1)%alpha_rho(1)": 0.8e00, |
| 96 | + "patch_icpp(1)%alpha(1)": 0.8e00, |
| 97 | + "patch_icpp(1)%alpha_rho(2)": 0.2e00, |
| 98 | + "patch_icpp(1)%alpha(2)": 0.2e00, |
| 99 | + # Patch: Sphere Immersed Boundary |
| 100 | + "patch_ib(1)%geometry": 8, |
| 101 | + "patch_ib(1)%x_centroid": 0.25, |
| 102 | + "patch_ib(1)%y_centroid": 0.25, |
| 103 | + "patch_ib(1)%z_centroid": 0.25, |
| 104 | + "patch_ib(1)%radius": 0.1, |
112 | 105 | # Fluids Physical Parameters |
113 | | - "fluid_pp(1)%gamma": 1.0e00 / (4.4e00 - 1.0e00), |
114 | | - "fluid_pp(1)%pi_inf": 4.4e00 * 6.0e08 / (4.4e00 - 1.0e00), |
115 | | - "fluid_pp(1)%G": 10e09, |
| 106 | + # Specify 2 fluids |
| 107 | + "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1.0e00), |
| 108 | + "fluid_pp(1)%pi_inf": 0, |
| 109 | + "fluid_pp(1)%Re(1)": 54000, |
| 110 | + "fluid_pp(2)%gamma": 1.0e00 / (1.4 - 1.0e00), |
| 111 | + "fluid_pp(2)%pi_inf": 0, |
| 112 | + "fluid_pp(2)%Re(1)": 54000, |
116 | 113 | } |
117 | 114 | ) |
118 | 115 | ) |
0 commit comments