Solution of Stefan problems by fully discrete linear schemes.
Caused by the problem of unilateral contact during vibrations of satellite solar arrays, the aim of this paper is to better understand such a phenomenon. Therefore, it is studied here a simplified model composed by a beam moving between rigid obstacles. Our purpose is to describe and compare some families of fully discretized approximations and their properties, in the case of non-penetration Signorini's conditions. For this, starting from the works of Dumont and Paoli, we adapt to our beam...
Caused by the problem of unilateral contact during vibrations of satellite solar arrays, the aim of this paper is to better understand such a phenomenon. Therefore, it is studied here a simplified model composed by a beam moving between rigid obstacles. Our purpose is to describe and compare some families of fully discretized approximations and their properties, in the case of non-penetration Signorini's conditions. For this, starting from the works of Dumont and Paoli, we adapt to our beam...
We consider a family of conforming finite element schemes with piecewise polynomial space of degree in space for solving the wave equation, as a model for second order hyperbolic equations. The discretization in time is performed using the Newmark method. A new a priori estimate is proved. Thanks to this new a priori estimate, it is proved that the convergence order of the error is in the discrete norms of and , where and are the mesh size of the spatial and temporal discretization, respectively....
There are many methods and approaches to solving convection--diffusion problems. For those who want to solve such problems the situation is very confusing and it is very difficult to choose the right method. The aim of this short overview is to provide basic guidelines and to mention the common features of different methods. We place particular emphasis on the concept of linear and non-linear stabilization and its implementation within different approaches.
The paper presents the theory of the discontinuous Galerkin finite element method for the space-time discretization of a linear nonstationary convection-diffusion-reaction initial-boundary value problem. The discontinuous Galerkin method is applied separately in space and time using, in general, different nonconforming space grids on different time levels and different polynomial degrees and in space and time discretization, respectively. In the space discretization the nonsymmetric interior...
This paper is concerned with the analysis and implementation of spectral Galerkin methods for a class of Fokker-Planck equations that arises from the kinetic theory of dilute polymers. A relevant feature of the class of equations under consideration from the viewpoint of mathematical analysis and numerical approximation is the presence of an unbounded drift coefficient, involving a smooth convex potential that is equal to along the boundary of the computational domain . Using a symmetrization...
This paper is concerned with the analysis and implementation of spectral Galerkin methods for a class of Fokker-Planck equations that arises from the kinetic theory of dilute polymers. A relevant feature of the class of equations under consideration from the viewpoint of mathematical analysis and numerical approximation is the presence of an unbounded drift coefficient, involving a smooth convex potential U that is equal to +∞ along the boundary ∂D of the computational domain D. Using a symmetrization...
We consider here the Interior Penalty Discontinuous Galerkin (IPDG) discretization of the wave equation. We show how to derive the optimal penalization parameter involved in this method in the case of regular meshes. Moreover, we provide necessary stability conditions of the global scheme when IPDG is coupled with the classical Leap–Frog scheme for the time discretization. Numerical experiments illustrate the fact that these conditions are also sufficient.
This paper is concerned with the stability analysis of the space-time discontinuous Galerkin method for the solution of nonstationary, nonlinear, convection-diffusion problems. In the formulation of the numerical scheme we use the nonsymmetric, symmetric and incomplete versions of the discretization of diffusion terms and interior and boundary penalty. Then error estimates are briefly characterized. The main attention is paid to the investigation of unconditional stability of the method. Theoretical...
We show stability and consistency of the linear semi-implicit complementary volume numerical scheme for solving the regularized, in the sense of Evans and Spruck, mean curvature flow equation in the level set formulation. The numerical method is based on the finite volume methodology using the so-called complementary volumes to a finite element triangulation. The scheme gives the solution in an efficient and unconditionally stable way.
We analyze two numerical schemes of Euler type in time and C0 finite-element type with -approximation in space for solving a phase-field model of a binary alloy with thermal properties. This model is written as a highly non-linear parabolic system with three unknowns: phase-field, solute concentration and temperature, where the diffusion for the temperature and solute concentration may degenerate. The first scheme is nonlinear, unconditionally stable and convergent. The other scheme is linear...