Skip to content

Commit f47784b

Browse files
ajnonakaajnonaka
committed
fix PML for both PEC and conductivity cases for inductor test
1 parent cb74d48 commit f47784b

2 files changed

Lines changed: 17 additions & 27 deletions

File tree

‎Examples/Tests/circuits/Candice/inputs_inductor_on_cap‎

Lines changed: 6 additions & 11 deletions
Original file line numberDiff line numberDiff line change
@@ -3,7 +3,7 @@
33
################################
44
####### GENERAL PARAMETERS ######
55
#################################
6-
max_step = 10000
6+
max_step = 200000
77

88
amr.n_cell = n_cellx n_celly n_cellz
99
amr.max_grid_size = max_grid_sizex max_grid_sizey max_grid_sizez
@@ -37,14 +37,14 @@ algo.macroscopic_sigma_method = laxwendroff
3737

3838
algo.lumped_inductor = on
3939

40-
inductor.inductor_x_function(x,y,z) = "4.e-10 * (x > -460.e-6) * (x < -395.e-6) * (y > -ddy) * (y < ddy) * (z > h_si - ddz) * (z < h_si + ddz)"
40+
inductor.inductor_x_function(x,y,z) = "4.e-13 * (x > -460.e-6) * (x < -395.e-6) * (y > -ddy) * (y < ddy) * (z > h_si - ddz) * (z < h_si + ddz)"
4141
inductor.inductor_y_function(x,y,z) = "0."
4242
inductor.inductor_z_function(x,y,z) = "0."
4343

