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

Shravan Veerapaneni

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

ResearcherUnclaimed source record
math.NANumerical Analysiscond-mat.softMachine Learningphysics.comp-phphysics.flu-dynmath.OCquant-phBiological PhysicsComputational Engineering, Finance, and Sciencecond-mat.dis-nnhep-lathep-thphysics.chem-phq-fin.MF

Catalog footprint

What is connected

14works
15topics
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

14 published item(s)

preprint2026arXiv

Convergence analysis for the Barrett--Garcke--Nurnberg method of transport type for evolving curves

In this paper, we propose a Barrett-Garcke-Nurnberg (BGN) method for evolving geometries under general flows and present the corresponding convergence analysis. Specifically, we examine the scenario where a closed curve evolves according to a prescribed background velocity field. Unlike mean curvature flow and surface diffusion, where the evolution velocities inherently exhibit parabolicity, this case is dominated by transport which poses a significant difficulty in establishing convergence proofs. To address the challenges imposed by this transport-dominant nature, we derive several discrete energy estimates of the transport type on discretized polynomial surfaces within the framework of the projection error. The use of the projection error is indispensable as it provides crucial additional stability through its orthogonality structure. We prove that the proposed method converges sub-optimally in the L2 norm, and this is the first convergence proof for a fully discrete numerical method solving the evolution of curves driven by general flows.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A fast direct solver for integral equations on locally refined boundary discretizations and its application to multiphase flow simulations

In transient simulations of particulate Stokes flow, to accurately capture the interaction between the constituent particles and the confining wall, the discretization of the wall often needs to be locally refined in the region approached by the particles. Consequently, standard fast direct solvers lose their efficiency since the linear system changes at each time step. This manuscript presents a new computational approach that avoids this issue by pre-constructing a fast direct solver for the wall ahead of time, computing a low-rank factorization to capture the changes due to the refinement, and solving the problem on the refined discretization via a Woodbury formula. Numerical results illustrate the efficiency of the solver in accelerating particulate Stokes simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Continuous-variable optimization with neural network quantum states

Inspired by proposals for continuous-variable quantum approximate optimization (CV-QAOA), we investigate the utility of continuous-variable neural network quantum states (CV-NQS) for performing continuous optimization, focusing on the ground state optimization of the classical antiferromagnetic rotor model. Numerical experiments conducted using variational Monte Carlo with CV-NQS indicate that although the non-local algorithm succeeds in finding ground states competitive with the local gradient search methods, the proposal suffers from unfavorable scaling. A number of proposed extensions are put forward which may help alleviate the scaling difficulty.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Numerical and geometrical aspects of flow-based variational quantum Monte Carlo

This article aims to summarize recent and ongoing efforts to simulate continuous-variable quantum systems using flow-based variational quantum Monte Carlo techniques, focusing for pedagogical purposes on the example of bosons in the field amplitude (quadrature) basis. Particular emphasis is placed on the variational real- and imaginary-time evolution problems, carefully reviewing the stochastic estimation of the time-dependent variational principles and their relationship with information geometry. Some practical instructions are provided to guide the implementation of a PyTorch code. The review is intended to be accessible to researchers interested in machine learning and quantum information science.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantum-inspired variational algorithms for partial differential equations: Application to financial derivative pricing

Variational quantum Monte Carlo (VMC) combined with neural-network quantum states offers a novel angle of attack on the curse-of-dimensionality encountered in a particular class of partial differential equations (PDEs); namely, the real- and imaginary time-dependent Schrödinger equation. In this paper, we present a simple generalization of VMC applicable to arbitrary time-dependent PDEs, showcasing the technique in the multi-asset Black-Scholes PDE for pricing European options contingent on many correlated underlying assets.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Scalable neural quantum states architecture for quantum chemistry

Variational optimization of neural-network representations of quantum states has been successfully applied to solve interacting fermionic problems. Despite rapid developments, significant scalability challenges arise when considering molecules of large scale, which correspond to non-locally interacting quantum spin Hamiltonians consisting of sums of thousands or even millions of Pauli operators. In this work, we introduce scalable parallelization strategies to improve neural-network-based variational quantum Monte Carlo calculations for ab-initio quantum chemistry applications. We establish GPU-supported local energy parallelism to compute the optimization objective for Hamiltonians of potentially complex molecules. Using autoregressive sampling techniques, we demonstrate systematic improvement in wall-clock timings required to achieve CCSD baseline target energies. The performance is further enhanced by accommodating the structure of resultant spin Hamiltonians into the autoregressive sampling ordering. The algorithm achieves promising performance in comparison with the classical approximate methods and exhibits both running time and scalability advantages over existing neural-network based methods.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A scalable computational platform for particulate Stokes suspensions

