A new stabilized method for quasi-Newtonian flows

For a generalized quasi-Newtonian flow, a new stabilized method focused on the low-order velocity-pressure pairs, （bi）linear/（bi）linear and （bi）linear/constant element, is presented. The pressure projection stabilized method is extended from Stokes problems to quasi-Newtonian flow problems. The theoretical framework developed here yields an estimate bound, which measures error in the approximate velocity in the W 1,r（Ω） norm and that of the pressure in the L r＇ （Ω） （1/r ＋ 1/r＇ = 1）. The power law model and the Carreau model are special ones of the quasi-Newtonian flow problem discussed in this paper. Moreover, a residual-based posterior bound is given. Numerical experiments are presented to confirm the theoretical results.

Two-level stabilized finite element method for Stokes eigenvalue problem

A two-level stabilized finite element method for the Stokes eigenvalue problem based on the local Gauss integration is considered. This method involves solving a Stokes eigenvalue problem on a coarse mesh with mesh size H and a Stokes problem on a fine mesh with mesh size h -- O（H2）, which can still maintain the asymptotically optimal accuracy. It provides an approximate solution with the convergence rate of the same order as the usual stabilized finite element solution, which involves solving a Stokes eigenvalue problem on a fine mesh with mesh size h. Hence, the two-level stabilized finite element method can save a large amount of computational time. Moreover, numerical tests confirm the theoretical results of the present method.

耦合拉格朗日-欧拉方法及其在海洋工程中的应用

Combining the strengths of Lagrangian and Eulerian descriptions,the coupled Lagrangian–Eulerian methods play an increasingly important role in various subjects.This work reviews their development and application in ocean engineering.Initially,we briefly outline the advantages and disadvantages of the Lagrangian and Eulerian descriptions and the main characteristics of the coupled Lagrangian–Eulerian approach.Then,following the developmental trajectory of these methods,the fundamental formulations and the frameworks of various approaches,including the arbitrary Lagrangian–Eulerian finite element method,the particle-in-cell method,the material point method,and the recently developed Lagrangian–Eulerian stabilized collocation method,are detailedly reviewed.In addition,the article reviews the research progress of these methods with applications in ocean hydrodynamics,focusing on free surface flows,numerical wave generation,wave overturning and breaking,interactions between waves and coastal structures,fluid–rigid body interactions,fluid–elastic body interactions,multiphase flow problems and visualization of ocean flows,etc.Furthermore,the latest research advancements in the numerical stability,accuracy,efficiency,and consistency of the coupled Lagrangian–Eulerian particle methods are reviewed;these advancements enable efficient and highly accurate simulation of complicated multiphysics problems in ocean and coastal engineering.By building on these works,the current challenges and future directions of the hybrid Lagrangian–Eulerian particle methods are summarized.

Direct discontinuous Galerkin method for the generalized Burgers-Fisher equation

In this study, we use the direct discontinuous Galerkin method to solve the generalized Burgers-Fisher equation. The method is based on the direct weak formulation of the Burgers-Fisher equation. The two adjacent cells are jointed by a numerical flux that includes the convection numerical flux and the diffusion numerical flux. We solve the ordinary differential equations arising in the direct Galerkin method by using the strong stability preserving Runge^Kutta method. Numerical results are compared with the exact solution and the other results to show the accuracy and reliability of the method.

3D stability analysis method of concave slope based on the Bishop method

In order to study the stability control mechanism of a concave slope with circular landslide, and remove the influence of differences in shape on slope stability, the limit analysis method of a simplified Bishop method was employed. The sliding body was divided into strips in a three-dimensional model, and the lateral earth pressure was put into mechanical analysis and the three-dimensional stability analysis methods applicable for circular sliding in concave slope were deduced. Based on geometric structure and the geological parameters of a concave slope, the influence rule of curvature radius and the top and bottom arch height on the concave slope stability were analyzed. The results show that the stability coefficient decreases after growth, first in the transition stage of slope shape from flat to concave, and it has been confirmed that there is a best size to make the slope stability factor reach a maximum. By contrast with average slope, the stability of a concave slope features a smaller range of ascension with slope height increase, which indicates that the enhancing effect of a concave slope is apparent only with lower slope heights.

Limit analysis method for active earth pressure on laggings between stabilizing piles

Stabilizing pile is a kind of earth shoring structure frequently used in slope engineering. When the piles have cantilever segments above the ground,laggings are usually installed to avoid collapse of soil between piles. Evaluating the earth pressure acting on laggings is of great importance in design process.Since laggings are usually less stiff than piles,the lateral pressure on lagging is much closer to active earth pressure. In order to estimate the lateral earth pressure on lagging more accurately,first,a model test of cantilever stabilizing pile and lagging systems was carried out. Then,basing the experimental results a three-dimensional sliding wedge model was established. Last,the calculation process of the total active force on lagging is presented based on the kinematic approach of limit analysis. A comparison is made between the total active force on lagging calculated by the formula presented in this study and the force on a same-size rigid retaining wall obtained from Rankine's theory. It is found that the proposed method fits well with the experimental results.Parametric studies show that the total active force on lagging increases with the growth of the lagging height and the lagging clear span; while decreases asthe soil internal friction angle and soil cohesion increase.

