## Dilated convolution neural operator for multiscale partial differential equations

Bo Xu , Xinliang Liu , Lei Zhang a b c a Shanghai Jiao Tong University, School of Mathematical Sciences, Shanghai, 200240, China b King Abdullah University of Science and Technology, Computer, Electrical and Mathematical Science and Engineering Division, King Abdullah University of Science and Technology, Thuwal, 23955, Saudi Arabia

- c Shanghai Jiao Tong University, School of Mathematical Sciences, Institute of Natural Sciences, and MOE-LSC, Shanghai, 200240, China

## Abstract

This paper introduces a data-driven operator learning method for multiscale partial differential equations, with a particular emphasis on preserving high-frequency information. Drawing inspiration from the representation of multiscale parameterized solutions as a combination of low-rank global bases (such as low-frequency Fourier modes) and localized bases over coarse patches (analogous to dilated convolution), we propose the Dilated Convolutional Neural Operator (DCNO). The DCNO architecture effectively captures both high-frequency and low-frequency features while maintaining a low computational cost through a combination of convolution and Fourier layers. We conduct experiments to evaluate the performance of DCNO on various datasets, including the multiscale elliptic equation, its inverse problem, Navier-Stokes equation, and Helmholtz equation. We show that DCNO strikes an optimal balance between accuracy and computational cost and offers a promising solution for multiscale operator learning.

Keywords: operator learning, multiscale PDE, dilated convolution, high frequency features, spectral bias

## 1. Introduction

In recent years, operator learning methods such as Fourier Neural Operator (FNO) (Li et al., 2020b), Galerkin Transformer (GT) (Cao, 2021), and Deep Operator Network (DeepONet) (Lu et al., 2021) have emerged as powerful tools for computing parameter-to-solution maps of partial differential equations (PDEs). In this paper, we focus on multiscale PDEs that encompass multiple temporal/spatial

scales. These multiscale PDE models are widely prevalent in physics, engineering, and other disciplines, playing a crucial role in addressing complex practical problems such as reservoir modeling, atmosphere and ocean circulation, and high-frequency scattering.

A well-known challenge with neural networks is their tendency to prioritize learning low-frequency components before high frequencies-a phenomenon referred to as 'spectral bias' or 'frequency principle' (Rahaman et al., 2019; Xu et al., 2019). This presents challenges when applying neural networks to functions characterized by multiscale or high-frequency properties, adapting neural network architectures (Cai and Xu, 2019; Wang et al., 2021) have been proposed to capture high-frequency details . In the context of operator learning, existing methods such as FNO and GT have shown spectral bias when applied to multiscale PDEs, as observed in Liu et al. (2023). To address this issue and recover high-frequency features, Liu et al. (2023) introduced an approach based on hierarchical attention and H 1 loss. However, despite providing improved accuracy, the method's high computational cost and implemational complexity to some extent counterbalances its strength.

In this paper, we present a novel method that strikes a balance among accuracy, computational cost, and the preservation of multiscale features. Our approach utilizes a carefully designed architecture that combines the strengths of dilated convolutions and Fourier layers. Dilated convolutions (Holschneider et al., 1990), also known as atrous convolutions, expand the kernel of a convolution layer in a convolutional neural network (CNN) by introducing gaps (holes) between the kernel elements. This technique allows for selectively skipping input values with specific step sizes, effectively covering a larger receptive field over the input feature map without introducing extra parameters or computational overhead. As a result, we can efficiently capture high-frequency local details. On the other hand, we leverage Fourier layers to capture the smooth global components of the data. DCNO achieves higher accuracy compared to existing models while maintaining lower computational costs by utilizing efficient implementations of both convolution and Fourier layers. This makes our approach well-suited for applications that require the preservation of multiscale features.

Figure 1.1: We demonstrate the effectiveness of the DCNO scheme using a challenging multiscale trigonometric benchmark. The coefficient and corresponding solution derivative are presented in (a) and (b), respectively (refer to Section 4.1.2) for a detailed description. We observe that DCNO accurately captures the solution derivatives. In (c) and (d), we analyze the (testing) dynamics for high-frequency ( &gt; 10 π ) and low-frequency ( ≤ 10 π ) errors, respectively. It is evident that DCNO achieves outstanding performance in terms of both high-frequency and low-frequency errors.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-6343804f-0342-4548-8ba7-7c052ae4f32f.jpg)

## 2. Background and Related Work

## 2.1. Multiscale PDEs

We briefly introduce some representative multiscale PDEs in this section. One notable example is the class of multiscale elliptic PDEs, which involve coefficients varying rapidly and are often encountered in heterogeneous and random media applications. For smooth coefficients, the coefficient to solution map can be effectively resolved by the FNO parameterization (Li et al., 2020b). However, when dealing with multiscale/rough coefficients, the presence of fast oscillation, high contrast ratios, and non-separable scales pose significant challenges from both scientific computing (Branets et al., 2009) and operator learning (Liu et al., 2023) perspectives. Other notable examples include the Navier-Stokes equation, which models fluid flow and exhibits turbulence behavior at high Reynolds numbers, and the Helmholtz equation, which models time-harmonic acoustic waves and is challenging to solve in the high wave number regime. In these multiscale PDEs, the accurate prediction of physical phenomena and properties necessitates the resolution of high-frequency components.

Numerical Methods for Multiscale PDEs. Multiscale PDEs, even with fixed parameters, present a challenge for classical numerical methods, as their computational cost typically scales inversely proportional to the finest scale ε of the problem. To overcome this issue, multiscale solvers have been developed by incorporating microscopic information to achieve computational cost independent of ε . One such technique

is numerical homogenization (Engquist and Souganidis, 2008), which identifies lowdimensional approximation spaces adapted to the corresponding multiscale operator. Similarly, fast solvers like multilevel/multigrid methods (Hackbusch, 2013; Xu and Zikatanov, 2017) and wavelet-based multiresolution methods (Brewster and Beylkin, 1995; Beylkin and Coult, 1998) may face limitations when applied to multiscale PDEs (Branets et al., 2009), while multilevel methods based on numerical homogenization techniques, such as Gamblets (Owhadi, 2017), have emerged as a way to discover scalable multilevel algorithms and operator-adapted wavelets for multiscale PDEs. In recent years, there has been increasing exploration of neural network methods for solving multiscale PDEs despite the spectral bias or frequency principle (Rahaman et al., 2019; Ronen et al., 2019; Xu et al., 2019) indicating that deep neural networks (DNNs) often struggle to effectively capture high-frequency components of functions. To address this limitation, specifically designed neural solvers (Li et al., 2020a; Wang et al., 2021; Li et al., 2023) have been developed to mitigate the spectral bias and accurately solve multiscale PDEs (with fixed parameters).

## 2.2. Neural Operator for Multiscale PDEs

