Accurate multigrid solution of the Euler equations on unstructured and adaptive meshes

A method for accurately solving inviscid compressible flow in the subcritical and supercritical regimes about complex configurations is presented. The method is based on the use of unstructured triangular meshes in two dimensions, and special emphasis is placed on the accuracy and efficiency of the solutions. High accuracy is achieved by careful scaling of the artificial dissipation terms, and by reformulating the inner and outer boundary conditions for both the convective and dissipative operators. An adaptive grid refinement strategy is presented which enhances the solution accuracy for complex flows. When coupled with an unstructured multigrid algorithm, this method is shown to produce an efficient solver for flows about arbitrary configurations

A C0 Finite Element Method For The Biharmonic Problem In A Polygonal Domain

This dissertation studies the biharmonic equation with Dirichlet boundary conditions in a polygonal domain. The biharmonic problem appears in various real-world applications, for example in plate problems, human face recognition, radar imaging, and hydrodynamics problems. There are three classical approaches to discretizing the biharmonic equation in the literature: conforming finite element methods, nonconforming finite element methods, and mixed finite element methods. We propose a mixed finite element method that effectively decouples the fourth-order problem into a system of one steady-state Stokes equation and one Poisson equation. As a generalization to the above-decoupled formulation, we propose another decoupled formulation using a system of two Poison equations and one steady-state Stokes equation. We solve Poisson equations using linear and quadratic Lagrange\u27s elements and the Stokes equation using Hood-Taylor elements and Mini finite elements. It is shown that the solution of each system is equivalent to that of the original fourth-order problem on both convex and non-convex polygonal domains. Two finite element algorithms are, in turn, proposed to solve the decoupled systems. Solving this problem in a non-convex domain is challenging due to the singularity occurring near re-entrant corners. We introduce a weighted Sobolev space and a graded mesh refine Algorithm to attack the singularity near re-entrant corners. We show the regularity results of each decoupled system in both Sobolev space and weighted Sobolev space. We derive the $H^1$ and $L^2$ error estimates for the numerical solutions on quasi-uniform and graded meshes. We present various numerical test results to justify the theoretical findings. Given the availability of fast Poisson solvers and Stokes solvers, our Algorithm is a relatively easy and cost-effective alternative to existing algorithms for solving the biharmonic equation