Source author record

Abhishek Saha

Abhishek Saha 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

21works
9topics
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

21 published item(s)

preprint2021arXiv

Neural ODE to model and prognose thermoacoustic instability

In reacting flow systems, thermoacoustic instability characterized by high amplitude pressure fluctuations, is driven by a positive coupling between the unsteady heat release rate and the acoustic field of the combustor. When the underlying flow is turbulent, as a control parameter of the system is varied and the system approach thermoacoustic instability, the acoustic pressure oscillations synchronize with heat release rate oscillations. Consequently, during the onset of thermoacoustic instability in turbulent combustors, the system dynamics transition from chaotic oscillations to periodic oscillations via a state of intermittency. Thermoacoustic systems are traditionally modeled by coupling the model for the unsteady heat source and the acoustic subsystem, each estimated independently. The response of the unsteady heat source, the flame, to acoustic fluctuations are characterized by introducing external unsteady forcing. This necessitates a powerful excitation module to obtain the nonlinear response of the flame to acoustic perturbations. Instead of characterizing individual subsystems, we introduce a neural ordinary differential equation (neural ODE) framework to model the thermoacoustic system as a whole. The neural ODE model for the thermoacoustic system uses time series of the heat release rate and the pressure fluctuations, measured simultaneously without introducing any external perturbations, to model their coupled interaction. Further, we use the parameters of neural ODE to define an anomaly measure that represents the proximity of system dynamics to limit cycle oscillations and thus provide an early warning signal for the onset of thermoacoustic instability.

preprint2020arXiv

Modeling ambient temperature and relative humidity sensitivity of respiratory droplets and their role in Covid-19 outbreaks

One of the many unresolved questions that revolves around the Covid-19 pandemic is whether local outbreaks can depend on ambient conditions like temperature and relative humidity. In this paper, we develop a model that tries to explain and describe the temperature and relative humidity sensitivity of respiratory droplets and their possible connection in determining viral outbreaks. The model has two parts. First, we model the growth rate of the infected population based on a reaction mechanism - the final equations of which are similar to the well-known SIR model. The advantage of modeling the pandemic using the reaction mechanism is that the rate constants have sound physical interpretation. The infection rate constant is derived using collision rate theory and shown to be a function of the respiratory droplet lifetime. In the second part, we have emulated the respiratory droplets responsible for disease transmission as salt solution droplets and computed their evaporation time accounting for droplet cooling, heat and mass transfer and finally crystallization of the salt. The model output favourably compares with the experimentally obtained evaporation characteristics of levitated droplets of pure water and salt solution, respectively, ensuring fidelity of the model. Droplet evaporation/desiccation time is indeed dependent on ambient temperature and relative humidity, considered at both outdoor and indoor conditions. Since the droplet evaporation time determines the infection rate constant, ambient temperature and relative humidity are shown to impact the outbreak growth rates.

preprint2020arXiv

Sup-norms of eigenfunctions in the level aspect for compact arithmetic surfaces, II: newforms and subconvexity

We improve upon the local bound in the depth aspect for sup-norms of newforms on $D^\times$ where $D$ is an indefinite quaternion division algebra over $\mathbb{Q}$. Our sup-norm bound implies a depth-aspect subconvexity bound for $L(1/2, f \times θ_χ)$, where $f$ is a (varying) newform on $D^\times$ of level $p^n$, and $θ_χ$ is an (essentially fixed) automorphic form on $\mathrm{GL}_2$ obtained as the theta lift of a Hecke character $χ$ on a quadratic field. For the proof, we augment the amplification method with a novel filtration argument and a recent counting result proved by the second-named author to reduce to showing strong quantitative decay of matrix coefficients of local newvectors along compact subsets, which we establish via $p$-adic stationary phase analysis. Furthermore, we prove a general upper bound in the level aspect for sup-norms of automorphic forms belonging to \emph{any} family whose associated matrix coefficients have such a decay property.

preprint2020arXiv

Synchronization framework for modeling transition to thermoacoustic instability in laminar combustors

