Simon Lemaire: Thursday 16 April at 3 pm, A415 Inria Paris. We are interested in physical settings presenting an interface between a classical (positive) material and a (negative) metamaterial, in such a way that the coefficients of the model change sign in the domain. We study, in the « elliptic » case, the numerical approximation of such sign-shifting problems. We introduce a new numerical method, based on domain decomposition and optimization, that we prove to be convergent, as soon as, for a given right-hand side, the problem admits a solution that is unique. The proof of convergence does not rely on any symmetry assumption on the mesh family with respect to the sign-changing interface. In that respect, it gives a more convenient alternative to T-coercivity based approximation in the situations when the latter is applicable, whereas it constitutes a new paradigm in the situations when the latter is not. We illustrate our findings on a comprehensive set of test-cases.
Thirupathi Gudi: Tuesday 20 February at 3 pm, A415 Inria Paris. In this talk, we review some approaches for formulating the Dirichlet boundary control problem and then we present a new energy space based approach. We show that this new approach allows high regularity for both optimal control and the optimal state. Using, the optimality conditions at continuous level, we propose a finite element method for numerical solution and derive subsequent error estimates. We show some numerical experiments to illustrate the method.
Franz Chouly: Thursday 15 February at 2 pm, A415 Inria Paris. In the first part of this talk, we will present a residual based a posteriori error estimate for contact problems in small strain elasticity, discretized with finite elements and Nitsche’s method. Upper and lower bounds are established under a saturation assumption. This theoretical results will be illustrated by some numerical experiments (joint work with Mathieu Fabre, Patrick Hild, Jérôme Pousin and Yves Renard). In the second part of this talk, we will present preliminary results on goal oriented error estimates for soft-tissue biomechanics, still under small strain assumptions. The performance of the Dual Weighted Residual method will be assessed for two simplified scenarios involving tongue muscular activation, and contraction of the arterial wall. Open mathematical questions and the potential interest of such a methodology for computational biomechanics will be discussed (joint work with Stéphane Bordas, Marek Bucki, Michel Duprez, Vanessa Lleras, Claudio Lobos, Alexei Lozinski, Pierre-Yves Rohan and Satyendra Tomar).
Sébastien Furic: Thursday 30 November at 14 pm, A415 Inria Paris. This presentation aims at introducing the concept of System-Level Physical Model, with some emphasis on model construction and simulation. In particular, soundness of models and meaning and solution of associated differential equations will be discussed.
Hend Benameur: Thursday 2 November at 11 am, A415 Inria Paris. We are interested in some inverse problems in porous media: parameter estimation, fracture identification and wells location. All these problems are formulated as optimization problems. The main and common tool in the developed techniques is “ the gradient” of a convenient function. An adaptive parameterization algorithm is developed, implemented and applied for the estimation of scalar and vector parameters in porous media. Values of parameters and shapes of hydrogeological zones are unknown. The main tool in the adaptive parameterization approach is a refinement indicator: Once the identification problem is set as a minimization of an objective function, the question is what is the effect on this function of allowing discontinuity of the parameter in some candidate location? Refinement indicators give the answer to this question . Since fractures are characterized by discontinuities, the idea is to extend previous indicators to locate fractures. We define fracture indicators and we proceed in an iterative way in order to identify fractures in porous media. The topological sensitivity analysis method has been recognized as a promising tool to solve topology optimization problems. It consists to provide an asymptotic expansion of a shape functional with respect to the size of a small hole created inside the domain. To solve the inverse problem where both parametrization and well’s location are unknown, we incorporate the topological gradient approach in the adaptive parametrization algorithm; results are promising.
Peter Minev: Tuesday 10 October at 11 am, A415 Inria Paris. The presentation will be focussed on two classes of recently developed splitting schemes for the Navier-Stokes equations. The first class is based on the classical artificial compressibility (AC) method. The original method proposed by J. Shen in 1995 reduces the solution of the incompressible Navier-Stokes equations to a set of two or three parabolic problems in 2D and 3D correspondingly. Unfortunately, its accuracy is limited to first order in time and can be extended further only if the resulting scheme involves an elliptic problem for the velocity vector. Recently, together with J.L. Guermond (Texas A&M University) we proposed a scheme that extends the AC method to any order in time using a bootstrapping approach to the incompressibility constraint that essentially requires to solve only a set of parabolic equations for the velocity. The conditioning of the corresponding linear systems is therefore O(Δth^-2). This is generally better than solving a parabolic equation for the velocity and an elliptic equation for the pressure required by the various projection schemes that are perhaps the most popular approach at present. Besides, the bootstrapping algorithm allows to achieve any order in time, subject to some initialization conditions, in contrast to the projection methods whose accuracy seems to be essentially limited to second order on the velocity in the L2 norm. The second class of methods is based on a novel approach to the Navier-Stokes equations that reformulates them in terms of stress variables. It was developed in a recent paper together with P. Vabishchevich (Russian Academy of Sciences). The main advantage of such an approach becomes clear when it is applied to fluid-structure interaction problems since in such case the problems for the fluid and the structure, both written in terms of stress variables,…
Théophile Chaumont: Monday 18 September at 1 pm, A415 Inria Paris. Time-harmonic wave propagation problems are costly to solve numerically since the corresponding PDE operators are not strongly elliptic, and as a result, discretization methods might become unstable. Specifically, the finite element solution is quasi-optimal (almost as good as the best approximation the finite element space can provide) only under restrictive assumptions on the mesh size. If the mesh size is too large, stability is lost, and the finite element solution can become completely inaccurate, even when the best approximation is. This phenomenon is called the « pollution effect » and becomes more important for larger frequencies. For the case of wave propagation problems in homogeneous media, it is known that high order finite element methods are less sensitive to the pollution effect. For this reason, they are employed in a wide range of applications, as the corresponding linear systems are smaller and easier to solve. In this talk, we investigate the use of high order finite element methods to solve wave propagation problems in highly heterogeneous media. Since the heterogeneities of the medium can exhibit small scale features, we consider « non-fitting » meshes, that are not aligned with the physical interfaces of the medium. Instead, the parameters defining the medium of propagation can be discontinuous inside each element. We propose a convergence analysis and draw two main conclusions: – the asymptotic convergence rate of the proposed finite element method is suboptimal due to the lack of regularity of the solution inside each cell – the pollution effect is greatly reduced by increasing the order of discretization. We illustrate our main conclusions with geophysical application benchmarks. These examples confirm that higher order methods are more efficient than linear finite elements.
Valentine Rey: Thursday 22 June at 3 pm, A415 Inria Paris. Domain decomposition methods are robust and efficient methods to solve mechanical problems with several million degrees of freedom. Taking advantage of increasing performances of computers, they exploit the clusters-parallel architecture and are numerically scalable. Verification has been widely developed since 1980’s and proposes tools to estimate the distance between the unknown exact solution of continuous problem and the computed solution. Among those techniques, estimators based on error in constitutive relation provide constant-free upper bounds and are available for varied range of problems. In this talk, we present techniques for steering parallel computation by objective of accuracy on quantities of interest. It relies on a parallel error estimator that provides strict guaranteed upper bound and separates the algebraic error (due to the use of iterative solver) from the discretization error (due to the finite element method). This estimator enables to adapt the solver’s stopping criterion to the discretization, which avoids over resolution and useless iterations. In [*], the estimator is used for goal-oriented error estimation and classical bounds for quantities of interest are rewritten in order to separate sources of error. Finally, we benefit the information provided by the error estimator and the Krylov subspaces built during the resolution to set an auto-adaptive strategy (adaptive remeshing and recycling search directions). *V. Rey, P. Gosselet, C. Rey, Strict bounding of quantities of interest in computations based on domain decomposition, Computer Methods in Applied Mechanics and Engineering. 2015 Apr 15;287:212-28
Gouranga Mallik: Thursday 29 June at 3pm, A415 Inria Paris. In this work we consider a priori and a posteriori error control for the nonsingular solution of von Karman plate bending problem. Conforming and nonconforming finite element methods are employed. Existence, uniqueness and error estimates for the discrete solution are presented. We discuss an abstract framework for a posteriori error control which includes conforming and nonconforming finite element methods. This allows us to compute reliable and efficient local estimators. The key ingredients in establishing well-posedness of the discrete problem rely on the linearization of the continuous problem and suitable enrichment operator. Numerical experiments are performed to justify the theoretical results.
Patrik Daniel: Thursday 15 June at 3:30pm, A415 Inria Paris. We propose a new practical adaptive refinement strategy for hp-finite element approximations of elliptic problems. Following some recent theoretical developments in polynomial-degree-robust a posteriori error analysis, we solve two complementary classes of discrete local problems on the vertex-based patches. The first class involves the solution on each patch of a mixed finite element problem with homogeneous Neumann boundary conditions, which leads to an H(div,Ω)-conforming equilibrated flux. This in turns yields a guaranteed upper bound on the error and serves to mark elements for refinement via a Dörfler bulk criterion. The second class of local problems involves the solution, on each marked patch only, of two separate primal finite element problems with homogeneous Dirichlet boundary conditions, which serve to decide between h-, p-, or hp-refinement. Altogether, we show that these ingredients lead to a computable error reduction factor; we guarantee that while performing the hp-adaptive refinement as suggested, the error will be reduced at least by this factor on the next hp-mesh. In a series of numerical experiments in two space dimensions, we first study the accuracy of our predicted reduction factor: in particular, we measure the ratio of the predicted reduction factor relative to the true error reduction, and we find that it is very close to the optimal value of one for both smooth and singular exact solutions. Finally, we study the overall performance of the proposed hp-refinement strategy on some test cases, for which we observe effectivity indices very close to one and exponential convergence rates.