112 research outputs found
Structures and waves in a nonlinear heat-conducting medium
The paper is an overview of the main contributions of a Bulgarian team of
researchers to the problem of finding the possible structures and waves in the
open nonlinear heat conducting medium, described by a reaction-diffusion
equation. Being posed and actively worked out by the Russian school of A. A.
Samarskii and S.P. Kurdyumov since the seventies of the last century, this
problem still contains open and challenging questions.Comment: 23 pages, 13 figures, the final publication will appear in Springer
Proceedings in Mathematics and Statistics, Numerical Methods for PDEs:
Theory, Algorithms and their Application
A linear domain decomposition method for partially saturated flow in porous media
The Richards equation is a nonlinear parabolic equation that is commonly used
for modelling saturated/unsaturated flow in porous media. We assume that the
medium occupies a bounded Lipschitz domain partitioned into two disjoint
subdomains separated by a fixed interface . This leads to two problems
defined on the subdomains which are coupled through conditions expressing flux
and pressure continuity at . After an Euler implicit discretisation of
the resulting nonlinear subproblems a linear iterative (-type) domain
decomposition scheme is proposed. The convergence of the scheme is proved
rigorously. In the last part we present numerical results that are in line with
the theoretical finding, in particular the unconditional convergence of the
scheme. We further compare the scheme to other approaches not making use of a
domain decomposition. Namely, we compare to a Newton and a Picard scheme. We
show that the proposed scheme is more stable than the Newton scheme while
remaining comparable in computational time, even if no parallelisation is being
adopted. Finally we present a parametric study that can be used to optimize the
proposed scheme.Comment: 34 pages, 13 figures, 7 table
Multidimensional Conservation Laws: Overview, Problems, and Perspective
Some of recent important developments are overviewed, several longstanding
open problems are discussed, and a perspective is presented for the
mathematical theory of multidimensional conservation laws. Some basic features
and phenomena of multidimensional hyperbolic conservation laws are revealed,
and some samples of multidimensional systems/models and related important
problems are presented and analyzed with emphasis on the prototypes that have
been solved or may be expected to be solved rigorously at least for some cases.
In particular, multidimensional steady supersonic problems and transonic
problems, shock reflection-diffraction problems, and related effective
nonlinear approaches are analyzed. A theory of divergence-measure vector fields
and related analytical frameworks for the analysis of entropy solutions are
discussed.Comment: 43 pages, 3 figure
An a posteriori error estimate for vertex-centered finite volume discretizations of immiscible incompressible two-phase flow
International audienceIn this paper we derive an a posteriori error estimate for the numerical approximation of the solution of a system modeling the flow of two incompressible and immiscible fluids in a porous medium. We take into account the capillary pressure, which leads to a coupled system of two equations: parabolic and elliptic. The parabolic equation may become degenerate, i.e., the nonlinear diffusion coefficient may vanish over regions that are not known a priori. We first show that, under appropriate assumptions, the energy-type-norm differences between the exact and the approximate nonwetting phase saturations, the global pressures, and the Kirchhoff transforms of the nonwetting phase saturations can be bounded by the dual norm of the residuals. We then bound the dual norm of the residuals by fully computable a posteriori estimators. Our analysis covers a large class of conforming, vertex-centered finite volume-type discretizations with fully implicit time stepping. As an example, we focus here on two approaches: a ''mathematical'' scheme derived from the weak formulation, and a phase-by-phase upstream weighting ''engineering'' scheme. Finally, we show how the different error components, namely the space discretization error, the time discretization error, the linearization error, the algebraic solver error, and the quadrature error can be distinguished and used for making the calculations efficient
Geometric partial differential equations: Theory, numerics and applications
This workshop concentrated on partial differential equations involving stationary and evolving surfaces in which geometric quantities play a major role. Mutual interest in this emerging field stimulated the interaction between analysis, numerical solution, and applications
An a posterior error estimate for vertex-centered finite volume discretizations of immiscible incompressible two-phase flow
In this paper we derive an a posteriori error estimate for the numerical approximation of the solution of a system modeling the flow of two incompressible and immiscible fluids in a porous medium. We take into account the capillary pressure, which leads to a coupled system of two equations: parabolic and elliptic. The parabolic equation may become degenerate, i.e., the nonlinear diffusion coefficient may vanish over regions that are not known a priori. We first show that, under appropriate assumptions, the energy-type-norm differences between the exact and the approximate nonwetting phase saturations, the global pressures, and the Kirchhoff transforms of the nonwetting phase saturations can be bounded by the dual norm of the residuals. We then bound the dual norm of the residuals by fully computable a posteriori estimators. Our analysis covers a large class of conforming, vertex-centered finite volume-type discretizations with fully implicit time stepping. As an example, we focus here on two approaches: a "mathematical" scheme derived from the weak formulation, and a phase-by-phase upstream weighting "engineering" scheme. Finally, we show how the different error components, namely the space discretization error, the time discretization error, the linearization error, the algebraic solver error, and the quadrature error can be distinguished and used for making the calculations efficient
Numerical analysis of a robust free energy diminishing Finite Volume scheme for parabolic equations with gradient structure
We present a numerical method for approximating the solutions of degenerate
parabolic equations with a formal gradient flow structure. The numerical method
we propose preserves at the discrete level the formal gradient flow structure,
allowing the use of some nonlinear test functions in the analysis. The
existence of a solution to and the convergence of the scheme are proved under
very general assumptions on the continuous problem (nonlinearities, anisotropy,
heterogeneity) and on the mesh. Moreover, we provide numerical evidences of the
efficiency and of the robustness of our approach
- …