The 3D Laplace equation is one of the important PDEs of physics and describes the phenomonology of electrostatics and magnetostatics. Frequently very precise solution of this PDE is required; but with conventional finite element or finite difference codes this is difficult to achieve because of the need for an exceedingly fine mesh which leads to often prohibitive CPU time. We present an altern...

The concept of topological sensitivity has been successfully employed as an imaging tool to obtain the correct initial topology and preliminary geometry of hidden obstacles for a variety of inverse scattering problems. In this paper, we extend these ideas to acoustic scattering involving transient waveforms and penetrable obstacles. Through a boundary integral equation framework, we present a d...

Finite element and finite difference methods have been widely used, among other methods, to numerically solve the Fokker-Planck equation for investigating the time history of the probability density function of linear and nonlinear 2d and 3d problems, and also the application to 4d problems has been addressed. However, due to the enormous increase of the computational costs, different strategie...

The thin plate is a key structure in various musical instruments, including many percussion instruments and the soundboard of the piano, and also is the mechanism underlying electromechanical plate reverberation. As such, it is a suitable candidate for physical modelling approaches to audio effects and sound synthesis, such as finite difference methods—though great attention must be paid to the...

Many applications based on finite element and finite difference methods include the solution of large sparse linear systems using preconditioned iterative methods. Matrix vector multiplication is one of the key operations that has a significant impact on the performance of any iterative solver. In this paper, recent developments in sparse storage formats on vector machines are reviewed. Then, s...

Most results related to discrete nonnegativity conservation principles (DNCP) for elliptic problems are limited to finite differences (FDM) and lowest-order finite element methods (FEM). In this paper we confirm that a straightforward extension to higher-order finite element methods (hpFEM) in the classical sense is not possible. We formulate a weaker DNCP for the Poisson equation in one spatia...

The prevalent model to describe propagation of neural membrane voltage responses to ion currents is the cable equation. The cable equation is a parabolic-type partial differential equation, which is frequently studied with the use of compartmental models. This paper explores other solution techniques for this type of partial differential equation, specifically the finite difference and finite e...