key: cord-0444715-jpi4aq37 authors: Samant, Omstavan; Alageshan, Jaya Kumar; Sharma, Sarveshwar; Kuley, Animesh title: Dynamic Mode Decomposition of inertial particle caustics in Taylor-Green flow date: 2021-02-09 journal: nan DOI: nan sha: 6441a9dcec6d1a4636ef3f7c6cf0daa9f25cecac doc_id: 444715 cord_uid: jpi4aq37 Inertial particles advected by a background flow can show complex structures. We consider inertial particles in a 2D Taylor-Green (TG) flow and characterize particle dynamics as a function of the particle's Stokes number using dynamic mode decomposition (DMD) method from particle image velocimetry (PIV) like-data. We observe the formation of caustic structures and analyze them using DMD to (a) determine the Stokes number of the particles, and (b) estimate the particle Stokes number composition. Our analysis in this idealized flow will provide useful insight to analyze inertial particles in more complex or turbulent flows. We propose that the DMD technique can be used to perform a similar analysis on an experimental system. experimentally setup using ion solutions in an array of magnets 29 . The TG flow is given by the vorticity field as ω(x, y) = ω 0 sin 2πx L sin 2πy L (1) and the corresponding velocity field is where x, y ∈ [0, L) and are periodic, u is the Eulerian velocity field such that ω = ∇ × u. We choose V 0 as the velocity scale and L as the length scale and write the system parameters in corresponding dimensionless form. We model the aerosol particles as small rigid spheres, which are effectively points, that have density different from the surrounding fluid. The equation of motion of the inertial particles in a background flow given by the simplified Maxey-Riley approximation 17 for small particles that are much denser than the fluid are where St is the Stokes number which captures the effect of particle inertia, x is the particle position and v is the particle velocity (see Appendix A for validity of the equations). The case when St → 0 the particles act as tracers that follow the velocity stream lines and the eq. 3 leads to v = u. We use RK4 numerical scheme to discretize and evolve the particle positions and velocities. Furthermore, we use periodic boundary conditions, such that the particles are reintroduced into the system when they exit the boundary. In our analysis we set the time step to ∆t = 0.01 (L/V 0 ). A typical feature of inertial particles in a background flow is that they are expelled from high vorticity regions. In a background flow with vortices, the inertial particles can form caustics 20 . So in the case of TG flow, the inertial particles tend to accumulate along the low vorticity regions that separate the vortices. We simulate the dynamics of inertial particles in a TG flow and observe spatial regions where the Lagrangian velocity field is multi-valued, which are referred to as caustics 30-32 (see Fig. 3(a) ). Figure. 1 shows the caustic structure formed by a mono-disperse inertial particle system, when we start the simulation with the particle positions initialized with uniform random distribution within the L × L box, and setting the initial velocities of the particles to zero. In the long time limit, the particles accumulate along a curve 33 . We observe that the caustics that form in the transient state are robust and stable to small perturbations (see Appendix A), whereas the steady state structures break up and lead to chaos 33 . Furthermore, the steady state behaviour strongly depends on the system size and the boundary conditions. We find that for a range of Stokes numbers the caustic structures preserve their shapes; and their sizes depend on St. In the following section we use dynamic mode decomposition (DMD) to extract features of the structures, in particular the caustic wavefront, and study its size dependence on the Stokes number. Furthermore, the sharpness of the caustic wavefront enables us to detect and extract their sizes even in presence of poly-disperse Stokes number systems. The caustics in Fig. 1 have a complex structure and in the presence of multiple Stokes number particles resolving these structures from a single snap-shot is hard. Therefore we employ the spatio-temporal data in the form of a video sequence that contains F frames of N × N pixel images and analyze them using the dynamic mode decomposition (DMD) method. DMD is a data analysis technique that has been used to extract coherent structures in fluid dynamic systems 34 , where it is able to extract different modes that are similar to normal modes in linear dynamical systems. The DMD is a data-driven technique introduced by Schmid as a numerical procedure for extracting dynamical features from flow data 35 . The DMD algorithm takes in a time-series data in the form of vectors { v 1 , v 2 , ... v T } and estimates a linear dynamical system that can generate a map where A is an N 2 × N 2 matrix and the eigenvectors of A form the DMD modes, with the corresponding eigenvalues. Finally, the eigenvectors are reshaped into N × N pixel image to obtain the modes. A Singular Value Decomposition (SVD) based algorithm for estimating the DMD modes is described in Appendix B. Let i stand for the iteration number such that the particles are in their stationary initial state and start their evolution at i = 0. Then for our data, the vectors v i are obtained by rearranging the N × N pixel images at instant i into N 2 × 1 vector. In our DMD analysis we employ {v 250 , v 251 , ...v 750 } (i.e. F = 500), as the modes obtained are the sharpest for this range. Let D (α) ( j, k) represent the ( j, k) th pixel of the α th eigenmode, ordered in terms of decreasing absolute eigenvalue. Since the caustics are localized around the central region of the domain, we use a zoomed-in region of size 512 × 512 pixels (i.e. N = 512) in our analysis, as shown in Fig. 1 . We find that the highest singular eigenvalue mode, namely D (1) , shown in Fig. 2 (a) highlights a straight-line caustic structure, which we refer to as the wavefront. The eigenvalues of other DMD modes decay exponentially. We observe that the position of the wavefront in D (1) has a systematic dependence on the Stokes number, and to extract this relation we detect the location of the wavefronts using edge detection techniques. Similarly, when we perform DMD analysis on a bi-disperse system D (1) ( j, k) shows two distinct wavefronts (see Fig. 2 (b)) corresponding to the two different Stokes numbers; and here DMD uses the velocity information to unambiguously extract the wave fronts. In particular, Fig. 3(a) shows the reduced phase space portrait of a typical particle which form the caustic 32 and Fig. 3(b) shows the particles overlaid on top of the DMD that demonstrates the DMD's ability to extract the caustic structures. Furthermore, we find that the intensity of each wavefront compared to the background, which we refer to as prominence 36 , depends on the corresponding initial particle concentrations in the system. We now prescribe a method to extract the position of the wavefront from DMD. We use a Sobel operator 37 such that the vertical gradient of the first DMD mode is given by where represents the 2D convolution operator 38 , and ∆y is the spacing in DMD along the y-axis. We then sum over the values in the x-direction to get a 1D function of y as where ∆x is the spacing in DMD along the x-axis, and we choose a square grid with ∆x = ∆y such that G x is by definition independent of the grid spacing and is dimensionless. In the next section we describe how the location and the value of the peaks in G x can be used to find the Stokes number of the particles and the relative initial concentrations in the case of a bi-disperse system. We extract the location of the peaks in G x for each St, as defined in eq. 7, to generate the plot in (b) and we find that the Y W F and St have a quadratic dependence. The G x is obtained from the DMD as described in Sec. 1.3 by simulating mono-disperse Stokes number particle systems to generate the plots in Fig. 4(a) which shows G x as a function of St. The alignment of the peaks in G x along a curve indicates a systematic dependence of the location of wavefront on the Stokes number. To extract this relation we first get the location of the caustic wavefront from the domain center using the position of the peaks in G x given by where arg max y gives the value of y for which G x is maximum. We then plot Y W F as a function of St as shown in Fig. 4(b) . Using a non-linear least squares fit method we find that the relation is of the form Y W F ∼ a St 2 + b St + c, with values of the parameters a, b, and c as indicated in Fig. 4(b) , where x represents St. Now, extrapolating the fit we find that Y W F = 0 at St = 0.3767 and becomes multi-valued for St > 3.0979, thus setting the limits on the validity of the relation. The three parameters in the relation can be estimated experimentally using calibrated measurements and the relation can be used to predict the St of new particle systems. In particular we demonstrate that the above method can be generalized to work in case of bi-disperse system. As shown in Fig. 2(b) , for a bi-disperse system the DMD has two sets of caustic wavefronts, corresponding to each Stokes number. Now we set one of the Stokes number fixed at one (St 1 = 1.0), vary the St 2 and find G x to generate the plots in Fig. 5 . The results in Fig. 5 show that even in the bi-disperse system the caustic wavefront has the same characteristic behaviour on the Stokes number as the mono-disperse system. In particular, the wavefront corresponding to St 1 has a fixed location and the wavefront due to St 2 preserves the dependence on Y W F of the mono-disperse system. Our studies with poly-disperse St systems show that the caustic wavefronts can be used to find the Stokes number of different particles in the system using the relation obtained from a mono-disperse system. Until now the bi-disperse systems that we considered had equal number of St 1 and St 2 particles, with uniform initial distribution in space. Now we study the variation in D (1) w.r.t. the change in relative number of particles. We observe that the intensity of the wavefront or the gradient in the DMD image depends on the number of particles or the initial uniform concentration, denoted by C (St) . The variable G x gives the gradient of D (1) along the vertical direction and the magnitude of the gradient indicates sharpness of the wavefronts (see Fig. 6(a) ). To measure the sharpness of the wavefront we define a "prominence" parameter, P, as the sum of the non zero values of| G x | in the neighbourhood of the wavefront, which takes into account multiple peaks in the vicinity of the wavefront. We find that the prominence of the wavefront has a systematic dependence on the initial concentration of the corresponding Stokes number particles and from Fig. 6(b) we find that on a log-log plot the relation is Notice that the peaks in G x corresponding to St 1 are aligned at the same location along y, whereas the peaks due to St 2 show similar trend as the plots for mono-disperse systems in Fig. 4(a) . ). Notice that the peaks corresponding to each wavefront is not unique and have a finite spread in y. In (b) the relation between the prominence corresponding to St 1 and St 2 are given by P 1 , P 2 respectively, as a function of the initial particle concentrations C is shown in a log-log plot. The linear fit shows that the ratios of the peaks of| G | and the ratios of the concentration are related by a power-law, with a power close to -1. linear, with a slope approximately equal to -1. This implies that in a bi-disperse system the ratio of the prominence is inversely related to the ratio of initial concentrations. We can use this relation to predict the concentration of various Stokes number particles in a system. We study the dynamics of inertial particles in a Taylor-Green flow with periodic boundary conditions in 2D. In a minimal model of inertial particles we observe that for a mono-disperse Stokes number system, starting from a uniform distribution of stationary particles, the particle distribution forms caustics in the strain dominated region of the flow. We use the DMD method to analyze the PIV-like time-series data of the spatio-temporal particle distribution and find that the largest absolute eigenvalue mode is effective in extracting the caustic wavefront-like structure. We notice that (a) the position of the wavefront depends on the particle Stokes number and employ standard image processing techniques to quantitatively extract a quadratic relation. Using this relation we can predict the Stokes number from the wavefront position. Furthermore, we find that for a bi-disperse system the DMD is able to extract two different wavefronts corresponding to each Stokes number and the positions of each wavefront follow the same quadratic relation as in the case of mono-disperse system. We also observe that (b) the sharpness of the wavefront in the DMD, measured in terms of prominence, depends on the initial particle concentration and find that for a bi-disperse Stokes number system the ratio of the wavefront prominence is inversely proportional to the corresponding Stokes number initial concentration. Hence the measurement of prominence can be used to estimate the concentration of the corresponding Stokes number particles. We propose that the DMD technique can be used to analyze real experimental PIV data of caustics and perform similar analysis to extract information about the Stokes numbers and concentrations of the particles. In future we will consider detailed Navier-Stokes equation for the self-consistent evolution of the velocities and analyze the caustic structures in turbulent flows. The particle dynamics we use in our study is a special case of 17 where η is the white noise and R = ρ f /(ρ p + ρ f /2) for density of particle ρ p and density of fluid ρ f . We present our results for η = 0 and R = 0, and use small values of these parameters to verify the stability of our results. The N × N pixel image at k th instant, I k (i, j), is rearranged into an N 2 × 1 vector X k (m) (note that we subtract the mean value from X k (m)). Now the F sequence of image frames are appended together to form N 2 × N T matrix Y . Let the N 2 × (F − 1) dimensional matrix formed from the first (F − 1) frames be Y 1 and the last (F − 1) frames be Y 2 , i.e. We find the singular value decomposition (SVD) of the matrix Y 1 , such that Y 1 = UΣV * , where U is a N 2 × N 2 complex unitary matrix, Σ is an N 2 × (F − 1) rectangular diagonal matrix with non-negative real numbers on the diagonal, and V is a (F − 1) × (F − 1) real or complex unitary matrix. Now choose a lower dimensional SVD matrices made up of first n T (<< F ) columns, represented byŨ,Ṽ , and the first n T × n T block of Σ asΣ. Define a matrix n T × n T matrix A as Then the dynamic modes are the N 2 × 1 eigenvectors of matrix A that is rearranged into N × N matrix. The movies referred to in the text are available in the additional materials. C.1 See "Vid1.mp4" C.2 See "Vid2.mp4" 7/9 Oceanic mesoscale eddies as revealed by Lagrangian coherent structures A new mixing diagnostic and Gulf oil spill movement Surface coastal circulation patterns by in-situ detection of Lagrangian coherent structures Forecasting sudden changes in environmental pollution patterns The geometry of inertial particle mixing in urban flows, from deterministic and random displacement models Toxic trace elements: preferential concentration in respirable particles Density-ratio effects on the capture of suspended particles in aquatic systems Preferential concentration of cloud droplets by turbulence: Effects on the early evolution of cumulus cloud droplet spectra Inertial particle dynamics in a hurricane Volcanic ash in the free atmosphere: A dynamical systems approach The flow physics of COVID-19 The dispersion of spherical droplets in source-sink flows and their relevance to the COVID-19 pandemic First-passage-time problem for tracers in turbulent flows applied to virus spreading Does Air Pollution Influence COVID-19 Outbreaks? Numerical investigation of aerosol transport in a classroom with relevance to COVID-19 Clustering criterion for inertial particles in two-dimensional time-periodic and threedimensional steady flows Equation of motion for a small rigid sphere in a nonuniform flow Fractal clustering of inertial particles in random flows Heavy particle concentration in turbulence at dissipative and inertial scales Caustics and clustering in the vicinity of a vortex Particle Image Velocimetry Full Lagrangian methods for calculating particle concentration fields in dilute gas-particle flows Applications of the dynamic mode decomposition Dynamic mode decomposition of numerical and experimental data Dynamic mode decomposition: data-driven modeling of complex systems On dynamic mode decomposition: Theory and applications Anomalous diffusion of tracer in convection rolls Nonlinear Response of Inertial Tracers in Steady Laminar Flows: Differential and Absolute Negative Mobility Two-dimensional turbulence: A physicist approach Caustics in turbulent aerosols Caustic activation of rain showers Inertial-particle dynamics in turbulent flows: Caustics, concentration fluctuations and random uncorrelated motion Chaotic dynamics of particle dispersion in fluids Modal analysis of fluid flows: An overview Dynamic mode decomposition of numerical and experimental data The Finest Peaks-Prominence and Other Mountain Measures Principles of filter design Computer Vision