Publications

2026

• Surface Plasmon Polariton Excitation in Time-modulated Media
  
  • Raziman T. V.
  • Touboul Marie
  • Sapienza Riccardo
  • Craster Richard V.
  • Rodríguez-Fortuño Francisco J.
  , 2025. Surface plasmon polaritons (SPPs) are central to application areas such as sensing, energy harvesting, and nanoscale optics, and are typically excited via spatial structuring -- an approach lacking dynamic control. We demonstrate that step-like time modulation allows the excitation and out-coupling of SPPs through modulating a dispersive Drude-like metal thereby modelling realistic transparent conducting oxides; this establishes a pathway for the active control and extraction of plasmons in experimentally viable time-varying systems. Using finite-difference time-domain simulations we show that time modulation facilitates both the launching and radiation of surface plasmons with frequencies governed by the dispersion of the bounding media. Our results also reveal the generation of time-reflected waves and the emergence of a magnetostatic mode required for matching boundary conditions at the temporal interface.
• Interacting particles in random irregular media
  
  • Bondi Luca
  • Issoglio Elena
  • Russo Francesco
  , 2026. In this paper, we first develop novel tools of independent interest in the field of stochastic calculus related to SDEs with distributional drift, namely the formulation of a quasi-Stroock-Varadhan martingale problem in the sense that the test functions are of class $C^{1,2}$ instead of solutions of Kolmogorov type PDE. This has several applications, among which, the time restriction and the notion of conditional rough martingale problem and related solution.
• Analysis of a two-level domain decomposition preconditioner for the time-harmonic Maxwell equations with heterogeneous and anisotropic coefficients
  
  • Bonazzoli Marcella
  • Ciarlet Patrick
  • Modave Axel
  • Rappaport Ari
  , 2026. We analyze a two-level additive Schwarz domain decomposition preconditioner for the time-harmonic Maxwell equations with heterogeneous and anisotropic material coefficients. The electric permittivity ε, magnetic permeability μ, and conductivity σ are spatially varying tensor fields, assumed to belong to PW¹,∞(Ω) relative to a fixed finite material partition and to be uniformly symmetric positive definite (USPD). Convergence estimates for the preconditioned GMRES solver are obtained through bounds on the norm and the field of values of the preconditioned operator. Our purpose is to extend the analysis of Bonazzoli et al. from homogeneous scalar material laws to spatially heterogeneous tensor-valued coefficients. While the overall argument follows the additive Schwarz framework therein, this extension requires substantial new ingredients. Among these are coefficient-weighted discrete decompositions and liftings, as well as regularity estimates adapted to material interfaces and piecewise-regular tensors. The resulting field-of-values and GMRES estimates depend explicitly on global spectral contrast factors that combine spatial heterogeneity and pointwise anisotropy. They cover the “proportional conductivity” baseline σ = cε for a scalar c > 0 and normalized deviations that are sufficiently small relative to the Schwarz resolution scales. Across the three tested frequencies in a resolution-scaled sweep, the two-level iteration counts are nearly constant at the highest absorption ratio, grow modestly at the intermediate ratio, and show clearer growth at the lowest ratio. The two-level method consistently requires fewer iterations than the one-level method. In the tested fixed-frequency problems, increasing either directional anisotropy or spatial heterogeneity degrades convergence, while reducing the coarse-mesh diameter largely offsets this deterioration. Numerically, non-proportional conductivity produces essentially the same iteration counts as a proportional control in the tested fixed-frequency configurations. These experiments provide evidence of convergence across a broader conductivity class than the sufficient field-of-values condition captures, without establishing parameter-robust behavior outside that condition.
• Circuit Design Problems in Binary Polynomial Quantum Optimization
  
  • Laplace Mermoud Dylan
  • Watel Dimitri
  • Elloumi Sourour
  • Simonetto Andrea
  , 2026. Solving combinatorial optimization problems on near-term quantum computers mainly relies on their quadratic unconstrained binary formulation. Lately, binary polynomial formulations start to take the center stage, having a better scaling with respect to the number of used qubits. However, while quadratic forms are straightforward to implement on quantum circuits, polynomial forms of degree d incur in circuit depth overheads, due to the lack of d-qubit gates, for d≥3, on current hardware. In this paper, we study the problem of designing the shallowest quantum circuit for a given polynomial form. We formulate it as a classical binary linear problem and we prove the NP-completeness of a variant of it. On the algorithmic side, we offer a scalable and performant heuristic method to tackle it, Tempura, as well as faster binary linear problem formulations, which offer computational feasibility while sacrificing optimality. Our performance evaluation shows that the algorithms yield far shorter circuits than the state-of-the-art approaches, and Tempura is mostly comparable with them in terms of runtime
