Related papers: A structure-preserving multiscale solver for parti…
Design of modern nanostructured semiconductor devices often calls for simulation tools capable of modeling arbitrarily-shaped multiscale geometries. In this work, to this end, a discontinuous Galerkin (DG) method-based framework is…
A variational method is used to derive a self-consistent macro-particle model for relativistic electromagnetic kinetic plasma simulations. Extending earlier work [E. G. Evstatiev and B. A. Shadwick, J. Comput. Phys., vol. 245, pp. 376-398,…
In high-temperature plasma physics, a strong magnetic field is usually used to confine charged particles. Therefore, for studying the classical mathematical models of the physical problems it is needed to consider the effect of external…
Numerical modeling of electromagnetic waves is an important tool for understanding the interaction of light and matter, and lies at the core of computational electromagnetics. Traditional approaches to injecting and evolving electromagnetic…
We investigate analytically and numerically the interaction of strong electromagnetic waves with unmagnetized pair plasmas. We show that the interaction is governed by a single nonlinearity parameter, $\varepsilon_{\rm p}$, defined as the…
Understanding fundamental kinetic processes is important for many problems, from plasma physics to gas dynamics. A first-principles approach to these problems requires a statistical description via the Boltzmann equation, coupled to…
The gas of the interacted electrons is usually described within Kohn-Sham approximation by the set of Poisson and Schr\"{o}dinger equations with an effective potential for the single-particle wave functions. The solution of these equations…
Numerical schemes that preserve the structure of the kinetic equations can provide stable simulation results over a long time. An electromagnetic particle-in-cell solver for the Vlasov-Maxwell equations that preserves at the discrete level…
We consider plasmonic metasurfaces constituted by an arbitrary periodic arrangement of spherical metallic nanoparticles. Each nanoparticle supports three degenerate dipolar localized surface plasmon (LSP) resonances. In the regime where the…
Electrically tunable graphene plasmons are anticipated to enable strong light-matter interactions with resonant quantum emitters. However, plasmon resonances in graphene are typically limited to infrared frequencies, below those of optical…
We introduce a new mixed discontinuous/continuous Galerkin finite element for solving the 2- and 3-dimensional wave equations and equations of incompressible flow. The element, which we refer to as P1dg-P2, uses discontinuous piecewise…
Chew, Goldberger & Low (CGL) equations are a set of hyperbolic PDEs with non-conservative products used to model the plasma flows, when the assumption of local thermodynamic equilibrium is not valid, and the pressure tensor is assumed to be…
This paper presents quasilinear theory (QLT) for classical plasma interacting with inhomogeneous turbulence. The particle Hamiltonian is kept general; for example, relativistic, electromagnetic, and gravitational effects are subsumed. A…
We discuss the development, analysis, implementation, and numerical assessment of a spectral method for the numerical simulation of the three-dimensional Vlasov-Maxwell equations. The method is based on a spectral expansion of the velocity…
This paper presents applications of weighted meshless scheme for conservation laws to the Euler equations and the equations of ideal magnetohydrodynamics. The divergence constraint of the latter is maintained to the truncation error by a…
Collisionless plasmas, mostly present in astrophysical and space environments, often require a kinetic treatment as given by the Vlasov equation. Unfortunately, the six-dimensional Vlasov equation can only be solved on very small parts of…
We develop a method for investigating the relationship between the shape of a 1-particle distribution and non-linear electrostatic oscillations in a collisionless plasma, incorporating transverse thermal motion. A general expression is…
In this work we (i) extend previous 1D studies of electromagnetic (EM) wave propagation in an over-dense plasma-metamaterial composite into two spatial dimensions and (ii) study blocking of EM waves by the composite 2D structures…
We develop a cut Discontinuous Galerkin method (cutDGM) for a diffusion-reaction equation in a bulk domain which is coupled to a corresponding equation on the boundary of the bulk domain. The bulk domain is embedded into a structured,…
Collisions are an important dissipation mechanism in plasmas. In one-dimensional modelling, a commonly used collision operator is the Lenard-Bernstein operator, or its modified energy- and momentum-conserving counterpart. When approximating…