首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 93 毫秒
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.  相似文献   

2.
Wave-tide-surge coupled simulation for typhoon Maemi   总被引:1,自引:0,他引:1  
The main task of this study focuses on studying the effect of wave-current interaction on currents, storm surge and wind wave as well as effects of current induced wave refraction and current on waves by using numerical models which consider the bottom boundary layer and sea surface roughness parameter for shallow and smooth bed area around Korean Peninsula. The coupled system (unstructured-mesh SWAN wave and ADCIRC) run on the same unstructured mesh. This identical and homogeneous mesh allows the physics of wave-circulation interactions to be correctly resolved in both models. The unstructured mesh can be applied to a large domain allowing all energy from deep to shallow waters to be seamlessly followed. There is no nesting or overlapping of structured wave meshes, and no interpolation is required. In response to typhoon Maemi (2003), all model components were validated independently, and shown to provide a faithful representation of the system’s response to this storm. The waves and storm surge were allowed to develop on the continental shelf and interact with the complex nearshore environment. The resulting modeling system can be used extensively for prediction of the typhoon surge. The result show that it is important to incorporate the wave-current interaction effect into coastal area in the wave-tide-surge coupled model. At the same time, it should consider effects of depth-induced wave breaking, wind field, currents and sea surface elevation in prediction of waves. Specially, we found that: (1) wave radiation stress enhanced the current and surge elevation otherwise wave enhanced nonlinear bottom boundary layer decreased that, (2) wind wave was significantly controlled by sea surface roughness thus we cautiously took the experimental expression. The resulting modeling system can be used for hindcasting (prediction) the wave-tide-surge coupled environments at complex coastline, shallow water and fine sediment area like areas around Korean Peninsula.  相似文献   

3.
This paper proposes an enhanced approach for evaluating the fatigue life of each metallic layer of unbonded flexible risers. Owing to the complex structure of unbonded flexible risers and the nonlinearity of the system, particularly in the critical touchdown zone, the traditional method is insufficient for accurately evaluating the fatigue life of these risers. The main challenge lies in the transposition from global to local analyses, which is a key stage for the fatigue analysis of flexible pipes owing to their complex structure. The new enhanced approach derives a multi-layer stress-decomposition method to meet this challenge. In this study, a numerical model validated experimentally is used to demonstrate the accuracy of the stress-decomposition method. And a numerical case is studied to validate the proposed approach. The results demonstrate that the multi-layer stress-decomposition method is accurate, and the fatigue lives of the metallic layers predicted by the enhanced multi-layer analysis approach are rational. The proposed fatigue-analysis approach provides a practical and reasonable method for predicting fatigue life in the design of unbonded flexible risers.  相似文献   

4.
This paper adopted a semi-analytical method based on eigenfunction matching to solve the problem of sharp resonance of cylindrical structures with a moonpool that has a restricted entrance. To eliminate the sharp resonance and to measure the viscous effect, a quadratic dissipation is introduced by assuming an additional dissipative disk at the moonpool entrance. The fluid domain is divided into five cylindrical subdomains, and the velocity potential in each subdomain is obtained by meeting the Laplace equation as well as the boundary conditions. The free-surface elevation at the center of the moonpool, along with the pressure and velocity at the restricted entrance for first-order wave are evaluated. By choosing appropriate dissipation coefficients, the free-surface elevation calculated at the center of the moonpool is in coincidence with the measurements in model tests both at the peak period and amplitude at resonance. It is shown that the sharp resonance in the potential flow theory can be eliminated and the viscous effect can be estimated with a simple method in some provided hydrodynamic models.  相似文献   

5.
吕彪 《中国海洋工程》2013,27(4):509-522
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.  相似文献   

6.
The characteristics of three-dimensional (3-D) tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, many islands, moveable boundary, and submerged dyke, therefore, a 3-D numerical model on an unstructured triangular grid has been developed. The σ coordinate transformation, the moveable boundary and submerged dyke treatment techniques were employed in the model so it is suitable for the tidal simu...  相似文献   

7.
—The Monte-Carlo method is used to simulate the surface fatigue crack growth rate foroffshore structural steel E36-Z35,and to determine the distributions and relevance of the parameters in theParis equation.By this method,the time and cost of fatigue crack propagation testing can be reduced.Theapplication of the method is demonstrated by use of four sets of fatigue crack propagation data foroffshore structural steel E36-Z35.A comparison of the test data with the theoretical prediction for surfacecrack growth rate shows the application of the simulation method to the fatigue crack propagation tests issuccessful.  相似文献   

