Researcher profile

Balasubramaniam Shanker

Balasubramaniam Shanker contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

9 published item(s)

preprint2020arXiv

High Performance Evaluation of Helmholtz Potentials using the Multi-Level Fast Multipole Algorithm

Evaluation of pair potentials is critical in a number of areas of physics. The classicalN-body problem has its root in evaluating the Laplace potential, and has spawned tree-algorithms, the fast multipole method (FMM), as well as kernel independent approaches. Over the years, FMM for Laplace potential has had a profound impact on a number of disciplines as it has been possible to develop highly scalable parallel algorithm for these potential evaluators. This is in stark contrast to parallel algorithms for the Helmholtz (oscillatory) potentials. The principal bottleneck to scalable parallelism are operations necessary to traverse up, across and down the tree, affecting both computation and communication. In this paper, we describe techniques to overcome bottlenecks and achieve high performance evaluation of the Helmholtz potential for a wide spectrum of geometries. We demonstrate that the resulting implementation has a load balancing effect that significantly reduces the time-to-solution and enhances the scale of problems that can be treated using full wave physics.

preprint2015arXiv

Computational dynamics of acoustically-driven microsphere systems

We propose a computational framework for the self-consistent dynamics of a microsphere system driven by a pulsed acoustic field in an ideal fluid. Our framework combines a molecular dynamics integrator describing the dynamics of the microsphere system with a time-dependent integral equation solver for the acoustic field that makes use of fields represented as surface expansions in spherical harmonic basis functions. The presented approach allows us to describe the inter-particle interaction induced by the field as well as the dynamics of trapping in counter-propagating acoustic pulses. The integral equation formulation leads to equations of motion for the microspheres describing the effect of non-dissipative drag forces. We show (1) that the field-induced interactions between the microspheres give rise to effective dipolar interactions, with effective dipoles defined by their velocities, and (2) that the dominant effect of an ultrasound pulse through a cloud of microspheres gives rise mainly to a translation of the system, though we also observe both expansion and contraction of the cloud determined by the initial system geometry.

preprint2015arXiv

Subdivision based Isogeometric Analysis technique for Electric Field Integral Equations for Simply Connected Structures

The analysis of electromagnetic scattering has long been performed on a discrete representation of the geometry. This representation is typically continuous but {\em not} differentiable. The need to define physical quantities on this geometric representation has led to development of sets of basis functions that need to satisfy constraints at the boundaries of the elements/tesselations (viz., continuity of normal or tangential components across element boundaries). For electromagnetics, these result in either curl/div-conforming basis sets. The geometric representation used for analysis is in stark contrast with that used for design, wherein the surface representation is higher order differentiable. Using this representation for {\em both} geometry and physics on geometry has several advantages, and is eludicated in Hughes et al., Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement, Computer Methods in Applied Mechanics and Engineering 194 (39-41) (2005). Until now, a bulk of the literature on isogeometric methods have been limited to solid mechanics, with some effort to create NURBS based basis functions for electromagnetic analysis. In this paper, we present the first complete isogeometry solution methodology for the electric field integral equation as applied to simply connected structures. This paper systematically proceeds through surface representation using subdivision, definition of vector basis functions on this surface, to fidelity in the solution of integral equations. We also present techniques to stabilize the solution at low frequencies, and impose a Calderón preconditioner. Several results presented serve to validate the proposed approach as well as demonstrate some of its capabilities.

preprint2014arXiv

A Discontinuous Galerkin Time Domain Framework for Periodic Structures Subject To Oblique Excitation

A nodal Discontinuous Galerkin (DG) method is derived for the analysis of time-domain (TD) scattering from doubly periodic PEC/dielectric structures under oblique interrogation. Field transformations are employed to elaborate a formalism that is free from any issues with causality that are common when applying spatial periodic boundary conditions simultaneously with incident fields at arbitrary angles of incidence. An upwind numerical flux is derived for the transformed variables, which retains the same form as it does in the original Maxwell problem for domains without explicitly imposed periodicity. This, in conjunction with the amenability of the DG framework to non-conformal meshes, provides a natural means of accurately solving the first order TD Maxwell equations for a number of periodic systems of engineering interest. Results are presented that substantiate the accuracy and utility of our method.

preprint2014arXiv

The Generalized Method of Moments for Electromagnetic Boundary Integral Equations