• Optimized Schwarz Methods in Time for Discrete Transport Control
  
  • Bui Duc-Quang
  • Delourme Bérangère
  • Halpern Laurence
  • Kwok Felix
  Numerical Algorithms, Springer Verlag, 2026. We investigate optimized Schwarz domain decomposition methods in time for the control of the 1D transport equation. In the case of an internal control over the whole domain, the optimization problem can be transformed into a system of two coupled PDEs. We then apply the time-domain decomposition (without overlap) strategy on this PDE system as well as on its discretized counterpart. Under Fourier analysis, we analyse three different iterations: the fixed point iteration, the relaxed iteration and the preconditioned GMRES method. For each case, we propose parameters for the transmission conditions that lead to fast convergence of the method. We illustrate our results by numerical examples. (10.1007/s11075-025-02270-y)
  DOI : 10.1007/s11075-025-02270-y
• Optimization-based numerical methods for the two-and three-dimensional time harmonic Maxwell problem with sign-changing coefficients
  
  • Chaaban Farah
  • Ciarlet Patrick
  • Rihani Mahran
  , 2026. We study the two-and three-dimensional time-harmonic Maxwell equations in a medium composed of a classical material and a metamaterial. The problem is modeled as an interface problem with sign-changing coefficients, with the electric field as the main vector-valued unknown. To carry out the theoretical and numerical analyses, we only require that the problem has a unique solution. We then construct a convergent numerical approximation, with the help of Nédélec finite elements and prove its convergence in the strong H(curl)-norm. At the discrete level, we consider two approaches. The first directly adapts the arguments used for the continuous problem, while the second uses the extra regularity of the tangential traces of discrete fields. Finally, numerical experiments illustrate the theoretical results.
• Tangential and normal traces for extension domains with non-Lipschitz boundaries
  
  • Cervera Romain
  • Ciarlet Patrick
  • Rozanova-Pierrat Anna
  • Teplyaev Alexander
  , 2026. We generalize the classical vector-valued tangential and normal trace theory on Lipschitz domains to the setting of non-Lipschitz $H^1$-extension domains, which includes domains with fractal boundaries such as the Koch snowflake. We define and study these operators based on the surjectivity of the trace operator for elements of $H^1(\Omega)$ and the Hilbert structure of the associated trace space. The normal trace operator is defined on $\boldsymbol{H}({\rm div},\Omega)$ in $\mathbb{R}^n$. A generalized Stokes formula allows us to introduce the tangential trace operator on $\boldsymbol{H}({\bf curl},\Omega)$ and $\boldsymbol{H}^1(\Omega)$ in two and three dimensions. Following the approach of Buffa, Costabel and Sheen (2002) for Lipschitz domains, we define two abstract tangential boundary spaces as images of these trace operators, establish their Hilbert structure, and construct an abstract rotation operator linking them, in place of the geometric rotation based on the normal vector. We also extend Costabel's (1991) coercive bilinear form approach for Maxwell's equations to non-Lipschitz domains, yielding a theory that supports the treatment of the Hodge-Dirac operator and Green's formulas and enables the solution of boundary-value problems for the $\operatorname{\textbf{curl}}\operatorname{\textbf{curl}}+\textbf{I}$ operator within the $H^1$-extension domains framework.
• Diffraction of waves by finite thin coatings : asymptotic expansion and effective models
  
  • Baudet Cédric
  , 2026. We consider the Helmholtz equation in a Dirichlet angular sector (which may be locally perturbed) covered on one side by a thin homogeneous layer. First, we propose an asymptotic expansion of the solution at any order and for any angle when the thickness of the layer tends to 0, using a matched asymptotic expansion. This consists in introducing different asymptotic expansions in three subdomains – the vicinity of the corner, the layer and the rest of the domain – that are linked through matching conditions. The presence of the corner makes these matching conditions delicate to derive because the fields have complex singular behaviors. We overcome this difficulty using a novel system of algebraic tools, including formal series and operators on these series. Moreover, we justify the asymptotic expansion with error estimates and provide details on the form of the expansion depending on the angle.We numerically validate this expansion up to order 5 for several angles. This requires computing the behavior at infinity of some corner profiles. To do this, we introduce a Galerkin-type method enriched by complements at infinity. These complements are known functions that appear in the behavior at infinity of the corner profiles with coefficients that become unknowns of the system.Finally, we propose effective models of order 3, resp. 5 where, far from the corner, the layer is replaced with an effective Robin, resp. Ventcel, boundary condition and, near the corner, an artificial boundary is introduced and a Dirichlet-to-Neumann (DtN) type condition is imposed. The DtN map relies on a modal decomposition of the truncated asymptotic expansion at order 3, resp. 5. The modes take into account the singular behavior of the expansion terms, the corner profiles and some algebraic matching coefficients. We show that our effective models are well-posed when the layer is thin enough, and rigorously justify these models with error estimates. These results are illustrated by numerical simulations.
