BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

B. Shanker

B. Shanker contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.comp-phmath.NAComputational Engineering, Finance, and Scienceeess.SPmath-phmath.MPNumerical Analysisphysics.plasm-ph

Trust snapshot

Quick read

Trust 21 - Emerging
8works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

8 published item(s)

preprint2022arXiv

An Envelope Tracking Approach for Particle in Cell Simulations

The state of the art in electromagnetic Finite Element Particle-in-Cell (EM-FEMPIC) has advanced significantly in the last few years. These have included understanding function spaces that must be used to represent sources and fields consistently, and how currents should be evolved in space and time. In concert, these achieve satisfaction of Gauss laws. All of these, were restricted to conditionally stable explicit time stepping. More recently, there has been advances to the state of art. It is now possible to use a implicit EM FEMPIC method while satisfying Gauss laws to machine precision. This enables choosing time step sizes dictated by physics as opposed to geometry. In this paper, we take this a step further. For devices characterized by a narrowband high frequency response, choosing a time-step size based on the highest frequency of interest is considerably expensive. In this paper, we use methods derived from envelope tracking to construct an EM-FEMPIC method that analytically provides for the high-frequency oscillations of the system, allowing for analysis at considerable coarser time-step sizes even in the presence of non-linear effects from active media such as plasmas. Consequentially, we demonstrate how the pointwise metric used for measuring satisfaction of Gauss Laws breaks down when prescribing analytical fast fields and provide a thorough analysis of how charge conservation can be measured. Through a number of examples, we demonstrate that the proposed approach retains the accuracy the regular scheme while requiring far fewer time steps.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Domain Decomposition Framework for Maxwell Finite Element Solvers and Application to PIC

The most popular methods for self-consistent simulation of fields interacting with charged species is using finite difference time domain (FDTD) methods together with Newton's laws of motion to evolve locations and velocities of particles. Despite their popularity, the limitation of FDTD particle in cell (EM-FDTDPIC) methods are well known. To address these, there has been significant interest over the past decade in exploring alternatives. In the past few years, the advances in electromagnetic finite element methods for particle in cell (EM-FEMPIC) has advanced by leaps and bounds. The mathematics necessary for implicit FEM methods that are unconditionally stable and charge conserving are now well understood. Some of these advances are more recent. The next bottleneck necessary to make EM-FEMPIC competitive with FDTD based scheme is overcoming computational cost. Our approach to resolving this challenge is develop two different finite element tearing and integration approaches, and using these to create domain decomposition schemes for EM-FEMPIC. Details of the proposed methodology are presented as well as a number of results that demonstrates charge conservation as well as amelioration of costs for a number of problems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Port Parameter Extraction Based Self Consistent Coupled EM-Circuit FEM Solvers

Self consistent solution to electromagnetic (EM)-circuit systems is of significant interest for a number of applications. This has resulted in exhaustive research on means to couple them. In time domain, this typically involves a tight integration (or coupling) with field and non-linear circuit solvers. This is in stark contrast to coupled analysis of linear/weakly non-linear circuits and EM systems in frequency domain. Here, one typically extracts equivalent port parameters that are then fed into the circuit solver. Such an approach has several advantages; (a) the number of ports is typically smaller than the number of degrees of freedom, resulting in cost savings; (b) is circuit agnostic; (c) can be integrated with a variety of device models. Port extraction is tantamount to obtaining impulse response of the linear EM system. In time domain, the deconvolution required to effect this is unstable. Recently, a novel approach was developed for time domain integral equations to overcome this bottleneck. We extend this approach to time domain finite element method, and demonstrate its utility via a number of examples; significantly, we demonstrate that self consistent solutions obtained using either a fully coupled or port extraction is identical to the desired precision for non-linear circuit systems. This is shown within a nodal network. We also demonstrate integration of port extracted data directly with drift diffusion equation to model device physics.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Laplace-Beltrami based Multi-Resolution Shape Reconstruction on Subdivision Surfaces

