Source author record

Jay Gopalakrishnan

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

14 published item(s)

preprint2026arXiv

Generalizing Riemann curvature to Regge metrics

In this paper, we propose a generalization of the Riemann curvature tensor on manifolds (of dimension two or higher) endowed with a Regge metric. Specifically, while all components of the metric tensor are assumed to be smooth within elements of a triangulation of the manifold, they need not be smooth across element interfaces, where only continuity of the tangential components is assumed. While linear derivatives of the metric can be generalized to Schwartz distributions, similarly generalizing the classical Riemann curvature tensor, a nonlinear second-order derivative of the metric, requires more care. We propose a generalization that combines the classical angle defect and jumps in the second fundamental form across element interfaces, and argue its correctness. Specifically, if a piecewise smooth metric approximates a globally smooth metric, then our generalized Riemann curvature tensor approximates the classical Riemann curvature tensor associated with the latter. Moreover, we show that if the metric approximation converges at some rate in a mesh-dependent norm equivalent to the $L^2$ norm, then the curvature approximation converges in the negative Sobolev space $H^{-2}$, the dual space of $H^2_0$, at the same rate, under additional assumptions. By appropriate contractions of the generalized Riemann curvature tensor, this work also provides generalizations of scalar curvature, the Ricci curvature tensor, and the Einstein tensor in any dimension.

preprint2020arXiv

A discrete elasticity complex on three-dimensional Alfeld splits

We construct conforming finite element elasticity complexes on the Alfeld splits of tetrahedra. The complex consists of vector fields and symmetric tensor fields, interlinked via the linearized deformation operator, the linearized curvature operator, and the divergence operator, respectively. The construction is based on an algebraic machinery that derives the elasticity complex from de~Rham complexes, and smoother finite element differential forms.

preprint2020arXiv

Structure aware Runge-Kutta time stepping for spacetime tents

We introduce a new class of Runge-Kutta type methods suitable for time stepping to propagate hyperbolic solutions within tent-shaped spacetime regions. Unlike standard Runge-Kutta methods, the new methods yield expected convergence properties when standard high order spatial (discontinuous Galerkin) discretizations are used. After presenting a derivation of nonstandard order conditions for these methods, we show numerical examples of nonlinear hyperbolic systems to demonstrate the optimal convergence rates. We also report on the discrete stability properties of these methods applied to linear hyperbolic equations.

preprint2019arXiv

Simulation of Optical Fiber Amplifier Gain Using Equivalent Short Fibers

Electromagnetic wave propagation in optical fiber amplifiers obeys Maxwell equations. Using coupled mode theory, the full Maxwell system within an optical fiber amplifier is reduced to a simpler model. The simpler model is made more efficient through a new scale model, referred to as an equivalent short fiber, which captures some of the essential characteristics of a longer fiber. The equivalent short fiber can be viewed as a fiber made using artificial (unphysical) material properties that in some sense compensates for its reduced length. The computations can be accelerated by a factor approximately equal to the ratio of the original length to the reduced length of the equivalent fiber. Computations using models of two commercially available fibers -- one doped with ytterbium, and the other with thulium -- show the practical utility of the concept. Extensive numerical studies are conducted to assess when the equivalent short fiber model is useful and when it is not.

preprint2018arXiv

The DPG-star method

This article introduces the DPG-star (from now on, denoted DPG$^*$) finite element method. It is a method that is in some sense dual to the discontinuous Petrov-Galerkin (DPG) method. The DPG methodology can be viewed as a means to solve an overdetermined discretization of a boundary value problem. In the same vein, the DPG$^*$ methodology is a means to solve an underdetermined discretization. These two viewpoints are developed by embedding the same operator equation into two different saddle-point problems. The analyses of the two problems have many common elements. Comparison to other methods in the literature round out the newly garnered perspective. Notably, DPG$^*$ and DPG methods can be seen as generalizations of $\mathcal{L}\mathcal{L}^\ast$ and least-squares methods, respectively. A priori error analysis and a posteriori error control for the DPG$^*$ method are considered in detail. Reports of several numerical experiments are provided which demonstrate the essential features of the new method. A notable difference between the results from the DPG$^*$ and DPG analyses is that the convergence rates of the former are limited by the regularity of an extraneous Lagrange multiplier variable.

preprint2015arXiv

A tent pitching scheme motivated by Friedrichs theory

Certain Friedrichs systems can be posed on Hilbert spaces normed with a graph norm. Functions in such spaces arising from advective problems are found to have traces with a weak continuity property at points where the inflow and outflow boundaries meet. Motivated by this continuity property, an explicit space-time finite element scheme of the tent pitching type, with spaces that conform to the continuity property, is designed. Numerical results for a model one-dimensional wave propagation problem are presented.

