The interactions between wave-currents and offshore structures with consideration of fluid viscosity

Study of the flow field around the large scale offshore structures under the action of waves and viscous currents is of primary importance for the scouring estimation and protection in the vicinity of the structures. But very little has been known in its mechanism when the viscous effects is taken into consideration. As a part of the efforts to tackle the problem, a numerical model is presented for the simulation of the flow field around a fixed vertical truncated circular cylinder subjected to waves and viscous currents based on the depth-averaged Reynolds equations and depth-averagedk-ɛ turbulence model. Finite difference method with a suitable iteration defect correct method and an artificial open boundary condition are adopted in the numerical process. Numerical results presented relate to the interactions of a pure incident viscous current with Reynolds numberRe=10⁵, a pure incident regular sinusoidal wave, and the coexisting of viscous current and wave with a circular cylinder, respectively. Flow fields associated with the hydrodynamic coefficients of the fixed cylinder, as well as corresponding free surface profiles and wave amplitudes, are discussed. The present method is found to be relatively straightforward, computationally effective and numerically stable for treating the problem of interactions among waves, viscous currents and bodies. The project supported by the National Natural Science Foundation of China and Foundation of State Key Laboratory of Ocean Engineering at Shanghai Jiao Tong University.     

Boundary element method for wave trapping by a multi-layered trapezoidal breakwater near a sloping rigid wall

This study examines the multiple layers in a rubble mound breakwater and their effect on reflection and dissipation of incoming ocean waves. The numerical model is developed using multi-domain boundary element method for oblique water wave trapping near a sloping wall by a multi-layered trapezoidal porous structure, which is utilized to model armour, filter and core layers while examining the hydrodynamics in different configurations. Both, the constant element and linear element approaches to boundary element method are discussed. The cases of bottom-standing porous structures as being submerged and fully extended are considered. The wave hydrodynamics over the structure is described by the reflection and dissipation coefficients along with the forces acting on the sloping wall, and is influenced by wave and structural parametrics of the system. The influence of armour layer in different configurations is highlighted for various structural and wave parameters.

     

Finite element method for shallow water equation including open boundary condition

A method to deal with an open boundary condition in the analysis of water surface waves, the tide, etc. by means of the finite element method is proposed in this paper. The present method has two important features relating to the treatment of the open boundary condition. The first feature is to consider the non-reflective virtual boundary condition which has been developed in the numerical wave analysis method. The incident wave conditions without spurious reflected waves can be imposed at the open boundary. The second feature is to identify the amplitude of the components of incident waves in terms of observed water elevations in the field of standing waves. This can be done as a parameter identification based on an optimization problem by applying the conjugate gradient method. The applicability of this method to wave propagation problems is verified by several numerical computations.     

直立码头前船波浪力耦合计算模型

建立了外域用Boussinesq方程、内域用刚体运动方程的直立码头前二维船剖面波浪力的时 域计算耦合模型，内域与外域在交界面的匹配条件是流量连续和压力相等. 进行了相关模型 实验，并把计算结果与实验结果进行了对比. 推导了船体与水底和直立码头之间间隙内流体 运动的自振频率，研究了间隙内流体运动的共振现象.     

港口非线性波浪耦合计算模型研究

建立了外域用差分法求解高阶Boussinesq方程、内域用边界元法求解Laplace方程的二维船 非线性波浪力时域计算的耦合模型. 研究了该类耦合模型的匹配条件、耦合求解过程和内域、 外域公共区域长度的确定. 该耦合模型计算结果与只用边界元求解Laplace方程模型的计算 结果和实验结果对比表明，该耦合模型不仅计算精度高，而且计算效率快，适用于研究较大 区域内波浪对物体的非线性作用.     

Variational approach to non-linear boundary value problems for elasto-plastic incompressible bending plate

Non-linear boundary value problems for inelastic isotropic homogeneous incompressible bending plate, within the range of J₂-deformation theory, are considered. An existence of the weak solution of the non-linear problem with clamped boundary condition is obtained in H²(Ω) by using monotone operator theory and Browder-Minty theorem. For linearization of the non-linear problem a monotone iteration scheme is constructed. It is shown that the sequence of potentials obtained from the sequence of approximate solutions (i.e. iterations), is a monotone decreasing one. Convergence of the iteration process in H²-norm is proved by using the convexity argument. Numerical solutions, based on finite-difference scheme, are given for linear bending problems with rigid clamped as well as simply supported boundary conditions. Further numerical examples are presented to illustrate the convergence of approximate solutions and monotonicity of the potentials as applied to the non-linear problems.     