Neural operators, as proposed by Li et al. (2020b); Gupta et al. (2021), have shown great promise in capturing the input-output relationship of parametric partial differential equations (PDEs). However, multiscale PDEs introduce new challenges for neural operators. Fourier or wavelet transforms, which are central to the construction of Li et al. (2020b); Gupta et al. (2021), may not always be effective, even for multiscale PDEs with fixed parameters. Moreover, while universal approximation theorems exist for FNO-type models (Kovachki et al., 2021), achieving a meaningful convergence rate often requires 'excessive smoothness' that may be absent for multiscale PDEs. Additionally, aliasing errors becomes significant in multiscale PDEs (Bartolucci et al., 2023), raising concerns about continuous-discrete equivalence. The work by Liu et al. (2023) addresses the issue of spectral bias in (multiscale) operator learning and highlights the challenges faced by existing neural operators in capturing high-frequency components of multiscale PDEs. These neural operators tend to prioritize the fitting of low-frequency components over high-frequency ones, limiting their ability to accurately capture fine details. To overcome this limitation, Liu et al. (2023) proposes a new architecture for multiscale operator learning that leverages hierarchical attention mechanisms and a tailored loss function. While these innovations help reduce the spectral bias and improve the prediction of multiscale solutions, it is worth noting that hierarchical attention induces a significant computational cost.

## 2.3. Dilated convolution

In this paper, we focus on utilizing dilated convolutions to capture the highresolution components of the data due to their simplicity and efficiency. Dilated convolution, also known as atrous convolution, was initially developed in the 'algorithme ` trous' for wavelet decomposition (Holschneider et al., 1990). a Its primary purpose was to increase image resolution and enable dense feature extraction without additional computational cost in deep convolutional neural networks (CNNs) by inserting 'holes' or zeros between pixels in convolutional kernels. By incorporating dilated convolution, networks can enlarge receptive fields, capture more global information, and gather contextual details, which are crucial for dense prediction tasks such as semantic segmentation. Various approaches have been proposed to leverage dilated convolution for this purpose (Yu and Koltun, 2015; Wang et al., 2018), and have demonstrated comparable results compared with U-Net and attention based models.

More recently, dilated convolution has also found applications in operator learning, such as the Dil-ResNet used for simulating turbulent flow (Stachenfeld et al., 2022). However, our work demonstrates that using dilated convolution alone is not sufficient to accurately capture the solution. Instead, we propose an interwoven global-local architecture of Fourier layers with dilated convolution layers. Furthermore, while Dil-ResNet requires up to 10 million training steps to achieve satisfactory results, our model offers a more efficient approach. It is worth noting that there are many alternative approaches to extract multiscale features, inspired by developments in numerical analysis and computer vision. These include hierarchical matrix methods (Fan et al., 2019), hierarchical attention (Liu et al., 2021, 2023), U-Net (Ronneberger et al., 2015) and U-NO (Rahman et al., 2022), wavelet-based methods (Gupta et al., 2021), among others. In Section 4, we will conduct a comprehensive benchmark of these different multiscale feature extraction techniques to evaluate their performance.

## 3. Methods

We adopt a data-driven approach to approximate the operator S : H ↦→ H 1 2 as in references Li et al. (2020b); Cao (2021); Liu et al. (2023). The operator S maps between two infinite-dimensional Banach spaces H 1 and H 2 , and aims to find the solution to the parametric partial differential equation (PDE) L a ( u ) = f , where the input/parameter a ∈ H 1 is drawn from a distribution µ , and the corresponding output/solution u ∈ H 2 .

To be specific, in this paper, our objective is to address the following operator learning problems:

↦

- · Approximating the nonlinear mapping S : a → u := S ( a ) from the varying parameter a to the solution u .

↦

- · Solving the inverse coefficient identification problem, which involves recovering the coefficient from a noisy measurement ˆ of the solution u u . In this scenario, we aim to approximate S -1 : ˆ u → a := S -1 ( ˆ). u

Figure 3.1: The architecture of the DCNO neural operator.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-f2b8df35-362c-4d98-8f98-d3287daf3a10.jpg)

## 3.1. Model architectures

In our model, we employ an Encode-Process-Decode architecture (Sanchez-Gonzalez et al., 2018, 2020), as shown in Figure 3.1.

- · The encoder incorporates a patch embedding function denoted as P, which utilizes a convolutional neural network (CNN). This step is performed to lift the input a x ( ) to a higher-dimensional channel (feature) space.
- · The processor part of the model comprises alternating Fourier layers (F layers) and Convolution layers (C layers). The role of the F layer is to approximate the low-frequency components, while the C layer is responsible for extracting high-frequency features. This alternating approach allows the DCNO model to effectively handle both low-frequency and high-frequency components present in the data. To gain further insights into the influence of the F and C layers, an ablation study is conducted, as described in Section 4.3.2. This study provides information about the impact of these layers on the model's performance and enhances our understanding of their significance within the overall architecture.

Figure 3.2: An example of a two-layer dilated convolution with dilation rates (1, 3).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-e3d08e37-af14-4375-b442-d6f24bac4fdd.jpg)

- · The decoder in our model adopts a structure similar to the one proposed in Cao (2021). Its primary role is to project high dimensional feature v x ( ) back to the target dimension of solutions. The decoder comprises two F layers, with the final F layer omitting an activation function, followed by a three-layer feedforward neural network ( FFN ).

F layers:. The Fourier layers, as proposed in Li et al. (2020b), consist of two main components that operate on the input v x ( ). In the first component, the input undergoes the Fourier transform, followed by a linear transform R acting only on the lower Fourier modes while filtering out the higher modes. The modified input is then obtained by applying the inverse Fourier transform. The first component of the Fourier layer aims to preserve low-frequency global information while reducing the influence of high-frequency components. The second component of the Fourier layers incorporates a local (pointwise) linear transform C to extract high-frequency information. Additionally, the outputs from both components are combined using the GELU activation function.

It is important to note that while the second component of the F layer helps capture some high-frequency details, relying solely on this part is not sufficient to accurately capture high-frequency information. This limitation is why the Fourier neural operator ( FNO ) approach may not perform well for multiscale PDE problems.

C layers:. Each convolution layer includes three convolutional neural networks, each utilizing a kernel of size 3 and followed by a GELU activation function. These convolutional neural networks employ dilation rates of (1 , 3 9 3 1). The dilation rate , , , determines the spacing between the points with which each point is convolved. A dilation rate of 1 corresponds to a regular convolution where each point is convolved with its immediate neighbors. Larger dilation rates, such as 3 and 9, expand the receptive field of each point to include more distant points. Figure 3.2 illustrates an example of a two-layer dilated convolution with dilation rates of (1 , 3). In the first layer, each point is convolved with its neighbors at a distance of 1, and in the second layer, each point is convolved with its neighbors at a distance of 3. As a result, the central red cell has a 9 × 9 receptive field. By incorporating multiple dilation rates, the model can capture long-range dependencies and maintain communication between distant points. This approach enhances the model's ability to capture both local and global information. Residual connections are applied to alleviate the vanishing gradient problem.

