We apply the ultraspherical spectral method to solving time-dependent PDEs by proposing two approaches to discretization based on the method of lines and show that these approaches produce approximately same results. We analyze the stability, the error, and the computational cost of the proposed method. In addition, we show how adaptivity can be incorporated to offer adequate spatial resolution efficiently. Both linear and nonlinear problems are considered. We also explore time integration using exponential integrators with the ultraspherical spatial discretization. Comparisons with the Chebyshev pseudospectral method are given along the discussion and they show that the ultraspherical spectral method is a competitive candidate for the spatial discretization of time-dependent PDEs.
We discuss applications of exact structures and relative homological algebra to the study of invariants of multiparameter persistence modules. This paper is mostly expository, but does contain a pair of novel results. Over finite posets, classical arguments about the relative projective modules of an exact structure make use of Auslander-Reiten theory. One of our results establishes a new adjunction which allows us to "lift" these arguments to certain infinite posets over which Auslander-Reiten theory is not available. We give several examples of this lifting, in particular highlighting the non-existence and existence of resolutions by upsets when working with finitely presentable representations of the plane and of the closure of the positive quadrant, respectively. We then restrict our attention to finite posets. In this setting, we discuss the relationship between the global dimension of an exact structure and the representation dimension of the incidence algebra of the poset. We conclude with our second novel contribution. This is an explicit description of the irreducible morphisms between relative projective modules for several exact structures which have appeared previously in the literature.
We study the problem of enumerating Tarski fixed points, focusing on the relational lattices of equivalences, quasiorders and binary relations. We present a polynomial space enumeration algorithm for Tarski fixed points on these lattices and other lattices of polynomial height. It achieves polynomial delay when enumerating fixed points of increasing isotone maps on all three lattices, as well as decreasing isotone maps on the lattice of binary relations. In those cases in which the enumeration algorithm does not guarantee polynomial delay on the three relational lattices on the other hand, we prove exponential lower bounds for deciding the existence of three fixed points when the isotone map is given as an oracle, and that it is NP-hard to find three or more Tarski fixed points. More generally, we show that any deterministic or bounded-error randomized algorithm must perform a number of queries asymptotically at least as large as the lattice width to decide the existence of three fixed points when the isotone map is given as an oracle. Finally, we demonstrate that our findings yield a polynomial delay and space algorithm for listing bisimulations and instances of some related models of behavioral or role equivalence.
This work is the first exploration of proof-theoretic semantics for a substructural logic. It focuses on the base-extension semantics (B-eS) for intuitionistic multiplicative linear logic (IMLL). The starting point is a review of Sandqvist's B-eS for intuitionistic propositional logic (IPL), for which we propose an alternative treatment of conjunction that takes the form of the generalized elimination rule for the connective. The resulting semantics is shown to be sound and complete. This motivates our main contribution, a B-eS for IMLL, in which the definitions of the logical constants all take the form of their elimination rule and for which soundness and completeness are established.
We develop the no-propagate algorithm for sampling the linear response of random dynamical systems, which are non-uniform hyperbolic deterministic systems perturbed by noise with smooth density. We first derive a Monte-Carlo type formula and then the algorithm, which is different from the ensemble (stochastic gradient) algorithms, finite-element algorithms, and fast-response algorithms; it does not involve the propagation of vectors or covectors, and only the density of the noise is differentiated, so the formula is not cursed by gradient explosion, dimensionality, or non-hyperbolicity. We demonstrate our algorithm on a tent map perturbed by noise and a chaotic neural network with 51 layers $\times$ 9 neurons. By itself, this algorithm approximates the linear response of non-hyperbolic deterministic systems, with an additional error proportional to the noise. We also discuss the potential of using this algorithm as a part of a bigger algorithm with smaller error.
In the context of adaptive remeshing, the virtual element method provides significant advantages over the finite element method. The attractive features of the virtual element method, such as the permission of arbitrary element geometries, and the seamless permission of 'hanging' nodes, have inspired many works concerning error estimation and adaptivity. However, these works have primarily focused on adaptive refinement techniques with little attention paid to adaptive coarsening (i.e. de-refinement) techniques that are required for the development of fully adaptive remeshing procedures. In this work novel indicators are proposed for the identification of patches/clusters of elements to be coarsened, along with a novel procedure to perform the coarsening. The indicators are computed over prospective patches of elements rather than on individual elements to identify the most suitable combinations of elements to coarsen. The coarsening procedure is robust and suitable for meshes of structured and unstructured/Voronoi elements. Numerical results demonstrate the high degree of efficacy of the proposed coarsening procedures and sensible mesh evolution during the coarsening process. It is demonstrated that critical mesh geometries, such as non-convex corners and holes, are preserved during coarsening, and that meshes remain fine in regions of interest to engineers, such as near singularities.
The use of numerical based multi-phase fluid flow simulation can significantly aid in the development of an effective remediation strategy for groundwater systems contaminated with Dense Non Aqueous Phase Liquid (DNAPL). Incorporating the lithological heterogeneities of the aquifer into the model domain is a crucial aspect in the development of robust numerical simulators. Previous research studies have attempted to incorporate lithological heterogeneities into the domain; however, most of these numerical simulators are based on Finite Volume Method (FVM) and Finite Difference Method (FDM) which have limited applicability in the field-scale aquifers. Finite Element Method (FEM) can be highly useful in developing the field-scale simulation of DNAPL infiltration due to its consistent accuracy on irregular study domain, and the availability of higher orders of basis functions. In this research work, FEM based model has been developed to simulate the DNAPL infiltration in a hypothetical field-scale aquifer. The model results demonstrate the effect of meso-scale heterogeneities, specifically clay lenses, on the migration and accumulation of Dense Non Aqueous Phase Liquid (DNAPL) within the aquifer. Furthermore, this research provides valuable insights for the development of an appropriate remediation strategy for a general contaminated aquifer.
It is a challenge to numerically solve nonlinear partial differential equations whose solution involves discontinuity. In the context of numerical simulators for multi-phase flow in porous media, there exists a long-standing issue known as Grid Orientation Effect (GOE), wherein different numerical solutions can be obtained when considering grids with different orientations under certain unfavorable conditions. Our perspective is that GOE arises due to numerical instability near displacement fronts, where spurious oscillations accompanied by sharp fronts, if not adequately suppressed, lead to GOE. To reduce or even eliminate GOE, we propose augmenting adaptive artificial viscosity when solving the saturation equation. It has been demonstrated that appropriate artificial viscosity can effectively reduce or even eliminate GOE. The proposed numerical method can be easily applied in practical engineering problems.
We propose a Hermite spectral method for the inelastic Boltzmann equation, which makes two-dimensional periodic problem computation affordable by the hardware nowadays. The new algorithm is based on a Hermite expansion, where the expansion coefficients for the VHS model are reduced into several summations and can be derived exactly. Moreover, a new collision model is built with a combination of the quadratic collision operator and a linearized collision operator, which helps us to balance the computational cost and the accuracy. Various numerical experiments, including spatially two-dimensional simulations, demonstrate the accuracy and efficiency of this numerical scheme.
This paper presents a robust numerical solution to the electromagnetic scattering problem involving multiple multi-layered cavities in both transverse magnetic and electric polarizations. A transparent boundary condition is introduced at the open aperture of the cavity to transform the problem from an unbounded domain into that of bounded cavities. By employing Fourier series expansion of the solution, we reduce the original boundary value problem to a two-point boundary value problem, represented as an ordinary differential equation for the Fourier coefficients. The analytical derivation of the connection formula for the solution enables us to construct a small-scale system that includes solely the Fourier coefficients on the aperture, streamlining the solving process. Furthermore, we propose accurate numerical quadrature formulas designed to efficiently handle the weakly singular integrals that arise in the transparent boundary conditions. To demonstrate the effectiveness and versatility of our proposed method, a series of numerical experiments are conducted.
The purpose of the paper is to provide a characterization of the error of the best polynomial approximation of composite functions in weighted spaces. Such a characterization is essential for the convergence analysis of numerical methods applied to non-linear problems or for numerical approaches that make use of regularization techniques to cure low smoothness of the solution. This result is obtained through an estimate of the derivatives of composite functions in weighted uniform norm.