Source author record

Manas Rachh

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

7works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

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

7 published item(s)

preprint2023arXiv

FMM-LU: A fast direct solver for multiscale boundary integral equations in three dimensions

We present a fast direct solver for boundary integral equations on complex surfaces in three dimensions using an extension of the recently introduced recursive strong skeletonization scheme. For problems that are not highly oscillatory, our algorithm computes an ${LU}$-like hierarchical factorization of the dense system matrix, permitting application of the inverse in $\mathcal O(n)$ time, where $n$ is the number of unknowns on the surface. The factorization itself also scales linearly with the system size, albeit with a somewhat larger constant. The scheme is built on a level-restricted adaptive octree data structure, and therefore it is compatible with highly nonuniform discretizations. Furthermore, the scheme is coupled with high-order accurate locally-corrected Nyström quadrature methods to integrate the singular and weakly-singular Green's functions used in the integral representations. Our method has immediate applications to a variety of problems in computational physics. We concentrate here on studying its performance in acoustic scattering (governed by the Helmholtz equation) at low to moderate frequencies, and provide rigorous justification for compression of submatrices via proxy surfaces.

preprint2021arXiv

Debye source representation for type-I superconductors, I

In this note, we analyze the classical magneto-static approach to the theory of type I superconductors, and a Debye source representation that can be used numerically to solve the resultant equations. We also prove that one of the fields, $\boldsymbol{J}^-$, found within the superconductor via the London equations, is the physical current in that the outgoing part of the magnetic field is given as the Biot-Savart integral of $\boldsymbol{J}^{-}$. Finally, we compute the static currents for moderate values of London penetration depth, $λ_L,$ for a sphere, a stellarator-like geometry and a two-holed torus.

preprint2020arXiv

On the discretization of Laplace's equation with Neumann boundary conditions on polygonal domains

In the present paper we describe a class of algorithms for the solution of Laplace's equation on polygonal domains with Neumann boundary conditions. It is well known that in such cases the solutions have singularities near the corners which poses a challenge for many existing methods. If the boundary data is smooth on each edge of the polygon, then in the vicinity of each corner the solution to the corresponding boundary integral equation has an expansion in terms of certain (analytically available) singular powers. Using the known behavior of the solution, universal discretizations have been constructed for the solution of the Dirichlet problem. However, the leading order behavior of solutions to the Neumann problem is $O(t^μ)$ for $μ\in (-1/2,0)$ depending on the angle at the corner (compared to $O(C+t^μ)$ with $μ>1/2$ for the Dirichlet problem); this presents a significant challenge in the design of universal discretizations. Our approach is based on using the discretization for the Dirichlet problem in order to compute a solution in the "weak sense" by solving an adjoint linear system; namely, it can be used to compute inner products with smooth functions accurately, but it cannot be interpolated. Furthermore we present a procedure to obtain accurate solutions arbitrarily close to the corner, by solving a sequence of small local subproblems in the vicinity of that corner. The results are illustrated with several numerical examples.

preprint2020arXiv

On the solution of Laplace's equation in the vicinity of triple-junctions

In this paper we characterize the behavior of solutions to systems of boundary integral equations associated with Laplace transmission problems in composite media consisting of regions with polygonal boundaries. In particular we consider triple junctions, i.e. points at which three distinct media meet. We show that, under suitable conditions, solutions to the boundary integral equations in the vicinity of a triple junction are well-approximated by linear combinations of functions of the form $t^β,$ where $t$ is the distance of the point from the junction and the powers $β$ depend only on the material properties of the media and the angles at which their boundaries meet. Moreover, we use this analysis to design efficient discretizations of boundary integral equations for Laplace transmission problems in regions with triple junctions and demonstrate the accuracy and efficiency of this algorithm with a number of examples.

preprint2020arXiv

Recovering missing data in coherent diffraction imaging

In coherent diffraction imaging (CDI) experiments, the intensity of the scattered wave impinging on an object is measured on an array of detectors. This signal can be interpreted as the square of the modulus of the Fourier transform of the unknown scattering density. A beam-stop obstructs the forward scattered wave and, hence, the modulus Fourier data from a neighborhood of k=0 cannot be measured. In this note, we describe a linear method for recovering this unmeasured modulus Fourier data from the measured values and an estimate of the support of the image's autocorrelation function without consideration of phase retrieval. We analyze the conditioning of this problem, which grows exponentially with the modulus of the maximum spatial frequency not measured, and the effects of noise.

preprint2017arXiv

On the location of maximal of solutions of Schrödinger's equation

We prove an inequality with applications to solutions of the Schrödinger equation. There is a universal constant $c>0$, such that if $Ω\subset \mathbb{R}^2$ is simply connected, $u:Ω\rightarrow \mathbb{R}$ vanishes on the boundary $\partial Ω$, and $|u|$ assumes a maximum in $x_0 \in Ω$, then $$ \inf_{y \in \partial Ω}{ \| x_0 - y\|} \geq c \left\| \frac{Δu}{u} \right\|^{-1/2}_{L^{\infty}(Ω)}.$$ It was conjectured by Pólya \& Szegő (and proven, independently, by Makai and Hayman) that a membrane vibrating at frequency $λ$ contains a disk of size $\sim λ^{-1/2}$. Our inequality implies a refined result: the point on the membrane that achieves the maximal amplitude is at distance $\sim λ^{-1/2}$ from the boundary. We also give an extension to higher dimensions (generalizing results of Lieb and Georgiev \& Mukherjee): if $u$ solves $-Δu = Vu$ on $Ω\subset \mathbb{R}^n$ with Dirichlet boundary conditions, then the ball $B$ with radius $\sim \|V\|_{L^{\infty}(Ω)}^{-1/2}$ centered at the point in which $|u|$ assumes a maximum is almost fully contained in $Ω$ in the sense that $|B \cap Ω| \geq 0.99 |B|.$

preprint2016arXiv

On the Diffusion Geometry of Graph Laplacians and Applications

We study directed, weighted graphs $G=(V,E)$ and consider the (not necessarily symmetric) averaging operator $$ (\mathcal{L}u)(i) = -\sum_{j \sim_{} i}{p_{ij} (u(j) - u(i))},$$ where $p_{ij}$ are normalized edge weights. Given a vertex $i \in V$, we define the diffusion distance to a set $B \subset V$ as the smallest number of steps $d_{B}(i) \in \mathbb{N}$ required for half of all random walks started in $i$ and moving randomly with respect to the weights $p_{ij}$ to visit $B$ within $d_{B}(i)$ steps. Our main result is that the eigenfunctions interact nicely with this notion of distance. In particular, if $u$ satisfies $\mathcal{L}u = λu$ on $V$ and $$ B = \left\{ i \in V: - \varepsilon \leq u(i) \leq \varepsilon \right\} \neq \emptyset,$$ then, for all $i \in V$, $$ d_{B}(i) \log{\left( \frac{1}{|1-λ|} \right) } \geq \log{\left( \frac{ |u(i)| }{\|u\|_{L^{\infty}}} \right)} - \log{\left(\frac{1}{2} + \varepsilon\right)}.$$ $d_B(i)$ is a remarkably good approximation of $|u|$ in the sense of having very high correlation. The result implies that the classical one-dimensional spectral embedding preserves particular aspects of geometry in the presence of clustered data. We also give a continuous variant of the result which has a connection to the hot spots conjecture.