Optimal control of information epidemics modeled as Maki Thompson rumors

We model the spread of information in a homogeneously mixed population using the Maki Thompson rumor model. We formulate an optimal control problem, from the perspective of single campaigner, to maximize the spread of information when the campaign budget is fixed. Control signals, such as advertising in the mass media, attempt to convert ignorants and stiflers into spreaders. We show the existence of a solution to the optimal control problem when the campaigning incurs non-linear costs under the isoperimetric budget constraint. The solution employs Pontryagin's Minimum Principle and a modified version of forward backward sweep technique for numerical computation to accommodate the isoperimetric budget constraint. The techniques developed in this paper are general and can be applied to similar optimal control problems in other areas. We have allowed the spreading rate of the information epidemic to vary over the campaign duration to model practical situations when the interest level of the population in the subject of the campaign changes with time. The shape of the optimal control signal is studied for different model parameters and spreading rate profiles. We have also studied the variation of the optimal campaigning costs with respect to various model parameters. Results indicate that, for some model parameters, significant improvements can be achieved by the optimal strategy compared to the static control strategy. The static strategy respects the same budget constraint as the optimal strategy and has a constant value throughout the campaign horizon. This work finds application in election and social awareness campaigns, product advertising, movie promotion and crowdfunding campaigns. (C) 2014 Elsevier B.V. All rights reserved.

Frustration in Optical Lattices of Interacting Bosons

The realization of optical lattices of cold atoms has opened up the possibility of engineering interacting lattice systems of bosons and fermions, stimulating a frenzy of research over the last decade. More recently, experimental techniques have been developed to apply synthetic gauge fields to these optical lattices. As a result, it has become possible to study quantum Hall physics and the effects of frustration in lattices of cold atoms. In this article we describe the combined effect of frustration and interactions on the superfluidity of bosons. By focussing on a frustrated ladder of interacting bosons, we show that the effect of frustration is for ``chiral'' order to develop, which manifests itself as an alternating pattern of circulating supercurrents. Remarkably, this order persists even when superfluidity is lost and the system enters a Mott phase giving rise to a novel chiral Mott insulator. We describe the combined physics of frustration and interactions by studying a fully frustrated one dimensional model of interacting bosons. The model is studied using mean-field theory, a direct quantum simulation and a higher dimensional classical theory in order to offer a full description of the different quantum phases contained in it and transitions between the different phases. In addition, we provide physical descriptions of the chiral Mott insulator as a vortex-anitvortex super solid and indirect excitonic condensate in addition to obtaining a variational wavefunction for it. We also briefly describe the chiral Mott states arising in other microscopic models.

Shear-imposed falling film

The study of a film falling down an inclined plane is revisited in the presence of imposed shear stress. Earlier studies regarding this topic (Smith, J. Fluid Mech., vol. 217, 1990, pp. 469-485; Wei, Phys. Fluids, vol. 17, 2005a, 012103), developed on the basis of a low Reynolds number, are extended up to moderate values of the Reynolds number. The mechanism of the primary instability is provided under the framework of a two-wave structure, which is normally a combination of kinematic and dynamic waves. In general, the primary instability appears when the kinematic wave speed exceeds the speed of dynamic waves. An equality criterion between their speeds yields the neutral stability condition. Similarly, it is revealed that the nonlinear travelling wave solutions also depend on the kinematic and dynamic wave speeds, and an equality criterion between the speeds leads to an analytical expression for the speed of a family of travelling waves as a function of the Froude number. This new analytical result is compared with numerical prediction, and an excellent agreement is achieved. Direct numerical simulations of the low-dimensional model have been performed in order to analyse the spatiotemporal behaviour of nonlinear waves by applying a constant shear stress in the upstream and downstream directions. It is noticed that the presence of imposed shear stress in the upstream (downstream) direction makes the evolution of spatially growing waves weaker (stronger).

Secondary instability as a possible mechanism for clear-air turbulence: a case study

