We give results for the approximation of a laminate withvarying volume fractions for multi-well energy minimizationproblems modeling martensitic crystals thatcan undergo either an orthorhombicto monoclinic or a cubic to tetragonal transformation. We construct energy minimizing sequences of deformations which satisfythe corresponding boundary condition, and we establish a series of error bounds in terms of the elastic energyfor the approximation of the limiting macroscopicdeformation and the simply laminated microstructure.Finally, we give results for the corresponding finite element approximation of the laminate with varying volume fractions.

In the present work, the symmetrized sequential-parallel decomposition method with the fourth order accuracy for the solution of Cauchy abstract problem with an operator under a split form is presented. The fourth order accuracy is reached by introducing a complex coefficient with the positive real part. For the considered scheme, the explicit a priori estimate is obtained.

We develop a Discrete Element Method (DEM) for elastodynamics using polyhedral elements. We show that for a given choice of forces and torques, we recover the equations of linear elastodynamics in small deformations. Furthermore, the torques and forces derive from a potential energy, and thus the global equation is an Hamiltonian dynamics. The use of an explicit symplectic time integration scheme allows us to recover conservation of energy, and thus stability over long time simulations. These theoretical results are illustrated by numerical simulations of test cases involving large displacements.

We consider models based on conservation laws. For the optimizationof such systems, a sensitivity analysis is essential to determinehow changes in the decision variables influence the objectivefunction. Here we study the sensitivity with respect to the initialdata of objective functions that depend upon the solution of Riemannproblems with piecewise linear flux functions. We presentrepresentations for the one–sided directional derivatives of theobjective functions. The results can be used in the numerical methodcalled Front-Tracking.

A hyperelastic constitutive law, for use in anatomically accurate finite element models ofliving structures, is suggested for the passive and the active mechanical properties of incompressiblebiological tissues. This law considers the passive and active states as a same hyperelastic continuummedium, and uses an activation function in order to describe the whole contraction phase.The variational and the FE formulations are also presented, and the FE code has been validatedand applied to describe the biomechanical behavior of a thick-walled anisotropic cylinder underdifferent active loading conditions.

We propose and analyse a abstract framework for augmented mixed formulations.We give a priori error estimate in the general case: conforming and nonconforming approximations with or without numerical integration.Finally, a posteriori error estimator is given. An example of stabilized formulation for Stokes problem is analysed.

We study order-adaptive implementations of Hermite methods for hyperbolic and singularly perturbed parabolic initial value problems. Exploiting the facts that Hermite methods allow the degree of the local polynomial representation to vary arbitrarily from cell to cell and that, for hyperbolic problems, each cell can be evolved independently over a time-step determined only by the cell size, a relatively straightforward method is proposed. Its utility is demonstrated on a number of model problems posed in 1+1 and 2+1 dimensions.

A numerically inexpensive globalization strategy of sequential quadratic programming methods (SQP-methods) for control of the instationary Navier Stokes equations is investigated. Based on the proper functional analytic setting a convergence analysis for the globalized methodis given. It is argued that the a priori formidable SQP-step can be decomposed into linear primal and linear adjoint systems, which is amenable for existing CFL-software. A report on a numerical test demonstrates the feasibility of the approach.

In this paper, a Dirichlet-Neumann substructuring domaindecomposition method is presented for a finite elementapproximation to the nonlinear Navier-Stokes equations. It isshown that the Dirichlet-Neumann domain decomposition sequenceconverges geometrically to the true solution provided the Reynoldsnumber is sufficiently small. In this method, subdomain problemsare linear. Other version where the subdomain problems are linearStokes problems is also presented.

In this paper, we investigate the coupling between operator splitting techniques and a timeparallelization scheme, the parareal algorithm,as a numericalstrategy for the simulation of reaction-diffusion equations modelling multi-scale reaction waves. This type of problems induces peculiar difficulties and potentially large stiffness which stem from the broad spectrum of temporal scales in the nonlinear chemical source term as well as from the presence of large spatial gradients in the reactive fronts,spatially very localized.In a series of previous studies, the numerical analysis of theoperator splitting as well as the parareal algorithmhas been conducted and such approaches have shown a great potential in the framework of reaction-diffusion and convection-diffusion-reaction systems. However, complementary studies are needed for a more complete characterizationof such techniques for these stiff configurations.Therefore,we conduct in this work a precise numerical analysis that considers thecombination of time operator splitting and the parareal algorithmin the context of stiff reaction fronts. The impact of the stiffness featured by these frontson the convergence of the method is thus quantified, and allows to conclude on an optimal strategy for the resolution ofsuch problems.We finally perform some numerical simulations in the field of nonlinear chemical dynamics that validate the theoretical estimatesand examine the performance of such strategiesin the context of academical one-dimensional test casesas well as multi-dimensional configurationssimulated on parallel architecture.

