Publications

2025

• An operator approach to the analysis of electromagnetic wave propagation in dispersive media. Part 2: transmission problems.
  
  • Cassier Maxence
  • Joly Patrick
  , 2025. In this second chapter, we analyse transmission problems between a dielectric and a dispersive negative material. In the first part, we consider a transmission problem between two half-spaces, filled respectively by the vacuum and a Drude material, and separated by a planar interface. In this setting, we answer the following question: does this medium satisfy a limiting amplitude principle? This principle defines the stationary regime as the large time asymptotic behavior of a system subject to a periodic excitation. In the second part, we consider the transmission problem of an infinite strip made of a Drude material embedded in the vacuum and analyse the existence and dispersive properties of guided waves. In both problems, our spectral analysis elucidates new and unusual physical phenomena for the considered transmission problems due to the presence of the dispersive negative material. In particular, we prove the existence of an interface resonance in the first part and the existence of slow light phenomena for guiding waves in the second part.
• Multistage stochastic optimization of a mono-site hydrogen infrastructure by decomposition techniques
  
  • Lefgoum Raian
  • Afsar Sezin
  • Carpentier Pierre
  • Chancelier Jean-Philippe
  • de Lara Michel
  Journal of Optimization Theory and Applications, Springer Verlag, 2025, 207 (3), pp.49. The development of hydrogen infrastructures requires to reduce their costs. In this paper, we develop a multistage stochastic optimization model for the management of a hydrogen infrastructure which consists of an electrolyser, a compressor and a storage to serve a transportation demand. This infrastructure is powered by three different sources: on-site photovoltaic panels (PV), renewable energy through a power purchase agreement (PPA) and the power grid. We consider uncertainties affecting on-site photovoltaic production and hydrogen demand. Renewable energy sources are emphasized in the hydrogen production process to ensure eligibility for a subsidy, which is awarded if the proportion of nonrenewable electricity usage stays under a predetermined threshold. We solve the multistage stochastic optimization problem using a decomposition method based on Lagrange duality. The numerical results indicate that the solution to this problem, formulated as a policy, achieves a small duality gap, thus proving the effectiveness of this approach. (10.1007/s10957-025-02795-1)
  DOI : 10.1007/s10957-025-02795-1
• On the constitutive behavior of linear viscoelastic solids under the plane stress condition
  
  • Guzina Bojan B
  • Bonnet Marc
  Journal of Elasticity, Springer Verlag, 2025, 157, pp.45. Motivated by the recent experimental and analytical developments enabling high-fidelity material characterization of (heterogeneous) sheet-like solid specimens, we seek to elucidate the constitutive behavior of linear viscoelastic solids under the plane stress condition. More specifically, our goal is to expose the relationship between the plane-stress viscoelastic constitutive parameters and their (native) ``bulk'' counterparts. To facilitate the sought reduction of the three-dimensional (3D) constitutive behavior, we deploy the concept of projection operators and focus on the frequency-domain behavior by resorting to the Fourier transform and the mathematical framework of tempered distributions, which extends the Fourier analysis to functions (common in linear viscoelasticity) for which Fourier integrals are not convergent. In the analysis, our primary focus is the on class of linear viscoelastic solids whose 3D rheological behavior is described by a set of constant-coefficient ordinary differential equations, each corresponding to a generic arrangement of ``springs'' and ``dashpots''. On reducing the general formulation to the isotropic case, we proceed with an in-depth investigation of viscoelastic solids whose bulk and shear modulus each derive from a suite of classical ``spring and dashpot'' configurations. To enable faithful reconstruction of the 3D constitutive parameters of natural and engineered solids via (i) thin-sheet testing and (ii) applications of the error-in-constitutive-relation approach to the inversion of (kinematic) sensory data, we also examine the reduction of thermodynamic potentials describing linear viscoelasticity under the plane stress condition. The analysis is complemented by a set of analytical and numerical examples, illustrating the effect on the plane stress condition on the behavior of isotropic and anisotropic viscoelastic solids. (10.1007/s10659-025-10136-6)
  DOI : 10.1007/s10659-025-10136-6
