Stability of finite-amplitude and overstable convection of a conducting fluid through fixed porous bed

The stability of infinitestimal steady and oscillatory motions and finite amplitude steady motions of a conducting fluid through porous media with free boundaries which is heated from below and cooled from above is investigated in the presence of a uniform magnetic field. Infinitesimal steady motions are investigated using Liapunov method and its is shown that the principle of exchange of stability is valid only when P_m/P_r≤1 with a restricted value of the Hartmann number. It is shown that overstable motions are due to the zonal current induced by the magnetic field. Finite amplitude steady motions are investigated using Veronis [1] analysis and it is shown that for a restricted range of Hartmann numbers and porous parameter P_l, steady finite-amplitude motions can exist for values of the Rayleigh number smaller than that value corresponding to oscillatory motions. Since the Busse number is greater than the wave number the horizontal scale of the steady finite-amplitude motions is larger than that of the overstable motions.     

Secondary convective motions in a plane vertical layer

Secondary plane-parallel motion in a vertical layer between isothermal planes heated to different temperatures is unstable at low and moderate values of the Prandtl number with respect to monotonically increasing disturbances [1]. The results of numerical experiments carried out by the method of networks [2, 3] indicate that this instability leads to the development of stationary secondary motions; the secondary motions have also been investigated in [4] by averaging the original equations. In the present paper we consider plane and three-dimensional stationary spatially periodic secondary motions near the threshold at which the motions develop. We make use of the methods of branching theory which were used earlier for the investigation of isothermal flows [5–9]. We determine the regions of “soft∝ and “hard∝ instability of the plane-parallel motion and the region of stability of the secondary motions. We give the results obtained by calculation of the basic characteristics of the secondary motions.     

Numerical study of convection of a liquid heated from below

The equilibrium of a liquid heated from below is stable only for small values of the vertical temperature gradient. With increase of the temperature gradient a critical equilibrium situation occurs, as a result of which convection develops. If the liquid fills a closed cavity, then there is a discrete sequence of critical temperature gradients (Rayleigh numbers) for which the equilibrium loses stability with respect to small characteristic disturbances. This sequence of critical gradients and motions may be found from the solution of the linear problem of equilibrium stability relative to small disturbances. If the temperature gradient exceeds the lower critical value, then (for steady-state heating conditions) there is established in the liquid a steady convective motion of a definite amplitude which depends on the magnitude of the temperature gradient. Naturally, the amplitude of the steady convective motion cannot be determined from linear stability theory; to find this amplitude we must solve the problem of convection with heating from below in the nonlinear formulation. A nonlinear study of the steady motion of a liquid in a closed cavity with heating from below was made in [1]. In that study it was shown that for Rayleigh numbers R which are less than the lower critical value R_c steady-state motions of the liquid are not possible. With R>R_c a steady convection arises, whose amplitude near the threshold is small and proportional to (R–R_c)^1/2 (the so-called soft instability)-this is in complete agreement with the results of the phenom-enological theory of Landau [2, 3].Primarily, various versions of the method of expansion in powers of the amplitude [4–8] have been used, and, consequently, the results obtained in those studies are valid only for values of R which are close to R_c, i. e., near the convection threshold.It is apparent that the study of developed convective motion far from the threshold can be carried out only numerically, with the use of digital computers. In [9, 10] the numerical methods have been successfully used for the study of developed convection in an infinite plane horizontal liquid layer.The present paper undertakes the numerical study of plane convective motions of a liquid in a closed cavity of square section. The complete nonlinear system of convection equations is solved by the method of finite differences on a digital computer for various values of the Rayleigh number, the maximal value exceeding by a factor of 40 the minimal critical value R_c. The numerical solution permits following the development of the steady motion which arises with R>R_c in the course of increase of the Rayleigh number and permits study of the oscillatory motions which occur at some value of the parameter R. The heat transfer through the cavity is studied. The corresponding linear problem on equilibrium stability is solved approximately by the Galerkin method.     

Effect of a latitudinal temperature gradient on convective motions arising in a rotating spherical layer

