On spatial instability of electrically forced axisymmetric jets with variable applied field

This paper considers the problem of instability of electrically forced axisymmetric jets with respect to spatially growing disturbances and in the presence of a variable applied electric field. A mathematical model, which is developed for the dependent variables of such disturbances, is based on the relevant approximated versions of the equations of the electrohydrodynamics for an electrically forced jet flow. The approximations include the assumptions that the length scale along the axial direction of the jet is much larger than that in the radial direction of the jet and the disturbances are axisymmetric and infinitesimal in amplitude. For neutral temporal stability boundary, we find, in particular, two new spatial modes of instabilities under certain conditions. Both modes are found to be enhanced with increasing the strength of the field. The more dominant instability mode is found to exist for a wider range of values of the wave number in the axial direction. The effect of variable applied electric field is found to increase the growth rates of the disturbances but operate over a more restricted domain in the axial wave number.     

On spatial instability of an electrically forced non-axisymmetric jet with curved centerline

This paper considers the problem of spatial instability of an electrically forced non-axisymmetric jet with curved centerline. A mathematical model, which is developed for the spatially growing oscillations of the centerline of the jet, is based on the relevant approximated versions of the equations of the motion for such electrically forced jet flow. The approximations include the assumptions that the radius of curvature of the centerline of the jet is much larger than the radius of the jet and the spatially growing disturbances are infinitesimal in amplitude. For the neutral temporal stability boundary, we identify, in particular, new spatial, conducting and viscous modes of instabilities which travel axially in the direction of increasing axial variable and are enhanced with increasing either the surface charge or the strength of the applied electric field. For given values of the parameters, there is a critical wave number at which the instability of these modes is maximized. The range of values of the wave number and the frequency of these instability modes increases with the externally imposed electric field.     

Homogenization of monotone parabolic problems with several temporal scales

In this paper we homogenize monotone parabolic problems with two spatial scales and any number of temporal scales. Under the assumption that the spatial and temporal scales are well-separated in the sense explained in the paper, we show that there is an H-limit defined by at most four distinct sets of local problems corresponding to slow temporal oscillations, slow resonant spatial and temporal oscillations (the ??slow?? self-similar case), rapid temporal oscillations, and rapid resonant spatial and temporal oscillations (the ??rapid?? self-similar case), respectively.     

Viscous linear instability of a sheared liquid-gas interface - Effect of a velocity deficit near the interface

The initial state of liquid atomization by a fast gas stream is considered by viscous linear spatial stability analysis for parallel two-fluid flow. The unbounded basic velocity profile is characterized by different asymptotic velocities and a velocity deficit near the interface. We examine the influence of the velocity deficit on the spatial growth rates of two competing modes originating from the Kelvin-Helmholtz and viscosity contrast mechanisms. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Resonant instability and nonlinear wave interactions in electrically forced jets

We investigate the problem of linear temporal instability of the modes that satisfy the dyad resonance conditions and the associated nonlinear wave interactions in jets driven by either a constant or a variable external electric field. A mathematical model, which is developed and used for the temporally growing modes with resonance and their nonlinear wave interactions in electrically driven jet flows, leads to equations for the unknown amplitudes of such waves. These equations are solved for both water and glycerol jet cases, and the expressions for the dependent variables of the corresponding modes are determined. The results of the generated data for these dependent variables versus time indicate, in particular, that the instability resulted from the nonlinear interactions of such modes is mostly quite strong but can also lead to significant reduction in the jet radius.     

On the resonant reflection of weak,nonlinear sound waves off an entropy wave

We analyze the resonant reflection of very weak, nonlinear sound waves off a weak sawtooth entropy wave for spatially periodic solutions of the one‐dimensional, nonisentropic gas dynamics equations. The case of an entropy wave with a sawtooth profile is of interest because the oscillations of the reflected sound waves are nondispersive with frequency independent of their wavenumber, leading to an unusual type of nonlinear dynamics. On an appropriate long time scale, we show that a complex amplitude function for the spatial profile of the sound waves satisfies a degenerate quasilinear Schrödinger equation. We present some numerical solutions of this equation that illustrate the generation of small spatial scales by a resonant four‐wave cascade and front propagation in compactly supported solutions.     