Development of a numerical wave tank for analysis of nonlinear and irregular wave field

A numerical wave-absorption filter has been developed for an open boundary condition in the analysis of nonlinear and irregular wave evolution. The filter is composed of a simulated sponge layer and Sommerfeld's radiation condition at the outer edge of the layer. The wave-absorption characteristics of the filter have been investigated by applying the linear potential theory and a two-dimensional nonlinear boundary element model. In both cases, the filter is found to he applicable for a wide range of wave parameters. In order to realize an idealized “numerical wave tank”, the present model also incorporates a nonreflective wave generator in the computational domain composed of a series of vertically aligned point sources. Numerous numerical experiments demonstrate that the present approach is effective in generating an arbitrary wave profile without reflection not only at the open boundaries but also at the wave generator.     

On flows behind shock waves reflected from a solid wall

The flow fields associated with the interaction of a normal shock wave with a plane wall kept at a constant temperature were studied based on kinetic theory which can describe appropriately the shock structure and its reflection process. With the use of a difference scheme, the time developments of the distributions of the fluid dynamic quantities (velocity, temperature, pressure and number density of the gas) were obtained numerically from the BGK model of the Boltzmann equation subject to the condition of diffusive-reflection at the wall for several cases of incident Mach number:M ₁=1.2, 1.5, 2.0, 3.0, 4.0, 5.0 and 6.0. The reflection process of the shocks is shown explicitly together with the resulting formation of the flow fields as time goes on. The nonzero uniform velocity toward the wall occurring between the viscous boundary layer and the reflected shock wave is found to be fairly large, the magnitude of which is of the order of several percent of the velocity induced behind the incident shock, decreasing as the incident Mach number increases. It is also seen that a region of positive velocity (away from the wall) within the viscous boundary layer manifests itself in the immediate vicinity of the wall, which is distinct for larger incident Mach numbers. Some of the calculated density profiles are compared with available experimental data and also with numerical results based on the Navier-Stokes equations. The agreement between the three results is fairly good except in the region close to the wall, where the difference in the conditions of these studies and the inappropriateness of the Navier-Stokes equations manifest themselves greatly in the gas behavior.This article was processed using Springer-Verlag TEX Shock Waves macro package 1990.     

推广的多次透射边界：标量波清形

利用空间外推将多次透射边界条件推广以提高这一边界在时域计算中的精度。对全波场仅由入射波构成及由入射波和充分发育的反射波构成两种情况导出了边界反射系数的频域公式。从时域计算精度角度分析了边界反射系数以阐明在透射大角度入射行进波动和消减运动方面推广边界的优点。     

A meshless numerical wave tank for simulation of nonlinear irregular waves in shallow water

Time domain simulation of the interaction between offshore structures and irregular waves in shallow water becomes a focus due to significant increase of liquefied natural gas (LNG) terminals. To obtain the time series of irregular waves in shallow water, a numerical wave tank is developed by using the meshless method for simulation of 2D nonlinear irregular waves propagating from deep water to shallow water. Using the fundamental solution of Laplace equation as the radial basis function (RBF) and locating the source points outside the computational domain, the problem of water wave propagation is solved by collocation of boundary points. In order to improve the computation stability, both the incident wave elevation and velocity potential are applied to the wave generation. A sponge damping layer combined with the Sommerfeld radiation condition is used on the radiation boundary. The present model is applied to simulate the propagation of regular and irregular waves. The numerical results are validated by analytical solutions and experimental data and good agreements are observed. Copyright © 2008 John Wiley & Sons, Ltd.     

Contact Interaction of the Faces of a Rectangular Crack under Normally Incident Tension–Compression Waves

A three-dimensional problem on the contact interaction between the faces of a rectangular crack under a normally incident harmonic tension–compression wave is considered. The problem is solved by using the method of boundary integral equations and an iterative algorithm. The contact forces and the discontinuity in the displacement of the crack faces are studied. The results obtained are compared with those for a finite plane crack.     

