Effect of Volume Energy Supply on the Stability of a Subsonic Vortex Flow

Calculations of the stability of an axisymmetric vortex flow of viscous heat-conducting gas with volume energy supply are presented. The unperturbed axisymmetric vortex flow was found numerically using a quasi-cylindrical approximation of the Navier-Stokes equations under the assumption of constant peripheral-velocity circulation in the ambient co-current flow. The volume energy supply in the viscous vortex core was modeled by an additional source term in the energy equation. The stability characteristics of the viscous vortex flow in a longitudinal vortex with respect to both axisymmetric and non-axisymmetric three-dimensional waves traveling along the vortex axis and corresponding to both positive and negative values of the azimuthal wave number were found using the time-dependent formulation of the linear stability theory for compressible three-dimensional plane-parallel flows.     

Asymptotic Behavior of Solutions to the Compressible Navier–Stokes Equation Around a Parallel Flow

The initial boundary value problem for the compressible Navier–Stokes equation is considered in an infinite layer of . It is proved that if the Reynolds and Mach numbers are sufficiently small, then strong solutions to the compressible Navier–Stokes equation around parallel flows exist globally in time for sufficiently small initial perturbations. The large time behavior of the solution is described by a solution of a one-dimensional viscous Burgers equation. The proof is given by a combination of spectral analysis of the linearized operator and a variant of the Matsumura–Nishida energy method.     

The Linear Stability of Shock Waves for the Nonlinear Schrödinger–Inviscid Burgers System

We investigate the coupling between the nonlinear Schrödinger equation and the inviscid Burgers equation, a system which models interactions between short and long waves, for instance in fluids. Well-posedness for the associated Cauchy problem remains a difficult open problem, and we tackle it here via a linearization technique. Namely, we establish a linearized stability theorem for the Schrödinger–Burgers system, when the reference solution is an entropy-satisfying shock wave to Burgers equation. Our proof is based on suitable energy estimates and on properties of hyperbolic equations with discontinuous coefficients. Numerical experiments support and expand our theoretical results.     

Stability of a compressible swirling flow in a circular pipe

The stability of an axisymmetric flow of viscous gas in a circular pipe, which models the Burgers vortex in the pipe axis neighborhood, is studied within the linear theory framework. Neutral curves for the most unstable disturbances are calculated. The influence of the characteristic Mach number on the flow stability is investigated. It is shown that for a given model velocity distribution the Mach number affects only the temperature and pressure profiles of the main undisturbed flow. In this case, for the disturbance types considered, as the Mach number increases, the critical Reynolds number corresponding to loss of stability decreases. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 35–41, January–February, 1999. The work received financial support from the Russian Foundation for Basic Research (project No. 96-01-00586).     

Evolution Equations for the Surface Concentration of an Insoluble Surfactant; Applications to the Stability of an Elongating Thread and a Stretched Interface

The general form of the convection–diffusion equation governing the evolution of the surface concentration of an insoluble surfactant over an evolving interface is reviewed and discussed for three-dimensional, axisymmetric, and two-dimensional configurations. The linearized form of the evolution equation is then derived around cylindrical and planar shapes in a framework that is suitable for carrying out a linear stability analysis for axisymmetric or two-dimensional perturbations. Particular attention is paid to the cases of quiescent unperturbed fluids, unidirectional shear flow, and elongational flow. By way of application, the linearized transport equations are combined with Stokes-flow hydrodynamics to investigate the stability of an elongating cylindrical viscous thread suspended in an ambient viscous fluid or in a vacuum, and the stability of a two-dimensional interface separating two semi-infinite fluids and stretched under the action of an orthogonal stagnation-point flow. The results illustrate the possibly important role of the surfactant on the linear growth of periodic waves on the cylindrical interface, and reveal that the surfactant has no effect on the stability of the planar interface.     

Stability of an axisymmetric swirled flow in a supersonic cocurrent stream for volume energy supply in the viscous vortex core