We, herein, present a new model based on the framework of synchronization to describe a thermoacoustic system and capture the multiple bifurcations that such a system undergoes. Instead of applying flame describing function to depict the unsteady heat release rate as the flame's response to acoustic perturbation, the new model considers the acoustic field and the unsteady heat release rate as a pair of nonlinearly coupled damped oscillators. By varying the coupling strength, multiple dynamical behaviors, including limit cycle oscillation, quasi-periodic oscillation, strange nonchaos, and chaos can be captured. Furthermore, the model was able to qualitatively replicate the different behaviors of a laminar thermoacoustic system observed in experiments by Kabiraj et al.~[Chaos 22, 023129 (2012)]. By analyzing the temporal variation of the phase difference between heat release rate oscillations and pressure oscillations under different dynamical states, we show that the characteristics of the dynamical states depend on the nature of synchronization between the two signals, which is consistent with previous experimental findings.

preprint2018arXiv

Some analytic aspects of automorphic forms on GL(2) of minimal type

Let $π$ be a cuspidal automorphic representation of $PGL_2(\mathbb{A}_\mathbb{Q})$ of arithmetic conductor $C$ and archimedean parameter $T$, and let $ϕ$ be an $L^2$-normalized automorphic form in the space of $π$. The sup-norm problem asks for bounds on $\| ϕ\|_\infty$ in terms of $C$ and $T$. The quantum unique ergodicity (QUE) problem concerns the limiting behavior of the $L^2$-mass $|ϕ|^2 (g) \, d g$ of $ϕ$. All previous work on these problems in the conductor-aspect has focused on the case that $ϕ$ is a newform. In this work, we study these problems for a class of automorphic forms that are not newforms. Precisely, we assume that for each prime divisor $p$ of $C$, the local component $π_p$ is supercuspidal (and satisfies some additional technical hypotheses), and consider automorphic forms $ϕ$ for which the local components $ϕ_p \in π_p$ are "minimal" vectors. Such vectors may be understood as non-archimedean analogues of lowest weight vectors in holomorphic discrete series representations of $PGL_2(\mathbb{R})$. For automorphic forms as above, we prove a sup-norm bound that is sharper than what is known in the newform case. In particular, if $π_\infty$ is a holomorphic discrete series of lowest weight $k$, we obtain the optimal bound $C^{1/8 -ε} k^{1/4 - ε} \ll_ε |ϕ|_\infty \ll_ε C^{1/8 + ε} k^{1/4+ε}$. We prove also that these forms give analytic test vectors for the QUE period, thereby demonstrating the equivalence between the strong QUE and the subconvexity problems for this class of vectors. This finding contrasts the known failure of this equivalence for newforms of powerful level.

preprint2016arXiv

A note on the growth of nearly holomorphic vector-valued Siegel modular forms

Let $F$ be a nearly holomorphic vector-valued Siegel modular form of weight $ρ$ with respect to some congruence subgroup of $\mathrm{Sp}_{2n}(\mathbb Q)$. In this note, we prove that the function on $\mathrm{Sp}_{2n}(\mathbb R)$ obtained by lifting $F$ has the moderate growth (or "slowly increasing") property. This is a consequence of the following bound that we prove: $\|ρ(Y^{1/2})F(Z) \| \ll \prod_{i=1}^n (μ_i(Y)^{λ_1/2} + μ_i(Y)^{-λ_1/2})$ where $ λ_1 \ge \ldots \ge λ_n$ is the highest weight of $ρ$ and $μ_i(Y)$ are the eigenvalues of the matrix $Y$.

preprint2015arXiv

Bayesian inference for Gaussian graphical models beyond decomposable graphs

Bayesian inference for graphical models has received much attention in the literature in recent years. It is well known that when the graph G is decomposable, Bayesian inference is significantly more tractable than in the general non-decomposable setting. Penalized likelihood inference on the other hand has made tremendous gains in the past few years in terms of scalability and tractability. Bayesian inference, however, has not had the same level of success, though a scalable Bayesian approach has its respective strengths, especially in terms of quantifying uncertainty. To address this gap, we propose a scalable and flexible novel Bayesian approach for estimation and model selection in Gaussian undirected graphical models. We first develop a class of generalized G-Wishart distributions with multiple shape parameters for an arbitrary underlying graph. This class contains the G-Wishart distribution as a special case. We then introduce the class of Generalized Bartlett (GB) graphs, and derive an efficient Gibbs sampling algorithm to obtain posterior draws from generalized G-Wishart distributions corresponding to a GB graph. The class of Generalized Bartlett graphs contains the class of decomposable graphs as a special case, but is substantially larger than the class of decomposable graphs. We proceed to derive theoretical properties of the proposed Gibbs sampler. We then demonstrate that the proposed Gibbs sampler is scalable to significantly higher dimensional problems as compared to using an accept-reject or a Metropolis-Hasting algorithm. Finally, we show the efficacy of the proposed approach on simulated and real data.

