In this paper, several modifications of the quasi-interpolation operator of Scott and Zhang [30] are discussed. The modified operators are defined for non-smooth functions and are suited for application on anisotropic meshes. The anisotropy of the elements is reflected in the local stability and approximation error estimates. As an application, an example is considered where anisotropic finite element meshes are appropriate, namely the Poisson problem in domains with edges.

Ce travail a pour objet l’étude d’une méthode de « discrétisation » du Laplacien dans le problème de Poisson à deux dimensions sur un rectangle, avec des conditions aux limites de Dirichlet. Nous approchons l’opérateur Laplacien par une matrice de Toeplitz à blocs, eux-mêmes de Toeplitz, et nous établissons une formule donnant les blocs de l’inverse de cette matrice. Nous donnons ensuite un développement asymptotique de la trace de la matrice inverse, et du déterminant de la matrice de Toeplitz....

In vivo visualization of cardiovascular structures is possible using medical images. However, one has to realize that the resulting 3D geometries correspond to in vivo conditions. This entails an internal stress state to be present in the in vivo measured geometry of e.g. a blood vessel due to the presence of the blood pressure. In order to correct for this in vivo stress, this paper presents an inverse method to restore the original zero-pressure geometry of a structure, and to recover the in vivo...

We consider the problem of determining the unknown source term $f=f(x,t)$ in a space fractional diffusion equation from the measured data at the final time $u(x,T)=\psi \left(x\right)$. In this way, a methodology involving minimization of the cost functional $J\left(f\right)={\int }_0^l{(u(x,t;f){|}_{t=T}-\psi \left(x\right))}^2\mathrm{d}x$ is applied and shown that this cost functional is Fréchet differentiable and its derivative can be formulated via the solution of an adjoint problem. In addition, Lipschitz continuity of the gradient is proved. These results help us to prove the monotonicity and convergence...

In this paper, we propose a new numerical method for solving elliptic equations in unbounded regions of ${ℝ}^n$. The method is based on the mapping of a part of the domain into a bounded region. An appropriate family of weighted spaces is used for describing the growth or the decay of functions at large distances. After exposing the main ideas of the method, we analyse carefully its convergence. Some 3D computational results are displayed to demonstrate its efficiency and its high performance.

In this paper, we propose a new numerical method for solving elliptic equations in unbounded regions of ${ℝ}^n$. The method is based on the mapping of a part of the domain into a bounded region. An appropriate family of weighted spaces is used for describing the growth or the decay of functions at large distances. After exposing the main ideas of the method, we analyse carefully its convergence. Some 3D computational results are displayed to demonstrate its efficiency and its high performance.

The article is devoted to the simulation of viscous incompressible fluid flow based on solving the Navier-Stokes equations. As a numerical model we chose isogeometrical approach. Primary goal of using isogemetric analysis is to be always geometrically exact, independently of the discretization, and to avoid a time-consuming generation of meshes of computational domains. For higher Reynolds numbers, we use stabilization techniques SUPG and PSPG. All methods mentioned in the paper are demonstrated...

Estimates for the combined effect of boundary approximation and numerical integration on the approximation of (simple) eigenvalues and eigenvectors of 4th order eigenvalue problems with variable/constant coefficients in convex domains with curved boundary by an isoparametric mixed finite element method, which, in the particular case of bending problems of aniso-/ortho-/isotropic plates with variable/constant thickness, gives a simultaneous approximation to bending moment tensor field $\Psi ={\left({\psi }_{ij}\right)}_{1\le i,j\le 2}$ and displacement...

Estimates for the combined effect of boundary approximation and numerical integration on the approximation of (simple) eigenvalues and eigenvectors of 4th order eigenvalue problems with variable/constant coefficients in convex domains with curved boundary by an isoparametric mixed finite element method, which, in the particular case of bending problems of aniso-/ortho-/isotropic plates with variable/constant thickness, gives a simultaneous approximation to bending moment tensor field $\Psi ={\left({\psi }_{ij}\right)}_{1\le i,j\le 2}$ and displacement...

We consider high order finite difference approximations to the Helmholtz equation in an exterior domain. We include a simplified absorbing boundary condition to approximate the Sommerfeld radiation condition. This yields a large, but sparse, complex system, which is not self-adjoint and not positive definite. We discretize the equation with a compact fourth or sixth order accurate scheme. We solve this large system of linear equations with a Krylov subspace iterative method. Since the method converges...

We consider high order finite difference approximations to the Helmholtz equation in an exterior domain. We include a simplified absorbing boundary condition to approximate the Sommerfeld radiation condition. This yields a large, but sparse, complex system, which is not self-adjoint and not positive definite. We discretize the equation with a compact fourth or sixth order accurate scheme. We solve this large system of linear equations with a Krylov subspace iterative method. Since the method converges...

In this paper, we are concerned with a kind of Signorini transmission problem in a unbounded domain. A variational inequality is derived when discretizing this problem by coupled FEM-BEM. To solve such variational inequality, an iterative method, which can be viewed as a variant of the D-N alternative method, will be introduced. In the iterative method, the finite element part and the boundary element part can be solved independently. It will be shown that the convergence speed of this iteration...