Analysis of a Non-Linear Difference scheme in Reaction-Diffusion.
We show that a non-standard mixed finite element method proposed by Barrios and Gatica in 2007, is a higher order perturbation of the least-squares mixed finite element method. Therefore, it is also superconvergent whenever the least-squares mixed finite element method is superconvergent. Superconvergence of the latter was earlier investigated by Brandts, Chen and Yang between 2004 and 2006. Since the new method leads to a non-symmetric system matrix, its application seems however more expensive...
We study the mathematical properties of a general model of cell division structured with several internal variables. We begin with a simpler and specific model with two variables, we solve the eigenvalue problem with strong or weak assumptions, and deduce from it the long-time convergence. The main difficulty comes from natural degeneracy of birth terms that we overcome with a regularization technique. We then extend the results to the case with several parameters and recall the link between this...
In this paper we propose a time discretization of a system of two parabolic equations describing diffusion-driven atom rearrangement in crystalline matter. The equations express the balances of microforces and microenergy; the two phase fields are the order parameter and the chemical potential. The initial and boundary-value problem for the evolutionary system is known to be well posed. Convergence of the discrete scheme to the solution of the continuous problem is proved by a careful development...
We consider a time optimal control problem arisen from the optimal management of a bioreactor devoted to the treatment of eutrophicated water. We formulate this realistic problem as a state-control constrained time optimal control problem. After analyzing the state system (a complex system of coupled partial differential equations with non-smooth coefficients for advection-diffusion-reaction with Michaelis-Menten kinetics, modelling the eutrophication processes) we demonstrate the existence of,...
We consider a time optimal control problem arisen from the optimal management of a bioreactor devoted to the treatment of eutrophicated water. We formulate this realistic problem as a state-control constrained time optimal control problem. After analyzing the state system (a complex system of coupled partial differential equations with non-smooth coefficients for advection-diffusion-reaction with Michaelis-Menten kinetics, modelling the eutrophication processes) we demonstrate the existence of,...
We consider an asymptotic preserving numerical scheme initially proposed by F. Filbet and S. Jin [J. Comput. Phys. 229 (2010)] and G. Dimarco and L. Pareschi [SIAM J. Numer. Anal. 49 (2011) 2057–2077] in the context of nonlinear and stiff kinetic equations. Here, we propose a convergence analysis of such a scheme for the approximation of a system of transport equations with a nonlinear source term, for which the asymptotic limit is given by a conservation law. We investigate the convergence of the...
Compatible schemes localize degrees of freedom according to the physical nature of the underlying fields and operate a clear distinction between topological laws and closure relations. For elliptic problems, the cornerstone in the scheme design is the discrete Hodge operator linking gradients to fluxes by means of a dual mesh, while a structure-preserving discretization is employed for the gradient and divergence operators. The discrete Hodge operator is sparse, symmetric positive definite and is...
So far optimal error estimates on Bakhvalov-type meshes are only known for finite difference and finite element methods solving linear convection-diffusion problems in the one-dimensional case. We prove (almost) optimal error estimates for problems with exponential boundary layers in two dimensions.
This paper studies the gradient flow of a regularized Mumford-Shah functional proposed by Ambrosio and Tortorelli (1990, 1992) for image segmentation, and adopted by Esedoglu and Shen (2002) for image inpainting. It is shown that the gradient flow with initial data possesses a global weak solution, and it has a unique global in time strong solution, which has at most finite number of point singularities in the space-time, when the initial data are in . A family of fully discrete approximation...
This paper studies the gradient flow of a regularized Mumford-Shah functional proposed by Ambrosio and Tortorelli (1990, 1992) for image segmentation, and adopted by Esedoglu and Shen (2002) for image inpainting. It is shown that the gradient flow with L2 x L∞ initial data possesses a global weak solution, and it has a unique global in time strong solution, which has at most finite number of point singularities in the space-time, when the initial data are in H1 x H1 ∩ L∞. A family of fully...
We study the partial differential equation max{Lu − f, H(Du)} = 0 where u is the unknown function, L is a second-order elliptic operator, f is a given smooth function and H is a convex function. This is a model equation for Hamilton-Jacobi-Bellman equations arising in stochastic singular control. We establish the existence of a unique viscosity solution of the Dirichlet problem that has a Hölder continuous gradient. We also show that if H is uniformly convex, the gradient of this solution...
This paper provides new results of consistence and convergence of the lumped parameters (ODE models) toward one-dimensional (hyperbolic or parabolic) models for blood flow. Indeed, lumped parameter models (exploiting the electric circuit analogy for the circulatory system) are shown to discretize continuous 1D models at first order in space. We derive the complete set of equations useful for the blood flow networks, new schemes for electric circuit analogy, the stability criteria that guarantee...
This paper provides new results of consistence and convergence of the lumped parameters (ODE models) toward one-dimensional (hyperbolic or parabolic) models for blood flow. Indeed, lumped parameter models (exploiting the electric circuit analogy for the circulatory system) are shown to discretize continuous 1D models at first order in space. We derive the complete set of equations useful for the blood flow networks, new schemes for electric circuit analogy, the stability criteria that...