• OPTEC (Optimized Parallelized Task Execution in C++)
  
  • Baudet Cédric
  , 2026. OPTEC (Optimized Parallelized Task Execution in C++) is a C++ library. Its main features are as follows: – Optimized parallel execution of a task list: each process is handed a new task as soon as it finishes the previous one, until the task list is exhausted. – Data communication between tasks through a simple send/request interface. – Estimation of the remaining computation time, based on the tasks' time complexity or their completion percentage (also useful for non-parallelized computations). – Memory (RAM) resource management, with automatic pausing of some processes to avoid saturating the RAM. More advanced features: task-order optimization, computation of undetermined constants appearing in the complexity formula, sending of custom types between processes, compatibility with the FEM-BEM library XLiFE++, etc.
• Space and Time Decompositions for LQG Problems in Decision-Hazard Information Structure
  
  • Carpentier Pierre
  , 2026. <div><p>In this paper, we implement spatial decomposition, temporal block decomposition, and the combination of these two decompositions for a Linear Quadratic Gaussian (LQG) problem that can be decomposed into subproblems linked by linear coupling constraints, with an optimization span comprising two time scales. This allows us to solve only quadratic problems under linear constraints during the decompositions, which we address by establishing the Riccati equation adapted to the different cases considered, and thus to be able to efficiently solve the different subproblems appearing in the decompositions. We can also compute the solution of the global problem using Riccati and thus measure the quality of the solutions obtained by the decompositions.</p></div>
• Pontryagin maximum principle for non-smooth state-constrained control problems over Wasserstein spaces
  
  • Urrea Fernanda
  • Treumún Ernesto
  , 2026. We study a state-constrained Bolza optimal control problem governed by a non-local continuity equation in the Wasserstein space $\sP_2(\R^d)$. We develop a Pontryagin maximum principle for these problems in which both the cost functionals and the constraints are only locally Lipschitz, and the initial measure does not necessarily have compact support. To derive the main result, we introduce a notion of a Clarke-type deterministic subdifferential over $\sP_2(\R^d)$. The proof of the main result combines a careful adaptation of the Lasry-Lions regularization to the Wasserstein space coupled with stability results of the Clarke subdifferential that allow us to recover the first-order information. We also show that, for functionals convex along plans, this subdifferential admits a global supporting characterization and provides a first-order characterization of global minimizers.
• Fast SDP certification of neural networks : towards large multi-class datasets
  
  • Boyer Margot
  • Rambour Clément
  • Alès Zacharie
  • Lambert Amélie
  , 2026. <div><p>We present a new quadratic model for the certification problem in adversarial robustness, which simultaneously accounts for all possible target classes. Building on this model, we propose a novel semidefinite programming (SDP) relaxation for incomplete verification. A key advantage of our approach is that it certifies robustness in a single optimization, avoiding the need for a separate resolution per class. This yields a significant computational speed-up and enables scalability to large datasets with many classes. To further improve efficiency, we also propose an effective pruning strategy of active neurons, thus reducing the problem dimensionality and accelerating convergence.</p></div>