In summary, the combination of convolution layers with multiple dilation rates effectively enlarges the receptive fields of the network, facilitating the aggregation of global information and leveraging the advantages of convolutions in extracting

localized features (see Section 4.5). However, in most cases in operator learning, relying solely on convolutional layers may not yield satisfactory results, and the combination with Fourier layers can boost performance. See the ablation results in Section 4.3.2 for more details.

## 4. Experiments

In this section, we present a series of numerical experiments comparing DCNO with different operator models based on several metrics, including relative L 2 error, parameter count, memory consumption, and training time per epoch. The goal of these experiments is to assess the performance of various operator models through different types of tasks. Specifically, we investigate the performance of these models in the context of multiscale elliptic equations, time-dependent Navier-Stokes equations, inverse coefficient identification for multiscale elliptic equations and the Helmholtz equation. The results consistently demonstrate that the DCNO model outperforms other operator models. It achieves superior accuracy, robustness, and cost-accuarcy trade-off in all the considered scenarios.

## 4.1. Data generation

Multiscale elliptic equations are a fundamental class of problems, exemplified by the following second-order elliptic equation in divergence form:

$$\begin{cases} - \nabla \cdot ( a ( x ) \nabla u ( x ) ) = f ( x ), & x \in D \\ u ( x ) = 0, & x \in \partial D \end{cases}$$

Here, the coefficient a x ( ) satisfies 0 &lt; a min ≤ a x ( ) ≤ a max for all x ∈ D , and f ∈ H -1 ( D ; R ) represents the forcing term. The coefficient-to-solution map is denoted as S : L ∞ ( D ; R + ) → H 1 0 ( D ; R ), where u = S ( a ). The coefficient a x ( ) may exhibit rapid oscillations (e.g., a x ( ) = a x/ε ( ) with ε ≪ 1), high contrast ratios with a max /a min ≫ 1, and even a continuum of non-separable scales. Handling rough coefficients poses significant challenges from both scientific computing (Branets et al., 2009) and operator learning perspectives. In following, we give the details of the two examples of multiscale elliptic equations benchmarked in this paper.

## 4.1.1. Darcy rough example

Darcy flow, originally introduced by Darcy (Darcy, 1856), has diverse applications in modeling subsurface flow pressure, linearly elastic material deformation, and electric potential in conductive materials. The two-phase coefficients and solutions are generated using the approach outlined in https://github.com/zongyi-li/

fourier\_neural\_operator/tree/master/data\_generation . Given the computational domain [0 , 1] 2 , the coefficients a x ( ) are generated according to

$$a \sim \mu \coloneqq \psi _ { \# } \mathcal { N } \left ( 0, ( - \Delta + c I ) ^ { - 2 } \right ),$$

where ∆ represents the Laplacian with zero Neumann boundary condition. The mapping ψ : R → R takes the value 12 for the positive part of the real line and 2 for the negative part, with a contrast of 6. The push-forward is defined in a pointwise manner. These generated datasets serve as benchmark examples for operator learning in various studies, including Li et al. (2020b), Gupta et al. (2021), and Cao (2021). The parameter c can be used to control the 'roughness' of the coefficient and corresponding solution. In the aforementioned references, the parameter c is set as c = 9, while in Liu et al. (2023), a value of c = 20 is used to generate a rougher coefficient. The forcing term is fixed as f ( x ) ≡ 1. Solutions u are obtained using a second-order finite difference scheme on a 512 × 512 grid. Lower-resolution datasets are created by sub-sampling from the original dataset through linear interpolation.

## 4.1.2. Multiscale trigonometric example

Multiscale trigonometric coefficients are described in Owhadi (2017), as an example of highly oscillatory coefficients. Given the domain D [ -1 1] , the coefficient , 2 a x ( ) is specified as follows:

$$a ( x ) = \prod _ { k = 1 } ^ { 6 } \left ( 1 + \frac { 1 } { 2 } \cos ( a _ { k } \pi ( x _ { 1 } + x _ { 2 } ) ) \right ) \left ( 1 + \frac { 1 } { 2 } \sin ( a _ { k } \pi ( x _ { 2 } - 3 x _ { 1 } ) ) \right )$$

Here, a k is uniformly distributed in [2 k -1 , 1 5 . × 2 k -1 ]. The forcing term is fixed at f ( x ) ≡ 1. To obtain the reference solutions, the P 1 Finite Element Method (FEM) is employed on a 1023 × 1023 grid. Lower-resolution datasets are generated by downsampling the higher-resolution dataset through linear interpolation.

Figure 4.1: (a) multiscale trigonometric coefficient, (b)comparison of predicted solutions on the slice x = 0, (c) comparison of predicted derivative ∂u ∂y on the slice x = 0.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-7893829e-993a-4e56-8043-f37341db337e.jpg)

We present the multiscale trigonometric coefficient, reference solution, and a comparison with other operator learning models in Figure 4.1. Among the models considered, DCNO demonstrates superior accuracy in predicting function values and, more importantly, accurately captures the fine-scale oscillations. This is evident in the predicted derivatives shown in (c) of Figure 4.1.

## 4.2. Training and evaluation setup

Unless stated otherwise, the train-val-test split datasets used consist of 1000, 100, and 100 samples, respectively, with a maximum of 500 training epochs and batch size of 8. The Adam optimizer is utilized with a decay of 1e -4 and a 1cycle learning rate scheduler (Smith and Topin, 2019). For the Navier-Stokes equations, the dataset resolution is set to 64 × 64, again with a maximum of 500 epochs and batch size of 32. All experiments are executed on an NVIDIA A100 GPU.

## 4.3. Multiscale Elliptic Equations

We examine the effectiveness of the DCNO model on the multiscale elliptic equation, given by the following second-order linear elliptic equation,

$$\begin{cases} - \nabla \cdot ( a ( x ) \nabla u ( x ) ) = f ( x ), & x \in D \\ u ( x ) = 0, & x \in \partial D \end{cases}$$

with rough coefficients and Dirichlet boundary conditions. Our goal is to approximate the operator S : L ∞ ( D ; R + ) → H 1 0 ( D ; R ), which maps the coefficient function

a x ( ) to the corresponding solution u . We assess the model on two-phase Darcy rough coefficients (Darcy rough) given in Liu et al. (2023) where the coefficients are significantly rougher compared to the well-known benchmark proposed in Li et al. (2020b). We also consider multiscale trigonometric coefficients with higher contrast, following the setup in Owhadi (2017); Liu et al. (2023). The coefficients and solutions are displayed in Figure 4.3.2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-8ecc1074-5044-4d8c-8cc1-f19b4f6cb761.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-bc207bf6-b38c-4f1f-9cb6-e4f590adfabc.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-0d1e55d6-f853-4de9-b917-feba67bd643f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-f637b3ae-bbe0-4345-9031-5cd1ddeb9673.jpg)

