BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Oscar P. Bruno

Oscar P. Bruno appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAphysics.comp-phNumerical Analysismath.APmath-phmath.MPphysics.flu-dynphysics.plasm-ph

Catalog footprint

What is connected

13works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2022arXiv

Massively Parallelized Interpolated Factored Green Function Method

This paper presents the first parallel implementation of the novel "Interpolated Factored Green Function" (IFGF) method introduced recently for the accelerated evaluation of discrete integral operators arising in wave scattering and other areas (Bauinger and Bruno, Jour. Computat. Phys., 2021). On the basis of the hierarchical IFGF interpolation strategy, the proposed (hybrid MPI-OpenMP) parallel implementation results in highly efficient data communication, and it exhibits in practice excellent parallel scaling up to large numbers of cores - without any hard limitations on the number of cores concurrently employed with high efficiency. Moreover, on any given number of cores, the proposed parallel approach preserves the O(N log N) computing cost inherent in the sequential version of the IFGF algorithm. Unlike other approaches, the IFGF method does not utilize the Fast Fourier Transform (FFT), and it is thus better suited than other methods for efficient parallelization in distributed-memory computer systems. In particular, the IFGF method relies on a "peer-to-peer" strategy wherein, at every level, field propagation is directly enacted via "exchanges" between "peer" polynomials of low and constant degree, without data accumulation in large-scale "telephone-central" mathematical constructs such as those in the Fast Multipole Method (FMM) or pure FFT-based approaches. A variety of numerical results presented in this paper illustrate the character of the proposed parallel algorithm, including excellent weak and strong parallel scaling properties in all cases considered - for problems of up to 4,096 wavelengths in acoustic size, and scaling tests spanning from 1 compute core to all 1,680 cores available in the High Performance Computing cluster used.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Vector potential-based MHD solver for non-periodic flows using Fourier continuation expansions

A high-order method to evolve in time electromagnetic and velocity fields in conducting fluids with non-periodic boundaries is presented. The method has a small overhead compared with fast FFT-based pseudospectral methods in periodic domains. It uses the magnetic vector potential formulation for accurately enforcing the null divergence of the magnetic field, and allowing for different boundary conditions including perfectly conducting walls or vacuum surroundings, two cases relevant for many astrophysical, geophysical, and industrial flows. A spectral Fourier continuation method is used to accurately represent all fields and their spatial derivatives, allowing also for efficient solution of Poisson equations with different boundaries. A study of conducting flows at different Reynolds and Hartmann numbers, and with different boundary conditions, is presented to study convergence of the method and the accuracy of the solenoidal and boundary conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

"Interpolated Factored Green Function" Method for accelerated solution of Scattering Problems

This paper presents a novel {\em Interpolated Factored Green Function} method (IFGF) for the accelerated evaluation of the integral operators in scattering theory and other areas. Like existing acceleration methods in these fields, the IFGF algorithm evaluates the action of Green function-based integral operators at a cost of $\mathcal{O}(N\log N)$ operations for an $N$-point surface mesh. The IFGF strategy, which leads to an extremely simple algorithm, capitalizes on slow variations inherent in a certain Green function {\em analytic factor}, which is analytic up to and including infinity, and which therefore allows for accelerated evaluation of fields produced by groups of sources on the basis of a recursive application of classical interpolation methods. Unlike other approaches, the IFGF method does not utilize the Fast Fourier Transform (FFT), and is thus better suited than other methods for efficient parallelization in distributed-memory computer systems. Only a serial implementation of the algorithm is considered in this paper, however, whose efficiency in terms of memory and speed is illustrated by means of a variety of numerical experiments -- including a 43 min., single-core operator evaluation (on 10 GB of peak memory), with a relative error of $1.5\times 10^{-2}$, for a problem of acoustic size of 512$λ$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Windowed Green Function method for elastic scattering problems on a half-space

This paper presents a windowed Green function (WGF) method for the numerical solution of problems of elastic scattering by "locally-rough surfaces" (i.e., local perturbations of a half space), under either Dirichlet or Neumann boundary conditions, and in both two and three spatial dimensions. The proposed WGF method relies on an integral-equation formulation based on the free-space Green function, together with smooth operator windowing (based on a "slow-rise" windowing function) and efficient high-order singular-integration methods. The approach avoids the evaluation of the expensive layer Green function for elastic problems on a half-space, and it yields uniformly fast convergence for all incident angles. Numerical experiments for both two and three dimensional problems are presented, demonstrating the accuracy and super-algebraically fast convergence of the proposed method as the window-size grows.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fourier continuation method for incompressible fluids with boundaries

