共查询到19条相似文献，搜索用时 125 毫秒
1.
建立了求解一维全非线性Green-Naghdi水波方程的中心有限体积/有限差分混合数值格式。采用结构化网格对守恒形式的控制方程进行离散和积分,界面数值通量采用有限体积法计算,剩余项则采用中心有限差分格式求解。其中,采用中心迎风有限体积格式计算控制体界面数值通量,并结合界面变量的线性重构方法,使其在空间上具有四阶精度,通过引入静压重构技术和波浪破碎指标使模型具备处理海岸水-陆动边界及波浪破碎的能力。时间积分则采用具有总时间变差减小(Total Variation Diminishing,TVD)性质的三阶龙格-库塔法进行。应用该模型对孤立波在常水深和斜坡海岸上的传播过程及规则波跨越潜堤传播的实验进行了数值模型研究,数值计算同解析解及实验数据吻合良好。 相似文献
2.
3.
4.
5.
6.
二维非线性浅水波的数值模拟 总被引：3，自引：0，他引：3
研究工作的目的在于建立一个能够模拟二维潮流、洪水波(长波)和浅水波浪(短波)的综合数学模型.基本模型建立在非线性的Boussinesq方程基础之上.本文主要讨论浅水波浪即短波的数值模拟.模型可以考虑必要的外力项,如柯氏力、风应力、大气压力和底摩阻力等.针对Boussinesq方程提出了一个全隐的二维差分格式,讨论了人工开边界的处理方法.模型被用来计算了突然扩张渠槽中的环流和单突堤后的水波绕射,取得了满意的结果. 相似文献
7.
Liu等给出的最高导数为2的双层Boussinesq水波方程具有较好的色散性和非线性,基于该方程建立了有限差分法的三维波浪数值模型。在矩形网格上对方程进行了空间离散,采用高阶导数近似方程中的时、空项,时间积分采用混合4阶Adams-Bashforth-Moulton的预报—校正格式。模拟了深水条件下的规则波传播过程,计算波面与解析结果吻合较好,反映出数值模型能很好地刻画波面过程及波面处的速度变化;在kh=2π条件下可较为准确获得沿水深分布的水平和垂向速度,这与理论分析结果一致。最后,利用数值模型计算了规则波在三维特征地形上的传播变形,数值结果和试验数据吻合较好;高阶非线性项会对波浪数值结果产生一定的影响,当波浪非线性增强,水深减少将产生更多的高次谐波。建立的双层Boussinesq模型对强非线性波浪的演化具有较好的模拟精度。 相似文献
8.
建立基于四阶完全非线性Boussinesq水波方程的二维波浪传播数值模型。采用Kennedy等提出的涡粘方法模拟波浪破碎。在矩形网格上对控制方程进行离散,采用高精度的数值格式对离散方程进行数值求解。对规则波在具有三维特征地形上的传播过程进行了数值模拟,通过数值模拟结果与实验结果的对比,对所建立的波浪传播模型进行了验证。同时,为了考察非线性对波浪传播的影响,给出和上述模型具有同阶色散性、变浅作用性能但仅具有二阶完全非线性特征的波浪模型的数值结果。通过对比两个模型的数值结果以及实验数据,讨论非线性在波浪传播过程中的作用。研究结果表明,所建立的Boussinesq水波方程在深水范围内不但具有较精确的色散性和变浅作用性能,而且具有四阶完全非线性特征,适合模拟波浪在近岸水域的非线性运动。 相似文献
9.
应用内外解匹配的方法和驻相法推导了柱坐标系下地震引起的水面波动方程的解,即推导了地震海啸波生成与传播的理论方程,并对越洋地震海啸理论方程进行了求解。基于Boussinesq方程出发,建立了二维越洋海啸波传播数学模型,并对越洋海啸进行了数值计算,计算方法采用有限差分方法,差分格式采用交替方向隐格式(Alternating Direction Implicit即ADI方法)。利用越洋海啸计算模式对发生在大洋中的地震海啸进行了模拟,将数值模拟结果与地震海啸波理论方程的计算结果进行了比较,两种计算结果吻合较好。 相似文献
10.
本文基于具备间断捕捉能力的二阶全非线性Boussinesq数值模型,对规则波和随机波在礁坪地形上的传播变形进行了数值模拟。该模型采用高阶有限体积法和有限差分方法求解守恒格式的控制方程,将波浪破碎视为间断,同时采用静态重构技术处理了海岸动边界问题。重点针对礁坪上波浪传播过程中的波高空间分布和沿程衰减,礁坪上的平均水位变化,以及波浪能量频谱的移动和空间差异等典型水动力现象开展数值计算。将数值结果与实验结果对比,两者吻合情况良好,验证了模型具有良好的稳定性,具备模拟破碎波浪和海-岸动边界的能力,能较为准确地模拟波浪在礁坪地形上的传播过程中发生的各种水动力现象。 相似文献
11.
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. 相似文献
12.
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. 相似文献
13.
This study aims to develop a robust, accurate and computationally efficient hydrodynamic and sediment transport model for dam break flows. The two dimensional shallow water equations are resolved based on the finite volume method with an unstructured quadtree mesh. The sediment transport and bed evolution modules are coupled with hydrodynamic module to predict simultaneously the hydrodynamics, sediment concentrations and morphological changes. The interface flux is computed by the HLL approximate Riemann solver with second order accuracy. The effects of pressure and gravity are included in source term in this model, which can simplify the computation and eliminate numerical imbalance between source and flux terms. For dam break flows occurring in complicated geometries, the quadtree rectangular mesh is used to refine the interesting area and important part. The model is first verified against results from laboratory experiments, existing numerical models and real life case. It is then used to simulate dam break flows over a mobile bed to investigate the bed evolution. The results are compared with experimental data and field data with good agreement. The method is simple, efficient, and conservative. It shows promise for handling hydrodynamic simulation and sediment transport for a wide range of dam break flows. 相似文献
14.
为完全拟合河口近海复杂岸线和工程结构以及有效局部加密,设计并建立了一个无结构三角形网格二维河口海岸水动力数值模式。空间离散主要基于有限体积法以保证守恒性,时间积分采用预估修正法以提高精度。水位在三角形网格中心通过连续方程求解;水平x方向和y方向的流速U和V均在网格边中点上通过动量方程求解。流速平流项的求解中采用了TVD格式。TVD流速平流通量为一个一阶迎风格式通量和一个二阶格式通量的组合,一阶格式通量和二阶格式通量根据流速的局部分布情况得出配比,最终组合得到TVD通量。TVD格式具有低耗散和无频散的优点,提高了模式的稳定性。应用实测资料验证建立的模式,结果显示水位、流速和流向的计算值与实测值均符合良好。 相似文献
15.
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. 相似文献
16.
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. 相似文献
17.
Presented here is a compact explicit difference scheme of high accuracy for solving the extended Boussinesq equations.For time discretization,a three-stage explicit Runge-Kutta method with TVD property is used at predicting stage,a cubic spline function is adopted at correcting stage,which made the time discretization accuracy up to fourth order;For spatial discretization,a three-point explicit compact difference scheme with arbitrary order accuracy is employed.The extended Boussinesq equations derived by Beji and Nadaoka are solved by the proposed scheme.The numerical results agree well with the experimental data.At the same time,the comparisons of the two numerical results between the present scheme and low accuracy difference method are made,which further show the necessity of using high accuracy scheme to solve the extended Boussinesq equations.As a valid sample,the wave propagation on the rectangular step is formulated by the present scheme,the modelled results are in better agreement with the experimental data than those of Kittitanasuan. 相似文献
18.
A Finite Volume Method with Unstructured Triangular Grids for Numerical Modeling of Tidal Current 总被引：1，自引：1，他引：1
The finite volume method （FVM） has many advantages in 2-D shallow water numerical simulation. In this study, the finite volume method is used with unstructured triangular grids to simulate the tidal currents. The Roe scheme is applied in the calculation of the intercell numerical flux, and the MUSCL method is introduced to improve its accuracy. The time integral is a two-step scheme of forecast and revision. For the verification of the present method, the Stoker＇s problem is calculated and the result is compared with the mathematically analytic solutions. The comparison indicates that the method is feasible. A sea area of a port is used as an example to test the method established here. The result shows that the present computational method is satisfactory, and it could be applied to the engineering fields. 相似文献
19.
Efficient Numerical Solution of the Modified Mild-Slope Equation 总被引：12，自引：1，他引：11
An efficient numerical model for wave refraction,diffraction and reflection is presented in thispaper.In the model,the modified time-dependent mild-slope equation is transformed into an evolutionequation and an improved ADI method involving a relaxation factor is adopted to solve it.The methodhas the advantage of improving the numerical stability and convergence rate by properly determining therelaxation factor.The range of the relaxation factor making the differential scheme unconditionally stableis determined by stability analysis.Several verifications are performed to examine the accuracy of the pres-ent model.The numerical results coincide with the analytic solutions or experimental data very well,andthe computer time is reduced. 相似文献