This paper deals with the numerical approximation of mild solutions of elliptic-parabolic equations, relying on the existence results of Bénilan and Wittbold (1996). We introduce a new and simple algorithm based on Halpern's iteration for nonexpansive operators (Bauschke, 1996; Halpern, 1967; Lions, 1977), which is shown to be convergent in the degenerate case, and compare it with existing schemes (Jäger and Kačur, 1995; Kačur, 1999).
The paper contributes to the problem of finding all possible structures and waves, which may arise and preserve themselves in the open nonlinear medium, described by the mathematical model of heat structures. A new class of self-similar blow-up solutions of this model is constructed numerically and their stability is investigated. An effective and reliable numerical approach is developed and implemented for solving the nonlinear elliptic self-similar problem and the parabolic problem. This approach...
In order to investigate effects of the dynamic capillary pressure-saturation relationship used in the modelling of a flow in porous media, a one-dimensional fully implicit numerical scheme is proposed. The numerical scheme is used to simulate an experimental procedure using a measured dataset for the sand and fluid properties. Results of simulations using different models for the dynamic effect term in capillary pressure-saturation relationship are presented and discussed.
Numerical schemes are presented for a class of fourth order diffusion problems. These problems arise in lubrication theory for thin films of viscous fluids on surfaces. The equations being in general fourth order degenerate parabolic, additional singular terms of second order may occur to model effects of gravity, molecular interactions or thermocapillarity. Furthermore, we incorporate nonlinear surface tension terms. Finally, in the case of a thin film flow driven by a surface active agent (surfactant),...
In this article, we investigate numerical schemes for solving a three component Cahn-Hilliard model. The space discretization is performed by using a Galerkin formulation and the finite element method. Concerning the time discretization, the main difficulty is to write a scheme ensuring, at the discrete level, the decrease of the free energy and thus the stability of the method. We study three different schemes and prove existence and convergence theorems. Theoretical results are illustrated by...
In this article, we investigate numerical schemes for solving a three component Cahn-Hilliard model. The space discretization is performed by using a Galerkin formulation and the finite element method. Concerning the time discretization, the main difficulty is to write a scheme ensuring, at the discrete level, the decrease of the free energy and thus the stability of the method. We study three different schemes and prove existence and convergence theorems. Theoretical results are illustrated by...
This article is devoted to the numerical study of a flame ball model, derived by Joulin, which obeys to a singular integro-differential equation. The numerical scheme that we analyze here, is based upon a one step method, and we are interested in its long-time behaviour. We recover the same dynamics as in the continuous case: quenching, or stabilization of the flame, depending on heat losses, and an energy input parameter.
This article is devoted to the numerical study of a flame ball model, derived by Joulin, which obeys to a singular integro-differential equation. The numerical scheme that we analyze here, is based upon a one step method, and we are interested in its long-time behaviour. We recover the same dynamics as in the continuous case: quenching, or stabilization of the flame, depending on heat losses, and an energy input parameter.