共查询到19条相似文献，搜索用时 0 毫秒
1.
建立了求解一维全非线性Green-Naghdi水波方程的中心有限体积/有限差分混合数值格式。采用结构化网格对守恒形式的控制方程进行离散和积分,界面数值通量采用有限体积法计算,剩余项则采用中心有限差分格式求解。其中,采用中心迎风有限体积格式计算控制体界面数值通量,并结合界面变量的线性重构方法,使其在空间上具有四阶精度,通过引入静压重构技术和波浪破碎指标使模型具备处理海岸水-陆动边界及波浪破碎的能力。时间积分则采用具有总时间变差减小(Total Variation Diminishing,TVD)性质的三阶龙格-库塔法进行。应用该模型对孤立波在常水深和斜坡海岸上的传播过程及规则波跨越潜堤传播的实验进行了数值模型研究,数值计算同解析解及实验数据吻合良好。 相似文献
2.
为提高复杂区域内风暴潮的预测能力,以上海黄浦江、长江口和沿岸近海为研究对象,构建了基于浅水方程的风暴潮模型。模型运用HLLC格式计算通量,具有高空间分辨率,能够模拟风暴潮、天文潮、洪水等多重因素作用下潮水位过程变化。它采用了非结构的四边形和三角形混合网格,适合应用于跨尺度河渠和河口等区域潮汐交互作用的分析模拟。为了验证模型,开展了0012号派比安台风和0509号麦莎台风影响下的潮汐变化过程计算。结果表明,该模型能够分析天文潮和台风组合作用下的潮水位变化,能够满足风暴潮预报要求。 相似文献
3.
海洋流体动力学的一种交替方向隐式二维数值模式 总被引：5，自引：2，他引：5
本文提出一种计算浅海潮汐、风暴潮及其他长周期波动的二维数值模式.动力学方程用建立在交错网格上的二步差分子方程近似,两步均为隐格式,其中一步用于计算水位和x方向速度分量,另一步计算水位和y方向速度分量.两步交替进行以获得每一时间步长的值,该格式具有二阶精度.文章用Von Neumann方法证明了该计算格式的线性差分方程是无条件稳定的.为检验数值模式的可靠性,本文作了一系列数值实验.在有一开边界的长方形平底和斜底水域中以开边界处周期性振动和海面处定常风应力驱动下,数值模拟所得结果与解析解完全一致;在水深不连续,边界封闭的矩形水域中,水位对定常风响应的数值解也与Leendertse模式计算结果完全一致.模式还用于对南海风暴潮进行后报,良好地再现了风暴增水过程. 相似文献
4.
5.
6.
流速逆变张量隐式求解方法及其在航道港池流场计算中的应用 总被引：6，自引：1，他引：5
运用高分辨率的边界适应网格进行流体动力学数值计算时,如何提高计算稳定性和减少计算量成为数值求解的关键性问题.在非正交的边界适应坐标系中,每个动量方程中同时出现了两个交叉方向的水位偏导数项,给隐式求解带来困难,而显式格式下的时间步长由于受与空间步长有关的Courant-Friedrichs-Lewy条件限制,计算量成倍增加.本文从广义曲线坐标系下浅海动力学方程组出发,导出了流速的逆变张量所满足的动量方程组,使方程中的水位偏导数项变成了沿某一协变基向量方向占优的形式,方便地采用了交替方向隐式差分格式,从而提高了计算稳定性并减小了计算量.本文通过对澳门海域航道和港池中流场的计算,证实了该模式是一种进行高分辩率数值计算的有效方法. 相似文献
7.
8.
海岸工程中实用三维流场数学模型研究 总被引：6，自引：0，他引：6
本文以完整的具有自由表面的三维非线性Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组为基础方程，经垂向σ坐标变换后得到模型方程。采用有限差分法进行数值求解；垂向上等间距分层，使用隐式中心差分格式；水平方向上给出两套网格－正方网格和不规则三角形网格离散形式。 相似文献
9.
贴体正交曲线坐标系下水流泥沙数学模型的构建与应用 总被引：1，自引：0，他引：1
为了研究深圳河河口福田保税区围垦填土占用部分河道对深圳河治理工程防洪能力和河床冲淤变化的影响情况，建立了贴体正交曲线坐标系下的平面二维的深圳河河道-河口-海湾整体水域的水流泥沙数学模型，模型中采用双步全隐有限差分法离散求解水流运动方程，悬沙不平衡输运方程采用三阶QUIKEST格式离散后用交错方向法求解。在经实测水文资料验证的基础上，应用数学模型计算了8种洪峰流量和河口高潮位组合4种围垦填土的挖除方案下深圳河沿程水位和河床冲淤变化，讨论了深圳河治理工程防洪能力降低的程度和近河口段洲滩冲刷以及所造成香港一侧红树林的破坏情况，提出了优化挖除方案。 相似文献
10.
潮汐河口断面悬沙通量组分模式及其在长江口的应用 总被引：2，自引：5，他引：2
断面泥沙通量估算的误差主要来源于计算方法、测点布局等,通量模式应该建立在通量估算统计误差最小的原则基础之上.在断面网格设计中采用统计误差最小的等面积单元网格,在泥沙通量估算中采用泥沙组分浓度,在流速变量插值上垂向采用对数函数插值,横向采用第一边界三次方样条函数插值.这样建立的潮汐河口悬沙断面通量组分模式较以往的任何模式更完善,断面通量估算的误差最小.将该模式应用于长江河口南港断面悬沙通量估算及其输移机制分析,断面泥沙通量表现为大潮期大进大出、大出大于大进;小潮期小进小出、小出大于小进;主要输移机制是拉格朗日输移和潮泵. 相似文献
11.
A two-dimensional vertical (2DV) non-hydrostatic boundary fitted model based on a Godunov-type shock-capturing scheme is introduced and applied to the simulation of waves from deep water up to the swash zone. The effects of shoaling, breaking, surf zone dissipation and swash motions are considered. The application of a Godunov-type shock-capturing algorithm together with an implicit solver on a standard staggered grid is proposed as a new approach in the 2DV simulation of large gradient problems such as wave breaking and hydraulic jumps. The complete form of conservative Reynolds averaged Navier–Stokes (RANS) equations are solved using an implicit finite volume method with a pressure correction technique. The horizontal advection of the horizontal velocity is solved by an explicit predictor–corrector method. Fluxes are predicted by an exact Riemann solver and corrected by a downwind scheme. A simple total variation diminishing (TVD) method with a monotonic upstream-centered scheme for conservation laws (MUSCL) limiter function is employed to eliminate undesirable oscillations across discontinuities. Validation of the model is carried out by comparing the results of the simulations with several experimental test cases of wave breaking and run-up and the analytical solution to linear short waves in deep water. Promising performance of the model has been observed. 相似文献
12.
Under conditions common in muddy coastal and estuarine environments, acoustic Doppler velocimeters (ADVs) can serve to estimate sediment settling velocity (w s) by assuming a balance between upward turbulent Reynolds flux and downward gravitational settling. Advantages of this method include simple instrument deployment, lack of flow disturbance, and relative insensitivity to biofouling and water column stratification. Although this method is being used with increasing frequency in coastal and estuarine environments, to date it has received little direct ground truthing. This study compared in situ estimates of w s inferred by a 5-MHz ADV to independent in situ observations from a high-definition video settling column over the course of a flood tide in the bottom boundary layer of the York River estuary, Virginia, USA. The ADV-based measurements were found to agree with those of the settling column when the current speed at about 40 cm above the bed was greater than about 20 cm/s. This corresponded to periods when the estimated magnitude of the settling term in the suspended sediment continuity equation was four or more times larger than the time rate of change of concentration. For ADV observations restricted to these conditions, ADV-based estimates of w s (mean 0.48±0.04 mm/s) were highly consistent with those observed by the settling column (mean 0.45±0.02 mm/s). However, the ADV-based method for estimating w s was sensitive to the prescribed concentration of the non-settling washload, C wash. In an objective operational definition, C wash can be set equal to the lowest suspended solids concentration observed around slack water. 相似文献
13.
ＡｎＡＤＩｔｅｃｈｎｉｑｕｅｆｏｒｓｏｌｖｉｎｇｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌｃｏａｓｔａｌｃｉｒｃｕｌａｔｉｏｎＹａｎｇＬｉａｎｗｕ１，＊ａｎｄＪｏｓｅＯｚｅｒ２（ＲｅｃｅｉｖｅｄＮｏｖｅｍｂｅｒ２７，１９９６；ａｃｃｅｐｔｅｄＤｅｃｅｍｂｅｒ... 相似文献
14.
We point out one problem of the grid advection schemes when used in wave models in coastal areas. The deficiency of the schemes is investigated by means of the ‘third' generation WAM wave model, in which the wave energy is advected by a first order upwind scheme. Two similar, alternative modifications of this scheme are analyzed, the second of which is shown to solve most of the problems encountered with advection along the co-ordinate axes. 相似文献
15.
四叉树网格下的椭圆型缓坡方程数值模型研究 总被引：1，自引：1，他引：0
波浪是近岸海域关键的水动力因素之一。考虑到近岸地形复杂、波浪演化显著的特点,建立了四叉树网格体系下的椭圆型缓坡方程数值模型,采用有限体积法对模型进行数值离散,应用GPBiCG(m, n)算法求解离散后的控制方程。模型中根据波浪波长布局计算网格,生成多层次四叉树网格,对复杂计算域有较好的适应性,并且在离散和方程求解中无需引入形函数、不产生复杂的交叉项,节约了存储空间和计算时间。将模型成功应用于物理模型实验及Acapulco海湾的波浪场数值模拟,结果表明该模型能够准确、高效地模拟近岸波浪场,可为近岸波浪场的模拟提供一定的理论和技术支持。 相似文献
16.
为建立高效的Boussinesq类水波数值模型,提出了一种新型的、基于有限差分和有限体积方法的混合数值格式。针对守恒形式的一维控制方程,在等间距矩形控制体内对其进行积分并离散,采用有限体积方法计算界面数值通量,剩余源项采用有限差分方法计算。其中,采用MUSTA格式并结合高精度状态插值方法计算控制体界面数值通量。时间积分则采用具有TVD性质的三阶龙格-库塔多步积分法进行。除验证模型外,重点对MUSTA格式和广泛使用的HLL格式进行了比较。结果表明,MUSTA格式可用于Boussinesq类水波方程数值求解,综合考虑数值精度、计算效率、程序编制和实际应用这几个方面,其较HLL格式更具有优势。 相似文献
17.
《Deep Sea Research Part II: Topical Studies in Oceanography》2009,56(24):2474-2486
The Coastal Gulf of Alaska (CGOA) is productive, with large populations of fish, seabirds, and marine mammals; yet it is subject to downwelling-favorable coastal winds. Downwelling regions in other parts of the world are typically much less productive than their upwelling counterparts. Alternate sources of nutrients to feed primary production in the topographically complex CGOA are poorly known and difficult to quantify. Here we diagnose the output from a spatially nested, coupled hydrodynamic and lower trophic level model of the CGOA, to quantify both horizontal and vertical nutrient fluxes into the euphotic zone. Our nested model includes both nitrogen and iron limitation of phytoplankton production, and is driven by a fine-scale atmospheric model that resolves the effects of local orography on the coastal winds. Results indicate significant “rivers” of cross-shelf nitrogen flux due to horizontal advection, as well as “fountains” of vertical transport over shallow banks due to tidal mixing. Using these results, we constructed a provisional budget of nutrient transport among subregions of the CGOA. Contrary to expectations, this budget reveals substantial upwelling of nutrients over major portions of the shelf, driven by local wind-stress curl. These effects are large enough to overwhelm the smaller downwelling flux at the coast throughout the growing season. Vertical mixing by winds and tides, and horizontal flux from the deep basin, are other substantial contributors to nutrients above the 15-m horizon. These findings help to explain the productivity of this coastal ecosystem. 相似文献
18.
《Coastal Engineering》2004,51(10):1021-1049
A numerical process-based model to forecast beach profile morphodynamics has been developed. In the present paper, an analysis of various modelling approaches and key parametrizations involved in the estimation of the wave-driven current and the suspended sediment concentration is carried out.Several resolution techniques for the 1DV horizontal (i.e., in the x-direction perpendicular to coastline) momentum equation governing the Mean Horizontal Velocity (MHV) are analysed. In the first kind of techniques, the mean horizontal velocity is computed from the momentum equation, whereas the Mean Water Level (MWL) is computed using a parametrization of the depth-averaged momentum equation. Two boundary or integral conditions are thus needed. In the second kind, both mean horizontal velocity and mean water level gradient in the x-direction are the unknowns of the momentum equation, thus, three boundary or integral conditions are needed. Various additional conditions are discussed. We show that using a technique of the first kind is equivalent to imposing the difference between the surface and the bottom shear stresses in the 1D vertical equation. Both techniques lead to results that are in good agreement with the Delta Flume experimental data, provided the Stokes drift flow discharge is imposed as an additional condition. The influence of the breaking roller model and of the turbulent viscosity parametrization are also analysed.Suspended sediment transport by the mean current and wave-induced bedload transport are taken into account in the sediment flux. Three turbulent diffusivity parametrizations are compared for suspended sediment concentration estimations. A linear profile for the turbulent diffusivity taking into account the wave bottom shear stress and the surface wave breaking turbulence production is shown to give the best results. Using experimental data, we put forward the poor estimation of the bottom sediment concentration given by the three implemented parametrizations. We thus propose a new parametrization relying on a Shields parameter based on the breaking roller induced surface shear stress. Using this new parametrization, the bottom profile used in the tests keeps its two bars which disappear otherwise. However, the morphodynamical model still overestimates the bars offshore motion, a bias already observed in other models. 相似文献
19.
A staggered finite-volume technique for non-hydrostatic, small amplitude free surface flow governed by the incompressible Navier-Stokes equations is presented there is a proper balance between accuracy and computing time. The advection and horizontal diffusion terms in the momentum equation are discretized by an integral interpolation method on the orthogonal unstructured staggered mesh and, while it has the attractive property of being conservative. The pressure-correction algorithm is employed for the non-hydrostatic pressure in order to achieve second-order temporal accuracy. A conservative scalar transport algorithm is also applied to discretize k-ε equations in this model. The eddy viscosity is calculated from the k-ε turbulent model. The resulting model is mass and momentum conservative. The model is verified by two examples to simulate unsteady small amplitude free surface flows where non-hydrostatic pressures have a considerable effect on the velocity field,and then applied to simulate the tidal flow in the Bohai Sea. 相似文献