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.

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.

THE SINGULAR PERTURBATION METHOD APPLIED TO THE NONLINEAR STABILITY PROBLEM OF A SHALLOW SPHERICAL SHELL(Ⅱ)

In this paper we consider the nonlinear stability of a thin elastic circular shallow spherical shell under the action of uniform normal pressure with a clamped edge. When the geometrical parameter k is large, the uniformly valid asymptotic solutions are obtained by means of the singular perturbation method. In addition, we give the analytic formula for determining the centre deflection and the critical load, and the stability curve is also derived. This paper is a continuation of the author's previous paper[11]

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.

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.

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.

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.

An Inf-Sup Stabilized Finite Element Method by Multiscale Functions for the Stokes Equations

In the paper,an inf-sup stabilized finite element method by multiscale functions for the Stokes equations is discussed.The key idea is to use a PetrovGalerkin approach based on the enrichment of the standard polynomial space for the velocity component with multiscale functions.The inf-sup condition for P_(1)-P_(0)triangular element(or Q_(1)-P_(0)quadrilateral element)is established.The optimal error estimates of the stabilized finite element method for the Stokes equations are obtained.

Design and operation problems related to water curtain system forunderground water-sealed oil storage caverns

The underground water-sealed storage technique is critically important and generally accepted for the national energy strategy in China. Although several small underground water-sealed oil storage caverns have been built in China since the 1970s, there is still a lack of experience for large-volume underground storage in complicated geological conditions. The current design concept of water curtain system and the technical instruction for system operation have limitations in maintaining the stability of surrounding rock mass during the construction of the main storage caverns, as well as the long-term stability. Although several large-scale underground oil storage projects are under construction at present in China, the design concepts and construction methods, especially for the water curtain system, are mainly based on the ideal porosity medium flow theory and the experiences gained from the similar projects overseas. The storage projects currently constructed in China have the specific features such as huge scale, large depth, multiple-level arrangement, high seepage pressure, complicated geological conditions, and high in situ stresses, which are the challenging issues for the stability of the storage caverns. Based on years’ experiences obtained from the first large-scale （millions of cubic meters） underground water-sealed oil storage project in China, some design and operation problems related to water curtain system during project construction are discussed. The drawbacks and merits of the water curtain system are also presented. As an example, the conventional concept of “filling joints with water” is widely used in many cases, as a basic concept for the design of the water curtain system, but it is immature. In this paper, the advantages and disadvantages of the conventional concept are pointed out, with respect to the long-term stability as well as the safety of construction of storage caverns. Finally, new concepts and principles for design and construction of the underground water-sealed oil storage caverns are proposed.

Neutrality Criteria for Stability Analysis of Dynamical Systems through Lorentz and Rossler Model Problems

Two methods of stability analysis of systems described by dynamical equations are being considered. They are based on an analysis of eigenvalues spectrum for the evolutionary matrix or the spectral equation and they allow determining the conditions of stability and instability, as well as the possibility of chaotic behavior of systems in case of a stability loss. The methods are illustrated for nonlinear Lorenz and Rossler model problems.

Evaluation of embankment foundation creep in conditions of deep bedding of roof in permafrost soils

Railroad operating experience in permafrost conditions has shown that deformations of embankments on thawing foun- dations last for a long time. After an initial period of heat settlement due to permafrost degradation, the determining factor is the plastic flow of seasonal thawed soils of the foundation upper layer under the embankment. This paper provides a method to evaluate these deformations, and calculation examples using data from line sections of the Chum-Labytnangi Railway in northwestern Russia. It also discusses several methods of embankment stabilization, including the use of ver- tical thermosiphons.

Analysis of nonlinear stability and post-buckling for Euler-type beam-column structure

Based on the assumption of finite deformation, the Hamilton variational principle is extended to a nonlinear elastic Euler-type beam-column structure located on a nonlinear elastic foundation. The corresponding three-dimensional （3D） mathematical model for analyzing the nonlinear mechanical behaviors of structures is established, in which the effects of the rotation inertia and the nonlinearity of material and geometry are considered. As an application, the nonlinear stability and the post-buckling for a linear elastic beam with the equal cross-section located on an elastic foundation are analyzed. One end of the beam is fully fixed, and the other end is partially fixed and subjected to an axial force. A new numerical technique is proposed to calculate the trivial solution, bifurcation points, and bifurcation solutions by the shooting method and the Newton- Raphson iterative method. The first and second bifurcation points and the corresponding bifurcation solutions are calculated successfully. The effects of the foundation resistances and the inertia moments on the bifurcation points are considered.

Large eddy simulation for wind field analysis based on stabilized finite element method

In this paper,a stabilized finite element technique,actualized by streamline upwind Petrov-Galerkin(SUPG) stabilized method and three-step finite element method(FEM),for large eddy simulation(LES) is developed to predict the wind flow with high Reynolds numbers.Weak form of LES motion equation is combined with the SUPG stabilized term for the spatial finite element discretization.An explicit three-step scheme is implemented for the temporal discretization.For the numerical example of 2D wind flow over a square rib at Re=4.2×105,the Smagorinsky's subgrid-scale(SSGS) model,the DSGS model,and the DSGS model with Cabot near-wall model are applied,and their results are analyzed and compared with experimental results.Furthermore,numerical examples of 3D wind flow around a surface-mounted cube with different Reynolds numbers are performed using DSGS model with Cabot near-wall model based on the present stabilized method to study the wind field and compared with experimental and numerical results.Finally,vortex structures for wind flow around a surface-mounted cube are studied by present numerical method.Stable and satisfactory results are obtained,which are consistent with most of the measurements even under coarse mesh.

ERROR ANALYSIS OF A STABILIZED FINITE ELEMENT METHOD FOR THE GENERALIZED STOKES PROBLEM

This paper is devoted to the establishment of sharper a priori stability and error estimates of a stabilized finite element method proposed by Barrenechea and Valentin for solving the generalized Stokes problem,which involves a viscosity v and a reaction consta nt a.With the establishmen t of sharper st ability estimates and the help of ad hoc finite element projections,we can explicitly establish the dependence of error bounds of velocity and pressure on the viscosity z/,the reaction constant cr,and the mesh size h.Our analysis reveals that the viscosity y and the reaction constant a respectively act in the numerator position and the denominator position in the error estimates of velocity and pressure in standard norms without any weights.Consequently,the stabilization method is indeed suitable for the generalized Stokes problem with a small viscosity u and a large reaction constant a.The sharper error estimates agree very well with the numerical results.