Stabilization meshless method for convection-dominated problems

It is weN-known that the standard Galerkin is not ideally suited to deal with the spatial discretization of convection-dominated problems. In this paper, several techniques are proposed to overcome the instabilitY issues in convection-dominated problems in the simulation with a meshless method. These stable techniques included nodal refinement, enlargement of the nodal influence domain, full upwind meshless technique and adaptive upwind meshless technique. Numerical results for sample problems show that these techniques are effective in solving convection-dominated problems, and the adaptive upwind meshless technique is the most effective method of all.

GENERALIZED TRANSFER FUNCTION OF CONTROL SYSTEM AND AN IMPROVEMENT ON THE DECISION METHOD OF MOVEMENT STABILITY

In this paper the definitions of generalized transfer functios of control system and itscontinuity are presented.Using generalized transfer function as a tool,a set of theorems fordeciding movement stability have been constructed.Thus basing understanding of thecharacteristics of a control dynamics system on its measured procedure will simplify thedecision method of movement stability problems.

Stability analysis of liquid filled spacecraft system with flexible attachment by using the energy–Casimir method

The stability of partly liquid filled spacecraft with flexible attachment was investigated in this paper. Liquid sloshing dynamics was simplified as the spring-mass model, and flexible attachment was modeled as the linear shearing beam. The dynamic equations and Hamiltonian of the coupled spacecraft system were given by analyzing the rigid body, liquid fuel, and flexible appendage. Nonlinear stability conditions of the coupled spacecraft system were derived by computing the variation of Casimir function which was added to the Hamiltonian. The stable region of the parameter space was given and validated by numerical computation. Related results suggest that the change of inertia matrix, the length of flexible attachment, spacecraft spinning rate, and filled ratio of liquid fuel tank have strong influence on the stability of the spacecraft system.

Reinforcement strength reduction in FEM for mechanically stabilized earth structures

The factor of safety of mechanically stabilized earth(MSE) structures can be analyzed either using limit equilibrium method(LEM) or strength reduction method(SRM) in finite element/difference method. In LEM, the strengths of the reinforcement members and soils are reduced with the same factor. While using the SRM, only soil strength is reduced during the calculation of the factor of safety. This causes inconsistence in calculating the factor of safety of the MSE structures. To overcome this, an iteration method is proposed to consider the strength reduction of the reinforcements in SRM. The method is demonstrated by using PLAXIS, a finite element software. The results show that the factor of safety converges after a few iterations. The reduction of strength has different effects on the factor of safety depending on the properties of the reinforcements and the soil, and failure modes.

THE STABILITY AND CONVERGENCE OF THE FINITE ANALYTIC METHOD FOR THE NUMERICAL SOLUTION OF CONVECTIVE DIFFUSION EQUATION

In this paper we make a close study of the finite analytic method by means of the maximum principles in differential equations and give the proof of the stability and convergence of the finite analytic method.

Error Analysis of the Nonconforming P1 Finite Element Method to the Sequential Regularization Formulation for Unsteady Navier-Stokes Equations

In this paper we investigate the nonconforming P_(1) finite element ap-proximation to the sequential regularization method for unsteady Navier-Stokes equations.We provide error estimates for a full discretization scheme.Typi-cally,conforming P_(1) finite element methods lead to error bounds that depend inversely on the penalty parameter ∈.We obtain an ∈-uniform error bound by utilizing the nonconforming P_(1) finite element method in this paper.Numerical examples are given to verify theoretical results.

Uniform analysis of a stabilized hybrid finite element method for Reissner-Mindlin plates

This paper presents a low order stabilized hybrid quadrilateral finite element method for ReissnerMindlin plates based on Hellinger-Reissner variational principle,which includes variables of displacements,shear stresses and bending moments.The approach uses continuous piecewise isoparametric bilinear interpolations for the approximations of the transverse displacement and rotation.The stabilization achieved by adding a stabilization term of least-squares to the original hybrid scheme,allows independent approximations of the stresses and moments.The stress approximation adopts a piecewise independent 4-parameter mode satisfying an accuracy-enhanced condition.The approximation of moments employs a piecewise-independent 5-parameter mode.This method can be viewed as a stabilized version of the hybrid finite element scheme proposed in [Carstensen C,Xie X,Yu G,et al.A priori and a posteriori analysis for a locking-free low order quadrilateral hybrid finite element for Reissner-Mindlin plates.Comput Methods Appl Mech Engrg,2011,200:1161-1175],where the approximations of stresses and moments are required to satisfy an equilibrium criterion.A priori error analysis shows that the method is uniform with respect to the plate thickness t.Numerical experiments confirm the theoretical results.

Stabilized Finite Element Methods for Biot’s Consolidation Problems Using Equal Order Elements

