The superconvergent property of the Discontinuous Galerkin (DG) method for linear hyperbolic systems of partial differential equations in one dimension is explained by relating the DG method to a particular continuous method, whose accuracy depends in part on a local analysis, and in part on information transferred from upwind elements.

We propose here a class of numerical schemes for the approximation of weak solutions to nonlinear hyperbolic systems in nonconservative form – the notion of solution being understood in the sense of Dal Maso, LeFloch, and Murat (DLM). The proposed numerical method falls within LeFloch-Mishra's framework of schemes with well-controlled dissipation (WCD), recently introduced for dealing with small-scale dependent shocks. We design WCD schemes which are consistent with a given nonconservative system at arbitrarily high-order and then analyze their linear stability. We then investigate several nonconservative hyperbolic models arising in complex fluid dynamics, and we numerically demonstrate the convergence of our schemes toward physically meaningful weak solutions.

Plasma waves with frequencies close to the particular gyrofrequencies of the charged particles in the plasma lose energy due to cyclotron damping. We briefly discuss the gyro-resonance of low frequency plasma waves and ions particularly with regard to particle-in-cell (PiC) simulations. A setup is outlined which uses artificially excited waves in the damped regime of the wave mode's dispersion relation to track the damping of the wave's electromagnetic fields. Extracting the damping rate directly from the field data in real or Fourier space is an intricate and non-trivial task. We therefore present a simple method of obtaining the damping rate Γ from the simulation data. This method is described in detail, focusing on a step-by-step explanation of the course of actions. In a first application to a test simulation we find that the damping rates obtained from this simulation generally are in good agreement with theoretical predictions. We then compare the results of one-, two- and three-dimensional simulation setups and simulations with different physical parameter sets.

Moment models are often used for the solution of kinetic equations such as the Boltzmann equation. Unfortunately, standard models like Grad's equations are not hyperbolic and can lead to nonphysical solutions. Newly derived moment models like the Hyperbolic Moment Equations and the Quadrature-Based Moment Equations yield globally hyperbolic equations but are given in partially conservative form that cannot be written as a conservative system.
In this paper we investigate the applicability of different dedicated numerical schemes to solve the partially conservative model equations. Caused by the non-conservative type of equation we obtain differences in the numerical solutions, but due to the structure of the moment systems we show that these effects are very small for standard simulation cases. After successful identification of useful numerical settings we show a convergence study for a shock tube problem and compare the results to a discrete velocity solution. The results are in good agreement with the reference solution and we see convergence considering an increasing number of moments.

A cheapest stable nonconforming finite element method is presented for solving the incompressible flow in a square cavity without smoothing the corner singularities. The stable cheapest nonconforming finite element pair based on P₁×P₀ on rectangular meshes [29] is employed with a minimal modification of the discontinuous Dirichlet data on the top boundary, where $\widetilde{\mathcal{P}^h_0}$ is the finite element space of piecewise constant pressures with the globally one-dimensional checker-board pattern subspace eliminated. The proposed Stokes elements have the least number of degrees of freedom compared to those of known stable Stokes elements. Three accuracy indications for our elements are analyzed and numerically verified. Also, various numerous computational results obtained by using our proposed element show excellent accuracy.

In this study an explicit Finite Difference Method (FDM) based scheme is developed to solve the Maxwell's equations in time domain for a lossless medium. This manuscript focuses on two unique aspects – the three dimensional time-accurate discretization of the hyperbolic system of Maxwell equations in three-point non-staggered grid stencil and it's application to parallel computing through the use of Graphics Processing Units (GPU). The proposed temporal scheme is symplectic, thus permitting conservation of all Hamiltonians in the Maxwell equation. Moreover, to enable accurate predictions over large time frames, a phase velocity preserving scheme is developed for treatment of the spatial derivative terms. As a result, the chosen time increment and grid spacing can be optimally coupled. An additional theoretical investigation into this pairing is also shown. Finally, the application of the proposed scheme to parallel computing using one Nvidia K20 Tesla GPU card is demonstrated. For the benchmarks performed, the parallel speedup when compared to a single core of an Intel i7-4820K CPU is approximately 190x.

