Dear all,
am still playing with this complex dielectric issue. I am running
two ctl files listed below. The first, sea_fake.ctl runs in about 75s
having only real dielectrics ie water with zero conductiivty. The second
which simulates the same situation but with complex dielectric set by
using a seawater conductivitty of 4S/m takes over 200000s to run!
I'm using meep which responds to the meep -V command with:
Meep 0.10, Copyright (C) 2005 Massachusetts Insitute of Technology.
Using libctl 3.0.2 and Guile 1.6.7.
Why is there this huge difference ? I've tried asking for complex fields
which doesn't make much difference in either case. I'm running for 30000
as I'm using rather low frequeencies but comapred to their free space
wavelengths also rather small volumes.
best regards
Chris
"sea_fake.ctl"
;"sea_fake.ctl" to compute fields below a an exciting current in air above a
1.4km depth of seawater without any conductivity - ie just dielectric
; scale here is 1x2km setting the scale unit to 1000m
; thus for 100Hz we have a freqeuency of 3.3e-4
; and time units are ????
(set! geometry-lattice (make lattice (size 1 2 no-size)))
(set! geometry (list
; block of dielectric of constant epsilon =81
(make block (center 0 0.33) (size 1 1.4 infinity)
(material (make dielectric (epsilon 81)
)))
))
(set! sources (list
(make source
(src (make continuous-src (frequency 3.3e-4)))
(component Ez)
(center 0 -0.5 0))))
(set! pml-layers (list (make pml (thickness 0.1))))
(set! resolution 100)
(run-until 3000 ; should run for 1 periods of the source if the scaling is right
(at-beginning output-epsilon)
; (to-appended "hy" (at-every 100 output-hfield-y))
(to-appended "hxm" (at-every 100 output-hfield-x))
(to-appended "ezm" (at-every 100 output-efield-z))
(at-end output-efield-z)
(at-end output-hfield-y)
(at-end output-hfield-x)
)
"sea_deep.ctl"
;"sea_deep.ctl" to compute fields below a an exciting current in air above a
1.4km depth of seawater
; scale here is 1x2km setting the scale unit to 1m
; thus for 100Hz we have a ffreqeuency of 3.3e-4
; and time units are anybodys guess.....
(set! geometry-lattice (make lattice (size 1 2 no-size)))
(set! geometry (list
; block of dielectric of constant with polarisation to model seawater with
conductivity of 4Sm-1 and epsilon =81
(make block (center 0 0.3) (size 1.0 1.4 infinity)
(material (make dielectric (epsilon 81)
(polarizations
(make polarizability
(omega 1e-20) (gamma 1e13) (delta-epsilon 4.5e65))
)
)))
))
(set! sources (list
(make source
(src (make continuous-src (frequency 3.3e-7)))
(component Ez)
(center 0 -0.5 0))))
(set! pml-layers (list (make pml (thickness 1.0))))
(set! resolution 100)
(run-until 30000
(at-beginning output-epsilon)
; (to-appended "hy" (at-every 100 output-hfield-y))
(to-appended "hxm" (at-every 100 output-hfield-x))
(to-appended "ezm" (at-every 100 output-efield-z))
(at-end output-efield-z)
(at-end output-hfield-y)
(at-end output-hfield-x)
)
--
..................................................................
. Dr.Chris Stevens . Oxford University .
. Tel: +44 1865 283272 . Department of Engineering Science .
. Fax: +44 1865 273905 . Parks ROad .
. [EMAIL PROTECTED] . Oxford, OX1 3PJ, UK .
..................................................................
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