(a)

coefficient

(b)

reference so-

lution

(c) abs. error of DCNO (in log 10 scale) (d) abs. error of FNO (in log 10 scale)

Figure 4.2: Top: Darcy rough example, (a) coefficient, (b) reference solution, (c) DCNO, absolute (abs.) error, (d) FNO, abs. error; Bottom: multiscale trignometric example, (a) coefficient (in log 10 scale), (b) reference solution, (c) DCNO, abs. error (in log 10 scale), (d) FNO, abs. error (in log 10 scale).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-31777d1e-918c-474a-a79d-5ad5f30979c9.jpg)

Figure 4.3: Costaccuracy (in log 2 scale) trade-off

The experimental results for multiscale elliptic equations conducted at different resolutions are presented in Table 1, and can be summarized below,

- · DCNO achieves the lowest relative error compared to other neural operators at various resolutions, and the errors remain approximately invariant with the resolution. Compared to FNO, DCNO has a remarkable accuracy improvement of 71% and 69% in the cases of Darcy rough and multiscale trigonometric. Furthermore, we observe that DCNO achieves the best cost-accuracy tradeoff among all the neural operators we tested in Figure 4.3, where the cost is measured by the training time per epoch for s = 256 (input resolution) of the Darcy rough example.
- · Dil-ResNet, with its convolutional architecture, has the fewest parameters among all models. However, it requires more memory and longer training time per epoch compared to DCNO and FNO. Furthermore, its accuracy is not as ideal as the other models.

These results highlight the superior performance of the DCNO model in terms of accuracy and efficiency, making it a favorable choice for solving multiscale elliptic equations.

| Model      | Parameters × 10 6   | Memory (GB)   | Time per epoch (s)   | Darcy rough       | Trigonometric     |
|------------|---------------------|---------------|----------------------|-------------------|-------------------|
| FNO        | 2 . 37              | 1 . 79        | 5 . 65               | 1 . 749 ± 0 . 029 | 1 . 744 ± 0 . 034 |
| MWT        | 9 . 81              | 2 . 54        | 17 . 94              | 1 . 289 ± 0 . 043 | 0 . 929 ± 0 . 029 |
| U-NO       | 16 . 39             | 1 . 45        | 8 . 71               | 1 . 159 ± 0 . 021 | 1 . 183 ± 0 . 018 |
| GT         | 2 . 22              | 9 . 32        | 35 . 25              | 1 . 351 ± 0 . 030 | 1 . 073 ± 0 . 013 |
| HANO       | 8 . 65              | 2 . 33        | 9 . 92               | 0 . 556 ± 0 . 031 | 0 . 594 ± 0 . 148 |
| Dil-ResNet | 0 . 58              | 5 . 71        | 10 . 69              | 4 . 794 ± 0 . 175 | 1 . 233 ± 0 . 050 |
| DCNO       | 2 . 51              | 3 . 94        | 10 . 14              | 0 . 531 ± 0 . 026 | 0 . 541 ± 0 . 051 |

Table 1: Benchmarks on multiscale elliptic equations at various input resolution s . Performance are measured with relative L 2 errors ( × 10 -2 ), number of parameters, memory consumption for a batch size of 8, and time per epoch for s = 256 of the Darcy rough example during the training process.

4.3.1. Hyperparameter Study(dilated convolution)

| Dilation                |   Parameters for C layers | Darcy ( × 10 - 2 )   | Trigonometric ( × 10 -   |
|-------------------------|---------------------------|----------------------|--------------------------|
| (1 , 3 , 9 , 3 , 1)     |                    138720 | 0 . 531              | 0 . 541                  |
| (1 , 1 , 1 , 1 , 1)     |                    138720 | 0 . 792              | 0 . 650                  |
| (1 , 3 , 9)             |                     83232 | 0 . 623              | 0 . 814                  |
| (1 , 3 , 1)             |                     83232 | 0 . 691              | 0 . 641                  |
| (1 , 1 , 1)             |                     83232 | 0 . 768              | 0 . 805                  |
| (1)                     |                     27744 | 1 . 065              | 0 . 976                  |
| (1) (kernel size 9 × 9) |                    248928 | 0 . 778              | 0 . 627                  |

Table 2: Hyperparameter study: impact of dilation rate and kernel size

We conducted a hyperparameter study to investigate the influence of different dilation rates and kernel sizes in the C layers of the DCNO model for multiscale elliptic PDEs. Those parameters will in turn determine the reception field and parameter count of the convolution. The results of this study are summarized in Table 2. The default size of the convolutional kernel is 3 × 3, unless otherwise specified. The dilation rates determine the configuration of the C layers in the DCNO model, where dilation rates of (1 , 3 9) means that the C layers consist of three dilated convolutions , with dilation factors of 1, 3, and 9. As expected, increasing the number of layers in the C layers leads to improved results. This can be attributed to the fact that additional layers can capture more complex features, thereby enhancing the model's

accuracy. To assess the impact of hierarchical dilated convolution, we compare the outcomes obtained with dilation rates (1 , 1 1 1 1) and (1 1 1) against those acquired , , , , , with dilation rates (1 , 3 9 3 1) and (1 3 9). , , , , , The results clearly demonstrate that hierarchical dilated convolution has a positive effect on the outcomes. This suggests that the ability to capture multiscale information through multiple dilation rates proves beneficial in enhancing the performance of the model. Additionally, we compared the performance of a single layer with a large 9 × 9 convolutional kernel against a layered dilated convolution with a 3 × 3 kernel size. Interestingly, we found that the latter yielded better results with fewer parameters. This highlights the effectiveness of layered dilated convolutions in capturing complex features while maintaining model efficiency.

## 4.3.2. Hyperparameter Study(architecture)

Table 3: Ablation Study of DCNO: different combinations of F- and C- layers.

| Model            | Darcy rough   | Trigonometric   |
|------------------|---------------|-----------------|
| DCNO(FCFCFCF)    | 0 . 531       | 0 . 541         |
| DCNO 1 (FFFFFFF) | 1 . 344       | 1 . 256         |
| DCNO 2 (CCCCCCC) | 4 . 324       | 1 . 045         |
| DCNO 3 (CCCCCCF) | 0 . 594       | 0 . 842         |
| DCNO 4 (FCCCCCC) | 1 . 842       | 0 . 923         |
| DCNO 5 (CFFFFFF) | 0 . 707       | 0 . 732         |
| DCNO 7 (FFFFFFC) | 1 . 180       | 0 . 864         |
| DCNO 8 (FFFFCCC) | 1 . 240       | 0 . 631         |
| DCNO 9 (CCCFFFF) | 0 . 676       | 0 . 872         |

To assess the impact of different components of the DCNO model on multiscale elliptic PDEs, we conducted experiments testing various combinations of F- and C-layers.

DCNO 1 is a variant of the DCNO model where all C layers are replaced with F layers, akin to an enhanced version of the FNO model. However, we observed that DCNO does not achieve comparable accuracy to the full DCNO model. In contrast, 1 DCNO 2 replaces all F layers with C layers within the DCNO framework. Interestingly, we found that DCNO 2 demonstrates significantly better performance in the multiscale trigonometric case compared to the Darcy rough case. This observation suggests that dilated convolutions in the C layers effectively enhance results when dealing with rough coefficients. The C layers' capability to capture local interactions and finer details appears particularly advantageous in scenarios with pronounced coefficient variations.

Figure 4.4: Top: Darcy rough example, Bottom: multiscale trignometric example.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-a9ef4b42-6500-40fb-a0d7-2d7b849f0c35.jpg)

