Maria Bayard Dühring - Solid Mechanics
Maria Bayard Dühring - Solid Mechanics
Maria Bayard Dühring - Solid Mechanics
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