Oblique water entry of a wedge into waves with gravity effect

The hydrodynamic problem of a two dimensional wedge entering waves with gravity effect is analysed based on the incompressible velocity potential theory. The problem is solved through the boundary element method in the time domain. The stretched coordinate system in the spatial domain, which is based on the ratio of the Cartesian system in the physic space to the vertical distance the wedge has travelled into the water, is adopted based on the consideration that the decay of the effect of the impact away from the body is proportional to this ratio. The solution is sought for the total potential which includes both the incident and disturbed potentials, and decays towards the incident potential away from the body. A separate treatment at initial stage is used, in which the solution for the disturbed potential is sought to avoid the very large incident potential amplified by dividing the small travelled vertical distance of the wedge. The auxiliary function method is used to calculate the pressure on the body surface. Detailed results through the free surface elevation and the pressure distribution are provided to show the effect of the gravity and the wave, and their physical implications are discussed.     

Fully nonlinear modeling of radiated waves generated by floating flared structures

The nonlinear radiated waves generated by a structure in forced motion, are simulated numerically based on the potential theory. A fully nonlinear numerical model is developed by using a higher-order boundary element method （HOBEM）. In this model, the instantaneous body position and the transient free surface are updated at each time step. A Lagrangian technique is employed as the time marching scheme on the free surface. The mesh regridding and interpolation methods are adopted to deal with the possible numerical instability. Several auxiliary functions are proposed to calculate the wave loads indirectly, instead of directly predicting the temporal derivative of the velocity potential. Numerical experiments are carried out to simulate the heave motions of a submerged sphere in infinite water depth, the heave and pitch motions of a truncated flared cylinder in finite depth. The results are verified against the published numerical results to ensure the effectiveness of the proposed model. Moreover, a series of higher harmonic waves and force components are obtained by the Fourier transformation to investigate the nonlinear effect of oscillation frequency. The difference among fully nonlinear, body-nonlinear and linear results is analyzed. It is found that the nonlinearity due to free surface and body surface has significant influences on the numerical results of the radiated waves and forces.     

三维非稳定渗流自由面边界积分项的精确数值计算

利用固定网格法分析三维非稳定渗流问题时,将要面对两项积分难题:以自由面及单元表面为边界的空间积分及以自由面为边界的曲面积分。针对常用的任意8结点6平面三维普通单元,提出采用坐标变换及等参变换技术求取空间积分项的精确数值解;至于曲面积分项,建议改用单元非饱和区部分表面作为积分边界,经过坐标变换及等参变换处理积分边界后,利用高斯数值积分可求出曲面积分项的精确数值解。通过一个普通单元及一项均质半无限边界堤坝的实例分析,表明此方法的精确性和稳定性良好。     

楔形体入水时域解的复边界元数值分析

基于速度势理论,利用复数变量边界元法对二维楔形体常速入水冲击的时域解进行了数值研究。 以相似解作为数值计算的初始条件,采用时域解射流线性近似处理方法,利用复数变量边界元法进行求解,以 减少计算量和数值误差。深入讨论了扩展坐标系求时域解、射流处理、网格划分和网格更新等关键技术。最 后数值计算了不同斜升角楔形体入水时的自由液面隆起、射流飞溅和压力分布,经与相似解结果作比较,自由 液面隆起轮廓基本吻合,而压力分布更符合实际情况,从而证明了模型及分析方法的正确性。     

Finite difference method for 2-D and 3-D nonlinear free surface wave problems

A finite difference method is developed for the numerical modelling of the 2-D and 3-D unsteady potential flow generated by transient disturbances on the free surface, on which the nonlinear boundary conditions are fully satisfied. The unknown function is computed with an iteration scheme processing in a transformed time-invariant space. After the velocity is calculated, the location of the free surface is renewed and so is the value of velocity on it. The boundary-value problem of the governing equation is then solved at the next time step. The present method incorporates the FFT. Consequently, a tri-diagonal equation system is obtained which could be readily solved. The feasibility of this method has been demonstrated by 2-D and 3-D examples corresponding to different initial disturbances. This work is supported by the science foundation of Academia Sinica. The paper had been accepted by the XVIth International Congress of IUTAM, Lyngby, Denmark, August, 1984.     

A computational method for simulating transient motions of an incompressible inviscid fluid with a free surface

