Organizations: Lanzhou University of Technology, School of Automation and Electrical Engineering, Lanzhou, 730050, Gansu, China · Hiroshima University, Graduate School of Advanced Science and Engineering, Higashihiroshima, 739-0046, Hiroshima, Japan · Nanjing University - Suzhou Campus, School of Robotics and Automation, Suzhou, 215163, Jiangsu, China · Princeton University, Department of Mechanical and Aerospace Engineering, Princeton, 08544, New Jersey, United States
This paper addresses the issue of self-intersecting trajectories (in phase space) in industrial reduced-order modeling and proposes the Latent-Augmented Neural Ordinary Differential Equations (LA-NODEs) framework. From the perspective of artificial intelligence, the proposed method augments conventional neural ordinary differential equations to enhance model expressiveness, enabling the representation of conflicting vector fields that may arise in reduced-order systems, thereby improving learning accuracy. Through theoretical analysis, the underlying mechanism of the framework is established, and a condition for determining the minimum required augmentation dimension is derived. From the perspective of engineering applications, the effectiveness of the proposed method is validated on the reduced-order system of two representative industrial models, namely an interior permanent magnet synchronous motor (IPMSM) drive and a distributed energy system (DES). Experimental results demonstrate that the proposed method can recover system features that are difficult to capture using conventional approaches and achieve superior performance in terms of prediction accuracy and modeling fidelity, thereby providing an effective approach for high-precision data-driven modeling of complex industrial systems.
Figures & tables
Figure 1: Comparison of learned and true trajectories for example system.
Figure 2: Architecture of the LA-NODE model, including primary neural dynamics (Block 1, Eq. ( 8a )), latent neural dynamics (Block 2, Eq. ( 8b )), adjoint primary neural dynamics (Block 3, Eq. ( 11a )), and adjoint latent neural dynamics (Block 4, Eq. ( 11b )). ayt and azt represent the adjoint state defined as Eq. ( 10 ). Blue lines denote forward propagation, red lines denote adjoint (loss-driven) backward propagation, and green lines denote parameter updates as shown in Eq. ( 12 ).
Parameter
Symbol
Value [Unit]
Rated power
P
150 [W]
Rated speed
Nr
1000 [rpm]
Rated load
TL
1.2 [N ⋅ m]
Rated voltage
Ur
24 [V]
Rated current
Ir
7.5 [A]
Number of pole pairs
np
4 [–]
Table 1: Key parameters of the IPMSM system
Parameter
Symbol
Value [Unit]
Maximum battery capacity
Ebat,max
5000 [kWh]
Charge/discharge efficiency
ηbat
0.95
Inertia constant
M
10 [s]
Damping coefficient
D
1 [p.u.]
Rated angular frequency
ω0
100π [rad/s]
PV dynamic time constant
TPV
5 [s]
Table 2: Key parameters of the DES
Figure 3: Comparison of trajectory-wise errors ( Error=(id,true(i)−id,pred(i))2+(iq,true(i)−iq,pred(i))2 ) across four frameworks over the full dataset for IPMSM. Index represents the data point position in the IPMSM dataset.
Figure 4: Performance comparison between CSODE and LA-NODE based on IPMSM trajectories: (a) Phase portraits, (b) Time-domain responses.
Figure 5: Comparison of MSE Loss across different frameworks with epochs for IPMSM: (a) Training set, (b) Testing set.
Figure 6: Comparison of state variables for DES (the color variation indicates the progression of time): (a) LA-NODE, (b) ANODE, (c) CSODE, (d) Transformer, for x3 ; (e) LA-NODE, (f) ANODE, (g) CSODE, (h) Transformer, for x4 .
Model
MAE
RMSE
Max Error
Std Error
LA-NODE (ours)
0.0075 ± 0.0037
0.0097 ± 0.0049
0.0270 ± 0.0138
0.0061 ± 0.0033
CSODE
0.1771 ± 0.0053
0.2266 ± 0.0046
0.5819 ± 0.0095
0.1404 ± 0.0008
Transformer
0.2625 ± 0.0404
0.3517 ± 0.0422
0.8586 ± 0.0498
0.2333 ± 0.0240
ANODE
0.0171 ± 0.0085
0.0236 ± 0.0114
0.0858 ± 0.0484
0.0161 ± 0.0083
Table 3: Comparison of four dynamic modeling frameworks for DES
Figure 7: Comparison of MSE Loss for the testing set across four frameworks for DES.
Figure 8: Comparison of MSE Loss (on logarithmic scale) for training set and testing set across different augmentation dimensions.
k
1
2
3
4
5
Time (sec)
637.2
638.2
648.9
684.0
740.2
Table 4: Average total training time (sec) over 5 runs for 700 epochs, varying k .
Figure 9: Trade-off analysis between computational cost and test loss under different augmentation dimensions.