The goal of this paper is to obtain a well-balanced, stable, fast, and robust HLLC-typeapproximate Riemann solver for a hyperbolic nonconservative PDE system arising in aturbidity current model. The main difficulties come from the nonconservative nature of thesystem. A general strategy to derive simple approximate Riemann solvers fornonconservative systems is introduced, which is applied to the turbidity current model toobtain two different HLLC solvers. Some results concerning the non-negativity preservingproperty of the corresponding numerical methods are presented. The numerical resultsprovided by the two HLLC solvers are compared between them and also with those obtainedwith a Roe-type method in a number of 1d and 2d test problems. This comparison shows that,while the quality of the numerical solutions is comparable, the computational cost of theHLLC solvers is lower, as only some partial information of the eigenstructure of thematrix system is needed.

We consider the Laplace equation posed in a three-dimensional axisymmetric domain. Wereduce the original problem by a Fourier expansion in the angular variable to a countablefamily of two-dimensional problems. We decompose the meridian domain, assumed polygonal,in a finite number of rectangles and we discretize by a spectral method. Then we describethe main features of the mortar method and use the algorithm Strang Fix to improve theaccuracy of our discretization.

In this paper we analyze the consistency, the accuracy and some entropy properties ofparticle methods with remeshing in the case of a scalar one-dimensional conservation law.As in [G.-H. Cottet and L. Weynans, C. R. Acad. Sci. Paris, Ser. I343 (2006) 51–56] we re-write particle methods with remeshing inthe finite-difference formalism. This allows us to prove the consistency of these methods,and accuracy properties related to the accuracy of interpolation kernels. Cottet and Magnidevised recently in [G.-H. Cottet and A. Magni, C. R. Acad. Sci. Paris, Ser. I347 (2009) 1367–1372] and [A. Magni and G.-H. Cottet, J.Comput. Phys. 231 (2012) 152–172] TVD remeshing schemes forparticle methods. We extend these results to the nonlinear case with arbitrary velocitysign. We present numerical results obtained with these new TVD particle methods for theEuler equations in the case of the Sod shock tube. Then we prove that with these new TVDremeshing schemes the particle methods converge toward the entropy solution of the scalarconservation law.

Searching for the optimal partitioning of a domain leads to the use of the adjoint methodin topological asymptotic expansions to know the influence of a domain perturbation on acost function. Our approach works by restricting to local subproblems containing theperturbation and outperforms the adjoint method by providing approximations of higherorder. It is a universal tool, easily adapted to different kinds of real problems and doesnot need the fundamental solution of the problem; furthermore our approach allows toconsider finite perturbations and not infinitesimal ones. This paper provides theoreticaljustifications in the linear case and presents some applications with topologicalperturbations, continuous perturbations and mesh perturbations. This proposed approach canalso be used to update the solution of singularly perturbed problems.

The surface Cauchy–Born (SCB) method is a computational multi-scale method for the simulation of surface-dominated crystalline materials. We present an error analysis of the SCB method, focused on the role of surface relaxation. In a linearized 1D model we show that the error committed by the SCB method is 𝒪(1) in the mesh size; however, we are able to identify an alternative “approximation parameter” – the stiffness of the interaction potential – with respect to which the relative error in the mean strain is exponentially small. Our analysis naturally suggests an improvement of the SCB model by enforcing atomistic mesh spacing in the normal direction at the free boundary. In this case we even obtain pointwise error estimates for the strain.

For the efficient numerical solution of indefinite linear systems arising from curl conforming edge element approximations of the time-harmonic Maxwell equation, we consider local multigrid methods (LMM) on adaptively refined meshes. The edge element discretization is done by the lowest order edge elements of Nédélec’s first family. The LMM features local hybrid Hiptmair smoothers of Jacobi and Gauss–Seidel type which are performed only on basis functions associated with newly created edges/nodal points or those edges/nodal points where the support of the corresponding basis function has changed during the refinement process. The adaptive mesh refinement is based on Dörfler marking for residual-type a posteriori error estimators and the newest vertex bisection strategy. Using the abstract Schwarz theory of multilevel iterative schemes, quasi-optimal convergence of the LMM is shown, i.e., the convergence rates are independent of mesh sizes and mesh levels provided the coarsest mesh is chosen sufficiently fine. The theoretical findings are illustrated by the results of some numerical examples.

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 prove the discrete compactness property of the edge elements of any order on a classof anisotropically refined meshes on polyhedral domains. The meshes, made up oftetrahedra, have been introduced in [Th. Apel and S. Nicaise, Math. Meth. Appl.Sci. 21 (1998) 519–549]. They are appropriately graded nearsingular corners and edges of the polyhedron.

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.

We propose a new reduced basis element-cum-component mode synthesis approach forparametrized elliptic coercive partial differential equations. In the Offline stage weconstruct a Library of interoperable parametrized reference componentsrelevant to some family of problems; in the Online stage we instantiate andconnect reference components (at ports) to rapidly form and query parametricsystems. The method is based on static condensation at the interdomainlevel, a conforming eigenfunction “port” representation at the interface level, andfinally Reduced Basis (RB) approximation of Finite Element (FE) bubble functions at theintradomain level. We show under suitable hypotheses that the RB Schur complement is closeto the FE Schur complement: we can thus demonstrate the stability of the discreteequations; furthermore, we can develop inexpensive and rigorous (system-level) aposteriori error bounds. We present numerical results for model many-parameterheat transfer and elasticity problems with particular emphasis on the Online stage; wediscuss flexibility, accuracy, computational performance, and also the effectivity of thea posteriori error bounds.

The numerical approximation of parametric partial differential equations is acomputational challenge, in particular when the number of involved parameter is large.This paper considers a model class of second order, linear, parametric, elliptic PDEs on abounded domain D with diffusion coefficients depending on the parametersin an affine manner. For such models, it was shown in [9, 10] that under very weak assumptionson the diffusion coefficients, the entire family of solutions to such equations can besimultaneously approximated in the Hilbert space V = H01(D) by multivariate sparse polynomials in the parametervector y with a controlled number N of terms. Theconvergence rate in terms of N does not depend on the number ofparameters in V, which may be arbitrarily large or countably infinite,thereby breaking the curse of dimensionality. However, these approximation results do notdescribe the concrete construction of these polynomial expansions, and should thereforerather be viewed as benchmark for the convergence analysis of numerical methods. Thepresent paper presents an adaptive numerical algorithm for constructing a sequence ofsparse polynomials that is proved to converge toward the solution with the optimalbenchmark rate. Numerical experiments are presented in large parameter dimension, whichconfirm the effectiveness of the adaptive approach.

Space-time approximations of the FitzHugh–Nagumo system of coupled semi-linear parabolicPDEs are examined. The schemes under consideration are discontinuous in time butconforming in space and of arbitrary order. Stability estimates are presented in thenatural energy norms and at arbitrary times, under minimal regularity assumptions.Space-time error estimates of arbitrary order are derived, provided that the naturalparabolic regularity is present. Various physical parameters appearing in the model aretracked and numerical examples are presented.

We consider a two-dimensional quantum waveguide composed of two semi-strips of width 1and 1 − ε, where ε > 0 is a small real parameter,i.e. the waveguide is gently converging. The width of the junction zonefor the semi-strips is 1 + O(√ε). We will present a sufficient condition for the existence of a weaklycoupled bound state below π2, the lower bound of thecontinuous spectrum. This eigenvalue in the discrete spectrum is unique and itsasymptotics is constructed and justified whenε → 0+.