Table of Contents
Fetching ...

Investigating the Surrogate Modeling Capabilities of Continuous Time Echo State Networks

Saakaar Bhatnagar

TL;DR

This work investigates CTESNs as surrogate models for stiff ODEs, focusing on how the reservoir-to-output projection affects generalization. It systematically compares linear and nonlinear projections, showing that NLPCTESN outperforms LPCTESN across benchmark problems (Robertson, Sliding Basepoint, POLLU) and across parameterizations of rates and initial conditions, with further gains from polynomial-augmented k-NN RBF interpolation. The CTESN surrogates deliver accurate predictions over challenging timescales and offer speedups of up to two orders of magnitude over traditional ODE solvers, enabling real-time prediction and potential digital-twin applications in engineering systems. The study provides practical guidance on algorithmic choices and hyper-parameter settings for CTESNs in stiff dynamical systems and motivates future physics-informed extensions for broader, high-dimensional problems.

Abstract

Continuous Time Echo State Networks (CTESNs) are a promising yet under-explored surrogate modeling technique for dynamical systems, particularly those governed by stiff Ordinary Differential Equations (ODEs). A key determinant of the generalization accuracy of a CTESN surrogate is the method of projecting the reservoir state to the output. This paper shows that of the two common projection methods (linear and nonlinear), the surrogates developed via the nonlinear projection consistently outperform those developed via the linear method. CTESN surrogates are developed for several challenging benchmark cases governed by stiff ODEs, and for each case, the performance of the linear and nonlinear projections is compared. The results of this paper demonstrate the applicability of CTESNs to a variety of problems while serving as a reference for important algorithmic and hyper-parameter choices for CTESNs

Investigating the Surrogate Modeling Capabilities of Continuous Time Echo State Networks

TL;DR

This work investigates CTESNs as surrogate models for stiff ODEs, focusing on how the reservoir-to-output projection affects generalization. It systematically compares linear and nonlinear projections, showing that NLPCTESN outperforms LPCTESN across benchmark problems (Robertson, Sliding Basepoint, POLLU) and across parameterizations of rates and initial conditions, with further gains from polynomial-augmented k-NN RBF interpolation. The CTESN surrogates deliver accurate predictions over challenging timescales and offer speedups of up to two orders of magnitude over traditional ODE solvers, enabling real-time prediction and potential digital-twin applications in engineering systems. The study provides practical guidance on algorithmic choices and hyper-parameter settings for CTESNs in stiff dynamical systems and motivates future physics-informed extensions for broader, high-dimensional problems.

Abstract

Continuous Time Echo State Networks (CTESNs) are a promising yet under-explored surrogate modeling technique for dynamical systems, particularly those governed by stiff Ordinary Differential Equations (ODEs). A key determinant of the generalization accuracy of a CTESN surrogate is the method of projecting the reservoir state to the output. This paper shows that of the two common projection methods (linear and nonlinear), the surrogates developed via the nonlinear projection consistently outperform those developed via the linear method. CTESN surrogates are developed for several challenging benchmark cases governed by stiff ODEs, and for each case, the performance of the linear and nonlinear projections is compared. The results of this paper demonstrate the applicability of CTESNs to a variety of problems while serving as a reference for important algorithmic and hyper-parameter choices for CTESNs
Paper Structure (14 sections, 24 equations, 6 figures, 11 tables, 2 algorithms)

This paper contains 14 sections, 24 equations, 6 figures, 11 tables, 2 algorithms.

Figures (6)

  • Figure 1: Depiction of a standard Echo State Network.
  • Figure 2: Figures show the time history for $y_{2}$ for a test parameter set (P2 from Table \ref{['table:rate_test_values']}). Each trial mentioned below is referenced from Table \ref{['table:robertson_rate_errors']}. The NLPCTESN prediction is the best out of all models. (\ref{['fig:nx_50_nd_50_no_poly_rate ']}) Trial 1 (\ref{['fig:nx_50_nd_50_L_no_neigh_rate ']}) Trial 2 (\ref{['fig:nx_50_nd_500_L_rate ']}) Trial 3 (\ref{['fig:nx_50_nd_50_L_rate ']}) Trial 4 (\ref{['fig:nx_50_nd_50_NL_rate ']}) Trial 5
  • Figure 3: Figures show time history for $y_{2}$ for a test parameter set (P4 from Table \ref{['table:init_test_values']}). Each trial mentioned below is referenced from Table \ref{['table:robertson_init_errors']}. The NLPCTESN performs the best out of all models. (\ref{['fig:nx_500_nd_50_L_no_neigh ']}) Trial 1 (\ref{['fig:nx_50_nd_50_L_no_neigh ']}) Trial 2 (\ref{['fig:nx_500_nd_50_L ']}) Trial 3 (\ref{['fig:nx_50_nd_50_L ']}) Trial 4 (\ref{['fig:nx_50_nd_50_NL ']}) Trial 5
  • Figure 4: Sliding Basepoint system for collision modeling.
  • Figure 5: Solution $v_{2}$ for several tests with different test parameters, the average error of which is shown in Table \ref{['table:crash_errors']}. The left and right columns show results from the linear and nonlinear projections respectively. Table \ref{['table: crash_test_values']} defines the values P3-P5 mentioned below. (\ref{['fig:v2_2_LPCTESN ']}) LPCTESN - P3 (\ref{['fig:v2_2_NLPCTESN ']}) NLPCTESN - P3 (\ref{['fig:v2_3_LPCTESN ']}) LPCTESN - P4 (\ref{['fig:v2_3_NLPCTESN ']}) NLPCTESN - P4 (\ref{['fig:v2_4_LPCTESN ']}) LPCTESN - P5 (\ref{['fig:v2_4_NLPCTESN ']}) NLPCTESN - P5
  • ...and 1 more figures