We present a high-order discontinuous Galerkin (DG) method for the time domain Maxwell's equations in three-dimensional heterogeneous media. New hierarchical orthonormal basis functions on unstructured tetrahedral meshes are used for spatial discretization while Runge-Kutta methods for time discretization. A uniaxial perfectly matched layer (UPML) is employed to terminate the computational domain. Exponential convergence with respect to the order of the basis functions is observed and large parallel speedup is obtained for a plane-wave scattering model. The rapid decay of the out-going wave in the UPML is shown in a dipole radiation simulation. Moreover, the low frequency electromagnetic fields excited by a horizontal electric dipole (HED) and a vertical magnetic dipole (VMD) over a layered conductive half-space and a high frequency ground penetrating radar (GPR) detection for an underground structure are investigated, showing the high accuracy and broadband simulation capability of the proposed method.

Pressure-correction projection finite element methods (FEMs) are proposed to solve nonstationary natural convection problems in this paper. The first-order and second-order backward difference formulas are applied for time derivative, the stability analysis and error estimates of the semi-discrete schemes are presented using energy method. Compared with characteristic variational multiscale FEM, pressure-correction projection FEMs are more efficient and unconditionally energy stable. Ample numerical results are presented to demonstrate the effectiveness of the pressure-correction projection FEMs for solving these problems.

The evolution of material points interacting via gravitational force in 3D space was investigated. At initial moment points with masses of 2.48 Sun masses are randomly distributed inside a cube with an edge of 5 light-years. The modeling was conducted at different initial distributions of velocities and different ratios between potential and kinetic energy of the points. As a result of modeling the time dependence of velocity distribution function of points was obtained. Dependence of particles fraction which had evaporated from initial cluster on time for different initial conditions is obtained. In particular, it was obtained that the fraction of evaporated particles varies between 0,45 and 0,63.
Mutual diffusion of two classes of particles at different initial conditions in the case when at initial moment of time both classes of particles occupy equal parts of cube was investigated.
The maximum Lyapunov exponent of the system with different initial conditions was calculated. The obtained value weakly depends on the ratio between initial kinetic and potential energies and amounts approximately 10⁻⁵. Corresponding time of the particle trajectories divergence turned out to be 40-50 thousand years.

We study a novel class of numerical integrators, the adapted nested force-gradient schemes, used within the molecular dynamics step of the Hybrid Monte Carlo (HMC) algorithm. We test these methods in the Schwinger model on the lattice, a well known benchmark problem. We derive the analytical basis of nested force-gradient type methods and demonstrate the advantage of the proposed approach, namely reduced computational costs compared with other numerical integration schemes in HMC. 

The method of analytic continuation in the coupling constant in combination with the use of statistical Padé approximation is applied to the determination of complex S-matrix poles, i.e. to the determination of resonance energy and widths. These parameters are of vital importance in many physical, chemical and biological processes. It is shown that an alternative to the method of analytic continuation in the coupling constant exists which in principle makes it possible to locate several resonances at once, in contrast to the original method which yields parameters of only one resonance. In addition, the new approach appears to be less sensitive to the choice of the perturbation interaction used for the analytical continuation than the original method. In this paper both approaches are compared and tested for model analytic separable potential. It is shown that the new variant of the method of analytic continuation in the coupling constant is more robust and efficient than the original method and yields reasonable results even for data of limited accuracy.

In this paper, a novel second-order two-scale (SOTS) computational method is developed for nonlinear dynamic thermo-mechanical problems of composites with cylindrical periodicity. The nonlinearities of these multi-scale problems were caused by the temperature-dependent properties of the composites. Firstly, the formal SOTS solutions for these problems are constructed by the multiscale asymptotic analysis. Then we theoretically explain the importance of the SOTS solutions by the error analysis in the pointwise sense. In addition, a SOTS numerical algorithm is proposed in detail to effectively solve these problems. Finally, some numerical examples verify the feasibility and effectiveness of the SOTS numerical algorithm we proposed.