A new numerical method has been developed for the analysis of unsteady free surface flow problems. The problem under consideration is formulated mathematically as a two-dimensional non-linear initial boundary value problem with unknown quantities of a velocity potential and a free surface profile. The basic equations are discretized spacewise with a boundary element method and timewise with a truncated forward-time Taylor series. The key feature of the present paper lies in the method used to compute the time derivatives of the unknown quantities in the Taylor series. The use of the Taylor series expansion has enabled us to employ a variable time-stepping method. The size of time increment is determined at each time step so that the remainders of the truncated Taylor series should be equal to a given small error limit. Such a variable time-stepping technique has made a great contribution to numerically stable computations. A wave-making problem in a two-dimensional rectangular water tank has been analysed. The computational accuracy has been verified by comparing the present numerical results with available experimental data. Good agreement is obtained.     

基于单元破裂的岩石裂纹扩展模拟方法

传统离散元方法在处理破裂问题时, 采用界面上的准则进行判断, 裂纹只能沿着单元边界扩展. 当物理问题存在宏观或微观裂隙时, 在界面上应用准则具有其合理性; 而裂纹沿着单元边界扩展, 使得裂纹路径受网格影响较大, 扩展方向受到限制. 针对上述情况, 可以基于单元破裂的方式, 构建连续- 非连续单元法, 并应用于岩石裂纹扩展问题的模拟. 该方法在连续计算时, 将单元离散为具有物理意义的弹簧系统, 在局部坐标系下由弹簧特征长度、面积求解单元变形和应力, 通过更新局部坐标系和弹簧特征量, 可进一步计算块体大位移、大转动, 连续问题计算结果与有限元一致, 同时提高了计算效率. 在此基础上, 引入最大拉应力与莫尔—库伦的复合准则, 判断单元破裂状态和破裂方向, 并采用局部块体切割的方式, 在单元内形成初始裂纹. 裂纹两侧相应增加新的计算节点, 同时引入内聚力模型描述裂纹两侧的法向、切向作用与张开度及滑移变形之间的关系. 按此方式, 裂纹尖端处的扩展路径可穿过单元内部和单元边界, 在扩展方向的选取上更为准确. 最后, 通过三点弯曲梁、单切口平板拉伸、双切口试样等典型数值试验, 模拟裂纹在拉伸、压剪等各种应力状态下的扩展问题, 并对岩石单轴压缩试验的破坏过程进行模拟, 分析裂纹形成与应力—应变曲线各阶段之间的对应关系. 结果表明: 连续—非连续单元法通过单元内部破裂的方式, 可以显示模拟裂纹萌生、扩展、贯通直至形成宏观裂缝的过程.      

自由活塞压缩管ALE方法数值模拟

当前国际上实现高焓气体流动的实验手段之一是自由活塞驱动类脉冲设备,包括自由活塞激波风洞和自由活塞膨胀管.采用自由活塞压缩管作为激波风洞和膨胀管的驱动段时,其驱动能力在很大程度上决定了该类设备的性能.本文采用计算流体力学中任意拉格朗日——欧拉方法(arbitrary Lagrangian Eulerian)数值模拟了压缩管内部的自由活塞运动和气体流动特征.采用移动网格技术来适应活塞运动边界,耦合求解网格运动和气体流动过程,并通过双时间步长方法进行流体运动的时间积分.为了满足几何守恒律(geometric conservation law),对移动网格的法向矢量和表面面积计算进行了修正.不同时刻的活塞位置试验测量结果及欧拉方法预测结果,以及基于简单波理论获得的运动活塞底部气体压力、活塞速度与活塞位置都与当前的ALE方法十分一致.该工作为下一步数值模拟自由活塞激波风洞和自由活塞膨胀管中包括压缩管、激波管和喷管等不同部位的耦合流动提供了基础.      

成层半空间出平面自由波场的一维化时域算法

提出了一种计算出平面SH波斜入射时弹性水平成层半空间中自由波场时域计算的一维化有 限元方法. 在进行有限元网格划分时,竖向单元取满足有限元模拟精度的任意尺寸,水平向 网格尺寸由时间离散步长和水平视波速确定,并自动进行虚拟网格划分. 基底设置人工边界, 并将波动输入转化为等效荷载施加在边界节点上. 然后将集中质量有限元法和中心差分法相 结合建立节点运动方程,并将水平方向相邻节点的运动用该节点相邻时刻的运动表示,从而 将求解节点运动的二维方程组转化为一维方程组. 求解此方程组,即得到自由场中竖向一列 节点的运动. 最后根据行波传播的特点,可方便地确定全部自由波场. 理论分析和数值算例 表明,该方法具有较高的精度和良好的稳定性.     

A finite element method for free surface flows of incompressible fluids in three dimensions. Part I. Boundary fitted mesh motion

