Codes, scripts & models

Codes and models

MATLAB scripts and models

• COMSOL Multiphysics — I have developed models of magnetic fields, flows, convection and aero-acoustics, including an axisymmetric ZoRo acoustics model (v.5.4) with non-zero azimuthal wavenumbers.

• FELINS (Free-surface Ellipsoids with LINear State) — a MATLAB code for unsteady uniform-vorticity flows in deforming, free-surface ellipsoids. Forcing can include self-gravity, tides, linear elasticity and the Lorentz force from a uniform imposed magnetic field in the inductionless limit ; for modest deformations, surface tension can be included. Benchmarks include rotating charged drops and tidally deformed self-gravitating fluids.

• SIF² (Spectral code for Inviscid Free Surface Irrotational Flow) — a MATLAB code for three-dimensional sloshing in a rectangular tank. It expands the velocity potential in the tank’s natural sloshing modes and solves the fully nonlinear free-surface problem without linearising its evolution. The spectral expansion follows Fenton and Rienecker (1982), while the implementation was inspired by D. Le Touzé’s thesis code (Nantes, 2003).

• FLIPPER (FLows at ImPosed Precession within Ellipsoids in Rotation) — a MATLAB script for uniform-vorticity flows in precessing, rotating ellipsoids. It was released with Cébron, Fluid Dynamics Research (2015). It implements several formulations from Busse, Noir and collaborators as well as more efficient recasts of the bulk-flow equations.

• LEAFS (Lagrangian Extended Analysis of Flow Stability) — a MATLAB script for short-wavelength stability analysis along the trajectories of a three-dimensional, time-dependent inviscid flow. Localised perturbations evolve through a WKB/Floquet system : three equations each for the trajectory and wavenumber, plus nine for the monodromy matrix. The resulting growth rate characterises local instability along a pathline.

• TREFLES (TRiaxial Ellipsoid Flows and LinEar Stability) — a MATLAB script for uniform-vorticity flows driven by imposed motions of a triaxial ellipsoid. It covers precession, nutation and libration, separately or together. The rotation axis and the ellipsoid figure axes can move independently. Polynomial perturbations up to degree five give linear growth rates and account for confinement by the ellipsoidal boundary.

Codes developed by collaborators

• ToCCo — Rémy Monville’s code for local topographic coupling in rotating magnetohydrodynamic flows. Related paper : Monville, Cébron and Jault (2025).
• SHINE, SWAN and SIREN — developed by Jérémie Vidal to study hydromagnetic modes and flow stability in ellipsoids.
• XSHELLS — Nathanaël Schaeffer’s spectral code for rotating flows and magnetic fields in spherical geometry.

Data and related material

• Data and code for Personnettaz et al. (2026), on viscous and ohmic damping of inner-core oscillations, by Paolo Personnettaz and coauthors.
• Numerical tools used in THEIA.