共查询到20条相似文献，搜索用时 21 毫秒
1.
A Vertical Two-Dimensional Model to Simulate Tidal Hydrodynamics in A Branched Estuary 总被引：1，自引：1，他引：1
A vertical (laterally averaged) two-dimensional hydrodynamic model is developed for tides, tidal current, and salinity in a branched estuarine system. The goveming equations are solved with the hydrostatic pressure distribution assumption and the Boussinesq approximation. An explicit scheme is employed to solve the continuity equations. The momentum and mass balance equations are solved implicitly in the Cartesian coordinate system. The tributaries are govemed by the same dynamic equations. A control volume at the junctions is designed to conserve mass and volume transport in the finite difference schemes, based on the physical principle of continuum medium of fluid. Predictions by the developed model are compared with the analytic solutions of steady wind-driven circulatory flow and tidal flow. The model results for the velocities and water surface elevations coincide with analytic results. The model is then applied to the Tanshui River estuarine system. Detailed model calibration and verification have been conducted with measured water surface elevations,tidal current, and salinity distributions. The overall performance of the model is in qualitative agreement with the available field data. The calibrated and verified numerical model has been used to quantify the tidal prism and flushing rate in the Tanshui River-Tahan Stream, Hsintien Stream, and Keelung River. 相似文献
2.
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. 相似文献
3.
The surface heat budget of the Arctic Ocean （SHEBA） project has shown that the study of the surface heat budget characteristics is crucial to understanding the interface process and environmental change in the polar region. An arctic single - column model （ARCSCM） of Colorado University is used to simulate the arctic surface radiation and energy budget during the summertime. The simulation results are analyzed and compared with the SHEBA measurements. Sensitivity analyses are performed to test microphysical and radiative parameterizations in this model. The results show that the ARCSCM model is able to simulate the surface radiation and energy budget in the arctic during the summertime, and the different parameterizations have a significant influence on the results. The combination of cloud microphysics and RRTM parameterizations can fairly derive the surface solar shortwave radiation and downwelling longwave radiation flux. But this cloud microphysics parameterization scheme deviates notably from the simulation of surface sensible and latent heat flux. Further improvement for the parameterization scheme applied to the Arctic Regions is necessary. 相似文献
4.
Wave Numerical Model for Shallow Water 总被引：4，自引：0，他引：4
The history of forecasting wind waves by wave energy conservation equation is briefly des-cribed.Several currently used wave numerical models for shallow water based on different wave theoriesare discussed.Wave energy conservation models for the simulation of shallow water waves are introduced,with emphasis placed on the SWAN model,which takes use of the most advanced wave research achieve-ments and has been applied to several theoretical and field conditions.The characteristics and applicabilityof the model,the finite difference numerical scheme of the action balance equation and its source termscomputing methods are described in detail.The model has been verified with the propagation refractionnumerical experiments for waves propagating in following and opposing currents;finally.the model is ap-plied to the Haian Gulf area to simulate the wave height and wave period field there,and the results arecompared with observed data. 相似文献
5.
Researches on breaking-induced currents by waves are summarized firstly in this paper. Then, a combined numerical model in orthogonal curvilinear coordinates is presented to simulate wave-induced current in areas with curved boundary or irregular coastline. The proposed wave-induced current model includes a nearshore current module established through orthogonal curvilinear transformation form of shallow water equations and a wave module based on the curvilinear parabolic approximation wave equation. The wave module actually serves as the driving force to provide the current module with required radiation stresses. The Crank-Nicolson finite difference scheme and the alternating directions implicit method are used to solve the wave and current module, respectively. The established surf zone currents model is validated by two numerical experiments about longshore currents and rip currents in basins with rip channel and breakwater. The numerical results are compared with the measured data and published numerical results. 相似文献
6.
Parabolic Approximation of the Weakly Nonlinear Mild Slope Equation with Bottom Friction 总被引：2，自引：0，他引：2
This paper presents a refined parabolic approximation model of the mild slope equation to simu-late the combination of water wave refraction and diffraction in the large coastal region.The bottom frictionand weakly nonlinear term are included in the model.The difference equation is established with the Crank-Nicolson scheme.The numerical test shows that some numerical prediction results will be inaccurate in com-plicated topography without considering weak nonlinearity;the bottom friction will make wave height damp-ing and it can not be neglected for calculation of wave field in large areas. 相似文献
7.
This study deals with the general numerical model to simulate the two-dimensional tidal flow, flooding wave (long wave) and shallow water waves (short wave). The foundational model is based on nonlinear Boussinesq equations. Numerical method for modelling the short waves is investigated in detail. The forces, such as Coriolis forces, wind stress, atmosphere and bottom friction, are considered. A two-dimensional implicit difference scheme of Boussinesq equations is proposed. The low-reflection outflow open boundary is suggested. By means of this model,both velocity fields of circulation current in a channel with step expansion and the wave diffraction behind a semi-infinite breakwater are computed, and the results are satisfactory. 相似文献
8.
To examine the circulation in shallow water with tidal flat, a finite element model for the numerical solution of the shallow water equations was developed by means of standard Galerkin's method. The domain computed was covered with triangular finite elements, and water elevation and velocity were approximated by linear interpolation functions, and the lumped coefficients were used to substitute for solving the high order algebraic equation system. The time-dependent land-water boundary changes are treated mathematically by interrelating the location of the land-water boundary with the instantaneous tidal level. The implicit scheme was adopted for the terms of the bottom friction and the Coriolis effect in the motion equation so that the numerical stability of the model has been improved.The model was applied to the tidal current on shoaling water with large tidal flat off Pikou, and a comparison between observed and calculated values showed good agreement, the flow pattern being reproduced. The result 相似文献
9.
Sediment transport in the Hangzhou Bay is extremely complicated due to its bathymetry and hydrodynamic conditions. The ECOMSED model is employed to simulate three-dimensional (3-D) cohesive sediment transport in Hangzhou Bay. Dynamical factors such as Coriolis force, tides, salinity, river discharges, and waves are considered in the model. The wave parameters, including the significant wave height, period, and direction, are calculated with the SWAN model. The Grant-Madsen model is introduced for the bed shear stress due to the combined effect of waves and currents. The formulation of bed shear stress used to calculate the sink/source terms is modified based on previous research that sufficiently validated the formulation with measurement data. The integrated model of the above-mentioned models is applied to simulate sediment transport in Hangzhou Bay. The results of the simulation agree well with field observations concerning the distribution of suspended sediment, indicating that the sediments are remarkably suspended in Hangzhou Bay under the action of waves and currents. 相似文献
10.
一个两时间层分裂显格式海洋环流模式(MASNUM)及其检验 总被引：1，自引：0，他引：1
HAN Lei 《海洋学报(英文版)》2014,33(11):11-35
A two-time-level, three-dimensional numerical ocean circulation model（named MASNUM） was established with a two-level, single-step Eulerian forward-backward time-differencing scheme. A mathematical model of large-scale oceanic motions was based on the terrain-following coordinated, Boussinesq, Reynolds-averaged primitive equations of ocean dynamics. A simple but very practical Eulerian forward-backward method was adopted to replace the most preferred leapfrog scheme as the time-differencing method for both barotropic and baroclinic modes. The forward-backward method is of second-order of accuracy, computationally efficient by requiring only one function evaluation per time step, and free of the computational mode inherent in the three-level schemes. This method is superior to the leapfrog scheme in that the maximum time step of stability is twice as large as that of the leapfrog scheme in staggered meshes thus the computational efficiency could be doubled. A spatial smoothing method was introduced to control the nonlinear instability in the numerical integration. An ideal numerical experiment simulating the propagation of the equatorial Rossby soliton was performed to test the amplitude and phase error of this new model. The performance of this circulation model was further verified with a regional（northwest Pacific） and a quasi-global（global ocean simulation with the Arctic Ocean excluded） simulation experiments. These two numerical experiments show fairly good agreement with the observations. The maximum time step of stability in these two experiments were also investigated and compared between this model and that model which adopts the leapfrog scheme. 相似文献
11.
The purpose of this paper is to extend the validity of Li's parabolic model (1994) by incorporating a combined energy factor in the mild-slope equation and by improving the traditional radiation boundary conditions. With wave breaking and energy dissipation expressed in a direct form in the equation, the proposed model could provide an efficient numerical scheme and accurate predictions of wave transformation across the surf zone. The radiation boundary conditions are iterated in the model without use of approximations. The numerical predictions for wave height distributions across the surf zone are compared with experimental data over typical beach profiles. In addition, tests of waves scattering around a circular pile show that the proposed model could also provide reasonable improvement on the radiation boundary conditions for large incident angles of waves. 相似文献
12.
Numerical Modelling of Standing Waves with Three-Dimensional Non-Linear Wave Propagation Model 总被引：1，自引：0，他引：1
ZHANG Hongsheng HONG Guangwen DING Pingxing 《中国海洋工程》2001,(4):521-530
Based on the theoretical high-order model with a dissipative term for non-linear and dispersive wave in water of varying depth, a 3-D mathematical model of non-linear wave propagation is presented. The model, which can be used to calculate the wave particle velocity and wave pressure, is suitable to the complicated topography whose relative depth ratio of the characteristic water depth to the characteristic wavelength in deep-water) is equal to or smaller than one. The governing equations are discretized with the improved 2-D Crank-Nicolson method in which the first-order derivatives are corrected by Taylor series expansion, .and the general boundary conditions with an arbitrary reflection coefficient and phase shift are adopted in the model. The surface elevation, horizontal and vertical velocity components and wave pressure of standing waves are numerically calculated. The results show that the numerical model can effectively simulate the complicated standing waves, and the general boundary conditions 相似文献
13.
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. 相似文献
14.
Two-Dimensional Mathematical Model of Tidal Current and Sediment for Oujiang Estnary and Wenzhou Bay 总被引：1，自引：0，他引：1
A 2-D mathematical model of tidal current and sediment has been developed for the Oujiang Estuary and the WenzhouBay. This model accomodates complicated features including multiple islands, existence of turbidity, and significant differ-ence in size distribution of bed material. The governing equations for non-uniform suspended load and bed load transport arepresented in a boundary-fitted orthogonal curvilinear coordinate system. The numerical solution procedures along with theirinitial conditions, boundary conditions, and movable boundary technique are presented. Strategies for computation of thecritical condition of deposition or erosion, sediment transport capacity, non-uniform bed load discharge, etc. are suggested.The model verification computation shows that, the tidal levels computed from the model are in good agreement with the fielddata at the 18 tidal gauge stations. The computed velocities and flow directions also agree well with the values measuredalong the totally 52 synchronously observed verticals distributed over 8 cross sections. The computed tidal water throughputsthrough the Huangda‘ao cross section are close to the measured data. And the computed values of bed deformation fromYangfushan to the estuary outfall and in the outer-sea area are in good agreement with the data observed from 1986 to 1992.The changes of tidal volumes through the estuary, velocities in different channels and the bed form due to the influence of thereclamation project on the Wenzhou shoal are predicted by means of this model. 相似文献
15.
Diapycnal mixing is important in oceanic circulation. An inverse method in which a semi-explicit scheme is applied to discretize the one-dimensional temperature diffusion equation is established to estimate the vertical temperature diffusion coefficient based on the observed temperature profiles. The sensitivity of the inverse model in the idealized and actual conditions is tested in detail. It can be found that this inverse model has high feasibility under multiple situations ensuring the stability of the inverse model, and can be considered as an efficient way to estimate the temperature diffusion coefficient in the weak current regions of the ocean. Here, the hydrographic profiles from Argo floats are used to estimate the temporal and spatial distribution of the vertical mixing in the north central Pacific based on this inverse method. It is further found that the vertical mixing in the upper ocean displays a distinct seasonal variation with the amplitude decreasing with depth, and the vertical mixing over rough topography is stronger than that over smooth topography It is suggested that the high-resolution profiles from Argo floats and a more reasonable design of the inverse scheme will serve to understand mixing processes. 相似文献
16.
In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split into diffusive and convective parts in each time step. The diffusive part is discretized by the backward difference method in time and discretized by the standard Galerkin method in space. The convective part is a first-order nonlinear equation.After the linearization of the nonlinear part by Newton's method, the convective part is also discretized by the backward difference method in time and discretized by least square scheme in space. C~0-type element can be used for interpolation of the velocity and pressure in the present model. Driven cavity flow and flow past a circular cylinder are conducted to validate the present model. Numerical results agree with previous numerical results, and the model has high accuracy and can be used to simulate problems with complex geometry. 相似文献
17.
The performance of a z-level ocean model, the Modular Ocean Model Version 4(MOM4), is evaluated in terms of simulating the global tide with different horizontal resolutions commonly used by climate models. The performance using various sets of model topography is evaluated. The results show that the optimum filter radius can improve the simulated co-tidal phase and that better topography quality can lead to smaller rootmean square(RMS) error in simulated tides. Sensitivity experiments are conducted to test the impact of spatial resolutions. It is shown that the model results are sensitive to horizontal resolutions. The calculated absolute mean errors of the co-tidal phase show that simulations with horizontal resolutions of 0.5° and 0.25° have about 35.5% higher performance compared that with 1° model resolution. An internal tide drag parameterization is adopted to reduce large system errors in the tidal amplitude. The RMS error of the best tuned 0.25° model compared with the satellite-altimetry-constrained model TPXO7.2 is 8.5 cm for M_2. The tidal energy fluxes of M_2 and K_1 are calculated and their patterns are in good agreement with those from the TPXO7.2. The correlation coefficients of the tidal energy fluxes can be used as an important index to evaluate a model skill. 相似文献
18.
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. 相似文献
19.
This paper presents a simulation model based on the finite element method. The method is used to analyze the motion response and mooring line tension of the flatfish cage system in waves. The cage system consists of top frames, netting, mooring lines, bottom frames, and floats. A series of scaled physical model tests in regular waves are conducted to verify the numerical model. The comparison results show that the simulated and the experimental results agree well under the wave conditions, and the maximum pitch of the bottom frame with two orientations is about 12o. The motion process of the whole cage system in the wave can be described with the computer visualized technology. Then, the mooring line tensions and the motion of the bottom frame with three kinds of weight are calculated under different wave conditions. According to the numerical results, the differences in mooring line tensions of flatfish cages with three weight modes are indistinct. The maximum pitch of the bottom frame decreases with the increase of the bottom weight. 相似文献
20.
2-D Current Field Numerical Simulation Integrating Yangtze Estuary with Hangzhou Bay 总被引：9，自引：2，他引：7
—In this paper,integrating the Yangtze Estuary with the Hangzhou Bay,a 2-D velocity fieldmodel is established.In the model,fine self-adaptive grids are employed to adapt to the complicated coast-al shape.The hydrodynamic equations satisfied by two contravariant components of velocity vector andsurface elevation in non-orthogonal curvilinear coordinates are used.In each momentum equation thecoefficients before the two partial derivatives of surface elevation with respect to variables of alternative di-rection coordinates have different orders of magnitude,i.e.,the derivative with the larger coefficient mayplay a more important role than that with the smaller one.With this advantage,the ADI scheme can beeasily employed.The hydrodynamic factors include tidal current,river runoff and wind-induced current.In terms of tidal current,seven main constituents in the area are considered in the open boundaries.Theverifications of surface elevation process and current velocity process in the spring tide and in the n 相似文献