Computational fluid mechanics techniques for examining free surface problems in two‐dimensional form are now well established. Extending these methods to three dimensions requires a reconsideration of some of the difficult issues from two‐dimensional problems as well as developing new formulations to handle added geometric complexity. This paper presents a new finite element formulation for handling three‐dimensional free surface problems with a boundary‐fitted mesh and full Newton iteration, which solves for velocity, pressure, and mesh variables simultaneously. A boundary‐fitted, pseudo‐solid approach is used for moving the mesh, which treats the interior of the mesh as a fictitious elastic solid that deforms in response to boundary motion. To minimize mesh distortion near free boundary under large deformations, the mesh motion equations are rotated into normal and tangential components prior to applying boundary conditions. The Navier–Stokes equations are discretized using a Galerkin–least square/pressure stabilization formulation, which provides good convergence properties with iterative solvers. The result is a method that can track large deformations and rotations of free surface boundaries in three dimensions. The method is applied to two sample problems: solid body rotation of a fluid and extrusion from a nozzle with a rectangular cross‐section. The extrusion example exhibits a variety of free surface shapes that arise from changing processing conditions. Copyright © 2000 John Wiley & Sons, Ltd.     

水泵叶片表面三维边界层计算

本文用有限差分法计算混流式可逆水力机械水泵工况叶片表面的三维边界层。水泵叶轮中主流区的三维势流由直接边界元法计算。对于叶片面附近的粘性流动。用三维半正交坐标系中的边界层方程表示。为了提高计算精度采用贴体坐标技术生成边界层区域的计算网格。并利用Cebeci等变换函数及Keller差分格式离散方程。用分块解法求解。计算叶轮叶片表面的压力分布与相应试验结果进行了对比。     

二维弹性结构入水冲击过程中的流固耦合效应

描述了一个研究弹性结构入水冲击过程中水弹性效应的数值方法,在弹性结构入水冲击过程中,流体域作用在结构上的水动力载荷由边界元法获得,而结构的弹性动力响应则由有限元方法求解,通过线性给离散Ｂｅｒｎｏｕｌｌｉ方程将有限元方程和边界元方程耦合到一起,从而获得了求解流场和结构动力响应的相互耦合的运动方程。在数值考虑了自由表面的非线性边界条件,通过引入射流单元以及最大射流厚度,较好地处理了冲击引起的射流问题。     

均匀流中双色波传播特性的模拟研究

针对双色波浪与均匀流相互作用问题,采用时域高阶边界元方法建立自由水面满足完全非线性边界条件的数学模型。求解中采用混合欧拉-拉格朗日方法追踪流体瞬时水面,运用四阶龙格库塔方法更新下一时间步的波面和速度势。通过与已发表试验结果对比,验证了本模型的准确性。通过数值计算研究了水流参数对各组成波及衍生的高阶波幅值、波浪和水流间能量交换的影响规律。     

Simulation of liquid sloshing in curved‐wall containers with arbitrary Lagrangian–Eulerian method

There are many challenges in the numerical simulation of liquid sloshing in horizontal cylinders and spherical containers using the finite element method of arbitrary Lagrangian–Eulerian (ALE) formulation: tracking the motion of the free surface with the contact points, defining the mesh velocity on the curved wall boundary and updating the computational mesh. In order to keep the contact points slipping along the curved side wall, the shape vector in each time advancement is defined to modify the kinematical boundary conditions on the free surface. A special function is introduced to automatically smooth the nodal velocities on the curved wall boundary based on the liquid nodal velocities. The elliptic partial differential equation with Dirichlet boundary conditions can directly rezone the inner nodal velocities in more than a single freedom. The incremental fractional step method is introduced to solve the finite element liquid equations. The numerical results that stemmed from the algorithm show good agreement with experimental phenomena, which demonstrates that the ALE method provides an efficient computing scheme in moving curved wall boundaries. This method can be extended to 3D cases by improving the technique to compute the shape vector. Copyright © 2007 John Wiley & Sons, Ltd.     

潜堤上波流传播的完全非线性数值模拟

利用时域高阶边界元方法建立模拟波流混合作用的完全非线性数值水槽模型, 其中自由水面满足完全非线性自由水面条件. 采用混合欧拉-拉格朗日方法追踪流体瞬时水面, 运用4阶Runge-Kutta方法更新下一时间步的波面和速度势. 为了减少计算域, 提高计算速度, 采用同时消除底面和侧面的镜像格林函数; 在每一时间步内, 对网格进行重新划分以避免由于网格运动变形而引起的数值不稳定问题. 对水流中淹没潜堤上的波浪变形在水槽中开展了物理模型试验, 并把试验结果和数值结果进行了对比, 吻合得很好. 进一步研究了水流及潜堤的存在对高阶谐波产生的影响.      