The results of numerical calculations of the stability of axisymmetric swirled flows in a viscous vortex embedded in a supersonic cocurrent stream with a constant circulation of the azimuthal velocity component are presented. The stability characteristics of the swirled three-dimensional viscous flow in the streamwise vortex are determined on the basis of the linearized system of Navier-Stokes equations for a viscous heat-conducting gas under the assumption that the basic undisturbed flow is locally plane-parallel. The disturbed flow stability is studied in the temporal formulation with respect to both symmetric and asymmetric three-dimensional waves traveling along the vortex axis and corresponding to both positive and negative values of the azimuthal wavenumber. It is shown that at external inviscid flow Mach numbers M = 2 and 3 thermal energy supply in a small region near the vortex axis leads to considerable restructuring of the basic undisturbed flow in the vicinity of the vortex core, the growth of the adverse pressure gradient along the vortex axis, and a significant change in the small perturbation stability and behavior.__________Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, 2005, pp. 71–80. Original Russian Text Copyright © 2005 by Kazakov.     

Stability of Gasless Combustion Fronts in One-Dimensional Solids

For gasless combustion in a one-dimensional solid, we show a type of nonlinear stability of the physical combustion front: if a perturbation of the front is small in both a spatially uniform norm and an exponentially weighted norm, then the perturbation stays small in the spatially uniform norm and decays in the exponentially weighted norm, provided the linearized operator has no eigenvalues in the right half-plane other than zero. Using the Evans function, we show that the zero eigenvalue must be simple. Factors that complicate the analysis are: (1) the linearized operator is not sectorial, and (2) the linearized operator has good spectral properties only when the weighted norm is used, but then the nonlinear term is not Lipschitz. The result is nevertheless physically natural. To prove it, we first show that when the weighted norm is used, the semigroup generated by the linearized operator decays on a subspace complementary to the operator’s kernel, by showing that it is a compact perturbation of the semigroup generated by a more easily analyzed triangular operator. We then use this result to help establish that solutions stay small in the spatially uniform norm, which in turn helps establish nonlinear convergence in the weighted norm.     

The Modulational Instability for a Generalized Korteweg–de Vries Equation

We study the spectral stability of a family of periodic standing wave solutions to the generalized Korteweg–de Vries in a neighborhood of the origin in the spectral plane using what amounts to a rigorous Whitham modulation theory calculation. In particular we are interested in understanding the role played by the null directions of the linearized operator in the stability of the traveling wave to perturbations of long wavelength. A study of the normal form of the characteristic polynomial of the monodromy map (the periodic Evans function) in a neighborhood of the origin in the spectral plane leads to two different instability indices. The first, an orientation index, counts modulo 2 the total number of periodic eigenvalues on the real axis. This index is a generalization of the one which governs the stability of the solitary wave. The second, a modulational instability index, provides a necessary and sufficient condition for the existence of a long-wavelength instability. This index is essentially the quantity calculated by Hǎrǎguş and Kapitula in the small amplitude limit. Both of these quantities can be expressed in terms of the map between the constants of integration for the ordinary differential equation defining the traveling waves and the conserved quantities of the partial differential equation. These two indices together provide a good deal of information about the spectrum of the linearized operator. We sketch the connection of this calculation to a study of the linearized operator—in particular we perform a perturbation calculation in terms of the Floquet parameter. This suggests a geometric interpretation attached to the vanishing of the orientation index previously mentioned.     

Dynamic stability of spinning viscoelastic cylinders at finite deformation

