A boundary integral method for modelling vibroacoustic energy distributions in uncertain built up structures

A phase-space boundary integral method is developed for modelling stochastic high-frequency acoustic and vibrational energy transport in both single and multi-domain problems. The numerical implementation is carried out using the collocation method in both the position and momentum phase-space variables. One of the major developments of this work is the systematic convergence study, which demonstrates that the proposed numerical schemes exhibit convergence rates that could be expected from theoretical estimates under the right conditions. For the discretisation with respect to the momentum variable, we employ spectrally convergent basis approximations using both Legendre polynomials and Gaussian radial basis functions. The former have the advantage of being simpler to apply in general without the need for preconditioning techniques. The Gaussian basis is introduced with the aim of achieving more efficient computations in the weak noise case with near-deterministic dynamics. Numerical results for a series of coupled domain problems are presented, and demonstrate the potential for future applications to larger scale problems from industry

Preconditioning for radial basis function partition of unity methods

Meshfree radial basis function (RBF) methods are of interest for solving partial differential equations due to attractive convergence properties, flexibility with respect to geometry, and ease of implementation. For global RBF methods, the computational cost grows rapidly with dimension and problem size, so localised approaches, such as partition of unity or stencil based RBF methods, are currently being developed. An RBF partition of unity method (RBF--PUM) approximates functions through a combination of local RBF approximations. The linear systems that arise are locally unstructured, but with a global structure due to the partitioning of the domain. Due to the sparsity of the matrices, for large scale problems, iterative solution methods are needed both for computational reasons and to reduce memory requirements. In this paper we implement and test different algebraic preconditioning strategies based on the structure of the matrix in combination with incomplete factorisations. We compare their performance for different orderings and problem settings and find that a no-fill incomplete factorisation of the central band of the original discretisation matrix provides a robust and efficient preconditioner

Sparse spectral-tau method for the three-dimensional helically reduced wave equation on two-center domains

We describe a multidomain spectral-tau method for solving the three-dimensional helically reduced wave equation on the type of two-center domain that arises when modeling compact binary objects in astrophysical applications. A global two-center domain may arise as the union of Cartesian blocks, cylindrical shells, and inner and outer spherical shells. For each such subdomain, our key objective is to realize certain (differential and multiplication) physical-space operators as matrices acting on the corresponding set of modal coefficients. We achieve sparse banded realizations through the integration "preconditioning" of Coutsias, Hagstrom, Hesthaven, and Torres. Since ours is the first three-dimensional multidomain implementation of the technique, we focus on the issue of convergence for the global solver, here the alternating Schwarz method accelerated by GMRES. Our methods may prove relevant for numerical solution of other mixed-type or elliptic problems, and in particular for the generation of initial data in general relativity.Comment: 37 pages, 3 figures, 12 table

Multidomain Spectral Method for the Helically Reduced Wave Equation

We consider the 2+1 and 3+1 scalar wave equations reduced via a helical Killing field, respectively referred to as the 2-dimensional and 3-dimensional helically reduced wave equation (HRWE). The HRWE serves as the fundamental model for the mixed-type PDE arising in the periodic standing wave (PSW) approximation to binary inspiral. We present a method for solving the equation based on domain decomposition and spectral approximation. Beyond describing such a numerical method for solving strictly linear HRWE, we also present results for a nonlinear scalar model of binary inspiral. The PSW approximation has already been theoretically and numerically studied in the context of the post-Minkowskian gravitational field, with numerical simulations carried out via the "eigenspectral method." Despite its name, the eigenspectral technique does feature a finite-difference component, and is lower-order accurate. We intend to apply the numerical method described here to the theoretically well-developed post-Minkowski PSW formalism with the twin goals of spectral accuracy and the coordinate flexibility afforded by global spectral interpolation.Comment: 57 pages, 11 figures, uses elsart.cls. Final version includes revisions based on referee reports and has two extra figure

hp-adaptive discontinuous Galerkin solver for elliptic equations in numerical relativity

