Instability investigation for rotating thin spherical membrane

A fluid-filled truncated spherical membrane fixed along its truncated edge to a horizontal, rigid and frictionless plane and spinning around a center axis was investigated. A two-parameter Mooney–Rivlin model was used to describe the material of the membrane. The truncated sphere was modeled in 3D using finite element meshes with different symmetry properties. A quadratic function was used for interpolating hydro-static pressure, giving a symmetric tangent stiffness matrix, thereby reproducing the conservative problem. Various problem settings were considered, related to the spinning, and different instability behaviors were observed. Multi-parametric problems were defined, generalized paths including primary and secondary paths were followed. Stability of the multi-parametric problem was evaluated using generalized eigenvalue analysis based on the total differential matrix for the constrained problem. Numerical results showed that mesh symmetry affected the simulated stability behavior. Fold line evaluations showed the parametric effects on critical solutions.     

A stability investigation for an incompressible simple fluid with fading memory

The nonlinear equations of motion for an incompressible simple fluid, occupying a fixed bounded container, are formulated on the basis of the finitelinear viscoelasticity theory for materials with fading memory; this formal boundary-initial value problem is then viewed as a nonlinear abstract evolution equation on a certain Hilbert space. It is shown that a linearized version of this evolution equation is associated with a linear dynamical system on this Hilbert space, and several results for stability and asymptotic behavior for this linearized problem are proved through the use of Liapunov stability methods. On the assumption that the original nonlinear evolution equation also is associated with some dynamical system on the same space, it is shown that the rest condition of the fluid is stable and all motions are bounded. The Liapunov function employed for this purpose can be interpreted as a mechanical energy function for the fluid.E. F. Infante's work was supported in part by the U.S. Office of Naval Research (grant N0014-76-C-0278P002), the U.S. National Science Foundation (grant MCS-76-07247 A03), and the U.S. Army Research Office (grant AROD 31-124-73-G-130); that of J. A. WALKER was supported in part by the U.S. National Science Foundation (grant ENG76-81570) and the U.S. Air Force (grant AFOSR-76-3063A).     

Motion stability dynamics for spacecraft coupled with partially filled liquid container

This paper presents a canonical Hamiltonian model of liquid sloshing for the container coupled with spacecraft. Elliptical shape of rigid body is considered as spacecraft structure. Hamiltonian system is an important form of mechanical system. It mostly used to stabilize the potential shaping of dynamical system. Free surface movement of liquid inside the container is called sloshing. If there is uncontrolled resonance between the motion of tank and liquid-frequency inside the tank then such sloshing can be a reason of attitude disturbance or structural damage of spacecraft. Equivalent mechanical model of simple pendulum or mass attached with spring for sloshing is used by many researchers. Mass attached with spring is used as an equivalent model of sloshing to derive the mathematical equations in terms of Hamiltonian model. Analytical method of Lyapunov function with Casimir energy function is used to find the stability for spacecraft dynamics. Vertical axial rotation is taken as the major axial steady rotation for the moving rigid body.     

Linear and non‐linear stability analysis for finite difference discretizations of high‐order Boussinesq equations

This paper considers a method of lines stability analysis for finite difference discretizations of a recently published Boussinesq method for the study of highly non‐linear and extremely dispersive water waves. The analysis demonstrates the near‐equivalence of classical linear Fourier (von Neumann) techniques with matrix‐based methods for formulations in both one and two horizontal dimensions. The matrix‐based method is also extended to show the local de‐stabilizing effects of the non‐linear terms, as well as the stabilizing effects of numerical dissipation. A comparison of the relative stability of rotational and irrotational formulations in two horizontal dimensions provides evidence that the irrotational formulation has significantly better stability properties when the deep‐water non‐linearity is high, particularly on refined grids. Computation of matrix pseudospectra shows that the system is only moderately non‐normal, suggesting that the eigenvalues are likely suitable for analysis purposes. Numerical experiments demonstrate excellent agreement with the linear analysis, and good qualitative agreement with the local non‐linear analysis. The various methods of analysis combine to provide significant insight into the numerical behaviour of this rather complicated system of non‐linear PDEs. Copyright © 2004 John Wiley & Sons, Ltd.     