While DCNO , DCNO , DCNO , and DCNO 3 4 5 6 each incorporate only one F layer or C layer compared to DCNO 1 and DCNO , their mixed utilization of both F and 2 C layers enables them to outperform DCNO 1 with only F layers and DCNO 2 with only C layers. This underscores the critical importance of combining F and C layers to achieve superior performance. F layers are adept at capturing global information and long-range dependencies, while C layers excel in capturing local interactions and finer details. Both components are essential for accurately modeling the complex dynamics of the system.

To underscore the importance of interleaved architecture, we scrutinized DCNO 7 and DCNO . DCNO 8 7 features three C layers following four F layers, while DCNO 8 places three C layers before four F layers. However, both DCNO 7 and DCNO 8 proved less effective than the complete DCNO model. This suggests that the specific interleaving of C and F layers within the DCNO model's architecture is critical for achieving optimal performance.

In summary, the comparison among the DCNO variants underscores the critical roles played by both F layers and C layers within the DCNO model. These components are essential and complement each other in achieving the superior performance demonstrated by the complete DCNO model. F layers excel in capturing global information and long-range dependencies, whereas C layers are adept at capturing local interactions and finer details. The interleaved architecture of the DCNO model, with a carefully ordered combination of F and C layers, is pivotal for achieving its outstanding performance.

## 4.3.3. Spectral bias

Figure 4.5: Error dynamics in the frequency domain for multiscale trigonometric example.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-44f194e1-2fce-4f19-bf60-3f1446a1bd4b.jpg)

The spectral bias, also known as the frequency principle, suggests that deep neural networks (DNNs) face challenges in effectively learning high-frequency components of functions that exhibit variations at multiple scales. This phenomenon has been extensively studied and discussed in the literature (Rahaman et al., 2019; Ronen et al., 2019; Xu et al., 2019) in the context of function approximation.

In the context of operator learning, we conducted an analysis on the relative error spectrum dynamics of a multiscale trigonometric example presented in Section 4.3. The analysis is illustrated in Figure 4.5. To begin, we computed the Fourier transform of the relative error in the frequency domain [ -128 π, 128 π ] 2 . Our examination focused on the error density ρ r ( ) within the annulus A r ( ), which satisfies the equation ∫ A r ( ) ρ r rdr ( ) = ∑ i ∈ A r ( ) ϵ i . The annulus A r ( ) := B r ( +1) \ B r ( ), where B r ( ) represents a sphere of radius r . The term ϵ i represents the Fourier-transformed relative error at a lattice point indexed by i in the discrete frequency domain.

In Figure 4.5, the x-axis represents dominant frequencies arranged from low to high frequency, while the y-axis represents the number of training epochs. The plot reveals several significant findings. Firstly, the DCNO model exhibits a faster error decay for higher frequencies, demonstrating its effectiveness in capturing highfrequency components. Additionally, the DCNO model shows a consistent reduction

in errors across all frequencies, indicating its capability to learn variations at multiple scales proficiently. Moreover, the DCNO model outperforms other methods by achieving lower testing errors, highlighting its effectiveness in addressing spectral bias and accurately predicting functions with diverse scale variations. This analysis provides compelling evidence of the DCNO model's advantages in operator learning and its ability to handle spectral bias while accurately modeling functions across various frequency ranges.

## 4.4. Navier-Stokes Equation

In this section, we focus on the following 2D Navier-Stokes equation in vorticity form on the unit torus T , as benchmarked in (Li et al., 2020b).

$$\begin{cases} \partial _ { t } w ( x, t ) + u ( x, t ) \cdot \nabla w ( x, t ) = \nu \Delta w ( x, t ) + f ( x ), & x \in ( 0, 1 ) ^ { 2 }, t \in ( 0, T ] \\ \nabla \cdot u ( x, t ) = 0, & x \in ( 0, 1 ) ^ { 2 }, t \in [ 0, T ] \\ w ( x, 0 ) = w _ { 0 } ( x ), & x \in ( 0, 1 ) ^ { 2 } \end{cases}$$

In the context of fluid dynamics, the variables used in the equations have the following interpretations:

- · The velocity field is represented by the symbol u .
- · The vorticity field is denoted as w , and it is defined as the curl of the velocity field, i.e., w = ∇× u .
- · The initial vorticity distribution is denoted by w 0 .
- · The viscosity coefficient is represented by ν .
- · The forcing term is given by f ( x ) = 0 1 (sin (2 . π x ( 1 + x 2 )) + cos (2 π x ( 1 + x 2 ))).

The Reynolds number, denoted as Re, is a dimensionless parameter defined as Re := ρuL ν , where ρ is the density (assumed to be 1), u is the fluid velocity, and L is the characteristic length scale of the fluid (set to 1). The Reynolds number is inversely proportional to the viscosity coefficient ν . An increase in the Reynolds number tends to promote the transition of the flow to turbulence. The initial vorticity distribution w 0 ( x ) is generated from a probability measure µ . Specifically, w 0 ∼ µ , where µ = N ( 0 7 , 3 2 / ( -∆+49 ) I -2 5 . ) , and periodic boundary conditions are applied.

The vorticity variable is denoted as ω x, t ( ), where x ∈ T represents the spatial domain and t ∈ [0 , T ] represents the time interval. Our goal is to learn the mapping from the vorticity at the previous time step, w ,t ( · -1), to the vorticity at the current time step, w ,t ( · ).

During the training process, we follow the approach used in Liu et al. (2023). The predicted vorticity at each time step, denoted as ˜ , w t is obtained using a recurrence relation: ˜ w t = G ( ˜ w t -1 ). Here, G represents an operator that approximates the mapping from the previous vorticity to the current vorticity. In the training process, we set ˜ w i = w i , where ˜ w i is the predicted vorticity and w i is the true vorticity at time step i .