4444
macroscopic.sigma_function(x,y,z) = "sigma_0 + (sigma_si - sigma_0) * (z <= h_si) + (sigma_metal - sigma_0) * (z >= h_si) * (z <= h_si + dz) *
4545
(1
46-
-(x > -1750.e-6 + ddx) * (x < -795.e-6 - ddx) * (y > -32.5e-6 + ddy) * (y < -27.5e-6 - ddy)
47-
-(x > -1750.e-6 + ddx) * (x < -795.e-6 - ddx) * (y > -20.e-6 + ddy) * (y < -15.e-6 - ddy)
46+
-(x < -795.e-6 - ddx) * (y > -32.5e-6 + ddy) * (y < -27.5e-6 - ddy)
47+
-(x < -795.e-6 - ddx) * (y > -20.e-6 + ddy) * (y < -15.e-6 - ddy)
4848
-(x > -800.e-6 + ddx) * (x < -795.e-6 - ddx) * (y > -27.5e-6 - ddy) * (y < -20.e-6 + ddy)
4949
-(x > -695.e-6 + ddx) * (x < -160.e-6 - ddx) * (y > -372.5e-6 + ddy) * (y < -272.5e-6 - ddy)
5050
-(x > -695.e-6 + ddx) * (x < -160.e-6 - ddx) * (y > 272.5e-6 + ddy) * (y < 372.5e-6 - ddy)
@@ -55,11 +55,6 @@ macroscopic.sigma_function(x,y,z) = "sigma_0 + (sigma_si - sigma_0) * (z <= h_si
5555

5656
macroscopic.epsilon_function(x,y,z) = "eps_0 + eps_0 * (eps_r_si - 1.) * (z <= h_si)"
5757

58-
# lines 1-4: qubit control
59-
# lines 5-9: readout resonator
60-
# lines 9-13: capacitor
61-
# lines 14-18: entrance to meanding line
62-
# lines 19-68: meandering line
6358
macroscopic.mu_function(x,y,z) = "mu_0 + mu_0 * (mu_r_si - 1.) * (z <= h_si)"
6459

6560
#################################
@@ -122,7 +117,7 @@ my_constants.w_port = 320.e-6
122117

123118
my_constants.pi = 3.14159265358979
124119

125-
my_constants.freq = 8.6e9
120+
my_constants.freq = 1.9e11
126121
my_constants.TP = 1./freq
127122

128123
# grid spacing
@@ -162,7 +157,7 @@ warpx.Ez_excitation_grid_function(x,y,z,t) = "0."
162157
diagnostics.diags_names = plt
163158
###############
164159
# full plotfiles
165-
plt.intervals = 100
160+
plt.intervals = 1000
166161
plt.fields_to_plot = Ex Ey Ez Bx By Bz sigma
167162
plt.diag_type = Full
168163
plt.file_min_digits = 7

‎Examples/Tests/circuits/Candice/inputs_inductor_on_cap_PEC‎

Lines changed: 11 additions & 16 deletions
Original file line numberDiff line numberDiff line change
@@ -3,7 +3,7 @@
33
################################
44
####### GENERAL PARAMETERS ######
55
#################################
6-
max_step = 10000
6+
max_step = 200000
77

88
amr.n_cell = n_cellx n_celly n_cellz
99
amr.max_grid_size = max_grid_sizex max_grid_sizey max_grid_sizez
@@ -37,7 +37,7 @@ algo.macroscopic_sigma_method = laxwendroff
3737

3838
algo.lumped_inductor = on
3939

40-
inductor.inductor_x_function(x,y,z) = "4.e-10 * (x > -460.e-6) * (x < -395.e-6) * (y > -ddy) * (y < ddy) * (z > h_si - ddz) * (z < h_si + ddz)"
40+
inductor.inductor_x_function(x,y,z) = "4.e-13 * (x > -460.e-6) * (x < -395.e-6) * (y > -ddy) * (y < ddy) * (z > h_si - ddz) * (z < h_si + ddz)"
4141
inductor.inductor_y_function(x,y,z) = "0."
4242
inductor.inductor_z_function(x,y,z) = "0."
4343

@@ -50,15 +50,10 @@ macroscopic.sigma_function(x,y,z) = "sigma_0 + (sigma_si - sigma_0) * (z <= h_si
5050

5151
macroscopic.epsilon_function(x,y,z) = "eps_0 + eps_0 * (eps_r_si - 1.) * (z <= h_si)"
5252

53-
# lines 1-4: qubit control
54-
# lines 5-9: readout resonator
55-
# lines 9-13: capacitor
56-
# lines 14-18: entrance to meanding line
57-
# lines 19-68: meandering line
5853
macroscopic.mu_function(x,y,z) = "mu_0 + mu_0 * (mu_r_test - 1.) * (z > h_si) * (z < h_si + dz) *
5954
(1
60-
-(x > -1750.e-6 + ddx) * (x < -795.e-6 - ddx) * (y > -32.5e-6 + ddy) * (y < -27.5e-6 - ddy)
61-
-(x > -1750.e-6 + ddx) * (x < -795.e-6 - ddx) * (y > -20.e-6 + ddy) * (y < -15.e-6 - ddy)
55+
-(x > -1750.e-6 - ddx) * (x < -795.e-6 - ddx) * (y > -32.5e-6 + ddy) * (y < -27.5e-6 - ddy)
56+
-(x > -1750.e-6 - ddx) * (x < -795.e-6 - ddx) * (y > -20.e-6 + ddy) * (y < -15.e-6 - ddy)
6257
-(x > -800.e-6 + ddx) * (x < -795.e-6 - ddx) * (y > -27.5e-6 - ddy) * (y < -20.e-6 + ddy)
6358
-(x > -695.e-6 + ddx) * (x < -160.e-6 - ddx) * (y > -372.5e-6 + ddy) * (y < -272.5e-6 - ddy)
6459
-(x > -695.e-6 + ddx) * (x < -160.e-6 - ddx) * (y > 272.5e-6 + ddy) * (y < 372.5e-6 - ddy)
@@ -127,7 +122,7 @@ my_constants.w_port = 320.e-6
127122

128123
my_constants.pi = 3.14159265358979
129124

130-
my_constants.freq = 8.6e9
125+
my_constants.freq = 1.9e11
131126
my_constants.TP = 1./freq
132127

133128
# grid spacing
@@ -151,8 +146,8 @@ warpx.E_excitation_on_grid_style = parse_E_excitation_grid_function
151146

152147
warpx.Ex_excitation_flag_function(x,y,z) = "flag_hs * (z > h_si - ddz) * (z < h_si + ddz) *
153148
(1
154-
-(x > -1750.e-6 + ddx) * (x < -795.e-6 - ddx) * (y > -32.5e-6 + ddy) * (y < -27.5e-6 - ddy)
155-
-(x > -1750.e-6 + ddx) * (x < -795.e-6 - ddx) * (y > -20.e-6 + ddy) * (y < -15.e-6 - ddy)
149+
-(x > -1750.e-6 - ddx) * (x < -795.e-6 - ddx) * (y > -32.5e-6 + ddy) * (y < -27.5e-6 - ddy)
150+
-(x > -1750.e-6 - ddx) * (x < -795.e-6 - ddx) * (y > -20.e-6 + ddy) * (y < -15.e-6 - ddy)
156151
-(x > -800.e-6 + ddx) * (x < -795.e-6 - ddx) * (y > -27.5e-6 - ddy) * (y < -20.e-6 + ddy)
157152
-(x > -695.e-6 + ddx) * (x < -160.e-6 - ddx) * (y > -372.5e-6 + ddy) * (y < -272.5e-6 - ddy)
158153
-(x > -695.e-6 + ddx) * (x < -160.e-6 - ddx) * (y > 272.5e-6 + ddy) * (y < 372.5e-6 - ddy)
@@ -163,8 +158,8 @@ warpx.Ex_excitation_flag_function(x,y,z) = "flag_hs * (z > h_si - ddz) * (z < h_
163158

164159
warpx.Ey_excitation_flag_function(x,y,z) = "flag_hs * (z > h_si - ddz) * (z < h_si + ddz) *
165160
(1
166-
-(x > -1750.e-6 + ddx) * (x < -795.e-6 - ddx) * (y > -32.5e-6 + ddy) * (y < -27.5e-6 - ddy)
167-
-(x > -1750.e-6 + ddx) * (x < -795.e-6 - ddx) * (y > -20.e-6 + ddy) * (y < -15.e-6 - ddy)
161+
-(x > -1750.e-6 - ddx) * (x < -795.e-6 - ddx) * (y > -32.5e-6 + ddy) * (y < -27.5e-6 - ddy)
162+
-(x > -1750.e-6 - ddx) * (x < -795.e-6 - ddx) * (y > -20.e-6 + ddy) * (y < -15.e-6 - ddy)
168163
-(x > -800.e-6 + ddx) * (x < -795.e-6 - ddx) * (y > -27.5e-6 - ddy) * (y < -20.e-6 + ddy)
169164
-(x > -695.e-6 + ddx) * (x < -160.e-6 - ddx) * (y > -372.5e-6 + ddy) * (y < -272.5e-6 - ddy)
170165
-(x > -695.e-6 + ddx) * (x < -160.e-6 - ddx) * (y > 272.5e-6 + ddy) * (y < 372.5e-6 - ddy)
@@ -189,7 +184,7 @@ warpx.Ez_excitation_grid_function(x,y,z,t) = "0."
189184
diagnostics.diags_names = plt
190185
###############
191186
# full plotfiles
192-
plt.intervals = 100
193-
plt.fields_to_plot = Ex Ey Ez Bx By Bz sigma
187+
plt.intervals = 1000
188+
plt.fields_to_plot = Ex Ey Ez Bx By Bz mu
194189
plt.diag_type = Full
195190
plt.file_min_digits = 7

0 commit comments

Comments
 (0)