Solving of heat shock on a hybrid system
The numerical approximation of parametric partial differential equations is a computational challenge, in particular when the number of involved parameter is large. This paper considers a model class of second order, linear, parametric, elliptic PDEs on a bounded domain D with diffusion coefficients depending on the parameters in an affine manner. For such models, it was shown in [9, 10] that under very weak assumptions on the diffusion coefficients, the entire family of solutions to such equations...
We are interested in the discretization of the heat equation with a diffusion coefficient depending on the space and time variables. The discretization relies on a spectral element method with respect to the space variables and Euler's implicit scheme with respect to the time variable. A detailed numerical analysis leads to optimal a priori error estimates.
We consider the Stokes problem provided with non standard boundary conditions which involve the normal component of the velocity and the tangential components of the vorticity. We write a variational formulation of this problem with three independent unknowns: the vorticity, the velocity and the pressure. Next we propose a discretization by spectral element methods which relies on this formulation. A detailed numerical analysis leads to optimal error estimates for the three unknowns and numerical...
We consider the Stokes problem provided with non standard boundary conditions which involve the normal component of the velocity and the tangential components of the vorticity. We write a variational formulation of this problem with three independent unknowns: the vorticity, the velocity and the pressure. Next we propose a discretization by spectral element methods which relies on this formulation. A detailed numerical analysis leads to optimal error estimates for the three unknowns and numerical...
We study spectral discretizations for singular perturbation problems. A special technique of stabilization for the spectral method is proposed. Boundary layer problems are accurately solved by a domain decomposition method. An effective iterative method for the solution of spectral systems is proposed. Suitable components for a multigrid method are presented.
Strong convergence estimates for pseudospectral methods applied to ordinary boundary value problems are derived. The results are also used for a convergence analysis of the Schwarz algorithm (a special domain decomposition technique). Different types of nodes (Chebyshev, Legendre nodes) are examined and compared.
We present a parallel preconditioning method for the iterative solution of the time-harmonic elastic wave equation which makes use of higher-order spectral elements to reduce pollution error. In particular, the method leverages perfectly matched layer boundary conditions to efficiently approximate the Schur complement matrices of a block LDLT factorization. Both sequential and parallel versions of the algorithm are discussed and results for large-scale problems from exploration geophysics are presented....
The author investigates time-periodic solutions of the quasilinear beam equation with the help of accelerated convergence methods. Using the Newton iteration scheme, the problem is approximated by a sequence of linear equations solved via the Galerkin method. The derivatiove loss inherent to this kind of problems is compensated by taking advantage of smoothing operators.
The method of transfer of boundary conditions yields a universal frame into which most methods for solving boundary value problems for ordinary differential equations can be included. The purpose of this paper is to show a possibility to extend the idea of transfer of conditions to a particular twodimensional problem.
We consider a family of quadrilateral or hexahedral mixed hp-finite elements for an incompressible flow problem with Qr-elements for the velocity and discontinuous -elements for the pressure where the order r can vary from element to element between 2 and an arbitrary bound. For multilevel adaptive grids with hanging nodes and a sufficiently small mesh size, we prove the inf-sup condition uniformly with respect to the mesh size and the polynomial degree.