An arbitrary Lagrangian Eulerian (ALE) formulation for free surface flows using the characteristic‐based split (CBS) scheme

An arbitrary Lagrangian Eulerian (ALE) method for non‐breaking free surface flow problems is presented. The characteristic‐based split (CBS) scheme has been employed to solve the ALE equations. A simple mesh smoothing procedure based on coordinate averaging (Laplacian smoothing) is employed in the calculations. The mesh velocity is calculated at each time step and incorporated as part of the scheme. Results presented show an excellent agreement with the available experimental data. Copyright © 2005 John Wiley & Sons, Ltd.     

Solution of incompressible flows with or without a free surface using the finite volume method on unstructured triangular meshes

An incompressible Navier–Stokes solver based on a cell‐centre finite volume formulation for unstructured triangular meshes is developed and tested. The solution methodology makes use of pseudocompressibility, whereby the convective terms are computed using a Godunov‐type second‐order upwind finite volume formulation. The evolution of the solution in time is obtained by subiterating the equations in pseudotime for each physical time step, with the pseudotime step set equal to infinity. For flows with a free surface the computational mesh is fitted to the free surface boundary at each time step, with the free surface elevation satisfying a kinematic boundary condition. A ‘leakage coefficient’, ε, is introduced for the calculation of flows with a free surface in order to control the leakage of flow through the free surface. This allows the assumption of stationarity of mesh points to be made during the course of pseudotime iteration. The solver is tested by comparing the output with a wide range of documented published results, both for flows with and without a free surface. The presented results show that the solver is robust. © 1999 John Wiley & Sons, Ltd.     

A finite element analysis of a free surface drainage problem of two immiscible fluids

A sharp interface problem arising in the flow of two immiscible fluids, slag and molten metal in a blast furnace, is formulated using a two-dimensional model and solved numerically. This problem is a transient two-phase free or moving boundary problem, the slag surface and the slag–metal interface being the free boundaries. At each time step the hydraulic potential of each fluid satisfies the Laplace equation which is solved by the finite element method. The ordinary differential equations determining the motion of the free boundaries are treated using an implicit time-stepping scheme. The systems of linear equations obtained by discretization of the Laplace equations and the equations of motion of the free boundaries are incorporated into a large system of linear equations. At each time step the hydraulic potential in the interior domain and its derivatives on the free boundaries are obtained simultaneously by solving this linear system of equations. In addition, this solution directly gives the shape of the free boundaries at the next time step. The implicit scheme mentioned above enables us to get the solution without handling normal derivatives, which results in a good numerical solution of the present problem. A numerical example that simulates the flow in a blast furnace is given.     

An assessment of a parallel,finite element method for three‐dimensional,moving‐boundary flows driven by capillarity for simulation of viscous sintering

A parallel, finite element method is presented for the computation of three‐dimensional, free‐surface flows where surface tension effects are significant. The method employs an unstructured tetrahedral mesh, a front‐tracking arbitrary Lagrangian–Eulerian formulation, and fully implicit time integration. Interior mesh motion is accomplished via pseudo‐solid mesh deformation. Surface tension effects are incorporated directly into the momentum equation boundary conditions using surface identities that circumvent the need to compute second derivatives of the surface shape, resulting in a robust representation of capillary phenomena. Sample results are shown for the viscous sintering of glassy ceramic particles. The most serious performance issue is error arising from mesh distortion when boundary motion is significant. This effect can be severe enough to stop the calculations; some simple strategies for improving performance are tested. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical simulation of the hydrodynamic interaction between a sedimenting particle and a neutrally buoyant particle

A new method for the simulation of the translational and rotational motions of a system containing a sedimenting particle interacting with a neutrally buoyant particle has been developed. The method is based on coupling the quasi-static Stokes equations for the fluid with the rigid body equations of motion for the particles. The Stokes equations are solved at each time step with the boundary element method. The stresses are then integrated over the surface of each particle to determine the resultant forces and moments. These forces and moments are inserted into the rigid body equations of motion to determine the translational and rotational motions of the particles. Unlike many other simulation techniques, no restrictions are placed on the shape of the particles. Superparametric boundary elements are employed to achieve accurate geometric representations of the particles. The simulation method is able to predict the local fluid velocity, resolve the forces and moments exerted on the particles, and track the particle trajectories and orientations.