Many theories and mechanisms have been proposed to explain the phenomenon of clear-air turbulence (CAT), and some of them have been successful in predicting light, moderate and, in some cases, severe turbulence. It is only recently that skill in the forecasting of the severe form of CAT, which could lead to injuries to passengers and damage to aircraft, has improved. Recent observations and simulations suggest that some severe to extreme turbulence could be caused by horizontal vortex tubes resulting from secondary instabilities of regions of high shear in the atmosphere. We have conducted direct numerical simulations to understand the scale relationship between primary structures (larger-scale structures related to one of the causes mentioned above) and secondary structures (smaller-sized, shear structures of the size of aircraft). From shear layer simulations, we find that the ratio of sizes of primary and secondary vortices is of the right order to generate aircraft-scale vortex tubes from typical atmospheric shear layers. We have also conducted simulations with a mesoscale atmospheric model, to understand possible causes of turbulence experienced by a flight off the west coast of India. Our simulations show the occurrence of primary flow structures related to synoptic conditions around the time of the incident. The evidence presented for this mechanism also has implications for possible methods of detection and avoidance of severe CAT.

Transition from dissipative to conservative dynamics in equations of hydrodynamics

We show, by using direct numerical simulations and theory, how, by increasing the order of dissipativity (alpha) in equations of hydrodynamics, there is a transition from a dissipative to a conservative system. This remarkable result, already conjectured for the asymptotic case alpha -> infinity U. Frisch et al., Phys. Rev. Lett. 101, 144501 (2008)], is now shown to be true for any large, but finite, value of alpha greater than a crossover value alpha(crossover). We thus provide a self-consistent picture of how dissipative systems, under certain conditions, start behaving like conservative systems and hence elucidate the subtle connection between equilibrium statistical mechanics and out-of-equilibrium turbulent flows.

Riemann shock tube: 1D normal shocks in air, simulations and experiments

This paper presents numerical simulation of the evolution of one-dimensional normal shocks, their propagation, reflection and interaction in air using a single diaphragm Riemann shock tube and validate them using experimental results. Mathematical model is derived for one-dimensional compressible flow of viscous and conducting medium. Dimensionless form of the mathematical model is used to construct space-time finite element processes based on minimization of the space-time residual functional. The space-time local approximation functions for space-time p-version hierarchical finite elements are considered in higher order GRAPHICS] spaces that permit desired order of global differentiability of local approximations in space and time. The resulting algebraic systems from this approach yield unconditionally positive-definite coefficient matrices, hence ensure unique numerical solution. The evolution is computed for a space-time strip corresponding to a time increment Delta t and then time march to obtain the evolution up to any desired value of time. Numerical studies are designed using recently invented hand-driven shock tube (Reddy tube) parameters, high/low side density and pressure values, high- and low-pressure side shock tube lengths, so that numerically computed results can be compared with actual experimental measurements.

Life of flame particles embedded in premixed flames interacting with near isotropic turbulence

Flame particles are surface points that always remain embedded on, by comoving with a given iso-scalar surface within a flame. Tracking flame particles allow us to study the fate of propagating surface locations uniquely identified throughout their evolution with time. In this work, using Direct Numerical Simulations we study the finite lifetime of such flame particles residing on iso-temperature surfaces of statistically planar H-2-air flames interacting with near-isotropic turbulence. We find that individual flame particles as well as their ensemble, experience progressively increasing tangential straining rate (K-t) and increasing negative curvature (kappa) near the end of their lifetime to finally get annihilated. By studying two different turbulent flow conditions, flame particle tracking shows that such tendency of local flame surfaces to be strained and cusped towards pinch-off from the main surface is a rather generic feature, independent of initial conditions, locations and ambient turbulence intensity levels. The evolution of the alignments between the flame surface normals and the principal components of the local straining rates are also tracked. We find that the surface normals initially aligned with the most extensive principal strain rate components, rotate near the end of flame particles' lifetime to enable preferential alignment between the surface tangent and the most extensive principal strain rate component. This could explain the persistently increasing tangential strain rate, sharp negative curvature formation and eventual detachment. (C) 2014 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Three Dimensional Impact Angle Constrained Guidance Law using Sliding Mode Control

In this paper, three dimensional impact angle control guidance laws are proposed for stationary targets. Unlike the usual approach of decoupling the engagement dynamics into two mutually orthogonal 2-dimensional planes, the guidance laws are derived using the coupled dynamics. These guidance laws are designed using principles of conventional as well as nonsingular terminal sliding mode control theory. The guidance law based on nonsingular terminal sliding mode guarantees finite time convergence of interceptor to the desired impact angle. In order to derive the guidance laws, multi-dimension switching surfaces are used. The stability of the system, with selected switching surfaces, is demonstrated using Lyapunov stability theory. Numerical simulation results are presented to validate the proposed guidance law.

