Über die instabilität konservativer systeme mit gyroskopischen kräften

Summary Routh's theorem states that a steady motion of a discrete, conservative mechanical system is stable if the dynamic potential W(q)=U(q)–T₀(q) assumes a minimum. This is a generalized version of the theorem on the stability of equilibrium at a minimum of the potential energy, which is due to Dirichlet. It is well known that a steady motion may also be stable if W(q) assumes a maximum instead of a minimum. The stability is then due to the gyroscopic terms in the equations of motion, without which the steady motion would be unstable. Here it is shown that the steady motion is always unstable if not only W(q) but also H ₀(q) assumes a maximum, H ₀(q) being the part of the Hamiltonian that does not depend on the momenta. It is astonishing that this unexpectedly simple criterion was not found before now. In the proof, a variational formulation is used for the problem, and the instability is shown directly from the existence of certain motions which diverge from the trivial solution.

---

---

Vorgelegt von C. Truesdell     

On pressure and thermal flux in gas-mixture flows and in two-phase flows

To begin with, two different definitions of pressure, thermal flux, etc. in the diffusion model and two-fluid model are given. Then the physical interpretations of the pressure and the thermal flux are provided by introducing the momentum and energy fluxes,M and ε, through a surface dS in the flow field. The quantities defined in the diffusion model are suggested when the motion of the mixture is studied as a whole, while the quantities defined in the two-fluid model are suggested when the motion of individual species is studied. The collision pressure and thermal flux in dense gas-mixtures are also discussed in detail, i.e. their origin, their expressions in the momentum and energy equations, and their distinctions from the normal partial pressure and thermal flux. A gas-particle flow can be treated as a flow of dense gas-mixtures. The long-standing controversy whether the “inertial coupling term” should exist in the momentum equation can be clarified by the two different definitions of pressure.     

Free Convection Boundary Layer Flow from a Heated Upward Facing Horizontal Flat Plate Embedded in a Porous Medium Filled by a Nanofluid with Convective Boundary Condition

The steady laminar incompressible free convective flow of a nanofluid over a permeable upward facing horizontal plate located in porous medium taking into account the thermal convective boundary condition is studied numerically. The nanofluid model used involves the effect of Brownian motion and the thermophoresis. Using similarity transformations the continuity, the momentum, the energy, and the nanoparticle volume fraction equations are transformed into a set of coupled similarity equations, before being solved numerically, by an implicit finite difference numerical method. Our analysis reveals that for a true similarity solution, the convective heat transfer coefficient related with the hot fluid and the mass transfer velocity must be proportional to x ^−2/3, where x is the horizontal distance along the plate from the origin. Effects of the various parameters on the dimensionless longitudinal velocity, the temperature, the nanoparticle volume fraction, as well as on the rate of heat transfer and the rate of nanoparticle volume fraction have been presented graphically and discussed. It is found that Lewis number, the Brownian motion, and the convective heat transfer parameters increase the heat transfer rate whilst the thermophoresis decreases the heat transfer rate. It is also found that Lewis number and the convective heat transfer parameter enhance the nanoparticle volume fraction rate whilst the thermophoresis parameter decreases nanoparticle volume fraction rate. A very good agreement is found between numerical results of the present article for special case and published results. This close agreement supports the validity of our analysis and the accuracy of the numerical computations.     

Model of steady motion of the interface in a layer of a strongly superheated liquid

Steady propagation of the boundary of a vapor cavity in a layer of a metastable liquid along the heater surface is considered. The temperature and velocity of interface propagation are determined from the equations of conservation of mass, momentum, and energy in the neighborhood of the stagnation point of the vapor cavity and the condition of stability of steady motion of the interface. It is shown that a solution of these equations exists only if the liquid is heated above a threshold value. The calculated velocity of interface motion and the threshold value of temperature are in reasonable agreement with available experimental data for various liquids within wide ranges of saturation pressures and temperatures of the superheated liquid. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 2, pp. 47–55, March–April, 2008.     

Global Nonlinear Stability for Steady Ideal Fluid Flow in Bounded Planar Domains

