深海采矿柔性立管系统非线性有限元模拟

在海洋环境载荷及集矿机牵引作用下,深海采矿柔性立管系统的动力学响应涉及几何非线性和非保守载荷的双重非线性源.基于三维固体有限变形理论建立数学模型,在完全拉格朗日格式下推导了系统运动平衡方程,针对非保守载荷的等效计算、非线性切线刚度矩阵及非线性方程的数值求解等关键问题提出了有效处理方案.根据处理方案开发了数值计算程序,并...     

Direct numerical simulation of particle–fluid systems by combining time-driven hard-sphere model and lattice Boltzmann method

A coupled numerical method for the direct numerical simulation of particle–fluid systems is formulated and implemented, resolving an order of magnitude smaller than particle size. The particle motion is described by the time-driven hard-sphere model, while the hydrodynamic equations governing fluid flow are solved by the lattice Boltzmann method (LBM). Particle–fluid coupling is realized by an immersed boundary method (IBM), which considers the effect of boundary on surrounding fluid as a restoring force added to the governing equations of the fluid. The proposed scheme is validated in the classical flow-around-cylinder simulations, and preliminary application of this scheme to fluidization is reported, demonstrating it to be a promising computational strategy for better understanding complex behavior in particle–fluid systems.     

Parametric vibrations and stability of an axially accelerating string guided by a non-linear elastic foundation

Parametric vibrations and stability of an axially accelerating string guided by a non-linear elastic foundation are studied analytically. The axial speed, as the source of parametric vibrations, is assumed to involve a mean speed, along with small harmonic variations. The method of multiple scales is applied to the governing non-linear equation of motion and then the natural frequencies and mode shape equations of the system are derived using the equation of order one, and satisfying the compatibility conditions. Using the equation of order epsilon, the solvability conditions are obtained for three distinct cases of axial acceleration frequency. For all cases, the stability areas of system are constructed analytically. Finally, some numerical simulations are presented to highlight the effects of system parameters on vibration, natural frequencies, frequency-response curves, stability, and bifurcation points of the system.     

剪切流场中含内流立管横向涡激振动特性

立管是海洋工程中输送油气或其他矿产资源的必备结构, 外部洋流引起的立管涡激振动影响着立管的疲劳寿命, 危害深海资源开发. 本文基于欧拉?伯努利梁方程, 结合半经验时域水动力模型, 建立剪切流与内流耦合作用下海洋立管涡激振动预报模型, 运用有限元方法和Newmark-β逐步积分法求解方程, 首先将数值模拟结果与实验数据进行对比, 验证模型正确性. 然后, 运用此模型, 对剪切流作用下含内流的顶张立管在不同内流速度和密度下的横向涡激振动响应特性进行研究, 主要分析了立管的横向振动模态、振动频率以及均方根位移等涡激振动参数随内流速度和密度等参数的变化规律. 结果表明, 在剪切流场中, 含内流海洋立管在横向上表现出多模态多频率的涡激振动;立管横向振动的最大均方根位移随内流速度和密度的增大而增大, 特别是当内流速度较大时, 横向最大均方根位移增大明显;立管横向振动的主导频率随内流速度和密度的增大而减小, 并且内流密度的增大同样会引起模态转换和频率转换.      

A harmonic balance technique for the reduced-order computation of vortex-induced vibration

We present a harmonic balance (HB) method to model frequency lock-in effect during vortex-induced vibration (VIV) of elastically mounted circular cylinder and a flexible riser section in a freestream uniform flow. The fluid flow and structure are coupled by a fixed-point iteration process through a frequency updating algorithm. By minimizing the structural residual in the standard least-square norm, the convergence of HB-based fixed-point algorithm is achieved for a range of reduced velocity. To begin with, the HB solver is first assessed for a periodic unsteady flow around a stationary circular cylinder. A freely vibrating circular cylinder is then adopted for the reduced-order computation of VIV at low Reynolds numbers of Re=100 and 180 with one- and two-degrees-of-freedom. The coupled VIV dynamics and the frequency lock-in phenomenon are accurately captured. The results show that the HB solver is able to predict the amplitude of vibration, frequency and forces comparable to its time domain counterpart, while providing a significant reduction with regard to overall computational cost. The proposed new scheme is then demonstrated for a fully-coupled three dimensional (3D) analysis of a linear-elastic riser section undergoing vortex-induced vibration in the lock-in range. The results reveal the 3D effects through isosurfaces of streamwise vorticity blobs distributed over the span of flexible riser section. In comparison to time domain results, the 3D flow-structure interactions are accurately predicted while providing a similar speed up rate that of 2D simulations. This further corroborates that the HB solver can be extended to 3D flow-structure dynamics without compromising efficiency and accuracy.     