Using the standard mixed Galerkin methods with equal order elements to solve Biot’s consolidation problems,the pressure close to the initial time produces large non-physical oscillations.In this paper,we propose a class of fully discrete stabilized methods using equal order elements to reduce the effects of non-physical oscillations.Optimal error estimates for the approximation of displacements and pressure at every time level are obtained,which are valid even close to the initial time.Numerical experiments illustrate and confirm our theoretical analysis.

Three-dimensional stability of landslides based on local safety factor

Unlike the limit equilibrium method(LEM), with which only the global safety factor of the landslide can be calculated, a local safety factor(LSF) method is proposed to evaluate the stability of different sections of a landslide in this paper. Based on three-dimensional(3D) numerical simulation results, the local safety factor is defined as the ratio of the shear strength of the soil at an element on the slip zone to the shear stress parallel to the sliding direction at that element. The global safety factor of the landslide is defined as the weighted average of all local safety factors based on the area of the slip surface. Some example analyses show that the results computed by the LSF method agree well with those calculated by the General Limit Equilibrium(GLE) method in two-dimensional(2D) models and the distribution of the LSF in the 3D slip zone is consistent with that indicated by the observed deformation pattern of an actual landslide in China.

Sensitivity analysis of influencing parameters in cavern stability

In order to analyze the stability of the underground rock structures,knowing the sensitivity of geomechanical parameters is important.To investigate the priority of these geomechanical properties in the stability of cavern,a sensitivity analysis has been performed on a single cavern in various rock mass qualities according to RMR using Phase 2.The stability of cavern has been studied by investigating the side wall deformation.Results showed that most sensitive properties are coefficient of lateral stress and modulus of deformation.Also parameters of Hoek-Brown criterion and r c have no sensitivity when cavern is in a perfect elastic state.But in an elasto-plastic state,parameters of Hoek-Brown criterion and r c affect the deformability;such effect becomes more remarkable with increasing plastic area.Other parameters have different sensitivities concerning rock mass quality(RMR).Results have been used to propose the best set of parameters for study on prediction of sidewall displacement.

Intra-layer synchronization in duplex networks

This paper explores the intra-layer synchronization in duplex networks with different topologies within layers and different inner coupling patterns between, within, and across layers. Based on the Lyapunov stability method, we prove theoretically that the duplex network can achieve intra-layer synchronization under some appropriate conditions, and give the thresholds of coupling strength within layers for different types of inner coupling matrices across layers. Interestingly,for a certain class of coupling matrices across layers, it needs larger coupling strength within layers to ensure the intra-layer synchronization when the coupling strength across layers become larger, intuitively opposing the fact that the intra-layer synchronization is seemly independent of the coupling strength across layers. Finally, numerical simulations further verify the theoretical results.

Slope stability analysis of Southern slope of Chengmenshan copper mine,China

The engineering geology and hydrogeology in the southern slope of Chengmenshan copper mine are very complicated,because there is a soft-weak layer between two kinds of sandstones.Field investigations demonstrate that some instability problems might occur in the slope.In this research,the southern slope,which is divided into six sections(I-0,I-1,I-2,II-0,II-1 and II-2),is selected for slope stability analysis using limit equilibrium and numerical method.Stability results show that the values of factor of safety(FOS) of sections I-0,I-1 and I-2 are very low and slope failure is likely to happen.Therefore reinforcement subjected to seismic,water and weak layer according to sections were carried out to increase the factor of safety of the three sections,two methods were used;grouting with hydration of cement and water to increase the cohesion(c) and pre-stressed anchor.Results of reinforcement showed that factor of safety increased more than 1.15.

Seismic stability evaluation of embankment slope based on catastrophe theory

An evaluation method for the seismic stability of embankment slope was presented based on catastrophe theory. Seven control factors, including internal frictional angle, cohesion force, slope height, slope angle, surface gradients, peak acceleration, and distance to fault were selected for analysis of multi-level objective decomposition. According to the normalization formula and the fuzzy subject function produced by combination of catastrophe theory and fuzzy math, a recursive calculation was carried out to obtain a catastrophic affiliated functional value, which can be used to evaluate the seismic stability of embankment slope. Fifteen samples were used to verify the effectiveness of this method. The results show that compared with the traditional quantitative method, the catastrophe progression owns higher accuracy and good application potential in predicting the seismic stability of embankment slope.

LOCALLY STABILIZED FINITE ELEMENT METHOD FOR STOKES PROBLEM WITH NONLINEAR SLIP BOUNDARY CONDITIONS

Based on the low-order conforming finite element subspace （Vh, Mh） such as the P1-P0 triangle element or the Q1-P0 quadrilateral element, the locally stabilized finite element method for the Stokes problem with nonlinear slip boundary conditions is investigated in this paper. For this class of nonlinear slip boundary conditions including the subdifferential property, the weak variational formulation associated with the Stokes problem is an variational inequality. Since （Vh, Mh） does not satisfy the discrete inf-sup conditions, a macroelement condition is introduced for constructing the locally stabilized formulation such that the stability of （Vh, Mh） is established. Under these conditions, we obtain the H1 and L2 error estimates for the numerical solutions.