共查询到20条相似文献，搜索用时 15 毫秒
1.
The present study develops a numerical model of the two-dimensional fully nonlinear shallow water equations （NSWE） for the wave run-up on a beach. The finite volume method （FVM） is used to solve the equations, and a second-order explicit scheme is developed to improve the computation efficiency. The numerical fluxes are obtained by the two dimensional Roe＇ s flux function to overcome the errors caused by the use of one dimensional fluxes in dimension splitting methods. The high-resolution Godunov-type TVD upwind scheme is employed and a second-order accuracy is achieved based on monotonic upstream schemes for conservation laws （MUSCL） variable extrapolation; a nonlinear limiter is applied to prevent unwanted spurious oscillation. A simple but efficient technique is adopted to deal with the moving shoreline boundary. The verification of the solution technique is carried out by comparing the model output with documented results and it shows that the solution technique is robust. 相似文献
2.
A hybrid finite-volume and finite-difference method is proposed for numerically solving the two-dimensional (2D) extended Boussinesq equations. The governing equations are written in such a way that the convective flux is approximated using finite volume (FV) method while the remaining terms are discretized using finite difference (FD) method. Multi-stage (MUSTA) scheme, instead of commonly used HLL or Roe schemes, is adopted to evaluate the convective flux as it has the simplicity of centred scheme and accuracy of upwind scheme. The third order Runge–Kutta method is used for time marching. Wave breaking and wet–dry interface are also treated in the model. In addition to model validation, the emphasis is given to compare the merits and limitations of using MUSTA scheme and HLL scheme in the model. The analytical and experimental data available in the literature have been used for the assessment. Numerical tests demonstrate that the developed model has the advantages of stability preserving, shock-capturing and numerical efficiency when applied in the complex nearshore region. Compared with that using HLL scheme, the proposed model has comparable numerical accuracy, but requires slightly less computation time and is much simpler to code. 相似文献
3.
为建立高效的Boussinesq类水波数值模型,提出了一种新型的、基于有限差分和有限体积方法的混合数值格式。针对守恒形式的一维控制方程,在等间距矩形控制体内对其进行积分并离散,采用有限体积方法计算界面数值通量,剩余源项采用有限差分方法计算。其中,采用MUSTA格式并结合高精度状态插值方法计算控制体界面数值通量。时间积分则采用具有TVD性质的三阶龙格-库塔多步积分法进行。除验证模型外,重点对MUSTA格式和广泛使用的HLL格式进行了比较。结果表明,MUSTA格式可用于Boussinesq类水波方程数值求解,综合考虑数值精度、计算效率、程序编制和实际应用这几个方面,其较HLL格式更具有优势。 相似文献
4.
In this paper, a hybrid scheme based on a set of 2DH extended Boussinesq equations for slowly varying bathymetries is introduced. The numerical code combines the finite volume technique, applied to solve the advective part of the equations, with the finite difference method, used to discretize dispersive and source terms. Time integration is performed using the fourth-order Adams–Bashforth–Moulton predictor–corrector method; the Riemann problem is solved employing an approximate HLL solver, a fourth-order MUSCL-TVD technique is applied. Five test cases, for non-breaking and breaking waves, are reproduced to verify the model comparing its results to laboratory data or analytical solutions. 相似文献
5.
A depth-integrated, non-hydrostatic model with hybrid finite difference and finite volume numerical algorithm is proposed in this paper. By utilizing a fraction step method, the governing equations are decomposed into hydrostatic and non-hydrostatic parts. The first part is solved by using the finite volume conservative discretization method, whilst the latter is considered by solving discretized Poisson-type equations with the finite difference method. The second-order accuracy, both in time and space, of the finite volume scheme is achieved by using an explicit predictor-correction step and linear construction of variable state in cells. The fluxes across the cell faces are computed in a Godunov-based manner by using MUSTA scheme. Slope and flux limiting technique is used to equip the algorithm with total variation dimensioning property for shock capturing purpose. Wave breaking is treated as a shock by switching off the non-hydrostatic pressure in the steep wave front locally. The model deals with moving wet/dry front in a simple way. Numerical experiments are conducted to verify the proposed model. 相似文献
6.
Vegetation damping effects on propagating water waves have been investigated by many researchers. This paper investigates the effects of damping due to vegetation on solitary water wave run-up via numerical simulation. The numerical model is based on an implementation of Morison's formulation for vegetation induced inertia and drag stresses in the nonlinear shallow water equations. The numerical model is solved via a finite volume method on a Cartesian cut cell mesh. The accuracy of the numerical scheme and the effects of the vegetation terms in the present model are validated by comparison with experiment results. The model is then applied to simulate a solitary wave propagating on a plane slope with vegetation. The sensitivity of solitary wave run-up to plant height, diameter and stem density is investigated by comparison of the numerical results for different patterns of vegetation. The numerical results show that vegetation can effectively reduce solitary wave propagation velocity and that solitary wave run-up is decreased with increase of plant height in water and also diameter and stem density. 相似文献
7.
一种基于Roe格式的有限体积法在二维溃坝问题中的应用 总被引：3，自引：0，他引：3
溃坝属于典型的非恒定含间断的浅水问题。应用有限体积法离散二维浅水控制方程的守恒型方程组,将基于近似黎曼解的Roe格式用于数值计算溃坝问题,并利用MUSCL方法构造二阶空间积分格式和预测-校正二步法构造二阶时间格式,从而使数值解的整体达到二阶,提高了精度。文中算法在一维溃坝的Stoker问题的数值结果与解析解进行对比,结果证明了此方法的可行性。应用此方法在二维溃坝问题上的结果,说明了此算法可有效模拟溃坝水流的演进过程。 相似文献
8.
数值模拟作为海啸预报的主要研究方法在海啸预警中起着关键作用。本文采用Godunov格式的有限体积方法,使用MUSCL-Hancock格式,并利用HLLC Riemann近似求解器计算单元界面上的流体通量,建立了球坐标系下二阶精度的海啸数值模型。模型所基于的全和谐型浅水方程保证了数值的稳定性,而地形重构方法实现了干湿边界的精准模拟。本文模拟了2015年9月16日智利Mw8.3级地震海啸,通过与智利近岸14个测站和环太平洋20个DART浮标实测数据比较,验证了模型对实际越洋海啸模拟预报的能力。 相似文献
9.
Impact problems associated with water entry have important applications in various aspects of naval architecture and ocean engineering. Estimation of hydrodynamic impact forces especially during the first instances after the impact is very important and is of interest. Since the estimation of hydrodynamic impact load plays an important role in safe design and also in evaluation of structural weight and costs, it is better to use a reliable and accurate prediction method instead of a simple estimation resulted by analyzing methods. In landing of flying boats, some phenomena such as weather conditions and strong winds can cause asymmetric instead of symmetric descent. In this paper, a numerical simulation of the asymmetric impact of a wedge, as the step of a flying boat, considering dynamic equations in two-phase flow is taken into account. The dynamic motion of the wedge in two-phase flow is solved based on finite volume method with volume of fluid (VOF) scheme considering dynamic equations. Then the effects of different angles of impact and water depth on the velocity change and slamming forces in an asymmetric impact are investigated. The comparison between the simulation results and experimental data verifies the accuracy of the method applied in the present study. 相似文献
10.
《Coastal Engineering》2004,51(1):1-15
The shallow water equations (SWE) have been used to model a series of experiments examining violent wave overtopping of a near-vertical sloping structure with impacting wave conditions. A finite volume scheme was used to solve the shallow water equations. A monotonic reconstruction method was applied to eliminate spurious oscillations and ensure proper treatment of bed slope terms. Both the numerical results and physical observations of the water surface closely followed the relevant Rayleigh probability distributions. However, the numerical model overestimated the wave heights and suffered from the lack of dispersion within the shallow water equations. Comparisons made on dimensionless parameters for the overtopping discharge and percentage of waves overtopping between the numerical model and the experimental observations indicated that for the lesser impacting waves, the shallow water equations perform satisfactorily and provide a good alternative to computationally more expensive methods. 相似文献
11.
The incompressible viscous uniform and shear flow past a circular cylinder is studied. The two-dimensional Navier-Stokes equations are solved by a finite element method. The governing equations are discretized by a weighted residual method in space. The stable three-step scheme is applied to the momentum equations in the time integration. The numerical model is firstly applied to the computation of the lid-driven cavity flow for its validation. The computed results agree well with the measured data and other numerical results. Then, it is used to simulate the viscous uniform and shear flow over a circular cylinder for Reynolds numbers from lO0 to lO00. The transient time interval before the vortex shedding occurs is shortened considerably by introduction of artificial perturbation. The computed Strouhal number, drag and lift coefficients agree well with the experimental data. The computation shows that the finite element model can be successfully applied to the viscous flow problem. 相似文献
12.
Finite volume scheme for the solution of 2D extended Boussinesq equations in the surf zone 总被引：2，自引：0，他引：2
In this paper, a hybrid finite volume-finite difference scheme is applied to study surf zone dynamics. The numerical model solves the 2DH extended Boussinesq equations proposed by Madsen and Sørensen (1992) where nonlinear and dispersive effects are both relevant whereas it solves NSWE equations where nonlinearity prevails. The shock-capturing features of the finite volume method allow an intrinsic representation of wave breaking and runup; therefore no empirical (calibration) parameters are necessary. Comparison with laboratory measurements demonstrates that the proposed model can accurately predict wave height decay and mean water level setup, for both regular and solitary wave breaking on a sloping beach. The model is also applied to reproduce two-dimensional wave transformation and breaking over a submerged circular shoal, showing good agreement with experimental data. 相似文献
13.
1 .IntroductionThe development and application of hydrodynamic models to predict circulation in estuarine ,coastal ,and shelf waters has been extremely active over the past three decades (Spauldinget al .,1992 ,1994) . However , most of these models are b… 相似文献
14.
An improved coupling of numerical and physical models for simulating 2D wave propagation is developed in this paper. In the proposed model, an unstructured finite element model (FEM) based Boussinesq equations is applied for the numerical wave simulation, and a 2D piston-type wavemaker is used for the physical wave generation. An innovative scheme combining fourth-order Lagrange interpolation and Runge-Kutta scheme is described for solving the coupling equation. A Transfer function modulation method is presented to minimize the errors induced from the hydrodynamic invalidity of the coupling model and/or the mechanical capability of the wavemaker in area where nonlinearities or dispersion predominate. The overall performance and applicability of the coupling model has been experimentally validated by accounting for both regular and irregular waves and varying bathymetry. Experimental results show that the proposed numerical scheme and transfer function modulation method are efficient for the data transfer from the numerical model to the physical model up to a deterministic level. 相似文献
15.
An improved coupling of numerical and physical models for simulating 2D wave propagation is developed in this paper. In the proposed model, an unstructured finite element model （FEM） based Boussinesq equations is applied for the numerical wave simulation, and a 2D piston-type wavemaker is used for the physical wave generation. An innovative scheme combining fourth-order Lagrange interpolation and Runge-Kutta scheme is described for solving the coupling equation. A Transfer function modulation method is presented to minimize the errors induced from the hydrodynamic invalidity of the coupling model and/or the mechanical capability of the wavemaker in area where nonlinearities or dispersion predominate. The overall performance and applicability of the coupling model has been experimentally validated by accounting for both regular and irregular waves and varying bathymetry. Experimental results show that the proposed numerical scheme and transfer function modulation method are efficient for the data transfer from the numerical model to the physical model up to a deterministic level. 相似文献
16.
17.
实际工程中存在大量的曲边界,因此在曲边界上的计算准确性可以考察出一个数值模型的实用价值。利用Beji的改进型Boussinesq方程建立了一个有限元方法的数值波浪模型。造波方面采用Fenton提出的非线性规则波浪解;在墙边界处,以求解法线方向和切线方向的速度和导数代替求解x、y方向的速度和导数,从而使边界条件直接适用、严格满足,保证了对曲边界计算的准确性。"重开始广义极小残量法"的使用保证了求解方程组的效率和精度,使造波和边界处理方法的有效性和准确性得到了合理地诠释。通过与试验数据、他人数值结果、解析解的比对,显示出该模型计算稳定、结果准确,真正体现出了有限元方法对曲边界适用的优势。 相似文献
18.
19.
《Ocean Modelling》2010,34(3-4):235-256
In this article we propose and implement a multilevel method to simulate the solution of the two-dimensional nonlinear shallow-water equations. The multilevel method is based on a central-upwind finite volume scheme and uses new incremental unknowns which enable to preserve the numerical conservation of the scheme. The method is tested and analyzed on two and three levels of discretization on different test cases and turns out to furnish a good solution of the problems while saving CPU time. 相似文献
20.
四叉树网格下的椭圆型缓坡方程数值模型研究 总被引：1，自引：1，他引：0
波浪是近岸海域关键的水动力因素之一。考虑到近岸地形复杂、波浪演化显著的特点,建立了四叉树网格体系下的椭圆型缓坡方程数值模型,采用有限体积法对模型进行数值离散,应用GPBiCG(m, n)算法求解离散后的控制方程。模型中根据波浪波长布局计算网格,生成多层次四叉树网格,对复杂计算域有较好的适应性,并且在离散和方程求解中无需引入形函数、不产生复杂的交叉项,节约了存储空间和计算时间。将模型成功应用于物理模型实验及Acapulco海湾的波浪场数值模拟,结果表明该模型能够准确、高效地模拟近岸波浪场,可为近岸波浪场的模拟提供一定的理论和技术支持。 相似文献