• Efficient boundary integral method to evaluate the acoustic scattering from coupled fluid-fluid problems excited by multiple sources
  
  • Pacaut Louise
  • Chaillat Stéphanie
  • Mercier Jean-François
  • Serre Gilles
  Journal of Computational Physics, Elsevier, 2025, 524, pp.113736. In the naval industry many applications require to study the behavior of a penetrable obstacle embedded in water, notably in presence of a turbulent flow. Such configuration is encountered in particular when noise is scattered by two-phase fluids, e.g., turbulent flows with air bubbles. Fast and efficient numerical methods are required to compute this scattering in the presence of realistic 3D geometries, such as bubble curtains. In [1], we have developed a very efficient approach in the case of a rigid obstacle of arbitrary shape, excited by a turbulent flow. It is based on the numerical evaluation of tailored Green's functions. Here we extend this fast method to the case of a penetrable obstacle. It is not a straightforward extension and we propose two main contributions. First, tailored Green's functions for a fluid-fluid coupled problem are derived theoretically and determined numerically. Second, we show the need of a regularized Boundary Integral formulation to obtain these Green's functions accurately in all configurations. Finally, we illustrate the efficiency of the method on various applications related to the scattering by multiple bubbles. (10.1016/j.jcp.2025.113736)
  DOI : 10.1016/j.jcp.2025.113736
• An operator approach to the analysis of electromagnetic wave propagation in dispersive media. Part 1: general results.
  
  • Cassier Maxence
  • Joly Patrick
  , 2025. In this chapter, we investigate the mathematical models for electromagnetic wave propagation in dispersive isotropic passive linear media for which the dielectric permittivity $\varepsilon$ and magnetic permeability $\mu$ depend on the frequency. We emphasize the link between physical requirements and mathematical properties of the models. A particular attention is devoted to the notions of causality and passivity and its connection to the existence of Herglotz functions that determine the dispersion of the material. We consider successively the cases of the general passive media and the so-called local media for which $\varepsilon$ and $\mu$ are rational functions of the frequency. This leads us to analyse the important classes of non-dissipative and dissipative generalized Lorentz models. In particular, we discuss the connection between mathematical and physical properties of models through the notions of stability, energy conservation, dispersion and modal analyses, group and phase velocities and energy decay in dissipative systems.
• A hybridizable discontinuous Galerkin method with transmission variables for time-harmonic acoustic problems in heterogeneous media
  
  • Pescuma Simone
  • Gabard Gwenael
  • Chaumont-Frelet Théophile
  • Modave Axel
  Journal of Computational Physics, Elsevier, 2025, 534, pp.114009. We consider the finite element solution of time-harmonic wave propagation problems in heterogeneous media with hybridizable discontinuous Galerkin (HDG) methods. In the case of homogeneous media, it has been observed that the iterative solution of the linear system can be accelerated by hybridizing with transmission variables instead of numerical traces, as performed in standard approaches. In this work, we extend the HDG method with transmission variables, which is called the CHDG method, to the heterogeneous case with piecewise constant physical coefficients. In particular, we consider formulations with standard upwind and general symmetric fluxes. The CHDG hybridized system can be written as a fixed-point problem, which can be solved with stationary iterative schemes for a class of symmetric fluxes. The standard HDG and CHDG methods are systematically studied with the different numerical fluxes by considering a series of 2D numerical benchmarks. The convergence of standard iterative schemes is always faster with the extended CHDG method than with the standard HDG methods, with upwind and scalar symmetric fluxes. (10.1016/j.jcp.2025.114009)
  DOI : 10.1016/j.jcp.2025.114009
• Global solution of quadratic problems using interval methods and convex relaxations
  
  • Elloumi Sourour
  • Lambert Amélie
  • Neveu Bertrand
  • Trombettoni Gilles
  Journal of Global Optimization, Springer Verlag, 2025, 91 (2), pp.331–353. Interval branch-and-bound solvers provide reliable algorithms for handling non-convex optimization problems by ensuring the feasibility and the optimality of the computed solutions, i.e. independently from the floating-point rounding errors. Moreover, these solvers deal with a wide variety of mathematical operators. However, these solvers are not dedicated to quadratic optimization and do not exploit nonlinear convex relaxations in their framework. We present an interval branch-andbound method that can efficiently solve quadratic optimization problems. At each node explored by the algorithm, our solver uses a quadratic convex relaxation which is as strong as a semi-definite programming relaxation, and a variable selection strategy dedicated to quadratic problems. The interval features can then propagate efficiently this information for contracting all variable domains. We also propose to make our algorithm rigorous by certifying firstly the convexity of the objective function of our relaxation, and secondly the validity of the lower bound calculated at each node. In the non-rigorous case, our experiments show significant speedups on general integer quadratic instances, and when reliability is required, our first results show that we are able to handle medium-sized instances in a reasonable running time. (10.1007/s10898-024-01370-8)
  DOI : 10.1007/s10898-024-01370-8
• A note on pliability and the openness of the multiexponential map in Carnot groups
  
  • Jean Frédéric
  • Sigalotti Mario
  • Socionovo Alessandro
  , 2025. In recent years, several notions of non-rigidity of horizontal vectors in Carnot groups have been proposed, motivated, in particular, by the characterization of monotone sets and Whitney extension properties. In this note we compare some of these notions.