Electronics World articles Popular Electronics articles QST articles Radio & TV News articles Radio-Craft articles Radio-Electronics articles Short Wave Craft articles Wireless World articles Google Search of RF Cafe website Sitemap Electronics Equations Mathematics Equations Equations physics Manufacturers & distributors Engineer Jobs LinkedIn Crosswords Engineering Humor Kirt's Cogitations Engineering Event Calendar RF Engineering Quizzes USAF radar shop Notable Quotes App Notes Calculators Education Engineering Magazines Engineering magazine articles Engineering software Engineering smorgasbord RF Cafe Archives RF Cascade Workbook 2018 RF Stencils for Visio RF & EE Shapes for Word Advertising RF Cafe Homepage Sudoku puzzles Thank you for visiting RF Cafe!
Exodus Advanced Communications

Trouble with UPML in FDTD - RF Cafe Forums

The original RF Cafe Forums were shut down in late 2012 due to maintenance issues - primarily having to spend time purging garbage posts from the board. At some point I might start the RF Cafe Forums again if the phpBB software gets better at filtering spam.

Below are the old forum threads, including responses to the original posts.

-- Amateur Radio
-- Anecdotes, Gripes & Humor
-- Antennas
-- CAE, CAD, & Software
-- Circuits & Components
-- Employment & Interviews
-- Miscellany
-- Swap Shop
-- Systems
-- Test & Measurement
-- Webmaster

metal
 Post subject: Trouble with UPML in FDTD
Posted: Thu Mar 09, 2006 10:32 am 
 
Lieutenant

Joined: Thu Mar 09, 2006 10:27 am
Posts: 1
Hi everyone,

I am currently working on a 3d FDTD code using UPML to truncate the physical domain. Nevertheless, when generating a plane wave (wave vector linear to z), it appears that the wavefront in the UPML becomes quickly angled, which results in large errors at late times. The algorithm I use for source generation is :

IF (Nt.LE.200) THEN
DO J=1,Ny
DO ii=1,Nx
Ex(ii,J,Nz/2)=Ex(ii,J,Nz/2)+(kappax(ii)-I*sigmax(ii)*30./(2.*Pi))/(kappay(j)-I*sigmay(j)*30./(2.*Pi))*(dcos(2.*pi*Nt/30.)+I*dsin(2.*pi*Nt/30.))*dexp(dble(-(Nt-200.)*(Nt-200.)/4000.))
ENDDO
ENDDO
ELSE
DO J=1,Ny
DO II=1,Nx
Ex(II,J,Nz/2)=Ex(II,J,Nz/2)+(kappax(ii)-I*sigmax(ii)*30./(2.*Pi))/(kappay(j)-I*sigmay(j)*30./(2.*Pi))*(dcos(2.*pi*Nt/30.)+I*dsin(2.*pi*Nt/30.))
ENDDO
ENDDO
ENDIF

You can notice that the electric field checks the conservation of the normal component of electric induction D at the interface with UPML. Nevertheless, this does not work well. Did someone experience similar problems with a FDTD code ? How did you solve them ? Thanks very much,

M.



Posted  11/12/2012
Windfreak Technologies PCB Directory (Assembly)
NorthEast RF Antenna Testing Services Berkeley Nucleonics Academy - RF Cafe
About RF Cafe
Kirt Blattenberger - RF Cafe Webmaster
Copyright: 1996 - 2024
Webmaster:
    Kirt Blattenberger,
    BSEE - KB3UON

RF Cafe began life in 1996 as "RF Tools" in an AOL screen name web space totaling 2 MB. Its primary purpose was to provide me with ready access to commonly needed formulas and reference material while performing my work as an RF system and circuit design engineer. The World Wide Web (Internet) was largely an unknown entity at the time and bandwidth was a scarce commodity. Dial-up modems blazed along at 14.4 kbps while typing up your telephone line, and a nice lady's voice announced "You've Got Mail" when a new message arrived...

All trademarks, copyrights, patents, and other rights of ownership to images and text used on the RF Cafe website are hereby acknowledged.

My Hobby Website:  AirplanesAndRockets.com

spacer

Please Support RF Cafe by purchasing my  ridiculously low−priced products, all of which I created.

These Are Available for Free