Source author record

C. P. Vyasarayani

C. P. Vyasarayani 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

5works
4topics
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

5 published item(s)

preprint2020arXiv

Complete dimensional collapse in the continuum limit of a delayed SEIQR network model with separable distributed infectivity

We take up a recently proposed compartmental SEIQR model with delays, ignore loss of immunity in the context of a fast pandemic, extend the model to a network structured on infectivity, and consider the continuum limit of the same with a simple separable interaction model for the infectivities $β$. Numerical simulations show that the evolving dynamics of the network is effectively captured by a single scalar function of time, regardless of the distribution of $β$ in the population. The continuum limit of the network model allows a simple derivation of the simpler model, which is a single scalar delay differential equation (DDE), wherein the variation in $β$ appears through an integral closely related to the moment generating function of $u=\sqrtβ$. If the first few moments of $u$ exist, the governing DDE can be expanded in a series that shows a direct correspondence with the original compartmental DDE with a single $β$. Even otherwise, the new scalar DDE can be solved using either numerical integration over $u$ at each time step, or with the analytical integral if available in some useful form. Our work provides a new academic example of complete dimensional collapse, ties up an underlying continuum model for a pandemic with a simpler-seeming compartmental model, and will hopefully lead to new analysis of continuum models for epidemics.

preprint2020arXiv

Fast generation of stability charts for time-delay systems using continuation of characteristic roots

Many dynamic processes involve time delays, thus their dynamics are governed by delay differential equations (DDEs). Studying the stability of dynamic systems is critical, but analyzing the stability of time-delay systems is challenging because DDEs are infinite-dimensional. We propose a new approach to quickly generate stability charts for DDEs using continuation of characteristic roots (CCR). In our CCR method, the roots of the characteristic equation of a DDE are written as implicit functions of the parameters of interest, and the continuation equations are derived in the form of ordinary differential equations (ODEs). Numerical continuation is then employed to determine the characteristic roots at all points in a parametric space; the stability of the original DDE can then be easily determined. A key advantage of the proposed method is that a system of linearly independent ODEs is solved rather than the typical strategy of solving a large eigenvalue problem at each grid point in the domain. Thus, the CCR method significantly reduces the computational effort required to determine the stability of DDEs. As we demonstrate with several examples, the CCR method generates highly accurate stability charts, and does so up to 10 times faster than the Galerkin approximation method.

preprint2020arXiv

Galerkin-Ivanov transformation for nonsmooth modeling of vibro-impacts in continuous structures

This work deals with the modeling of nonsmooth vibro-impact motion of a continuous structure against a rigid distributed obstacle. Galerkin's approach is used to approximate the solutions of the governing partial differential equations of the structure, which results in a system of ordinary differential equations (ODEs). When these ODEs are subjected to unilateral constraints and velocity jump conditions, one must use an event detection algorithm to calculate the time of impact accurately. Event detection in the presence of multiple simultaneous impacts is a computationally demanding task. Ivanov proposed a nonsmooth transformation for a vibro-impacting multi-degree-of-freedom system subjected to a single unilateral constraint. This transformation eliminates the unilateral constraints from the problem and, therefore, no event detection is required during numerical integration. Ivanov used his transformation to make analytical calculations for the stability and bifurcations of vibro-impacting motions; however, he did not explore its application for simulating distributed collisions in spatially continuous structures. We adopt Ivanov's transformation to deal with multiple unilateral constraints in spatially continuous structures. Also, imposing the velocity jump conditions exactly in the modal coordinates is nontrivial and challenging. Therefore, in this work we use a modal-physical transformation to convert the system from modal to physical coordinates on a spatially discretized grid. We then apply Ivanov's transformation on the physical system to simulate the vibro-impact motion of the structure. The developed method is demonstrated by modeling the distributed collision of a nonlinear string against a rigid distributed surface. For validation, we compare our results with the well-known penalty approach.

preprint2020arXiv

New approximations, and policy implications, from a delayed dynamic model of a fast pandemic

We study an SEIQR (Susceptible-Exposed-Infectious-Quarantined-Recovered) model for an infectious disease, with time delays for latency and an asymptomatic phase. For fast pandemics where nobody has prior immunity and everyone has immunity after recovery, the SEIQR model decouples into two nonlinear delay differential equations (DDEs) with five parameters. One parameter is set to unity by scaling time. The subcase of perfect quarantining and zero self-recovery before quarantine, with two free parameters, is examined first. The method of multiple scales yields a hyperbolic tangent solution; and a long-wave approximation yields a first order ordinary differential equation (ODE). With imperfect quarantining and nonzero self-recovery, the long-wave approximation is a second order ODE. These three approximations each capture the full outbreak, from infinitesimal initiation to final saturation. Low-dimensional dynamics in the DDEs is demonstrated using a six state non-delayed reduced order model obtained by Galerkin projection. Numerical solutions from the reduced order model match the DDE over a range of parameter choices and initial conditions. Finally, stability analysis and numerics show how correctly executed time-varying social distancing, within the present model, can cut the number of affected people by almost half. Alternatively, faster detection followed by near-certain quarantining can potentially be even more effective.

preprint2020arXiv

Performance limit for base-excited energy harvesting, and comparison with experiments

We consider the theoretical maximum extractable average power from an energy harvesting device attached to a vibrating table which provides a unidirectional displacement $A\sin(ωt)$. The total mass of moving components in the device is $m$. The device is assembled in a container of dimension $L$, limiting the displacements and deformations of components within. The masses in the device may be interconnected in arbitrary ways. The maximum extractable average power is bounded by $\frac{mLAω^3}π$, for motions in 1, 2, or 3 dimensions; with both rectilinear and rotary motions as special cases; and with either single or multiple degrees of freedom. The limiting displacement profile of the moving masses for extracting maximum power is discontinuous, and not physically realizable. But smooth approximations can be nearly as good: with $15$ terms in a Fourier approximation, the upper limit is $99$\% of the theoretical maximum. Purely sinusoidal solutions are limited to $\fracπ{4}$ times the theoretical maximum. For both single-degree-of-freedom linear resonant devices and nonresonant whirling devices where the energy extraction mimics a linear torsional damper, the maximum average power output is $\frac{mLAω^3}{4}$. Thirty-six experimental energy harvesting devices in the literature are found to extract power amounts ranging from $0.0036$\% to $29$\% of the theoretical maximum. Of these thirty-six, twenty achieve less than 2\% and three achieve more than 20\%. We suggest, as a figure of merit, that energy extraction above $\frac{0.2 mLAω^3}π$ may be considered excellent, and extraction above $\frac{0.3 mLAω^3}π$ may be considered challenging.