Detection of a quantum particle on a lattice under repeated projective measurements

We consider a quantum particle, moving on a lattice with a tight-binding Hamiltonian, which is subjected to measurements to detect its arrival at a particular chosen set of sites. The projective measurements are made at regular time intervals tau, and we consider the evolution of the wave function until the time a detection occurs. We study the probabilities of its first detection at some time and, conversely, the probability of it not being detected (i.e., surviving) up to that time. We propose a general perturbative approach for understanding the dynamics which maps the evolution operator, which consists of unitary transformations followed by projections, to one described by a non-Hermitian Hamiltonian. For some examples of a particle moving on one-and two-dimensional lattices with one or more detection sites, we use this approach to find exact expressions for the survival probability and find excellent agreement with direct numerical results. A mean-field model with hopping between all pairs of sites and detection at one site is solved exactly. For the one-and two-dimensional systems, the survival probability is shown to have a power-law decay with time, where the power depends on the initial position of the particle. Finally, we show an interesting and nontrivial connection between the dynamics of the particle in our model and the evolution of a particle under a non-Hermitian Hamiltonian with a large absorbing potential at some sites.

TURBULENCE-TRANSPORT-CHEMISTRY INTERACTION IN STATISTICALLY PLANAR PREMIXED FLAMES AND IGNITION KERNELS IN NEAR ISOTROPIC TURBULENCE

