The aim of this work is to compare a new uncoupled solver for the cardiac Bidomain model with a usual coupled solver. The Bidomain model describes the bioelectric activity of the cardiac tissue and consists of a system of a non-linear parabolic reaction-diffusion partial differential equation (PDE) and an elliptic linear PDE. This system models at macroscopic level the evolution of the transmembrane and extracellular electric potentials of the anisotropic cardiac tissue. The evolution equation is coupled through the non-linear reaction term with a stiff system of ordinary differential equations (ODEs), the so-called membrane model, describing the ionic currents through the cellular membrane. A novel uncoupled solver for the Bidomain system is here introduced, based on solving twice the parabolic PDE and once the elliptic PDE at each time step, and it is compared with a usual coupled solver. Three-dimensional numerical tests have been performed in order to show that the proposed uncoupled method has the same accuracy of the coupled strategy. Parallel numerical tests on structured meshes have also shown that the uncoupled technique is as scalable as the coupled one. Moreover, the conjugate gradient method preconditioned by Multilevel Hybrid Schwarz preconditioners converges faster for the linear systems deriving from the uncoupled method than from the coupled one. Finally, in all parallel numerical tests considered, the uncoupled technique proposed is always about two or three times faster than the coupled approach.

One of the current debate about simulating the electrical activity in the heart is thefollowing: Using a realistic anatomical setting, i.e. realisticgeometries, fibres orientations, etc., is it enough to use a simplified 2-variablephenomenological model to reproduce the main characteristics of the cardiac actionpotential propagation, and in what sense is it sufficient? Using a combination ofdimensional and asymptotic analysis, together with the well-known Mitchell − Schaeffermodel, it is shown that it is possible to accurately control (at least locally) thesolution while spatial propagation is involved. In particular, we reduce the set ofparameters by writing the bidomain model in a new nondimensional form. The parameters ofthe bidomain model with Mitchell − Schaeffer ion kinetics are then set and shown to be inone-to-one relation with the main characteristics of the four phases of a propagatedaction potential. Explicit relations are derived using a combination of asymptotic methodsand ansatz. These relations are tested against numerical results. We illustrate how theserelations can be used to recover the time/space scales and speed of the action potentialin various regions of the heart.

This paper is devoted to the general problem of reconstructing the cost from theobservation of trajectories, in a problem of optimal control. It is motivated by thefollowing applied problem, concerning HALE drones: one would like them to decide bythemselves for their trajectories, and to behave at least as a good human pilot. Thisapplied question is very similar to the problem of determining what is minimized in humanlocomotion. These starting points are the reasons for the particular classes of controlsystems and of costs under consideration. To summarize, our conclusion is that in general,inside these classes, three experiments visiting the same values of the control are neededto reconstruct the cost, and two experiments are in general not enough. The method isconstructive.

The proof of these results is mostly based upon the Thom’s transversality theory.

This study is partly supported by FUI AAP9 project SHARE, and by ANR Project GCM, program“blanche”, project number NT09-504490.

In this paper, we consider the approximation of second order evolution equations. It is well known that the approximated system by finite element or finite difference is not uniformly exponentially or polynomially stable with respect to the discretization parameter, even if the continuous system has this property. Our goal is to damp the spurious high frequency modes by introducing numerical viscosity terms in the approximation scheme. With these viscosity terms, we show the exponential or polynomial decay of the discrete scheme when the continuous problem has such a decay and when the spectrum of the spatial operator associated with the undamped problem satisfies the generalized gap condition. By using the Trotter–Kato Theorem, we further show the convergence of the discrete solution to the continuous one. Some illustrative examples are also presented.