We present a Fourier Continuation-based parallel pseudospectral method for incompressible fluids in cuboid non-periodic domains. The method produces dispersionless and dissipationless derivatives with fast spectral convergence inside the domain, and with very high order convergence at the boundaries. Incompressibility is imposed by solving a Poisson equation for the pressure. Being Fourier-based, the method allows for fast computation of spectral transforms. It is compatible with uniform grids (although refined or nested meshes can also be implemented), which in turn allows for explicit time integration at sufficiently high Reynolds numbers. Using a new parallel code named SPECTER we illustrate the method with two problems: channel flow, and plane Rayleigh-Bénard convection under the Boussinesq approximation. In both cases the method yields results compatible with previous studies using other high-order numerical methods, with mild requirements on the time step for stability.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

High-order, Dispersionless "Fast-Hybrid" Wave Equation Solver. Part I: $\mathcal{O}(1)$ Sampling Cost via Incident-Field Windowing and Recentering

This paper proposes a frequency/time hybrid integral-equation method for the time dependent wave equation in two and three-dimensional spatial domains. Relying on Fourier Transformation in time, the method utilizes a fixed (time-independent) number of frequency-domain integral-equation solutions to evaluate, with superalgebraically-small errors, time domain solutions for arbitrarily long times. The approach relies on two main elements, namely, 1) A smooth time-windowing methodology that enables accurate band-limited representations for arbitrarily-long time signals, and 2) A novel Fourier transform approach which, in a time-parallel manner and without causing spurious periodicity effects, delivers numerically dispersionless spectrally-accurate solutions. A similar hybrid technique can be obtained on the basis of Laplace transforms instead of Fourier transforms, but we do not consider the Laplace-based method in the present contribution. The algorithm can handle dispersive media, it can tackle complex physical structures, it enables parallelization in time in a straightforward manner, and it allows for time leaping---that is, solution sampling at any given time $T$ at $\mathcal{O}(1)$-bounded sampling cost, for arbitrarily large values of $T$, and without requirement of evaluation of the solution at intermediate times. The proposed frequency-time hybridization strategy, which generalizes to any linear partial differential equation in the time domain for which frequency-domain solutions can be obtained (including e.g. the time-domain Maxwell equations), and which is applicable in a wide range of scientific and engineering contexts, provides significant advantages over other available alternatives such as volumetric discretization, time-domain integral equations, and convolution-quadrature approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Regularized integral equation methods for elastic scattering problems in three dimensions

This paper presents novel methodologies for the numerical simulation of scattering of elastic waves by both closed and open surfaces in three-dimensional space. The proposed approach utilizes new integral formulations as well as an extension to the elastic context of the efficient high-order singular-integration methods~\cite{BG18} introduced recently for the acoustic case. In order to obtain formulations leading to iterative solvers (GMRES) which converge in small numbers of iterations we investigate, theoretically and computationally, the character of the spectra of various operators associated with the elastic-wave Calderón relation---including some of their possible compositions and combinations. In particular, by relying on the fact that the eigenvalues of the composite operator $NS$ are bounded away from zero and infinity, new uniquely-solvable, low-GMRES-iteration integral formulation for the closed-surface case are presented. The introduction of corresponding low-GMRES-iteration equations for the open-surface equations additionally requires, for both spectral quality as well as accuracy and efficiency, use of weighted versions of the classical integral operators to match the singularity of the unknown density at edges. Several numerical examples demonstrate the accuracy and efficiency of the proposed methodology.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Superalgebraically Convergent Smoothly-Windowed Lattice Sums for Doubly Periodic Green Functions in Three-Dimensional Space

This paper, Part I in a two-part series, presents (i) A simple and highly efficient algorithm for evaluation of quasi-periodic Green functions, as well as (ii) An associated boundary-integral equation method for the numerical solution of problems of scattering of waves by doubly periodic arrays of scatterers in three-dimensional space. Except for certain "Wood frequencies" at which the quasi-periodic Green function ceases to exist, the proposed approach, which is based on use of smooth windowing functions, gives rise to lattice sums which converge superalgebraically fast--that is, faster than any power of the number of terms used--in sharp contrast with the extremely slow convergence exhibited by the corresponding sums in absence of smooth windowing. (The Wood-frequency problem is treated in Part II.) A proof presented in this paper establishes rigorously the superalgebraic convergence of the windowed lattice sums. A variety of numerical results demonstrate the practical efficiency of the proposed approach.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-order integral equation methods for problems of scattering by bumps and cavities on half-planes