A considerable amount of attention has been given to discontinuous Galerkin methods for hyperbolic problems in numerical relativity, showing potential advantages of the methods in dealing with hydrodynamical shocks and other discontinuities. This paper investigates discontinuous Galerkin methods for the solution of elliptic problems in numerical relativity. We present a novel hp-adaptive numerical scheme for curvilinear and non-conforming meshes. It uses a multigrid preconditioner with a Chebyshev or Schwarz smoother to create a very scalable discontinuous Galerkin code on generic domains. The code employs compactification to move the outer boundary near spatial infinity. We explore the properties of the code on some test problems, including one mimicking Neutron stars with phase transitions. We also apply it to construct initial data for two or three black holes

Fast Mesh Refinement in Pseudospectral Optimal Control

Mesh refinement in pseudospectral (PS) optimal control is embarrassingly easy --- simply increase the order $N$ of the Lagrange interpolating polynomial and the mathematics of convergence automates the distribution of the grid points. Unfortunately, as $N$ increases, the condition number of the resulting linear algebra increases as $N^2$; hence, spectral efficiency and accuracy are lost in practice. In this paper, we advance Birkhoff interpolation concepts over an arbitrary grid to generate well-conditioned PS optimal control discretizations. We show that the condition number increases only as $\sqrt{N}$ in general, but is independent of $N$ for the special case of one of the boundary points being fixed. Hence, spectral accuracy and efficiency are maintained as $N$ increases. The effectiveness of the resulting fast mesh refinement strategy is demonstrated by using \underline{polynomials of over a thousandth order} to solve a low-thrust, long-duration orbit transfer problem.Comment: 27 pages, 12 figures, JGCD April 201

Approximation of the Neutron Diffusion Equation on Hexagonal Geometries

La ecuación de la difusión neutrónica describe la población de neutrones de un reactor nuclear. Este trabajo trata con este modelo para reactores nucleares con geometría hexagonal. En primer lugar se estudia la ecuación de la difusión neutrónica. Este es un problema diferencial de valores propios, llamado problema de los modos Lambda. Para resolver el problema de los modos Lambda se han comparado diferentes métodos en geometrías unidimensionales, resultando como el mejor el método de elementos espectrales. Usando este método discretizamos los operadores en geometrías bidimensiones y tridimensionales, resolviendo el problema algebraica de valores propios resultante con el método de Arnoldi. La distribución de neutrones estado estacionario se utiliza como condición inicial para la integración de la ecuación de la difusión neutrónica dependiente del tiempo. Se utiliza un método de Euler implícito para integrar en el tiempo. Cuando un nodo está parcialmente insertado aparece un comportamiento no físico de la solución, el efecto ``rod cusping'', que se corrige mediante la ponderación de las secciones eficaces con el flujo del paso de tiempo anterior. Cuando la solución de los sistemas algebraicos que surgen en el método hacia atrás, un método de Krylov se utiliza para resolver los sistemas resultantes, y diferentes estrategias de precondicionamiento se evalúan se. La primera consiste en el uso de la estructura de bloque obtenido por los grupos de energía para resolver el sistema por bloques, y diferentes técnicas de aceleración para el esquema iterativo de bloques y un precondicionador utilizando esta estructura de bloque se proponen. Además se estudia un precondicionador espectral, que hace uso de la información en un subespacio de Krylov para precondicionar el siguiente sistema. También se proponen métodos exponenciales de segundo y cuarto orden integrar la ecuación de difusión neutrónica dependiente del tiempo, donde la exponencial de la matriz del sistema tiene quGonzález Pintor, S. (2012). Approximation of the Neutron Diffusion Equation on Hexagonal Geometries [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/17829Palanci

On certain (block) Toeplitz matrices related to radial functions

AbstractInterpolation of smooth functions and the discretization of elliptic PDEs by means of radial functions lead to structured linear systems which, for equidistant grid points, have almost the (block) Toeplitz structure. We prove upper bounds for the condition numbers of the n×n Toeplitz matrices which discretize the model problem u″(x)=f(x), x∈(0,1), u(0)=a, u(1)=b over an equally spaced grid of n+2 points in [0,1] by means of the collocation method based on radial functions of the multiquadric, inverse multiquadric and Gaussian type. These bounds are asymptotically sharp