共查询到20条相似文献，搜索用时 78 毫秒
1.
A three-dimensional non-hydrostatic numerical model for simulation of the free-surface stratified flows is presented. The model is a non-hydrostatic extension of free-surface primitive equation model with a general vertical coordinate and horizontal orthogonal curvilinear coordinates. The model equations are integrated with mode-splitting technique and decomposition of pressure and velocity fields on hydrostatic and non-hydrostatic components. The model was tested against laboratory experiments on the steep wave transformation over the longshore bar, solitary wave impact on the vertical wall, the collapse of the mixed region in the thin pycnocline, mixing in the lock-exchange flows and water exchange through the sea strait. The agreement is generally fair.Responsible Editor: Hans Burchard 相似文献
2.
A 3D non-hydrostatic model is developed to compute internal waves. A novel grid arrangement is incorporated in the model. This not only ensures the homogenous Dirichlet boundary condition for the non-hydrostatic pressure can be precisely and easily imposed but also renders the model relatively simple in its discretized form. The Perot scheme is employed to discretize horizontal advection terms in the horizontal momentum equations, which is based on staggered grids and has the conservative property. Based on previous water wave models, the main works of the present paper are to (1) utilize a semi-implicit, fractional step algorithm to solve the Navier-Stokes equations (NSE); (2) develop a second-order flux-limiter method satisfying the max–min property; (3) incorporate a density equation, which is solved by a high-resolution finite volume method ensuring mass conservation and max–min property based on a vertical boundary-fitted coordinate system; and (4) validate the developed model by using four tests including two internal seiche waves, lock-exchange flow, and internal solitary wave breaking. Comparisons of numerical results with analytical solutions or experimental data or other model results show reasonably good agreement, demonstrating the model’s capability to resolve internal waves relating to complex non-hydrostatic phenomena. 相似文献
3.
A wetting and drying method for free-surface problems for the three-dimensional, non-hydrostatic Navier–Stokes equations is proposed. The key idea is to use a horizontally fixed mesh and to apply different boundary conditions on the free-surface in wet and dry zones. In wet areas a combined pressure/free-surface kinematic boundary condition is applied, while in dry areas a positive water level and a no-normal flow boundary condition are enforced. In addition, vertical mesh movement is performed to accurately represent the free-surface motion. Non-physical flow in the remaining thin layer in dry areas is naturally prevented if a Manning–Strickler bottom drag is used. The treatment of the wetting and drying processes applied through the boundary condition yields great flexibility to the discretisation used. Specifically, a fully unstructured mesh with any finite element choice and implicit time discretisation method can be applied. The resulting method is mass conservative, stable and accurate. It is implemented within Fluidity-ICOM [1] and verified against several idealized test cases and a laboratory experiment of the Okushiri tsunami. 相似文献
4.
A numerical experiment is carried out to investigate the suitability of a Boussinesq-type momentum model for simulating transcritical flows at short length transitions in open channel flow measuring structures. Two one-dimensional Boussinesq-type equation models, which incorporate different degrees of dynamic pressure corrections, are considered for this purpose. A finite difference method is employed to discretise and solve the equations. The models are then applied to simulate different test cases for flows in such channels with predominant non-hydrostatic pressure distribution effects. A comparison of the computed results with the corresponding experimental data is presented. Results of this study reveal that the proposed model, which includes a higher-order correction for the effect of the centrifugal pressure, describes well even relatively abrupt changes from sub- to super-critical flow state. 相似文献
5.
Three-dimensional numerical modelling of wind-driven circulation in a homogeneous lake 总被引：1，自引：0，他引：1
A three-dimensional numerical model has been developed to study wind-induced circulation patterns in a shallow homogeneous lake with a complex bathymetry. The governing equations are the unsteady Reynolds-Averaged Navier–Stokes equations in which the non-hydrostatic pressure distribution has been included. The model was tested against analytical solutions and laboratory data for wind-induced currents and then applied to Esthwaite Water, a small lake in Cumbria, UK. The model was used to study the main model parameters and to generate typical circulation patterns for a variety of conditions in the lake. Simulations showed that a non-hydrostatic pressure distribution did not have any noticeable influence on the overall circulation pattern in the lake. However, comparisons with field data at some measurement stations in the near-shore region with sharply varying bottom topography showed that the hydrodynamic pressure component had some influence on the vertical velocity profile. 相似文献
6.
在空间域偶层位法的基础上,研究了完整的频率域偶层位曲面位场处理和转换方法.该法可应用于平面或曲面、规则网或非规则网的位场数据处理和转换.通过对偶层面z坐标和计算面z坐标平移不同的量来加速正演快速收敛和保证反演稳定、快速收敛;提出了适合于不规则网曲面处理和转换的核心算法——单点快速Fourier变换;提出了频率域不规则网曲面处理和转换方法技术.通过以上技术措施解决了大数据量特别是曲面不规则网的位场处理和转换问题,模型试算以及实际资料处理验证了该方法的应用效果. 相似文献
7.
A robust well-balanced finite volume model for shallow water flows with wetting and drying over irregular terrain 总被引：2，自引：0，他引：2
An unstructured Godunov-type finite volume model is developed for the numerical simulation of geometrically challenging two-dimensional shallow water flows with wetting and drying over convoluted topography. In the framework of sloping bottom model, a modified formulation of shallow water equations is used to preserve mass conservation during flooding and recession. The key ingredient of the model is the use of this combination of the sloping bottom model and the modified shallow water equations to provide a robust technique for wet/dry fronts tracking and, together with centered discretization of the bed slope source term, to exactly preserve the static flow on irregular topographies. The variable reconstruction technique ensures nonnegative reconstructed water depth and reasonable reconstructed velocity, and the friction terms are solved by semi-implicit scheme that does not invert the direction of velocity components. The robustness and accuracy of the proposed model are assessed by comparing numerical and reference results of extensive test cases. Moreover, the results of a dam-break flooding over real topography are presented to show the capability of the model on field-scale application. 相似文献
8.
9.
Secondary flows induced by the blocking effect of a river plume on a transverse upwelling are investigated in a microtidal region of freshwater influence (ROFI). A nested version of the SYMPHONIE primitive-equation free-surface model for 3-D baroclinic coastal flows has been developed for the Rhône ROFI. The main characteristics of the model are a generalized sigma coordinate system in finite differences, using a time splitting for external and internal modes and high-order numerical advection schemes for density fields in combination with an modified turbulence closure scheme. The nesting system consists of two grids forced by the high-resolution ALADIN model atmospheric data. The coarse grid of 3 km resolution for the whole Gulf of Lions allows the forcing of the Liguro-Provençal large-scale current when the fine mesh of 1-km resolution is centred on the river mouth of the Grand Rhône. Documented field experiments from the Biodypar 3 field campaign performed during March 1999 are used for validation. Numerical results, CTD profiles and a SPOT TSM visible image are in good agreement concerning the shape and structure of the river plume. Other coastal flow features can be observed from satellite imagery. Computations of realistic situations recover these main secondary structures. Complementary process-oriented runs give an explanation of how the coastal upwelling induced by an inhomogeneous offshore wind is destabilized by the combination of the river plume and along-shelf current-blocking effects. In the end, a factor-separation analysis provides evidence that the locally non-linear effects in momentum contribute to the occurrence of secondary vortices.Responsible Editor: Phil Dyke 相似文献
10.
Margarita Markina Alexander Gavrikov Sergey Gulev Bernard Barnier 《Ocean Dynamics》2018,68(11):1593-1604
The spatial resolution of wind forcing fields is critical for modeling ocean surface waves. We analyze here the performance of the non-hydrostatic numerical weather prediction system WRF-ARW (Weather Research and Forecasting) run with a 14-km resolution for hindcasting wind waves in the North Atlantic. The regional atmospheric model was run in the domain from 20° N to 70° N in the North Atlantic and was forced with ERA-Interim reanalysis as initial and boundary conditions in a spectral nudging mode. Here, we present the analysis of the impact of spectral nudging formulation (cutoff wavelengths and depth through which full weighting from reanalysis data is applied) onto the performance of the modeled 10-m wind speed and wind wave fields for 1 year (2010). For modeling waves, we use the third-generation spectral wave model WAVEWATCH III. The sensitivity of the atmospheric and wave models to the spectral nudging formulation is investigated via the comparison with reanalysis and observational data. The results reveal strong and persistent agreement with reanalysis data during all seasons within the year with well-simulated annual cycle and regional patterns independently of the nudging parameters that were tested. Thus, the proposed formulation of the nudging provides a reliable framework for future long-term experiments aiming at hindcasting climate variability in the North Atlantic wave field. At the same time, dynamical downscaling allows for simulation of higher waves in coastal regions, specifically near the Greenland east coast likely due to a better representation of the mesoscale atmospheric dynamics in this area. 相似文献
11.
The main goal of this work is to appraise the finite element method in the way it represents barotropic instabilities. To that end, three different formulations are employed. The free-surface formulation solves the primitive shallow-water equations and is of predominant use for ocean modeling. The vorticity–stream function and velocity–pressure formulations resort to the rigid-lid approximation and are presented because theoretical results are based on the same approximation. The growth rates for all three formulations are compared for hyperbolic tangent and piecewise linear shear flows. Structured and unstructured meshes are utilized. The investigation is also extended to time scales that allow for instability meanders to unfold, permitting the formation of eddies. We find that all three finite element formulations accurately represent barotropic instablities. In particular, convergence of growth rates toward theoretical ones is observed in all cases. It is also shown that the use of unstructured meshes allows for decreasing the computational cost while achieving greater accuracy. Overall, we find that the finite element method for free-surface models is effective at representing barotropic instabilities when it is combined with an appropriate advection scheme and, most importantly, adapted meshes. 相似文献
12.
On the solution of transient free-surface flow problems in porous media by the finite element method
A numerical procedure is presented to deal with solution of transient free-surface flows in porous media. The governing boundary-value problem for the piezometric potential is solved by the finite element method. The initial-value problem which describes the transient motion of the free-surface is solved by the method of quasi-linearization. The numerical scheme has been applied to isotropic and anisotropic earth dam problem and also to a ditch drainage problem. Excellent agreements have been reached when compared with known solutions. This computational procedure is shown to be stable and suitable for this class of problems with the aid of a digital computer. 相似文献
13.
On the importance of non-hydrostatic processes in determining tidally induced mixing in sill regions
The importance of using a non-hydrostatic model to compute tidally induced mixing and flow in the region of a sill is examined using idealized topography representing the sill at the entrance to Loch Etive. This site is chosen since detailed measurements were recently made there. Calculations are performed with and without the inclusion of non-hydrostatic dynamics using a vertical slice model for a range of sill widths corresponding to typical sill regions. Initial non-hydrostatic calculations showed that the model could reproduce the observed flow characteristics in the region. However, when calculations were performed using the model in hydrostatic form, the significant artificial convective mixing that occurred in order to remove density inversions led to excessively high vertical mixing. This influenced the computed temperature field and the intensity of the current jet that separated from the sill on its lee side. In addition it affected the magnitude and spatial characteristics of the lee waves generated on the lee side of the sill. Calculations with a range of sill widths, showed that as the sill width decreased the difference between the solution computed with the non-hydrostatic and hydrostatic model increased. 相似文献
14.
在深水陡坡带和复杂海底构造区域,即使采用长电缆,由绕射产生的自由表面相关多次波因为传播方向的异变,仍然难以避免在自由表面发生下行反射的反射点位置位于电缆长度范围之外,造成预测多次波的近偏移距信息不足或者负偏移距信息缺失等问题,以致不能准确地预测和压制陡坡带的绕射多次波.本文在常规反馈迭代自由表面多次波压制方法的基础上做出相应改进,将邻炮的炮记录依照炮点、检波点互换原理引入到反馈迭代循环当中,补充绕射多次波预测所需要的近偏移距和负偏移距信息.经过Sigsbee2b模型和实际资料测试,达到较好压制深水陡坡带的自由表面相关绕射多次波的目地. 相似文献
15.
The effect of random recharge on uniform steady free-surface flow in heterogeneous porous formations
The effect of parametric uncertainty in recharge rate and spatial variability of hydraulic conductivity upon free-surface flow is investigated in a stochastic framework. We examine the three-dimensional free-surface gravitational flow problem for sloped mean uniform flow in a randomly heterogeneous porous medium under the influence of random recharge. We develop analytic solutions for the variance of free-surface position, head, and specific discharge on the free surface. Additionally, we obtain semi-analytic solutions for the statistical moments of head and specific discharge beneath the free-surface. Statistical moments are derived using a first-order approximation and then compared with their parallel in an unbounded medium. The effect of recharge mean and variability on the statistical moments is analyzed. Results can be applied to more complex flows, slowly varying in the mean. 相似文献
16.
J. Wolf 《Continental Shelf Research》1982,2(4):215-229
As the capacity of computers increases, the size and resolution of numerical models can be increased. In tidal models, however, using an explicit scheme together with high spatial resolution results in an unreasonably small time-step demanded by the CFL condition for stability. This condition is usually the most restrictive and applies to the propagation of the gravity wave.A fully nonlinear three-dimensional model has been developed, using the Galerkin method in the vertical dimension, in which the gravity wave terms are treated by an alternating-direction implicit scheme, while the friction, viscous and advective terms are treated explicitly. This permits a stable solution with a longer time-step than that required in an equivalent explicit scheme, while not needing as much computational effort as a fully implicit scheme. This semi-implicit model is compared with an explicit model in terms of efficiency, accuracy, and stability. Tidal and wind-driven flows and free oscillations in a rectangular sea model are examined, using various boundary conditions. The semi-implicit scheme takes approximately 1.1 times as long to run (per time-step) on the CRAY-1 computer as the explicit scheme. 相似文献
17.
Sheet flows occur widely in natural free-surface flows including rivers in flood, tidal estuaries and coastal waters in storm conditions when bed shear stress becomes sufficiently high. Particle volumetric concentration in sheet flows normally follows a linear distribution with the Rouse [Rouse H. Modern conceptions of the mechanics of fluid turbulence. Trans ASCE, 1937;102:463–543] distribution applicable in the dilute water column above the sheet-flow layer. However, a well-verified formula for determining particle velocity distribution in sheet flows is still lacking. This paper presents formulas to describe the particle velocity profile in steady or oscillatory sheet flows. They compare well with measured data. In particular, the novel formula for determining the particle velocity at the top of bedload–sediment-dominated sublayer in sheet flows is also well verified with measured data. 相似文献
18.
Irrespective of their spatial extent, free-surface shallow flows are challenging measurement environments for most instruments due to the relatively small depths and velocities typically associated with these flows. A promising candidate for enabling measurements in such conditions is Large-scale Particle Image Velocimetry (LSPIV). This technique uses a non-intrusive approach to measure two-dimensional surface velocity fields with high spatial and temporal resolutions. Although there are many publications documenting the successful use of LSPIV in various laboratory and field open-channel flow situations, its performance has not been equally substantiated for measurement in shallow flows. This paper aims at filling in this gap by demonstrating the capabilities of LSPIV to: (a) accurately evaluate complex flow patterns in shallow channel flows; and (b) estimate depth in shallow flows using exclusively LSPIV measurements. The demonstration is provided by LSPIV measurements in three shallow flow laboratory situations with flow depths ranging from 0.05 to 0.31 m. The obtained measurements illustrate the LSPIV flexibility and reliability in measuring velocities in shallow and low-velocity (near-zero) flows. Moreover, the technique is capable to evaluate and map velocity-derived quantities that are difficult to document with alternative measurement techniques (e.g. vorticity and shear stress distributions and mapping of large-scale structure in the body of water). 相似文献
19.
Numerical studies of flow over a sill: sensitivity of the non-hydrostatic effects to the grid size 总被引：1，自引：1，他引：0
A non-hydrostatic terrain-following model in cross sectional form is applied to study the processes in the lee of a sill in
an idealized stratified fjord during super-critical tidal inflow. A sequence of numerical studies with horizontal grid sizes
in the range from 100 to 1.5625 m are performed. All experiments are repeated using both hydrostatic and non-hydrostatic versions
of the model, allowing a systematic study of possible non-hydrostatic pressure effects and also of the sensitivity of these
effects to the horizontal grid size. The length scales and periods of the internal waves in the lee of the sill are gradually
reduced and the amplitudes of these waves are increased as the grid size is reduced from 100 down to 12.5 m. With a further
reduction in grid size, more short time and space scale motions become superimposed on the internal waves. Associated with
the internal wave activity, there is a deeper separation point that is fairly robust to all parameters investigated. Another
separation point nearer to the top of the sill appears in the numerical results from the high-resolution studies with the
non-hydrostatic model. Associated with this shallower separation point, an overturning vortex appears in the same set of numerical
solutions. This vortex grows in strength with reduced grid size in the non-hydrostatic experiments. The effects of the non-hydrostatic
pressure on the velocity and temperature fields grow with reduced grid size. In the experiments with horizontal grid sizes
equal to 100 or 50 m, the non-hydrostatic pressure effects are small. For smaller grid sizes, the time mean velocity and temperature
fields are also clearly affected by the non-hydrostatic pressure adjustments. 相似文献
20.
During seismic wave propagation on a free surface, a strong material contrast boundary develops in response to interference by P- and S- waves to create a surfacewave phenomenon. To accurately determine the effects of this interface on surface-wave propagation, the boundary conditions must be accurately modeled. In this paper, we present a numerical approach based on the dynamic poroelasticity for a space–time-domain staggeredgrid finite-difference simulation in porous media that contain a free-surface boundary. We propose a generalized stess mirror formulation of the free-surface boundary for solids and fluids in porous media for the grid mesh on which lays the free-surface plane. Its analog is that used for elastic media, which is suitable for precise and stable Rayleigh-type surface-wave modeling. The results of our analysis of first kind of Rayleigh (R1) waves obtained by this model demonstrate that the discretization of the mesh in a similar way to that for elastic media can realize stable numerical solutions with acceptable precision. We present numerical examples demonstrating the efficiency and accuracy of our proposed method. 相似文献