[論文レビュー] Fast Jacobi Spectral Methods and Closure Approximations for the Homogeneous FENE Model of Complex Fluids

Runkai Feng, Jie Shen+1 more|arXiv (Cornell University)|Feb 8, 2026

Model Reduction and Neural Networks被引用数 0

ひとこと要約

要約: 論文は同質 FENE Fokker–Planck 方程式の2つの高速 Jacobi-球面調和関数スペクトル-Galerkin法を開発し、BDF2 時間安定性を証明し、閉包モデル（FENE-P、FENE-QE、FENE-QE-NN）をニューラルネットワーク強化と共に比較する。QE閉包の近似を効率化するニューラルネットワーク閉包も導入し、従来法との精度と効率を検証。

ABSTRACT

The Finitely Extensible Nonlinear Elastic (FENE) dumbbell model is a widely used mathematical model for complex fluids. Direct simulation of the FENE Fokker--Planck equation is computationally challenging due to high dimensionality and singularity of its potential. In this paper, we develop two fast Jacobi-Spherical Harmonic spectral methods for the spatially homogeneous FENE Fokker--Planck equation. These methods effectively resolve the singularity near the boundary by combining properly designed Jacobi polynomials with a weighted variational formulation. A semi-implicit backward differentiation formula of second-order (BDF2) is employed for time marching, and its energy stability is rigorously proved. The resulting linear algebraic system possesses a sparse structure and can be efficiently solved. Numerical results verify the spectral convergence and efficiency of the direct spectral solvers, establishing them as a reliable tool for generating reference solutions for challenging benchmark problems. Furthermore, to achieve an optimal trade-off between accuracy and efficiency, we compare several closure approximation models, including the industry workhorse Peterlin approximation (FENE-P), the quasi-equilibrium approximation (FENE-QE), and a novel neural network implementation for FENE-QE proposed in this paper (FENE-QE-NN). Numerical experiments in extensional and shear flows demonstrate the superior accuracy and efficiency of the proposed methods compared to traditional approaches.

研究の動機と目的

• 境界特異性と高次元性に起因する FENE Fokker–Planck 方程式の数値的課題に対処。
• 0+3 次元同質 FENE ダイナミクスの二つの高速 Jacobi-球面調和関数スペクトル-Galerkin法を開発。
• 半古典的な BDF2 スキームによる安定した時間積分と厳密な収束検証を提供。
• 閉包モデル（FENE-P、FENE-QE）を比較・実装し、ニューラルネットワーク支援の QE 閉包を導入して精度と効率を向上。

提案手法

• 境界特異性の処理を伴う FENE Fokker–Planck 方程式の加重弱形を定式化し、(1-|q|^2)^s 変換で境界近傍を扱う。
• 配置空間における高速 Jacobi-Galerkin スペクトル離散化を適用し、角度変数には球面調和関数、半径方向には Jacobi 多項式を用いる（JG1 および JGinf 変種）。
• 対流項を明示的な線形外挿と組み合わせた半陰的な BDF2 時間ステップを用いて解法をデカップリング・安定化。
• 時間離散化のコercivity とエネルギー安定性を、(h, q·∇h/(1-|q|^2))_s をノルム項と relating する重要な等式を用いて証明し、安定性条件を提示。
• モーメント閉包モデル（FENE-P、FENE-QE）を導出・実装し、構成テンソル固有値をラグランジュ乗数へ写像する PLA およびニューラルネットワーク（FENE-QE-NN）アプローチを開発。
• スペクトル解法と閉包の精度と効率を評価するためのひずみ・せん断流における数値ベンチマークを提供。

Figure 1 : Spectral convergence of the relative $L_{\omega}^{2}$ error with respect to the polynomial degree $N$ .

実験結果

リサーチクエスチョン

• RQ1 境界特異性がある場合でも、速い Jacobi-球面調和関数スペクトル法は 0+3 FENE Fokker–Planck 方程式を正確に解けるのか。
• RQ2 拡張およびせん断流に対する FENE-P、FENE-QE、およびニューラルネットワーク強化 QE 閉包の精度と計算効率はどう比較されるのか。

主な発見

• 二つの高速 Jacobi スペクトル法はスペクトル収束を達成し、同質 FENE Fokker–Planck 方程式の効率的な直接解法を可能にする。
• BDF2 時間離散化は指定条件下でエネルギー安定性を証明され、頑健な時間積分を保証。
• FENE-P および FENE-QE 閉包を分析・対比し、 QE は高忠実度を提供する一方、PLA および NN 増強 QE は効率的な代替案となる。
• ニューラルネットワーク（FENE-QE-NN）は構成テンソル固有値から QE ラグランジュ乗数への写像を学習し、テーブルルックアップなしで閉包の効率を向上させる。
• 伸長流およびせん断流のテストは、提案手法が従来法よりも精度と計算コストの点で優れていることを示す。

Figure 2 : Training diagnostics of the FENE-QE-NN constitutive model. (a) Loss history (MSE) showing the transition from Adam to L-BFGS. (b) Parity plot of the predicted vs. true Lagrange multipliers $\lambda$ on the validation set.