Search

Neural operators have been validated as promising deep surrogate models for solving partial differential equations (PDEs). Despite the critical role of boundary conditions in PDEs, however, only a limited number of neural operators robustly enforce these conditions. In this paper we introduce semi-periodic Fourier neural operator (SPFNO), a novel spectral operator learning method, to learn the target operators of PDEs with non-periodic BCs. This method extends our previous work (arXiv:2206.12698), which showed significant improvements by employing enhanced neural operators that precisely satisfy the boundary conditions. However, the previous work is associated with Gaussian grids, restricting comprehensive comparisons across most public datasets. Additionally, we present numerical results for various PDEs such as the viscous Burgers' equation, Darcy flow, incompressible pipe flow, and coupled reactiondiffusion equations. These results demonstrate the computational efficiency, resolution invariant property, and BC-satisfaction behavior of proposed model. An accuracy improvement of approximately 1.7X-4.7X over the non-BC-satisfying baselines is also achieved. Furthermore, our studies on SOL underscore the significance of satisfying BCs as a criterion for deep surrogate models of PDEs.

By learning the mappings between infinite function spaces using carefully designed neural networks, the operator learning methodology has exhibited significantly more efficiency than traditional methods in solving complex problems such as differential equations, but faces concerns about their accuracy and reliability. To overcomes these limitations, combined with the structures of the spectral numerical method, a general neural architecture named spectral operator learning (SOL) is introduced, and one variant called the orthogonal polynomial neural operator (OPNO), developed for PDEs with Dirichlet, Neumann and Robin boundary conditions (BCs), is proposed later. The strict BC satisfaction properties and the universal approximation capacity of the OPNO are theoretically proven. A variety of numerical experiments with physical backgrounds show that the OPNO outperforms other existing deep learning methodologies, as well as the traditional 2nd-order finite difference method (FDM) with a considerably fine mesh (with the relative errors reaching the order of 1e-6), and is up to almost 5 magnitudes faster than the traditional method.

In this paper, we are concerned with arbitrarily high-order momentum-preserving and energy-preserving schemes for solving the generalized Rosenau-type equation, respectively. The derivation of the momentum-preserving schemes is made within the symplectic Runge-Kutta method, coupled with the standard Fourier pseudo-spectral method in space. Then, combined with the quadratic auxiliary variable approach and the symplectic Runge-Kutta method, together with the standard Fourier pseudo-spectral method, we present a class of high-order mass- and energy-preserving schemes for the Rosenau equation. Finally, extensive numerical tests and comparisons are also addressed to illustrate the performance of the proposed schemes., Comment: 23 pages, 31 figures

Much recent work has addressed the solution of a family of partial differential equations by computing the inverse operator map between the input and solution space. Toward this end, we incorporate function-valued reproducing kernel Hilbert spaces in our operator learning model. We use neural networks to parameterize Hilbert-Schmidt integral operator and propose an architecture. Experiments including several typical datasets show that the proposed architecture has desirable accuracy on linear and nonlinear partial differential equations even with a small amount of data. By learning the mappings between function spaces, the proposed method can find the solution of a high-resolution input after learning from lower-resolution data., Comment: 14 pages, 8 figures, 4 tables

In this paper, a class of arbitrarily high-order linear momentum-preserving and energy-preserving schemes are proposed, respectively, for solving the regularized long-wave equation. For the momentum-preserving scheme, the key idea is based on the extrapolation/prediction-correction technique and the symplectic Runge-Kutta method in time, together with the standard Fourier pseudo-spectral method in space. We show that the scheme is linear, high-order, unconditionally stable and preserves the discrete momentum of the system. For the energy-preserving scheme, it is mainly based on the energy quadratization approach and the analogous linearized strategy used in the construction of the linear momentum-preserving scheme. The proposed scheme is linear, high-order and can preserve a discrete quadratic energy exactly. Numerical results are addressed to demonstrate the accuracy and efficiency of the proposed scheme., Comment: 26 pages, 52 figures

