We present a discontinuous Galerkin scheme for the numerical approximation of the one-dimensional periodic Vlasov-Poisson equation. The scheme is based on a Galerkin-characteristics method in which the distribution function is projected onto a space of discontinuous functions. We present comparisons with a semi-Lagrangian method to emphasize the good behavior of this scheme when applied to Vlasov-Poisson test cases.
In this paper we prove the discrete compactness property for a discontinuous Galerkin approximation of Maxwell's system on quite general tetrahedral meshes. As a consequence, a discrete Friedrichs inequality is obtained and the convergence of the discrete eigenvalues to the continuous ones is deduced using the theory of collectively compact operators. Some numerical experiments confirm the theoretical predictions.
Our studies are motivated by a desire to model long-time simulations of possible scenarios for a waste disposal. Numerical methods are developed for solving the arising systems of convection-diffusion-dispersion-reaction equations, and the received results of several discretization methods are presented. We concentrate on linear reaction systems, which can be solved analytically. In the numerical methods, we use large time-steps to achieve long simulation times of about 10 000 years. We propose...
Two new time-dependent versions of div-curl results in a bounded domain are presented. We study a limit of the product , where the sequences and belong to . In Theorem 2.1 we assume that is bounded in the -norm and is controlled in the -norm. In Theorem 2.2 we suppose that is bounded in the -norm and is controlled in the -norm. The time derivative of is bounded in both cases in the norm of . The convergence (in the sense of distributions) of to the product of weak limits...
This article presents some results of numerical tests of solving the two-dimensional non-linear unsteady viscous Burgers equation. We have compared the known convergence and parallel performance properties of the additive Schwarz domain decomposition method with or without a coarse grid for the model Poisson problem with those obtained by experiments for the Burgers problem.
A mixed finite element method for the Navier–Stokes equations is introduced in which the stress is a primary variable. The variational formulation retains the mathematical structure of the Navier–Stokes equations and the classical theory extends naturally to this setting. Finite element spaces satisfying the associated inf–sup conditions are developed.
We derive a biomembrane model consisting of a fluid enclosed by a lipid membrane. The membrane is characterized by its Canham-Helfrich energy (Willmore energy with area constraint) and acts as a boundary force on the Navier-Stokes system modeling an incompressible fluid. We give a concise description of the model and of the associated numerical scheme. We provide numerical simulations with emphasis on the comparisons between different types of flow:...
The Maxwell type elastic-plastic solids are characterized by decaying the absolute values of the principal components of the deviatoric part of the stress tensor during the plastic relaxation step. We propose a mathematical formulation of such a model which is compatible with the von Mises criterion of plasticity. Numerical examples show the ability of the model to deal with complex physical phenomena.
We introduce a new way of computation of time dependent partial differential equations using hybrid method FEM in space and FDM in time domain and explicit computational scheme. The key idea is quick transformation of standard basis functions into new simple basis functions. This new way is used for better computational efficiency. We explain this way of computation on an example of elastodynamic equation using quadrilateral elements. However, the method can be used for more types of elements and...
We introduce a new efficient way of computation of partial differential equations using a hybrid method composed from FEM in space and FDM in time domain. The overall computational scheme is explicit in time. The key idea of the suggested way is based on a transformation of standard basis functions into new basis functions. The results of this matrix transformation are e-invariants (effective invariants) with such suitable properties which save the number of arithmetical operations needed for a...
In this paper, we extend the reduced-basis approximations developed earlier for linear elliptic and parabolic partial differential equations with affine parameter dependence to problems involving (a) nonaffine dependence on the parameter, and (b) nonlinear dependence on the field variable. The method replaces the nonaffine and nonlinear terms with a coefficient function approximation which then permits an efficient offline-online computational decomposition. We first review the coefficient function...
The use of parallel computers makes it feasible to simulate elastic waves throughout large heterogeneous structures, and new domain decomposition methods can be used to increase their efficiency and decrease the computing time spent in the simulation. In this paper we introduce a simple parallel algorithm for the propagation of elastic waves in complex heterogeneous media after a finite element discretization. This method performs more efficiently than classic domain decomposition techniques based...