28.08.2013 Views

Maria Bayard Dühring - Solid Mechanics

Maria Bayard Dühring - Solid Mechanics

Maria Bayard Dühring - Solid Mechanics

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.

28 Chapter 5 Design of photonic-crystal fibers by topology optimization [P2]<br />

Figure 5.1 Geometry of the photonic crystal fiber used in the optimization. Black<br />

indicates silica and white is air. The design domain is Ωd and Ωc is the core region<br />

where the objective function is optimized.<br />

radius of the air core in the center is R = Λ. The boundary condition along all<br />

the borders is the perfect electric conductor. The chosen optical wavelength in free<br />

space is λ0 = 2 µm and the refractive indexes for this wavelength are na = 1 for<br />

air and ns = 1.43791 + 0.0001i for silica. The imaginary part of ns indicates the<br />

absorbing coefficient αs. It is assumed that the propagating optical modes have<br />

harmonic solutions on the form<br />

Hp,ν(x1, x2, x3) = Hp,ν(x1, x2)e −iβνx3 , (5.1)<br />

where Hp,ν is the magnetic field of the optical wave and βν is the propagation<br />

constant for a given optical mode ν. This is entered into the time-harmonic wave<br />

equation<br />

<br />

<br />

∂<br />

∂Hp,ν<br />

eijk bklblmemnp − k<br />

∂xj<br />

∂xn<br />

2 0Hi,ν = 0, (5.2)<br />

where k0 = 2π/λ0 is the free space propagation constant, eijk is the alternating<br />

symbol and biknkj = δij. For the problem considered in this chapter the materials<br />

are assumed to be isotropic and therefore bklblm reduces to n −2 . For the given value of<br />

k0 the propagation constant βν for the possible modes are found by solving the wave<br />

equation as an eigenvalue problem. In order to define the problem for the topology<br />

optimization, the propagation constant for the guided mode β1 is calculated. The<br />

wave equation is then solved in a way where both k0 and β1 are fixed, and the guided<br />

mode is excited by applying a magnetic source term in the center of the fiber, which<br />

is directed in the horizontal direction. The design variable ξ takes the value 0 for<br />

air and 1 for silica and the refractive index is interpolated linearly between the two<br />

material phases. The objective function Φ is an energy measure and is defined as

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

Saved successfully!

Ooh no, something went wrong!