Related papers: Staggered grid leap-frog scheme for the (2+1)D Dir…
The usual treatment of a (first order) classical field theory such as electromagnetism has a little drawback: It has a primary constraint submanifold that arise from the fact that the dynamics is governed by the antisymmetric part of the…
Mathematical models for flow and reactive transport in porous media often involve non-linear, degenerate parabolic equations. Their solutions have low regularity, and therefore lower order schemes are used for the numerical approximation.…
We propose finite-volume schemes for the Cahn-Hilliard equation which unconditionally and discretely preserve the boundedness of the phase field and the dissipation of the free energy. Our numerical framework is applicable to a variety of…
We show that spin-flip probabilities emerge in the relativistic regime for scalar potentials, absent in the standard Dirac representation. We examine 1D scattering for the Dirac equation employing an alternate matrix representation…
A methodology on making the variational principle well-posed in degenerate systems is constructed. In the systems including higher-order time derivative terms being compatible with Newtonian dynamics, we show that a set of position…
The analyses of interior penalty discontinuous Galerkin methods of any order k for solving elliptic and parabolic problems with Dirac line sources are presented. For the steady state case, we prove convergence of the method by deriving a…
Traditional representations of graphs and their duals suggest the requirement that the dual vertices be placed inside their corresponding primal faces, and the edges of the dual graph cross only their corresponding primal edges. We consider…
Interest on 2 + 1 dimensional electron systems has increased considerably after the realization of novel properties of graphene sheets, in which the behaviour of electrons is effectively described by relativistic equations. Having this fact…
How to develop efficient numerical schemes while preserving the energy stability at the discrete level is a challenging issue for the three component Cahn-Hilliard phase-field model. In this paper, we develop first and second order temporal…
Non-commutative electrodynamics obtained through the Seiberg-Witten map ceases to have equivalent action-level and equation-level realizations once fixed external currents are introduced, and in the action-level construction associated with…
We present an extended range of stable flux reconstruction (FR) methods on triangles through the development and application of the summation-by-parts framework in two-dimensions. This extended range of stable schemes is then shown to…
We present an approach to solving problems in micromechanics that is amenable to massively parallel calculations through the use of graphical processing units and other accelerators. The problems lead to nonlinear differential equations…
It is well-known that the tight-binding Hamiltonian of graphene describes the low-energy excitations that appear to be massless chiral Dirac fermions. Thus, in the continuum limit one can analyze the crystal properties using the formalism…
The shifted-i expansion technique (SLET) is extended to solve for Dirac particle trapped in spherically symmetric scalar and/or 4-vector potentials. A parameter {\lambda}=0,1 is introduced in such a way that one can obtain the Klein-Gordon…
The Dirac equation for a massive spin-1/2 field in a central potential V in three dimensions is studied without fixing a priori the functional form of V. The second-order equations for the radial parts of the spinor wave function are shown…
The aim of this paper is the derivation of structure preserving schemes for the solution of the EPDiff equation, with particular emphasis on the two dimensional case. We develop three different schemes based on the Discrete Variational…
In this paper, we develop sparse grid central discontinuous Galerkin (CDG) scheme for linear hyperbolic systems with variable coefficients in high dimensions. The scheme combines the CDG framework with the sparse grid approach, with the aim…
The stationary Dirac equation $(p\cdot\sigma)\psi=E\psi$, confined to a two-dimensional (2D) region, supports states propagating along the boundary and decaying exponentially away from the boundary. These edge states appear on the 2D…
A boundary integral equation formulation is presented for the electromagnetic transmission problem where an incident electromagnetic wave is scattered from a bounded dielectric object. The formulation provides unique solutions for all…
In this paper, we present a systematic framework to derive a Lagrangian scheme for porous medium type generalized diffusion equations by employing a discrete energetic variational approach. Such discrete energetic variational approaches are…