We prove stability of steady flows of an ideal fluid in a bounded, simply connected, planar region, that are strict maximisers or minimisers of kinetic energy on an isovortical surface. The proof uses conservation of energy and transport of vorticity for solutions of the vorticity equation with initial data in L^p for p>4/3. A related stability theorem using conservation of angular momentum in a circular domain is also proved.     

A momentum coupling method for the unsteady incompressible Navier–Stokes equations on the staggered grid

A new numerical method is developed to efficiently solve the unsteady incompressible Navier–Stokes equations with second-order accuracy in time and space. In contrast to the SIMPLE algorithms, the present formulation directly solves the discrete x- and y-momentum equations in a coupled form. It is found that the present implicit formulation retrieves some cross convection terms overlooked by the conventional iterative methods, which contribute to accuracy and fast convergence. The finite volume method is applied on the fully staggered grid to solve the vector-form momentum equations. The preconditioned conjugate gradient squared method (PCGS) has proved very efficient in solving the associate linearized large, sparse block-matrix system. Comparison with the SIMPLE algorithm has indicated that the present momentum coupling method is fast and robust in solving unsteady as well as steady viscous flow problems. © 1998 John Wiley & Sons, Ltd.     

Stability analysis of steady rotations of a rigid body bearing two-degree-of-freedom control moment gyros with dissipation in gimbal suspension axes

To study the stability of steady rotations of a control moment gyro system with internal dissipation, we use the Barbashin-Krasovskii theorem and the relation, established in [1], between the Lyapunov function and steady motions. Taking into account the special properties of the original problem, we reduce it to a lower-dimensional problem.We give a detailed presentation of an algorithm for analyzing the stability of steady motions of a gyrostat and use this algorithm to perform a complete study for two systems consisting, respectively, of one and two gyros whose gimbal axes are parallel to the principal axis of inertia of the system. Each steady motion is identified as either asymptotically stable or unstable. We find periodic motions that exist only in the presence of dynamic symmetry and which are regular precessions. For the system with two gyros, we prove the asymptotic stability of quiescent states and prove that in the angular momentum range where these states are defined the system does not have any other stable motions.     

Performance of a multigrid calculation procedure in three-dimensional sudden expansion flows

The performance of a recently developed calculation procedure for steady incompressible flows is assessed in a variety of three-dimensional sudden expansion type flows representative of those encountered in several types of industrial equipment. The calculation procedure, called here BLIMM (for block-implicit multigrid method), is based on a coupled solution of the three-dimensional momentum and continuity equations in primitive variables, using the multigrid technique. Different Reynolds numbers and finite difference grids are considered for each flow situation. The rates of convergence and the computational times are reported for each case.     

A solar system model with dissipation

A problem of motion for an arbitrary number of planets is discussed with consideration of the forces of gravitational interaction according to the law of universal gravitation. The planets are assumed to be homogeneous viscoelastic spheres. In the process of motion, the planets are deformed and the dissipation of energy takes place due to internal viscous forces. On the basis of the motion separation method, an approximate system of equations is obtained to describe the motion of planet centers of mass and the variation of planet angular momenta with respect to the centers of mass. The equations of motion contain small conservative corrections to the law of universal gravitation and small dissipative forces whose influence causes a decrease of the total mechanical energy. The motion under consideration admits the following first integral: the law of angular momentum conservation for the system with respect to the centers of mass. When the system executes the steady motion corresponding to its rotation with a constant angular velocity as a rigid body, the dissipative forces do not perform work, since the deformed planets have no time-dependent deformations.     

Inertial steady flow of a bed of granular material down a surface with microrelief with allowance for finite time of intergranular contacting

An inertial flow of a granular material can be described by the laws of conservation of mass, momentum, and energy of random motion of solid particles by invoking some closing relations. In this work, these closing relations are inferred from the dimensional theory. The system of equations obtained is used to determine characteristics of a steady flow of a bed of a granular material down an inclined surface with a microrelief for various Richardson numbers and finite contact times of the particles during their collisions. Novosibirsk Military Institute, Novosibirsk 630103. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 6, pp. 128–132, November–December, 1999.     

Effect of interface deformability on thermocapillary motion of a drop in a tube

