Approximation of a Nondifferentiable Nonlinear Problem Related to MHD Equilibria.
The main goal of this article is to establish a priori and a posteriori error estimates for the numerical approximation of some non linear elliptic problems arising in glaciology. The stationary motion of a glacier is given by a non-newtonian fluid flow model which becomes, in a first two-dimensional approximation, the so-called infinite parallel sided slab model. The approximation of this model is made by a finite element method with piecewise polynomial functions of degree 1. Numerical results...
In this paper, a nonlinear problem corresponding to a simplified Oldroyd-B model without convective terms is considered. Assuming the domain to be a convex polygon, existence of a solution is proved for small relaxation times. Continuous piecewise linear finite elements together with a Galerkin Least Square (GLS) method are studied for solving this problem. Existence and a priori error estimates are established using a Newton-chord fixed point theorem, a posteriori error estimates are also derived....
In this article, we present a numerical scheme based on a finite element method in order to solve a time-dependent convection-diffusion equation problem and satisfy some conservation properties. In particular, our scheme is able to conserve the total energy for a heat equation or the total mass of a solute in a fluid for a concentration equation, even if the approximation of the velocity field is not completely divergence-free. We establish a priori errror estimates for this scheme and we give some...
In this paper, a nonlinear problem corresponding to a simplified Oldroyd-B model without convective terms is considered. Assuming the domain to be a convex polygon, existence of a solution is proved for small relaxation times. Continuous piecewise linear finite elements together with a Galerkin Least Square (GLS) method are studied for solving this problem. Existence and error estimates are established using a Newton-chord fixed point theorem, error estimates are also derived. An Elastic Viscous...
The computation of glacier movements leads to a system of nonlinear partial differential equations. The existence and uniqueness of a weak solution is established by using the calculus of variations. A discretization by the finite element method is done. The solution of the discrete problem is proved to be convergent to the exact solution. A first simple numerical algorithm is proposed and its convergence numerically studied.
The main goal of this article is to establish and error estimates for the numerical approximation of some non linear elliptic problems arising in glaciology. The stationary motion of a glacier is given by a non-Newtonian fluid flow model which becomes, in a first two-dimensional approximation, the so-called infinite parallel sided slab model. The approximation of this model is made by a finite element method with piecewise polynomial functions of degree 1. Numerical results show that the theoretical...
We formulate a finite element method for the computation of solutions to an anisotropic phase-field model for a binary alloy. Convergence is proved in the -norm. The convergence result holds for anisotropy below a certain threshold value. We present some numerical experiments verifying the theoretical results. For anisotropy below the threshold value we observe optimal order convergence, whereas in the case where the anisotropy is strong the numerical solution to the phase-field equation does not...
The motion of a three-dimensional glacier is considered. Ice is modeled as an incompressible non-Newtonian fluid. At each time step, given the shape of the glacier, a nonlinear elliptic system has to be solved in order to obtain the two components of the horizontal velocity field. Then, the shape of the glacier is updated by solving a transport equation. Finite element techniques are used to compute the velocity field and to solve the transport equation. Numerical results are compared to experiments...
Page 1