This paper presents high-order integral equation methods for evaluation of electromagnetic wave scattering by dielectric bumps and dielectric cavities on perfectly conducting or dielectric half-planes. In detail, the algorithms introduced in this paper apply to eight classical scattering problems, namely: scattering by a dielectric bump on a perfectly conducting or a dielectric half-plane, and scattering by a filled, overfilled or void dielectric cavity on a perfectly conducting or a dielectric half-plane. In all cases field representations based on single-layer potentials for appropriately chosen Green functions are used. The numerical far fields and near fields exhibit excellent convergence as discretizations are refined--even at and around points where singular fields and infinite currents exist.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A generalized Calderon Formula for open-arc diffraction problems: theoretical considerations

We deal with the general problem of scattering by open-arcs in two-dimensional space. We show that this problem can be solved by means of certain second-kind integral equations of the form $\tilde{N} \tilde{S}[φ] = f$, where $\tilde{N}$ and $\tilde{S}$ are first-kind integral operators whose composition gives rise to a generalized Calderón formula of the form $\tilde{N} \tilde{S} = \tilde{J}_0^τ+ \tilde{K}$ in a {\em weighted, periodized} Sobolev space. The $\tilde{N} \tilde{S}$ formulation provides, for the first time, a second-kind integral equation for the open-arc scattering problem with Neumann boundary conditions. Numerical experiments show that, for both the Dirichlet and Neumann boundary conditions, our second-kind integral equations have spectra that are bounded away from zero and infinity as $k\to \infty$; to the authors' knowledge these are the first integral equations for these problems that possess this desirable property. Our proofs rely on three main elements: 1) Algebraic manipulations enabled by the presence of integral weights; 2) Use of the classical result of continuity of the Cesàro operator; and 3) Explicit characterization of the point spectrum of $\tilde{J}^τ_0$, which, interestingly, can be decomposed into the union of a countable set and an open set, both tightly clustered around -1/4. As shown in a separate contribution, the new approach can be used to construct simple spectrally-accurate numerical solvers and, when used in conjunction with Krylov-subspace solvers such as GMRES, gives rise to dramatic reductions of Krylov-subspace iteration numbers vs. those required by other approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A high-order integral solver for scalar problems of diffraction by screens and apertures in three dimensional space

We present a novel methodology for the numerical solution of problems of diffraction by infinitely thin screens in three dimensional space. Our approach relies on new integral formulations as well as associated high-order quadrature rules. The new integral formulations involve weighted versions of the classical integral operators associated with the thin-screen Dirichlet and Neumann problems as well as a generalization to the open surface problem of the classical Calderon formulae. The high-order quadrature rules we introduce for these operators, in turn, resolve the multiple Green function and edge singularities (which occur at arbitrarily close distances from each other, and which include weakly singular as well as hypersingular kernels) and thus give rise to super-algebraically fast convergence as the discretization sizes are increased. When used in conjunction with Krylov-subspace linear algebra solvers such as GMRES, the resulting solvers produce results of high accuracy in small numbers of iterations for low and high frequencies alike. We demonstrate our methodology with a variety of numerical results for screen and aperture problems at high frequencies---including simulation of classical experiments such as the diffraction by a circular disc (including observation of the famous Poisson spot), interference fringes resulting from diffraction across two nearby circular apertures, as well as more complex geometries consisting of multiple scatterers and cavities.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Convergence analysis of a high-order Nystrom integral-equation method for surface scattering problems

In this paper we present a convergence analysis for the Nystrom method proposed in [Jour. Comput. Phys. 169 pp. 2921-2934, 2001] for the solution of the combined boundary integral equation formulations of sound-soft acoustic scattering problems in three-dimensional space. This fast and efficient scheme combines FFT techniques and a polar change of variables that cancels out the kernel singularity. We establish the stability of the algorithms in the $L^2$ norm and we derive convergence estimates in both the $L^2$ and $L^\infty$ norms. In particular, our analysis establishes theoretically the previously observed super-algebraic convergence of the method in cases in which the right-hand side is smooth.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Second-Kind integral solvers for TE and TM problems of diffraction by open-arcs

We present a novel approach for the numerical solution of problems of diffraction by open arcs in two dimensional space. Our methodology relies on composition of {\em weighted versions} of the classical integral operators associated with the Dirichlet and Neumann problems (TE and TM polarizations, respectively) together with a generalization to the open-arc case of the well known closed-surface Calderón formulae. When used in conjunction with spectrally accurate discretization rules and Krylov-subspace linear algebra solvers such as GMRES, the new second-kind TE and TM formulations for open arcs produce results of high accuracy in small numbers of iterations and short computing times---for low and high frequencies alike.

0 citations0 reviewsNo code linkedOpen access