We address in this article the computation of the convex solutions of the Dirichletproblem for the real elliptic Monge − Ampère equation for general convex domains in twodimensions. The method we discuss combines a least-squares formulation with a relaxationmethod. This approach leads to a sequence of Poisson − Dirichlet problems and anothersequence of low dimensional algebraic eigenvalue problems of a new type. Mixed finiteelement approximations with a smoothing procedure are used for the computer implementationof our least-squares/relaxation methodology. Domains with curved boundaries are easilyaccommodated. Numerical experiments show the convergence of the computed solutions totheir continuous counterparts when such solutions exist. On the other hand, when classicalsolutions do not exist, our methodology produces solutions in a least-squares sense.

The partially observed optimal control problem is considered for forward-backward doubly stochastic systems with controls entering into the diffusion and the observation. The maximum principle is proven for the partially observable optimal control problems. A probabilistic approach is used, and the adjoint processes are characterized as solutions of related forward-backward doubly stochastic differential equations in finite-dimensional spaces. Then, our theoretical result is applied to study a partially-observed linear-quadratic optimal control problem for a fully coupled forward-backward doubly stochastic system.

We establish some new results about the Γ-limit, with respect to the L1-topology, of two different (but related) phase-field approximations \hbox{$\{\mathcal E_\eps\}_\eps,\,\{\widetilde{\mathcalE}_\eps\}_\eps$}$\mathrm{\{}{\mathrm{ℰ}}_{\mathit{\epsilon }}{\mathrm{\}}}_{\mathit{\epsilon }}\mathit{,}\hspace{0.17em}\mathrm{\{}\begin{array}{}{}_{\mathrm{􏽥}}\\ {\mathrm{ℰ}}_{\mathit{\epsilon }}\end{array}{\mathrm{\}}}_{\mathit{\epsilon }}$ of the so-called Euler’s Elastica Bending Energy for curves in the plane. In particular we characterize the Γ-limit as ε → 0 of ℰε, and show that in general the Γ-limits of ℰε and \hbox{$\widetilde{\mathcal E}_\eps$}$\begin{array}{}{\mathrm{􏽥}}_{}\\ {\mathrm{ℰ}}_{\mathit{\epsilon }}\end{array}$ do not coincide on indicator functions of sets with non-smooth boundary. More precisely we show that the domain of the Γ-limit of \hbox{$\widetilde{\mathcal E}_\eps$}$\begin{array}{}{\mathrm{􏽥}}_{}\\ {\mathrm{ℰ}}_{\mathit{\epsilon }}\end{array}$ strictly contains the domain of the Γ-limit of ℰε.

This article is the starting point of a series of works whose aim is the study of deterministic control problems where the dynamic and the running cost can be completely different in two (or more) complementary domains of the space ℝN. As a consequence, the dynamic and running cost present discontinuities at the boundary of these domains and this is the main difficulty of this type of problems. We address these questions by using a Bellman approach: our aim is to investigate how to define properly the value function(s), to deduce what is (are) the right Bellman Equation(s) associated to this problem (in particular what are the conditions on the set where the dynamic and running cost are discontinuous) and to study the uniqueness properties for this Bellman equation. In this work, we provide rather complete answers to these questions in the case of a simple geometry, namely when we only consider two different domains which are half spaces: we properly define the control problem, identify the different conditions on the hyperplane where the dynamic and the running cost are discontinuous and discuss the uniqueness properties of the Bellman problem by either providing explicitly the minimal and maximal solution or by giving explicit conditions to have uniqueness.

The adjoint method, recently introduced by Evans, is used to study obstacle problems,weakly coupled systems, cell problems for weakly coupled systems of Hamilton − Jacobiequations, and weakly coupled systems of obstacle type. In particular, new results aboutthe speed of convergence of some approximation procedures are derived.

By disintegration of transport plans it is introduced the notion of transport class. This allows to consider the Monge problem as a particular case of the Kantorovich transport problem, once a transport class is fixed. The transport problem constrained to a fixed transport class is equivalent to an abstract Monge problem over a Wasserstein space of probability measures. Concerning solvability of this kind of constrained problems, it turns out that in some sense the Monge problem corresponds to a lucky case.

