A Numerical Study on Water Waves Generated by A Submerged Moving Body in A Two-Layer Fluid System

This is a numerical study on the time development of surface waves generated by a submerged body moving steadily in a two-layer fluid system, in which a layer of water is underlain by a layer of viscous mud. The fully nonlinear Navier-Stokes equations are solved on FLUENT with the Volume-of-Fluid (VOF) multiphase scheme in order to simulate the free surface waves as well as the water-mud interface waves as functions of time. The numerical model is validated by mimicking a reported experiment in a one-layer ...     

A Numerical Study on Water Waves Generated by a A Submerged Moving Body in a A Two-Llayer Fluid System

This is a numerical study on the time development of surface waves generated by a submerged body moving steadily in a two-layer fluid system, in which a layer of water is underlain by a layer of viscous mud. The fully nonlinear Navier–Stokes equations are solved on FLUENT with the Volume-of-Fluid (VOF) multiphase scheme in order to simulate the free surface waves as well as the water–mud interface waves as functions of time. The numerical model is validated by mimicking a reported experiment in a one-layer system before it is applied to a two-layer system. It is found that the presence of bottom mud in a water layer can lead to large viscous damping of the surface waves. To investigate the problem systematically, the effects of the Froude number and the mud layer thickness, density and viscosity relative to those of water are evaluated and discussed in detail.     

Numerical Solution of Membrane Forces for A Free-Free Floating Plate with Large Deflection

1 .IntroductionVerylargefloatingstructures (VLFS)haveattractedconsiderableattentionrecentlyduetoitspo tentialuseintheexploitationofoceanresourcesandintheutilizationofoceanspaces (Cui,2 0 0 2 ) .Owingtoitslargesizeandrelativelysmallthickness ,theinfluenceofelasticdeformationonthehydro dynamicresponseisimportant.Thisisknownashydroelasticity .Threedimensionallinear (Wu ,1 984;PriceandWu ,1 985)andnonlinear (Chenetal.,2 0 0 3)hydroelasticityiswellappliedtotheanalysisofthreedimensionalshipsando…     

主动式声纳列阵拖曳系统姿态数值计算

主动式声纳列阵拖曳系统是用于探测潜艇的新型声纳系统，为了准确探测潜艇的位置，必须首先预报声纳列阵的瓷态，本文通过对其三维力学模型的分析，得到该系统的运动微分方程，其中缆索的力学方程是基于Ablow和Milinazzo的模型，而对于拖体则运用六自由度空间运动方程模拟，结合边界条件，用有限差分法求解，通过对拖船的不同运动状态如匀速，变速和回转的计算，证明本文的方法对于预报声纳列阵的姿态是有效的。     

各种间隙比下海底管线绕流二维数值模拟与实验研究

通过数值模拟方法与物理模型实验,对不同间隙比情况下的海底管线绕流特性进行研究。数值模拟采用有限差分法对标准k-ε湍流方程和控制方程进行离散,利用GMRES算法求解离散方程,模拟绕流流场。物理模型实验中采用超声测速仪ADV系统测取不同断面的流速分布情况,分别分析不同间隙比情况下的漩涡运动特性。同时分析各垂直断面不同间隙比时尾流流速变化情况,数值模拟结果与实验结果吻合较好,本文成果为进一步研究海底管线的防护措施打下基础。     

长兴岛海区波流相互作用数值模拟研究

波和流是近岸海区的主要动力因素。应用二维潮流数学模型和最新第三代近岸海浪模式SWAN,建立了非结构网格下二维情况近岸波流耦合作用数学模型。时间离散采用欧拉向前格式,空间离散采用有限体积法显式格式。通过将波浪场及潮流场进行迭代耦合计算,实现了对波流共同作用下波浪场和潮流场的数值模拟。将模型应用于矩形海湾实验和李孟国数模实验等理想地形以及大连长兴岛海区实际复杂地形算例,并用现场实测资料对计算结果进行验证,结果表明:耦合结果与实测结果吻合良好,并且要优于未耦合的结果。     

Dynamic Finite Element Analysis for Interaction between Two Phase Saturated Soil Foundation and Platform

In this paper, the foundation soil of offshore structure is simulated as a two phase saturated porous medium. The dynamic equations of porous medium and finite element formulation are given. For structural analysis, the technique of multilevel substructure is used, and the saturated soil analysis is set in the highest level substructure model. Based on these theories a dynamic finite element analysis program DIASS for the analysis of interaction between two phase ocean soil foundation and platform structures has been developed. A numerical example is given here to illustrate the influence of the pore water in soil on the structural response of an ocean platform.     

Numerical Study on Plume Interaction Above An Alternating Diffuser in Stagnant Water

1.IntroductionOcean outfall diffusers are used in most sewage disposal schemes to discharge the wastewater asmultiple buoyant jets which can be rapidly mixed with and diluted by the ambient current.Diffuserstake ona varietyof forms,rangingfromone withjets…     

