Hi,
I have a question in code mentioned in step - 81 tutorial.
I see the discretized form is written as (attached image).
Aij=∫Ω(μ−1r∇×φj)⋅(∇×φ¯i)dx−∫Ωεrφj⋅φ¯idx−i∫Σ(σΣrφjT)⋅φ¯iTdox−i∫∂Ω(μ−1rε−−−−√φ
jT)⋅(∇×φ¯iT)dox,
Fi=i∫ΩJa⋅φi¯dx−∫Ωμ−1rMa⋅(∇×φi¯)dx.
But in the code,
the magnetization term in the cell_rhs is absent.
This part of the rhs accounts for interior domain.
My question is: The reason of omitting that term is PML boundary condition?
or assumption Ma = 0?
I am not able to understand the flow here.
I would request some help with this.
const auto rhs_value =
(imag * scalar_product(J_a, phi_i)) *
fe_values.JxW(q_point);
cell_rhs(i) += rhs_value.real();
Thanks & Regards,
Pratyusa Kar
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion visit
https://groups.google.com/d/msgid/dealii/9391199d-bb98-4595-a11f-eefd0783c545n%40googlegroups.com.
