On disturbance evolution and nonlinear stability of the maltilayer quasi-geostrophic basic flows

The status of disturbances of both initial values and parameters in the model is further investigated, and the exact explicit estimates on the disturbance energy and disturbance potential enstrophy are given, while the initial disturbance fields depend only on the initial disturbance potential enstrophy, initial disturbance velocity circulation along the boundary, and disturbance parameters. And based on the generalized nonlinear stability concept in the sense of Lyapunov, the nonlinear stability criteria paralleling to Arnold’s second theorem are obtained. Project supported by the National Natural Science Foundation of China (Grant No. 49776286).     

Baroclinic nonlinear exchanges of energy and potential enstrophy in the quasi-geostrophic two-layer model

Abstract

Merilees and Warn's (1975) nonlinear interaction analysis of two-dimensional nondivergent flow is extended to examine the quasi-geostrophic two-layer model. Two sets of triads exist in this model (Salmon, 1978). The purely barotropic triads are the same as the triads examined by Merilees and Warn. Baroclinic-barotropic triads are found to exchange more energy or potential enstrophy with smaller or larger scales depending on the scale of motion as compared with the internal Rossby deformation radius and the relative wavenumber position of baroclinic and barotropic components.     

Nonlinear stability of baroclinic flows over topography

Abstract

A new nonlinear stability criterion is derived for baroclinic flows over topography in spherical geometry. The stability of a wide class of exact three-dimensional nonlinear steady state solutions subject to arbitrary disturbances is established. The resonance condition, at the highest total wavenumber, for the steady state solutions and the stability criteria for baroclinic flow in the absence of topography provide the boundaries of the regions of stability in the presence of topography. The analogous results for flow on periodic or infinite beta planes incorporating non-orthogonal function large scale flows are also discussed.     

A one-dimensional turbulence model; Enstrophy cascades and small-scale intermittency in decaying turbulence

Abstract

Severe unidirectional Fourier truncation of the equations for 2-D incompressible flow leads to a system of three coupled PDEs in one space dimension with the same quadratic invariants as the original set (i.e. energy and enstrophy). Numerically generated equilibria for inviscid, truncated versions of the reduced system are well approximated by Kraichnan's energy-enstrophy equipartition spectra. Viscous calculations for decaying turbulence at moderate resolution (1024 degrees of freedom) also appear to be consistent with a direct, k ^?3, enstrophy cascading inertial range when the dissipation is small. Dissipation range intermittency in the form of spatially intermittent enstrophy dissipation with occasional strong bursts producing linear phase locking is also observed. In contrast to full 2-D simulations, no tendency towards the emergence of isolated, coherent vorticity structures is observed. The model consequently mimics some, but not all, of the properties of the full 2-D set.     

The amplitude of turbulent shear flow

A new formulation of the problem of the statistical stability of fully turbulent shear flow is proposed, in which one seeks mean fields that bound the observed flow from the stable side. In the spirit of maximum transport theory, this formulation admits a larger set of “flows” than are dynamically possible. A sequence of constraints derived from the equations of motion can narrow this set, permitting at each step the determination of a “most stable” field free of any empirical elements. Turbulent channel flow is proposed as the first application and test of this quantitative theory. Past deductive theories for this flow, from “mean field” to “transport upper bounds,” are assessed. It is shown why these theories do not retain the significant destabilizing mechanisms of the actual flow. The implications for turbulent flow of recent work on the nonlinear and three-dimensional instability of laminar shearing flow are described. In first exploration of the “decoupled mean” stability theory proposed here, approximate analytical and numerical stability methods are used to find an amplitude and structure for the averaged flow propoerties. The quantitative results differ by considerably less than two from the observed values, providing an incentive for a more complete numerical study and for further constraints on the admitted class of flows. In the language now current for nonlinear stability theory, evidence is advanced here that anN-dimensional central manifold is adjacent to the realized turbulent flow, whereN has the largest possible value compatible with the dynamical relations.     

Seismic stability of concrete gravity dams strengthened by rockfill buttressing

Rockfill buttressing resting on the downstream face of masonry or concrete gravity dam is often considered as a strengthening method to improve the stability of existing dam for hydrostatic and seismic loads. Simplified methods for seismic stability analysis of composite concrete-rockfill dams are discussed. Numerical analyses are performed using a nonlinear rockfill model and nonlinear dam-rockfill interface behavior to investigate the effects of backfill on dynamic response of composite dams. A typical 35 m concrete gravity dam, strengthened by rockfill buttressing is considered. The results of analyses confirm that backfill can improve the seismic stability of gravity dams by exerting pressure on the dam in opposition to hydrostatic loads. According to numerical analyses results, the backfill pressures vary during earthquake base excitations and the inertia forces of the backfill are the main source for those variations. It is also shown that significant passive (or active) pressure cannot develop in composite dams with a finite backfill width. A simplified model is also proposed for dynamic analysis of composite dam by replacing the backfill with by a series of vertical cantilever shear beams connected to each other and to the dam by flexible links.     