During testing, we employ a 'rollout' strategy, commonly utilized in fluid dynamics studies (Li et al., 2020b; Rahman et al., 2022). This strategy entails predicting the vorticity at each time step using the previously mentioned recurrence relation, with only the initial vorticity ˜ w 0 = w 0 known. Our objective is to evaluate the DCNO model and other benchmark operators on their capability to accurately forecast the intricate dynamics of fluid flow over time by assessing their performance in predicting vorticity dynamics.

Table 4 presents a comprehensive summary of the Navier-Stokes experiment results. For lower Reynolds numbers, FNO, MWT, U-NO, HANO and DCNO achieve similar levels of accuracy. However, as the Reynolds number increases, the DCNO models demonstrate a notable advantage in accuracy. This highlights the effectiveness of the DCNO models in capturing the complex dynamics of fluid flow and accurately predicting vorticity, especially in scenarios with higher Reynolds numbers. Overall, the results demonstrate that DCNO consistently provides more accurate predictions compared to other methods while maintaining a reasonable computational cost.

Table 4: Benchmarks on Navier Stokes equation.Performance are measured with relative L 2 errors , number of parameters, memory consumption for a batch size of 16, and time per epoch for ( ν = 1e -5 , T 0 = 10 s, T = 20 ) during the training process. s

| Model      | Parameters × 10 6   | Memory Requirement (GB)   | Time per epoch (s)   | ν = 1e - 3 T 0 = 1 s T = 19 s   | ν = 1e - 4 T 0 = 1 s T = 19 s   | ν = 1e - 5 T 0 = 1 s T = 19 s   |
|------------|---------------------|---------------------------|----------------------|---------------------------------|---------------------------------|---------------------------------|
| FNO        | 2 . 37              | 0 . 22                    | 4 . 20               | 4 . 489 × 10 - 4                | 1 . 004 × 10 - 2                | 6 . 288 × 10 - 2                |
| FNO †      | 25 . 97             | 0 . 99                    | 10 . 26              | 6 . 095 × 10 - 4                | 5 . 100 × 10 - 3                | 3 . 555 × 10 - 2                |
| MWT        | 9 . 81              | 0 . 32                    | 27 . 42              | 7 . 075 × 10 - 4                | 1 . 015 × 10 - 2                | 5 . 288 × 10 - 2                |
| U-NO       | 11 . 91             | 0 . 22                    | 36 . 97              | 4 . 532 × 10 - 4                | 8 . 192 × 10 - 3                | 6 . 743 × 10 - 2                |
| HANO       | 15 . 37             | 0 . 82                    | 18 . 68              | 5 . 734 × 10 - 4                | 3 . 759 × 10 - 3                | 2 . 697 × 10 - 2                |
| Dil-ResNet | 0 . 58              | 0 . 71                    | 16 . 42              | 2 . 786 × 10 - 2                | 1 . 259 × 10 - 1                | 2 . 319 × 10 - 1                |
| DCNO       | 14 . 40             | 1 . 47                    | 16 . 21              | 4 . 056 × 10 - 4                | 3 . 764 × 10 - 3                | 2 . 227 × 10 - 2                |

Figure 4.6: Comparison of Navier Stokes solutions ( ν = 1 e -5) and DCNO predictions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-bcff6bb2-9415-4081-b0af-d7d9ea3fe102.jpg)

## 4.5. Inverse coefficient identification for multiscale elliptic PDEs

↦

In this section, we address an inverse coefficient identification problem using the same data as the previous example of multiscale elliptic PDE in section 4.3. Inverse problems play a crucial role in various scientific fields, including geological sciences and medical imaging. In this example, our objective is to learn an approximation to an ill-posed operator S -1 : H 1 0 ( D ) → L ∞ ( D ), where ˆ = u u + ϵN ( u ) → a . Here, ϵ represents the level of Gaussian noise added to both the training and evaluation data. The term N ( u ) accounts for the data-related noise. This task is challenging due to the ill-posed nature of the problem and the presence of noise.

↦

In Figure 4.7, we display the solution and predicted coefficients for the inverse coefficient identification problem at various levels of noise. Notably, even with a noise level of 10%, the predicted coefficient successfully recovers the interface present in the ground truth. This resilience to noise highlights the robustness and effectiveness of the coefficient prediction in capturing the underlying structure accurately.

The results of the inverse coefficient identification problem with noise are presented in Table 5. It is worth noting that the memory consumption and training time per epoch remain the same as reported in Table 1. The DCNO model outperforms other methods in this example, which highlights DCNO's ability to effectively address the challenges posed by this ill-posed inverse problem with noisy data. Interestingly, FNO and U-NO, known for their effectiveness in smoothing and filtering

high-frequency modes, encounter difficulties in recovering targets that exhibit highfrequency characteristics, such as irregular interfaces, highly oscillatory coefficients, and the presence of Gaussian noise. In contrast, Dil-ResNet (and MWT to a lesser extent) performs significantly better in this specific problem, particularly when the noise level is higher. The use of dilated convolutions in Dil-ResNet proves advantageous in capturing high-frequency features. While HANO achieves the second-best overall accuracy, it comes at a considerable computational cost.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-51217e89-3ce5-4f4b-8667-bcaba277927e.jpg)

Figure 4.7: These figures demonstrate the DCNO's prediction of coefficients given solutions with varying noise levels, ϵ : (a) Solution, slice at x=0.5. (b) Solution in 2D for ϵ = 0. (c) Solution for ϵ = 0 01. (d) Solution for . ϵ = 0 1. (e) Coefficient. . (f) Predicted coefficients for ϵ = 0. (g) Predicted coefficients for ϵ = 0 01. (h) Predicted coefficients for . ϵ = 0 1. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-671d3172-5418-4ef1-9b1c-353ad73d8991.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-1506bc58-e098-4260-91b0-c41ff299a559.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-789073fe-7ce3-4b85-b490-02bd1dba4cd3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-2fbd1b91-6044-4e5c-bb8e-7b1627e4a0cd.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-5faaef8f-d243-42f1-af13-9dd6c1de5b92.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-771f314d-db9b-4488-999b-02702b9c561d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-1f0e101c-1e33-4f49-8ed5-1563d0c77ef6.jpg)

Table 5: Relative error of the inverse coefficient identification. Also see Figure 4.7 for solutions and predicted coefficients at various noise levels.