SPH Numerical Modeling for the Wave-Thin Structure Interaction

In this paper, a numerical model of 2D weakly compressible smoothed particle hydrodynamics (WCSPH) is developed to simulate the interaction between waves and thin structures. A new color domain particle (CDP) technique is proposed to overcome difficulties of applying the ghost particle method to thin structures in dealing with solid boundaries. The new technique can deal with zero-thickness structures. To apply this enforcing technique, the computational fluid domain is divided into sub domains, i.e., boundary domains and internal domains. A color value is assigned to each particle, and contains the information of the domains in which the particle belongs to and the particles can interact with. A particle, nearby a thin boundary, is prevented from interacting with particles, which should not interact with on the other side of the structure. It is possible to model thin structures, or the structures with the thickness negligible with this technique. The proposed WCSPH module is validated for a still water tank, divided by a thin plate at the middle section, with different water levels in the subdomains, and is applied to simulate the interaction between regular waves and a perforated vertical plate. Finally, the computation is carried out for waves and submerged twin-horizontal plate interaction. It is shown that the numerical results agree well with experimental data in terms of the pressure distribution, pressure time series and wave transmission.     

Numerical Simulation of Wave Field near Submerged Bars by PLIC-VOF Model

Investigating the wave field near structures in coastal and offshore engineering is of increasing significance. In the present study, simulation is done of the wave profile and flow field for waves propagating over submerged bars, using PLIC-VOF （Pieeewise Linear Interface Construction） to trace the free surface of wave and finite difference method to solve vertical 2D Navier-Stokes （N-S） equations. A comparison of the numerical results for two kinds of submerged bars with the experimental ones shows that the PLIC-VOF model used in this study is effective and can compute the wave field precisely.     

印度-欧亚板块碰撞期间红河断裂带活动性的数值模拟研究

红河断裂带是一巨型走滑断裂带,在印度板块与欧亚板块碰撞过程中成为印支地块和华南地块之间相互运动的主位移带。通过有限差分数值模拟试验,结果表明,在碰撞期间,红河断裂带首先经历了左旋走滑运动,红河断裂带中央部分的左旋偏移量最大,达到约325km,印支地块向南东挤出,在5Ma BP后表现为右旋走滑,说明沿红河断裂带的剪切运动是红河断裂带西端与印度-欧亚碰撞边界之间距离的函数,随时间而减小并控制了印支地块与华南地块之间的相对运动。     

Numerical Simulation of Waves Interaction with A Submerged Horizontal Twin-Plate Breakwater

The numerical investigation of regular waves interacting with a submerged horizontal twin-plate breakwater is presented in this paper. A numerical model with an absorbing wave-maker is established based on the VOF method. The validity of the model is verified by experimental results. Comparisons between the numerical and experimental results show that both the water surface profiles and the wave-induced pressures can be modeled accurately. Wave deformation over the breakwater,water particle velocities aroun...     

拖曳线列阵阵形与姿态数值计算

本文提出的计算程序ＣＡＢＬＥ能预报拖缆系统中拖曳线列阵的三维运动,同时可计算各时刻缆索及列阵上相应的应力。系统的运动方程与Ａｂｌｏｗ［１］提出的相似,采用有限差分法计算求解。程序能够计算不同航速定常直航状态下列阵的阵形与姿态,可解决不定常运动（例如回转运动）、考虑流速以及船舶有升沉等多种运动状况的阵形预测。所得结果与其它计算及试验所提出的结果显示出良好的一致性。程序所花计算时间少,计算可靠。程序可运用于实际操作中。     

Numerical storm surge model with higher order finite difference method of lines for the coast of Bangladesh

In this study, the method of lines (MOLs) with higher order central difference approximation method coupled with the classical fourth order Runge-Kutta (RK(4,4)) method is used in solving shallow water equations (SWEs) in Cartesian coordinates to foresee water levels associated with a storm accurately along the coast of Bangladesh. In doing so, the partial derivatives of the SWEs with respect to the space variables were discretized with 5-point central difference, as a test case, to obtain a system of ordinary differential equations with time as an independent variable for every spatial grid point, which with initial conditions were solved by the RK(4,4) method. The complex land-sea interface and bottom topographic details were incorporated closely using nested schemes. The coastal and island boundaries were rectangularized through proper stair step representation, and the storing positions of the scalar and momentum variables were specified according to the rules of structured C-grid. A stable tidal regime was made over the model domain considering the effect of the major tidal constituent, M₂ along the southern open boundary of the outermost parent scheme. The Meghna River fresh water discharge was taken into account for the inner most child scheme. To take into account the dynamic interaction of tide and surge, the generated tidal regime was introduced as the initial state of the sea, and the surge was then made to come over it through computer simulation. Numerical experiments were performed with the cyclone April 1991 to simulate water levels due to tide, surge, and their interaction at different stations along the coast of Bangladesh. Our computed results were found to compare reasonable well with the limited observed data obtained from Bangladesh Inland Water Transport Authority (BIWTA) and were found to be better in comparison with the results obtained through the regular finite difference method and the 3-point central difference MOLs coupled with the RK(4,4) method with regard to the root mean square error values.     

