Source author record

Cody Lorton

Cody Lorton appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA

Catalog footprint

What is connected

4works

1topics

2close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2016arXiv

An efficient Monte Carlo interior penalty discontinuous Galerkin method for elastic wave scattering in random media

This paper develops and analyzes an efficient Monte Carlo interior penalty discontinuous Galerkin (MCIP-DG) method for elastic wave scattering in random media. The method is constructed based on a multi-modes expansion of the solution of the governing random partial differential equations. It is proved that the mode functions satisfy a three-term recurrence system of partial differential equations (PDEs) which are nearly deterministic in the sense that the randomness only appears in the right-hand side source terms, not in the coefficients of the PDEs. Moreover, the same differential operator applies to all mode functions. A proven unconditionally stable and optimally convergent IP-DG method is used to discretize the deterministic PDE operator, an efficient numerical algorithm is proposed based on combining the Monte Carlo method and the IP-DG method with the $LU$ direct linear solver. It is shown that the algorithm converges optimally with respect to both the mesh size $h$ and the sampling number $M$, and practically its total computational complexity is only amount to solving very few deterministic elastic Helmholtz equations using the $LU$ direct linear solver. Numerically experiments are also presented to demonstrate the performance and key features of the proposed MCIP-DG method.

preprint2015arXiv

An unconditionally stable discontinuous Galerkin method for the elastic Helmholtz equations with large frequency

In this paper we propose and analyze an interior penalty discontinuous Galerkin (IP-DG) method using piecewise linear polynomials for the elastic Helmholtz equations with the first order absorbing boundary condition. It is proved that the sesquilinear form for the problem satisfies a generalized weak coercivity property, which immediately infers a stability estimate for the solution of the differential problem in all frequency regimes. It is also proved that the proposed IP-DG method is unconditionally stable with respect to both frequency $ω$ and mesh size $h$. Sub-optimal order (with respect to $h$) error estimates in the broken $H^1$-norm and in the $L^2$-norm are obtained in all mesh regimes. These estimate improve to optimal order when the mesh size $h$ is restricted to the pre-asymptotic regime (i.e., $ω^βh =O(1)$ for some $1\leq β<2$). The novelties of the proposed IP-DG method include: first, the method penalizes not only the jumps of the function values across the element edges but also the jumps of the normal derivatives; second, the penalty parameters are taken as complex numbers with positive imaginary parts. In order to establish the desired unconditional stability estimate for the numerical solution, the main idea is to exploit a (simple) property of linear functions to overcome the main difficulty caused by non-Hermitian nature and strong indefiniteness of the Helmholtz-type problem. The error estimate is then derived using a nonstandard technique adapted from \cite{Feng_Wu09}. Numerical experiments are also presented to validate the theoretical results and to numerically examine the pollution effect (with respect to $ω$) in the error bounds.

preprint2014arXiv

An efficient numerical method for acoustic wave scattering in random media

This paper is concerned with developing efficient numerical methods for acoustic wave scattering in random media which can be expressed as random perturbations of homogeneous media. We first analyze the random Helmholtz problem by deriving some wave-number-explicit solution estimates. We then establish a multi-modes representation of the solution as a power series of the perturbation parameter and analyze its finite modes approximations. Based on this multi-modes representation, we develop a Monte Carlo interior penalty discontinuous Galerkin (MCIP-DG) method for approximating the mode functions, which are governed by recursively defined nearly deterministic Helmholtz equations. Optimal order error estimates are derived for the method and an efficient algorithm, which is based on the LU direct solver, is also designed for efficiently implementing the proposed multi-modes MCIP-DG method. It is proved that the computational complexity of the whole algorithm is comparable to that of solving one deterministic Helmholtz problem using the LU director solver. Numerical experiments are provided to validate the theoretical results and to gauge the performance of the proposed numerical method and algorithm.

preprint2013arXiv

On Schwarz Methods for Nonsymmetric and Indefinite Problems

In this paper we introduce a new Schwarz framework and theory, based on the well-known idea of space decomposition, for nonsymmetric and indefinite linear systems arising from continuous and discontinuous Galerkin approximations of general nonsymmetric and indefinite elliptic partial differential equations. The proposed Schwarz framework and theory are presented in a variational setting in Banach spaces instead of Hilbert spaces which is the case for the well-known symmetric and positive definite (SPD) Schwarz framework and theory. Condition number estimates for the additive and hybrid Schwarz preconditioners are established. The main idea of our nonsymmetric and indefinite Schwarz framework and theory is to use weak coercivity (satisfied by the nonsymmetric and indefinite bilinear form) induced norms to replace the standard bilinear form induced norm in the SPD Schwarz framework and theory. Applications of the proposed nonsymmetric and indefinite Schwarz framework and theory. Applications of the proposed nonsymmetric and indefinite Schwarz framework to solutions of discontinuous Galerkin approximations of convection-diffusion problems are also discussed. Extensive 1-D numerical experiments are also provided to gauge the performance of the proposed Schwarz methods.