A novel class of high-order linearly implicit energy-preserving integrating factor Runge-Kutta methods are proposed for the nonlinear Schr\"odinger equation. Based on the idea of the scalar auxiliary variable approach, the original equation is first reformulated into an equivalent form which satisfies a quadratic energy. The spatial derivatives of the system are then approximated with the standard Fourier pseudo-spectral method. Subsequently, we apply the extrapolation technique/prediction-correction strategy to the nonlinear terms of the semi-discretized system and a linearized energy-conserving system is obtained. A fully discrete scheme is gained by further using the integrating factor Runge-Kutta method to the resulting system. We show that, under certain circumstances for the coefficients of a Runge-Kutta method, the proposed scheme can produce numerical solutions along which the modified energy is precisely conserved, as is the case with the analytical solution and is extremely efficient in the sense that only linear equations with constant coefficients need to be solved at every time step. Numerical results are addressed to demonstrate the remarkable superiority of the proposed schemes in comparison with other existing structure-preserving schemes., Comment: 28 pages, 49 figures

We propose and analyze two regularized finite difference methods for the logarithmic Klein-Gordon equation (LogKGE). Due to the blowup phenomena caused by the logarithmic nonlinearity of the LogKGE, it is difficult to construct numerical schemes and establish their error bounds. In order to avoid singularity, we present a regularized logarithmic Klein-Gordon equation (RLogKGE) with a small regularized parameter $0<\varepsilon\ll1$. Besides, two finite difference methods are adopted to solve the regularized logarithmic Klein-Gordon equation (RLogKGE) and rigorous error bounds are estimated in terms of the mesh size $h$, time step $\tau$, and the small regularized parameter $\varepsilon$. Finally, numerical experiments are carried out to verify our error estimates of the two numerical methods and the convergence results from the LogKGE to the RLogKGE with the linear convergence order $O(\varepsilon)$.

We present and analyze two regularized finite difference methods which preserve energy of the logarithmic Klein-Gordon equation (LogKGE). In order to avoid singularity caused by the logarithmic nonlinearity of the LogKGE, we propose a regularized logarithmic Klein-Gordon equation (RLogKGE) with a small regulation parameter $0<\varepsilon\ll1$ to approximate the LogKGE with the convergence order $O(\varepsilon)$. By adopting the energy method, the inverse inequality, and the cut-off technique of the nonlinearity to bound the numerical solution, the error bound $O(h^{2}+\frac{\tau^{2}}{\varepsilon^{2}})$ of the two schemes with the mesh size $h$, the time step $\tau$ and the parameter $\varepsilon$. Numerical results are reported to support our conclusions., Comment: arXiv admin note: text overlap with arXiv:2006.08079

We present and analyze a series of conservative diagonally implicit Runge--Kutta schemes for the nonlinear Schr\"odiner equation. With the application of the newly developed invariant energy quadratization approach, these schemes possess not only high accuracy , high order convergence (up to fifth order) and efficiency due to diagonally implicity but also mass and energy conservative properties. Both theoretical analysis and numerical experiments of one- and two-dimensional dynamics are carried out to verify the invariant conservative properties, convergence orders and longtime simulation stability.

In this paper, we construct stochastic symplectic Runge--Kutta (SSRK) methods of high strong order for Hamiltonian systems with additive noise. By means of colored rooted tree theory, we combine conditions of mean-square order 1.5 and symplectic conditions to get totally derivative-free schemes. We also achieve mean-square order 2.0 symplectic schemes for a class of second-order Hamiltonian systems with additive noise by similar analysis. Finally, linear and non-linear systems are solved numerically, which verifies the theoretical analysis on convergence order. Especially for the stochastic harmonic oscillator with additive noise, the linear growth property can be preserved exactly over long-time simulation., Comment: 23 pages, 5 figures

In this paper, we consider the numerical methods preserving single or multiple conserved quantities, and these methods are able to reach high order of strong convergence simultaneously based on some kinds of projection methods. The mean-square convergence orders of these methods under certain conditions are given, which can reach order 1.5 or even 2 according to the supporting methods embedded in the projection step. Finally, three numerical experiments are taken into account to show the superiority of the projection methods.