This paper is concerned with the construction of local observers for nonlinear systems without inputs, satisfying an observability rank condition. The aim of this study is, first, to define an homogeneous approximation that keeps the observability property unchanged at the origin. This approximation is further used in the synthesis of a local observer which is proven to be locally convergent for Lyapunov-stable systems. We compare the performance of the homogeneous approximation observer with the classical linear approximation observer on an example.

The notion of adequate (resp. strongly adequate) function has been recently introduced tocharacterize the essentially strictly convex (resp. essentially firmly subdifferentiable)functions among the weakly lower semicontinuous (resp. lower semicontinuous) ones. In thispaper we provide various necessary and sufficient conditions in order that the lowersemicontinuous hull of an extended real-valued function on a reflexive Banach space isessentially strictly convex. Some new results on nearest (farthest) points are derivedfrom this approach.

The paper studies discrete/finite-difference approximations of optimal control problemsgoverned by continuous-time dynamical systems with endpoint constraints. Finite-differencesystems, considered as parametric control problems with the decreasing step ofdiscretization, occupy an intermediate position between continuous-time and discrete-time(with fixed steps) control processes and play a significant role in both qualitative andnumerical aspects of optimal control. In this paper we derive an enhanced version of theApproximate Maximum Principle for finite-difference control systems, which is new even forproblems with smooth endpoint constraints on trajectories and occurs to be the firstresult in the literature that holds for nonsmooth objectives and endpoint constraints. Theresults obtained establish necessary optimality conditions for constrained nonconvexfinite-difference control systems and justify stability of the Pontryagin MaximumPrinciple for continuous-time systems under discrete approximations.

We characterize generalized Young measures, the so-called DiPerna–Majda measures whichare generated by sequences of gradients. In particular, we precisely describe thesemeasures at the boundary of the domain in the case of the compactification of ℝm × n by the sphere. We show that this characterization is closely related to the notion of quasiconvexity at the boundary introduced by Ball and Marsden [J.M. Ball and J. Marsden, Arch. Ration. Mech.Anal. 86 (1984) 251–277]. As a consequence we get new results onweak W1,2(Ω; ℝ3) sequentialcontinuity ofu → a· [Cof∇u] ϱ,where Ω ⊂ ℝ3 has a smooth boundary and a,ϱare certain smooth mappings.

Nonlocal calculus is often overlooked in the mathematics curriculum. In this paper we present an interesting new class of nonlocal problems that arise from modelling the growth and division of cells, especially cancer cells, as they progress through the cell cycle. The cellular biomass is assumed to be unstructured in size or position, and its evolution governed by a time-dependent system of ordinary differential equations with multiple time delays. The system is linear and taken to be autonomous. As a result, it is possible to reduce its solution to that of a nonlinear matrix eigenvalue problem. This method is illustrated by considering case studies, including a model of the cell cycle developed recently by Simms, Bean and Koeber. The paper concludes by explaining how asymptotic expressions for the distribution of cells across the compartments can be determined and used to assess the impact of different chemotherapeutic agents.

In 1991, McNabb introduced the concept of mean action time (MAT) as a finite measure of the time required for a diffusive process to effectively reach steady state. Although this concept was initially adopted by others within the Australian and New Zealand applied mathematics community, it appears to have had little use outside this region until very recently, when in 2010 Berezhkovskii and co-workers [A. M. Berezhkovskii, C. Sample and S. Y. Shvartsman, “How long does it take to establish a morphogen gradient?” Biophys. J. 99 (2010) L59–L61] rediscovered the concept of MAT in their study of morphogen gradient formation. All previous work in this area has been limited to studying single-species differential equations, such as the linear advection–diffusion–reaction equation. Here we generalize the concept of MAT by showing how the theory can be applied to coupled linear processes. We begin by studying coupled ordinary differential equations and extend our approach to coupled partial differential equations. Our new results have broad applications, for example the analysis of models describing coupled chemical decay and cell differentiation processes.