The hydrodynamics of planetary atmospheres and Interiors are frequently directly or indirectly connected with convective motions taking place in rotating liquid spherical layers in the field of a central force. Convective stability in a spherical layer at rest, in a central gravity field, was first discussed in [1, 2]. It was shown that the critical Rayleigh number Rao at which convective instability sets in and the wave number of the critical perturbations depend essentially on the thickness of the layer. As in the plane case, the problem of the convective stability of a spherical layer is found to be degenerate, and the form of the critical perturbations cannot be determined from the linear problem. In actuality, minimization of the Rayleigh number permits establishing only the wave numberl for the spherical harmonic Y _l ^m (θ, ?), realized at the limit of stability; the parameter m remains indeterminate and thus 2l+1 independent convective modes correspond to Rao. In [3] a study was made of the convective stability of a liquid in a slowly rotating thin spherical layer. It was shown that the presence of rotation eliminates the degeneracy; at the limit of stability there arise motions corresponding to the Y _l l (θ, ?) -harmonic with a degenerate maximum at the equator, and propagating in a wave manner toward the side opposite to the rotation. In the present work a study is made of the convective stability of a flow of liquid, arising in a rotating spherical layer due to a nonuniform distribution of the temperatures at one of the boundaries of the layer. In such a statement of the problem it is possible to model large-scale motions in the atmospheres of large planets having internal sources of heat and absorbing solar radiation near the cloud cover of the atmosphere. It is established that, depending on the relationships between the parameters imparting the rotation and the inhomogeneous distribution of the temperature, there is either stabilization or destabilization of the layer in comparison with a fixed layer of the same thickness and with the same, but uniformly distributed heat flux supplied to the layer. A study is made of the form of the corresponding critical perturbations.     

Nonlinear interaction of convective waves and the onset of turbulence in a rotating horizontal layer

Nonlinear development and interaction of disturbances in the classical Bernard-Rayleigh problem has been considered in a number of works [1, 2]. As a rule, the investigation leads to the construction and study of model equations of Landau's type [3–6]. In the present work the convection equation in the Boussinesq approximation is solved with the aid of direct methods [7–10]. We investigate the evolution of the disturbances which are a superposition of two waves with different wave numbers ₁ and ₂. We consider the appearance of harmonics of the form n₁ ± m₂, where n and m are integers. The main attention is given to problems of the onset of turbulence [11]. The numerical experments carried out showed that the consideration of the interaction of a limited number of harmonics (from 23 to 500) allows one to reproduce some characteristic features of the motion during the onset of turbulence.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 9–15, March–April, 1977.     

Bifurcation of the solution to the problem of steady traveling waves in a layer of viscous liquid on an inclined plane

Traveling waves in a viscous liquid flowing down an inclined plane can be described at small and moderate Reynolds numbers by an ordinary differential equation in the thickness of the layer [1, 2]. As the Reynolds number tends to zero, this equation goes over into an equation of third order with quadratic nonlinearity [3]. Periodic solutions of this last equation bifurcating from the plane-parallel solution have been investigated by Nepomnyashchii and Tsvelodub [3–6]. In the present paper, a study is made of the bifurcation of periodic solutions from periodic solutions, namely, an investigation is made of the values of the wave number for which a periodic solution is not unique; a bifurcation equation is derived, the number of bifurcating solutions is found, and their behavior near a bifurcation point is considered; and the bifurcating solutions are continued numerically with respect to a parameter (the wave number) from the neighborhoods of the bifurcation points.     

Instability of Sedimenting Bidisperse Particle Gas Suspensions

A linear stability analysis for a sedimenting bidisperse gas-solid suspension (or gas fluidized bed) is performed. Mass, momentum and energy conservation equations for each of the two species are derived using constitutive equations that are valid at high Stokes numbers, (St₁ 1). The homogeneous suspension becomes unstable at sufficiently large St₁ to waves of particle volume fraction with the wave number in the vertical direction. Numerical calculations of the growth rate in an unbounded suspension indicate that the marginal stability limits are controlled by the small wave number (k 1) behavior. Depending on the Stokes number and the volume fractions ₁ and ₂ of the two species, the suspension becomes unstable due to O(k) or O(k²) contributions to the growth rate. The O(k) term corresponds to an instability due to kinematic waves similar to that predicted for bidisperse suspensions of particles in viscous liquids [22]. The O(k²) contribution represents an instability to dynamic waves similar to that obtained from an analysis of averaged equations for monodisperse fluidized beds [4].     