|            | Darcy rough   | Darcy rough   | Darcy rough   | Darcy rough   | Darcy rough   | Darcy rough   |
|------------|---------------|---------------|---------------|---------------|---------------|---------------|
| Model      | ϵ =0          | ϵ =0.01       | ϵ =0.03       | ϵ =0.05       | ϵ =0.07       | ϵ =0.1        |
| FNO        | 0 . 327       | 0 . 328       | 0 . 330       | 0 . 334       | 0 . 336       | 0 . 339       |
| MWT        | 0 . 103       | 0 . 129       | 0 . 162       | 0 . 183       | 0 . 198       | 0 . 214       |
| U-NO       | 0 . 235       | 0 . 236       | 0 . 239       | 0 . 245       | 0 . 251       | 0 . 261       |
| GT         | 0 . 146       | 0 . 172       | 0 . 178       | 0 . 199       | 0 . 214       | 0 . 230       |
| HANO       | 0 . 095       | 0 . 118       | 0 . 149       | 0 . 167       | 0 . 179       | 0 . 197       |
| Dil-RESNET | 0 . 086       | 0 . 122       | 0 . 150       | 0 . 168       | 0 . 185       | 0 . 201       |
| DCNO       | 0 . 061       | 0 . 115       | 0 . 147       | 0 . 166       | 0 . 180       | 0 . 196       |

## 4.6. Helmholtz equations

We evaluate the performance of DCNO for the acoustic Helmholtz equation in highly heterogeneous media, as an illustrative example of multiscale wave phenomena. Solving this equation is particularly challenging and computationally expensive, especially for complex material and large geological models. We adopt the setup described in Freese et al. (2021),

$$\begin{cases} - \text{div} ( a ( x ) \nabla u ( x ) ) - \kappa ^ { 2 } u = f ( x ), & x \in D, \\ u ( x ) = 0, & x \in \partial D. \end{cases}$$

where the coefficient a x ( ) takes the value 1 or ε as shown in Figure 4.6 with ε -1 ∈ rand(128 256) , , κ = 9, and

$$f ( x _ { 1 }, x _ { 2 } ) = \begin{cases} 1 0 ^ { 4 } \exp \left ( \frac { - 1 } { 1 - \frac { ( x _ { 1 } - 0. 1 2 5 ) ^ { 2 } + ( x _ { 2 } - 0. 5 ) ^ { 2 } } { 0. 0 5 ^ { 2 } } } \right ), & ( x _ { 1 } - 0. 1 2 5 ) ^ { 2 } + ( x _ { 2 } - 0. 5 ) ^ { 2 } < 0. 0 5 ^ { 2 }, \\ 0, & \text{else.} \end{cases}$$

|       | Parameters 6   | Memory   | Time per epoch (s)   | Relative L 2 error   | Relative L 2 error   | Relative L 2 error   |
|-------|----------------|----------|----------------------|----------------------|----------------------|----------------------|
| Model | × 10           | (GB)     |                      | s=128                | s=256                | s=512                |
| FNO   | 2 . 37         | 1 . 79   | 5 . 61               | 5 . 033              | 5 . 405              | 6 . 295              |
| MWT   | 9 . 81         | 2 . 54   | 17 . 74              | 2 . 653              | 2 . 731              | 2 . 875              |
| U-NO  | 16 . 39        | 1 . 57   | 9 . 78               | 2 . 917              | 2 . 852              | 3 . 028              |
| GT    | 2 . 22         | 9 . 32   | 35 . 36              | 13 . 505             | 13 . 828             | 16 . 601             |
| HANO  | 8 . 65         | 2 . 33   | 9 . 96               | 2 . 472              | 3 . 286              | 5 . 189              |
| DCNO  | 2 . 51         | 3 . 94   | 10 . 14              | 2 . 571              | 2 . 682              | 2 . 748              |

Table 6: Benchmarks on Helmholtz equations at various input resolution s . Performance are measured with relative L 2 errors ( × 10 -2 ), number of parameters, memory consumption for a batch size of 8, and time per epoch for s = 256 during the training process.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00775v1-figure-d5dd9276-dfb0-4d17-84eb-8b72340a5882.jpg)

↦

Figure 4.8: The mapping a x ( ) → u . (a) Heterogeneous coefficient a x ( ), (b) the reference solution for ε -1 = 237 3, which is solved by . P 1 FEM implemented in FreeFEM++ (Hecht, 2012), (c)DCNO predicted solution

The Helmholtz equation presents a significant challenge due to the highly oscillatory nature of its solution, illustrated in Figure 4.6, and the presence of the κ -dependent pollution effect. However, the DCNO model has emerged as a notable solution, consistently outperforming other methods in effectively tackling this problem.

## 5. Conclusion

In this paper, we introduce DCNO (Dilated Convolution Neural Operator) as a novel and effective method for learning operators in multiscale PDEs. DCNO leverages the strengths of both Fourier layers, adept at representing low-frequency global components, and convolution layers with multiple dilation rates, capable of

capturing high-resolution local details. This hybrid architecture enables DCNO to surpass existing operator methods, providing a highly accurate and computationally efficient approach for learning operators in multiscale settings. Through extensive experiments, we demonstrate the effectiveness of DCNO in addressing multiscale PDEs, highlighting its superior performance and potential for diverse applications.

## References

Bartolucci, F., de B´ ezenac, E., Raoni´ c, B., Molinaro, R., Mishra, S., Alaifari, R., 2023. Are neural operators really neural operators? frame theory meets operator learning. arXiv preprint arXiv:2305.19913 .

Beylkin, G., Coult, N., 1998. A multiresolution strategy for reduction of elliptic pdes and eigenvalue problems. Applied and Computational Harmonic Analysis 5, 129-155.

Branets, L.V., Ghai, S.S., Lyons, S.L., Wu, X.H., 2009. Challenges and technologies in reservoir modeling. Communications in Computational Physics 6, 1.

Brewster, M.E., Beylkin, G., 1995. A multiresolution strategy for numerical homogenization. Applied and Computational Harmonic Analysis 2, 327-349.

Cai, W., Xu, Z.Q.J., 2019. Multi-scale deep neural networks for solving high dimensional pdes. arXiv preprint arXiv:1910.11710 .

Cao, S., 2021. Choose a transformer: Fourier or galerkin. Advances in neural information processing systems 34, 24924-24940.

Darcy, H., 1856. Les fontaines publiques de la ville de Dijon: exposition et application des principes ` suivre et des formules ` employer dans les questions de distribution a a d'eau. volume 1. Victor dalmont.

Engquist, B., Souganidis, P.E., 2008. Asymptotic and numerical homogenization. Acta Numerica 17, 147-190.

Fan, Y., Lin, L., Ying, L., Zepeda-N´ unez, L., 2019. A multiscale neural network based on hierarchical matrices. Multiscale Modeling &amp; Simulation 17, 1189-1213.

Freese, P., Hauck, M., Peterseim, D., 2021. Super-localized orthogonal decomposition for high-frequency helmholtz problems. arXiv preprint arXiv:2112.11368 .