The study of spinning axisymmetric cylinders undergoing finite deformation is a classic problem in several industrial settings – the tire industry in particular. We present a stability analysis of spinning elastic and viscoelastic cylinders using ARPACK to compute eigenvalues and eigenfunctions of finite element discretizations of the linearized evolution operator. We show that the eigenmodes correspond to N-peak standing or traveling waves for the linearized problem with an additional index describing the number of oscillations in the radial direction. We find a second hierarchy of bifurcations to standing waves where these eigenvalues cross zero, and confirm numerically the existence of finite-amplitude standing waves for the nonlinear problem on one of the new branches. In the viscoelastic case, this analysis permits us to study the validity of two popular models of finite viscoelasticity. We show that a commonly used finite deformation linear convolution model results in non-physical energy growth and finite-time blow-up when the system is perturbed in a linearly unstable direction and followed nonlinearly in time. On the other hand, Sidoroff-style viscoelastic models are seen to be linearly and nonlinearly stable, as is physically required.     

Asymptotic Stability of N-Solitary Waves of the FPU Lattices

We study stability of N-solitary wave solutions of the Fermi-Pasta-Ulam (FPU) lattice equation. Solitary wave solutions of the FPU lattice equation cannot be characterized as critical points of conservation laws due to the lack of infinitesimal invariance in the spatial variable. In place of standard variational arguments for Hamiltonian systems, we use an exponential stability property of the linearized FPU equation in a weighted space which is biased in the direction of motion. The dispersion of the linearized FPU equation balances the potential term for low frequencies, whereas the dispersion is superior for high frequencies.We approximate the low frequency part of a solution of the linearized FPU equation by a solution to the linearized Korteweg-de Vries (KdV) equation around an N-soliton solution. We prove an exponential stability property of the linearized KdV equation around N-solitons by using the linearized Bäcklund transformation and use the result to analyze the linearized FPU equation.     

轴对称螺旋流解析解的探讨

针对一类轴对称螺旋流问题进行了解析分析与讨论．发现若干N-S方程的叠加精确解,其中包括Burgers涡与Hagen－Poiseuille流的叠加     

Properties of vortices in the self-similar turbulent jet

Instantaneous, two-dimensional velocity measurements were conducted in the axial plane of a self-similar turbulent axisymmetric jet. The velocity fields were high-pass filtered to expose the vortical structures. An automated method was used to identify the radial and axial coordinates of the vortex centers and rotational sense, and to measure their size, circulation, vorticity, and energy. New insights into turbulent jets are obtained by plotting statistical distributions for vortex properties as functions of Reynolds number and radial position. While the probability of finding a vortex is uniform up to the edge of the jet, the strongest eddies in the high-pass filtered field occur near the jet axis. The average circulation is directly proportional to the vortex size. The Reynolds number strongly affects the average vorticity, circulation, and energy of the eddies. However, the normalized curves show a good collapse implying that the jet is indeed self-similar. Results for the left and right half-planes of the jet are also presented. Interestingly, we find that contrary to customary drawings of jet flows, a substantial number of both clockwise and counter-clockwise rotating eddies exist on both sides of the jet axis, with almost equal numbers of oppositely rotating vortices close to the jet axis. Further, the disparity in the number of oppositely rotating eddies in each half-plane increases with the eddy size. Nevertheless, these results are consistent with the well-known radial vorticity distribution of axisymmetric jets.     

Stability of a viscous subsonic swirl flow

The results of calculating the stability of a three-dimensional swirl flow of a viscous heat-conducting gas are presented. The stability characteristics are determined using the linear time-dependent theory of plane-parallel flow stability. The main undisturbed axisymmetric vortex flow was determined numerically using a quasi-cylindrical approximation for the complete set of Navier-Stokes equations. The circulation of the peripheral velocity in the cocurrent flow surrounding the viscous vortex core was assumed to be constant. In analyzing the stability, nonaxisymmetric perturbations in the shape of waves traveling along the vortex axis with both positive and negative wavenumbers were considered; in these two cases the perturbation rotation is either the same or opposite in sense to the rotation in the vortex core. Neutral stability curves are determined for various values of the swirling parameter and the cocurrent flow Mach number. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 50–59, May–June, 1998.     

Deformation and Filtration Near a Well in a Porous Stratum with Linear Memory

