Surface effects of internal wave generated by a moving source in a two-layer fluid of finite depth

IntroductionInternalwaves,whichcanbeexcitedbymanykindsofdisturbancesfromthesurface,bottomorinteriorofastratifiedocean ,suchasthewindstrain ,theflowoverunevenbottomsandthemovingbodyatthesurfaceorunderwater,areremarkablyobviousasthereexistsapycnoclineorathermoclineintheocean[1,2 ] .ThepropagationofgravitysurfacewaveindeepwatergeneratedbyasteadilymovingdisturbanceisrestrictedinaV_shapedregionwithahalfangleθc =19.5°,whichiswell_knownKelvinshipwave[3 ] .Thewaveinducedbyamovingunderwaterobjectisa…     

Experiments on the interaction of internal waves with the semi-submersible platform in a stratified fluid

对密度分层流体中内波与半潜平台的相互作用问题 进行了模型试验研究. 采用摇板方法进行了内波造波试验,对内波波长、周期和波高进行了 测量分析,获得了内波波高和波长与周期之间的相关关系. 利用激光和倾角仪方法对半潜平 台的纵荡和纵摇运动响应进行了测量分析,获得了平台纵荡运动及纵摇角幅值与内波周期之 间的相关关系. 结果表明,在半潜平台的设计与应用中,内波对其运动响应的影响是不可忽 视的. 特别地,发现了在内波周期的某个范围内, 半潜平台的纵荡和纵摇都会出现 倍频周期运动响应的现象.     

分层流体中内孤立波在台阶上的反射和透射

基于匹配渐近展开和格林函数的方法,研究了两层流体系统中内孤立波在台阶地形上透射、 反射及其分裂的演化特征. 通过保角变换和求解奇异Fredholm积分方程,获得了反映地形 效应对Boussinesq方程影响的约化边界条件,藉此建立了KdV演化方程的``初值'问题, 根据散射反演理论获得了反射波和透射波的解析表达式. 分析结果表明：上下流体层的厚度 比、密度比以及台阶高度对于反射和透射波振幅及其分裂具有显著的影响. 尤其当上层流体 厚度小于下层厚度时,由于存在临界点,在其附近反射波的幅值随台阶高度的演化由单调增 变为单调减,透射波的幅值由单调减变为单调增;上台阶的反射波与入射波反相,其最大幅 值可达到入射波的数倍;此外,下台阶反射波也可发展为单支孤立波,它区别于单层流体中 反射波仅为衰减的振荡波列.     

Exact pressure integrations on submerged bodies in waves using a quadtree adaptive mesh algorithm

The development of an adaptive free surface, mesh cutting, methodology, in order to analytically integrate pressures on varying wet parts of partially submerged surfaces in the presence of waves, is presented. Given a function of free‐surface elevation, the algorithm checks for the intersection of the body with the free surface and, based on user‐defined parameters, modifies the initial mesh, by subdividing the elements where necessary and eliminating others, via a quadtree approach. Redundant sub‐divisions, generated in the quad‐division process, are partially eliminated, but the quadrilateral nature of the elements is always kept. The free‐surface function must be single‐valued and its definition domain simply connected. Hydrostatic and Froude–Krylov forces are computed exactly on each panel by means of analytical formulations, which are derived and presented, based on the theory of linear gravity waves and from applying Green's theorem. Copyright © 2014 John Wiley & Sons, Ltd.     

连续分层海洋中内波传播的一种数值模式

基于Euler方程,使用有限体积法建立了一种密度为连续分层情况下、适应水深变化的水域中内波传播的数值模式.为了使计算格式能够达到二阶精度,对流项的处理使用了TVD (total variation diminishing)格式.将SIMPLE算法引入连续分层海洋中内波的数值计算,为了简化计算并方便地适应多种TVD格式,在计算预估速度场时采用了显式格式,而没有采用传统的隐式格式;鉴于在原始的SIMPLE算法中没有涉及到由于密度扰动而引起的静水压力场的改变问题,给出了该问题的计算方法.因此改进了SIMPLE算法.出流边界的处理采用阻尼消波和Sommerfeld辐射条件相结合的方式,以使内波得到有效的衰减和释放.将等水深水域的数值解和理论解进行了比较,两者吻合较好;并对存在潜堤时数值计算的不同时刻密度变化的空间分布进行了详细的定性分析.计算结果表明,所建立的数值模式能有效地模拟内波的传播和变形.      