Numerical Tidal Current Modeling Using Fine Boundary Fitted Grids for Tidal Flats

Fine grids with small spacing in boundary-fitted coordinates are normally used to treat thecomputation of fluid dynamics for estuarine areas and tidal flats.However,the adoption of Cartesian com-ponents of velocity vectors in this kind of non-orthogonal coordinates will definitely result in a diffculty insolving implicitly the transformed momentum equations,and also complicate the wet-dry point judgementused for flood areas.To solve this problem,equations in terms of generalized contravariant velocityvectors in curvilinear coordinates are derived in the present study,by virtue of which,anAlternative-Direction-Implicit numerical scheme in non-orthogonal grids would then be easily obtained,and wet-dry point judgement would as well be largely simplified.A comparison is made between the explic-it scheme and implicit scheme,showing that the present model is accurate and numerically stable for com-putations of fluid dynamics for estuarine areas and tidal flats.     

A Numerical and Experimental Study on the Hull-Propeller Interaction of A Long Range Autonomous Underwater Vehicle

Range is an important factor to the design of autonomous underwater vehicles(AUVs), while drag reduction efforts are pursued, the investigation of body-propeller interaction is another vital consideration. We present a numerical and experimental study of the hull-propeller interaction for deeply submerged underwater vehicles, using a proportionalintegral-derivative(PID) controller method to estimate self-propulsion point in CFD environment. The hydrodynamic performance of hull and propeller at the balance state when the AUV sails at a fixed depth is investigated, using steady RANS solver of Star-CCM+. The proposed steady RANS solver takes only hours to reach a reasonable solution. It is more time efficient than unsteady simulations which takes days or weeks, as well as huge consumption of computing resources. Explorer 1000, a long range AUV developed by Shenyang Institute of Automation, Chinese Academy of Sciences, was studied as an object, and self-propulsion point, thrust deduction,wake fraction and hull efficiency were analyzed by using the proposed RANS method. Behind-hull performance of the selected propeller MAU4-40, as well as the hull-propeller interaction, was obtained from the computed hydrodynamic forces. The numerical results are in good qualitative and quantitative agreement with the experimental results obtained in the Qiandao Lake of Zhejiang province, China.     

基于双层Boussinesq方程的三维波浪数值模型及其验证

Liu等给出的最高导数为2的双层Boussinesq水波方程具有较好的色散性和非线性,基于该方程建立了有限差分法的三维波浪数值模型。在矩形网格上对方程进行了空间离散,采用高阶导数近似方程中的时、空项,时间积分采用混合4阶Adams-Bashforth-Moulton的预报—校正格式。模拟了深水条件下的规则波传播过程,计算波面与解析结果吻合较好,反映出数值模型能很好地刻画波面过程及波面处的速度变化;在kh=2π条件下可较为准确获得沿水深分布的水平和垂向速度,这与理论分析结果一致。最后,利用数值模型计算了规则波在三维特征地形上的传播变形,数值结果和试验数据吻合较好;高阶非线性项会对波浪数值结果产生一定的影响,当波浪非线性增强,水深减少将产生更多的高次谐波。建立的双层Boussinesq模型对强非线性波浪的演化具有较好的模拟精度。     

The Boundary Element Method for Pile-Soil Interaction

In this paper, several mathmatical models for the pile- soil interaction are outlined. The Boundary Element Method is one of the very effective methods for the reasonable models of elasticity and elastoplasticity. The major of this paper is concerned with the Boundary Element Method for the pile-soil interaction, including general methods and calculating formulation of static and dynamic analysis of the pile and pile groups. Some results of analysis are also given.     

Study on the Interaction of Water Waves with Semi-Circular Breakwater

—The present study investigates the interaction of steep waves with semi-circular breakwaterwith the complex plane＇s Cauchy boundary integral theorem.The boundary integral method is used totransform the calculation in fluid domain into its boundary alone.In the calculation the computation do-main is moved with the propagation of waves.A numerical solution is obtained for incident Stokes wavespassing the submerged obstacles.This method has been extended to the calculation of wave run-up on aslope for estimating wave overtopping.     

一种基于表面状态参数的三维海洋动力学模型数值求解方法*

大面积的海洋表面实时运动状态可以通过遥感技术获取, 而海洋内部的运动状态只能进行定点观测, 无法达到大面积的实时监测。基于海洋动力学基本原理, 在正压浅海大陆架模式下, 在三维空间构建海洋表面与内部运动状态的关系模型; 利用遥感探测的海洋表面流速与流向数据, 结合海域的浪高、风速状态参数, 运用有限差分法反演出深层海流的流速与流向信息。反演结果符合海洋动力学规律, 反映出了深层流整体分布状态, 扩展了雷达遥感应用范围。