8.
A numerical model for wave propagation in a harbour is verified by use of physical models.The extended time-dependent mild slope equation is employed as the governing equation,and the model is solved by use of ADI method containing the relaxation factor.Firstly,the reflection coefficient of waves in front of rubble-mound breakwaters under oblique incident waves is determined through physical model tests,and it is regarded as the basis for simulating partial reflection boundaries of the numerical model.Then model tests on refraction,diffraction and reflection of waves in a harbour are performed to measure wave height distribution.Comparative results between physical and numerical model tests show that the present numerical model can satisfactorily simulate the propagation of regular and irregular waves in a harbour with complex topography and boundary conditions.  相似文献   

9.
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.  相似文献   

10.
Smoothed Particle Hydrodynamics method (SPH) has a good adaptability for simulating of free surface flow problems. However, there are some shortcomings of SPH which are still in open discussion. This paper presents a corrected solid boundary handling method for weakly compressible SPH. This improved method is very helpful for numerical stability and pressure distribution. Compared with other solid boundary handling methods, this corrected method is simpler for virtual ghost particle interpolation and the ghost particle evaluation relationship is clearer. Several numerical tests are given, like dam breaking, solitary wave impact and sloshing tank waves. The results show that the corrected solid boundary processing method can recover the spurious oscillations of pressure distribution when simulating the problems with complex geometry boundary.  相似文献   

11.
适于模拟不规则水域波浪的缓坡方程两种数值模型比较   总被引:1,自引:1,他引:0  
本文分析比较了适于不规则水域波浪模拟的椭圆型缓坡方程两种数值模型。两种数值模型均采用有限体积法离散,分别基于四叉树网格和非结构化三角形网格建立。首先结合近岸缓坡地形上波浪传播的经典物理模型实验对两种数值模型分别进行了验证,并结合计算结果对比分析了两种模型的计算精度和效率。计算结果表明,两种数值模型均可有效地模拟近岸波浪的传播变形;相对非结构化三角形网格下的模型,基于四叉树网格建立的数值模型在数值离散和求解过程中无需引入形函数、不产生复杂的交叉项,离散简单,易于程序实现,且节约计算存储空间,计算效率高。  相似文献   

12.
文章基于Delaunay三角化的思想实现了非结构化数值网格生成技术,并LOP(Local Optimization Procedure)算法对网格进行了优化。用有限体积法在非结构三角形网格上对平面二维潮流方程进行了数值求解。实例应用分析表明,文中的数值网格生成技术和有限体积方法相结合用于求解具有复杂几何边界区域内的潮流问题能取得很好的结果。  相似文献   

13.
基于2011年兴化湾及周边海域的流速、流向的观测资料,应用短期资料的潮流准调和分析方法,对海流周日连续观测数据进行分析,探讨海域潮流特征。最后,针对兴化湾岸线复杂、岛屿众多的特点,基于FVCOM数值模式,采用非结构三角形网格,建立兴化湾及周边海域潮流动力数值模型,经实测站潮汐、潮流数据对模拟结果验证,模型能够较好的刻画各测站潮汐、潮流过程特点,反映研究海域流场特征规律,可为该海域海洋环境保护和海洋开发利用提供科学论证依据。  相似文献   

14.
强潮河口三维无结构网格盐度数学模型   总被引:3,自引:0,他引:3  
采用平面无结构网格和垂向直接分层,建立了强潮河口盐度三维斜压力数学模型。无结构网格使模型能够适应河口复杂的边界,并可以根据需要进行局部加密。Casulli的半隐式法离散模型水位梯度和垂项紊动扩散项,克服了水位梯度和垂项紊动扩散对模型稳定性的影响,Semi-Lagrangian法用于模型对流项的离散,使模型具有"无条件"稳定性。采用水池风生流和盐度异重流算例检验了模型,并应用模型模拟了强潮河口(瓯江口)的潮流盐度运动,计算和实测结果进行比较,表明该模型是准确、可靠的。  相似文献   

