Date of Award

Summer 2011

Document Type


Degree Name

Doctor of Philosophy (PhD)


Mathematics and Statistics


Computational and Applied Mathematics

Committee Director

Fang Q. Hu

Committee Member

Chester Grosch

Committee Member

Richard D. Noren

Committee Member

Yan Peng

Committee Member

Hideaki Kaneko


Perfectly Matched Layer (PML) absorbing boundary conditions were first proposed by Berenger in 1994 for the Maxwell's equations of electromagnetics. Since Hu first applied the method to Euler's equations in 1996, progress made in the application of PML to Computational Aeroacoustics (CAA) includes linearized Euler equations with non-uniform mean flow, non-linear Euler equations, flows with an arbitrary mean flow direction, and non-linear clavier-Stokes equations. Although Boltzmann-BGK methods have appeared in the literature and have been shown capable of simulating aeroacoustics phenomena, very little has been done to develop absorbing boundary conditions for these methods. The purpose of this work was to extend the PML methodology to the discrete velocity Boltzmann-BGK equation (DVBE) for the case of a horizontal mean flow in two and three dimensions. The proposed extension of the PML has been accomplished in this dissertation. Both split and unsplit PML absorbing boundary conditions are presented in two and three dimensions. A finite difference and a lattice model are considered for the solution of the PML equations. The linear stability of the PML equations is investigated for both models. The small relaxation time needed for the discrete velocity Boltzmann-BC4K model to solve the Euler equations renders the explicit Runge-Kutta schemes impractical. Alternatively, implicit-explicit Runge-Kutta (IMEX) schemes are used in the finite difference model and are implemented explicitly by exploiting the special structure of the Boltzmann-BGK equation. This yields a numerically stable solution by the finite difference schemes. As the lattice model proves to be unstable, a coupled model consisting of a lattice Boltzmann (LB) method for the Ulterior domain and an IMEX finite difference method for the PML domains is proposed and investigated. Numerical examples of acoustic and vorticity waves are included to support the validity of the PML equations. In each example, accurate solutions are obtained, supporting the conclusion that PML is an effective absorbing boundary condition.