EXPERIMENTAL INVESTIGATION ON THE EVOLUTION OF INTERNAL SOLITARY WAVE OVER A RIDGE AND ITS ACTING FORCES ON A SUBMERGED SLENDER BODY

内孤立波沿山脊地形传播过程中的结构变化使得影响海洋结构物及水下航行器安全的不确定因素增加,用大型重力式分层流试验水槽系统研究内孤立波过山脊地形及其对细长潜体作用力特性,可有效提高对其危害性机理的认识.为此,利用电导率探头阵列,结合染色标识方法,测量内孤立波演变特性,同时利用三分量传感器测量水下细长体模型受力特性.研究结果表明：山脊地形显著改变下凹型内孤立波结构,表现为坡前波幅增大,坡顶背风波面抬升,坡后波长变长;内波结构变化直接影响细长体受力特性,表现为坡前向下作用力明显增强,坡顶向上作用力突显,坡后作用力持续时间变长;细长体受力特性变化影响其运动趋势,坡前下沉运动增强、朝山脊方向纵荡运动显著.     

Comparisons of CIP,compact and CIP‐CSL2 schemes for reproducing internal solitary waves

Six different models were evaluated for reproducing internal solitary waves which occur and propagate in a stratified flow field with a sharp interface. Three stages were used to compute internal solitary waves in a stratified field: (1) first‐phase computation of momentum equations, (2) second‐phase computation of momentum equations, which corresponds to computing the Poisson's equation, and (3) density computation. The six models discussed in this paper consisted of combinations of four different schemes, a three‐point combined compact difference scheme (CCD), a normal central difference scheme (CDS), a cubic‐polynomial interpolation (CIP), and an exactly conservative semi‐Lagrangian scheme (CIP‐CSL2). The residual cutting method was used to solve the Poisson's equation. Three tests were used to confirm the validity of the computations using KdV theory; i.e. the incremental wave speed and amplitude of internal solitary waves, the maximum horizontal velocity and amplitude, and the wave form. In terms of the shape of an internal solitary wave, using CIP for momentum equations was found to provide better performance than CCD. These results suggest one of the most appropriate scheme for reproducing internal solitary waves may be one in which CIP is used for momentum equations and CCD to solve the Poisson's equation. Copyright © 2005 John Wiley & Sons, Ltd.     

Wave drag of an ellipsoid of revolution rapidly moving in a horizontal direction in a uniformly stratified fluid

The problem of the wave drag of ellipsoids moving in a uniformly stratified ideal fluid is considered by means of modeling the bodies by surface distributions of mass sources. Analytical results are obtained using the distributions known from the theory of a uniform fluid, which make it possible to describe the emission of internal waves by rapidly moving ellipsoids of revolution (spheroids) in the limit of large Froude numbers. An asymptotically simplified form of the dependence of the wave drag on the Froude number and the spheroid axis ratio is found. In the particular case of a sphere, the result confirmed earlier by Greenslade by making comparisons with a numerical calculation and experimental data is obtained.     

两层流体中水波在垂直薄板上的反射与透射

研究在两层流体中表面波模态和内波模态的波浪与半潜式刚性垂直薄板 相互作用的问题. 基于特征函数展开理论,建立了两种模态入射波作用下,半潜式刚性垂直 薄板的反射与透射能量的计算方法,证明了对每一种模态的入射波,另一种模态波浪的反射 与透射能量是相等的. 对水面漂浮和座底半潜式薄板的反射与透射能量,以及作用在其上的 水平波浪力进行了数值计算分析,表明在某个频率范围内,流体的分层效应对这些水动力 量的影响是不可忽视的. 特别地,当薄板的一端位于两层流体的内界面上时,两种模态波浪 的能量转化是最大的.     

An SPH model for free surface flows with moving rigid objects

This paper presents a computational model for free surface flows interacting with moving rigid bodies. The model is based on the SPH method, which is a popular meshfree, Lagrangian particle method and can naturally treat large flow deformation and moving features without any interface/surface capture or tracking algorithm. Fluid particles are used to model the free surface flows which are governed by Navier–Stokes equations, and solid particles are used to model the dynamic movement (translation and rotation) of moving rigid objects. The interaction of the neighboring fluid and solid particles renders the fluid–solid interaction and the non‐slip solid boundary conditions. The SPH method is improved with corrections on the SPH kernel and kernel gradients, enhancement of solid boundary condition, and implementation of Reynolds‐averaged Navier–Stokes turbulence model. Three numerical examples including the water exit of a cylinder, the sinking of a submerged cylinder and the complicated motion of an elliptical cylinder near free surface are provided. The obtained numerical results show good agreement with results from other sources and clearly demonstrate the effectiveness of the presented meshfree particle model in modeling free surface flows with moving objects. Copyright © 2013 John Wiley & Sons, Ltd.     

