In the present thesis, computational methods for theanalysis and design of photonic bandgap structure areconsidered. Many numerical methods have been used to study suchstructures. Among them, the plane wave expansion method is veryoften used. Using this method, we show that inclusions ofelliptic air holes can be used effectively to obtain a largercomplete band gap for two-dimensional (2D) photonic crystals.An optimal design of a 2D photonic crystal is also consideredin the thesis using a combination of the plane wave expansionmethod and the conjugate gradient method. We find that amaximum complete 2D band gap can be obtained by connectingdielectric rods with veins for a photonic crystal with a squarelattice of air holes in GaAs. For some problems, such as defect modes, the plane waveexpansion method is extremely time-consuming. It seems that thefinite-difference time-domain (FDTD) method is promising, sincethe computational time is proportional to the number of thediscretization points in the computation domain (i.e., it is oforderN). A FDTD scheme in a nonorthogonal coordinate systemis presented in the thesis to calculate the band structure of a2D photonic crystal consisting of askew lattice. The algorithmcan easily be used for any complicated inclusion configuration,which can have both the dielectric and metallic constituents.The FDTD method is also applied to calculate the off-plane bandstructures of 2D photonic crystals in the present thesis. Wealso propose a numerical method for computing defect modes in2D crystals (with dielectric or metallic inclusions). Comparedto the FDTD transmission spectra method, our method reduces thecomputation time and memory significantly, and finds as manydefect modes as possible, including those that are not excitedby an incident plane wave in the FDTD transmission spectramethod. The FDTD method has also been applied to calculateguided modes and surface modes in 2D photonic crystals using acombination of the periodic boundary condition and theperfectly matched layer for the boundary treatment. Anefficient FDTD method, in which only real variables are used,is also proposed for the full-wave analysis of guided modes inphotonic crystal fibers. / QC 20100629

Although it is possible to directly solve an entire flow-acoustics problem in one computation, this approach remains prohibitively large in terms of the computational resource required for most practical applications. Aeroacoustic problems are therefore usually split into two parts; one consisting of the source computation and one of the source propagation. Although both these parts entail great challenges on the computational method, in terms of accuracy and efficiency, it is still better than the direct solution alternative. The source usually consists of highly turbulent flows, which for most cases will need to be, at least partly, resolved. Then, acoustic waves generated by these sources often have to be propagated for long distances compared to the wavelength and might be subjected to scattering by solid objects or convective effects by the flow. Numerical methods used solve these problems therefore have to possess low dispersion and dissipation error qualities for the solution to be accurate and resource efficient. The wave expansion method (WEM) is an efficient discretization technique, which is used for wave propagation problems. The method uses fundamental solutions to the wave operator in the discretization procedure and will thus produce accurate results at two to three points per wavelength. This thesis presents a method that uses the WEM in an aeroacoustic context. Addressing the propagation of acoustic waves and transfer of sources from flow to acoustic simulations. The proposed computational procedure is applied to a co-rotating vortex pair and a cylinder in cross-flow. Overall, the computed results agree well with analytical solutions. Although the WEM is efficient in terms of the spatial discretization, the procedure requires that a Moore-Penrose pseudo-inverse is evaluated at each unique node-neighbour stencil in the grid. This evaluation significantly slows the procedure. In this thesis, a method with a regular grid is explored to speed-up this process. / <p>QC 20161121</p>