• Accelerating the Method of Reflections with Domain Decomposition techniques for Boundary Integral Equations in Multiple Scattering
  
  • Chaillat Stéphanie
  • Darbas Marion
  • Gander Martin J
  • Halpern Laurence
  BIT Numerical Mathematics, Springer Verlag, 2026. The Method of Reflections was historically introduced to obtain approximate solu-tions as series expansions for the motion of particles in suspension. It can however equally well be used for solving multiple scattering problems numerically. We show for Helmholtz multiple scattering problems that the Method of Reflections, whether applied in its alternating or parallel version, suffers from convergence problems when scatterers are close. We use boundary integral equations to formulate the methods, and then identify them as algebraic Schwarz methods, thereby interpreting them as boundary domain decomposition techniques. This connection allows us to introduce remedies such as overlap (which can be partial, covering only the illuminating region of the obstacles) and coarse spaces from domain decomposition into the Method of Reflections. This leads to substantially accelerated variants, and also naturally makes them suitable preconditioners for GMRES. These new approaches are particularly efficient for closeby obstacles. Moreover, numerical experiments show that the number of iterations remains robust with respect to the wavenumber. (10.1007/s10543-026-01151-7)
  DOI : 10.1007/s10543-026-01151-7
• A note on pliability and the openness of the multiexponential map in Carnot groups
  
  • Jean Frédéric
  • Sigalotti Mario
  • Socionovo Alessandro
  Symmetry, Integrability and Geometry : Methods and Applications, National Academy of Science of Ukraine, 2026. 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. (10.3842/SIGMA.2026.078)
  DOI : 10.3842/SIGMA.2026.078
• Local Multiple Traces Formulation for Transmission Problems in Linear Elasticity
  
  • Chaillat Stéphanie
  • Darbas Marion
  • Escapil-Inchauspé Paul
  • Jerez-Hanckes Carlos
  Computers & Mathematics with Applications, Elsevier, 2026, 220 (15). We investigate the use of the local Multiple Trace Formulation (MTF) in solving time-harmonic elastic wave transmission problems. Originally devised for heterogeneous acoustic media, MTF recasts the boundary value problem as a well-posed system of first-kind boundary integral equations, naturally amenable to parallelization and preconditioning. The formulation takes independent displacement and traction unknowns per subdomain, enforces Calderón identities locally, and imposes transmission conditions weakly across interfaces. We restrict our analysis to single homogeneous scatterers, representing a foundational step toward the application of MTF techniques in heterogeneous elasticity transmission problems. For the sake of clarity, all derivations in the one-dimensional setting are carried out explicitly, with illustrative examples provided in two dimensions. We analyze the effects of frequency and material contrast on the convergence of the GMRES iterative solver. Finally, we present preliminary results for an elastic Calderón preconditioner and discuss its potential to further accelerate iterative solvers. (10.1016/j.camwa.2026.08.005)
  DOI : 10.1016/j.camwa.2026.08.005
• A hybridizable discontinuous Galerkin method with transmission variables for time-harmonic electromagnetic problems
  
  • Rappaport Ari
  • Chaumont-Frelet Théophile
  • Modave Axel
  SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2026, 48 (4), pp.B629-B650. The CHDG method is a hybridizable discontinuous Galerkin (HDG) finite element method suitable for the iterative solution of time-harmonic wave propagation problems. Hybrid unknowns corresponding to transmission variables are introduced at the element interfaces and the physical unknowns inside the elements are eliminated, resulting in a hybridized system with favorable properties for fast iterative solution. In this paper, we extend the CHDG method, initially studied for the Helmholtz equation, to the time-harmonic Maxwell equations. We prove that the local problems stemming from hybridization are well-posed and that the fixed-point iteration naturally associated to the hybridized system is contractive. We propose a 3D implementation with a discrete scheme based on nodal basis functions. The resulting solver and different iterative strategies are studied with several numerical examples using a high-performance parallel C++ code. (10.1137/25M1756818)
  DOI : 10.1137/25M1756818
• Exponential twist of probability measures: drift correction in term of a generalized gradient
  
  • Bourdais Thibaut
  • Oudjane Nadia
  • Russo Francesco
  Stochastic Processes and their Applications, Elsevier, 2026. In this paper we study the exponential twist, i.e. a path-integral exponential change of measure, of a Markovian reference probability measure $\P$. This type of transformation naturally appears in variational representation formulae originating from the theory of large deviations and can be interpreted in some cases, as the solution of a specific stochastic control problem. Under a very general Markovian assumption on $\P$, we fully characterize the exponential twist probability measure as the solution of a martingale problem and prove that it inherits the Markov property of the reference measure. The ''generator'' of the martingale problem shows a drift depending on a {\it generalized gradient} of some suitable {\it value function} $v$. The analysis focuses on the fact that any Markovian probability measure fulfills an {\it intrinsic martingale problem} for which no uniqueness is required. (10.1016/j.spa.2026.105064)
  DOI : 10.1016/j.spa.2026.105064
