Newer
Older
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
-- Resonant cavity test. See Pozar, section 6.3 for theory.
-----------------------
--[[ Problem setup ]]--
sim.name = "test_maxwell_resonator" -- simulation name
sim.dt = 1e-12 -- timestep size
sim.timesteps = 10001 -- num of iterations
sim.gmsh_model = "resonator.geo" -- gmsh model filename
sim.use_gpu = 1 -- 0: cpu, 1: gpu
sim.approx_order = 1 -- approximation order
sim.time_integrator = "leapfrog"
postpro.silo_output_rate = 100
postpro.cycle_print_rate = 100 -- console print rate
postpro["J"].silo_mode = "none"
local epsr = 1
local mur = 1
materials[1] = {}
materials[1].epsilon = epsr
materials[1].mu = mur
materials[1].sigma = 0
function electric_initial_condition(x, y, z)
local Ex = 0
local Ey = math.sin(math.pi*x) * math.sin(math.pi*z)
local Ez = 0
return Ex, Ey, Ez
end
--------------------------
--[[ Validation stuff ]]--
debug = {}
local c0 = 1/math.sqrt(const.eps0*const.mu0)
-- Mode
local m = 1 -- along x
local n = 0 -- along y
local l = 1 -- along z
-- Cavity dimensions (must match sim.gmsh_model)
local a = 1 -- along x
local b = 0.1 -- along y
local d = 1 -- along z
local u = m*math.pi/a
local v = n*math.pi/b
local w = l*math.pi/d
-- Compute resonant frequency
local omega0 = c0*math.sqrt(u*u + v*v + w*w)/math.sqrt(epsr*mur)
-- Compute impedance
local eps = materials[1].epsilon * const.eps0
local mu = materials[1].mu * const.mu0
local Y = math.sqrt(eps/mu)
function ansol(tag, x, y, z, t)
local Ex = 0.0
local Ey = math.cos(omega0*t)*math.sin(math.pi*x)*math.sin(math.pi*z)
local Ez = 0.0
local Hx = -Y*math.sin(omega0*t)*math.sin(math.pi*x)*math.cos(math.pi*z)
local Hy = 0.0
local Hz = Y*math.sin(omega0*t)*math.cos(math.pi*x)*math.sin(math.pi*z)
return Ex, Ey, Ez, Hx, Hy, Hz
end
--debug.analytical_solution = ansol