Investigation of the hypersonic viscous shock layer around blunt bodies in a nonuniform flow

Unseparated viscous gas flow past a body is numerically investigated within the framework of the theory of a thin viscous shock layer [13–15]. The equations of the hypersonic viscous shock layer with generalized Rankine-Hugoniot conditions at the shock wave are solved by a finite-difference method [16] over a broad interval of Reynolds numbers and values of the temperature factor and nonuniformity parameters. Calculation results characterizing the effect of free-stream nonuniformity on the velocity and temperature profiles across the shock layer, the friction and heat transfer coefficients and the shock wave standoff distance are presented. The unseparated flow conditions are investigated and the critical values of the nonuniformity parameter a_k [10] at which reverse-circulatory zones develop on the front of the body are obtained as a function of the Reynolds number. The calculations are compared with the asymptotic solutions [10, 12].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 154–159, May–June, 1987.     

Heat transfer enhancement in grooved channels due to flow bifurcations

The flow bifurcation scenario and heat transfer characteristics in grooved channels, are investigated by direct numerical simulations of the mass, momentum and energy equations, using the spectral element methods for increasing Reynolds numbers in the laminar and transitional regimes. The Eulerian flow characteristics show a transition scenario of two Hopf bifurcations when the flow evolves from a laminar to a time-dependent periodic and then to a quasi-periodic flow. The first Hopf bifurcation occurs to a critical Reynolds number Rec₁ that is significantly lower than the critical Reynolds number for a plane-channel flow. The periodic and quasi-periodic flows are characterized by fundamental frequencies ω₁ and m· ω₁+n·ω₂, respectively, with m and n integers. Friction factor and pumping power evaluations demonstrate that these parameters are much higher than the plane channel values. The time-average mean Nusselt number remains mostly constant in the laminar regime and continuously increases in the transitional regime. The rate of increase of this Nusselt number is higher for a quasi-periodic than for a periodic flow regime. This higher rate originates because better flow mixing develops in quasi-periodic flow regimes. The flow bifurcation scenario occurring in grooved channels is similar to the Ruelle-Takens-Newhouse transition scenario of Eulerian chaos, observed in symmetric and asymmetric wavy channels.     

On periodic motions of an orbital dumbbell-shaped body with a cabin-elevator

The motion of a dumbbell-shaped body (a pair of massive points connected with each other by a weightless rod along which the elevator, i.e., a third point, is moving according to a given law) in an attractive Newtonian central field is considered. In particular, such a mechanical system can be considered as a simplified model of an orbital cable system equipped with an elevator. The practically most interesting case where the cabin performs periodic ??shuttle??motions is studied. Under the assumption that the elevator mass is small compared with the dumbbell mass, the Poincaré theory is used to determine the conditions for the existence of families of system periodic motions analytically depending on the arising small parameter and passing into some stable radial steady-state motion of the unperturbed problem as the small parameter tends to zero. It is also proved that, for sufficiently small parameter values, each of the radial relative equilibria generates exactly one family of such periodic motions. The stability of the obtained periodic solutions is studied in the linear approximation, and these solutions themselves are calculated up to terms of the firstorder in the small parameter. The contemporary studies of the motion of orbital dumbbell systems apparently originated in Okunev??s papers [1, 2]. These studies were continued in [3], where plane motions of an orbit tether (represented as a dumbbell-shaped satellite) in a circular orbit were considered in the satellite approximation. In [4], in the case of equal masses and in the unbounded statement, the energy-momentum method was used to perform the dynamic reduction of the problem and analyze the stability of relative equilibria. A similar technique was used in [5], where, in contrast to the above-mentioned problems, the massive points were connected by an elastic spring resisting to compression and forming a dumbbell with elastic properties. Under such assumptions, the stability of radial configurations was investigated in that paper. The bifurcations and stability of steady-state configurations of a deformable elastic dumbbell were also studied in [6]. Various obstacles arising in the construction of orbital cable systems, in particular, the strong deformability of known materials, were discussed in [7]. In [8], the problem of orbital motion of a pair of massive points connected by an inextensible weightless cable was considered in the exact statement. In other words, it was assumed that a unilateral constraint is imposed on themassive points. The conditions of stability of vertical positions of the relative equilibria of the cable system, which were obtained in [8], can be used for any ratio of the subsatellite and station masses. In turn, these results agree well with the results obtained earlier in the studies of stability of vertical configurations in the case of equal masses of the system end bodies [3, 4]. One of the basic papers in the dynamics of three-body orbital cable systems is the paper [9]. The steady-state motions and their bifurcations and stability were studied depending on the elevator cabin position in [10].     