• Multi-domain FEM-BEM coupling with several impenetrable obstacles
  
  • Boisneault Antonin
  • Bonazzoli Marcella
  • Claeys Xavier
  • Marchand Pierre
  , 2026. In a recent paper we have analyzed a new formulation of the coupling of finite and boundary element methods (FEM-BEM) for Helmholtz problems, involving several heterogeneous bounded subdomains and one homogeneous unbounded subdomain. This formulation, called Generalized Optimized Schwarz Method (GOSM), is substructured, that is, its unknowns are associated with the subdomains interfaces. To derive the GOSM, the first step was to prove that a solution to the Helmholtz problem satisfies a specific multi-domain variational formulation, which involves one operator for each subdomain, each operator being independent of the others. In the present contribution, we design a variational formulation of that type for a more general geometrical and material configuration: several heterogeneous bounded subdomains, impenetrable obstacles and homogeneous subdomains are allowed. Note that, like in our recent paper, we assume that only one subdomain is unbounded, and that its boundary is bounded. The domain partition can have cross-points, that is, points where at least three subdomains are adjacent. We also prove that a solution to the initial Helmholtz problem can be recovered from a solution to the multi-domain variational formulation. This shows that the GOSM is a general and flexible framework to model acoustic wave propagation, as it can handle both multi-domain FEM-BEM coupling and (weakly imposed) boundary conditions on several obstacles.
• Analysis of the geometric spectral properties of electromagnetic waveguides
  
  • Briet Philippe
  • Cassier Maxence
  • Ourmières-Bonafos Thomas
  • Zaccaron Michele
  , 2026. <div><p>We consider a homogeneous isotropic electromagnetic waveguide with simply connected crosssection, embedded in a perfect conductor. We study how bending and twisting affect the real spectrum of the associated self-adjoint Maxwell operator. We first give geometric conditions ensuring preservation of the essential spectrum, which is symmetric, has a gap around the origin, and coincides with the spectrum of a reference waveguide with constant twist and zero curvature. We then provide in the case of purely bended waveguides sufficient conditions for the appearance of discrete eigenvalues in the gap (associated with trapped modes) and analyze them both theoretically and numerically.</p></div>
• High-order numerical integration on self-affine sets
  
  • Joly Patrick
  • Kachanovska Maryna
  • Moitier Zoïs
  SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2026, 48 (4), pp.A1870-A1897. We construct an interpolatory high-order cubature rule to compute integrals of smooth functions over self-affine sets with respect to an invariant measure. The main difficulty is the computation of the cubature weights, which we characterize algebraically, by exploiting a self-similarity property of the integral. We propose an \( h \)-version and a \( p \)-version of the cubature, present an error analysis and conduct numerical experiments. (10.1137/24M1697141)
  DOI : 10.1137/24M1697141
• Douglas-Rachford splitting method for nonconvex stochastic optimal control
  
  • Caffier François
  • Malisani Paul
  • Pfeiffer Laurent
  , 2026. Motivated by the numerical resolution of stochastic optimal control problems, particularly those involving nonlinear dynamics, we investigate the convergence properties of the Douglas-Rachford splitting algorithm in a nonconvex and infinite-dimensional setting. We introduce an extension of the method that applies when one component of the objective function is continuously differentiable with a locally Lipschitz gradient, while the other is weakly lower semicontinuous with a bounded domain.
• High-Resolution Inertial Dynamics with Time-Rescaled Gradients for Nonsmooth Convex Optimization
  
  • Le Manh Hung
  • Simonetto Andrea
  Computational Optimization and Applications, Springer Verlag, 2026. We study nonsmooth convex minimization through a continuous-time dynamical system that can be seen as a high-resolution ODE of Nesterov Accelerated Gradient (NAG) adapted to the nonsmooth case. We apply a time-varying Moreau envelope smoothing to a proper convex lower semicontinuous objective function and introduce a controlled time-rescaling of the gradient, coupled with a Hessian-driven damping term, leading to our proposed inertial dynamic. We provide a well-posedness result for this dynamical system, and construct a Lyapunov energy function capturing the combined effects of inertia, damping, and smoothing. For an appropriate scaling, the energy dissipates and yields fast decay of the objective function and gradient, stabilization of velocities, and weak convergence of trajectories to minimizers under mild assumptions. Conceptually, the system is a nonsmooth high-resolution model of Nesterov's method that clarifies how time-varying smoothing and time rescaling jointly govern acceleration and stability. We further extend the framework to the setting of maximally monotone operators, for which we propose and analyze a corresponding dynamical system and establish analogous convergence results. We also present numerical experiments illustrating the effect of the main parameters and comparing the proposed system with several benchmark dynamics. (10.1007/s10589-026-00805-0)
  DOI : 10.1007/s10589-026-00805-0
