Physical Sciences and Mathematics Commons^™

Simulating Dislocation Densities With Finite Element Analysis, Ja'nya Breeden, Dow Drake, Saurabh Puri

REU Final Reports

A one-dimensional set of nonlinear time-dependent partial differential equations developed by Acharya (2010) is studied to observe how differing levels of applied strain affect dislocation walls. The framework of this model consists of a convective and diffusive term which is used to develop a linear system of equations to test two methods of the finite element method. The linear system of partial differential equations is used to determine whether the standard or Discontinuous Galerkin method will be used. The Discontinuous Galerkin method is implemented to discretize the continuum model and the results of simulations involving zero and non-zero applied strain …

Filtered Subspace Iteration For Selfadjoint Operators, Jay Gopalakrishnan, Luka Grubišić, Jeffrey S. Ovall

Portland Institute for Computational Science Publications

We consider the problem of computing a cluster of eigenvalues (and its associated eigenspace) of a (possibly unbounded) selfadjoint operator in a Hilbert space. A rational function of the operator is constructed such that the eigenspace of interest is its dominant eigenspace, and a subspace iteration procedure is used to approximate this eigenspace. The computed space is then used to obtain approximations of the eigenvalues of interest. An eigenvalue and eigenspace convergence analysis that considers both iteration error and dis- cretization error is provided. A realization of the proposed approach for a model second-order elliptic operator is based on a …

Accuracy Of Wave Speeds Computed From The Dpg And Hdg Methods For Electromagnetic And Acoustic Waves, Nicole Michelle Olivares

Dissertations and Theses

We study two finite element methods for solving time-harmonic electromagnetic and acoustic problems: the discontinuous Petrov-Galerkin (DPG) method and the hybrid discontinuous Galerkin (HDG) method.

The DPG method for the Helmholtz equation is studied using a test space normed by a modified graph norm. The modification scales one of the terms in the graph norm by an arbitrary positive scaling parameter. We find that, as the parameter approaches zero, better results are obtained, under some circumstances. A dispersion analysis on the multiple interacting stencils that form the DPG method shows that the discrete wavenumbers of the method are complex, explaining …

Discrete Stability Of Dpg Methods, Ammar Harb

Dissertations and Theses

This dissertation presents a duality theorem of the Aubin-Nitsche type for discontinuous Petrov Galerkin (DPG) methods. This explains the numerically observed higher convergence rates in weaker norms. Considering the specific example of the mild-weak (or primal) DPG method for the Laplace equation, two further results are obtained. First, for triangular meshes, the DPG method continues to be solvable even when the test space degree is reduced, provided it is odd. Second, a non-conforming method of analysis is developed to explain the numerically observed convergence rates for a test space of reduced degree. Finally, for rectangular meshes, the test space is …

Comparison Of Mesh And Meshless Methods For Partial Differential Equations Of Galerkin Form, Wallace F. Atterberry

UNLV Theses, Dissertations, Professional Papers, and Capstones

There are two purposes of this research project. The first purpose is to compare two types of Galerkin methods: The finite element mesh method and moving least sqaures meshless Galerkin (EFG) method. The second purpose of this project is to determine if a hybrid between the mesh and meshless method is beneficial.

This manuscript will be divided into three main parts. The first part is chapter one which develops the finite element method. The second part (Chapter two) will be developing the meshless method. The last part will provide a method for combining the mesh and meshless methods for a …

Analysis Of Hdg Methods For Stokes Flow, Bernardo Cockburn, Jay Gopalakrishnan, Ngoc Cuong Nguyen, Jaume Peraire, Francisco-Javier Sayas

Mathematics and Statistics Faculty Publications and Presentations

In this paper, we analyze a hybridizable discontinuous Galerkin method for numerically solving the Stokes equations. The method uses polynomials of degree $ k$ for all the components of the approximate solution of the gradient-velocity-pressure formulation. The novelty of the analysis is the use of a new projection tailored to the very structure of the numerical traces of the method. It renders the analysis of the projection of the errors very concise and allows us to see that the projection of the error in the velocity superconverges. As a consequence, we prove that the approximations of the velocity gradient, the …

Analysis Of The Dpg Method For The Poisson Equation, Leszek Demkowicz, Jay Gopalakrishnan

Mathematics and Statistics Faculty Publications and Presentations

We give an error analysis of the recently developed DPG method applied to solve the Poisson equation and a convection-dffusion problem. We prove that the method is quasioptimal. Error estimates in terms of both the mesh size h and the polynomial degree p (for various element shapes) can be derived from our results. Results of extensive numerical experiments are also presented.

A Class Of Discontinuous Petrov–Galerkin Methods. Part Iv: The Optimal Test Norm And Time-Harmonic Wave Propagation In 1d., Jeffrey Zitelli, Leszek Demkowicz, Jay Gopalakrishnan, D. Pardo, V. M. Calo

Mathematics and Statistics Faculty Publications and Presentations

