Hi all!
What I'm trying to do is simulate a tx/rx system as
follow:
2 identical patch antennas (one transmitting and the
other receiving) with dimension about 41*0.15*88
(expressed in mm). The distance between the two
antennas is 60 cm (in Y direction).
To contain the structure below I made a computational
cell of dimension 100*650*150 (I choose pml
thickness=10). Then I set resolution=1.
What I want to do is compute the field, saying Ey, at
each unit of time.
With these conditions my simulation takes too much
time...I would know if it is possible to make
subgridding for example for the vacuum space between
the two antennas...Otherwise what can I do to speed up
my simulation?
thanks in advance.
Federica
;***********************************************
; *
; Planar Differential Elliptical UWB Antenna *
; *
;***********************************************
;------------------------------
; define parameters |
;------------------------------
;fattore di scala a=1e-3 ---> tutto espresso in unitÃ
di mm
;dimensioni lungo X
(define-param Rx1 20.320) ; raggio ellisse maggiore
lungo X
(define-param Rx2 9.906) ; raggio ellisse media lungo
X
(define-param Rx3 2.8448) ; raggio ellisse piccola
lungo X
(define-param s1 0.3302) ; distanza ellisse maggiore
da bordo piastra lungo X
;dimensioni lungo Z
(define-param Rz1 15.240) ; raggio ellisse maggiore
lungo Z
(define-param Rz2 9.8806) ; raggio ellisse media lungo
Z
(define-param Rz3 4.2164) ; raggio ellisse piccola
lungo Z
(define-param h1 5.080) ; semi-distanza tra le 2
ellissi maggiori
(define-param h2 0.4318) ; distanza minima tra ellisse
media ed ellisse piccola
(define-param h3 8.5598) ; distanza ellisse maggiore
da bordo piastra lungo Z
(define-param w 0.254) ; feeding gap
;parametri substrato a 5GHZ
(define-param h 0.1016) ; spessore del substrato
dielettrico
(define-param eps_r 3.36) ; costante dielettrica
(define-param tan_delta 0.0037) ; loss tangent
;spessore piano metallico
(define-param hm 0.035052)
;spessore PML
(define-param dpml 10)
(define D 600) ; distanza antenne
;dimensioni computational cell
(define Cx 100) ; dimensione cella lungo X
(define Cy 650) ; dimensione cella lungo Y
(define Cz 150) ; dimensione cella lungo Z
;definizione dimensioni geometria
;blocco dielettrico
(define Ld_x (* 2 (+ Rx1 s1)))
(define-param Ld_y h)
(define Ld_z (* 2 (+ h3 h1 (* 2 Rz1))))
(define Cd_y (+ (/ D 2) hm (/ h 2)))
(define minus_Cd_y (- 0 Cd_y))
;blocco metallo
(define Lm_x (* 2 (+ Rx1 s1)))
(define-param Lm_y hm)
(define Lm_z (* 2 (+ h3 h1 (* 2 Rz1))))
(define Cm_y (+ (/ D 2) (/ hm 2)))
(define minus_Cm_y (- 0 Cm_y))
;ellisse maggiore
(define LE_x (* 2 Rx1)); asse ellisse maggiore lungo x
(define LE_z (* 2 Rz1)); asse ellisse maggiore lungo z
(define CE_z (+ h1 Rz1))
;ellisse minore
(define Le_x (* 2 Rx3)); asse ellisse minore lungo x
(define Le_z (* 2 Rz3)); asse ellisse minore lungo z
(define Ce_z (+ h1 w h2 Rz3))
;------------------------------
; define geometry |
;------------------------------
(set! geometry-lattice (make lattice (size Cx Cy Cz)))
(set! geometry (list
; antenna TX
; blocco dielettrico
(make block (center 0 minus_Cd_y 0)
(size Ld_x Ld_y Ld_z)
(material (make dielectric
(epsilon eps_r))))
; blocco metallo
(make block (center 0 minus_Cm_y 0)
(size Lm_x Lm_y Lm_z)
(material (make
perfect-metal)))
; ellisse maggiore 1-2 quadrante
(make ellipsoid (center 0 minus_Cm_y
CE_z)
(size LE_x hm LE_z)
(material (make
dielectric (epsilon eps_r))))
; ellisse minore 1-2 quadrante
(make ellipsoid (center 0 minus_Cm_y
Ce_z)
(size Le_x hm Le_z)
(material (make
dielectric (epsilon eps_r))))
; ellisse maggiore 3-4 quadrante
(make ellipsoid (center 0 minus_Cm_y
(- 0 CE_z))
(size LE_x hm LE_z)
(material (make
dielectric(epsilon eps_r))))
; ellisse minore 3-4 quadrante
(make ellipsoid (center 0 minus_Cm_y
(- 0 Ce_z))
(size Le_x hm Le_z)
(material
(make dielectric (epsilon eps_r))))
; antenna RX
; blocco dielettrico
(make block (center 0 Cd_y 0)
(size Ld_x Ld_y Ld_z)
(material (make dielectric
(epsilon eps_r))))
; blocco metallo
(make block (center 0 Cm_y 0)
(size Lm_x Lm_y Lm_z)
(material (make
perfect-metal)))
; ellisse maggiore 1-2 quadrante
(make ellipsoid (center 0 Cm_y CE_z)
(size LE_x hm LE_z)
(material (make
dielectric (epsilon eps_r))))
; ellisse minore 1-2 quadrante
(make ellipsoid (center 0 Cm_y Ce_z)
(size Le_x hm Le_z)
(material (make
dielectric (epsilon eps_r))))
; ellisse maggiore 3-4 quadrante
(make ellipsoid (center 0 Cm_y (- 0
CE_z))
(size LE_x hm LE_z)
(material (make
dielectric(epsilon eps_r))))
; ellisse minore 3-4 quadrante
(make ellipsoid (center 0 Cm_y (- 0
Ce_z))
(size Le_x hm Le_z)
(material (make
dielectric (epsilon eps_r))))
));end geometry
;------------------------------
; define source |
;------------------------------
(define-param 2PIsqr (sqrt (* 2 pi)))
(define-param f 0.15)
(define-param w 15.3)
(define Cr_y (/ D 2)); coord. y punto ricezione
(define Cs_y (- 0 (/ D 2))); coordinata y sorgente
; derivata quinta dellimpulso gaussiano
(define ((GM_5 sigma on_time) t1)
(begin;
(let ((t (- t1 (/ on_time 2))))
(if (and (<= t1 on_time) (>= t1 0))
(* (- (/ (* 10 ( expt t 3)) (* 2PIsqr (expt sigma
9)))
(/ ( expt t 5) (* 2PIsqr (expt sigma
11)))
(/ (* 15 t) (* 2PIsqr (expt sigma
7))))
(exp (/ (- 0 (expt t 2)) (* 2 (expt
sigma 2)))))
(+ 0 0))
) ))
(set! sources (list (make source (src (make
custom-src (src-func (GM_5 w 300)) (end-time 1000)))
(component Ey) (center 0 Cs_y 0))))
;inseriamo un PML per ridurre le riflessini
(set! pml-layers (list (make pml (thickness dpml))))
;definiamo la simmetria
(set! symmetries (list (make mirror-sym (direction X))
(make mirror-sym (direction
Z))))
(set! resolution 1)
(define (my_func)
(print (get-field-point Ey (vector3 0 0 0)) "\n"))
(run_until 600 (at-every 1 my_func))
___________________________________
L'email della prossima generazione? Puoi averla con la nuova Yahoo! Mail:
http://it.docs.yahoo.com/nowyoucan.html
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