The nonlinear development of disturbances in a zonal shear flow

Abstract

A study is made of the nonlinear stability of a weakly supercritical zonal shear flow in the β-plane approximation. The dynamics of initially small disturbances are examined. The main nonlinear effects are associated with the rearrangement of the critical layer. It is shown that as the wave grows in amplitude, linear regimes of the critical layer (viscous and nonstationary) change over to a nonlinear regime while the exponential law of disturbance growth becomes a power-law.     

Chaotic mixing of tracer and vorticity by modulated travelling Rossby waves

Abstract

We consider the mixing of passive tracers and vorticity by temporally fluctuating large scale flows in two dimensions. In analyzing this problem, we employ modern developments stemming from properties of Hamiltonian chaos in the particle trajectories; these developments generally come under the heading “chaotic advection” or “Lagrangian turbulence.” A review of the salient properties of this kind of mixing, and the mathematics used to analyze it, is presented in the context of passive tracer mixing by a vacillating barotropic Rossby wave. We then take up the characterization of subtler aspects of the mixing. It is shown the chaotic advection produces very nonlocal mixing which cannot be represented by eddy diffusivity. Also, the power spectrum of the tracer field is found to be k ^{? l} at shortwaves—precisely as for mixing by homogeneous, isotropic two dimensional turbulence,—even though the physics of the present case is very different. We have produced two independent arguments accounting for this behavior.

We then examine integrations of the unforced barotropic vorticity equation with initial conditions chosen to give a large scale streamline geometry similar to that analyzed in the passive case. It is found that vorticity mixing proceeds along lines similar to passive tracer mixing. Broad regions of homogenized vorticity ultimately surround the separatrices of the large scale streamline pattern, with vorticity gradients limited to nonchaotic regions (regions of tori) in the corresponding passive problem.

Vorticity in the chaotic zone takes the form of an arrangement of strands which become progressively finer in scale and progressively more densely packed; this process transfers enstrophy to small scales. Although the enstrophy cascade is entirely controlled by the large scale wave, the shortwave enstrophy spectrum ultimately takes on the classical k ^{? l} form. If one accepts that the enstrophy cascade is indeed mediated by chaotic advection, this is the expected behavior. The extreme form of nonlocality (in wavenumber space) manifest in this example casts some doubt on the traditional picture of enstrophy cascade in the Atmosphere, which is based on homogeneous two dimensional turbulence theory. We advance the conjecture that these transfers are in large measure attributable to large scale, low frequency, planetary waves.

Upscale energy transfers amplifying the large scale wave do indeed occur in the course of the above-described process. However, the energy transfer is complete long before vorticity mixing has gotten very far, and therefore has little to do with chaotic advection. In this sense, the vorticity involved in the enstrophy cascade is “fossil vorticity,” which has already given up its energy to the large scale.

We conclude with some speculations concerning statistical mechanics of two dimensional flow, prompted by our finding that flows with identical initial energy and enstrophy can culminate in very different final states. We also outline prospects for further applications of chaotic mixing in atmospheric problems.     

Conservation laws of wave action and potential enstrophy for rossby waves in a stratified atmosphere

The purpose of this article is to discuss the evolution of wave energy, enstrophy and action for atmospheric Rossby waves in a variable mean flow. The presentation is theoretical, but does not represent original research; rather, it is pedagogic in nature. The work of a number of people has been drawn together into a unified account, with much of the algebra implicit in previous work made explicit here. The central results are that wave energy is conserved only when there are no spatial variations in the mean flow, and wave action is conserved even in the presence of such variations as long as they are not in the longitudinal direction. Finally, wave enstrophy is conserved in the presence of arbitrary (slow) mean flow variations.     

A nonlinear stability analysis of convection in a generalized incompressible fluid

Abstract

We establish a nonlinear stability result for convection in a generalized incompressible fluid. Both numerical calculations and an asymptotic analysis are carried out. The linear and nonlinear results are shown to be very close in both cases, implying that the region of possible subcritical instabilities is very small.

