Hello everyone, I am trying to build a cross device of YIG like this in the figure (capture.png) image .
This is my script for my mumax3 code on this device
T = 220e-9
#dt = 1e-11
dt = 20e-12
dx = 500e-9
dy = 500e-9
Nx = 600
Ny = 900
A = 3.5e-12
Ms = 140e3
alpha = 0.0005
gamma = 1.758e11
sizeX = 300e-6
sizeY = 450e-6
sizeZ = 90e-9
script = f"""
//Specify output format
OutputFormat = OVF2_TEXT
//mesh
dx := {dx}
dy := {dy}
dz := 90e-9
Nx := {Nx}
Ny := {Ny}
Nz := 1
sizeX:= {sizeX}
sizeY:= {sizeY}
sizeZ:= {sizeZ}
a:=Cuboid(300e-6, 65e-6, {sizeZ})
//b:=Cuboid(40e-6, 450e-6, {sizeZ})
b:=Cuboid(40e-6, 192.5e-6, {sizeZ}).transl(0,128.70e-6, 0)
c:=Cuboid(40e-6, 192.5e-6, {sizeZ}).transl(0,-128.70e-6, 0)
d:=a.add(b).add(c)
SetMesh(Nx, Ny, Nz, dx, dy, dz, 0, 0, 0)
SetGeom(d)
Snapshot(geom)
EdgeSmooth=8
//External parameters
f := 1.73e9
Hac :=5e-5
Hdc := 0.018
Msat = {Ms}
Aex = {A}
alpha = {alpha}
tableadd(B_ext)
//Static field
B_ext = vector(Hdc, 0, 0)
defregion(2,xrange(-50e-6,-40e-6))
defregion(3,yrange(-97.5e-6, -87.5e-6))
defregion(4,yrange( 87.5e-6, 97.5e-6))
//Initial state
m=uniform(1,0,0)
//Solver
relax()
saveas(m, "relaxed")
//minimize()
//Output
B_ext.setregion(2, vector(Hdc, 0, Hac*sin(2*pi*f*t)))
B_ext.setregion(3, vector(Hdc, Hac*sin(2*pi*f*t ), 0))
//running
defregion(27,xrange(-Inf,-100e-6))
alpha.setregion(27,1)
defregion(28,xrange( 148e-6, Inf))
alpha.setregion(28,1)
defregion(17,yrange(223e-6,Inf))
alpha.setregion(17,1)
defregion(18,yrange(-Inf,-130e-6))
alpha.setregion(18,1)
autosave(m,200*{dt})
tableAdd(m.Region(5));
tableautosave({dt})
tableadd(B_ext)
run({T})
saveas(m, "final")
"""
Here along Horizontal axis BVMSW and along verticle axis MSSW is happening and different wavelength is observed for same signal and I have got the same figure for 0 degree phase and pi shifted signal the below one (output.png) . I mean proper constructive or destructive interference is not happening. What steps should I take to solve this issue and make a proper constructive and destructive interference?