The effect of an externally imposed axial temperature gradient on the mobility and deformation of a drop in an otherwise stagnant liquid within an insulated cylindrical tube is investigated. In the absence of bulk transport of momentum and energy, the boundary integral technique is used to obtain the flow and temperature fields inside and outside the deformable drop. The steady drop shapes and the corresponding migration velocities are examined over a wide range of the dimensionless parameters. The steady drop shape is nearly spherical for dimensionless drop sizes <0.5, but becomes slightly elongated in the axial direction for drop sizes comparable to tube diameter. The adverse effect of drop deformation on the effective temperature gradient driving the motion is slightly more pronounced than its favorable effect of reducing drag, thereby leading to a slight reduction in drop mobility with increasing drop deformation. Increasing the viscosity ratio reduces drop deformation and leads to a slight enhancement in the relative mobility (with respect to free thermocapillary motion) of confined drops. When the drop fluid has a lower thermal conductivity than the exterior phase, the presence of the thermally-insulating wall increases the thermal driving force for drop motion (compared to that for the same drop in unbounded domain) by causing more pronounced bending of the isotherms toward the drop. However, the favorable thermal effect of the confining wall is overwhelmed by its retarding hydrodynamic effect, causing the confined drop to always move slower than its unbounded counterpart regardless of the value of the thermal conductivity ratio.     

Some Primitive Concepts in Continuum Mechanics Regarded in Terms of Objective Space-Time Molecular Averaging: The Key Rôle Played by Inertial Observers

In Continuum Mechanics the notions of body, material point, and motion, are primitive. Here these concepts are derived for any (possibly time-dependent) material system via mass and momentum densities whose values are local spacetime averages of molecular quantities. The averaging procedure necessary to ensure molecular-based densities can be agreed upon by all observers (that is, are objective) has implications for constitutive relations. Specifically, such relations should first be expressed in terms of Galilean-invariant functions of the motion relative to an inertial frame. Thereafter such relations can be re-phrased for general observers, thereby yielding general-frame constitutive relations compatible with material frame-indifference. Two postulates concerning observer agreement (which together constitute a statement of material frame-indifference) are shown to imply that any stress response function which is assumed to depend upon the motion in an inertial (general) frame must be Galilean-invariant (invariant under superposed rigid body motions). Accordingly, invariance under superposed rigid body motions is not a fundamental tenet of continuum physics, but rather a consequence of material frame-indifference whenever constitutive dependence upon motion in a general observer frame is postulated.     

A numerical study of turbulent line puffs via the renormalization group (RNG) k–ϵ model

The time evolution of a line puff, a turbulent non-buoyant element with significant momentum, is studied using the renormalization group (RNG) k–ϵ model. The numerical results show that the puff motion is characterized by a vortex pair flow; the computed flow details and scalar mixing characteristics can be described by self-similar relations beyond a dimensionless time of around 30. The added mass coefficient of the puff motion is found to be approximately unity. The predicted puff flow and mixing rate are substantially similar to those obtained from the standard k–ϵ model and are well supported by experimental data. The computed scalar field reveals significant secondary concentration peaks trailing behind in the wake of the puff. The present results suggest that the overall mixing rate of a puff is primarily determined by the large-scale motion and that streamline curvature probably plays a minor role. © 1998 John Wiley & Sons, Ltd.     

一类刚-梁耦合系统平面稳态运动的稳定性

研究了中心力场中的一类刚-弹耦合系统的平面运动动力学，模型是带有一悬臂 梁的刚体. 综合考虑了系统轨道运动与姿态运动，在Lagrange力学体系下给出了系统的运 动方程，在保守系统和考虑梁的材料黏滞阻尼两种情况下，利用能量-动量方法给出了一类 相对平衡点稳定性的充分条件.     

A property of the set of Rivlin-Ericksen tensors

It is shown that if a motion is such that the first two Rivlin-Ericksen tensors A ₁ and A ₂ have the same unit proper vectors, and A₁ has distinct proper numbers, then all the Rivlin-Ericksen tensors A ₁, A ₂, ..., A _n, ... have the same proper vectors. If further this motion is steady and isochoric, then the product of velocity components corresponding to the unit proper vectors, assumed to be non-vanishing, and the product of the abnormalities of the vector fields of the proper vectors, also assumed to be non-vanishing, each bear a constant value along a streamline. If A ₁ has a pair of equal proper numbers the first statement remains true, but the results for steady isochoric motions will not necessarily hold.     