| Gupta, G., Xiao, X., Bogdan, P., 2021. Multiwavelet-based operator learning for differential equations. Advances in Neural Information Processing Systems 34, 24048-24062.                                                                                                                                                         |
|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Hackbusch, W., 2013. Multi-grid methods and applications. volume 4. Springer Science & Business Media.                                                                                                                                                                                                                             |
| Hecht, F., 2012. FreeFem++. https://freefem.org/ . Version 4.13.                                                                                                                                                                                                                                                                   |
| Holschneider, M., Kronland-Martinet, R., Morlet, J., Tchamitchian, P., 1990. A real-time algorithm for signal analysis with the help of the wavelet transform, in: Wavelets: Time-Frequency Methods and Phase Space Proceedings of the In- ternational Conference, Marseille, France, December 14-18, 1987, Springer. pp. 286-297. |
| Kovachki, N., Lanthaler, S., Mishra, S., 2021. On universal approximation and error bounds for fourier neural operators. The Journal of Machine Learning Research 22, 13237-13312.                                                                                                                                                 |
| Li, X.A., Xu, Z.Q.J., Zhang, L., 2020a. A multi-scale dnn algorithm for nonlinear elliptic equations with multiple scales. arXiv preprint arXiv:2009.14597 .                                                                                                                                                                       |
| Li, X.A., Xu, Z.Q.J., Zhang, L., 2023. Subspace decomposition based dnn algorithm for elliptic type multi-scale pdes. Journal of Computational Physics 488, 112242.                                                                                                                                                                |
| Li, Z., Kovachki, N., Azizzadenesheli, K., Liu, B., Bhattacharya, K., Stuart, A., Anandkumar, A., 2020b. Fourier neural operator for parametric partial differential equations. arXiv preprint arXiv:2010.08895 .                                                                                                                  |
| Liu, X., Xu, B., Zhang, L., 2023. Mitigating spectral bias for the multiscale operator learning with hierarchical attention. arXiv preprint arXiv:2210.10890 .                                                                                                                                                                     |
| Liu, Z., Lin, Y., Cao, Y., Hu, H., Wei, Y., Zhang, Z., Lin, S., Guo, B., 2021. Swin transformer: Hierarchical vision transformer using shifted windows, in: Proceed- ings of the IEEE/CVF international conference on computer vision, pp. 10012- 10022.                                                                           |
| Lu, L., Jin, P., Pang, G., Zhang, Z., Karniadakis, G.E., 2021. Learning nonlinear operators via deeponet based on the universal approximation theorem of operators. Nature machine intelligence 3, 218-229.                                                                                                                        |

| Owhadi, H., 2017. Multigrid with rough coefficients and multiresolution operator decomposition from hierarchical information games. Siam Review 59, 99-149.                                                                                                                                                       |
|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Rahaman, N., Baratin, A., Arpit, D., Draxler, F., Lin, M., Hamprecht, F., Bengio, Y., Courville, A., 2019. On the spectral bias of neural networks, in: International Conference on Machine Learning, PMLR. pp. 5301-5310.                                                                                        |
| Rahman, M.A., Ross, Z.E., Azizzadenesheli, K., 2022. U-no: U-shaped neural oper- ators. arXiv preprint arXiv:2204.11127 .                                                                                                                                                                                         |
| Ronen, B., Jacobs, D., Kasten, Y., Kritchman, S., 2019. The convergence rate of neural networks for learned functions of different frequencies. Advances in Neural Information Processing Systems 32.                                                                                                             |
| Ronneberger, O., Fischer, P., Brox, T., 2015. U-net: Convolutional networks for biomedical image segmentation, in: Medical Image Computing and Computer- Assisted Intervention-MICCAI 2015: 18th International Conference, Munich, Ger- many, October 5-9, 2015, Proceedings, Part III 18, Springer. pp. 234-241. |
| Sanchez-Gonzalez, A., Godwin, J., Pfaff, T., Ying, R., Leskovec, J., Battaglia, P., 2020. Learning to simulate complex physics with graph networks, in: International conference on machine learning, PMLR. pp. 8459-8468.                                                                                        |
| Sanchez-Gonzalez, A., Heess, N., Springenberg, J.T., Merel, J., Riedmiller, M., Had- sell, R., Battaglia, P., 2018. Graph networks as learnable physics engines for inference and control, in: International Conference on Machine Learning, PMLR. pp. 4470-4479.                                                 |
| Smith, L.N., Topin, N., 2019. Super-convergence: Very fast training of neural net- works using large learning rates, in: Artificial intelligence and machine learning for multi-domain operations applications, SPIE. pp. 369-386.                                                                                |
| Stachenfeld, K., Fielding, D.B., Kochkov, D., Cranmer, M., Pfaff, T., Godwin, J., Cui, C., Ho, S., Battaglia, P., Sanchez-Gonzalez, A., 2022. Learned simulators for turbulence, in: International Conference on Learning Representations.                                                                        |
| Wang, P., Chen, P., Yuan, Y., Liu, D., Huang, Z., Hou, X., Cottrell, G., 2018. Under- standing convolution for semantic segmentation, in: 2018 IEEE winter conference on applications of computer vision (WACV), Ieee. pp. 1451-1460.                                                                             |

Wang, S., Wang, H., Perdikaris, P., 2021. On the eigenvector bias of fourier feature networks: From regression to solving multi-scale pdes with physics-informed neural networks. Computer Methods in Applied Mechanics and Engineering 384, 113938.

Xu, J., Zikatanov, L., 2017. Algebraic multigrid methods. Acta Numerica 26, 591721.

Xu, Z.Q.J., Zhang, Y., Luo, T., Xiao, Y., Ma, Z., 2019. Frequency principle: Fourier analysis sheds light on deep neural networks. arXiv preprint arXiv:1901.06523 .

Yu, F., Koltun, V., 2015. Multi-scale context aggregation by dilated convolutions. arXiv preprint arXiv:1511.07122 .
Bo Xu Xinliang Liu Lei Zhang
Dilated convolution neural operator for multiscale partial differential equations
This paper introduces a data-driven operator learning method for multiscale
partial differential equations, with a particular emphasis on preserving
high-frequency information. Drawing inspiration from the representation of
multiscale parameterized solutions as a combination of low-rank global bases
(such as low-frequency Fourier modes) and localized bases over coarse patches
(analogous to dilated convolution), we propose the Dilated Convolutional Neural
Operator (DCNO). The DCNO architecture effectively captures both high-frequency
and low-frequency features while maintaining a low computational cost through a
combination of convolution and Fourier layers. We conduct experiments to
evaluate the performance of DCNO on various datasets, including the multiscale
elliptic equation, its inverse problem, Navier-Stokes equation, and Helmholtz
equation. We show that DCNO strikes an optimal balance between accuracy and
computational cost and offers a promising solution for multiscale operator
learning.