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

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.

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.

Application of Arnoldi method to boundary layer instability

The Arnoldi method is applied to boundary layer instability, and a finite difference method is employed to avoid the limit of the finite element method. This modus operandi is verified by three comparison cases, i.e., comparison with linear stability theory（LST） for two-dimensional（2D） disturbance on one-dimensional（1D） basic flow, comparison with LST for three-dimensional（3D） disturbance on 1D basic flow, and comparison with Floquet theory for 3D disturbance on 2D basic flow. Then it is applied to secondary instability analysis on the streaky boundary layer under spanwise-localized free-stream turbulence（FST）. Three unstable modes are found, i.e., an inner mode at a high-speed center streak, a sinuous type outer mode at a low-speed center streak, and a sinuous type outer mode at low-speed side streaks. All these modes are much more unstable than Tollmien–Schlichting（TS） waves, implying the dominant contribution of secondary instability in bypass transition. The modes at strong center streak are more unstable than those at weak side streaks, so the center streak is ‘dangerous＇ in secondary instability.

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.

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.

Multiple solutions and stability of the steady transonic small-disturbance equation

Numerical solutions of the steady transonic small-disturbance（TSD） potential equation are computed using the conservative Murman-Cole scheme. Multiple solutions are discovered and mapped out for the Mach number range at zero angle of attack and the angle of attack range at Mach number 0.85 for the NACA 0012 airfoil. We present a linear stability analysis method by directly assembling and evaluating the Jacobian matrix of the nonlinear finite-difference equation of the TSD equation. The stability of all the discovered multiple solutions are then determined by the proposed eigen analysis. The relation of stability to convergence of the iterative method for solving the TSD equation is discussed. Computations and the stability analysis demonstrate the possibility of eliminating the multiple solutions and stabilizing the remaining unique solution by adding a sufficiently long splitter plate downstream the airfoil trailing edge. Finally, instability of the solution of the TSD equation is shown to be closely connected to the onset of transonic buffet by comparing with experimental data.

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.

Tip-splitting instability in directional solidification based on bias field method

Tip splitting instability of cellular interface morphology in directional solidification is analyzed based on the bias field method proposed recently by Glicksman. The physical mechanism of tip instability is explained by analyzing the interface potential, the tangential energy flux, and the normal energy flux. A rigorous criterion for tip-splitting instability is established analytically, i.e., the ratio of the cellular tip radius to the cellular width α 〉3/2/π≈ 0.3899, which is in good agreement with simulation results. This study also reveals that the cellular tip splitting instability is attributable to weak Gibbs–Thomson energy acting on the interface.

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.

Stability Criteria for Large－Scale Linear Systems with Structured Uncertainties

The robust stability analysis for large scale linear systems with structured time varying uncertainties is investigated in this paper.By using the scalar Lyapunov functions and the properties of M matrix and nonnegative matrix,stability robustness measures are proposed.The robust stability criteria obtained are applied to derive an algebric criterion which is expressed directly in terms of plant parameters and is shown to be less conservative than the existing ones.A numerical example is given to demonstrate the stability criteria obtained and to compare them with the previous ones.

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.

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.

Unstable and exact periodic solutions of three-particles time-dependent FPU chains

For lower dimensional Fermi–Pasta–Ulam（FPU） chains, the α-chain is completely integrable and the Hamiltonian of the β-chain can be identified with the H′enon–Heiles Hamiltonian. When the strengths α, β of the nonlinearities depend on time periodically with the same frequencies as the natural angular frequencies, the resonance phenomenon is inevitable. In this paper, for certain periodic functions α（t） and β（t） with resonance frequencies, we give the existence and stability of some nontrivial exact periodic solutions for a one-dimensional αβ-FPU model composed of three particles with periodic boundary conditions.

Solution to multiple attribute group decision making problems with two decision makers

A kind of multiple attribute group decision making （MAGDM） problem is discussed from the perspective of statistic decision-making. Firstly, on the basis of the stability theory, a new idea is proposed to solve this kind of problem. Secondly, a con- crete method corresponding to this kind of problem is proposed. The main tool of our research is the technique o~ the jackknife method. The main advantage of the new method is that it can identify and determine the reliability degree of the existed decision making information. Finally, a traffic engineering example is given to show the effectiveness of the new method.

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.

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.

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.

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.