The Generalized Method of Moments (GMM) is a partition of unity based technique for solving electromagnetic and acoustic boundary integral equations. Past work on the GMM for electromagnetics was confined to geometries modeled by piecewise flat tessellations and suffered from spurious internal line charges. In the present article, we redesign the GMM scheme and demonstrate its ability to model scattering from PEC scatterers composed of mixtures of smooth and non-smooth geometrical features. Furthermore, we demonstrate that because the partition of unity provides function and effective geometrical continuity between patches, the GMM permits mixtures of local geometry descriptions and approximation function spaces with significantly more freedom than traditional moment methods.

preprint2014arXiv

Time-dependent Debye-Mie Series Solutions for Electromagnetic Scattering

Frequency domain Mie solutions to scattering from spheres have been used for a long time. However, deriving their transient analogue is a challenge as it involves an inverse Fourier transform of the spherical Hankel functions (and their derivatives) that are convolved with inverse Fourier transforms of spherical Bessel functions (and their derivatives). Series expansion of these convolutions are highly oscillatory (therefore, poorly convergent) and unstable. Indeed, the literature on numerical computation of this convolution is very sparse. In this paper, we present a novel quasi-analytical approach to computing transient Mie scattering that is both stable and rapidly convergent. The approach espoused here is to use vector tesseral harmonics as basis function for the currents in time domain integral equations together with a novel addition theorem for the Green's functions that renders these expansions stable. This procedure results in an orthogonal, spatially-meshfree and singularity-free system, giving us a set of one dimensional Volterra Integral equations. Time-dependent multipole coefficients for each mode are obtained via a time marching procedure. Finally, several numerical examples are presented to show the accuracy and stability the proposed method.

preprint2011arXiv

An O(N) Method for Rapidly Computing Periodic Potentials Using Accelerated Cartesian Expansions

The evaluation of long-range potentials in periodic, many-body systems arises as a necessary step in the numerical modeling of a multitude of interesting physical problems. Direct evaluation of these potentials requires O(N^2) operations and O(N^2) storage, where N is the number of interacting bodies. In this work, we present a method, which requires O(N) operations and O(N) storage, for the evaluation of periodic Helmholtz, Coulomb, and Yukawa potentials with periodicity in 1-, 2-, and 3-dimensions, using the method of Accelerated Cartesian Expansions (ACE). We present all aspects necessary to effect this acceleration within the framework of ACE including the necessary translation operators, and appropriately modifying the hierarchical computational algorithm. We also present several results that validate the efficacy of this method with respect to both error convergence and cost scaling, and derive error bounds for one exemplary potential.

preprint2011arXiv

Generalized Method of Moments: A novel framework for analyzing acoustic scattering from complex objects using a locally smooth surface parametrization and adaptive basis spaces

The analysis of scattering from complex objects using surface integral equations is a challenging problem. Its resolution has wide ranging applications- from crack propagation to diagnostic medicine. The two ingredients of any integral equation methodology is the representation of the domain and the design of approximation spaces to represent physical quantities on the domain. The order of convergence depends on both the surface and geometry representation. For instance, most surface models are restricted to piecewise at or second order tessellations. Similarly, the most commonly known basis spaces for acoustics are piecewise constant functions. What is desirable is a framework that permits adaptivity (of size and order) in both geometry and function representations. Unlike volumetric, differential equation solvers, such as the finite element method, developing an hpadaptive framework for surface integral equations is very difficult. This papers proposes a resolution to this problem by developing a novel framework that relies on reconstruction of the surface using locally smooth parameterizations, and defining partition of unity functions and higher order basis spaces on overlapping domains. This permits easy refinement of both the geometry and function representation. This capabilities of the proposed framework are shown via a number of numerical examples

preprint2011arXiv

The Rapid Analysis of Scattering from Periodic Dielectric Structures Using Accelerated Cartesian Expansions (ACE)

The analysis of fields in periodic dielectric structures arise in numerous applications of recent interest, ranging from photonic bandgap (PBG) structures and plasmonically active nanostructures to metamaterials. To achieve an accurate representation of the fields in these structures using numerical methods, dense spatial discretization is required. This, in turn, affects the cost of analysis, particularly for integral equation based methods, for which traditional iterative methods require O(N^2) operations, N being the number of spatial degrees of freedom. In this paper, we introduce a method for the rapid solution of volumetric electric field integral equations used in the analysis of doubly periodic dielectric structures. The crux of our method is the ACE algorithm, which is used to evaluate the requisite potentials in O(N) cost. Results are provided that corroborate our claims of acceleration without compromising accuracy, as well as the application of our method to a number of compelling photonics applications.