The phase error, or the pollution effect in the finite element solution of wave propagation problems, is a well known phenomenon that must be confronted when solving problems in the high-frequency range. This paper presents a new method with no phase errors for one-dimensional (1D) time-harmonic wave propagation problems using new ideas that hold promise for the multidimensional case. The method is constructed within the framework of the discontinuous Petrov–Galerkin (DPG) method with optimal test functions. We have previously shown that such methods select solutions that are the best possible approximations in an energy norm dual to any selected test …

A Class Of Discontinuous Petrov–Galerkin Methods. Ii. Optimal Test Functions, Leszek Demkowicz, Jay Gopalakrishnan

Mathematics and Statistics Faculty Publications and Presentations

We lay out a program for constructing discontinuous Petrov–Galerkin (DPG) schemes having test function spaces that are automatically computable to guarantee stability. Given a trial space, a DPG discretization using its optimal test space counterpart inherits stability from the well posedness of the undiscretized problem. Although the question of stable test space choice had attracted the attention of many previous authors, the novelty in our approach lies in the fact we identify a discontinuous Galerkin (DG) framework wherein test functions, arbitrarily close to the optimal ones, can be locally computed. The idea is presented abstractly and its feasibility illustrated through …

A Class Of Discontinuous Petrov–Galerkin Methods. Part I: The Transport Equation, Leszek Demkowicz, Jay Gopalakrishnan

Mathematics and Statistics Faculty Publications and Presentations

Considering a simple model transport problem, we present a new finite element method. While the new method fits in the class of discontinuous Galerkin (DG) methods, it differs from standard DG and streamline diffusion methods, in that it uses a space of discontinuous trial functions tailored for stability. The new method, unlike the older approaches, yields optimal estimates for the primal variable in both the element size h and polynomial degree p, and outperforms the standard upwind DG method.

A Projection-Based Error Analysis Of Hdg Methods, Jay Gopalakrishnan, Bernardo Cockburn, Francisco-Javier Sayas

Mathematics and Statistics Faculty Publications and Presentations

We introduce a new technique for the error analysis of hybridizable discontinuous Galerkin (HDG) methods. The technique relies on the use of a new projection whose design is inspired by the form of the numerical traces of the methods. This renders the analysis of the projections of the discretization errors simple and concise. By showing that these projections of the errors are bounded in terms of the distance between the solution and its projection, our studies of influence of the stabilization parameter are reduced to local analyses of approximation by the projection. We illustrate the technique on a specific HDG …

The Derivation Of Hybridizable Discontinuous Galerkin Methods For Stokes Flow, Bernardo Cockburn, Jay Gopalakrishnan

Mathematics and Statistics Faculty Publications and Presentations

In this paper, we introduce a new class of discontinuous Galerkin methods for the Stokes equations. The main feature of these methods is that they can be implemented in an efficient way through a hybridization procedure which reduces the globally coupled unknowns to certain approximations on the element boundaries. We present four ways of hybridizing the methods, which differ by the choice of the globally coupled unknowns. Classical methods for the Stokes equations can be thought of as limiting cases of these new methods.

Unified Hybridization Of Discontinuous Galerkin, Mixed, And Continuous Galerkin Methods For Second Order Elliptic Problems, Bernardo Cockburn, Jay Gopalakrishnan, Raytcho Lazarov

Mathematics and Statistics Faculty Publications and Presentations

We introduce a unifying framework for hybridization of finite element methods for second order elliptic problems. The methods fitting in the framework are a general class of mixed-dual finite element methods including hybridized mixed, continu- ous Galerkin, non-conforming and a new, wide class of hybridizable discontinuous Galerkin methods. The distinctive feature of the methods in this framework is that the only globally coupled degrees of freedom are those of an approximation of the solution defined only on the boundaries of the elements. Since the associated matrix is sparse, symmetric and positive definite, these methods can be efficiently implemented. Moreover, the …

Locally Conservative Fluxes For The Continuous Galerkin Method, Bernardo Cockburn, Jay Gopalakrishnan, Haiying Wang

Mathematics and Statistics Faculty Publications and Presentations

The standard continuous Galerkin (CG) finite element method for second order elliptic problems suffers from its inability to provide conservative flux approximations, a much needed quantity in many applications. We show how to overcome this shortcoming by using a two step postprocessing. The first step is the computation of a numerical flux trace defined on element inter- faces and is motivated by the structure of the numerical traces of discontinuous Galerkin methods. This computation is non-local in that it requires the solution of a symmetric positive definite system, but the system is well conditioned independently of mesh size, so it …

A Multilevel Discontinuous Galerkin Method, Jay Gopalakrishnan, Guido Kanschat

Mathematics and Statistics Faculty Publications and Presentations

A variable V-cycle preconditioner for an interior penalty finite element discretization for elliptic problems is presented. An analysis under a mild regularity assumption shows that the preconditioner is uniform. The interior penalty method is then combined with a discontinuous Galerkin scheme to arrive at a discretization scheme for an advection-diffusion problem, for which an error estimate is proved. A multigrid algorithm for this method is presented, and numerical experiments indicating its robustness with respect to diffusion coefficient are reported.