Stability of a nonuniformly heated fluid in a porous horizontal layer

We investigate the stability of a nonuniformly heated fluid in the gravitational field in a plane horizontal porous layer through which vertical forced motion is effected. A similar system was studied in [1, 2]. In the present paper, the nonuniformity of the permeability of the porous layer with respect to the depth and the dependence of the viscosity of the saturating fluid on the temperature are taken into account in addition. As a result of the application of the linear stability theory, an eigenvalue problem arises, which is solved numerically. A family of curves representing the dependence of the critical modified Rayleigh number Ra _k ^⋆ on the injection parameter (the Péclet number Pe) for different degrees of inhomogeneity of the permeability and the viscosity is obtained. It is found that although Pe=0 corresponds to Ra _k ^⋆ for uniform permeability and viscosity and the stability increases monotonically as Pe increases, the presence of nonuniformity of the permeability or the viscosity leads to the appearance of a stability minimum in the region Pe≈1, while under the simultaneous influence of these two factors, the minimum is shifted into the region Pe≈2. The results of the paper can be used, for example, in the investigation of heat transfer in the case of forced fluid motion in the fissures of a permeable rock mass, when, in the case of pumping through a horizontal fissure, the fluid penetrates vertically across its permeable walls into the stratum. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 3–7, November–December, 1986.     

Self-similar magnetogasdynamic flows accompanied by detonation and combustion waves

Magnetogasdynamic (MGD) flows with detonation waves and combustion fronts have attracted more and more attention in recent years. Intensive heat supply assures such a significant increase in the temperature and pressure behind the heat liberation fronts that the gaseous combustion products become conductive so that the flow map in the electric and magnetic fields can vary substantially as compared with ordinary gasdynamics. In the case of finite gas conductivity, when the magnetic Reynolds numbers R_m are low, the asymptotic laws of detonation wave propagation which either go over into the Chapman-Jouguet (CJ) mode (in a number of cases at a finite distance from the initiation source) or remain overcompressed, have been studied [1]. Stationary flow modes behind detonation waves have been investigated in [2] and the problem of the detonation wave originating at the closed end of the tube emerging in the stationary mode in crossed homogeneous magnetic and electric fields has been examined. Results are presented in this paper of an investigation of one-dimensional self-similar flows caused by piston motion in a hot gas mixture in which a detonation wave or combustion front is propagated. The motion is realized in external electric and magnetic fields which exert a substantial effect on the flow of the conductive combustion products. Domains of application of the governing parameters in which the various flow modes are realized are found by using a qualitative and numerical analysis. The results obtained are used to solve problems about the hypersonic gas flow around a thin wedge in an axial magnetic field.     

Instability of small-amplitude convective flows in a rotating layer with free boundaries

The stability of steady convective flows in a horizontal layer with free boundaries, heated from below and rotating about a vertical axis, is studied in the Boussinesq approximation (Rayleigh-Bénard convection). The flows considered are convective rolls or square cells that are sums of two perpendicular rolls with equal wave numbers k. It is assumed that the Rayleigh number is almost critical in order for convective flows with a wave number k: R = R _c (k) + ε² to arise, the amplitude of the supercritical states being of the order of ε. It is shown that the flows are always unstable relative to perturbations that are the sum of one long-and two short-wave modes corresponding to linear rolls turned through small angles in opposite directions.     

Free vibration analysis of thin rotating cylindrical shells using wave propagation approach

