Mathematics to keep the power grid running: François Pacaud’s research

Research Science and society Decoding
Published on 14 March 2026
Delivering electricity to the right place at the right time and in complete safety is a real balancing act on a national scale. Every day, decisions must be made between production, transmission, and consumption on networks with tens of thousands of connection points, while anticipating unforeseen events and peaks in demand.
To meet this challenge, operators rely on highly sophisticated mathematical models. On International Mathematics Day, Mines Paris – PSL is highlighting the work of François Pacaud, a professor and researcher at the Centre Automatique et Systèmes (CAS), who is developing new algorithms to optimize large electrical networks through mathematical optimization and high-performance computing.
 

Mathematics driving electrical networks

Modern electrical networks are among the most complex infrastructures ever built. They must constantly maintain a balance between what is produced and what is consumed, integrate inherently fluctuating renewable energies, and ensure the stability of the system in the face of climatic or technical uncertainties. To achieve this, operators rely not only on physical equipment, but also on mathematical models capable of describing the grid as a whole and optimizing its operation.

At Mines Paris – PSL, François Pacaud works specifically on these mathematical tools. A professor and researcher at CAS, he develops optimization methods based on applied mathematics and scientific computing. His research aims to make optimization problems that have long been considered too large or too complex to be solved quickly calculable in practice.

Optimizing an electrical network

A huge challenge

Managing an electrical network is not just a matter of moving electricity from point A to point B. It involves constantly deciding how to distribute electricity flows while complying with numerous physical constraints (Kirchhoff’s laws) and technical constraints (line capacities, system security, and economic costs). These issues are formalized in the form of nonlinear optimization problems that calculate the optimal power flow (OPF) to minimize losses within the grid.

On the scale of large interconnected grids, such as the European grid, these models reach dizzying proportions, with several hundred thousand variables and constraints. They must also take into account uncertainties, such as those related to weather or variations in wind and solar production. Solving these problems accurately and in a timely manner is therefore a key challenge, both for improving grid stability and for the gradual integration of new energy sources, in a context where the energy transition is becoming essential.

Map of the European electricity transmission network (ENTSO-E), showing the infrastructure operated by ENTSO-E member operators, both existing and under construction (power plants, converters, substations, and high-voltage lines or cables).

Large-scale nonlinear optimization

Equations to maintain equilibrium

To solve these problems, François Pacaud relies on methods derived from mathematical optimization, in particular methods known as “interior point methods.” The principle behind these methods is to gradually approach the best possible solution by performing a series of successive calculations. At each stage, the algorithm must solve large systems of equations, known as Karush–Kuhn–Tucker (KKT) conditions, which mathematically translate all of the network’s physical and technical constraints.

The main obstacle is not so much the formulation of the problem as its numerical solution. KKT systems are huge, mostly composed of zeros (they are said to be “hollow”) and often difficult to manipulate numerically. Their size and structure make calculations long and memory-intensive. Scaling, i.e., the ability to process entire networks rather than simplified models, is one of the central challenges in this field.

Accelerating calculations with GPUs

A major advance in François Pacaud’s recent work is to exploit the computing power of graphics processing units (GPUs), which are massively parallel architectures, to solve these optimization problems. GPUs work differently from conventional processors: they are particularly good at performing thousands of small calculations in parallel, rather than solving them step by step on a single processor. GPUs thus contributed greatly to the deep learning revolution in the 2010s.

In collaboration with researchers at MIT and Argonne National Laboratory, he contributed to the development of new methods known as “condensed space” methods, which consist of reformulating KKT equation systems in a mathematical form better suited to parallel computing. These reformulations make it possible to use specialized numerical libraries capable of efficiently factoring and solving very large sparse systems on GPUs. Thanks to these approaches, it is now possible to solve problems representing an entire large national power grid in less than twenty seconds, where hours of computation were previously required.

Beyond the performance gains, these results make it possible to optimize very large networks in near real time, which is essential for dynamic and secure management of energy infrastructures.

Nvidia GeForce 6600GT graphics processing unit (GPU)

A coherent set of optimization projects

This research is not limited to electrical networks. It focuses more broadly on large-scale nonlinear and stochastic optimization, i.e., finding the best possible solution when mathematical relationships are nonlinear and certain data are uncertain or variable, for example due to weather or user behavior.

To this end, François Pacaud develops and improves solvers, software programs capable of automatically solving very large mathematical problems. Among them, MadNLP and MadNCL are optimization solvers designed to handle thousands or even millions of variables, numerous physical constraints, and phenomena of uncertainty. These tools have a variety of applications: energy networks, signal processing, and robotics. The challenge is not only to find a solution, but to ensure that it is reliable, stable, and calculable within a reasonable time frame, which until now has been beyond the capabilities of conventional algorithms.

Between theoretical mathematics and industrial applications

Trained in applied mathematics, François Pacaud worked as an optimization engineer in industry before pursuing his research in major international laboratories. Today at Mines Paris – PSL, he develops advanced computational mathematics methods while collaborating with players in the energy sector. His career path reflects the link between theory and practice that is so important to the School, with a dual focus on science and industry. Mathematics thus becomes a decision-making tool, capable of transforming vast sets of data and physical constraints into concrete strategies for infrastructure management.

Mathematics plays a central role in modern energy infrastructure. Optimizing an electrical grid is not just about solving an equation: it is about ensuring the stability of a vital system, supporting the integration of renewable energies, and securing the electricity supply. On this International Mathematics Day, the research conducted at Mines Paris – PSL illustrates the power of mathematical tools when combined with scientific computing and engineering. Far from being abstract, equations can make a very concrete contribution to meeting today’s energy challenges.


To go further :

  • François Pacaud, Sungho Shin, Alexis Montoison, Michel Schanen, Mihai Anitescu. Condensed-space methods for nonlinear programming on GPUs. 2024. ⟨hal-04875904⟩ 
  • Wei Kuang, Alexis Montoison, Vishwas Rao, François Pacaud, Mihai Anitescu. Recovering sparse DFT from missing signals via interior point method on GPU. 2025. ⟨hal-05434700⟩ 
  • François Pacaud, Armin Nurkanović, Anton Pozharskiy, Alexis Montoison, Sungho Shin. An Augmented Lagrangian Method on GPU for Security-Constrained AC Optimal Power Flow. 2025. ⟨hal-05434667⟩ 
  • François Pacaud. Sensitivity analysis for parametric nonlinear programming: A tutorial. 2025. ⟨hal-05434699⟩ 
  • François Pacaud, Sungho Shin. GPU-accelerated dynamic nonlinear optimization with ExaModels and MadNLP François Pacaud and Sungho Shin. Conference on Decision and Control, IEEE, Dec 2024, Milano, Italy. ⟨hal-04875892⟩ 
  • Alexis Montoison, François Pacaud, Michael Saunders, Sungho Shin, Dominique Orban. MadNCL: A GPU Implementation of Algorithm NCL for Large-Scale, Degenerate Nonlinear Programs. 2025. ⟨hal-05434698⟩https://hal.science/hal-05403642v1

Also to be discovered