The linearized equations for saturated elastic porous media and for surrounding elastic rock are solved simultaneously; and the Volterra principle is used to derive an integro-differential filtration equation for a homogeneous weakly compressible fluid in an axisymmetric stratum with linear memory and central well. An analytical expression for porosity variation is obtained and then used to determine the permeability coefficient. The solutions are analyzed for the case where the stratum exhibits memory described by regular and singular kernels of the integral operator     

Existence and Stability of Compressible Current-Vortex Sheets in Three-Dimensional Magnetohydrodynamics

Compressible vortex sheets are fundamental waves, along with shocks and rarefaction waves, in entropy solutions to multidimensional hyperbolic systems of conservation laws. Understanding the behavior of compressible vortex sheets is an important step towards our full understanding of fluid motions and the behavior of entropy solutions. For the Euler equations in two-dimensional gas dynamics, the classical linearized stability analysis on compressible vortex sheets predicts stability when the Mach number \(M > \sqrt{2}\) and instability when \(M < \sqrt{2}\) ; and Artola and Majda’s analysis reveals that the nonlinear instability may occur if planar vortex sheets are perturbed by highly oscillatory waves even when \(M > \sqrt{2}\) . For the Euler equations in three dimensions, every compressible vortex sheet is violently unstable and this instability is the analogue of the Kelvin–Helmholtz instability for incompressible fluids. The purpose of this paper is to understand whether compressible vortex sheets in three dimensions, which are unstable in the regime of pure gas dynamics, become stable under the magnetic effect in three-dimensional magnetohydrodynamics (MHD). One of the main features is that the stability problem is equivalent to a free-boundary problem whose free boundary is a characteristic surface, which is more delicate than noncharacteristic free-boundary problems. Another feature is that the linearized problem for current-vortex sheets in MHD does not meet the uniform Kreiss–Lopatinskii condition. These features cause additional analytical difficulties and especially prevent a direct use of the standard Picard iteration to the nonlinear problem. In this paper, we develop a nonlinear approach to deal with these difficulties in three-dimensional MHD. We first carefully formulate the linearized problem for the current-vortex sheets to show rigorously that the magnetic effect makes the problem weakly stable and establish energy estimates, especially high-order energy estimates, in terms of the nonhomogeneous terms and variable coefficients. Then we exploit these results to develop a suitable iteration scheme of the Nash–Moser–Hörmander type to deal with the loss of the order of derivative in the nonlinear level and establish its convergence, which leads to the existence and stability of compressible current-vortex sheets, locally in time, in three-dimensional MHD.     

Singularities and spectral asymptotics of a random nonlinear wave in a nondispersive system

In this paper, we study the propagation of high-intensity acoustic noise in free space and in waveguide systems. A mathematical model generalizing the Burgers equation is used. It describes the nonlinear wave evolution inside tubes of variable cross-section, as well as in ray tubes, if the geometric approximation for heterogeneous media is used. The generalized equation transforms to the common Burgers equation with a dissipative parameter, known as the “Reynolds–Goldberg number”. In our model, this number depends on the distance travelled by the wave. With a zero “viscous” dissipative term, the model reduces to the Riemann (or Hopf) equation. Its solution presents the field by an implicit function. The spectral form of this solution makes it possible to derive explicit expressions for both dynamic and statistical characteristics of intense waves. The use of a spectral approach allowed us to describe the high-intensity noise in media with zero and finite viscosity. Applicability conditions of these solutions are defined. Since the phase matching is fulfilled for any triplet of interacting spectral components, there is an avalanche-like increase in the number of harmonics and the formation of shocks. The relationship between these discontinuities and other singularities and the high-frequency asymptotic of intense noise is studied. The possibility is shown to enhance nonlinear effects in waveguide systems during the evolution of noise.     

Numerical simulation of entry flow of the IUPAC-LDPE melt