An adaptive algorithm of precise integration for transient analysis

This paper presents an improved precise integration algorithm for transient analysis of heat transfer and some other problems. The original precise integration method is improved by means of the inverse accuracy analysis so that the parameterN, which has been taken as a constant and an independent parameter without consideration of the problems in the original method, can be generated automatically by the algorithm itself. Thus, the improved algorithm is adaptive and the accuracy of the algorithm is not dependent on the length of the time step in the integration process. It is shown that the numerical results obtained by the method proposed are more accurate than those obtained by the conventional time integration methods such as the difference method and others. Four examples are given to demonstrate the validity, accuracy and efficiency of the new method. Project supported by the National Natural Science Foundation of China (No. 19872016, 19872017), the National Key Basic Research Special Foundation (G1999032805) and the Foundation for University Key Teachers by the Ministry of Education of China.     

An adaptive quenching algorithm for a nonlinear single-degree-of-freedom system

An indirect adaptive quenching algorithm for a nonlinear single-degree-of-freedom system with unknown constant system parameters is presented. The system is subject to external or parametric sinusoidal disturbances and the resulting control signal is also sinusoidal. The quenching algorithm provides a reduction in the control effort required compared to direct disturbance cancellation. The disturbance sinusoid and the unknown parameters are incorporated into the system model and an extended Kalman filter (EKF) with modified update equations is used to estimate the system state and parameters. The estimates are then used to form the quenching signal. The adaptive quenching algorithm is found to work well inside a quenching region defined by the separatrices and suggests the use of a hybrid control law. The algorithm was verified by implementing it on an analog computer.     

An adaptive artificial viscosity for the displacement shallow water wave equation

The numerical oscillation problem is a difficulty for the simulation of rapidly varying shallow water surfaces which are often caused by the unsmooth uneven bottom,the moving wet-dry interface, and so on. In this paper, an adaptive artificial viscosity(AAV) is proposed and combined with the displacement shallow water wave equation(DSWWE) to establish an effective model which can accurately predict the evolution of multiple shocks effected by the uneven bottom and the wet-dry interface. The effec...     

An adaptive reanalysis method for genetic algorithm with application to fast truss optimization

Although the genetic algorithm （GA） for structural optimization is very robust, it is very computationally intensive and hence slower than optimality criteria and mathematical programming methods. To speed up the design process, the authors present an adaptive reanalysis method for GA and its applications in the optimal design of trusses. This reanalysis technique is primarily derived from the Kirsch＇s combined approximations method. An iteration scheme is adopted to adaptively determine the number of basis vectors at every generation. In order to illustrate this method, three classical examples of optimal truss design are used to validate the proposed reanalysis-based design procedure. The presented numerical results demonstrate that the adaptive reanalysis technique affects very slightly the accuracy of the optimal solutions and does accelerate the design process, especially for large-scale structures.     

An adaptive genetic algorithm for solving bilevel linear programming problem

Bilevel linear programming,which consists of the objective functions of the upper level and lower level,is a useful tool for modeling decentralized decision problems. Various methods are proposed for solving this problem.Of all the algorithms,the ge- netic algorithm is an alternative to conventional approaches to find the solution of the bilevel linear programming.In this paper,we describe an adaptive genetic algorithm for solving the bilevel linear programming problem to overcome the difficulty of determining the probabilities of crossover and mutation.In addition,some techniques are adopted not only to deal with the difficulty that most of the chromosomes may be infeasible in solving constrained optimization problem with genetic algorithm but also to improve the efficiency of the algorithm.The performance of this proposed algorithm is illustrated by the examples from references.     

具有乘性噪声和随机量测时滞的目标跟踪算法

