Related papers: Accurate and Efficient Nystrom Volume Integral Equ…
This paper presents a direct numerical scheme to approximate the solution of all classes of nonlinear Volterra integral equations of the first kind. This computational method is based on operational matrices and vectors. The operational…
This paper is the direct-formulation companion to [Burbano-Gallegos, P\'erez-Arancibia, and Turc, ESAIM: M2AN, 60(1):273--315, 2026], which developed indirect combined-field-only boundary integral equations (BIEs) for time-harmonic…
Validity of fluid models breaks down for non-thermal or weakly collisional plasmas which often occur e.g. in the solar wind. In these regimes one has to resort to modelling through the first-principle Vlasov-Maxwell system, but its…
The paper describes a numerical method for solving acoustic multibody scattering problems in two and three dimensions. The idea is to compute a highly accurate approximation to the scattering operator for each body through a local…
A numerical algorithm for solving mantle convection problems with strongly variable viscosity is presented. Equations for conservation of mass and momentum for highly viscous and incompressible fluids are solved iteratively by a multigrid…
The virtual element method (VEM) allows discretization of elasticity and plasticity problems with polygons in 2D and polyhedrals in 3D. The polygons (and polyhedrals) can have an arbitrary number of sides and can be concave or convex. These…
Multiscale finite element methods for 2D/1D problems have been studied in this work to demonstrate their excellent ability to solve real-world problems. These methods are much more efficient than conventional 3D finite element methods and…
A paramount goal in the field of nuclear physics is to unify ab-initio treatments of bound and unbound states. The position-space quantum Monte Carlo (QMC) methods have a long history of successful bound state calculations in light systems…
In the present paper we are concerned with the Novikov--Veselov equation at negative energy, i.e. with the $ (2 + 1) $--dimensional analog of the KdV equation integrable by the method of inverse scattering for the two--dimensional…
This hybrid method (FE-DVR), introduced by Resigno and McCurdy, Phys. Rev. A 62, 032706 (2000), uses Lagrange polynomials in each partition, rather than "hat" functions or Gaussian functions. These polynomials are discrete variable…
Recently, the volume integral equation (VIE) approach has been proposed as an efficient simulation tool for silicon photonics applications [J. Lightw. Technol. 36, 3765 (2018)]. However, for the high-frequency and strong contrast problems…
We present a software suite for the analysis and optimization of ideal convex polyhedra in hyperbolic 3-space $\mathbb{H}^3$. Using Rivin's variational characterization of ideal polyhedra, we develop efficient algorithms for checking…
High frequency integral equation methodologies display the capability of reproducing single-scattering returns in frequency-independent computational times and employ a Neumann series formulation to handle multiple-scattering effects. This…
This paper introduces a novel method for the efficient and accurate computation of volume fractions on unstructured polyhedral meshes, where the phase boundary is an orientable hypersurface, implicitly given as the iso-contour of a…
We consider the numerical algorithm for the two-dimensional time-harmonic elastic wave scattering by unbounded rough surfaces with Dirichlet boundary condition. A Nystr\"{o}m method is proposed for the scattering problem based on the…
A mathematical model for description of the viscous fingering induced by a chemical reaction is under study. This complicated five-component model is reduced to a three-component diffusive Lotka-Volterra system with convection by…
We develop a new, efficient, and accurate method to simulate frequency-domain borehole electromagnetic (EM) measurements acquired in the presence of three-dimensional (3D) variations of the anisotropic subsurface conductivity. The method is…
Surface integral equation (SIE) methods are of great interest for the numerical solution of Maxwell's equations in the presence of homogeneous objects. However, existing SIE algorithms have limitations, either in terms of scalability,…
A butterfly-based fast direct integral equation solver for analyzing high-frequency scattering from two-dimensional objects is presented. The solver leverages a randomized butterfly scheme to compress blocks corresponding to near- and…
Multiple scattering methods are widely used to reduce the computational complexity of acoustic or electromagnetic scattering problems when waves propagate through media containing many identical inclusions. Historically, this numerical…