The chemical master equation is a fundamental equation in chemical kinetics. It underliesthe classical reaction-rate equations and takes stochastic effects into account. In thispaper we give a simple argument showing that the solutions of a large class of chemicalmaster equations are bounded in weighted ℓ1-spaces and possess high-ordermoments. This class includes all equations in which no reactions between two or morealready present molecules and further external reactants occur that add mass to thesystem. As an illustration for the implications of this kind of regularity, we analyze theeffect of truncating the state space. This leads to an error analysis for the finite stateprojections of the chemical master equation, an approximation that forms the basis of manynumerical methods.

Using the min-plus version of the spectral radius formula, one proves: 1) that the unique eigenvalue of a min-plus eigenvalue problem depends continuously on parameters involved in the kernel defining the problem; 2) that the numerical method introduced by Chou and Griffiths to compute this eigenvalue converges.A toolbox recently developed at I.n.r.i.a. helps to illustrate these results. Frenkel-Kontorova models serve as example. The analogy with homogenization of Hamilton-Jacobi equations is emphasized.

We present in this paper a pressure correction scheme for the drift-flux model combining finite element and finite volume discretizations, which is shown to enjoy essential stability features of the continuous problem: the scheme is conservative, the unknowns are kept within their physical bounds and, in the homogeneous case (i.e. when the drift velocity vanishes), the discrete entropy of the system decreases; in addition, when using for the drift velocity a closure law which takes the form of a Darcy-like relation, the drift term becomes dissipative.Finally, the present algorithm preserves a constant pressure and a constant velocity through moving interfaces between phases.To ensure the stability as well as to obtain this latter property, a key ingredient is to couple the mass balance and the transport equation for the dispersed phase in an original pressure correction step.The existence of a solution to each step of the algorithm is proven; in particular, the existence of a solution to the pressure correction step is derived as a consequence of a more general existence result for discrete problems associated to the drift-flux model.Numerical tests show a near-first-order convergence rate for the scheme, both in time and space, and confirm its stability.

In the present work the symmetrized sequential-parallel decomposition methodof the third degree precision for the solution of Cauchy abstract problemwith an operator under a split form, is presented. The third degreeprecision is reached by introducing a complex coefficient with the positivereal part. For the considered schema the explicit a priori estimation isobtained.

We propose, analyze, and compare several numerical methods for the computation of the deformation of a pressurized martensitic thin film. Numerical results have been obtained for the hysteresis of the deformation as the film transforms reversibly from austenite to martensite.

From Kantorovich’s theory we present a semilocal convergence result for Newton’s methodwhich is based mainly on a modification of the condition required to the second derivativeof the operator involved. In particular, instead of requiring that the second derivativeis bounded, we demand that it is centered. As a consequence, we obtain a modification ofthe starting points for Newton’s method. We illustrate this study with applications tononlinear integral equations of mixed Hammerstein type.

We present a numerical algorithm to solve the micromagnetic equations based on tangential-plane minimization for the magnetization update and a homothethic-layer decomposition of outer space for the computation of the demagnetization field. As a first application, detailed results on the flower-vortex transition in the cube of Micromagnetic Standard Problem number 3 are obtained, which confirm, with a different method, those already present in the literature, and validate our method and code. We then turn to switching of small cubic or almost-cubic particles, in the single-domain limit. Our data show systematic deviations from the Stoner-Wohlfarth model due to the non-ellipsoidal shape of the particle, and in particular a non-monotone dependence on the particle size.

We propose a numerical scheme to compute the motion of a two-dimensional rigid body in a viscous fluid. Our method combines the method of characteristics with a finite element approximation to solve an ALE formulation of the problem. We derive error estimates implying the convergence of the scheme.

We prove error estimates for the ultra weak variational formulation (UWVF) in 3D linearelasticity. We show that the UWVF of Navier’s equation can be derived as an upwinddiscontinuous Galerkin method. Using this observation, error estimates are investigatedapplying techniques from the theory of discontinuous Galerkin methods. In particular, wederive a basic error estimate for the UWVF in a discontinuous Galerkin type norm and thenan error estimate in the L2(Ω) norm in terms of the bestapproximation error. Our final result is an L2(Ω) norm errorestimate using approximation properties of plane waves to give an estimate for the orderof convergence. Numerical examples are presented.

The propagation of the action potential in the heart chambers is accurately described bythe Bidomain model, which is commonly accepted and used in the specialistic literature.However, its mathematical structure of a degenerate parabolic system entails highcomputational costs in the numerical solution of the associated linear system. Domaindecomposition methods are a natural way to reduce computational costs, and OptimizedSchwarz Methods have proven in the recent years their effectiveness in accelerating theconvergence of such algorithms. The latter are based on interface matching conditions moreefficient than the classical Dirichlet or Neumann ones. In this paper we analyze anOptimized Schwarz approach for the numerical solution of the Bidomain problem. We assessthe convergence of the iterative method by means of Fourier analysis, and we investigatethe parameter optimization in the interface conditions. Numerical results in 2D and 3D aregiven to show the effectiveness of the method.