preprint2015arXiv

Large values of newforms on GL(2) with highly ramified central character

We give a lower bound for the sup-norm of an $L^2$-normalized newform in an irreducible, unitary, cuspidal representation $π$ of $GL_2$ over a number field. When the central character of $π$ is sufficiently ramified, this bound improves upon the trivial bound by a positive power of $N$ where $N$ is the norm of the conductor of $π$. This generalizes a result of Templier, who dealt with the special case when the conductor of the central character equals the conductor of the representation. We also make a conjecture about the true size of the sup-norm in the $N$-aspect that takes into account this central character phenomenon. Our results depend upon some explicit formulas and bounds for the Whittaker newvector over a non-archimedean local field, which may be of independent interest.

preprint2015arXiv

Representations of SL_2(R) and nearly holomorphic modular forms

In this semi-expository note, we give a new proof of a structure theorem due to Shimura for nearly holomorphic modular forms on the complex upper half plane. Roughly speaking, the theorem says that the space of all nearly holomorphic modular forms is the direct sum of the subspaces obtained by applying appropriate weight-raising operators on the spaces of holomorphic modular forms and on the one-dimensional space spanned by the weight 2 nearly holomorphic Eisenstein series. While Shimura's proof was classical, ours is representation-theoretic. We deduce the structure theorem from a decomposition for the space of n-finite automorphic forms on SL_2(R). To prove this decomposition, we use the mechanism of category O and a careful analysis of the various possible indecomposable submodules. It is possible to achieve the same end by more direct methods, but we prefer this approach as it generalizes to other groups. This note may be viewed as the toy case of our paper ["Lowest weight modules of Sp_4(R) and nearly holomorphic Siegel modular forms"], where we prove an analogous structure theorem for vector-valued nearly holomorphic Siegel modular forms of degree two.

preprint2014arXiv

Local and global Maass relations (expanded version)

We characterize the irreducible, admissible, spherical representations of GSp(4,F) (where F is a p-adic field) that occur in certain CAP representations in terms of relations satisfied by their spherical vector in a special Bessel model. These local relations are analogous to the Maass relations satisfied by the Fourier coefficients of Siegel modular forms of degree 2 in the image of the Saito-Kurokawa lifting. We show how the classical Maass relations can be deduced from the local relations in a representation theoretic way, without recourse to the construction of Saito-Kurokawa lifts in terms of Fourier coefficients of half-integral weight modular forms or Jacobi forms. As an additional application of our methods, we give a new characterization of Saito-Kurokawa lifts involving a certain average of Fourier coefficients.

preprint2013arXiv

Bounds for Rankin--Selberg integrals and quantum unique ergodicity for powerful levels

Let f be a classical holomorphic newform of level q and even weight k. We show that the pushforward to the full level modular curve of the mass of f equidistributes as qk -> infinity. This generalizes known results in the case that q is squarefree. We obtain a power savings in the rate of equidistribution as q becomes sufficiently "powerful" (far away from being squarefree), and in particular in the "depth aspect" as q traverses the powers of a fixed prime. We compare the difficulty of such equidistribution problems to that of corresponding subconvexity problems by deriving explicit extensions of Watson's formula to certain triple product integrals involving forms of non-squarefree level. By a theorem of Ichino and a lemma of Michel--Venkatesh, this amounts to a detailed study of Rankin--Selberg integrals int|f|^2 E attached to newforms f of arbitrary level and Eisenstein series E of full level. We find that the local factors of such integrals participate in many amusing analogies with global L-functions. For instance, we observe that the mass equidistribution conjecture with a power savings in the depth aspect is equivalent to the union of a global subconvexity bound and what we call a "local subconvexity bound"; a consequence of our local calculations is what we call a "local Lindelof hypothesis".

preprint2013arXiv

Transfer of Siegel cusp forms of degree 2