Reflection of longitudinal displacement wave at the viscoelastically supported boundary of micropolar half-space

The reflection of elastic waves at the boundary of micropolar half-space with the viscoelastic support is studied in this paper. The spring-dashpot model is used to model the viscoelastic support. The boundary condition includes the force stress, couple stress, the displacements and the micro-rotation and the contribution from the spring and the dashpot on them. The amplitude ratios and phase shifts of reflection waves with respect to the incident wave are obtained from the visco-elastic boundary conditions. Further, the energy flux ratios of the reflection waves to the incident wave are estimated and the energy flux conservation with consideration of the energy dissipation of visco-elastic boundary is used to validate the numerical results. Based on the numerical results, the influences of spring and dashpot are studied respectively. It is found that the elastically supported boundary and the viscously supported boundary have evident different influences on the amplitude ratio and the phase shift. The causes resulting in these deviations are related with the instantaneous elasticity of elastic boundary and the time-delay effects of viscous boundary.     

轴承非线性油膜力的一种变分近似解

基于自由边值理论和凸集上的变分方法，提出一种求解当轴颈大扰动时实际轴承瞬态油膜力的近似公式。公式中引入一个参数来模拟油膜破裂自由边界（雷诺边界条件），则凸集上的泛函极值问题就转换为求此参数的代数极值问题。对椭圆轴承在轴颈作大范围扰动情况下的计算结果表明，这一方法达到了很高的精度，可用于转子－轴承系统的非线性动力分析，能大量降低数值计算瞬态油膜力所需的计算量。     

Wave induced thermal boundary layers in a compressible fluid: analysis and numerical simulations

The response of a semi-infinite compressible fluid to a step-wise change in temperature of its boundary is investigated analytically and numerically. Numerical results of the boundary layer structure are compared with Clarke’s analytical solution for a gas with thermal conductivity proportional to temperature. To avoid unwanted numerical dissipation in the numerical analysis, the space-time conservation element and solution element (CESE) method has been adopted to solve the unsteady 1-D Navier-Stokes equations. Good agreement between analytical and numerical results has been found for the development of the thermal boundary layer on a long time scale. Weak shock waves and expansion waves induced by the thermal boundary layer due to its compressibility, are observed in the numerical simulation. Finally, the numerical method has been applied to the reflection of a non-linear expansion wave and to a shock wave from an isothermal wall, thereby illustrating the effect of the boundary layer on the external flow field.     

Formation and evolution of vortex rings induced by interactions between shock waves and a low-density bubble

The deformation and instability of a low-density spherical bubble induced by an incident and its reflected shock waves are studied by using the large eddy simulation method. The computational model is firstly validated by experimental results from the literature and is further used to examine the effect of incident shock wave strength on the formations and three-dimensional evolutions of the vortex rings. For the weak shock wave case (Ma?=?1.24), the baroclinic effect induced by the reflected shock wave is the key mechanism for the formation of new vortex rings. The vortex rings not only move due to the self-induced effect and the flow field velocity, but also generate azimuthal instability due to the pressure disturbance. For the strong shock wave case (Ma?=?2.2), a boundary layer is formed adjacent to the end wall owing to the approach of vortex ring, and unsteady separation of the boundary layer near the wall results in the ejection and formation of new vortex rings. These vortex rings interact in the vicinity of the end wall and finally collapse to a complicated vortex structure via azimuthal instability. For both shock wave strength cases, the evolutions of vortex rings due to the instability lead to the formation of the complicated structure dominated by the small-scale streamwise vortices.     

Theory of diffraction of weak shock waves

A numerical solution is considered to the universal nonlinear boundary-value diffraction problem which occurs in various problems of weak interaction [1, 2] in the asymptotic analysis of the flow in a region with large gradients of the parameters near the point of intersection of the incident, diffracted, and reflected waves. The analytical solutions to this type of problem usually approximately satisfy the conditions on the diffracted front, the position of which is not known beforehand, but is found along with the solution. In the present paper, the problem is solved by the numerical method of [3], which reduces the initial boundary-value problem for the system of short-wave equations with an unknown boundary to the solution of a series of boundary-value problems with a fixed boundary. The problem of the diffraction of a weak shock wave on a wedge with a finite apex angle is considered as an application of the solution. The data calculated by the asymptotic theory agree significantly better with the experimental data [5] than the theoretical data of [4].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6. pp. 176–178, November–December, 1984.     

