A new numerical approach, called the “subdomain Chebyshev spectral method” is presented for calculation of the spatial derivatives in a curved coordinate system, which may be employed for numerical solutions of part...A new numerical approach, called the “subdomain Chebyshev spectral method” is presented for calculation of the spatial derivatives in a curved coordinate system, which may be employed for numerical solutions of partial differential equations defined in a 2D or 3D geological model. The new approach refers to a “strong version” against the “weak version” of the subspace spectral method based on the variational principle or Galerkin’s weighting scheme. We incorporate local nonlinear transformations and global spline interpolations in a curved coordinate system and make the discrete grid exactly matches geometry of the model so that it is achieved to convert the global domain into subdomains and apply Chebyshev points to locally sampling physical quantities and globally computing the spatial derivatives. This new approach not only remains exponential convergence of the standard spectral method in subdomains, but also yields a sparse assembled matrix when applied for the global domain simulations. We conducted 2D and 3D synthetic experiments and compared accuracies of the numerical differentiations with traditional finite difference approaches. The results show that as the points of differentiation vector are larger than five, the subdomain Chebyshev spectral method significantly improve the accuracies of the finite difference approaches.展开更多
This paper derives two new integrated and explicit boundary conditions, named the “explicit?normal version” and “explicit tangential versions” respectively for electromagnetic fields at an arbitrary interface betw...This paper derives two new integrated and explicit boundary conditions, named the “explicit?normal version” and “explicit tangential versions” respectively for electromagnetic fields at an arbitrary interface between two anisotropic media. The new versions combine two implicit boundary equations into a single explicit matrix formula and reveal the boundary values linked by a 3 × 3 matrix, which depends on the interface topography and model property tensors. We analytically demonstrate the new versions equivalent to the common implicit boundary conditions and their application to transformation of the boundary values in the boundary integral equations. We also give two synthetic examples that show recovery of the boundary values on a hill and a ridge, and highlight the advantage of the new versions of being a simpler and more straightforward method to compute the electromagnetic boundary values.展开更多
文摘A new numerical approach, called the “subdomain Chebyshev spectral method” is presented for calculation of the spatial derivatives in a curved coordinate system, which may be employed for numerical solutions of partial differential equations defined in a 2D or 3D geological model. The new approach refers to a “strong version” against the “weak version” of the subspace spectral method based on the variational principle or Galerkin’s weighting scheme. We incorporate local nonlinear transformations and global spline interpolations in a curved coordinate system and make the discrete grid exactly matches geometry of the model so that it is achieved to convert the global domain into subdomains and apply Chebyshev points to locally sampling physical quantities and globally computing the spatial derivatives. This new approach not only remains exponential convergence of the standard spectral method in subdomains, but also yields a sparse assembled matrix when applied for the global domain simulations. We conducted 2D and 3D synthetic experiments and compared accuracies of the numerical differentiations with traditional finite difference approaches. The results show that as the points of differentiation vector are larger than five, the subdomain Chebyshev spectral method significantly improve the accuracies of the finite difference approaches.
文摘This paper derives two new integrated and explicit boundary conditions, named the “explicit?normal version” and “explicit tangential versions” respectively for electromagnetic fields at an arbitrary interface between two anisotropic media. The new versions combine two implicit boundary equations into a single explicit matrix formula and reveal the boundary values linked by a 3 × 3 matrix, which depends on the interface topography and model property tensors. We analytically demonstrate the new versions equivalent to the common implicit boundary conditions and their application to transformation of the boundary values in the boundary integral equations. We also give two synthetic examples that show recovery of the boundary values on a hill and a ridge, and highlight the advantage of the new versions of being a simpler and more straightforward method to compute the electromagnetic boundary values.