Nonlinear instability of a viscous liquid jet

We investigate the weakly nonlinear temporal instability of an axisymmetric Newtonian liquid jet. Early nonlinear studies on the capillary instability of inviscid liquid jets were carried up to the third order contributions to the jet deformation and showed the nonlinear interaction between different modes. A recent study on the weakly nonlinear instability of planar Newtonian liquid sheets revealed the role of the liquid viscosity in the sheet stability behavior and showed a complicated influence [1]. Here, the instability of a liquid jet is examined as the axisymmetric counterpart of the sheet, in search for corresponding insight into the role of the liquid viscosity in the jet instability mechanism. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Permanent Wave Structures and Resonant Triads in a Layered Fluid

A closed three layer fluid with small density differences between the layers has two closely related modes of gravity wave propagation. The nonlinear interactions between the wave modes are investigated, particularly the nearly resonant or significant interactions. Permanent wave solutions are calculated, and it is shown that a permanent wave of the slower mode can generate resonantly a wave harmonic of the faster mode. The equations governing resonant triads of the two modes are derived, and solutions having a permanent structure are calculated from them. It is found that some resonant triad solutions vanish when the triad is embedded in the set of all harmonics with wavenumbers in its neighborhood     

Spatio-temporal symmetries and bifurcations via bi-orthogonal decompositions

Summary A tool for analyzing spatio-temporal complex physical phenomena was recently proposed by the authors, Aubry et al. [5]. This tool consists in decomposing a spatially and temporally evolving signal into orthogonal temporal modes (temporal “structures”) and orthogonal spatial modes (spatial “structures”) which are coupled. This allows the introduction of a temporal configuration space and a spatial one which are related to each other by an isomorphism. In this paper, we show how such a tool can be used to analyze space-time bifurcations, that is, qualitative changes in both space and time as a parameter varies. The Hopf bifurcation and various spatio-temporal symmetry related bifurcations, such as bifurcations to traveling waves, are studied in detail. In particular, it is shown that symmetry-breaking bifurcations can be detected by analyzing the temporal and spatial eigenspaces of the decomposition which then lose their degeneracy, namely their invariance under the symmetry. Furthermore, we show how an extension of the theory to “quasi-symmetries” permits the treatment of nondegenerate signals and leads to an exponentially decreasing law of the energy spectrum. Examples extracted from numerically obtained solutions of the Kuramoto-Sivashinsky equation, a coupled map lattice, and fully developed turbulence are given to illustrate the theory.     

Variants of Dynamic Mode Decomposition: Boundary Condition, Koopman, and Fourier Analyses

Dynamic mode decomposition (DMD) is an Arnoldi-like method based on the Koopman operator. It analyzes empirical data, typically generated by nonlinear dynamics, and computes eigenvalues and eigenmodes of an approximate linear model. Without explicit knowledge of the dynamical operator, it extracts frequencies, growth rates, and spatial structures for each mode. We show that expansion in DMD modes is unique under certain conditions. When constructing mode-based reduced-order models of partial differential equations, subtracting a mean from the data set is typically necessary to satisfy boundary conditions. Subtracting the mean of the data exactly reduces DMD to the temporal discrete Fourier transform (DFT); this is restrictive and generally undesirable. On the other hand, subtracting an equilibrium point generally preserves the DMD spectrum and modes. Next, we introduce an ??optimized?? DMD that computes an arbitrary number of dynamical modes from a data set. Compared to DMD, optimized DMD is superior at calculating physically relevant frequencies, and is less numerically sensitive. We test these decomposition methods on data from a two-dimensional cylinder fluid flow at a Reynolds number of?60. Time-varying modes computed from the DMD variants yield low projection errors.     

Dynamic Mode Decomposition for Swirling Jet Flow Undergoing Vortex Breakdown

