How to analytically (and numerically) define a spatially varying deltaH for a ferromagnet in HFSS?
I'm using Ansys HFSS 2019 R1 and trying to do magnetostatic wave simulations with a ferromagnetic film (this is a microwave simulation, so Maxwell 3D won't work). The problem is that I need to manually create an absorbing boundary condition (ABC) for these waves (radiation and PML don't work for my structure and these types of waves). I can only model a small extent of this ferromagnetic film due to mesh and computational constraints (magnetostatic wavelengths are typically five orders of magnitude less than EM waves at a given frequency).
To create an ABC, I'd like to define the spatially varying magnetic line width (deltaH), preferably with an analytical expression. Defining it numerically (and choosing the axis along which it varies) would work too. I know how to reference intrinsic variables like X, Y, and Z, but when I define the material with an analytical expression, the results don't seem to make sense. For an exponentially graded ABC (more than a wavelength at the lowest frequency), the results seem to yield significant reflections, so I'm not convinced HFSS is interpreting my expression as I intend. I would also like to plot this expression as a function of space so I see exactly what HFSS is receiving as input.
If I define the spatially varying line width numerically, I'd like to select along which axis, and this is not clear from the help documentation. I do not want to have to do a coordinate transformation, since this complicates defining the magnetic bias field.
I create two project variables, $aa and $continuous_alpha:
$aa = 2*ln($alpha_max/$alpha_min)/($SubsY)
$continuous_alpha = $alpha_min*exp(abs($aa*Y))
Then, under View/Edit Material, I enter in the following expression which converts alpha (the damping factor) to the linewidth (deltaH):
Delta H = -($continuous_alpha)*4*pi*5.714e9/(4*pi*1e-7*$gyro)
I have no units in my expression except for the intrinsic variable Y, which is defined with mm units. $SubsY is 0.2 (without units, but is the number if it were defined in mm units). $alpha_max is 0.01 and $alpha_min is 1e-5. So in 100 um I grade the damping constant alpha 3 orders of magnitude. Below is a picture of the setup and of a linecut along the y-axis of the wave I get.
I've verified this expression in Matlab. This is one specific example, but I've tried various other alpha values, substrate widths, and CPW dimensions and it's hard to consistently get a good ABC. I'd first like to convince myself that HFSS is receiving the profile I think it is, but I have no idea how to do that.
Thanks for your help!
P.S. I have made my ABC using discrete layers (40), and that worked acceptably well for one CPW size but not others. For other sizes, I still see some reflections. When I simulate over various frequencies, this shows up as a frequency shift in the s-parameters. So I believe my best bet is to specify a continuous ABC and refine the mesh as needed.