Numerics of boundary-domain integral and integro-differential equations for BVP with variable coefficient in 3D

This is the post-print version of the article. The official published version can be accessed from the links below - Copyright @ 2013 Springer-VerlagA numerical implementation of the direct boundary-domain integral and integro-differential equations, BDIDEs, for treatment of the Dirichlet problem for a scalar elliptic PDE with variable coefficient in a three-dimensional domain is discussed. The mesh-based discretisation of the BDIEs with tetrahedron domain elements in conjunction with collocation method leads to a system of linear algebraic equations (discretised BDIE). The involved fully populated matrices are approximated by means of the H-Matrix/adaptive cross approximation technique. Convergence of the method is investigated.This study is partially supported by the EPSRC grant EP/H020497/1:"Mathematical Analysis of Localised-Boundary-Domain Integral Equations for Variable-Coefficients Boundary Value Problems"

Complexity Analysis of a Fast Directional Matrix-Vector Multiplication

We consider a fast, data-sparse directional method to realize matrix-vector products related to point evaluations of the Helmholtz kernel. The method is based on a hierarchical partitioning of the point sets and the matrix. The considered directional multi-level approximation of the Helmholtz kernel can be applied even on high-frequency levels efficiently. We provide a detailed analysis of the almost linear asymptotic complexity of the presented method. Our numerical experiments are in good agreement with the provided theory.Comment: 20 pages, 2 figures, 1 tabl

Deviational simulation of phonon transport in graphene ribbons with ab initio scattering

We present a deviational Monte Carlo method for solving the Boltzmann-Peierls equation with ab initio 3-phonon scattering, for temporally and spatially dependent thermal transport problems in arbitrary geometries. Phonon dispersion relations and transition rates for graphene are obtained from density functional theory calculations. The ab initio scattering operator is simulated by an energy-conserving stochastic algorithm embedded within a deviational, low-variance Monte Carlo formulation. The deviational formulation ensures that simulations are computationally feasible for arbitrarily small temperature differences, while the stochastic treatment of the scattering operator is both efficient and exhibits no timestep error. The proposed method, in which geometry and phonon-boundary scattering are explicitly treated, is extensively validated by comparison to analytical results, previous numerical solutions and experiments. It is subsequently used to generate solutions for heat transport in graphene ribbons of various geometries and evaluate the validity of some common approximations found in the literature. Our results show that modeling transport in long ribbons of finite width using the homogeneous Boltzmann equation and approximating phonon-boundary scattering using an additional homogeneous scattering rate introduces an error on the order of 10% at room temperature, with the maximum deviation reaching 30% in the middle of the transition regime.Singapore-MIT Alliance for Research and TechnologyAmerican Society for Engineering Education. National Defense Science and Engineering Graduate FellowshipNational Science Foundation (U.S.). Graduate Research Fellowshi

Current large deviations in a driven dissipative model

We consider lattice gas diffusive dynamics with creation-annihilation in the bulk and maintained out of equilibrium by two reservoirs at the boundaries. This stochastic particle system can be viewed as a toy model for granular gases where the energy is injected at the boundary and dissipated in the bulk. The large deviation functional for the particle currents flowing through the system is computed and some physical consequences are discussed: the mechanism for local current fluctuations, dynamical phase transitions, the fluctuation-relation

A posteriori error estimates for the virtual element method

An a posteriori error analysis for the virtual element method (VEM) applied to general elliptic problems is presented. The resulting error estimator is of residual-type and applies on very general polygonal/polyhedral meshes. The estimator is fully computable as it relies only on quantities available from the VEM solution, namely its degrees of freedom and element-wise polynomial projection. Upper and lower bounds of the error estimator with respect to the VEM approximation error are proven. The error estimator is used to drive adaptive mesh refinement in a number of test problems. Mesh adaptation is particularly simple to implement since elements with consecutive co-planar edges/faces are allowed and, therefore, locally adapted meshes do not require any local mesh post-processing