Let $π$ be the automorphic representation of $\GSp_4(\A)$ generated by a full level cuspidal Siegel eigenform that is not a Saito-Kurokawa lift, and $τ$ be an arbitrary cuspidal, automorphic representation of $\GL_2(\A)$. Using Furusawa's integral representation for $\GSp_4\times\GL_2$ combined with a pullback formula involving the unitary group $\GU(3,3)$, we prove that the $L$-functions $L(s,π\timesτ)$ are "nice". The converse theorem of Cogdell and Piatetski-Shapiro then implies that such representations $π$ have a functorial lifting to a cuspidal representation of $\GL_4(\A)$. Combined with the exterior-square lifting of Kim, this also leads to a functorial lifting of $π$ to a cuspidal representation of $\GL_5(\A)$. As an application, we obtain analytic properties of various $L$-functions related to full level Siegel cusp forms. We also obtain special value results for $\GSp_4\times\GL_1$ and $\GSp_4\times\GL_2$.

preprint2013arXiv

Yoshida lifts and simultaneous non-vanishing of dihedral twists of modular L-functions

Given elliptic modular forms f and g satisfying certain conditions on their weights and levels, we prove (a quantitative version of the statement) that there exist infinitely many imaginary quadratic fields K and characters chi of the ideal class group Cl_K such that L(1/2, f_K \times chi) \neq 0 and L(1/2, g_K \times chi) \neq 0. The proof is based on a non-vanishing result for Fourier coefficients of Siegel modular forms combined with the theory of Yoshida liftings.

preprint2012arXiv

Determination of modular forms by fundamental Fourier coefficients

This article is a research exposition based on the author's talk at the International Colloquium on Automorphic Representations and L-Functions, 2012, held at TIFR, Mumbai. We consider some special cases of the following question: when is a natural subset of the Fourier coefficients sufficient to uniquely determine a "modular form"? Two kinds of modular forms are considered in this article: a) classical modular forms of half-integral weight, and b) Siegel modular forms of genus 2 and integral weight. These two apparently different scenarios turn out to be closely related. Our results were motivated by, and have several interesting connections to, automorphic L-functions and Bessel models.

preprint2012arXiv

Siegel cusp forms of degree 2 are determined by their fundamental Fourier coefficients

We prove that a Siegel cusp form of degree 2 for the full modular group is determined by its set of Fourier coefficients a(S) with 4 det(S) ranging over odd squarefree integers. As a key step to our result, we also prove that a classical cusp form of half-integral weight and level 4N, with N odd and squarefree, is determined by its set of Fourier coefficients a(d) with d ranging over odd squarefree integers, a result that was previously known only for Hecke eigenforms.

preprint2010arXiv

Prime density results for Hecke eigenvalues of a Siegel cusp form

Let F in S_k(Sp(2g, Z)) be a cuspidal Siegel eigenform of genus g with normalized Hecke eigenvalues mu_F(n). Suppose that the associated automorphic representation pi_F is locally tempered everywhere. For each c>0 we consider the set of primes p for which |mu_F(p)| >= c and we provide an explicit upper bound on the density of this set. In the case g=2, we also provide an explicit upper bound on the density of the set of primes p for which mu_F(p) >= c.

preprint2010arXiv

Structure formation by nanosilica particles suspended in levitated droplet

Vaporization of liquid droplets containing particles has been studied extensively for its applications in combustion, thermal coating, ink-jet printing, spray cooling, drug delivery, and surface patterning. Droplets containing solid particles show a preferential solute-migration during drying process. Recently we carried out experiments with vaporizing droplet suspended in an acoustic levitator. In this work, we present detailed study of a laser irradiated droplet containing nanosilica particles. Infrared and High speed imaging of the heating process for different concentrations of nanosilica revealed an interesting solute migration pattern. Further investigation with Particle Image Velocimetry shows presence of strong recirculation within the levitated droplet. It also reveals that with increasing viscosity of the liquid the strength of this recirculation decreases. Due to the droplets rotation about the levitator axis, a centrifugal force also dominated the flow field within the droplet. High speed imaging shows that initial stage of heating is dominated by sharper diameter reduction. At the end of this stage the droplet takes a bowl shape. Next stage, which is termed as structure formation shows different characteristics for higher and lower concentrations. For nanosilica concentration of less than 1.3%, the bowl shape remains unaltered with further heating. On the other hand, for concentration greater than 1.9%, due to stronger centrifugal effect than recirculation, the particles migrate towards the equator plane creating horizontal ring structure. The horizontal ring very soon starts oscillating about the levitator axis. This creates an imbalance of forces resulting in reorientation of the ring. Thus a vertical ring is created.