A manufactured solution for a two-dimensional steady wall-bounded incompressible turbulent flow

This paper presents a manufactured solution (MS), resembling a two-dimensional, steady, wall-bounded, incompressible, turbulent flow for RANS codes verification. The specified flow field satisfies mass conservation, but requires additional source terms in the momentum equations. To also allow verification of the correct implementation of the turbulence models transport equations, the proposed MS exhibits most features of a true near-wall turbulent flow. The model is suited for testing six eddy-viscosity turbulence models: the one-equation models of Spalart and Allmaras and Menter; the standard two-equation k–ε model and the low-Reynolds version proposed by Chien; the TNT and BSL versions of the k–ω model.     

Mass,momentum and total excess energy transported by a weak planeN wave

A weakly nonlinear plane acoustic wave is emitted into an ideal gas of semi-infinite extent from an infinite plate by its sinusoidal motion of single period. The wave develops into anN wave in the far field, as long as the energy dissipation is negligible everywhere except for discontinuous shock fronts. The third-order effects at shock fronts are evaluated, i.e., the generation of reflected acoustic wave as a result of the interaction of shock and expansion wave and the production of entropy by the energy dissipation at shock fronts. Consideration of these effects enables one to estimate the whole mass, momentum and total excess energy (sum of the kinetic energy and excess of internal energy over an initial undisturbed value) transported by theN wave to the accuracy of third order of wave amplitude. It is shown that the mass and total excess energy transported by theN wave increase and the momentum decreases to asymptotic limits as the wave propagates. The result shows good agreement with a numerical result obtained by solving the Euler equations with a high-resolution TVD upwind scheme.     

Thermomechanical instability in steady operating regime of continuous-flow chemical reactor with fixed catalyzer bed

A continuous-flow chemical reactor with fixed bed is a vessel filled with a porous catalyzer through which a liquid or gaseous reagent mixture is filtered. The motion of the mixture in this system is maintained by the pressure differential which compensates for the hydraulic resistance.The complete system of equations describing the mass-and heat-transfer processes in the continuous-flow chemical reactor includes the mass, momentum, and energy conservation equations [1], Usually, in the study of the existence and stability of steady reactor operating regimes, very simple reactor models are analyzed, in which only the mass-conservation equation in the form of the diffusion equation and the energy-conservation equation in the form of the heat-conduction equation are taken into account. The equation of motion drops out of consideration, since the velocity of the reagent mixture with the reaction products through the reactor is considered given and unvarying for disturbances of the steady operating regime.The purpose of the present study is to indicate the possibility of instability in the steady-process regime of a chemical reactor, which is associated with the fact that any temperature (and also concentration) variation within the reactor with disturbances of the steady regime will, generally speaking, lead to variation of the mixture viscosity and, consequently, to variation of the hydraulic resistance and the filtration velocity. This interconnection is an additional factor which may have a significant effect on the behavior of disturbances of the steady regime.In contrast with thermal, kinetic, and thermokinetic instability [2], it is natural to refer to the instability being investigated as thermomechanical instability. Let us consider an example.     

3D pressure imaging of an aircraft propeller blade-tip flow by phase-locked stereoscopic PIV

The flow field at the tip region of a scaled DHC Beaver aircraft propeller, running at transonic speed, has been investigated by means of a multi-plane stereoscopic particle image velocimetry setup. Velocity fields, phase-locked with the blade rotational motion, are acquired across several planes perpendicular to the blade axis and merged to form a 3D measurement volume. Transonic conditions have been reached at the tip region, with a revolution frequency of 19,800 rpm and a relative free-stream Mach number of 0.73 at the tip. The pressure field and the surface pressure distribution are inferred from the 3D velocity data through integration of the momentum Navier-Stokes equation in differential form, allowing for the simultaneous flow visualization and the aerodynamic loads computation, with respect to a reference frame moving with the blade. The momentum and pressure data are further integrated by means of a contour-approach to yield the aerodynamic sectional force components as well as the blade torsional moment. A steady Reynolds averaged Navier-Stokes numerical simulation of the entire propeller model has been used for comparison to the measurement data.