Dear MEEP user
I am very sorry for non-sense question. Actually, the x-axis is a number of
points which is surely increase by increasing the resolution. However, the
distance along the line doesn't change.
Best regard,
P. Tapsanit
________________________________
From: "[email protected]"
<[email protected]>
To: [email protected]
Sent: Friday, March 15, 2013 1:00 AM
Subject: meep-discuss Digest, Vol 85, Issue 8
Send meep-discuss mailing list submissions to
[email protected]
To subscribe or unsubscribe via the World Wide Web, visit
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
or, via email, send a message with subject or body 'help' to
[email protected]
You can reach the person managing the list at
[email protected]
When replying, please edit your Subject line so it is more specific
than "Re: Contents of meep-discuss digest..."
Today's Topics:
1. Re: output-epsilon problem (Tim Saucer)
2. Why the distance along a line change by chaning the
resolution? (Epsilon Piyawath)
----------------------------------------------------------------------
Message: 1
Date: Wed, 13 Mar 2013 16:23:40 -0400
From: Tim Saucer <[email protected]>
To: [email protected]
Subject: Re: [Meep-discuss] output-epsilon problem
Message-ID: <[email protected]>
Content-Type: text/plain; charset="us-ascii"
Hi Alex,
I believe you need to change your size in the y direction to infinity rather
than no-size.
(set! geometry (list
(make block
(center comb_width no-size (- 1 (/ comb_dept 2)) )
(size comb_width infinity comb_dept)
(material (make dielectric (epsilon 1.8))) )
(make block
(center 0.5 no-size (/ (- 1 comb_dept) 2) )
(size 1.0 infinity (- 1 comb_dept) )
(material (make dielectric (epsilon 1.8))) )
))
Attached is a .png of the permittivity with this change.
Best wishes,
Tim
Tim Saucer
Physics 236/241/261 Lead GSI
Sih Research Group PhD Candidate
Office: 4409 Randall, Lab: 1277 Randall
Department of Physics, University of Michigan
Alex <to.sanja@...> writes:
>
> Hello.
> I try to define geomentry of 2D structure with relative permittivity
> 1.8. But then I look h5 output file by h5totxt I see one "1,1,1,1".
> Ctl file:
> (set! geometry-lattice (make lattice (size 1 no-size 1))); Specifies
> the size of the unit cell
> (define comb_dept (/ 9 10))
> (define comb_width (/ 5 10) )
>
> (set! geometry (list
> (make block (center comb_width no-size (- 1 (/ comb_dept 2)) ) (size
> comb_width no-size comb_dept) (material (make dielectric (epsilon
> 1.8))) )
> (make block (center 0.5 no-size (/ (- 1 comb_dept) 2) ) (size 1.0
> no-size (- 1 comb_dept) ) (material (make dielectric (epsilon 1.8))) )
> ))
> (display geometry)
> (set! resolution 32)
> (set! k-point (vector3 1 0 0))
> (set! ensure-periodicity true)
> (define-param fcen 0.5)
> (define-param df 3.0)
> (set! sources (list (make source (src (make gaussian-src (frequency
> fcen) (fwidth df))) (component Hz) (center 0.1234 0))))
>
> (set-param! k-point (vector3 0.4 0))
> (run-sources+ 300 (after-sources (harminv Hz (vector3 0.1234) fcen df)))
> (define-param k-interp 50)
> ;(run-k-points 300 (interpolate k-interp (list (vector3 0) (vector3 0.5))))
> (run-k-points 300 (append
> (interpolate 15 (list (vector3 0) (vector3 0.1)))
> (interpolate 5 (list (vector3 0.1) (vector3 0.4)))
> (interpolate 15 (list (vector3 0.4) (vector3 0.5)))
> ) )
> (run-until 200 (at-beginning output-epsilon) (at-end
> output-efield-z))
>
> OUTPUT:
> [alex <at> nout meep]$ h5totxt 1d_walls-eps-000000.00.h5
> 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1
> ....
>
> Alex Friman
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL:
<http://ab-initio.mit.edu/pipermail/meep-discuss/attachments/20130313/2bc1a791/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: temp-eps-000000.00.png
Type: image/png
Size: 231 bytes
Desc: not available
URL:
<http://ab-initio.mit.edu/pipermail/meep-discuss/attachments/20130313/2bc1a791/attachment-0001.png>
------------------------------
Message: 2
Date: Wed, 13 Mar 2013 22:05:27 -0700 (PDT)
From: Epsilon Piyawath <[email protected]>
To: "[email protected]" <[email protected]>
Subject: [Meep-discuss] Why the distance along a line change by
chaning the resolution?
Message-ID:
<[email protected]>
Content-Type: text/plain; charset="iso-8859-1"
Dear MEEP user,
I have a short question. I've been trying to get the focal length of a len. So,
I calculate the intensity along the particular line (e.g., (0,sy,0)), and then
I use h5totxt to convert .h5 file to .dat file.? I want the converged solution
of the intensity along this line, so I change the resolution and then obtain
the intensity of each resolution.? I use gnuplot to plot the intensity and I
found that the distance of this line increases(x-axis) by increasing the
resolution even the computational cell is fixed. The intensity curved is
shifted to the right according to the increase of the distance. The distance
should be constant. Do you know how to fix this problem and what is its origin.
I give you the ctl file as follows:
;length scale a = 100 nm
(define pi 3.141593)
(define-param lam 3.0)
;lens geometry
(define-param lenlen (* 2 lam))
(define-param lenwidth (* 2 lam))
;lens slab geometry
(define-param dlenslab (* 1.0 lam))
;define the computational cell
(define-param dpml lam) ; thickness of PML layers/shold be larger than
wavelength
(define-param pad lam)
(define-param sx (* 2 (+ lam pad)))
(define-param sy (* 2 (+ lam pad)))
(define-param sx0 (+ sx (* 2 dpml)))
(define-param sy0 (+ sy (* 2 dpml)))
(set! geometry-lattice (make lattice (size sx0 sy0 no-size)))?? ;2D
computational cell
;line source parameters
(define-param llen (* 1.0 lenlen))??????????????? ;distance between point
sources
(define-param lx 0)
(define-param ly (- (- (* 0.5 sy) pad)))
(define e0 1e+5)
(define epsquartz (make dielectric (epsilon 2.174)))
(set! eps-averaging? true)
(set! output-single-precision? true)
(set! force-complex-fields? true)
(set! Courant 0.6)
(set! geometry
? (if structure?
??? (list
?????? ;cell
?????? (make block (center 0 0 0) (size sx sy infinity) (material air))
?????? ;lens
?????? (make ellipsoid (center 0 0 0) (size lenlen lenwidth infinity) (e1 1 0
0) (e2 0 1 0) (e3 0 0 1) (material epsquartz))
?????? ;lens mask
?????? (make block (center 0 (* 0.25 sy) 0) (size sx (* 0.5 sy) infinity)
(material air))
?????? ;lens slab
?????? (make block (center 0 (* 0.25 dlenslab) 0) (size lenlen (* 0.5 dlenslab)
infinity) (material epsquartz))
??? )
??? (list
????? (make block (center 0 0 0) (size infinity infinity infinity) (material
air))?
??? )
? )
)
;define PML layers in all directions
(set! pml-layers (list (make pml (thickness dpml) )))?
;set symmetry associated with the source
(set! symmetries
????? (list
????? ??? ??? (make mirror-sym (direction X) (phase -1))???????????????????
;odd along x-axis
????? )
)
(set! sources
? (list
??? (make source (src (make continuous-src (wavelength lam)))
???????????????? (component Ex)
???????????????? (amplitude e0)
???????????????? (center lx ly 0)
???????????????? (size llen 0 0)?????
???? )
? )
)
;define the resolution
(set! resolution resl)
(define (esq)
? (output-real-field-function "esq" (list Ex Ey Ez)
??? (lambda (r ex ey ez) (expt (/ (sqrt (+ (* ex (conj ex)) (* ey (conj ey)) (*
ez (conj ez)) )) e0) 2) )
? )
)
(run-until 200
? (at-beginning
??? (in-volume (volume (center 0 0 0) (size sx sy 0)) output-epsilon)
? )
? ;fiels distribution
? (at-end
??? (to-appended "esq-line-vertical"
????? (in-volume (volume (center 0 (* 0.25 sy) 0) (size 0 (* 0.5 sy) 0) ) esq)
??? )
)??
Best regard,
P. Tapsanit
-------------- next part --------------
An HTML attachment was scrubbed...
URL:
<http://ab-initio.mit.edu/pipermail/meep-discuss/attachments/20130313/d202791c/attachment-0001.html>
------------------------------
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
End of meep-discuss Digest, Vol 85, Issue 8
*******************************************
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
