The Shapes a Machine Draws

Recalibration-invariant relational atoms of multi-physics operation for industrial telemetry
Sehaj Singh
Independent Researcher

A machine under load moves along a coupled thermal, mechanical, and electrical trajectory. We propose nine scalar descriptors — atoms — computed from the joint trajectory of every channel pair. Eight are built from within-unit ranks and are exactly invariant to any strictly monotone recalibration of each channel separately. The ninth, the signed Lévy area, captures the arrow of time and is the unique atom that changes sign under time reversal. On seven real industrial systems (309 records) the atlas identifies the source system at 99.4% accuracy (CI [97.7, 100.0]) against 96.1% for per-channel marginals. Under re-instrumentation the atlas holds at 99.0% while marginals fall to 80.3%. In controlled simulated fleets of 48 machines, an autoencoder trained on raw telemetry beats the atlas on clean data by 1.5× but collapses to chance under re-instrumentation; the atlas is numerically unchanged. The method requires a minimum of approximately 10³ samples per window and does not improve short-horizon forward prediction, where the identity is redundant with the state history. Its value is in identification, auditing, and transfer rather than in accuracy.

The idea: route on geometry, not on raw state

A servo joint accelerating a payload dissipates ohmic heat into a winding whose rising temperature derates the torque available to the next command — on a thermal time constant three decades slower than the control loop that issued it. Watch any one channel and you see a line wandering. Watch two against each other and you see a shape, and the shape is where the coupling lives.

Each channel pair gets nine atoms. The whole array over all pairs is one unit's atlas. Nine quantiles of each atom across pairs gives a signature of fixed length whatever the channel count — which is how a 12-channel motor and a 24-channel turbofan can be compared at all.

Even · rank-based

Spearman ρ

Monotone association. How strongly one channel tracks another, regardless of scale.

Even · rank-based

Distance correlation

Functional strength. Detects any dependence, linear or not. Zero iff independent.

Even · rank-based

Hoeffding D

Nonlinearity gap. How far the joint distribution deviates from independence.

Even · rank-based

Copula asymmetry

Whether dependence is stronger in one tail than the other.

Odd · time-ordered

Signed Lévy area

The arrow of time. Flips sign under time reversal. The one atom a correlation matrix cannot reach.

Even · rank-based

Jump share

Fraction of trajectory spent in jumps rather than smooth motion.

Even · rank-based

Timescale ratio

Power distribution across frequencies. Separates slow thermal drift from fast oscillation.

Even · rank-based

Support occupancy

How much of the rank-plane the trajectory actually visits.

Even · rank-based

Scaling exponent

Hurst-like. Whether the trajectory is persistent, anti-persistent, or a random walk.

Recalibration invariance is exact

Replace a sensor, rescale a channel, apply any strictly increasing warp to each channel independently: the atlas does not move. This is not a trained property — it follows from the construction. Eight of nine atoms are computed on within-unit ranks, and monotone transforms preserve ranks exactly. Verified to machine precision (0.00×10⁰ change) on all three robotic platforms.

Nine atoms and recalibration invariance
Fig 1 — The atoms

Four exemplar channel pairs showing different coupling geometries. The full atlas as a heatmap. Under monotone warp: unchanged to machine precision.

Time reversal parity
Fig ED1 — Parity under time reversal

Eight atoms overlap perfectly after reversal (even sector). The Lévy area flips sign (odd sector) to 4.7×10⁻¹⁴.

Seven real industrial systems, 309 records

Nothing simulated except the recalibration warp. The atlas identifies which system a record came from at 99.4% against 96.1% for marginals and 86.4% for Spearman correlation. Under re-instrumentation: atlas 99.0% vs marginals 80.3%.

Turbomachinery
Gas turbine
Channels11
Years5
Real driftYes
Fluid + thermal
Hydraulic rig
Channels17
Fault labels4
PLC-controlledYes
Power grid
Synchrophasor PMU
Channels~24
Grid infraYes
Power grid
Transformer (ETT)
Channels7
Oil temp + loadReal
Turbomachinery
Wind farm SCADA
Channels4
ThinnestSystem
Motor bench
PMSM 52 kW
Channels12
Sessions40/69
Open loopRich
Pneumatic + pumping
SKAB + MetroPT
Channels8–15
AnomalyBenchmark
Seven-system benchmark
Fig 2 — Seven-system benchmark

Atlas 99.4% (CI [97.7, 100.0]) vs marginals 96.1%. Under re-instrumentation: atlas 99.0% vs marginals 80.3%.

Main results
Fig 3 — Atlas vs learned baselines

AE wins on clean data by 1.5×. Under warp: AE collapses to chance, atlas unchanged.