Swirling jets undergoing vortex breakdown occur in many technical applications, e.g. vortex burners, turbines and jet engines. At the stage of vortex breakdown the flow is dominated by a conical shear layer and a large recirculation zone around the jet axis. We performed Large-Eddy Simulations (LES) of compressible swirling jet flows at Re=5000, Ma=0.6 in the high swirl number regime (S=1). A nozzle is included in our computational setup to account for more realistic inflow conditions. The obtained velocity fields are analyzed by means of temporal and spatial dynamic mode decomposition (DMD) to get further insight into the characteristic structures dominating the flow. We present eigenvalue spectra for the case under consideration and discuss the stability behaviour in time and space. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modeling the disintegration of modulated liquid jets using volume-of-fluid (VOF) methodology

In this study, we present the numerical investigations on the effect of finite velocity modulations imposed on an otherwise unperturbed cylindrical liquid jet issuing into stagnant gas. Sinusoidal velocity fluctuations of finite frequency and amplitude are imposed at the liquid jet inlet and the resulting liquid jet surface deformation is captured using a volume of fluid (VOF) methodology, utilizing compressive interface capturing scheme for arbitrary meshes (CICSAM) scheme. Variation of the simulation parameters, comprising of the mean liquid jet velocity, modulation amplitude and frequency grouped together using a set of non-dimensional parameters, leads to the formation of a wide gamut of reproducible liquid structures such as waves, upstream/downstream directed bells, chains of droplets similar to those observed in experiments. Elaborate tests on the effect of injection velocity and inlet jet diameter are investigated to characterize the breakup process. The computations efficiently capture the diverse flow structures generated by the evolving modulated liquid jet inclusive of several non-linear dynamics such as growth of surface waves, ligament interaction with shear vortices and its subsequent thinning process. The simulations identify the deterministic behavior of modulated liquid jets to predict liquid disintegration modes under given set of non-dimensional parameters.     

On the stability of resonant rotation of a symmetric satellite in an elliptical orbit

We deal with the stability problem of resonant rotation of a symmetric rigid body about its center of mass in an elliptical orbit. The resonant rotation is a planar motion such that the body completes one rotation in absolute space during two orbital revolutions of its center of mass. In [1–3] the stability analysis of the above resonant rotation with respect to planar perturbations has been performed in detail.In this paper we study the stability of the resonant rotation in an extended formulation taking into account both planar and spatial perturbations. By analyzing linearized equations of perturbed motion, we found eccentricity intervals, where the resonant rotation is unstable. Outside of these intervals a nonlinear stability study has been performed and subintervals of formal stability and stability for most initial data have been found. In addition, the instability of the resonant rotation was established at several eccentricity values corresponding to the third and fourth order resonances.Our study has also shown that in linear approximation the spatial perturbations have no effect on the stability of the resonant rotation, whereas in a nonlinear system they can lead to its instability at some resonant values of the eccentricity.     

Abstract principal component analysis

We present the basic idea of abstract principal component analysis(APCA)as a general approach that extends various popular data analysis techniques such as PCA and GPCA.We describe the mathematical theory behind APCA and focus on a particular application to mode extractions from a data set of mixed temporal and spatial signals.For illustration,algorithmic implementation details and numerical examples are presented for the extraction of a number of basic types of wave modes including,in particular,dynamic modes involving spatial shifts.     

The trajectory and stability of a spiralling liquid jet: Viscous theory

We examine a spiralling slender viscous jet emerging from a rapidly rotating orifice, extending Wallwork et al. [I.M. Wallwork, S.P. Decent, A.C. King, R.M.S.M. Schulkes, The trajectory and stability of a spiralling liquid jet. Part 1. Inviscid theory, J. Fluid Mech. 459 (2002) 43–65] by incorporating viscosity. The effects of viscosity on the trajectory of the jet and its linear instability are determined using a mixture of computational and asymptotic methods, and verified using experiments. A non-monotonic relationship between break-up length and rotation rate is demonstrated with the trend varying with viscosity. The sizes of the droplets produced by this instability are determined by considering the most unstable wave mode. It is also found that there is a non-monotonic relationship between droplet size and viscosity. Satellite droplet formation is also considered by analysing very short wavelength modes. The effects of long wavelength modes are examined, and a wave which propagates down the trajectory of the jet is identified for the highly viscous case. A comparison between theoretical and experimental results is made, with favourable agreement. In particular, a quantitative comparison is made between droplet sizes predicted from the theory with experimental observations, with encouraging agreement obtained. Four different types of break-up are identified in our experiments. The experimentally observed break-up mechanisms are discussed in light of our theory.     

