Score
Design, implement, and evaluate compact and efficient Fourier Neural Operator (FNO) architectures that minimize parameter count and inference latency while preserving predictive accuracy and the ability to generalize across finer mesh or resolution grids. Analyze and optimize trade-offs among model size, computational cost, and accuracy by applying parameter-reduction and efficiency techniques (e.g., streamlined spectral convolutions and architecture modifications) to produce deployable FNO models.
This work addresses aliasing errors introduced in the discrete implementation of Fourier Neural Operators (FNOs), which remain theoretically unquantified despite their practical significance. Specifically, the discrepancy between the continuous FNO formulation and its grid-based discretization has not been systematically characterized, and isolating this discretization error from other sources—such as approximation and optimization errors—remains challenging. Method: Leveraging tools from Fourier analysis, numerical functional analysis, and FFT theory, we derive an explicit algebraic convergence rate for the discretization error with respect to grid resolution and establish its quantitative dependence on the Sobolev regularity of the input function. Contribution/Results: We provide the first verifiable upper bound on this error, revealing intrinsic trade-offs among resolution, input smoothness, and model stability. Numerical experiments confirm the theoretical prediction: for inputs with higher Sobolev regularity, the error decays significantly under mesh refinement—thereby furnishing a rigorous foundation for reliable discrete FNO design.
Fourier Neural Operators (FNOs) show promise for PDE modeling but suffer from insufficient theoretical grounding, implementation ambiguities, and sensitivity to discretization—undermining reliability in practice. Method: We systematically formalize FNO’s operator-theoretic foundations and clarify its intrinsic frequency-domain modeling mechanism, correcting the common misconception that “frequency domain” equates to mere Fourier transforms. We propose a discretization-invariant implementation paradigm based on learnable frequency-domain filters. Leveraging NeuralOperator 2.0.0, we deliver a modular, theory-aligned open-source implementation supporting multi-resolution generalization and global correlation modeling. Contribution/Results: Experiments demonstrate substantial improvements in stability, generalization, and deployment efficiency across challenging PDE systems—including Navier–Stokes and Darcy flow—establishing a robust, accessible infrastructure for operator learning in scientific machine learning.
研究通过区分稀疏表示、存储参数、理论操作数和测量运行时间,探索了稀疏Fourier神经算子的不同实现路径,以解决模型大小和推理成本问题。
Fourier Neural Operators (FNOs) excel at modeling global frequency-domain patterns but struggle to capture critical local spatial structures in partial differential equations (PDEs), limiting accuracy. To address this, we propose Conv-FNO—a hybrid architecture that integrates a lightweight CNN preprocessor to explicitly extract local features, coupled with frequency-domain parameterized convolutions and an adaptive grid scaling mechanism to ensure resolution-invariant feature alignment. Theoretical analysis demonstrates that local feature enhancement improves generalization. Evaluated on multiple PDE benchmarks, Conv-FNO reduces mean error by 27–41% over standard FNOs and state-of-the-art hybrid models, while retaining FNO-level inference efficiency. Our key contribution is the first seamless, scale-robust coupling of CNNs and FNOs—establishing a new paradigm for PDE surrogates that jointly achieves local sensitivity and global spectral modeling.
This work addresses the limitations of Fourier neural operators in terms of generalization accuracy and computational efficiency for mappings in function spaces. To overcome these challenges, the authors propose an efficient architecture that uniquely integrates rank-1 lattice sampling with hyperbolic cross truncation in the frequency domain. This combination reduces high-dimensional Fourier transforms to one-dimensional fast Fourier transforms and introduces a structured lattice-based training set in parameter space. The resulting method substantially decreases the number of model parameters, spatial sampling points, and required training samples. Evaluated on elliptic partial differential equations over toroidal domains, the approach achieves higher approximation accuracy while maintaining lower computational complexity compared to existing methods.
Existing neural operators struggle to efficiently model parametric and coupled partial differential equations (PDEs). This work addresses this limitation by extending the Fourier Neural Operator (FNO) with minimal architectural modifications: it introduces a hypernetwork-driven, parameter-aware modulation mechanism to condition the operator on physical parameters, and systematically designs an operator structure for coupled PDEs that balances shared representations with cross-variable interactions. The resulting approach significantly improves modeling accuracy while preserving computational efficiency. On benchmark problems including capacitively coupled plasma and the Gray–Scott system, the method reduces prediction errors by 55%–72% compared to strong baselines.
This study addresses the limitations of Fourier Neural Operators (FNOs) in learning high-frequency components due to frequency truncation, as well as the challenges in coordinate encoding design. To this end, we propose the CAFE+FNO framework, which introduces a novel CAFE+ mechanism that fuses Fourier-Chebyshev features via the Hadamard product to enhance frequency-domain interactions. By integrating implicit neural representations, parallel affine branches, and shared-kernel MLPs, the method generates Fourier kernels through explicit feature composition, enabling efficient PDE solving with model parameters independent of the number of modes. Experiments across five PDE benchmarks demonstrate that the proposed approach significantly outperforms existing FNO variants, fully validating its advantages in bandwidth learnability.
This study addresses the limitation of Fourier Neural Operators in modeling non-periodic and local fine-scale structures due to their fixed spectral bases and hard truncation. To overcome this, we propose SAFDNO, which introduces stochastic adaptive Fourier decomposition theory to construct input-adaptive Takenaka-Malmquist orthogonal systems that replace predefined Fourier modes. Furthermore, a neural pole predictor is designed to amortize the stochastic selection process, enabling learned filtering on the analytic branch coefficients within Hardy spaces. Experimental results demonstrate that SAFDNO achieves superior performance across nine PDE benchmarks, significantly improving solution accuracy for equations such as Darcy flow and enhancing zero-shot super-resolution capabilities.
This work addresses the challenge of balancing accuracy and computational efficiency in modeling two-dimensional Rayleigh–Bénard convection by proposing a lightweight Fourier Neural Operator (FNO) architecture based on temporal increment prediction. Instead of predicting the full solution directly, the model forecasts the state increment, substantially reducing its complexity. The resulting network contains only 314k parameters (1.26 MB) and achieves a single inference time of 7 milliseconds. While maintaining accuracy comparable to existing methods, the approach significantly lowers computational overhead, offering a novel pathway toward efficient, high-fidelity simulations in fluid dynamics.
This work addresses the severe communication bottleneck in distributed Fourier Neural Operators (FNOs) at high resolutions, where spectral layers require frequent all-to-all communication. To mitigate this, the authors propose Distributed Truncated Spectral Transform (DTST), which computes only a small set of critical frequency modes locally on each GPU and aggregates them with minimal communication, achieving equivalence to truncated FFT while drastically reducing communication overhead during both training and inference. DTST uniquely integrates local discrete Fourier transforms with efficient collective communication, unifying spatial data parallelism and spectral weight model parallelism. Experiments demonstrate that, across 4–32 GPUs (up to 8 nodes), DTST accelerates forward propagation by 38–64× and training by 37×, reducing communication time from 97% to under 6% of total execution time, with scalability improving further as resolution increases.