• Risk-Averse Control for Continuous-Time Stochastic System Under Signal Temporal Logic Constraints
  
  • Lai En
  • Bonalli Riccardo
  • Girard Antoine
  • Jean Frédéric
  , 2026. Signal Temporal Logic (STL) has become a powerful formalism for specifying complex temporal-spatial behaviors in autonomous systems. Handling STL constraints within stochastic setting has received increasing research interest but still poses challenges. This paper proposes a general framework to efficiently solve continuous-time nonlinear stochastic optimal control problems under chance STL constraints. The STL formulae are implemented through extended dynamics, yielding a more classical chance constraint on the terminal state uniquely that we reliably relax via Conditional Value-at-Risk. The resulting new optimal control problem is then solved using established algorithms from risk--averse control. The efficiency and feasibility of the proposed approach are demonstrated through numerical simulations.
• Transport of dissolved and particulate arsenic in the Orbiel River during minor flood events downstream the historic Salsigne gold mining district, Southern France
  
  • Carrière Lali
  • Resongles Eléonore
  • Heydon Marie
  • Roux Hélène
  • Viers Jérôme
  • Schreck Eva
  • Freydier Rémi
  • Marchand Pierre
  • Domeau Aurélien
  • Horgue Pierre
  • Behra Philippe
  • Casiot Corinne
  Journal of Contaminant Hydrology, Elsevier, 2026, 281, pp.104955. Understanding and quantifying the transport of toxic metals and metalloids during flood events downstream of former mining sites is challenging as these events are transient, highly variable and unpredictable, resulting in a scarcity of observational data across watersheds worldwide. In this study, the concentrations and fluxes of dissolved and particulate arsenic (As) were determined to investigate the dynamic of As transfer during flood events in the Orbiel River catchment, in Southern France. This river drains the former Salsigne mining district, where gold-bearing arsenopyrite mineralization was historically exploited. In this Mediterranean catchment, characterized by long dry periods and episodic rainfall, high-frequency automatic sampling of river water was conducted at four stations during five minor to moderate flood events between 2021 and 2023. Arsenic dynamics during each flood were event-specific and varied according to the antecedent hydrological conditions of the basin and the respective contribution of mine-impacted and non-impacted tributaries, which depend on localized rainfall patterns. Dissolved As concentrations ranged from 10 to 35 μg L-1 at the outlet of the catchment, comparable to baseflow levels, while particulate As varied from 20 to 500 mg kg-1, with enrichment factors up to 29 relative to the local geochemical background. Although these recurring, low-intensity floods did not generate strong As concentration peaks, the total As load (i.e., the sum of dissolved and particulate As) transported during the quick-flow phases of individual flood events ranged from 7.5 to 93.2 kg at the catchment outlet, reaching up to 3.9 times that measured during a 50-day low-flow summer period (24 kg). These results demonstrate that recurring minor floods, which occur several times each year, substantially contribute to As export from the Orbiel River catchment despite their moderate magnitude. These findings highlight the need to account for frequent low-intensity floods in contaminant transport assessments and management of legacy mining areas under climate change. (10.1016/j.jconhyd.2026.104955)
  DOI : 10.1016/j.jconhyd.2026.104955
• Low-thrust Interplanetary Trajectories with Missed Thrust Events: a Numerical Approach
  
  • Carpentier Pierre
  • Chancelier Jean-Philippe
  • Cohen Guy
  • Dargent Thierry
  • Epenoy Richard
  Journal of Guidance, Control, and Dynamics, American Institute of Aeronautics and Astronautics, 2026, pp.1-18. The problem under consideration is to drive a spatial vehicle to a target at a given final time while minimizing fuel consumption. This is a classical optimal control problem in a deterministic setting. However temporary stochastic failures of the engine may prevent reaching the target after the engine usage is recovered. Therefore, a stochastic optimal control problem is formulated under the constraint of ensuring a minimal probability of hitting the target. This problem is modeled, improved and finally solved by dualizing the probability constraint and using an Arrow-Hurwicz stochastic algorithm. Numerical results concerning an interplanetary mission are presented. (10.2514/1.G009777)
  DOI : 10.2514/1.G009777