Method of large-scale system and analysis of three-axis stability of system partially filled with liquid

In the paper the method of large-scale system is investigated. The criteria of three-axis stability of a partially liquid-filled system are obtained by means of the method of large-scale system. Numerical results are given for a liquid-filled system with three-flywheels.     

Nonlinear axisymmetric bending and stability of thin spherical shallow shell with variable thickness under uniformly distributed loads

In this paper, the nonlinear bending and stability of thin spherical shallow shell with variable thickness under uniformly distributed loads are investigated by a new modified iteration method proposed by Prof. Yeh Kai-yuan and the author^[1]. Deflections and critical loads have been calculated and the numerical results obtained have been given in figures and tabular forms. It is shown that the final equation determining the central deflection and the load obtained coincides with the cusp catastrophe manifold. Projects Supported by the Science Fund of the Chinese Academy of Sciences     

变截面高层筒体结构的有限元线法整体稳定和二阶分析

本文用有限元线法对变截面的高层简体结构进行空间整体稳定和二阶分析。先把实际框筒结构分段连续化为正交各向异性折板结构；用有限元线法，通过有限条元半离散化，取结线上位移为基本未知函数，由势能驻值原理建立稳定和二阶分析的常微分方程组；由常微分方程求解器直接求解。     

Nonlinear stability of Poiseuille flow of a Bingham fluid: theoretical results and comparison with phenomenological criteria

We present new results on the nonlinear stability of Bingham fluid Poiseuille flows in pipes and plane channels. These results show that the critical Reynolds number for transition, Re_c, increases with Bingham number, B, at least as fast as Re_cB^1/2 as B→∞. Estimates for the rate of increase are also provided. We compare these bounds and existing linear stability bounds with predictions from a series of phenomenological criteria for transition, as B→∞, concluding that only Hanks [AIChE J. 9 (1963) 306; 15 (1) (1963) 25] criteria can possibly be compatible with the theoretical criteria as B→∞. In the more practical range of application, 0≤B≤50, we show that there exists a large disparity between the different phenomenological criteria that have been proposed.     

Linear stability diagrams for the shear flow of an Oldroyd-B fluid with slip along the fixed wall

We consider the time-dependent shear flow of an Oldroyd-B fluid with slip along the fixed wall. Slip is allowed by means of a generic slip equation predicting that the shear stress is a non-monotonic function of the velocity at the wall. The complete one-dimensional stability analysis to one-dimensional disturbances is carried out and the corresponding neutral stability diagrams are constructed. Asymptotic results for large values of the elasticity number and finite element calculations are also presented. The instability regimes are within or coincide with the negative-slope regime of the slip equation. The numerical calculations agree with the linear stability results when the size of the initial perturbation is small. Large perturbations may destabilize a linearly stable steady state, leading to a periodic solution. The period and the amplitude of the periodic solutions increase with elasticity. Received: 19 June 1997 Accepted: 22 September 1997     

Analytical investigation of the fluid flow in the interface region between a porous medium and a clear fluid in channels partially filled with a porous medium

In this paper analytical solutions for the steady fully developed laminar fluid flow in the parallel-plate and cylindrical channels partially filled with a porous medium and partially with a clear fluid are presented. The Brinkman-extended Darcy equation is utilized to model the flow in a porous region. The solutions account for the boundary effects and for the stress jump boundary condition at the interface recently suggested by Ochoa-Tapia and Whitaker. The dependence of the velocity on the Darcy number and on the adjustable coefficient in the stress jump boundary condition is investigated. It is shown that accounting for a jump in the shear stress at the interface essentially influences velocity profiles.     

Dynamic stability analysis and DQM for beams with variable cross-section

