Modeling the kinetics of a precipitation dissolution reaction occurring in a porous medium where diffusion alsotakes place leads to a system of two parabolic equations and one ordinary differential equation coupled with a stiff reaction term. This system is discretized by a finite volume scheme which is suitable for the approximation of thediscontinuous reaction term of unknown sign. Discrete solutions are shown to exist and converge towards a weak solution of the continuous problem. Uniqueness is proved under a Lipschitz condition on the equilibrium gap function. Numerical tests are shown which prove the efficiency of the scheme.

Low order edge elements are widely used for electromagnetic field problems. Higher order edge approximations are receiving increasing interest but their definition become rather complex.In this paper we propose a simple definition for Whitney edge elements of polynomial degree higher than one. We give a geometrical localization of all degrees of freedom over particular edges and provide a basis for these elements on simplicial meshes.As for Whitney edge elements of degree one, the basis is expressedonly in terms of the barycentric coordinates of the simplex.

The superconsistent collocation method, which is based on acollocation grid different from the one used to represent thesolution, has proven to be very accurate in the resolution ofvarious functional equations. Excellent results can be alsoobtained for what concerns preconditioning. Some analysis andnumerous experiments, regarding the use of finite-differencespreconditioners, for matrices arising from pseudospectralapproximations of advection-diffusion boundary value problems, arepresented and discussed, both in the case of Legendre andChebyshev representation nodes.

We introduce and analyse a mixed formulation of theMonge-Kantorovich equations, which express optimality conditions forthe mass transportation problem with cost proportional to distance.Furthermore, we introduce and analyse the finite elementapproximation of this formulation using the lowest orderRaviart-Thomas element. Finally, we present some numericalexperiments, where both the optimal transport density and theassociated Kantorovich potential are computed for a coupling problemand problems involving obstacles and regions of cheaptransportation.

In this paper, we compare a biomechanics empirical model of the heart fibrous structure to two models obtained by a non-periodic homogenization process. To this end, the two homogenized models are simplified using the small amplitude homogenization procedure of Tartar, both in conduction and in elasticity. A new small amplitude homogenization expansion formula for a mixture of anisotropic elastic materials is also derived and allows us to obtain a third simplified model.


Systems of mixed hyperbolic-elliptic conservation laws can serve as models for the evolution of a liquid-vapor fluid with possible sharp dynamical phase changes. We focus on the equations of ideal hydrodynamics in the isothermal case and introduce a thermodynamically consistent solution of the Riemann problem in one space dimension. This result is the basis for an algorithm of ghost fluid type to solve the sharp-interface model numerically. In particular the approach allows to resolve phase transitions sharply, i.e., without artificial smearing in the physically irrelevant elliptic region. Numerical experiments demonstrate the reliability of the method.