A high-order approximation method for semilinear parabolic equations on spheres

We describe a novel discretisation method for numerically solving (systems of) semilinear parabolic equations on Euclidean spheres. The new approximation method is based upon a discretisation in space using spherical basis functions and can be of arbitrary order. This, together with the fact that the solutions of semilinear parabolic problems are known to be infinitely smooth, at least locally in time, allows us to prove stability and convergence of the discretisation in a straight-forward way

Solving the Helmholtz Equation for the Neumann Boundary Condition for the Pseudosphere by the Galerkin Method

In this paper, the Helmholtz equation for the exterior Neumann boundary condition for the pseudosphere in three dimensions using the global Galerkin method is studied. The Galerkin method will be used to solve Jones’ modified integral equation approach (modified as a series of radiating waves will be added to the fundamental solution) for the Neumann problem for the Helmholtz equation, which uses a series of double sums to approximate the integral. A Fortran 77 program is used and some required subroutines from the Naval Warfare Center are called to help increase ouraccuracy since these boundary integrals are difficult to solve. The solutions obtained arecompared to the true solution for the Neumann problem to understand how well the method converges. The lower errors obtained show that the method for complete reflection of the sound waves off of the pseudosphere is accurate and successful. Also presented in this paper are both computational and theoretical details of the method ofdifferent values of k for the pseudosphere

Efficient Resolution of Anisotropic Structures

We highlight some recent new delevelopments concerning the sparse representation of possibly high-dimensional functions exhibiting strong anisotropic features and low regularity in isotropic Sobolev or Besov scales. Specifically, we focus on the solution of transport equations which exhibit propagation of singularities where, additionally, high-dimensionality enters when the convection field, and hence the solutions, depend on parameters varying over some compact set. Important constituents of our approach are directionally adaptive discretization concepts motivated by compactly supported shearlet systems, and well-conditioned stable variational formulations that support trial spaces with anisotropic refinements with arbitrary directionalities. We prove that they provide tight error-residual relations which are used to contrive rigorously founded adaptive refinement schemes which converge in $L_2$. Moreover, in the context of parameter dependent problems we discuss two approaches serving different purposes and working under different regularity assumptions. For frequent query problems, making essential use of the novel well-conditioned variational formulations, a new Reduced Basis Method is outlined which exhibits a certain rate-optimal performance for indefinite, unsymmetric or singularly perturbed problems. For the radiative transfer problem with scattering a sparse tensor method is presented which mitigates or even overcomes the curse of dimensionality under suitable (so far still isotropic) regularity assumptions. Numerical examples for both methods illustrate the theoretical findings

Higher-order compatible finite element schemes for the nonlinear rotating shallow water equations on the sphere

We describe a compatible finite element discretisation for the shallow water equations on the rotating sphere, concentrating on integrating consistent upwind stabilisation into the framework. Although the prognostic variables are velocity and layer depth, the discretisation has a diagnostic potential vorticity that satisfies a stable upwinded advection equation through a Taylor-Galerkin scheme; this provides a mechanism for dissipating enstrophy at the gridscale whilst retaining optimal order consistency. We also use upwind discontinuous Galerkin schemes for the transport of layer depth. These transport schemes are incorporated into a semi-implicit formulation that is facilitated by a hybridisation method for solving the resulting mixed Helmholtz equation. We illustrate our discretisation with some standard rotating sphere test problems.Comment: accepted versio

Superconvergence of Iterated Solutions for Linear and Nonlinear Integral Equations: Wavelet Applications

In this dissertation, we develop the Petrov-Galerkin method and the iterated Petrov-Galerkin method for a class of nonlinear Hammerstein equation. We also investigate the superconvergence phenomenon of the iterated Petrov-Galerkin and degenerate kernel numerical solutions of linear and nonlinear integral equations with a class of wavelet basis. The Fredholm integral equations and the Hammerstein equations are considered in linear and nonlinear cases respectively. Alpert demonstrated that an application of a class of wavelet basis elements in the Galerkin approximation of the Fredholm equation of the second kind leads to a system of linear equations which is sparse. The main concern of this dissertation is to address the issue of how this sparsity manifests itself in the setting of nonlinear equations, particularly for Hammerstein equations. We demonstrate that sparsity appears in the Jacobian matrix when one attempts to solve the system of nonlinear equations by the Newton\u27s iterative method. Overall, the dissertation generalizes the results of Alpert to nonlinear equations setting as well as the results of Chen and Xu, who discussed the Petrov-Galerkin method for Fredholm equation, to nonlinear equations setting