A robust high-resolution method for the time–frequency analysis of vortex-induced-vibration signals

The dynamic interaction between ocean current and marine riser is complex in nature, and the riser׳s vortex-induced vibration (VIV) due to the current often strongly exhibits a non-stationary phenomenon. For investigating the time-varying dominant frequencies of the VIV motion, a joint time–frequency analysis is necessary. Traditional methods for time–frequency analysis include the Short Time Fourier Transform (STFT) and Wavelet Transform (WT) methods, though both methods have significant drawbacks. Specifically, the STFT method suffers frequency resolution and leakage problems, while the WT method is highly sensitive to its basic wavelet selection. This paper newly introduces a robust high-resolution method, named the STPT-SS method, which is the Short Time Prony Transform (STPT) using a State-Space (SS) model. In particular, the STPT algorithm contributes to the high-resolution feature of the proposed method, and the SS model to the robustness. Using test VIV data that include a synthesized signal and measurements from laboratory and field experiments, the STPT-SS method is found to significantly outperform the STFT and WT methods in the time–frequency analysis.     

Numerical simulations of falling leaves using a pseudo-spectral method with volume penalization

The dynamics of falling leaves is studied by means of numerical simulations. The two-dimensional incompressible Navier–Stokes equations, coupled with the equations governing solid body dynamics, are solved using a Fourier pseudo-spectral method with volume penalization to impose no-slip boundary conditions. Comparison with other numerical methods is made. Simulations performed for different values of the Reynolds number show that its decrease stabilizes the free fall motion.     

A Newton's method scheme for solving free-surface flow problems

A Newton's method scheme is described for solving the system of non-linear algebraic equations arising when finite difference approximations are applied to the Navier–Stokes equations and their associated boundary conditions. The problem studied here is the steady, buoyancy-driven motion of a deformable bubble, assumed to consist of an inviscid, incompressible gas. The linear Newton system is solved using both direct and iterative equation solvers. The numerical results are in excellent agreement with previous work, and the method achieves quadratic convergence.     

A numerical method to simulate turbulent cavitating flows

The objective of this paper is to develop a numerical method for simulating multiphase cavitating flows on unstructured grids. The multiphase medium is represented using a homogeneous mixture model that assumes thermal equilibrium between the liquid and vapor phases. We develop a predictor–corrector approach to solve the governing Navier–Stokes equations for the liquid/vapor mixture, together with the transport equation for the vapor mass fraction. While a non-dissipative and symmetric scheme is used in the predictor step, a novel characteristic-based filtering scheme with a second order TVD filter is developed for the corrector step to handle shocks and material discontinuities in non-ideal gases and mixtures. Additionally, a sensor based on vapor volume fraction is proposed to localize dissipation to the vicinity of discontinuities. The scheme is first validated for simple one dimensional canonical problems to verify its accuracy in predicting jump conditions across material discontinuities and shocks. It is then applied to two turbulent cavitating flow problems – over a hydrofoil using RANS and over a wedge using LES. Our results show that the simulations are in good agreement with experimental data for the above tested cases, and that the scheme can be successfully applied to both RANS and LES methodologies.     

A numerical continuous model for the hydrodynamics of fluid particle systems

In order to understand the hydrodynamic interactions that can appear in a fluid particle motion, an original method based on the equations governing the motion of two immiscible fluids has been developed. These momentum equations are solved for both the fluid and solid phases. The solid phase is assumed to be a fluid phase with physical properties, such as its behaviour can be assimilated to that of pseudo‐rigid particles. The only unknowns are the velocity and the pressure defined in both phases. The unsteady two‐dimensional momentum equations are solved by using a staggered finite volume formulation and a projection method. The transport of each particle is solved by using a second‐order explicit scheme. The physical model and the numerical method are presented, and the method is validated through experimental measurements and numerical results concerning the flow around a circular cylinder. Good agreement is observed in most cases. The method is then applied to study the trajectory of one settling particle initially off‐centred between two parallel walls and the corresponding wake effects. Different particle trajectories related to particulate Reynolds numbers are presented and commented. A two‐body interaction problem is investigated too. This method allows the simulation of the transport of particles in a dilute suspension in reasonable time. One of the important features of this method is the computational cost that scales linearly with the number of particles. Copyright © 1999 John Wiley & Sons, Ltd.     

Free vibration and critical speed of moderately thick rotating laminated composite conical shell using generalized differential quadrature method