Designing a contact process: the piecewise-homogeneous process on a finite set with applications

We consider how to choose the reproduction rates in a one-dimensional contact process on a finite set to maximize the growth rate of the extinction time with the population size. The constraints are an upper bound on the average reproduction rate, and that the rate profile must be piecewise constant. We show that the optimum growth rate is achieved by a rate profile with at most two rates, and we characterize the solution in terms of a “spatial correlation length” of the supercritical process. We examine the analogous problem for the simpler biased voter model, for which we completely characterize the optimum profile. The contact process proofs make use of a planar-graph duality in the graphical representation, due to Durrett and Schonmann.     

Multiple-scale Homogenization for Weakly Nonlinear Conservation Laws with Rapid Spatial Fluctuations

We consider hyperbolic conservation laws with rapid periodic spatial fluctuations and study initial value problems that correspond to small perturbations about a steady state. Weakly nonlinear solutions are computed asymptotically using multiple spatial and temporal scales to capture the homogenized solution as well as its long-term behavior. We show that the linear problem may be destabilized through interactions between two solution modes and the periodic structure. We also show that a discontinuity, either in the initial data or due to shock formation, introduces rapid spatial and temporal fluctuations to leading order in its zone of influence. The evolution equations we derive for the homogenized leading-order solution are more general than their counterparts for conservation laws having no rapid spatial variations. In particular, these equations may be diffusive for certain general flux vectors. Selected examples are solved numerically to substantiate the asymptotic results.     

Modeling and Anomalous Cluster Detection for Point Processes Using Process Convolutions

We present a model using process convolutions, which describes spatial and temporal variations of the intensity of events that occur at random geographical locations. An inhomogeneous Poisson process is used to model the intensity over a spatial region with multiplicative spatial and temporal covariate effects. Temporal variation in the structure of the intensity is obtained by employing a time-varying process for the convolution. Use of a compactly supported kernel in the convolution improves the computational efficiency. Additionally, anomalous cluster detection in the event rates is developed based on exceedance probabilities. The methods are demonstrated on data of major crimes in Cincinnati during 2006. Supplementary materials for this article are available online.     

Numerical spectral analysis of temporal stability of laminar duct flows with constant cross sections

Problems related to the temporal stability of laminar viscous incompressible flows in ducts with a constant cross section are formulated, justified, and numerically solved. For the systems of ordinary differential and algebraic equations obtained by a spatial approximation, a new dimension reduction technique is proposed and substantiated. The solutions to the reduced systems are decomposed over subspaces of modes, which considerably improves the computational stability of the method and reduces the computational costs as compared with the usual decompositions over individual modes. The optimal disturbance problem is considered as an example. Numerical results for Poiseuille flows in a square duct are presented and discussed.     

Multiscale analysis for diffusion-driven neutrally stable states

The Turing instabilities for reaction–diffusion systems are studied from the Fourier normal modes which appear by searching the solution obtained from linearization of the reaction–diffusion system at the spatially homogeneous steady state. The linear stability analysis is only appropriate when the temporal eigenvalues associated to every given spatial eigenvalue have non-zero real part. If the real part of the temporal eigenvalue in a normal mode is equal to zero there is no enough information coming from the linearized system. Given an arbitrary spatial eigenvalue, by equating to zero the real part of the corresponding temporal eigenvalue will lead to a neutral stability manifold in the parameter space. If for a given spatial eigenvalue the other parameters in the reaction–diffusion process drive the system to the neutral manifold, then neither stability nor instability can be warranted by the usual linear analysis. In order to give a sketch of the nonlinear analysis we use a multiple scales method. As an application, we analyze the behavior of solutions to the Schnakenberg trimolecular reaction kinetics in the presence of diffusion.