Cryogenics · NMPC · state estimation
Real-time NMPC for the LHC superfluid-helium circuit
- Context
- CERN, Geneva — researcher and doctoral student (2007–2013); PhD, University of Valladolid (2015).
- Problem
- Keep the superconducting magnets of the LHC at their operating temperature of 1.9 K. Above 2.1 K an interlock de-energizes the magnets; above 2.16 K the helium leaves the superfluid state and a magnet can quench.
- Dynamics and constraints
- Strongly nonlinear, distributed thermohydraulics with dead time, inverse response, actuator saturation and strong coupling — and only a few measurements along a long cryogenic circuit.
- Existing approach
- Conventional PI control, tuned conservatively to stay clear of the temperature limits.
- Solution
- A first-principles thermohydraulic model simplified for real-time use; nonlinear moving-horizon estimation of the unmeasured states (5 thermodynamic states from 3 pressure sensors and temperature, at 1 Hz); output-feedback economic NMPC with the temperature limits as hard constraints.
- Implementation and validation
- Semi-automatic code generation for the NMPC of a stiff, distributed-parameter system. About 1 s per estimation and 7–14 s per optimization. The controller was tested at CERN.
- Result
- After the same perturbation, the NMPC stabilized the bath temperature about ten times faster than the legacy PI control, without violating the 2.1 K limit.
- Lessons
- The estimator decides how good the controller can be. A simplified model that runs in real time beats a complete one that does not. Hard limits belong in the optimization problem, not in conservative tuning.
- Evidence
-
- Noga et al., IEEE CDC 2014 — moving-horizon estimation and control
- Noga et al., IFAC NMPC 2015 — NMPC for superfluid helium cryogenics
- Noga, PhD thesis, University of Valladolid, 2015
- Noga & Ohtsuka, ICPC 2011 — code generation for NMPC of a stiff distributed-parameter system
- CERN recommendation letter (Geneva, 10.03.2009): “I recommend Mr. Noga without any reservation.”
Have a similar control or optimization problem? Book a 30-min technical fit call →