共查询到20条相似文献，搜索用时 0 毫秒
1.
Zang Yuelong 《Acta Mechanica Solida Sinica》1997,10(1):66-75
A boundary element method for three-dimensional steady ship wave-making potentialproblems is established with the Rankine souree function as its fundamental solution.In the treatmentof the linearized free surface condition,one sided.upstream finite difference operator(FDO)is usedto suppress the upstream waves,and the equation of the disturbance velocity is established so that thefirst order FDO can be used in place of the second order FDO.Compared with the method with the sec-ond order FDO,the current method gives better precision and stability.Numerical examples are pre-sented for verification. 相似文献
2.
利用满足Laplace方程,线性化自由面条件及无穷远处条件的Havelock兴波源涵数,建立了关于常航速稳态船波势函数的边界积分方程.针对这个积分方程,建立了相应的数值计算方法,编制了一般三维问题的边界元法计算机程序,可用来计算全潜和半潜物体的稳态绕流场及船舶兴波阻力. 相似文献
3.
An accurate three‐dimensional numerical model, applicable to strongly non‐linear waves, is proposed. The model solves fully non‐linear potential flow equations with a free surface using a higher‐order three‐dimensional boundary element method (BEM) and a mixed Eulerian–Lagrangian time updating, based on second‐order explicit Taylor series expansions with adaptive time steps. The model is applicable to non‐linear wave transformations from deep to shallow water over complex bottom topography up to overturning and breaking. Arbitrary waves can be generated in the model, and reflective or absorbing boundary conditions specified on lateral boundaries. In the BEM, boundary geometry and field variables are represented by 16‐node cubic ‘sliding’ quadrilateral elements, providing local inter‐element continuity of the first and second derivatives. Accurate and efficient numerical integrations are developed for these elements. Discretized boundary conditions at intersections (corner/edges) between the free surface or the bottom and lateral boundaries are well‐posed in all cases of mixed boundary conditions. Higher‐order tangential derivatives, required for the time updating, are calculated in a local curvilinear co‐ordinate system, using 25‐node ‘sliding’ fourth‐order quadrilateral elements. Very high accuracy is achieved in the model for mass and energy conservation. No smoothing of the solution is required, but regridding to a higher resolution can be specified at any time over selected areas of the free surface. Applications are presented for the propagation of numerically exact solitary waves. Model properties of accuracy and convergence with a refined spatio‐temporal discretization are assessed by propagating such a wave over constant depth. The shoaling of solitary waves up to overturning is then calculated over a 1:15 plane slope, and results show good agreement with a two‐dimensional solution proposed earlier. Finally, three‐dimensional overturning waves are generated over a 1:15 sloping bottom having a ridge in the middle, thus focusing wave energy. The node regridding method is used to refine the discretization around the overturning wave. Convergence of the solution with grid size is also verified for this case. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
4.
采用常数边界元对船舶与流体界面进行离散，求解船舶兴波势及船舶兴波阻力。这种方法可避免在船舶与流体自由面交线上安置节点，因而避免了这些节点建立补充方程。因为满足自由面条件的Ｈａｖｅｌｏｃｋ源函数的源点和场点不能同时在自由面上，使得自由面上的节点无法用Ｈａｖｅｌｏｃｋ源函数的建立方程。如对自由面交线上的节点建立补充方程，则要对线性自由面条件中包含的未知势函数的二阶导数用差分形式表达，引入较大误差。 相似文献
5.
6.
A methodology for computing three‐dimensional interaction between waves and fixed bodies is developed based on a fully non‐linear potential flow theory. The associated boundary value problem is solved using a finite element method (FEM). A recovery technique has been implemented to improve the FEM solution. The velocity is calculated by a numerical differentiation technique. The corresponding algebraic equations are solved by the conjugate gradient method with a symmetric successive overrelaxation (SSOR) preconditioner. The radiation condition at a truncated boundary is imposed based on the combination of a damping zone and the Sommerfeld condition. This paper (Part 1) focuses on the technical procedure, while Part 2 [Finite element simulation of fully non‐linear interaction between vertical cylinders and steep waves. Part 2. Numerical results and validation. International Journal for Numerical Methods in Fluids 2001] gives detailed numerical results, including validation, for the cases of steep waves interacting with one or two vertical cylinders. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
7.
The present paper makes use of a wave equation formulation of the primitive shallow water equations to simulate one-dimensional free surface flow. A numerical formulation of the boundary element method is then developed to solve the wave continuity equation using a time-dependent fundamental solution, while an explicit finite difference scheme is used to derive velocities from the primitive momentum equation. One-dimensional free surface flows in open channels are treated and the results compared with analytical and numerical solutions. © 1997 John Wiley & Sons, Ltd. 相似文献
8.
A fully nonlinear irregular wave tank has been developed using a three‐dimensional higher‐order boundary element method (HOBEM) in the time domain. The Laplace equation is solved at each time step by an integral equation method. Based on image theory, a new Green function is applied in the whole fluid domain so that only the incident surface and free surface are discretized for the integral equation. The fully nonlinear free surface boundary conditions are integrated with time to update the wave profile and boundary values on it by a semi‐mixed Eulerian–Lagrangian time marching scheme. The incident waves are generated by feeding analytic forms on the input boundary and a ramp function is introduced at the start of simulation to avoid the initial transient disturbance. The outgoing waves are sufficiently dissipated by using a spatially varying artificial damping on the free surface before they reach the downstream boundary. Numerous numerical simulations of linear and nonlinear waves are performed and the simulated results are compared with the theoretical input waves. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
9.
A new fully non‐hydrostatic model is presented by simulating three‐dimensional free surface flow on a vertical boundary‐fitted coordinate system. A projection method, known as pressure correction technique, is employed to solve the incompressible Euler equations. A new grid arrangement is proposed under a horizontal Cartesian grid framework and vertical boundary‐fitted coordinate system. The resulting model is relatively simple. Moreover, the discretized Poisson equation for pressure correction is symmetric and positive definite, and thus it can be solved effectively by the preconditioned conjugate gradient method. Several test cases of surface wave motion are used to demonstrate the capabilities and numerical stability of the model. Comparisons between numerical results and analytical or experimental data are presented. It is shown that the proposed model could accurately and effectively resolve the motion of short waves with only two layers, where wave shoaling, nonlinearity, dispersion, refraction, and diffraction phenomena occur. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
10.
Min-Joon Kim Hie-Tae Moon Yong-Bum Lee Seok-Ki Choi Yong-Kyun Kim Ho-Yun Nam Mann Cho 《国际流体数值方法杂志》1998,28(6):887-902
In this paper, an efficient numerical method for unsteady free surface motions, with simple geometries, has been devised. Under the potential flow assumption, the governing equation of free surface flows becomes a Laplace equation, which is treated here by means of a series expansions of the velocity potential. The free surface is represented with a height function. The present method is applied to surface gravity waves to test the stability and accuracy of the method. To show the versatility of the method, a model for a dip formation is considered. © 1998 John Wiley & Sons, Ltd. 相似文献
11.
In this paper, the authors treat the free‐surface waves generated by a moving disturbance with a constant speed in water of finite and constant depth. Specifically, the case when the disturbance is moving with the critical speed is investigated. The water is assumed inviscid and its motion irrotational. The surface tension is neglected. It is well‐known that the linear theory breaks down when a disturbance is moving with the critical speed. As a remedy to overcome the invalid linear theory, approximate non‐linear theories have been applied with success in the past, i.e. Boussinesq and Korteweg de Vries equations, for example. In the present paper, the authors describe a finite element method applied to the non‐linear water‐wave problems in two dimensions. The present numerical method solves the exact non‐linear formulation in the scope of potential theory without any additional assumptions on the magnitude of the disturbances. The present numerical results are compared with those obtained by other approximate non‐linear theories. Also presented are the discussions on the validity of the existing approximate theories applied to two types of the disturbances, i.e. the bottom bump and the pressure patch on the free‐surface at the critical speed. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
12.
In this study, a method is developed to simulate the interaction between free surface flows and moving or deforming boundaries using the flux‐difference splitting scheme on the hybrid Cartesian/immersed boundary method. At each physical time step, the boundary is defined by an unstructured triangular surface grid. Immersed boundary (IB) nodes are distributed inside an instantaneous fluid domain based on edges crossing the boundary. At an IB node, dependent variables are reconstructed along the local normal line to the boundary. Inviscid fluxes are computed using Roe's flux‐difference splitting scheme for immiscible and incompressible fluids. The free surface is considered as a contact discontinuity in the density field. The motion of free surface is captured without any additional treatment along the fluid interface. The developed code is validated by comparisons with other experimental and computational results for a piston‐type wave maker, impulsive motion of a submerged circular cylinder, flow around a submerged hydrofoil, and Rayleigh–Taylor instability. The developed code is applied to simulate wave generation due to a continuously deforming bed beneath the free surface. The violent motion of a free surface caused by sloshing in a spherical tank is simulated. In this case, the free surface undergoes breakup and reconnection. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
13.
Mirmosadegh Jamali 《国际流体数值方法杂志》2006,51(4):371-383
An effective numerical technique is presented to model turbulent motion of a standing surface wave in a tank. The equations of motion for turbulent boundary layers at the solid surfaces are coupled with the potential flow in the bulk of the fluid, and a mixed BEM–finite difference technique is used to model the wave motion and the corresponding boundary layer flow. A mixing‐length theory is used for turbulence modelling. The model results are in good agreement with previous physical and numerical experiments. Although the technique is presented for a standing surface wave, it can be easily applied to other free surface problems. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
14.
Roy A. Walters 《国际流体数值方法杂志》2005,49(7):721-737
The objective of this research is to develop a model that will adequately simulate the dynamics of tsunami propagating across the continental shelf. In practical terms, a large spatial domain with high resolution is required so that source areas and runup areas are adequately resolved. Hence efficiency of the model is a major issue. The three‐dimensional Reynolds averaged Navier–Stokes equations are depth‐averaged to yield a set of equations that are similar to the shallow water equations but retain the non‐hydrostatic pressure terms. This approach differs from the development of the Boussinesq equations where pressure is eliminated in favour of high‐order velocity and geometry terms. The model gives good results for several test problems including an oscillating basin, propagation of a solitary wave, and a wave transformation over a bar. The hydrostatic and non‐hydrostatic versions of the model are compared for a large‐scale problem where a fault rupture generates a tsunami on the New Zealand continental shelf. The model efficiency is also very good and execution times are about a factor of 1.8 to 5 slower than the standard shallow water model, depending on problem size. Moreover, there are at least two methods to increase model accuracy when warranted: choosing a more optimal vertical interpolation function, and dividing the problem into layers. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
15.
16.
17.
The paper describes the implementation of moving‐mesh and free‐surface capabilities within a 3‐d finite‐volume Reynolds‐averaged‐Navier–Stokes solver, using surface‐conforming multi‐block structured meshes. The free‐surface kinematic condition can be applied in two ways: enforcing zero net mass flux or solving the kinematic equation by a finite‐difference method. The free surface is best defined by intermediate control points rather than the mesh vertices. Application of the dynamic boundary condition to the piezometric pressure at these points provides a hydrostatic restoring force which helps to eliminate any unnatural free‐surface undulations. The implementation of time‐marching methods on moving grids are described in some detail and it is shown that a second‐order scheme must be applied in both scalar‐transport and free‐surface equations if flows driven by free‐surface height variations are to be computed without significant wave attenuation using a modest number of time steps. Computations of five flows of theoretical and practical interest—forced motion in a pump, linear waves in a tank, quasi‐1d flow over a ramp, solitary wave interaction with a submerged obstacle and 3‐d flow about a surface‐penetrating cylinder—are described to illustrate the capabilities of our code and methods. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
18.
A level set formulation in a generalized curvilinear coordinate is developed to simulate the free surface waves generated by moving bodies or the sloshing of fluid in a container. The Reynolds‐averaged Navier–Stokes (RANS) equations are modified to account for variable density and viscosity in two‐phase (i.e. water–air) fluid flow systems. A local level set method is used to update the level set function and a least square technique adopted to re‐initialize it at each time step. To assess the developed algorithm and its versatility, a selection of different fluid–structure interaction problems are examined, i.e. an oscillating flow in a two‐dimensional square tank, a breaking dam involving different density fluids, sloshing in a two‐dimensional rectangular tank and a Wigley ship hull travelling in calm water. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
19.
The lubrication theory is extended for transient free‐surface flow of a viscous fluid inside a three‐dimensional thin cavity. The problem is closely related to the filling stage during the injection molding process. The pressure, which in this case is governed by the Laplace's equation, is determined using the boundary element method. A fully Lagrangian approach is implemented for the tracking of the evolving free surface. The domain of computation is the projection of the physical domain onto the (x, y) plane. This approach is valid for simple and complex cavities as illustrated for the cases of a flat plate and a curved plate. It is found that the flow behavior is strongly influenced by the shape of the initial fluid domain, the shape of the cavity, and inlet flow pressure. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
20.
An unsteady single‐phase level set RANS method is used to resolve and investigate bow wave breaking around a surface combatant advancing in calm water, including induced vortices and free surface scars. A level set free surface capturing approach was extended and combined with local overset grid refinement for resolution of complex interfacial topologies and small‐scale free surface features. Although the focus of the paper is on wave breaking at Fr=0.35, results over three speeds (Fr=0.28, 0.35, and 0.41) show that the method can accurately predict the changes in resistance and free surface topology, with the two highest speeds showing bow wave breaking. For the Fr=0.35 case, comparison of wave elevation results shows good agreement with the data, including the development and thickening of the bow wave sheet, sequential formation of two overturning plungers with reconnections, and the formation of two free surface scars at the reconnection sites. The computational fluid dynamics (CFD) solution shows a steep shoulder wave, similar to the experiment, but does not predict the experimentally observed weak spilling breaking shoulder wave. Although the current predictions converge to steady state, the region of unsteady free surface measured experimentally can be reasonably well predicted from the region of the simulation where the wave slope exceeds 17°. Comparisons of velocity components and axial vorticity at four cross planes show that the method can accurately predict the wake of low axial velocity and vortical cross flow associated with the breaking bow wave. In addition, the simulation is used to explain the initial development of the overturning bow wave, induced vortices and scars and to fill in the relatively sparse experimental data set by providing a global picture of the axial vortex structure near the free surface. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献