During this work I was supported by a research studentship awarded by the Science and Engineering Council of the United Kingdom.     

Energy-consistent integration method and its application to hybrid testing

The time-integration algorithm is an indispensable element to determine response of the boundary of the numerical as well as physical parts in a hybrid test. Instability of the time-integration algorithm may directly lead to failure of the test, so stability of an integration algorithm is particularly important for hybrid testing. The explicit algorithms are very popular in hybrid testing, because iteration is not needed. Many unconditionally stable explicit-algorithms have been proposed for hybrid testing. However, the stability analysis approaches used in all these methods are valid only for linear systems. In this paper, a uniform formulation for energy-consistent time integrations, which are unconditionally stable, is proposed for nonlinear systems. The solvability and accuracy are analyzed for typical energy-consistent algorithms. Some numerical examples and the results of a hybrid test are provided to validate the effectiveness of energy-consistent algorithms.     

Nonlinear evaluations of unconditionally stable explicit algorithms

Two explicit integration algorithms with unconditional stability for linear elastic systems have been successfully developed for pseudodynamic testing.Their numerical properties in the solution of a linear elastic system have been well explored and their applications to the pseudodynamic testing of a nonlinear system have been shown to be feasible. However,their numerical properties in the solution of a nonlinear system are not apparent.Therefore,the performance of both algorithms for use in the solution...     

Extensions of nonlinear error propagation analysis for explicit pseudodynamic testing

Two important extensions of a technique to perform a nonlinear error propagation analysis for an explicit pseudodynamic algorithm (Chang, 2003) are presented. One extends the stability study from a given time step to a complete step-by-step integration procedure. It is analytically proven that ensuring stability conditions in each time step leads to a stable computation of the entire step-by-step integration procedure. The other extension shows that the nonlinear error propagation results, which are derived for a nonlinear single degree of freedom (SDOF) system, can be applied to a nonlinear multiple degree of freedom (MDOF) system. This application is dependent upon the determination of the natural frequencies of the system in each time step, since all the numerical properties and error propagation properties in the time step are closely related to these frequencies. The results are derived from the step degree of nonlinearity. An instantaneous degree of nonlinearity is introduced to replace the step degree of nonlinearity and is shown to be easier to use in practice. The extensions can be also applied to the results derived from a SDOF system based on the instantaneous degree of nonlinearity, and hence a time step might be appropriately chosen to perform a pseudodynamic test prior to testing.     

Hydrodynamical Modeling Of Oceanic Vortices

Mesoscale coherent vortices are numerous in the ocean.Though they possess various structures in temperature and salinity,they are all long-lived, fairly intense and mostly circular. Thephysical variable which best describes the rotation and the density anomaly associated with coherent vortices is potential vorticity. It is diagnostically related to velocity and pressure, when the vortex is stationary. Stationary vortices can be monopolar (circular or elliptical) or multipolar; their stability analysis shows thattransitions between the various stationary shapes are possible when they become unstable. But stable vortices can also undergo unsteady evolutions when perturbed by environmental effects, likelarge-scale shear or strain fields, -effect or topography. Changes in vortex shapes can also result from vortex interactions. such as the pairing, merger or vertical alignment of two vortices, which depend on their relative polarities and depths. Such interactions transfer energy and enstrophy between scales, and are essential in two-dimensional and in geostrophic turbulence. Finally, in relation with the observations, we describe a few mechanisms of vortex generation.     

Analysis of implicit HHT‐α integration algorithm for real‐time hybrid simulation

Real‐time hybrid simulation is a viable experiment technique to evaluate the performance of structures equipped with rate‐dependent seismic devices when subject to dynamic loading. The integration algorithm used to solve the equations of motion has to be stable and accurate to achieve a successful real‐time hybrid simulation. The implicit HHT α‐algorithm is a popular integration algorithm for conducting structural dynamic time history analysis because of its desirable properties of unconditional stability for linear elastic structures and controllable numerical damping for high frequencies. The implicit form of the algorithm, however, requires iterations for nonlinear structures, which is undesirable for real‐time hybrid simulation. Consequently, the HHT α‐algorithm has been implemented for real‐time hybrid simulation using a fixed number of substep iterations. The resulting HHT α‐algorithm with a fixed number of substep iterations is believed to be unconditionally stable for linear elastic structures, but research on its stability and accuracy for nonlinear structures is quite limited. In this paper, a discrete transfer function approach is utilized to analyze the HHT α‐algorithm with a fixed number of substep iterations. The algorithm is shown to be unconditionally stable for linear elastic structures, but only conditionally stable for nonlinear softening or hardening structures. The equivalent damping of the algorithm is shown to be almost the same as that of the original HHT α‐algorithm, while the period elongation varies depending on the structural nonlinearity and the size of the integration time‐step. A modified form of the algorithm is proposed to improve its stability for use in nonlinear structures. The stability of the modified algorithm is demonstrated to be enhanced and have an accuracy that is comparable to that of the existing HHT α‐algorithm with a fixed number of substep iterations. Both numerical and real‐time hybrid simulations are conducted to verify the modified algorithm. The experimental results demonstrate the effectiveness of the modified algorithm for real‐time testing. Copyright © 2011 John Wiley & Sons, Ltd.     