The generalized differential quadrature method (GDQM) is employed to consider the free vibration and critical speed of moderately thick rotating laminated composite conical shells with different boundary conditions developed from the first-order shear deformation theory (FSDT). The equations of motion are obtained applying Hamilton’s concept, which contain the influence of the centrifugal force, the Coriolis acceleration, and the preliminary hoop stress. In addition, the axial load is applied to the conical shell as a ratio of the global critical buckling load. The governing partial differential equations are given in the expressions of five components of displacement related to the points lying on the reference surface of the shell. Afterward, the governing differential equations are converted into a group of algebraic equations by using the GDQM. The outcomes are achieved considering the effects of stacking sequences, thickness of the shell, rotating velocities, half-vertex cone angle, and boundary conditions. Furthermore, the outcomes indicate that the rate of the convergence of frequencies is swift, and the numerical technique is superior stable. Three comparisons between the selected outcomes and those of other research are accomplished, and excellent agreement is achieved.     

海洋柔性结构涡激振动的流固耦合机理和响应

对近几十年来国内外在涡激振动的基础研究包括机理认识和动响应分析等方面的进展进行了论述,尤其针对海洋油气平台中的立管、隔水管等细长柔性结构的涡激振动.描述了涡激振动这种典型的非线性流固耦合现象所具有的特征,包括自激、自限制、展向相关、尾迹水动力与结构动力的流固耦合等及其主要影响参数.介绍了目前常用的结构响应预测方法和相关实验.通过讨论当前理论研究和实际工程中的热点问题,诸如多模态宽带振动、浮体运动与水下立管的耦合、响应抑制措施、双向振动、高雷诺数下的大尺度物理实验等,对今后该领域的研究方向进行了力所能及的展望.     

Nonlinear flutter behavior of a plate with motion constraints in subsonic flow

This paper is aimed at presenting the nonlinear flutter peculiarities of a cantilevered plate with motion-limiting constraints in subsonic flow. A non-smooth free-play structural nonlinearity is considered to model the motion constraints. The governing nonlinear partial differential equation is discretized in space and time domains by using the Galerkin method. The equilibrium points and their stabilities are presented based on qualitative analysis and numerical studies. The system loses its stability by flutter and undergoes the limit cycle oscillations (LCOs) due to the nonlinearity. A heuristic analysis scheme based on the equivalent linearization method is applied to theoretical analysis of the LCOs. The Hopf and two-multiple semi-stable limit cycle bifurcation bifurcations are supercritical or subcritical, which is dependent on the location of the motion constraints. For some special cases the bifurcations are, interestingly, both supercritical and subcritical. The influence of varying parameters on the dynamics is discussed in detail. The results predicted by the analysis scheme are in good agreement with the numerical ones.     

Fluid-solid coupling and dynamic response of vortex-induced vibration of slender ocean cylinders

对近几十年来国内外在涡激振动的基础研究包括机理认识和动响应分析等方面的进展进行了论述,尤其针对海洋油气平台中的立管、隔水管等细长柔性结构的涡激振动.描述了涡激振动这种典型的非线性流固耦合现象所具有的特征,包括自激、自限制、展向相关、尾迹水动力与结构动力的流固耦合等及其主要影响参数.介绍了目前常用的结构响应预测方法和相关实验.通过讨论当前理论研究和实际工程中的热点问题,诸如多模态宽带振动、浮体运动与水下立管的耦合、响应抑制措施、双向振动、高雷诺数下的大尺度物理实验等,对今后该领域的研究方向进行了力所能及的展望.     

Vortex-induced vibrations and waves under shear flow with a wake oscillator model

This paper addresses the interaction of a slender structure and a sheared incident cross-flow. The oscillating wake of the structure is modeled using a distribution of van der Pol oscillators. Elementary configurations are first considered to assess the basic phenomena and structure dynamics allowing to reliably investigate more complex and realistic cases. The scope of the study thus ranges from a forced oscillating cylinder to a tensioned cable. This last configuration is found to experience wave-packets of vortex-induced motion due to a series of local lock-ins. A theoretical analysis is carried-out to predict the wave-packets amplitude and distribution. It is shown to be in reasonable agreement with the results of numerical simulations. This indicates that the system behavior can be described in terms of the local interactions between the wake and the structure only.     

Further application of hybrid solution to another form of Boussinesq equations and comparisons

Recently, a new hybrid scheme is introduced for the solution of the Boussinesq equations. In this study, the hybrid scheme is used to solve another form of the Boussinesq equations. The hybrid solution is composed of finite‐volume and finite difference method. The finite‐volume method is applied to conservative part of the governing equations, whereas the higher order Boussinesq terms are discretized using the finite‐difference scheme. Fourth‐order accuracy is provided in both time and space. The solution is then applied to several test cases, which are taken from the previous studies. The results of this study are compared with experimental and theoretical results as well as those of the previous ones. The comparisons indicate that the Boussinesq equations solved here and in the previous study produce quite similar results. Copyright © 2006 John Wiley & Sons, Ltd.     

