Fast Lattice Green's Function Methods for Viscous Incompressible Flows on Unbounded Domains

In this thesis, a collection of novel numerical techniques culminating in a fast, parallel method for the direct numerical simulation of incompressible viscous flows around surfaces immersed in unbounded fluid domains is presented. At the core of all these techniques is the use of the fundamental solutions, or lattice Green’s functions, of discrete operators to solve inhomogeneous elliptic difference equations arising in the discretization of the three-dimensional incompressible Navier-Stokes equations on unbounded regular grids. In addition to automatically enforcing the natural free-space boundary conditions, these new lattice Green’s function techniques facilitate the implementation of robust staggered-Cartesian-grid flow solvers with efficient nodal distributions and fast multipole methods. The provable conservation and stability properties of the appropriately combined discretization and solution techniques ensure robust numerical solutions. Numerical experiments on thin vortex rings, low-aspect-ratio flat plates, and spheres are used verify the accuracy, physical fidelity, and computational efficiency of the present formulations.

Scattering by two degenerate anisotropic modes in square-lattice dielectric photonic crystals

We systematically investigate the square-lattice dielectric photonic crystals that have been used to demonstrate flat slab imaging experimentally. A right-handed Bloch mode is found in the left-handed frequency region by using the plane wave expansion method to analyze the photonic band structure and equifrequency contours. Using the multiple scattering theory, numerical simulations demonstrate that the left-handed mode and the right-handed mode are excited simultaneously by a point source and result in two kinds of transmitted waves. Impacted by the evanescent waves, superposition of these transmitted waves brings on complicated near field distributions such as the so-called imaging and its disappearance.

Optical properties of a square-lattice photonic crystal within the partial bandgap

Optical properties of a two-dimensional square-lattice photonic crystal are systematically investigated within the partial bandgap through anisotropic characteristics analysis and numerical simulation of field pattern. Using the plane-wave expansion method and Hellmann-Feynman theorem, the relationships between the incident and refracted angles for both phase and group velocities are calculated to analyze light propagation from air to photonic crystals. Three kinds of flat slab focusing are summarized and demonstrated by numerical simulations using the multiple scattering method. (c) 2007 Optical Society of America

Negative refraction in a honeycomb-lattice photonic crystal

In this paper we theoretically investigate a photonic crystal with dielectric rods in a honeycomb lattice. Two left-handed frequency regions are found in the second and third photonic band by using the plane wave expansion method to analyze the photonic band structure and equifrequency contours. Subwavelength imaging by the photonic crystal flat lens are systematically studied by numerical simulations using the multiple scattering method. Different from the photonic crystals with noncircular dielectric rods in air, this structure is almost isotropic at the optimal frequency for superlensing. As a comparison, flat slab focusing is also demonstrated at other frequencies in the two left-handed regions. (c) 2006 Elsevier Ltd. All rights reserved.

Band structure of photonic crystal fibers with silica rings in triangular lattice

The characteristics of the cladding band structure of air-core photonic crystal fibers with silica rings in triangular lattice are investigated by using a standard plane wave method. The numerical results show that light can be localized in the air core by the photonic band gaps of the fiber. By increasing the air-filling fraction, the band gap edges of the low frequency photonic band gaps shift to shorter wavelength.. whereas the band gap width decreases linearly. In order to make a specified light fall in the low frequency band gaps of the fiber, the interplay of the silica ring spacing and the air-filling fraction is also analyzed. It shows that the silica ring spacing increases monotonously when the air-filling fraction is increased, and the spacing range increases exponentially. This type fiber might have potential in infrared light transmission. (c) 2006 Elsevier B.V. All rights reserved.

Convergence Properties and Fixed Points of Two General Iterative Schemes with Composed Maps in Banach Spaces with Applications to Guaranteed Global Stability

This paper investigates the boundedness and convergence properties of two general iterative processes which involve sequences of self-mappings on either complete metric or Banach spaces. The sequences of self-mappings considered in the first iterative scheme are constructed by linear combinations of a set of self-mappings, each of them being a weighted version of a certain primary self-mapping on the same space. The sequences of self-mappings of the second iterative scheme are powers of an iteration-dependent scaled version of the primary self-mapping. Some applications are also given to the important problem of global stability of a class of extended nonlinear polytopic-type parameterizations of certain dynamic systems.