Wave-in-deck experiments with focussed waves into a solid deck

This paper examines impact forces resulting from wave-in-deck processes from two separate series of experiments: one with a generic solid deck model, and the other with a combined jacket and deck model, both were conducted with and without an I-beam grillage in-place below the solid deck. A range of inundation levels from 2.1 to 7.1 cm at 1:80 scale is considered. The focus is on global impact forces, which are considered more relevant for integrity assessment of overall bottom-founded structures when survivability is in question and local slamming is not addressed. The objective is to characterise the resultant impact forces as well as to investigate whether there is any interaction between the flows through the jacket and hitting the deck. Focussed wave groups were generated to impinge on the models which were suspended from a carriage over a towing tank. The motion of the support carriage mimics uniform current in-line with the incident waves. Both undisturbed surface elevations as well as impact force time histories were measured. From the first series of tests, a large increase in peak forces as well as high frequency oscillations (force spikes) is observed with the grillage in-place. As soon as the jacket model is in-place for the second series of tests, albeit with a different mounting support arrangement, such a large difference vanishes, which could likely be due to the effect of frequency-dependent transfer functions. We provide experimental evidence of the presence of the jacket in modifying the wave impact on the deck through a significant reduction in the total horizontal impulse. The effects of current on the wave impact forces are also investigated. A simple analytical model based on a momentum argument is used to describe the scaling of horizontal peak force with currents and inundation levels. Finally, the importance of the short duration force spikes as well as vertical impact loads on a real structure at full-scale is discussed based on the same analytical model.     

The non-linear,dynamic response of an impulsively loaded circular plate with a central hole

A numerical method, called Direct Analysis, is described and applied to solve the problem of a plate undergoing a large impulsive load. For generality, an expanded, non-linear form of the equations of motion is used and shear correction and rotatory inertia are considered. The wave speeds are calculated from the non-linear equations and appropriate boundary conditions are applied so that reflected waves are included. The results for two types of step loading pulses are presented and compared with previously presented solutions. The response of the plate is discussed and conclusions as to the effects of the non-linearities are given.     

Variable density bore interaction with block obstacles

A horizontally variable density flow model is used to simulate hydraulic bore interactions with idealised urban obstacles. The 2D non-linear shallow water equations are solved using a second-order Monotonic Upstream-centered Schemes for Conservation Laws-Hancock Godunov-type HLLC approximate Riemann scheme. Validation test results are reported for wave propagation over a hump, a constant-density circular dam break and two 1D dam breaks involving different spatial distributions of solute concentration. Detailed parameter studies are then considered for hydraulic bore interactions with single and multiple-square obstacles under subcritical, critical and supercritical flow conditions. In all cases, reflected and diffracted wave patterns are generated immediately after the bore impacts the obstacle(s). Later, the incident bore reconstitutes itself downstream of the obstacle(s). Variable density flows are also considered, with the upstream volumetric concentrations set to values corresponding to water–sediment mixture densities of 1165 and 1495 kg/m³. It is found that the upstream Froude number, gap spacing between obstacles and upstream to downstream density difference influence the strength of the bore–structure interaction, run-up at the front face of the obstacle(s), and subsequent wave–wave interactions.     

Relection of elastic waves at the elastically supported boundary of a couple stress elastic half-space

The relection elastic waves at the elastically supported boundary of a couple stress elastic half-space are studied in this paper. Different from the classical elastic solid, there are three kinds of elastic waves in the couple stress elastic solid, and two of them are dispersive. The boundary conditions of a couple stress elastic half-space include the couple stress vector and the rotation vector which disappear in the classical elastic solids. These boundary conditions are used to obtain a linear algebraic equation set, from which the amplitude ratios of relection waves to the incident wave can be determined. Then, the relection coeficients in terms of energy lux ratios are calculated numerically, and the normal energy lux conservation is used to validate the numerical results. Based on these numerical results,the inluences of the boundary parameters, which relect the mechanical behavior of elastic support, on the relection energy partition are discussed. Both the incident longitudinal wave(the P wave) and incident transverse wave(the SV wave) are considered.