The wave propagation approach is extended to study the frequency characteristics of thin rotating cylindrical shells. Based on Sanders’ shell theory, the governing equations of motion, which take into account the effects of centrifugal and Coriolis forces as well as the initial hoop tension due to rotation, are derived. And, the displacement field is expressed in the form of wave propagation associated with an axial wavenumber k _m and circumferential wavenumber n. Using the wavenumber of an equivalent beam with similar boundary conditions as the cylindrical shell, the axial wavenumber k _m is determined approximately. Then, the relation between the natural frequency with the axial wavenumber and circumferential wavenumber is established, and the traveling wave frequencies corresponding to a certain rotating speed are calculated numerically. To validate the results, comparisons are carried out with some available results of previous studies, and good agreements are observed. Finally, the relative errors induced by the approximation using the axial wavenumber of an equivalent beam are evaluated with respect to different circumferential wavenumbers, length-to-radius ratios as well as thickness-to-radius ratios, and the conditions under which the analysis presented in this paper will be accurate are discussed.     

The near wake of sinusoidal wavy cylinders: Three-dimensional POD analyses

Three-dimensional (3D) proper orthogonal decomposition (POD) analyses are conducted to investigate the near wake of sinusoidal wavy cylinders. For a wave amplitude a/D_m = 0.152, three typical spanwise wavelengths (λ_z) of the wavy cylinder are taken into account, i.e., λ_z/D_m = 1.89, 3.79 and 6.06, where D_m is the mean diameter of the wavy cylinder, among which λ_z/D_m = 1.89 and 6.06 are the optimum wavelengths corresponding to the largest reduction/suppression of fluid forces acting on the wavy cylinder. Time- and space-resolved three-component velocities of the near wake flow, obtained from large eddy simulation (LES) at a subcritical Reynolds number Re = 3 × 10³, are used in the 3D POD analyses. Comparison is made among the wavy cylinders of the three λ_z/D_m values as well as between them and a smooth cylinder, in terms of POD modes, mode energy, mode coefficients, as well as reconstructed flow structures by lower modes. For the optimum λ_z/D_m = 1.89 and 6.06, energy associated with the first two POD modes is significantly reduced compared with that for λ_z/D_m = 3.79 and the smooth cylinder. Distinct characteristics are observed on the lower POD modes for the wavy cylinders. It is found that the first two POD modes for λ_z/D_m = 1.89 and 6.06 are linked to large-scale streamwise vortices that are additionally introduced into the near wake due to the wavy geometry. Meanwhile, POD mode 3 suggests that the wavy cylinder with the larger optimum λ_z/D_m (= 6.06) generates dominant hairpin-like and spanwise coherent structures (CSs) shedding from the saddle at a different frequency from those shedding from the node. Evolutionary development of these CSs is discussed based on reconstructed flows.     

Stability of plane oscillations of a satellite in a circular orbit

We study motions of a rigid body (a satellite) about the center of mass in a central Newtonian gravitational field in a circular orbit. There is a known particular motion of the satellite in which one of its principal central axes of inertia is perpendicular to the orbit plane and the satellite itself exhibits plane pendulum-like oscillations about this axis. Under the assumption that the satellite principal central moments of inertia A, B, and C satisfy the relation B = A + C corresponding to the case of a thin plate, we perform rigorous nonlinear analysis of the orbital stability of this motion.In the plane of the problem parameters, namely, the oscillation amplitude ε and the inertial parameter, there exist countably many domains of orbital stability of the satellite oscillations in the linear approximation. Nonlinear orbital stability analysis was carried out in thirteen of these domains. Isoenergetic reduction of the system of equations of the perturbed motion is performed at the energy level corresponding to the unperturbed periodic motion. Further, using the algorithm developed in [1], we construct the symplectic mapping generated by the equations of the reduced system, normalize it, and analyze the stability. We consider resonance and nonresonance cases. For small values of the oscillation amplitude, we perform analytic investigations; for arbitrary values of ε, numerical analysis is used.Earlier, numerical analysis of stability of plane pendulum-like motions of a satellite in a circular orbit was performed in several special cases in [1–4].     

Secondary vibrational convective motions in a flat vertical layer of liquid