Turbulence-transport-chemistry interaction plays a crucial role on the flame surface geometry, local and global reactionrates, and therefore, on the propagation and extinction characteristics of intensely turbulent, premixed flames encountered in LPP gas-turbine combustors. The aim of the present work is to understand these interaction effects on the flame surface annihilation and extinction of lean premixed flames, interacting with near isotropic turbulence. As an example case, lean premixed H-2-air mixture is considered so as to enable inclusion of detailed chemistry effects in Direct Numerical Simulations (DNS). The work is carried out in two phases namely, statistically planar flames and ignition kernel, both interacting with near isotropic turbulence, using the recently proposed Flame Particle Tracking (FPT) technique. Flame particles are surface points residing and commoving with an iso-scalar surface within a premixed flame. Tracking flame particles allows us to study the evolution of propagating surface locations uniquely identified with time. In this work, using DNS and FPT we study the flame speed, reaction rate and transport histories of such flame particles residing on iso-scalar surfaces. An analytical expression for the local displacement flame speed (SO is derived, and the contribution of transport and chemistry on the displacement flame speed is identified. An examination of the results of the planar case leads to a conclusion that the cause of variation in S-d may be attributed to the effects of turbulent transport and heat release rate. In the second phase of this work, the sustenance of an ignition kernel is examined in light of the S-curve. A newly proposed Damkohler number accounting for local turbulent transport and reaction rates is found to explain either the sustenance or otherwise propagation of flame kernels in near isotropic turbulence.

A flame particle tracking analysis of turbulence-chemistry interaction in hydrogen-air premixed flames

Interactions of turbulence, molecular transport, and energy transport, coupled with chemistry play a crucial role in the evolution of flame surface geometry, propagation, annihilation, and local extinction/re-ignition characteristics of intensely turbulent premixed flames. This study seeks to understand how these interactions affect flame surface annihilation of lean hydrogen-air premixed turbulent flames. Direct numerical simulations (DNSs) are conducted at different parametric conditions with a detailed reaction mechanism and transport properties for hydrogen-air flames. Flame particle tracking (FPT) technique is used to follow specific flame surface segments. An analytical expression for the local displacement flame speed (S-d) of a temperature isosurface is considered, and the contributions of transport, chemistry, and kinematics on the displacement flame speed at different turbulence-flame interaction conditions are identified. In general, the displacement flame speed for the flame particles is found to increase with time for all conditions considered. This is because, eventually all flame surfaces and their resident flame particles approach annihilation by reactant island formation at the end of stretching and folding processes induced by turbulence. Statistics of principal curvature evolving in time, obtained using FPT, suggest that these islands are ellipsoidal on average enclosing fresh reactants. Further examinations show that the increase in S-d is caused by the increased negative curvature of the flame surface and eventual homogenization of temperature gradients as these reactant islands shrink due to flame propagation and turbulent mixing. Finally, the evolution of the normalized, averaged, displacement flame speed vs. stretch Karlovitz number are found to collapse on a narrow band, suggesting that a unified description of flame speed dependence on stretch rate may be possible in the Lagrangian description. (C) 2015 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Autoregressive cascades on random networks

A network cascade model that captures many real-life correlated node failures in large networks via load redistribution is studied. The considered model is well suited for networks where physical quantities are transmitted, e.g., studying large scale outages in electrical power grids, gridlocks in road networks, and connectivity breakdown in communication networks, etc. For this model, a phase transition is established, i.e., existence of critical thresholds above or below which a small number of node failures lead to a global cascade of network failures or not. Theoretical bounds are obtained for the phase transition on the critical capacity parameter that determines the threshold above and below which cascade appears or disappears, respectively, that are shown to closely follow numerical simulation results. (C) 2015 Elsevier B.V. All rights reserved.

Autoregressive cascades on random networks

A network cascade model that captures many real-life correlated node failures in large networks via load redistribution is studied. The considered model is well suited for networks where physical quantities are transmitted, e.g., studying large scale outages in electrical power grids, gridlocks in road networks, and connectivity breakdown in communication networks, etc. For this model, a phase transition is established, i.e., existence of critical thresholds above or below which a small number of node failures lead to a global cascade of network failures or not. Theoretical bounds are obtained for the phase transition on the critical capacity parameter that determines the threshold above and below which cascade appears or disappears, respectively, that are shown to closely follow numerical simulation results. (C) 2015 Elsevier B.V. All rights reserved.

Depletion of nonlinearity in magnetohydrodynamic turbulence: Insights from analysis and simulations

It is shown how suitably scaled, order-m moments, D-m(+/-), of the Elsasser vorticity fields in three-dimensional magnetohydrodynamics (MHD) can be used to identify three possible regimes for solutions of the MHD equations with magnetic Prandtl number P-M = 1. These vorticity fields are defined by omega(+/-) = curl z(+/-) = omega +/- j, where z(+/-) are Elsasser variables, and where omega and j are, respectively, the fluid vorticity and current density. This study follows recent developments in the study of three-dimensional Navier-Stokes fluid turbulence Gibbon et al., Nonlinearity 27, 2605 (2014)]. Our mathematical results are then compared with those from a variety of direct numerical simulations, which demonstrate that all solutions that have been investigated remain in only one of these regimes which has depleted nonlinearity. The exponents q(+/-) that characterize the inertial range power-law dependencies of the z(+/-) energy spectra, epsilon(+/-)(k), are then examined, and bounds are obtained. Comments are also made on (a) the generalization of our results to the case P-M not equal 1 and (b) the relation between D-m(+/-) and the order-m moments of gradients of magnetohydrodynamic fields, which are used to characterize intermittency in turbulent flows.

Traveling waves in trimer granular lattice I: Bifurcation structure of traveling waves in the unit-cell model

Present paper is the first one in the series devoted to the dynamics of traveling waves emerging in the uncompressed, tri-atomic granular crystals. This work is primarily concerned with the dynamics of one-dimensional periodic granular trimer (tri-atomic) chains in the state of acoustic vacuum. Each unit cell consists of three spherical particles of different masses subject to periodic boundary conditions. Hertzian interaction law governs the mutual interaction of these particles. Under the assumption of zero pre-compression, this interaction is modeled as purely nonlinear, which means the absence of linear force component. The dynamics of such chains is governed by the two system parameters that scale the mass ratios between the particles of the unit cell. Such a system supports two different classes of periodic solutions namely the traveling and standing waves. The primary objective of the present study is the numerical analysis of the bifurcation structure of these solutions with emphasis on the dynamics of traveling waves. In fact, understanding of the bifurcation structure of the traveling wave solutions emerging in the unit-cell granular trimer is rather important and can shed light on the more complex nonlinear wave phenomena emerging in semi-infinite trimer chains. (c) 2016 Elsevier B.V. All rights reserved.