共查询到20条相似文献，搜索用时 15 毫秒
1.
In this paper, a numerical model is established for simulating the wave forces on a submarine pipeline. A set of two-dimensional Navier–Stokes equations is discretized numerically with a finite volume method in a moving mesh system. After each time step, the mesh is modified according to the changed wave surface boundary. The deffered correction second-order upwind scheme (SUDC) is adopted here to discretize the convective fluxes. The effects of the clearance between the pipeline and the seabed, water depth and wave height on wave forces are studied, respectively. The results by the numerical simulation agree well with the experimental data and theory value. 相似文献
2.
数值模拟作为海啸预报的主要研究方法在海啸预警中起着关键作用。本文采用Godunov格式的有限体积方法,使用MUSCL-Hancock格式,并利用HLLC Riemann近似求解器计算单元界面上的流体通量,建立了球坐标系下二阶精度的海啸数值模型。模型所基于的全和谐型浅水方程保证了数值的稳定性,而地形重构方法实现了干湿边界的精准模拟。本文模拟了2015年9月16日智利Mw8.3级地震海啸,通过与智利近岸14个测站和环太平洋20个DART浮标实测数据比较,验证了模型对实际越洋海啸模拟预报的能力。 相似文献
3.
《Ocean Modelling》2007,16(1-2):47-60
The symmetry and stability properties of two unstructured C-grid discretisations of the shallow water equations are discussed. We establish that a scheme in which the circumcentres of the mesh triangles are used as the surface elevation points has advantageous symmetry properties and derive a Coriolis discretisation which preserves these properties. It is shown that the resulting scheme is conservative in a discretised energy norm. We then establish that schemes in which the water surface elevations are stored at the mesh triangle centroids do not share these advantageous symmetry properties. Finally we show examples which demonstrate that the centroid based scheme is subject to unstable growing modes, particularly in long timescale, Coriolis dominated problems; while the energy conservative circumcentre based scheme suffers from no such limitation. We conclude that unstructured C-grid methods using the triangle circumcentres and the conservative Coriolis scheme derived here therefore have advantages for this sort of problem over those schemes based on centroids. 相似文献
4.
无结构网格二维河口海岸水动力数值模式的建立及其应用 总被引：2，自引：0，他引：2
为完全拟合河口近海复杂岸线和工程结构以及有效局部加密,设计并建立了一个无结构三角形网格二维河口海岸水动力数值模式。空间离散主要基于有限体积法以保证守恒性,时间积分采用预估修正法以提高精度。水位在三角形网格中心通过连续方程求解;水平x方向和y方向的流速U和V均在网格边中点上通过动量方程求解。流速平流项的求解中采用了TVD格式。TVD流速平流通量为一个一阶迎风格式通量和一个二阶格式通量的组合,一阶格式通量和二阶格式通量根据流速的局部分布情况得出配比,最终组合得到TVD通量。TVD格式具有低耗散和无频散的优点,提高了模式的稳定性。应用实测资料验证建立的模式,结果显示水位、流速和流向的计算值与实测值均符合良好。 相似文献
5.
In this study, 1D and 2D shallow-water models were coupled to simulate unsteady flow in channel networks and embayment. The 1D model solved the 1D shallow-water equations (St. Venant) using the Preissmann box method and targeted long narrow reaches of the river networks, while the 2D model targeted broad channels and embayment and solved the 2D shallow-water equations using a semi-implicit scheme applied to an unstructured grid of triangular cells. The 1D and 2D models were solved simultaneously by building a matrix for the free surface elevation at every 1D junction and 2D cell center. Velocities were then computed explicitly based on the results at the previous time step and the updated water level. The originality of the scheme arose from a novel coupling method. The results showed that the coupled 1D/2D model produced identical results as the full 2D model in classical to benchmark problems with considerable savings in computational effort. Application of the model to the Pearl River Estuary in southern China showed that complex patterns of tidal wave propagation could be efficiently modeled. 相似文献
6.
7.
This paper addresses the issue of the non conservation of tracer that occurs when a leapfrog time stepping scheme is used in association with a Robert-Asselin time filter within the framework of a time varying vertical coordinate system. A re-formulation of both the forcing and filtering terms that provides a globally conservative scheme is presented. Local conservation is then achieved when the model levels and subsequently the sea surface height satisfy a local compatibility condition. The properties of the resulting scheme are illustrated by a global coupled ocean–sea ice model based on the Nucleus for European Modelling of the Ocean (NEMO). 相似文献
8.
Free vibrations of three-dimensional extensible marine cables with specified top tension via a variational method 总被引：1，自引：0，他引：1
This paper presents a model formulation that can be used for analyzing the three-dimensional vibration behaviours of an inclined extensible marine cable. The virtual work-energy functional, which involves strain energy due to axial stretching of the cable and virtual work done by external hydrostatic forces is formulated. The coupled equations of motion in the Cartesian coordinates of global systems are obtained by taking into account the difference between Euler’s equations and equilibrium equations. The method of Galerkin finite element is used to obtain the mass and stiffness matrices which are transformed into the local coordinate systems. Then the eigenvalue problem is solved to determine its natural frequencies and corresponding mode shapes. The model formulation developed herein is conveniently applied for the cases of specified top tension. The numerical investigations are carried out to demonstrate the validity of the model and to explore in details the influence of various parameters on the behaviours of marine cables. Results for the frequency avoidance phenomenon, maximum dynamic tension and coupled transverse mode shapes are presented and discussed. 相似文献
9.
A time-domain simulation method based on potential flow model has been developed to investigate the berthing problem between two floating bodies in wave. The boundary value problem is formulated with respect to an earth-fixed coordinate system because the relative positions of the two vessels continuously change during the berthing operation. The classical finite element method is used to solve the Laplace equation in the fluid domain with moving boundary. The linearized free-surface boundary conditions are integrated in time by applying 4th-order Adams–Bashforth–Moulton method. A simple re-mesh algorithm with local and global mesh systems is introduced to update mesh by considering large horizontal movement of the berthing vessel. The developed numerical method is used to investigate the berthing problem between a FPSO and shuttle tanker in waves. The focus is on the wave-induced motion response during the berthing process. The characteristics of the motion responses in berthing operation are examined with various wave frequencies, berthing speeds and wave headings. 相似文献
10.
In this paper we propose an integral form of the fully non-linear Boussinesq equations in contravariant formulation, in which Christoffel symbols are avoided, in order to simulate wave transformation phenomena, wave breaking and nearshore currents in computational domains representing the complex morphology of real coastal regions. Following the approach proposed by Chen (2006), the motion equations retain the term related to the approximation to the second order of the vertical vorticity. A new Upwind Weighted Essentially Non-Oscillatory scheme for the solution of the fully non-linear Boussinesq equations on generalised curvilinear coordinate systems is proposed. The equations are rearranged in order to solve them by a high resolution hybrid finite volume–finite difference scheme. The conservative part of the above-mentioned equations, consisting of the convective terms and the terms related to the free surface elevation, is discretised by a high-order shock-capturing finite volume scheme in which an exact Riemann solver is involved; dispersive terms and the term related to the approximation to the second order of the vertical vorticity are discretised by a cell-centred finite difference scheme. The shock-capturing method makes it possible to intrinsically model the wave breaking, therefore no additional terms are needed to take into account the breaking related energy dissipation in the surf zone. The model is verified against several benchmark tests, and the results are compared with experimental, theoretical and alternative numerical solutions. 相似文献
11.
An unstructured-grid procedure for SWAN is presented. It is a vertex-based, fully implicit, finite difference method which can accommodate unstructured meshes with a high variability in geographic resolution suitable for representing complicated bottom topography in shallow areas and irregular shoreline. The numerical solution is found by means of a point-to-point multi-directional Gauss–Seidel iteration method requiring a number of sweeps through the grid. The approach is stable for any time step while permitting local mesh refinements in areas of interest. A number of applications are shown to verify the correctness and numerical accuracy of the unstructured version of SWAN. 相似文献
12.
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. 相似文献
13.
《Coastal Engineering》2005,52(1):43-62
A vertical two-dimensional (2D) numerical model for time dependent local scour below offshore pipelines subject to unidirectional steady flow is developed. The governing equations for the flow and sediment transport are solved by using finite difference method in a general curvilinear coordinate system. The performance of two turbulence models, the standard k–ɛ model and Smagorinsky subgrid scale (SGS) model, on modeling time dependent scour processes is examined. Both suspended load and bed load are considered in the scour model. The suspended-load model is verified against two channel sediment transport cases. The change of bed level is calculated from the continuity equation of total sediment transport. A new time marching scheme and a sand slide scheme are proposed for the scour calculation. It is found that the proposed time marching scheme and sand slide model work well for both clear-water and live-bed scour situations and the standard k–ɛ turbulence closure is more preferable than the SGS model in the 2D scour model developed in this study. 相似文献
14.
In this paper, the generalised two-dimensional differential transform method (DTM) of solving the time-fractional coupled KdV equations is proposed. The fractional derivative is described in the Caputo sense. The presented method is a numerical method based on the generalised Taylor series expansion which constructs an analytical solution in the form of a polynomial. An illustrative example shows that the generalised two-dimensional DTM is effective for the coupled equations. 相似文献
15.
Parametric vertical coordinate formulation for multiscale, Boussinesq, and non-Boussinesq ocean modeling 总被引：2，自引：1，他引：2
Two physical parameters are introduced into the basic ocean equations to generalize numerical ocean models for various vertical coordinate systems and their hybrid features. The two parameters are formulated by combining three techniques: the arbitrary vertical coordinate system of Kasahara [Kasahara, A., 1974. Various vertical coordinate systems used for numerical weather prediction. Mon. Weather Rev. 102, 509–522], the Jacobian pressure gradient formulation of Song [Song, Y.T., 1998. A general pressure gradient formation for ocean models. Part I: Scheme design and diagnostic analysis. Mon. Weather Rev. 126 (12), 3213–3230], and a newly introduced parametric function that permits both Boussinesq (volume-conserving) and non-Boussinesq (mass-conserving) conditions. Based on this new formulation, a generalized modeling approach is proposed. Several representative oceanographic problems with different scales and characteristics––coastal canyon, seamount topography, non-Boussinesq Pacific Ocean with nested eastern Tropics, and a global ocean model––have been used to demonstrate the model’s capabilities for multiscale applications. The inclusion of non-Boussinesq physics in the topography-following ocean model does not incur computational expense, but more faithfully represents satellite-observed ocean-bottom-pressure data. Such a generalized modeling approach is expected to benefit oceanographers in solving multiscale ocean-related problems by using various coordinate systems on the same numerical platform. 相似文献
16.
A full time-domain analysis program is developed for the coupled dynamic analysis of offshore structures. For the hydrodynamic loads, a time domain second order method is developed. In this approach, Taylor series expansions are applied to the body surface and free-surface boundary conditions, and the Stokes perturbation procedure is then used to establish the corresponding boundary value problems with time-independent boundaries. A higher-order boundary element method (HOBEM) is developed to calculate the velocity potential of the resulting flow field at each time step. The free-surface boundary condition is satisfied to the second order by fourth order Adams–Bashforth–Moultn method. An artificial damping layer is adopted on the free surface to avoid the wave reflection. The mooring-line/tendon/riser dynamics are based on the rod theory and the finite element method (FEM), with the governing equations described in a global coordinate system. In the coupled dynamic analysis, the motion equation for the hull and dynamic equations for mooring-lines/tendons/risers are solved simultaneously using the Newmark method. The coupled analysis program is applied for a truss Spar motion response simulation. Numerical results including motions and tensions at the top of mooring-lines/risers are presented, and some significant conclusions are derived. 相似文献
17.
《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. 相似文献
18.
A novel computational approach is introduced for the efficient computation of equilibrium solutions of seasonally forced ocean biogeochemical models. The essential idea is to formulate the problem as a large system of nonlinear algebraic equations to be solved with a class of methods known as matrix-free Newton–Krylov (MFNK). MFNK is a combination of Newton-type methods for superlinearly convergent solution of nonlinear equations, and Krylov subspace methods for solving the Newton correction equations. The basic link between the two methods is the Jacobian-vector product, which may be probed approximately without forming and storing the elements of the true Jacobian. To render this approach practical for global models with O(106) degrees of freedom, a flexible preconditioning strategy is developed. The result is an essentially “black-box” numerical scheme than can be applied to most existing biogeochemical models. The method is illustrated by applying it to find the equilibrium solutions of two realistic biogeochemical problems. Compared with the conventional approach of direct time integration, the preconditioned-MFNK scheme is shown to be roughly two orders of magnitude more efficient. Several potential refinements of the basic algorithm that may yield further performance gains are discussed. The numerical scheme described here addresses a fundamental challenge to using ocean biogeochemical models more effectively. 相似文献
19.
The three-dimensional (3-D) eddy-current transient field problem is formulated first using theu-v method. This method breaks the vector Helmholtz equation into two scalar Helmholtz equations. Null-field integral equations and the appropriate boundary conditions germane to the problem are used to set up an identification matrix which is independent of null-field point locations. Embedded in the identification matrix are the unknown eigenvalues of the problem representing its impulse response in time. These eigenvalues are found by equating the determinant of the identification matrix to zero. The eigenvalues, which can be equated with temporal response, are found to be intimately linked to the initial forcing function which triggers the transient in question. When this initial forcing function is Fourier decomposed into its respective spatial harmonics, it is possible to associate with each Fourier component a unique eigenvalue by this technique. The true transient solution comes through a convolution of the impulse response so obtained with the particular imposed external field governing the problem at hand. The technique is applied to the FELIX medium cylinder (a conducting cylinder placed in a collapsing external field) and compared to data. A pseudoanalytic confirmation of the eigenvalues so obtained is formulated to validate the procedure. The technique proposed is applied in the low-frequency regime where the near-field effects must be considered. Application of the technique to a high frequency follows directly if the Coulomb gauge is adopted to represent the vector potential. 相似文献
20.
Numerical modeling of nonlinear resonance of semi-enclosed water bodies: Description and experimental validation 总被引：1，自引：0，他引：1
Inigo J. Losada Jose M. Gonzalez-Ondina Gabriel Diaz-HernandezErnesto M. Gonzalez 《Coastal Engineering》2008
The purpose of this work is to a present a numerical model to solve a set of modified Boussinesq equations to analyse nonlinear resonance of semi-enclosed water bodies. The equations are solved on a finite element unstructured grid in order to achieve an optimal mesh resolution with the local geometry. The model is able to simulate long time lapses and realistic forcing in real bathymetries with a reasonable computational cost. To validate the numerical results, a set of experiments was carried out in a physical model of two adjacent elongated basins. Comparisons between numerical and experimental results for different geometries and nonlinear conditions show that the model is able to simulate with an excellent agreement the transient nonlinear resonant process. 相似文献