针对在机动目标跟踪过程中产生的乘性量测噪声与随机量测时滞问题,提出了一种改进的高斯滤波(GF)算法,并给出了该算法的一种具体实现形式——随机时滞和乘性噪声容积卡尔曼滤波器(CKF-RDMN)。首先,利用一组满足伯努利分布的随机序列描述随机出现的量测时滞现象。其次,利用乘性噪声满足高斯条件分布的特性,改进量测更新过程。最后依据三阶球径容积法则,对高斯积分进行求解。由于该算法是对经典GF算法的改进,因此,也可作为通用的滤波框架解决上述问题。通过与推广无迹卡尔曼滤波和推广扩展卡尔曼滤波对比,仿真结果表明,CKF-RDMN在解决乘性量测噪声与随机量测时滞问题时,具有更高的估计精度。     

Multiphase computational method for thermal interactions between compressible fluid and arbitrarily shaped solids

A numerical method was investigated for multiphase fields consisting of compressible gas and arbitrarily shaped solids. Since the proposed model is based on a one‐fluid model in which variables are averaged according to the phase fractions in the computational cells; it enables us to estimate gas‐solid momentum and thermal interactions without setting up adapting grids even if the solids have extremely complicated shapes. The governing equations are derived with the characteristics of an ideal gas assuming the specific heat to be uniform in the multiphase field. The derived equations in conservative form are discretized with a finite volume method. In addition, the pressure is calculated implicitly in a similar way to incompressible flow solvers. Because of these improvements, the proposed method allows us to calculate low Mach number compressible flows free from the Courant‐Friedrichs‐Lewy condition based on the speed of sound and to conserve the mass more accurately. To confirm the validity of the proposed method, it was applied to natural convection around an isothermal cylinder and a heat‐conducting pipe. In comparison with previous studies, it was confirmed that the gas flows and temperature distributions are predicted reasonably. In addition, a numerical experiment was conducted under more complicated conditions, namely, gas leaking from a container including heat sources. As a result, it was demonstrated that the proposed method enables us to predict unsteady variations of pressure and temperature distributions in the container due to the leakage while still conserving mass accurately.     

具有小密度差的两层流体中运动点源的二阶内波解

在具有自由面的两层流体中，运动点源生成的Kelvin船波存在两种模式，即表面波模式和内波模式。当上、下层流体密度比趋于1时，由内波模式计算的界面波幅趋于无穷大，这与实验事实相违背。为克服此困难，在自由面和界面作小波幅运动的假设，引入一个小密度差参数。研究了运动点源在无粘、不可压且具有小密度差的两层有限深流体中生成的高阶波动。首先利用摄动方法推导了各阶小参数满足的边值问题；其次，给出了小密度差情形下的可解性条件。证明了在密度比趋于1的极限情形，不存在导致界面波幅无穷大的内波模式；最后，利用Phillips的非线性共振相互作用理论，构造了具有自由面的两层有限深流体中Kelvin船波系的二阶一致有效波动解，并证明了该解在深水情形下退化为Newman关于均匀流体中自由面的二阶波动解。     

一种冲击波作用下结构毁伤算法研究

针对爆炸冲击波与建筑物结构相互作用过程,分析了冲击波与结构碎块作用机理,发展了一种能够模拟建筑物结构破坏及冲击波传播过程的计算模型和方法。采用建筑物结构工程毁伤载荷作为判据,处理结构在冲击波作用下的破坏问题;利用流固耦合界面算法处理结构运动引起的泄压效应,利用“虚拟网格通气技术”处理结构碎块对冲击波的阻碍作用,模拟了冲击波作用下典型建筑物的毁伤过程及冲击波传播过程。结果表明,该模型在模拟冲击波与结构的作用过程中,压力计算结果与非结构动网格模拟结果符合较好;在典型建筑物毁伤过程的数值模拟中,计算得到的建筑物毁伤效果和冲击波超压分布与建筑物物理毁伤特点符合。

     

二维矩形腔体中混合流体行进波对流

Rayleigh-Benard模型是研究对流稳定性,时空结构和非线性特性的典型模型之一。本文的兴趣集中在二维矩形腔体中混合流体对流场的结构方面。利用SIMPLE算法数值求解流体力学方程组,模拟了充分发展的二维矩形腔体中混合流体对流。结果说明偏离传导失去稳定的系统经过亚临界分叉产生了振动对流。进一步,我们给出了分叉曲线及其沿分叉曲线的上部分支三个Rayleigh数对应的对流图案的垂直速度场,流线,温度场,浓度场和Shadowgraph强度的等值线图。所有场的结构分析表明浓度场及Shadowgraph强度的等值线图可以很好的特征行进波的运动特性。