The eigenfunctions of the Laplace-Beltrami operator have widespread applications in a number of disciplines of engineering, computer vision/graphics, machine learning, etc. These eigenfunctions or manifold harmonics, provide the means to smoothly interpolate data on a manifold. They are highly effective, specifically as it relates to geometry representation and editing; manifold harmonics form a natural basis for multi-resolution representation (and editing) of complex surfaces and functioned defined therein. In this paper, we seek to develop the framework to exploit the benefits of manifold harmonics for shape reconstruction. To this end, we develop a highly compressible, multi-resolution shape reconstruction scheme using manifold harmonics. The method relies on subdivision basis sets to construct both boundary element isogeometric methods for analysis and surface finite elements to construct manifold harmonics. We pair this technique with the volumetric source reconstruction method to determine an initial starting point. Examples presented highlight efficacy of the approach in the presence of noisy data, including significant reduction in the number of degrees of freedom for complex objects, the accuracy of reconstruction, and multi-resolution capabilities.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Quasi-Helmholtz Decomposition, Gauss' Laws and Charge Conservation for Finite Element Particle-in-Cell

Development of particle in cell methods using finite element based methods (FEMs) have been a topic of renewed interest; this has largely been driven by (a) the ability of finite element methods to better model geometry, (b) better understanding of function spaces that are necessary to represent all Maxwell quantities, and (c) more recently, the fundamental rubrics that should be obeyed in space and time so as to satisfy Gauss' laws and the equation of continuity. In that vein, methods have been developed recently that satisfy these equations and are agnostic to time stepping methods. While is development is indeed a significant advance, it should be noted that implicit FEM transient solvers support an underlying null space that corresponds to a gradient of a scalar potential $\nabla Φ(\textbf{r})$ (or $t \nabla Φ(\textbf{r})$ in the case of wave equation solvers). While explicit schemes do not suffer from this drawback, they are only conditionally stable, time step sizes are mesh dependent, and very small. A way to overcome this bottleneck, and indeed, satisfy all four Maxwell's equation is to use a quasi-Helmholtz formulation on a tessellation. In the re-formulation presented, we strictly satisfy the equation of continuity and Gauss' laws for both the electric and magnetic flux densities. Results demonstrating the efficacy of this scheme will be presented.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Stable Higher Order Space-Time Galerkin Scheme for Time Domain Integral Equations

Stability of time domain integral equation (TDIE) solvers has remained an elusive goal for many years. Advancement of this research has largely progressed on four fronts: (1) Exact integration, (2) Lubich quadrature, (3) smooth temporal basis functions, and (4) Space-time separation of convolutions with the retarded potential. The latter method was explored in [Pray et al. IEEE TAP 2012]. This method's efficacy in stabilizing solutions to the time domain electric field integral equation (TD-EFIE) was demonstrated on first order surface descriptions (flat elements) in tandem with 0th order functions as the temporal basis. In this work, we develop the methodology necessary to extend to higher order surface descriptions as well as to enable its use with higher order temporal basis functions. These higher order temporal basis functions are used in a Galerkin framework. A number of results that demonstrate convergence, stability, and applicability are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On late-time stability of time domain integral equations for electromagnetics

The problem of late time instability in time domain integral equations for electromagnetics is longstanding. While several techniques have been suggested for addressing this problem, they either require impractically high degrees of freedom in the basis function or an analytical computation of matrix elements. The authors recently proposed a method that demonstrates stability without requiring either of these. The paper, however, does not present theoretical foundations for the choice of, or a rigorous error bounds on the approximation of, the bilinear form. This paper complements the authors' previous work by presenting a construction of the bilinear form based on the minimization of the energy in the system and a proof for the bounds on the approximation. We present results on the bounds developed and few sample scattering results that demonstrate the stability of the proposed scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Stability Properties of the Time Domain Electric Field Integral Equation Using a Separable Approximation for the Convolution with the Retarded Potential

The state of art of time domain integral equation (TDIE) solvers has grown by leaps and bounds over the past decade. During this time, advances have been made in (i) the development of accelerators that can be retrofitted with these solvers and (ii) understanding the stability properties of the electric field integral equation. As is well known, time domain electric field integral equation solvers have been notoriously difficult to stabilize. Research into methods for understanding and prescribing remedies have been on the uptick. The most recent of these efforts are (i) Lubich quadrature and (ii) exact integration. In this paper, we re-examine the solution to this equation using (i) the undifferentiated form of the TD-EFIE and (ii) a separable approximation to the spatio-temporal convolution. The proposed scheme can be constructed such that the spatial integrand over the source and observer domains is smooth and integrable. As several numerical results will demonstrate, the proposed scheme yields stable results for long simulation times and a variety of targets, both of which have proven extremely challenging in the past.

0 citations0 reviewsNo code linkedOpen access