The present paper deals with the dynamic behaviour of a clamped beam subjected to a sub-tangential follower force at the free end. The aim of this work is to obtain the frequency–axial load relationship for a beam with a variable circular cross-section. In this way, one can identify both divergence critical loads – where the frequency goes to zero – and the flutter critical load – in correspondence with two frequencies coalescence. The numerical approach adopted for solving the partial differential equation of motion is the differential quadrature method (henceforth DQM). This method was proposed by Bellmann and Casti [Bellmann, R.E., Casti, J., 1971. Differential quadrature and long-term integration. J. Math. Anal. 34, 235–238] and has been employed recently in the solution of solid mechanics problems by Bert and Malik [Bert, C.W., Malik, M., 1996. Differential quadrature method in computational mechanics: a review. Appl. Mech. Rev., ASME, 49 (1), 1–28] and Chen et al. [Chen, W., Stritz, A.G., Bert, C.W., 1997. A new approach to the differential quadrature method for fourth-order equations. Int. J. Numer. Method Eng. 40, 1941–1956]. More precisely, a modified version of this method has been used, as proposed by De Rosa and Franciosi [De Rosa, M.A., Franciosi, C., 1998a. On natural boundary conditions and DQM. Mech. Res. Commun. 25 (3), 279–286; De Rosa, M.A., Franciosi, C., 1998b. Non classical boundary conditions and DQM. J. Sound Vibrat. 212(4), 743–748] to satisfy all the boundary conditions.Some frequencies–axial loads relationships are reported in order to show the influence of tapering on the critical loads.     

ON THE BOUNDEDNESS AND THE STABILITY RESULTS FOR THE SOLUTION OF CERTAIN FOURTH ORDER DIFFERENTIAL EQUATIONS VIA THE INTRINSIC METHOD

ＯＮＴＨＥＢＯＵＮＤＥＤＮＥＳＳＡＮＤＴＨＥＳＴＡＢＩＬＩＴＹＲＥＳＵＬＴＳＦＯＲＴＨＥＳＯＬＵＴＩＯＮＯＦＣＥＲＴＡＩＮＦＯＵＲＴＨＯＲＤＥＲＤＩＦＦＥＲＥＮＴＩＡＬＥＱＵＡＴＩＯＮＳＶＩＡＴＨＥＩＮＴＲＩＮＳＩＣＭＥＴＨＯＤＣｅｍｉｌＴＵＮＣ;Ａ...     

The p-moment exponential robust stability for stochastic systems with distributed delays and interval parameters

The p-moment exponential robust stability for stochastic systems with distributed delays and interval parameters is studied.By constructing the LyapunovKrasovskii functional and employing the decomposition technique of interval matrix and Ito's formula,the delay-dependent criteria for the p-moment exponential robust stability are obtained.Numerical examples show the validity and practicality of the presented criteria.     

Analysis of hybrid algorithms for the Navier–Stokes equations with respect to hydrodynamic stability theory

Hybrid three‐dimensional algorithms for the numerical integration of the incompressible Navier–Stokes equations are analyzed with respect to hydrodynamic stability in both linear and nonlinear fields. The computational schemes are mixed—spectral and finite differences—and are applied to the case of the channel flow driven by constant pressure gradient; time marching is handled with the fractional step method. Different formulations—fully explicit convective term, partially and fully implicit viscous term combined with uniform, stretched, staggered and non‐staggered meshes, x‐velocity splitted and non‐splitted in average and perturbation component – are analyzed by monitoring the evolution in time of both small and finite amplitude perturbations of the mean flow. The results in the linear field are compared with correspondent solutions of the Orr–Sommerfeld equation; in the nonlinear field, the comparison is made with results obtained by other authors. Copyright © 2002 John Wiley & Sons, Ltd.     

Nonlinear evolution analysis and stability for convection in a mushy layer with permeable interface

We consider the problem of two- and three-dimensional nonlinear buoyant flows in horizontal mushy layers during the solidification of binary alloys. We study the nonlinear evolution of such flow based on a recently developed realistic model for the mushy layer with permeable interface. The evolution approach is based on a Landau type equation for the amplitude of the secondary nonlinear solution, which can be in the form of rolls, squares, rectangles or hexagons. Using both analytical and computational methods, we calculate the solutions to the evolution equation near the onset of motion for both subcritical and supercritical regimes and determine the stable solutions. We find, in particular, that for several investigated cases with different parameter regimes, secondary solution in the form of subcritical down-hexagons or supercritical up-hexagons can be stable. However, the preferred solution for smallest values of the Rayleigh number and the amplitude of motion is in the form of subcritical down-hexagons. This result appears to agree with the experimental observation on the form of the convective flow near the onset of motion.     