Dynamic response of a deepwater riser subjected to combined axial and transverse excitation by the nonlinear coupled model

In offshore engineering long slender risers are simultaneously subjected to both axial and transverse excitations. The axial load is the fluctuating top tension which is induced by the floater’s heave motion, while the transverse excitation comes from environmental loads such as waves. As the time-varying axial load may trigger classical parametric resonance, dynamic analysis of a deepwater riser with combined axial and transverse excitations becomes more complex. In this study, to fully capture the coupling effect between the planar axial and transverse vibrations, the nonlinear coupled equations of a riser’s dynamic motion are formulated and then solved by the central difference method in the time domain. For comparison, numerical simulations are carried out for both linear and nonlinear models. The results show that the transverse displacements predicted by both models are similar to each other when only the random transverse excitation is applied. However, when the combined axial dynamic tension and transverse wave forces are both considered, the linear model underestimates the response because it ignores the coupling effect. Thus the coupled model is more appropriate for deep water. It is also found that the axial excitation can significantly increase the riser’s transverse response and hence the bending stress, especially for cases when the time-varying tension is located at the classical parametric resonance region. Such time-varying effects should be taken into account in fatigue safety assessment.     

外激励作用下亚音速二维壁板的复杂响应研究

研究了亚音速流中二维壁板在外激励作用下的复杂响应问题。采用迦辽金方法将非线性运动控制方程离散为常微分方程组,采用数值方法进行计算,研究了壁板系统的复杂响应。应用最大李亚普诺夫指数和庞加莱截面方法对系统的运动性质进行了判定。结果表明,系统随着参数的变化呈现出复杂的响应,系统的周期运动与混沌运动会相间出现;系统由周期运动进...     

Numerical simulation of flows over 2D and 3D backward-facing inclined steps

The work deals with the numerical solution of incompressible turbulent flow in a channel with a backward-facing step having various inclination angles. Also, the inclination of upper wall is considered. The mathematical model is based on the Reynolds averaged Navier–Stokes equations. The governing equations are closed by the explicit algebraic Reynolds stress (EARSM) model according to Wallin and Johansson or by linear eddy viscosity models (SST, TNT k–ω). The numerical solution is carried out by the implicit finite-volume method based on the artificial compressibility and by the finite-element method amd both approaches compared. The numerical simulations use as reference the experimental data by Makiola and Driver and Seegmiller in large aspect ratio channels. In these cases, the results are obtained by 2D and 3D simulations. Further narrow channel PIV experimental data are used as reference for 3D simulations.     

Modal analysis of measurements from a large-scale VIV model test of a riser in linearly sheared flow

Large-scale model testing of a tensioned steel riser in well-defined sheared current was performed at Hanøytangen outside Bergen, Norway in 1997. The length of the model was 90 m and the diameter was 3 cm. The aim of the present work is to look into this information and try to improve the understanding of vortex-induced vibrations (VIV) for cases with very high order of responding modes, and in particular to study if and under which circumstances the riser motions would be single-mode or multi-mode. The measurement system consisted of 29 biaxial gauges for bending moment. The signals are processed to yield curvature and displacement and further to identify modes of vibration. A modal approach is used successfully employing a combination of signal filtering and least-squares fitting of precalculated mode-shapes. As a part of the modal analysis, it is demonstrated that the equally spaced instrumentation limited the maximum mode number to be extracted to be equal to the number of instrumentation locations. This imposed a constraint on the analysis of in-line (IL) vibration, which occurs at higher frequencies and involves higher modes than cross-flow (CF). The analysis has shown that in general the riser response was irregular (i.e. broad-banded) and that the degree of irregularity increases with the flow speed. In some tests distinct spectral peaks could be seen, corresponding to a dominating mode. No occurrences of single-mode (lock-in) were seen. The IL response is more broad-banded than the CF response and contains higher frequencies. The average value of the displacement r.m.s over the length of the riser is computed to indicate the magnitude of VIV motion during one test. In the CF direction the average displacement is typically 1/4 of the diameter, almost independent of the flow speed. For the IL direction the values are in the range 0.05–0.08 of the diameter. The peak frequency taken from the spectra of the CF displacement at riser midpoint show approximately to be equal to the Strouhal frequency. The peak frequency in IL direction was typically twice the Strouhal frequency.