Dear everyone
I have writen 3D PML (based on Barenger) code in matlab but I have confused when I joint between pml and main together at interface region , So I have changed to write in UPML . I would like to ask for suggestion in this method ,I am not sure that in main region also must use UPML formula too or not, or use general formula(in free space)?
Could anyone summarize the method of 3D UPML ?
I summarize this PML in Chapter 4 of my dissertation. You can certainly use the PML equations in your problem space as long as your "lossy" parameters are zero. The lossy terms should slowly increase as you progress into any of your absorbing boundaries. Since implementing the PML involves integrations and hence extra calculations, it is most efficient to only perform these computations in the PML sections.
It's my 3d upml code Reference by formula in Taflove's book.
You should have "Computational Electrodynamics FDTD" it's Taflove's book.
It'll be helpful for do your project and write a code.
Hope this help
regards
CPML is better than UPML because CPML can treats dispersive case. Moreover, the codes are also clearly stated in Taflove's book. Please refer to CHAPTER 7