Boundedness and persistence and global asymptotic stability for a class of delay difference equations with higher order

ＩｎｔｒｏｄｕｃｔｉｏｎＣｏｎｓｉｄｅｒｔｈｅｆｏｌｌｏｗｉｎｇｄｅｌａｙｄｉｆｆｅｒｅｎｃｅｅｑｕａｔｉｏｎｗｉｔｈｈｉｇｈｅｒｏｒｄｅｒｘｎ+1=Ａ0ｘｐ0ｎ+ Ａ1ｘｐ1 ｎ- 1+… + Ａｋｘｐｋｎ-ｋ　　 (ｎ=0 ,1 ,2 ,… ) ,( 1 )ｗｈｅｒｅＡｉ,ｐｉ ∈ [0 ,∞ )　　 (ｉ=0 ,1 ,2 ,…ｋ;ｋ∈ 1 ,2 ,… ) ( 2 )ａｎｄｔｈｅｉｎｉｔｉａｌｃｏｎｄｉｔｉｏｎｓｘ-ｋ,ｘ-ｋ+1,… ,ｘ- 1,ｘ0 ａｒｅａｒｂｉｔｒａｒｙｐｏｓｉｔｉｖｅｎｕｍｂｅｒｓ.Ｓｏｍｅｓｐｅｃｉａｌｃａｓｅｓｏｆｋ=1ｆｏｒＥｑ .( 1 )ｈａｖｅｂｅｅｎ…     

A limiter-based well-balanced discontinuous Galerkin method for shallow-water flows with wetting and drying: Triangular grids

A novel wetting and drying treatment for second-order Runge-Kutta discontinuous Galerkin methods solving the nonlinear shallow-water equations is proposed. It is developed for general conforming two-dimensional triangular meshes and utilizes a slope limiting strategy to accurately model inundation. The method features a nondestructive limiter, which concurrently meets the requirements for linear stability and wetting and drying. It further combines existing approaches for positivity preservation and well balancing with an innovative velocity-based limiting of the momentum. This limiting controls spurious velocities in the vicinity of the wet/dry interface. It leads to a computationally stable and robust scheme, even on unstructured grids, and allows for large time steps in combination with explicit time integrators. The scheme comprises only one free parameter, to which it is not sensitive in terms of stability. A number of numerical test cases, ranging from analytical tests to near-realistic laboratory benchmarks, demonstrate the performance of the method for inundation applications. In particular, superlinear convergence, mass conservation, well balancedness, and stability are verified.     

The activation method for discretized conservative nonlinear stability problems with multiple parameter and state variables

For nonlinear stability problems of discretized conservative systems with multiple parameter variables and multiple state variables,the activation method is put forward,by which activated potential functions and activated equilibrium equations are derived.The activation method is the improvement and enhancement of Liapunov-Schmidt method in elastic stability theory.It is more generalized and more normalized than conventional perturbation methods.The activated potential functions may be transformed into normalized catastrophe potential functions.The activated equilibrium equations may be treated as bifurcation equations.The researches in this paper will motivate the combination of elastic stability theory with catastrophe theory and bifurcation theory.     

开口薄壁杆件结构稳定分析的精确单元和两步求解算法

从控制微分方程的通解出发,构造受偏心压力作用开口薄壁杆件的精确形函数,建立用于开口薄壁杆件结构稳定性分析的精确有限元,得到了单元刚度矩阵和几何刚度矩阵的显式表达,提出了计算给定区间内各阶临界荷载以及相应失稳模态的两步计算方法。计算结果表明,与常规单元相比,采用精确单元无需进行网格细分就可以获得精确的数值结果,结合本文的两步求解算法,可以准确获得给定区间内全部临界荷载和失稳模态。