Numerical simulations have been undertaken for the creeping entry flow of a well-characterized polymer melt (IUPAC-LDPE) in a 4:1 axisymmetric and a 14:1 planar contraction. The fluid has been modeled using an integral constitutive equation of the K-BKZ type with a spectrum of relaxation times (Papanastasiou–Scriven–Macosko or PSM model). Numerical values for the constants appearing in the equation have been obtained from fitting shear viscosity and normal stress data as measured in shear and elongational data from uniaxial elongation experiments. The numerical solutions show that in the axisymmetric contraction the vortex in the reservoir first increases with increasing flow rate (or apparent shear rate), goes through a maximum and then decreases following the behavior of the uniaxial elongational viscosity. For the planar contraction, the vortex diminishes monotonically with increasing flow rate following the planar extensional viscosity. This kinematic behavior is not in agreement with recent experiments. The PSM strain-memory function of the model is then modified to account for strain-hardening in planar extension. Then the vortex pattern shows an increase in both axisymmetric and planar flows. The results for planar flow are compared with recent experiments showing the correct trend.     

Numerical comparison of hybridized discontinuous Galerkin and finite volume methods for incompressible flow

In this study, a hybridizable discontinuous Galerkin method is presented for solving the incompressible Navier–Stokes equation. In our formulation, the convective part is linearized using a Picard iteration, for which there exists a necessary criterion for convergence. We show that our novel hybridized implementation can be used as an alternative method for solving a range of problems in the field of incompressible fluid dynamics. We demonstrate this by comparing the performance of our method with standard finite volume solvers, specifically the well‐established finite volume method of second order in space, such as the icoFoam and simpleFoam of the OpenFOAM package for three typical fluid problems. These are the Taylor–Green vortex, the 180‐degree fence case and the DFG benchmark. Our careful comparison yields convincing evidence for the use of hybridizable discontinuous Galerkin method as a competitive alternative because of their high accuracy and better stability properties. Copyright © 2014 John Wiley & Sons, Ltd.     

Free surface waves in equilibrium with a vortex

Finite amplitude solitary waves of uniform depth which interact with a stationary point vortex are considered. Waves both with and without a submerged obstacle are computed. The method of solution is collocation of Bernoulli's equation at a finite number of points on the free surface coupled with equations for equilibrium of a point vortex. The stream function and vortex location are found by computing a conformal map of the flow domain to an infinite strip. For a given obstacle the solutions are parametrized with respect to Froude number and vortex circulation. When no obstacle is present there are two families of solutions, in one of which the amplitude of the wave increases by increasing the circulation, while in the other amplitude increases by decreasing the circulation. Beyond a certain critical Froude number the maximum amplitude wave has a sharp crest with an angle of 120 degrees. Similar behavior is observed for the flow past a submerged obstacle except that there is a critical Froude number below which there is no solution at all.     

The Linearized Crocco Equation

In this paper, we study the existence and uniqueness of a degenerate parabolic equation, with nonhomogeneous boundary conditions, coming from the linearization of the Crocco equation [12]. The Crocco equation is a nonlinear degenerate parabolic equation obtained from the Prandtl equations with the so-called Crocco transformation. The linearized Crocco equation plays a major role in stabilization problems of fluid flows described by the Prandtl equations [5]. To study the infinitesimal generator associated with the adjoint linearized Crocco equation – with homogeneous boundary conditions – we first study degenerate parabolic equations in which the x-variable plays the role of a time variable. This equation is doubly degenerate: the coefficient in front of ∂_x vanishes on a part of the boundary, and the coefficient of the elliptic operator vanishes in another part of the boundary. This makes very delicate the proof of uniqueness of solution. To overcome this difficulty, a uniqueness result is first obtained for an equation in which the elliptic operator is symmetric, and it is next extended to the original equation by combining an iterative process and a fixed point argument (see Th. 4.9). This kind of argument is also used to prove estimates, which cannot be obtained in a classical way.