30.07.2013 Views

WAVES AND VIBRATIONS IN INHOMOGENEOUS STRUCTURES ...

WAVES AND VIBRATIONS IN INHOMOGENEOUS STRUCTURES ...

WAVES AND VIBRATIONS IN INHOMOGENEOUS STRUCTURES ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

5. Photonic Wire Waveguides<br />

Figure 4 shows a schematic setup of a PhW waveguide T-splitter. A straight waveguide with width w<br />

and of dielectric constant ε = n 2 , carries an incoming wave that propagates from left towards right. The<br />

input signal is to be split into two equal waves so that 50% of the incoming wave energy is transmitted<br />

through both the upper and lower output ports. We use a two-dimensional model of plane polarized<br />

light that for time-harmonic waves is governed by the Helmholtz equation:<br />

∇ · (A∇u) + Bω 2 u = 0, (1)<br />

where A = ε −1 and B = c −2 for TE-polarized waves, and A = 1 and B = εc −2 for TM-polarized waves.<br />

The frequency of the wave is denoted ω and the unknown magnetic or electric field is u.<br />

The boundary condition specifies an incoming wave at the left waveguide port:<br />

n · (A∇u) = 2iω √ ABu0, (2)<br />

where u0 is the specified amplitude of the incoming wave. In addition to the wave source, the computational<br />

domain in Figure 4 is embedded in a perfectly matched layer (PML) [15] that ensures a minimized<br />

reflection of outgoing waves. The governing equation in the PML is:<br />

∂<br />

∂x (sy A<br />

sx<br />

∂u ∂<br />

) +<br />

∂x<br />

∂y (sx<br />

sy<br />

A ∂u<br />

∂y ) + Bω2 sxsyu = 0, (3)<br />

where the two complex functions sx and sy specify the absorbing properties in the x- and y−direction,<br />

respectively. At the interface to the computational domain, sx = sy = 1 which transforms Eq. (3) into<br />

the original Eq. (1), thus creating a perfectly matching interface.<br />

5.1. Topology optimization<br />

We use the method of topology optimization [4] to find a suitable material distribution in the design<br />

domain (the area indicated with a question mark in Figure 4).<br />

The optimization problem is formulated as:<br />

⎧<br />

⎨<br />

⎩<br />

max min Φ1, Φ2<br />

subject to : S(̺, ω)p = f(ω)<br />

0 ≤ ̺e ≤ 1 e = 1, . . . , Nd.<br />

w<br />

n=1<br />

Figure 4: Design problem for a photonic wire waveguide T-splitter. The material distribution in the area<br />

indicated with a question mark is found using topology optimization. The objective is to distribute 50%<br />

of the input energy into the top and bottom output ports and thus eliminate reflection and radiation at<br />

the waveguide junction.<br />

4<br />

?<br />

n<br />

(4)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!