Numerical manifold method for dynamic nonlinear analysis of saturated porous media

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.     

Restricted Three-Body Problem: An Approximation of its General Solution – Part One – the Manifold of Symmetric Periodic Solutions

This paper gives the results of a programme attempting to exploit ‘la seule bréche’ (Poincaré, 1892, p. 82) of non-integrable systems, namely to develop an approximate general solution for the three out of its four component-solutions of the planar restricted three-body problem. This is accomplished by computing a large number of families of ‘solutions précieuses’ (periodic solutions) covering densely the space of initial conditions of this problem. More specifically, we calculated numerically and only for μ = 0.4, all families of symmetric periodic solutions (1st component of the general solution) existing in the domain D:(x ₀ ∊ [−2,2],C ∊ [−2,5]) of the (x ₀, C) space and consisting of symmetric solutions re-entering after 1 up to 50 revolutions (see graph in Fig. 4). Then we tested the parts of the domain D that is void of such families and established that they belong to the category of escape motions (2nd component of the general solution). The approximation of the 3rd component (asymmetric solutions) we shall present in a future publication. The 4th component of the general solution of the problem, namely the one consisting of the bounded non-periodic solutions, is considered as approximated by those of the 1st or the 2nd component on account of the `Last Geometric Theorem of Poincaré' (Birkhoff, 1913). The results obtained provoked interest to repeat the same work inside the larger closed domain D:(x ₀ ∊ [−6,2], C ∊ [−5,5]) and the results are presented in Fig. 15. A test run of the programme developed led to reproduction of the results presented by Hénon (1965) with better accuracy and many additional families not included in the sited paper. Pointer directions construed from the main body of results led to the definition of useful concepts of the basic family of order n, n = 1, 2,… and the completeness criterion of the solution inside a compact sub-domain of the (x ₀, C) space. The same results inspired the ‘partition theorem’, which conjectures the possibility of partitioning an initial conditions domain D into a finite set of sub-domains D i that fulfill the completeness criterion and allow complete approximation of the general solution of this problem by computing a relatively small number of family curves. The numerical results of this project include a large number of families that were computed in detail covering their natural termination, the morphology, and stability of their member solutions. Zooming into sub-domains of D permitted clear presentation of the families of symmetric solutions contained in them. Such zooming was made for various values of the parameter N, which defines the re-entrance revolutions number, which was selected to be from 50 to 500. The areas generating escape solutions have being investigated. In Appendix A we present families of symmetric solutions terminating at asymptotic solutions, and in Appendix B the morphology of large period symmetric solutions though examples of orbits that re-enter after from 8 to 500 revolutions. The paper concludes that approximations of the general solution of the planar restricted problem is possible and presents such approximations, only for some sub-domains that fulfill the completeness criterion, on the basis of sufficiently large number of families.     

Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

伶仃洋沉积动力特点的研究

位于珠江三角洲东侧的伶仃洋，因径流下泄与潮流进退的流向不一，使它各分流口门的出口水道都有主槽和支槽之分，即都有主干水道和分汊水道。陆架高盐海水入侵又使伶仃洋内沉积动力过程在空间分布上发生差异，如沉积物分布有粗－细－稍粗之分；而水体中的密度、速度差异，常常产生锋带，对水下地形的发展有不可忽视的影响。因而全面认识发生在伶仃洋内的沉积动力作用，对深水航道的选线极为重要。     

多传感器观测下带乘性噪声系统的逆向滤波与反褶积融合算法

针对多传感器观测环境下带乘性噪声系统的逆向最优滤波与反褶积融合估计问题 ,本文提出了 1种基于极大似然准则的最优融合算法。该算法中各单传感器间并行计算 ,并且融合中心与单传感器处理中心间无反向通讯 ,因而执行效率较高。仿真表明 ,该融合算法产生的逆向滤波与反褶积比单传感器处理结果有较明显提高     

钦州湾潮流特征分析

根据1994年5月和11—12月两个航次的调查及近年有关海湾的海流观测资料，分析了钦州湾的潮流特征，观测资料来自分布于该海湾东中西三个航道及相关区域的8个测站。钦州湾涨落潮流特征：落潮流速大于涨潮流速：东部最大涨落潮流速小于西部：夏季落潮流总是大于冬季：龙门水道附近的流速最大，而其余区域的流速相对较小。外湾(钦州湾)余流，是气旋式环流：水体东进而西出。从湾内来的泥沙和污染物质，主要从西部进入外海；中、东槽是主要“进水”通道，外海低泥沙含量、少污染物质的清洁水从这里潮流而上。因此，与西部相比，东部盐度高，水质清洁，底质重金属含量少。     

Semi-Probabilistic Evaluation of System Redundancy for Existing Offshore Structures

In this paper, a new definition of structure system redundancy is proposed in view of the various measures for structure redundancy. By introducing the terms of structure system failure at the mechanism level and equivalent reliability index, the safety for existing offshore platforms can be evaluated by the semi-probabilistic method presented in this paper. Some numerical examples are given and satisfactory results have been obtained.     

连云港近岸海域物理自净能力研究及水质预测——Ⅲ.西大堤建成后对物质输运的影响

在连云港近岸海域计算潮流场基础上建立拉格朗日余流模型,并对连云港市两大堤建成前后的拉格朗日余流变化进行了分析,且选择有代表性的排污口进行了数值跟踪。     

Numerical study of enhanced energy dissipation near a seamount

Using a two-dimensional primitive equation model, we examine nonlinear responses of a semidiurnal tidal flow impinging on a seamount with a background Garrett-Munk-like (GM-like) internal wavefield. It is found that horizontally elongated pancake-like structures of high vertical wavenumber near-inertial current shear are created both in the near-field (the region over the slope of the seamount) and far-field (the region over the flat bottom of the ocean). An important distinction is that the high vertical wavenumber near-inertial current shear is amplified only at mid-latitudes in the far-field (owing to a parametric subharmonic instability (PSI)), whereas it is amplified both at mid-and high-latitudes (above the latitude where PSI can occur) in the near-field. In order to clarify the generating mechanism for the strong shear in the near-field, additional numerical experiments are carried out with the GM-like background internal waves removed. The experiments show that the strong shear is also created, indicating that it is not caused by the interaction between the background GM-like internal waves and the semidiurnal internal tides. One possible explanation is proposed for the amplification of high vertical wavenumber near-inertial current shear in the near-field where tide residual flow resulting from tide-topography interaction plays an important role in transferring energy from high-mode internal tides to near-inertial internal waves.