Investigations of the stability of steady-state plane-parallel convective motion between vertical planes heated to different temperatures [1–5] have shown that this motion, depending on the value of the Prandtl number P, exhibits instability of two types. With small and moderate Prandtl numbers, the instability is of a hydrodynamic nature. It is brought about by monotonic perturbations which, in the supercritical region, develop into a periodic, with respect to the vertical, system of steady-state vortices at the interface between the opposing convective flows. Articles [6, 7] are devoted to the numerical investigation of nonlinear secondary steady-state flows. If P>11.4, there appears a new mode of instability, i.e., running thermal waves increasing in the flow; with P>12, this mode becomes more dangerous [4]. This instability is connected with the development of vibrational perturbations, and it can be considered that in the supercritical region the perturbations lead to the establishment of steady-state vibrations. Linear theory has made it possible to determine the boundaries of the regions of stability. In the present article a numerical investigation is made of nonlinear supercritical conditions developing as a result of a loss of stability of the steady-state flow with respect to vibrational perturbations.     

Slip motion and stability of a single degree of freedom elastic system with rate and state dependent friction

We consider quasistatic motion and stability of a single degree of freedom elastic system undergoing frictional slip. The system is represented by a block (slider) slipping at speed V and connected by a spring of stiffness k to a point at which motion is enforced at speed V₀ We adopt rate and state dependent frictional constitutive relations for the slider which describe approximately experimental results of Dieterich and Ruina over a range of slip speeds V. In the simplest relation the friction stress depends additively on a term A In V and a state variable θ; the state variable θ evolves, with a characteristic slip distance, to the value ? B In V, where the constants A, B are assumed to satisfy B > A > 0. Limited results are presented based on a similar friction law using two state variables.Linearized stability analysis predicts constant slip rate motion at V₀ to change from stable to unstable with a decrease in the spring stiffness k below a critical value k_cr. At neutral stability oscillations in slip rate are predicted. A nonlinear analysis of slip motions given here uses the Hopf bifurcation technique, direct determination of phase plane trajectories, Liapunov methods and numerical integration of the equations of motion. Small but finite amplitude limit cycles exist for one value of k, if one state variable is used. With two state variables oscillations exist for a small range of k which undergo period doubling and then lead to apparently chaotic motions as k is decreased.Perturbations from steady sliding are imposed by step changes in the imposed load point motion. Three cases are considered: (1) the load point speed V₀ is suddenly increased; (2) the load point is stopped for some time and then moved again at a constant rate; and (3) the load point displacement suddenly jumps and then stops. In all cases, for all values of k:, sufficiently large perturbations lead to instability. Primary conclusions are: (1) ‘stick-slip’ instability is possible in systems for which steady sliding is stable, and (2) physical manifestation of quasistatic oscillations is sensitive to material properties, stiffness, and the nature and magnitude of load perturbations.     

Wave jumps in media described by a modified Schrödinger equation

The nonlinear Schrödinger equationA _t ±iA _xx+i‖A‖² A=0 describes an envelope of periodic waves with slowly varying parameters on water, in plasmas, and in nonlinear optics [1]. This equation can also be applied to steady periodic waves (the wave amplitude and wave number do not depend on time, the variablest andx are replaced by the variables of a horizontal coordinate system on the surface of the fluid [2]). In the present paper the properties of a modified Schrödinger equation involving the third and higher derivatives are studied. Solutions describing transition regions between uniform wave states are obtained numerically. If the structure of the transition region whose extent increases with time is not considered, these solutions may be interpreted as jumps.     

Unsteady motion of a profile near a wavy shield

The investigation of the motion of a profile near a plane shield and near a wavy shield is one of the ways of approximately taking into account the effect of the surface of a fluid on the characteristics of surface wings. The steady-state characteristics of the profiles near a plane shield have been investigated in a number of theoretical and experimental studies [1, 2]. In the present study we consider the unsteady motion of a thin profile near plane and wavy shields in an ideal incompressible fluid at rest at infinity, and when there is a fluid flow.Translated from Izvestiya Akadeinii Nauk SSSR. Mekhanika Zhidkosti i Gaza, No. 3, pp. 10–16, May–June, 1973.