15.
A two-dimensional scour model based on coupled system of shallow water equations (SWEs) and sediment transport on unstructured mesh is developed. The coupled system of hydrodynamic and morphodynamic equations is solved by finite volume method using Godunov scheme. Roe's approximate Riemann solver is used to calculate the inviscid fluxes. The use of unstructured mesh makes the model applicable to complex domains. However, it is difficult to evaluate the eigenvalues and eigenvectors of the Jacobian matrix in the global coordinate. The method proposed herein to deal with this difficulty is to transform the system into the local coordinate with one of the axes in the same direction as the interface outward normal vector. In the local coordinate system, the Jacobian matrix is simplified and the eigenvalues are analyzed using asymptotic method. Regular expansion breaks down when the flow is near critical. Uniformity of the expansion is achieved by changing the scales. Rotational invariance theorem is used to relate the interfacial fluxes in the global and local coordinate systems. Special treatment of the source term on unstructured grid makes the scheme stable and physically balanced (both mass and momentum). The method proposed in this paper for the eigen-system is very efficient comparing to iterative numerical methods. Results from the test cases show good agreement with the experiments.  相似文献   

16.
为完全拟合河口近海复杂岸线和工程结构以及有效局部加密,设计并建立了一个无结构三角形网格二维河口海岸水动力数值模式。空间离散主要基于有限体积法以保证守恒性,时间积分采用预估修正法以提高精度。水位在三角形网格中心通过连续方程求解;水平x方向和y方向的流速U和V均在网格边中点上通过动量方程求解。流速平流项的求解中采用了TVD格式。TVD流速平流通量为一个一阶迎风格式通量和一个二阶格式通量的组合,一阶格式通量和二阶格式通量根据流速的局部分布情况得出配比,最终组合得到TVD通量。TVD格式具有低耗散和无频散的优点,提高了模式的稳定性。应用实测资料验证建立的模式,结果显示水位、流速和流向的计算值与实测值均符合良好。  相似文献   

17.
A method for generating high quality unstructured computational meshes for applications in geophysical, and specifically ocean, modelling is presented. The method generates a three-dimensional mesh which is optimised to approximate the bathymetry to a user-specified accuracy, which itself may be defined as a function of space and bathymetry. To do this a mesh optimisation method is employed which allows the formation of highly anisotropic elements aligned with bathymetric variations of high curvature. In addition a coastal recovery algorithm, based upon feature decimation and edge recovery, is described.  相似文献   

18.
《Ocean Modelling》2011,40(3-4):248-261
Accurate representation of geostrophic and hydrostatic balance is an essential requirement for numerical modelling of geophysical flows. Potentially, unstructured mesh numerical methods offer significant benefits over conventional structured meshes, including the ability to conform to arbitrary bounding topography in a natural manner and the ability to apply dynamic mesh adaptivity. However, there is a need to develop robust schemes with accurate representation of physical balance on arbitrary unstructured meshes. We discuss the origin of physical balance errors in a finite element discretisation of the Navier–Stokes equations using the fractional timestep pressure projection method. By considering the Helmholtz decomposition of forcing terms in the momentum equation, it is shown that the components of the buoyancy and Coriolis accelerations that project onto the non-divergent velocity tendency are the small residuals between two terms of comparable magnitude. Hence there is a potential for significant injection of imbalance by a numerical method that does not compute these residuals accurately. This observation is used to motivate a balanced pressure decomposition method whereby an additional “balanced pressure” field, associated with buoyancy and Coriolis accelerations, is solved for at increased accuracy and used to precondition the solution for the dynamical pressure. The utility of this approach is quantified in a fully non-linear system in exact geostrophic balance. The approach is further tested via quantitative comparison of unstructured mesh simulations of the thermally driven rotating annulus against laboratory data. Using a piecewise linear discretisation for velocity and pressure (a stabilised P1P1 discretisation), it is demonstrated that the balanced pressure decomposition method is required for a physically realistic representation of the system.  相似文献   

19.
In this work, a combined immersed boundary (IB) and volume of fluid (VOF) methodology is developed to simulate the interactions of free-surface waves and submerged solid bodies. The IB method is used to account for the no-slip boundary condition at solid interfaces and the VOF method, utilizing a piecewise linear interface calculation, is employed to track free surfaces. The combined model is applied in several case studies, including the propagation of small-amplitude progressive waves over a submerged trapezoidal dike, a solitary wave traveling over a submerged rectangular object, and wave generation induced by a moving bed. Numerical results depicting the free-surface evolutions and velocity fields are in good agreement with either experimental data or numerical results obtained by other researchers. In addition, the simplification of the initial free-surface deformation used in most tsunami earthquake source study is justified by the present model application. The methodology presented in the paper serves as a good tool for solving many practical problems involving free surfaces and complex boundaries.  相似文献   

20.
设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号