Stability and identification for rational approximation of frequency response function of unbounded soil

Exact representation of unbounded soil contains the single output–single input relationship between force and displacement in the physical or transformed space. This relationship is a global convolution integral in the time domain. Rational approximation to its frequency response function (frequency‐domain convolution kernel) in the frequency domain, which is then realized into the time domain as a lumped‐parameter model or recursive formula, is an effective method to obtain the temporally local representation of unbounded soil. Stability and identification for the rational approximation are studied in this paper. A necessary and sufficient stability condition is presented based on the stability theory of linear system. A parameter identification method is further developed by directly solving a nonlinear least‐squares fitting problem using the hybrid genetic‐simplex optimization algorithm, in which the proposed stability condition as constraint is enforced by the penalty function method. The stability is thus guaranteed a priori. The infrequent and undesirable resonance phenomenon in stable system is also discussed. The proposed stability condition and identification method are verified by several dynamic soil–structure‐interaction examples. Copyright © 2009 John Wiley & Sons, Ltd.     

几何缺陷对拱结构动力稳定性的影响

分析了外激励下几何缺陷对拱结构动力稳定性的影响。推导了拱结构边界确定而结构本身节点坐标偏差随机且指数相关时的条件相关矩阵,分解得到几何缺陷的分布方式和大小。从非线性运动方程出发,分别得出了周期荷载作用下非线性刚度矩阵可线性化,非周期荷载作用下同时考虑几何、材料非线性的Lyapunov指数计算方法。最后以一圆弧拱为例分别对周期荷载、阶跃荷载、脉冲荷载及地震荷载作用下几何缺陷的影响进行了数值分析。结果表明周期激励作用下拱结构存在动力失稳频域;在不同分布方式几何缺陷中动力稳定性对与屈曲模态相似的缺陷最为敏感。     

A Rosenbrock‐W method for real‐time dynamic substructuring and pseudo‐dynamic testing

A variant of the Rosenbrock‐W integration method is proposed for real‐time dynamic substructuring and pseudo‐dynamic testing. In this variant, an approximation of the Jacobian matrix that accounts for the properties of both the physical and numerical substructures is used throughout the analysis process. Only an initial estimate of the stiffness and damping properties of the physical components is required. It is demonstrated that the method is unconditionally stable provided that specific conditions are fulfilled and that the order accuracy can be maintained in the nonlinear regime without involving any matrix inversion while testing. The method also features controllable numerical energy dissipation characteristics and explicit expression of the target displacement and velocity vectors. The stability and accuracy of the proposed integration scheme are examined in the paper. The method has also been verified through hybrid testing performed of SDOF and MDOF structures with linear and highly nonlinear physical substructures. The results are compared with those obtained from the operator splitting method. An approach based on the modal decomposition principle is presented to predict the potential effect of experimental errors on the overall response during testing. Copyright © 2009 John Wiley & Sons, Ltd.     

壳体稳定性问题的研究进展

为探索大地震成因的物理机理,根据板壳理论的研究结果,综述多年来柔韧壳体非线性大挠度变形研究中与地震成因相关的壳体结构变形问题的研究进展、壳体结构变形过程中的稳定性问题。同时,介绍壳体曲率变化对其变形过程中稳定性的影响,期望能为探索大地震成因机理的研究提供某些有用的物理依据与理论基础。     

Stability of average acceleration method for structures with nonlinear damping

The energy approach is used to theoretically verify that the average acceleration method （AAM）, which is unconditionally stable for linear dynamic systems, is also unconditionally stable for structures with typical nonlinear damping, including the special case of velocity power type damping with a bilinear restoring force model. Based on the energy approach, the stability of the AAM is proven for SDOF structures using the mathematical features of the velocity power function and for MDOF structures by applying the virtual displacement theorem. Finally, numerical examples are given to demonstrate the accuracy of the theoretical analysis.