We describe a computational framework for simulating suspensions of rigid particles in Newtonian Stokes flow. One central building block is a collision-resolution algorithm that overcomes the numerical constraints arising from particle collisions. This algorithm extends the well-known complementarity method for non-smooth multi-body dynamics to resolve collisions in dense rigid body suspensions. This approach formulates the collision resolution problem as a linear complementarity problem with geometric `non-overlapping' constraints imposed at each timestep. It is then reformulated as a constrained quadratic programming problem and the Barzilai-Borwein projected gradient descent method is applied for its solution. This framework is designed to be applicable for any convex particle shape, e.g., spheres and spherocylinders, and applicable to any Stokes mobility solver, including the Rotne-Prager-Yamakawa approximation, Stokesian Dynamics, and PDE solvers (e.g., boundary integral and immersed boundary methods). In particular, this method imposes Newton's Third Law and records the entire contact network. Further, we describe a fast, parallel, and spectrally-accurate boundary integral method tailored for spherical particles, capable of resolving lubrication effects. We show weak and strong parallel scalings up to $8\times 10^4$ particles with approximately $4\times 10^7$ degrees of freedom on $1792$ cores. We demonstrate the versatility of this framework with several examples, including sedimentation of particle clusters, and active matter systems composed of ensembles of particles driven to rotate.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An Active Learning Framework for Constructing High-fidelity Mobility Maps

A mobility map, which provides maximum achievable speed on a given terrain, is essential for path planning of autonomous ground vehicles in off-road settings. While physics-based simulations play a central role in creating next-generation, high-fidelity mobility maps, they are cumbersome and expensive. For instance, a typical simulation can take weeks to run on a supercomputer and each map requires thousands of such simulations. Recent work at the U.S. Army CCDC Ground Vehicle Systems Center has shown that trained machine learning classifiers can greatly improve the efficiency of this process. However, deciding which simulations to run in order to train the classifier efficiently is still an open problem. According to PAC learning theory, data that can be separated by a classifier is expected to require $\mathcal{O}(1/ε)$ randomly selected points (simulations) to train the classifier with error less than $ε$. In this paper, building on existing algorithms, we introduce an active learning paradigm that substantially reduces the number of simulations needed to train a machine learning classifier without sacrificing accuracy. Experimental results suggest that our sampling algorithm can train a neural network, with higher accuracy, using less than half the number of simulations when compared to random sampling.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Optimal slip velocities of micro-swimmers with arbitrary axisymmetric shapes

This article presents a computational approach for determining the optimal slip velocities on any given shape of an axisymmetric micro-swimmer suspended in a viscous fluid. The objective is to minimize the power loss to maintain a target swimming speed, or equivalently to maximize the efficiency of the micro-swimmer. Owing to the linearity of the Stokes equations governing the fluid motion, we show that this PDE-constrained optimization problem reduces to a simpler quadratic optimization problem, whose solution is found using a high-order accurate boundary integral method. We consider various families of shapes parameterized by the reduced volume and compute their swimming efficiency. {Among those, prolate spheroids were found to be the most efficient micro-swimmer shapes for a given reduced volume. We propose a simple shape-based scalar metric that can determine whether the optimal slip on a given shape makes it a pusher, a puller, or a neutral swimmer.}

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Simulating cilia-driven mixing and transport in complex geometries

Cilia and flagella are self-actuated microtubule-based structures that are present on many cell surfaces, ranging from the outer surface of single-cell organisms to the internal epithelial surfaces in larger animals. A fast and robust numerical method that simulates the coupled hydrodynamics of cilia and the constituent particles in the fluid such as rigid particles, drops or cells would be useful to not only understand several disease and developmental pathologies due to ciliary dysfunction but also to design microfluidic chips with ciliated cultures for some targeted functionality---e.g., maximizing fluid transport or particle mixing. In this paper, we develop a hybrid numerical method that employs a boundary integral method for handling the confining geometries and the constituent rigid particles and the method of regularized Stokeslets for handling the cilia. We provide several examples demonstrating the effects of confining geometries on cilia-generated fluid mixing as well as the cilia-particle hydrodynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Solution of Stokes flow in complex nonsmooth 2D geometries via a linear-scaling high-order adaptive integral equation scheme

We present a fast, high-order accurate and adaptive boundary integral scheme for solving the Stokes equations in complex---possibly nonsmooth---geometries in two dimensions. The key ingredient is a set of panel quadrature rules capable of evaluating weakly-singular, nearly-singular and hyper-singular integrals to high accuracy. Near-singular integral evaluation, in particular, is done using an extension of the scheme developed in J.~Helsing and R.~Ojala, {\it J. Comput. Phys.} {\bf 227} (2008) 2899--2921. The boundary of the given geometry is ``panelized'' automatically to achieve user-prescribed precision. We show that this adaptive panel refinement procedure works well in practice even in the case of complex geometries with large number of corners. In one example, for instance, a model 2D vascular network with 378 corners required less than 200K discretization points to obtain a 9-digit solution accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A unified integral equation scheme for doubly-periodic Laplace and Stokes boundary value problems in two dimensions

We present a spectrally-accurate scheme to turn a boundary integral formulation for an elliptic PDE on a single unit cell geometry into one for the fully periodic problem. Applications include computing the effective permeability of composite media (homogenization), and microfluidic chip design. Our basic idea is to exploit a small least squares solve to apply periodicity without ever handling periodic Green's functions. We exhibit fast solvers for the two-dimensional (2D) doubly-periodic Neumann Laplace problem (flow around insulators), and Stokes non-slip fluid flow problem, that for inclusions with smooth boundaries achieve 12-digit accuracy, and can handle thousands of inclusions per unit cell. We split the infinite sum over the lattice of images into a directly-summed "near" part plus a small number of auxiliary sources which represent the (smooth) remaining "far" contribution. Applying physical boundary conditions on the unit cell walls gives an expanded linear system, which, after a rank-1 or rank-3 correction and a Schur complement, leaves a well-conditioned square system which can be solved iteratively using fast multipole acceleration plus a low-rank term. We are rather explicit about the consistency and nullspaces of both the continuous and discretized problems. The scheme is simple (no lattice sums, Ewald methods, nor particle meshes are required), allows adaptivity, and is essentially dimension- and PDE-independent, so would generalize without fuss to 3D and to other non-oscillatory elliptic problems such as elastostatics. We incorporate recently developed spectral quadratures that accurately handle close-to-touching geometries. We include many numerical examples, and provide a software implementation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamics of a multicomponent vesicle in shear flow

We study the fully nonlinear, nonlocal dynamics of two-dimensional multicomponent vesicles in a shear flow with matched viscosity of the inner and outer fluids. Using a nonstiff, pseudo-spectral boundary integral method, we investigate dynamical patterns induced by inhomogeneous bending for a two phase system. Numerical results reveal that there exist novel phase-treading and tumbling mechanisms that cannot be observed for a homogeneous vesicle. In particular, unlike the well-known steady tank-treading dynamics characterized by a fixed inclination angle, here the phase-treading mechanism leads to unsteady periodic dynamics with an oscillatory inclination angle. When the average phase concentration is around 1/2, we observe tumbling dynamics even for very low shear rate, and the excess length required for tumbling is significantly smaller than the value for the single phase case. We summarize our results in phase diagrams in terms of the excess length, shear rate, and concentration of the soft phase. These findings go beyond the well known dynamical regimes of a homogeneous vesicle and highlight the level of complexity of vesicle dynamics in a fluid due to heterogeneous material properties.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A fast algorithm for simulating multiphase flows through periodic geometries of arbitrary shape

This paper presents a new boundary integral equation (BIE) method for simulating particulate and multiphase flows through periodic channels of arbitrary smooth shape in two dimensions. The authors consider a particular system---multiple vesicles suspended in a periodic channel of arbitrary shape---to describe the numerical method and test its performance. Rather than relying on the periodic Green's function as classical BIE methods do, the method combines the free-space Green's function with a small auxiliary basis, and imposes periodicity as an extra linear condition. As a result, we can exploit existing free-space solver libraries, quadratures, and fast algorithms, and handle a large number of vesicles in a geometrically complex channel. Spectral accuracy in space is achieved using the periodic trapezoid rule and product quadratures, while a first-order semi-implicit scheme evolves particles by treating the vesicle-channel interactions explicitly. New constraint-correction formulas are introduced that preserve reduced areas of vesicles, independent of the number of time steps taken. By using two types of fast algorithms, (i) the fast multipole method (FMM) for the computation of the vesicle-vesicle and the vesicle-channel hydrodynamic interaction, and (ii) a fast direct solver for the BIE on the fixed channel geometry, the computational cost is reduced to $O(N)$ per time step where $N$ is the spatial discretization size. Moreover, the direct solver inverts the wall BIE operator at $t = 0$, stores its compressed representation and applies it at every time step to evolve the vesicle positions, leading to dramatic cost savings compared to classical approaches. Numerical experiments illustrate that a simulation with $N=128, 000$ can be evolved in less than a minute per time step on a laptop.

0 citations0 reviewsNo code linkedOpen access