preprint2013arXiv

Dispersive and dissipative errors in the DPG method with scaled norms for Helmholtz equation

We consider the discontinuous Petrov-Galerkin (DPG) method, wher the test space is normed by a modified graph norm. The modificatio scales one of the terms in the graph norm by an arbitrary positive scaling parameter. Studying the application of the method to the Helmholtz equation, we find that better results are obtained, under some circumstances, as the scaling parameter approaches a limiting value. We perform a dispersion analysis on the multiple interacting stencils that form the DPG method. The analysis shows that the discrete wavenumbers of the method are complex, explaining the numerically observed artificial dissipation in the computed wave approximations. Since the DPG method is a nonstandard least-squares Galerkin method, we compare its performance with a standard least-squares method.

preprint2012arXiv

A locking free hp DPG method for linear elasticity with symmetric stresses

We present two new methods for linear elasticity with simultaneously yield stress and displacement approximations of optimal accuracy in both the mesh size h and polynomial degree p. This is achieved within the recently developed discontinuous Petrov-Galerkin (DPG) framework. In this framework, both the stress and the displacement approximations are discontinuous across element interfaces. We study locking-free convergence properties and the interrelationships between the two DPG methods.

preprint2012arXiv

An analysis of the practical DPG method

In this work we give a complete error analysis of the Discontinuous Petrov Galerkin (DPG) method, accounting for all the approximations made in its practical implementation. Specifically, we consider the DPG method that uses a trial space consisting of polynomials of degree $p$ on each mesh element. Earlier works showed that there is a "trial-to-test" operator $T$, which when applied to the trial space, defines a test space that guarantees stability. In DPG formulations, this operator $T$ is local: it can be applied element-by-element. However, an infinite dimensional problem on each mesh element needed to be solved to apply $T$. In practical computations, $T$ is approximated using polynomials of some degree $r > p$ on each mesh element. We show that this approximation maintains optimal convergence rates, provided that $r\ge p+N$, where $N$ is the space dimension (two or more), for the Laplace equation. We also prove a similar result for the DPG method for linear elasticity. Remarks on the conditioning of the stiffness matrix in DPG methods are also included.

preprint2012arXiv

Convergence analysis of a multigrid algorithm for the acoustic single layer equation

We present and analyze a multigrid algorithm for the acoustic single layer equation in two dimensions. The boundary element formulation of the equation is based on piecewise constant test functions and we make use of a weak inner product in the multigrid scheme as proposed in \cite{BLP94}. A full error analysis of the algorithm is presented. We also conduct a numerical study of the effect of the weak inner product on the oscillatory behavior of the eigenfunctions for the Laplace single layer operator.

preprint2012arXiv

Instability in a generalized Keller-Segel model

We present a generalized Keller-Segel model where an arbitrary number of chemical compounds react, some of which are produced by a species, and one of which is a chemoattractant for the species. To investigate the stability of homogeneous stationary states of this generalized model, we consider the eigenvalues of a linearized system. We are able to reduce this infinite dimensional eigenproblem to a parametrized finite dimensional eigenproblem. By matrix theoretic tools, we then provide easily verifiable sufficient conditions for destabilizing the homogeneous stationary states. In particular, one of the sufficient conditions is that the chemotactic feedback is sufficiently strong. Although this mechanism was already known to exist in the original Keller-Segel model, here we show that it is more generally applicable by significantly enlarging the class of models exhibiting this instability phenomenon which may lead to pattern formation.

preprint2011arXiv

Mixed finite element approximation of the vector Laplacian with Dirichlet boundary conditions

We consider the finite element solution of the vector Laplace equation on a domain in two dimensions. For various choices of boundary conditions, it is known that a mixed finite element method, in which the rotation of the solution is introduced as a second unknown, is advantageous, and appropriate choices of mixed finite element spaces lead to a stable, optimally convergent discretization. However, the theory that leads to these conclusions does not apply to the case of Dirichlet boundary conditions, in which both components of the solution vanish on the boundary. We show, by computational example, that indeed such mixed finite elements do not perform optimally in this case, and we analyze the suboptimal convergence that does occur. As we indicate, these results have implications for the solution of the biharmonic equation and of the Stokes equations using a mixed formulation involving the vorticity.

preprint2011arXiv

Partial expansion of a Lipschitz domain and some applications

We show that a Lipschitz domain can be expanded solely near a part of its boundary, assuming that the part is enclosed by a piecewise C1 curve. The expanded domain as well as the extended part are both Lipschitz. We apply this result to prove a regular decomposition of standard vector Sobolev spaces with vanishing traces only on part of the boundary. Another application in the construction of low-regularity projectors into finite element spaces with partial boundary conditions is also indicated.