ComponentAtlasAtlas CIMarginalsWarp (atlas)Warp (marginals)
Cooler96.4%[89.7, 98.7]94.7%96.4%73.6%
Valve97.8%[93.2, 99.5]98.9%97.8%82.1%
Pump73.6%[62.8, 82.4]81.6%73.6%55.5%
Accumulator89.5%[80.3, 95.1]90.4%89.5%68.4%
Hydraulic rig fault diagnosis. For relational faults (cooler, valve, accumulator): atlas matches marginals and is unmoved by warp. For level faults (pump): marginals stronger — the correct design rule.

The controlled fleet experiment

Simulated fleets of 48 machines per platform (UR5e, Panda, iiwa14), each with independent physical parameters. Three fleet types isolate the two axes of machine identity: level (magnitude only — atlas should fail), shape (coupling only — atlas should succeed), both.

Fleet control
Fig 4 — Fleet control, 3 platforms × 3 types

Atlas at chance on every level fleet. AE wins clean, collapses under warp. Atlas numerically identical before and after.

PlatformFleetAtlasAtlas warpedAE cleanAE warped
UR5elevel6.2% (ns)6.2%8.3%2.1%
UR5eshape10.4% (p=.003)10.4%4.2%2.1%
UR5eboth14.6% (p<.001)14.6%37.5%2.1%
Pandalevel2.1% (ns)2.1%10.4%2.1%
Pandashape20.8% (p<.001)20.8%41.7%0.0%
Pandaboth39.6% (p<.001)39.6%58.3%2.1%
iiwa14level0.0% (ns)0.0%0.0%0.0%
iiwa14shape6.2% (ns)6.2%4.2%2.1%
iiwa14both12.5% (p<.001)12.5%4.2%2.1%
Chance 2.1%. Bold = significant (p<0.05). The AE wins on clean data (58.3% vs 39.6% on Panda all) but collapses to 2.1% under warp. The atlas is numerically identical before and after.

The parity argument: why correlation cannot substitute

Reverse a record end to end. A lag becomes a lead. Any statistic built from the joint distribution of simultaneous values is unchanged — reversal only permutes the samples. The Lévy area satisfies \(\mathcal{R}(\mathrm{levy}) = -\mathrm{levy}\) to \(4.7 \times 10^{-14}\). The vocabulary splits into an even sector (8 atoms) and a one-dimensional odd sector (1 atom).

A machine's defining multi-physics behaviour is a lag with a direction — torque heating a winding that then derates the torque. That lag is precisely the content of the odd sector. Every descriptor in the standard relational toolkit lives in the even sector.

What the atlas carries — and what it does not

TargetAtlas R²Marginals R²Reading
Payload0.610.99Level property — atlas sees it weakly
Bearing friction0.36−5.49Atlas carries it; marginals diverge
Ambient temperature−0.050.99Atlas is blind to workload
Duty cycle0.040.82Atlas is blind to workload
The workload control: duty cycle and ambient temperature are properties of the task, not the machine. A representation that encodes them is a duty detector wearing the costume of a machine fingerprint.

What the atlas cannot do

Real drift: 200× less movement

The PMSM bench has 69 sessions from one machine. Splitting into early (1–34) and late (35–69) gives a real sensor-ageing proxy without synthetic warps. The atlas centroid shifts 0.63 while marginals shift 126.3 — a 200× difference — yet both identify early vs late at the same accuracy (6.1%).

The atlas tracks relational physics while being immune to the magnitude drift that dominates the marginals.

Intrinsic dimension survey

SystemRegimeChannelsLive dimd/n
Panda robotclosed, repeated143.190.23
UR5e robotclosed, repeated102.280.23
SKABclosed, free102.280.23
Hydraulic rigclosed, repeated143.550.25
Transformerclosed, free145.360.38
Gas turbineclosed, free114.260.39
Steel plantclosed, free72.730.39
Turbofan (sim)closed, free249.450.39
PMSM benchopen, rich125.720.48
Wind farmclosed, free42.630.66
Ratio groups by drive regime: 0.23 for closed-loop repeated, 0.39 for closed-loop free, 0.60 for open-loop rich. Below d/n ≈ 0.35, geometric methods have something to work with.

Reproduce

GitHub repo · Kaggle kernels · HF Space

git clone https://github.com/sehajr-singhs/operating-atlas
cd operating-atlas
pip install -r requirements.txt

# Generate atoms from a motor session
python -m atlas.atoms --input data/robot_ur5e.npz

# Run the 7-system benchmark on Kaggle
kaggle kernels push -p kaggle_kernel/ioo_real2

# Build the findings page
python atlas/build_page.py

All figures generated by matplotlib with Okabe-Ito colour-blind-safe palette at 300 DPI. Paper compiled with pdflatex (22 pages, 44 references, zero errors).