共查询到20条相似文献，搜索用时 219 毫秒
1.
The Incremental Differential Quadrature Method (IDQM) was applied to a tidal and surge model of the Bristol Channel, UK. The method is considered as an alternative new numerical technique in the field of marine hydraulics and its performance was examined through this case study. For validation of the simulated results, tide gauge data along the Bristol Channel was used. Another well known 1D model (MIKE11) and a quasi-3D model (POLCOMS) provided more insight into the flow pattern of the estuary and additional validation data. Similar to MIKE11, IDQM is unconditionally stable and so time steps of around 45 min achieved good results for the Bristol Channel, whereas for methods which are restricted to the CFL criterion (e.g. explicit finite differencing schemes), the time step is limited to around 3 min. Since there is no stability constraint in IDQM, the time step must be selected with reference to accuracy rather than stability. The usefulness of IDQM was also demonstrated by using a small number of grid points (11 along the 110 km length of the Bristol Channel) to produce accurate results. Based on the results of this case study, it is concluded that IDQM can be successfully implemented for 1D modelling of tidal elevations and surges in non-prismatic irregular channels. 相似文献
2.
A one-dimensional high-resolution finite volume model capable of simulating storm waves propagating in the coastal surf zone and overtopping a sea wall is presented. The model (AMAZON) is based on solving the non-linear shallow water (NLSW) equations. A modern upwind scheme of the Godunov-type using an HLL approximate Riemann solver is described which captures bore waves in both transcritical and supercritical flows. By employing a finite volume formulation, the method can be implemented on an irregular, structured, boundary-fitted computational mesh. The use of the NLSW equations to model wave overtopping is computationally efficient and practically flexible, though the detailed structure of wave breaking is of course ignored. It is shown that wave overtopping at a vertical wall may also be approximately modelled by representing the wall as a steep bed slope. The AMAZON model solutions have been compared with analytical solutions and laboratory data for wave overtopping at sloping and vertical seawalls and good agreement has been found. The model requires more verification tests for irregular waves before its application as a generic design tool. 相似文献
3.
A spatial fixed σ-coordinate is used to transform the Navier–Stokes equations from the sea bed to the still water level. In the fixed σ-coordinate system only a very small number of vertical grid points are required for the numerical model. The time step for using the spatial fixed σ-coordinate is efficiently larger than that of using a time dependent σ-coordinate, as there is substantial truncation error involved in the time dependent σ-coordinate transformation. There is no need to carry out the σ-coordinate transformation at each time step, which can reduce computational times. It is important that wave breaking can be potentially modeled in the fixed σ-coordinate system, but in a time-dependent σ-coordinate system the wave breaking cannot be modeled. A projection method is used to separate advection and diffusion terms from the pressure terms in Navier–Stokes equations. The pressure variable is further separated into hydrostatic and hydrodynamic pressures so that the computer rounding errors can be largely avoided. In order to reduce computational time of solving the hydrodynamic pressure equation, at every time step the initial pressure is extrapolated in time domain using computed pressures from previous time steps, and then corrected in spatial domain using a multigrid method. For each time step, only a few of iterations (typically six iterations) are required for solving the pressure equation. The model is tested against available experimental data for regular and irregular waves and good agreement between calculation results and the measured data has been achieved. 相似文献
4.
The transformation of nonlinear water waves over porous beds is studied by applying a numerical model based on Chen's [2006. Fully nonlinear Boussinesq-type equations for waves and currents over porous beds. Journal of Engineering Mechanics, 132:2, 220–230] Boussinesq-type equations for highly nonlinear waves on permeable beds. The numerical model uses a high-order time-marching solution and fourth-order finite-difference schemes for discretization of first-order spatial derivatives to obtain a computational accuracy consistent with the model equations. By forcing the wave celerity and spatial porous-damping rate of the linearized model to match the exact linear theory for horizontal porous bed over a prescribed range of relative depths, the values of the model parameters are optimally determined. Numerical simulations of the damped wave propagation over finite-thickness porous layer demonstrate the accuracy of both the numerical model and governing equations, which have been shown by prior theoretical analyses to be accurate for both nominal and thick porous layers. These simulations also elucidate on the significance of the higher-order porous-damping terms and the influence of the hydraulic parameters. Application of the model to the simulation of the wave field around a laboratory-scale submerged porous mound provides a measure of its capability, as well as useful insight into the scaling of the porous-resistance coefficients. For application to heterogeneous porous beds, the assumption of weak spatial variation of the porous resistance is examined using truncated forms of the governing equations. The results indicate that the complete set of Boussinesq-type equations is applicable to porous beds of nonhomogeneous makeup. 相似文献
5.
Unsteady nonlinear wave motions on the free surface in shallow water and over slopes of various geometries are numerically simulated using a finite difference method in rectangular grid system. Two-dimensional Navier–Stokes equations and the continuity equation are used for the computations. Irregular leg lengths and stars are employed near the boundaries of body and free surface to satisfy the boundary conditions. Also, the free surface which consists of markers or segments is determined every time step with the satisfaction of kinematic and dynamic free surface conditions. Moreover, marker-density method is also adopted to allow plunging jets impinging on the free surface. Either linear or Stokes wave theory is employed for the generation of waves on the inflow boundary. For the simulation of wave breaking phenomena, the computations are carried out with various wave periods and sea bottom slopes in surf zone. The results are compared with other existing computational and experimental results. Agreement between the experimental data and the computation results is good. 相似文献
6.
一个两时间层分裂显格式海洋环流模式(MASNUM)及其检验 总被引：1，自引：0，他引：1
HAN Lei 《海洋学报(英文版)》2014,33(11):11-35
A two-time-level, three-dimensional numerical ocean circulation model（named MASNUM） was established with a two-level, single-step Eulerian forward-backward time-differencing scheme. A mathematical model of large-scale oceanic motions was based on the terrain-following coordinated, Boussinesq, Reynolds-averaged primitive equations of ocean dynamics. A simple but very practical Eulerian forward-backward method was adopted to replace the most preferred leapfrog scheme as the time-differencing method for both barotropic and baroclinic modes. The forward-backward method is of second-order of accuracy, computationally efficient by requiring only one function evaluation per time step, and free of the computational mode inherent in the three-level schemes. This method is superior to the leapfrog scheme in that the maximum time step of stability is twice as large as that of the leapfrog scheme in staggered meshes thus the computational efficiency could be doubled. A spatial smoothing method was introduced to control the nonlinear instability in the numerical integration. An ideal numerical experiment simulating the propagation of the equatorial Rossby soliton was performed to test the amplitude and phase error of this new model. The performance of this circulation model was further verified with a regional（northwest Pacific） and a quasi-global（global ocean simulation with the Arctic Ocean excluded） simulation experiments. These two numerical experiments show fairly good agreement with the observations. The maximum time step of stability in these two experiments were also investigated and compared between this model and that model which adopts the leapfrog scheme. 相似文献
7.
O. SaoutP. Ananthakrishnan 《Applied Ocean Research》2011,33(2):158-167
A theoretical methodology to determine the open-loop directional stability of a near-surface underwater vehicle is presented. It involves a solution of coupled sway and yaw equations of motion in a manner similar to that carried out for surface ships. The stability derivatives are obtained numerically through simulation of motions corresponding to planar motion mechanism (PMM) model tests. For the numerical simulation, a boundary-integral method based on the mixed Lagrangian-Eulerian formulation is developed. The free-surface effect on the vehicle stability is determined by comparing the results with that obtained for vehicle motion in infinite fluid. The methodology was used to determine the stability of the Florida Atlantic University’s Ocean EXplorer (OEX) AUV. The presence of the free surface, through radiation damping, is found to suppress unsteady oscillations and thereby enhance the directional stability of the vehicle. With effects of free surface, forward speed, location and geometry of rudders, location of the center of gravity etc. all being significant factors affecting stability, a general conclusion cannot be drawn on their combined effect on the vehicle stability. The present computational methodology is therefore a useful tool to determine an underwater vehicle’s stability for a given configuration and thus the viability of an intended mission a priori. 相似文献
8.
Non-linear evolution of uni-directional focussed wave-groups on a deep water: A comparison of models
Up until the point at which ocean waves break, their dynamics are generally assumed to be accurately modelled by potential flow theory. For practical and computational reasons it is often useful to approximate the full potential flow solution with bandwidth and amplitude limited equations. A approximation used for waves on deep water is the Broad-banded Modified Non-linear Schrödinger equation (also known as the modified Dysthe equation). In this paper we compare this approximate model with potential flow simulations of focussing uni-directional wave-groups. We find that for moderate non-linearity the approximate model predicts very similar changes to the potential flow model. However, one of the dominant non-linear changes to the wave-group is a localised increase in the bandwidth and contraction in physical length, and beyond a certain point the approximate model fails to accurately reproduce this causing other elements, such as the maximum wave amplitude, to be poorly modelled. This modelling inaccuracy occurs in cases where, based on the initial conditions of the simulation, the approximate model would be expected to be accurate. 相似文献
9.
The paper presents an Incompressible Smoothed Particle Hydrodynamics (ISPH) method to simulate wave interactions with a porous medium. The SPH method is a mesh free particle modeling approach that is capable of tracking the large deformation of free surfaces in an easy and accurate manner. The ISPH method employs a strict incompressible hydrodynamic formulation to solve the fluid pressure and the numerical solution is obtained by using a two-step semi-implicit scheme. The ISPH flow model solves the unsteady 2D Navier–Stokes (NS) equations for the flows outside the porous media and the NS type model equations for the flows inside the porous media. The presence of porous media is considered by including additional friction forces into the equations. The developed ISPH model is first validated by the solitary and regular waves damping over a porous bed and the solitary wave interacting with a submerged porous breakwater. The convergence of the method and the sensitivity of relevant model parameters are discussed. Then the model is applied to the breaking wave interacting with a breakwater covered with a layer of porous materials. The computational results demonstrate that the ISPH flow model could provide a promising simulation tool in coastal hydrodynamic applications. 相似文献
10.
The regional oceanic modeling system (ROMS): a split-explicit, free-surface, topography-following-coordinate oceanic model 总被引：54，自引：1，他引：54
The purpose of this study is to find a combination of optimal numerical algorithms for time-stepping and mode-splitting suitable for a high-resolution, free-surface, terrain-following coordinate oceanic model. Due to mathematical feedback between the baroclinic momentum and tracer equations and, similarly, between the barotropic momentum and continuity equations, it is advantageous to treat both modes so that, after a time step for the momentum equation, the computed velocities participate immediately in the computation of tracers and continuity, and vice versa, rather than advancing all equations for one time step simultaneously. This leads to a new family of time-stepping algorithms that combine forward–backward feedback with the best known synchronous algorithms, allowing an increased time step due to the enhanced internal stability without sacrificing its accuracy. Based on these algorithms we design a split-explicit hydrodynamic kernel for a realistic oceanic model, which addresses multiple numerical issues associated with mode splitting. This kernel utilizes consistent temporal averaging of the barotropic mode via a specially designed filter function to guarantee both exact conservation and constancy preservation properties for tracers and yields more accurate (up to second-order), resolved barotropic processes, while preventing aliasing of unresolved barotropic signals into the slow baroclinic motions. It has a more accurate mode-splitting due to redefined barotropic pressure-gradient terms to account for the local variations in density field, while maintaining the computational efficiency of a split model. It is naturally compatible with a variety of centered and upstream-biased high-order advection algorithms, and helps to mitigate computational cost of expensive physical parameterization of mixing processes and submodels. 相似文献
11.
The three-dimensional numerical model with σ-coordinate transformation in the vertical direction is applied to the simulation of surface water waves and wave-induced laminar boundary layers. Unlike most of the previous investigations that solved the simplified one-dimensional boundary layer equation of motion and neglected the interaction between boundary layer and outside flow, the present model solves the full Navier–Stokes equations (NSE) in the entire domain from bottom to free surface. A non-uniform mesh system is used in the vertical direction to resolve the thin boundary layer. Linear wave, Stokes wave, cnoidal wave and solitary wave are considered. The numerical results are compared to analytical solutions and available experimental data. The numerical results agree favorably to all of the experimental data. It is found that the analytical solutions are accurate for both linear wave and Stokes wave but inadequate for cnoidal wave or solitary wave. The possible reason is that the existing analytical solutions for cnoidal and solitary waves adopt the first-order approximation for free stream velocity and thus overestimate the near bottom velocity. Besides velocity, the present model also provides accurate results for wave-induced bed shear stress. 相似文献
12.
波浪作用下缆船拖带系统非线性运动数值模拟 总被引：1，自引：0，他引：1
基于船舶操纵性运动方程和拖缆的三维动力学运动方程,提出了被拖点位置匹配的方法,建立了拖船—拖缆—被拖船系统整体非线性拖带动力学模型。为了考察被拖船航向稳定性与横向稳性的关系以及波浪载荷作用的影响,被拖船采用水平面四自由度运动方程,并引入了波浪的作用力和力矩。拖船采用PD控制方法较真实地模拟了拖船航向改变的运动过程。对一个拖船—拖缆—被拖船系统(5 000 t的拖船和3 000 t的被拖船)在时域内进行了规则波浪作用下拖带运动的模拟,计算结果表明被拖带船舶在波浪中运动呈现运动稳定、不稳定和临界状态3种可能的特性。根据模拟计算结果,认为波浪中拖带航向稳定是被拖带船舶保持稳性的必要条件。 相似文献
13.
《Ocean Modelling》1999,1(2-4):71-80
Ocean general circulation models (OGCMs) which represent the governing equations on a finite difference grid require shorter time steps with increasing resolution. Thus, until now, in the absence of filtering, the time step length has been determined by the smallest grid spacing within the model domain. Here we present a method for reducing the time step length (and increasing the number of time steps taken) at selected points in the grid, so as to minimise the computational cost of integrating the OGCM, whilst achieving numerical stability throughout the model domain without filtering. This variable time stepping method can be used to overcome numerical constraints associated with the convergence of longitude–latitude grids at the poles, and also to allow efficient integration of model domains with variable resolution. Examples of the computational saving are given. 相似文献
14.
New laboratory and field data are presented on fluid advection into the swash zone. The data illustrate the region of the inner surf zone from which sediment can be directly advected into the swash zone during a single uprush, which is termed the advection length. Experiments were conducted by particle tracking in a Lagrangian reference frame, and were performed for monochromatic breaking waves, solitary bores, non-breaking solitary waves and field conditions. The advection length is normalised by the run-up length to give an advection ratio, A, and different advection ratios are identified on the basis of the experimental data. The data show that fluid enters the swash zone from a region of the inner surf zone that can extend a distance seaward of the bore collapse location that is approximately equal to half of the run-up length. This region is about eight times wider than the region predicted by the classical swash solution of Shen and Meyer [Shen, M.C., Meyer, R.E., 1963. Climb of a bore on a beach. Part 3. Runup. Journal of Fluid Mechanics 16, 113–125], as illustrated by Pritchard and Hogg [Pritchard, D., Hogg, A.J., 2005. On the transport of suspended sediment by a swash event on a plane beach. Coastal Engineering 52, 1–23]. Measured advection ratios for periodic waves show no significant trend with Iribarren number, consistent with self-similarity in typical swash flows. The data are compared to recent characteristic solutions of the non-linear shallow water wave (NLSW) equations and both finite difference and finite volume solutions of the NLSW equations. 相似文献
15.
16.
In this paper, a numerical model is established for estimating the wave forces on a submerged horizontal circular cylinder. For predicting the wave motion, a set of two-dimensional Navier-Stokes equations is solved numerically with a finite element method. In order to track the moving non-linear wave surface boundary, the Navier-Stokes equations are discretized in a moving mesh system. After each computational time step, the mesh is modified according to the changed wave surface boundary. In order to stabilize the numerical procedure, a three-step finite element method is applied in the time integration. The water sloshing in a tank and wave propagation over a submerged bar are simulated for the first time to validate the present model. The computational results agree well with the analytical solution and the experimental data.Finally, the model is applied to the simulation of interaction between waves and a submerged horizontal circular cylinder.The effects of the KC number and the cylinder depth on the wave forces are studied. 相似文献
17.
18.
多向不规则波浪的确定性模拟 总被引：1，自引：0，他引：1
波浪波动时间过程及波列的模拟，对于开展实际波浪对于工程建筑物的作用具有重要的意义。本文采用线性叠加的单叠加模型，建立了多向不规则波浪的确定性模拟方法。基于理论模拟的规则波、单向不规则波和多向不规则波，验证了波浪确定性模拟方法的有效性。定性地对比分析了模拟波列和已知波列的一致性；定量地研究了模拟波浪在空间范围rr/Ls的误差分布情况(rr表示指定位置与给定位置的空间距离，Ls为有效波长)。并且建议，采用本文方法进行波浪确定性模拟时，最佳的浪高仪间距应小于0.12Ls。 相似文献
19.
非线性波浪时域计算的三维耦合模型 总被引：3，自引：1，他引：2
将计算区域Ω划分为内域Ω1和外域Ω2(Ω2=Ω-Ω1),外域控制方程采用改进线性频散特性的二维Boussinesq方程,用预报一校正法数值求解;结构物附近的内域控制方程为三维Navier-Stokes方程,由VOF方法数值求解。通过在外域和内域相匹配的交界面上设置合适的速度和波面边界条件,建立了三维非线性波浪时域计算的耦合模型。模拟试验表明:(1)耦合模型数值波浪水池可以产生稳定的、重复性较好的波动过程;(2)用耦合模型数值波浪水池求解较大浅水区域上的非线性波浪数值计算问题可以取得较高的计算效率,同时又能得出结构物附近的复杂流场。 相似文献
20.
本文基于具备间断捕捉能力的二阶全非线性Boussinesq数值模型,对规则波和随机波在礁坪地形上的传播变形进行了数值模拟。该模型采用高阶有限体积法和有限差分方法求解守恒格式的控制方程,将波浪破碎视为间断,同时采用静态重构技术处理了海岸动边界问题。重点针对礁坪上波浪传播过程中的波高空间分布和沿程衰减,礁坪上的平均水位变化,以及波浪能量频谱的移动和空间差异等典型水动力现象开展数值计算。将数值结果与实验结果对比,两者吻合情况良好,验证了模型具有良好的稳定性,具备模拟破碎波浪和海-岸动边界的能力,能较为准确地模拟波浪在礁坪地形上的传播过程中发生的各种水动力现象。 相似文献