共查询到20条相似文献，搜索用时 15 毫秒
1.
A higher-efficient three-dimensional numerical model for small amplitude free surface flows 总被引：1，自引：0，他引：1
A higher-efficient three-dimensional non-hydrostatic model is developed to simulate small amplitude free surface flows based on a staggered unstructured grid. In this model, a fractional step algorithm is adopted to solve the Navier-Stokes equations in two major steps. A top-layer pressure method is proposed to minimize the number of vertical layers and subsequently the computational cost. Three classical examples of small amplitude free surface flows are used to demonstrate the capability and efficiency of the model. The satisfactory results demonstrated the capability and efficiency of modelling a range of small amplitude free surface flows with only a small number of vertical layers. 相似文献
2.
This paper presents CCHE2D-NHWAVE, a depth-integrated non-hydrostatic finite element model for simulating nearshore wave processes. The governing equations are a depth-integrated vertical momentum equation and the shallow water equations including extra non-hydrostatic pressure terms, which enable the model to simulate relatively short wave motions, where both frequency dispersion and nonlinear effects play important roles. A special type of finite element method, which was previously developed for a well-validated depth-integrated free surface flow model CCHE2D, is used to solve the governing equations on a partially staggered grid using a pressure projection method. To resolve discontinuous flows, involving breaking waves and hydraulic jumps, a momentum conservation advection scheme is developed based on the partially staggered grid. In addition, a simple and efficient wetting and drying algorithm is implemented to deal with the moving shoreline. The model is first verified by analytical solutions, and then validated by a series of laboratory experiments. The comparison shows that the developed wave model without the use of any empirical parameters is capable of accurately simulating a wide range of nearshore wave processes, including propagation, breaking, and run-up of nonlinear dispersive waves and transformation and inundation of tsunami waves. 相似文献
3.
A higher-order non-hydrostatic model in a σ-coordinate system is developed. The model uses an implicit finite difference scheme on a staggered grid to simultaneously solve the unsteady Navier-Stokes equations (NSE) with the free-surface boundary conditions. An integral method is applied to resolve the top-layer non-hydrostatic pressure, allowing for accurately resolving free-surface wave propagation. In contrast to the previous work, a higher-order spatial discretization is utilized to approximate the large horizontal pressure gradient due to steep surface waves or rapidly varying topographies. An efficient direct solver is developed to solve the resulting block hepta-diagonal matrix system. Accuracy of the new model is validated by linear and nonlinear standing waves and progressive waves. The model is then used to examine freak (extreme) waves. Features of downshifting focusing location and wave asymmetry characteristics are predicted on the temporal and spatial domains of a freak wave. 相似文献
4.
A new three-dimensional, non-hydrostatic free surface flow model is presented. For simulating water wave motions over uneven bottoms, the model employs an explicit project method on a Cartesian the staggered gird system to solve the complete three-dimensional Navier–Stokes equations. A bi-conjugated gradient method with a pre-conditioning procedure is used to solve the resulting matrix system. The model is capable of resolving non-hydrostatic pressure by incorporating the integral method of the top-layer pressure treatment, and predicting wave propagation and interaction over irregular bottom by including a partial bottom-cell treatment. Four examples of surface wave propagation are used to demonstrate the capability of the model. Using a small of vertical layers (e.g. 2–3 layers), it is shown that the model could effectively and accurately resolve wave shoaling, non-linearity, dispersion, fission, refraction, and diffraction phenomena. 相似文献
5.
Hydrodynamic forces exerting on a pipeline partially buried in a permeable seabed subjected to combined oscillatory flow and steady current are investigated numerically. Two-dimensional Reynolds-Averaged Navier-Stokes equations with a k−ω turbulent model closure are solved to simulate the flow around the pipeline. The Laplace equation is solved to calculate the pore pressure below the seabed with the simulated seabed hydrodynamic pressure as boundary conditions. The numerical model is validated against the experimental data of a fully exposed pipeline resting on a plane boundary under various flow conditions. Then the flow with different embedment depths, steady current ratios and KC numbers is simulated. The amplitude of seepage velocity is much smaller than the amplitude of free stream velocity as expected. The normalized Morison inertia, drag and lift coefficients based on the corresponding force coefficients of a fully exposed pipeline are investigated. The normalized Morison force coefficients reduce almost linearly with the increase of embedment depth and that the KC only has minor effect on the normalized Morison coefficients. It is also found that the permeable seabed condition causes a slight increase on the inline force and has a little effect on the lift force, compared with corresponding conditions in an impermeable bed. 相似文献
6.
《Coastal Engineering》2004,51(7):557-579
In this paper, a Reynolds Averaged Navier–Stokes (RANS) model was developed to simulate the vortex generation and dissipation caused by progressive waves passing over impermeable submerged double breakwaters. The dynamics of the turbulence are described by introducing a k–ɛ model with Boussinesq closure. The Height Function (HF) is implemented to define the free-surface configuration. The governing equations are discretized by means of a finite volume method based on a staggered grid system with variable width and height. The feasibility of the numerical model was verified through a series of comparisons of numerical results with the existing analytical solutions and the experimental data. The good agreements demonstrate the satisfactory performance of the developed numerical model. The flow separation mechanism both near the upstream and the downstream edges of the obstacles demonstrates the physical and expected nature of development of the flow. The present model provides an accurate and efficient tool for the simulation of flow field and wave transformation near coastal structures without breaking. 相似文献
7.
A vertical two-dimensional numerical model has been applied to solving the Reynolds Averaged Navier- Stokes （RANS} equations in the simulation of current and wave propagation through vegetated and non- vegetated waters. The k-e model is used for turbulence closure of RANS equations. The effect of vegeta- tion is simulated by adding the drag force of vegetation in the flow momentum equations and turbulence model. To solve the modified N-S equations, the finite difference method is used with the staggered grid system to solver equations. The Youngs＇ fractional volume of fluid （VOF） is applied tracking the free sur- face with second-order accuracy. The model has been tested by simulating dam break wave, pure current with vegetation, solitary wave runup on vegetated and non-vegetated channel, regular and random waves over a vegetated field. The model reasonably well reproduces these experimental observations, the model- ing approach presented herein should be useful in simulating nearshore processes in coastal domains with vegetation effects. 相似文献
8.
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. 相似文献
9.
In this study an Euler-Euler two-phase model was developed to investigate the tunnel erosion beneath a submarine pipeline exposed to unidirectional flow. Both of the fluid and sediment phases were described via the Navier-Stokes equations, i.e. the model was implemented using time-averaged continuity and momentum equations for the fluid and sediment phases and a modified k−ε turbulence closure for the fluid phase. The fluid and sediment phases were coupled by considering the drag and lift interaction forces. The model was employed to simulate the tunnel erosion around the pipeline laid on an erodible bed. Comparison between the numerical result and experimental measurement confirms that the numerical model successfully predicts the bed profile and velocity field during the tunnel erosion. It is evident that the sediments are transported as the sheet-flow mode in the tunnel erosion stage. Also the transport rate under the pipe increases rapidly at the early stage and then reduces gradually at the end of the tunnel erosion beneath pipelines. 相似文献
10.
A new three-dimensional semi-implicit finite-volume ocean model has been developed for simulating the coastal ocean circulation, which is based on the staggered C -unstructured non-orthogonal grid in the horizontal direction and z -level grid in the vertical direction. The three-dimensional model is discretized by the semi-implicit finite-volume method, in that the free-surface and the vertical diffusion are semi-implicit, thereby removing stability limitations associated with the surface gravity wave and vertical diffusion terms. The remaining terms in the momentum equations are discretized explicitly by an integral method. The partial cell method is used for resolving topography, which enables the model to better represent irregular topography. The model has been tested against analytical cases for wind and tidal oscillation circulation, and is applied to simulating the tidal flow in the Bohai Sea. The results are in good agreement both with the analytical solutions and measurement results. 相似文献
11.
A higher-order non-hydrostatic σ model is developed to simulate non-linear refraction–diffraction of water waves. To capture non-linear (or steep) waves, a 4th-order spatial discretization is utilized to approximate the large horizontal pressure gradient. A higher-order top-layer pressure treatment is further implemented to resolve wave propagation. The model's characteristics including linear wave dispersion and non-linearity are carefully examined. The accuracy of the present model using only two vertical layers is validated by laboratory data and the available results predicted by the non-linear Schrödinger equation, Boussinesq-type equations, the non-linear mild slope equation, and the Laplace equation. Features of harmonic generation as well as the influences of dispersion and non-linearity on wave energy transfer processes are discussed. 相似文献
12.
A numerical model is developed by combining a porous flow model and a two-phase flow model to simulate wave transformation in porous structure and hydraulic performances of a composite type low-crest seawall. The structure consists of a wide submerged reef, a porous terrace at the top and an impermeable rear wall. The porous flow model is based on the extended Navier-Stokes equations for wave motion in porous media and k−ε turbulence equations. The two-phase flow model combines the water domain with the air zone of finite thickness above water surface. A unique solution domain is established by satisfying kinematic boundary condition at the interface of air and water. The free surface advection of water wave is modeled by the volume of fluid method with newly developed fluid advection algorithm. Comparison of computed and measured wave properties shows reasonably good agreement. The influence of terrace width and structure porosity is investigated based on numerical results. It is concluded that there exist optimum value of terrace width and porosity that can maximize hydraulic performances. The velocity distributions inside and in front of the structure are also investigated. 相似文献
13.
S. N. Moshonkin V. B. Zalesny A. V. Gusev 《Izvestiya Atmospheric and Oceanic Physics》2018,54(5):495-506
The algorithm for splitting k–ω turbulence equations is used to parameterize viscosity and diffusion coefficients in the ocean general circulation model. The k–ω equations are split into stages describing the transport-diffusion and generation-dissipation of the turbulent kinetic energy and frequency function ω. At the generation-dissipation stage, the equations are solved analytically. Calculations of circulation in the North Atlantic–Arctic Ocean for 1948–2009 have been carried out. The experiments demonstrate an adequate reproduction of hydrophysical characteristics and high efficiency of the algorithm. It is shown that considering the climatic annual mean buoyancy frequency in the turbulence equations at the generation-dissipation stage is an important factor in improving the accuracy of simulated fields. 相似文献
14.
Some changes in the eddy mixing in the atmospheric boundary layer (ABL) are investigated with the use of the mesoscale RANS turbulence model. It is found that the behavior of parameters of the eddy turbulence mixing is in compliance with the recently obtained data of laboratory and atmospheric measurements. In particular, the flow Richardson number (Ri f ) during the transient flow to a strongly stable state can behave nonmonotonically, growing with the increasing gradient Richardson number (Ri g ) to the state of saturation at a certain gradient Richardson number (Ri g ? 1), which separates two different turbulent regimes: the regimes of strong mixing and weak mixing. An analysis of the energetics based on the balance equations of kinetic and potential turbulence energies shows, in particular, that the weak mixing (Ri g > 1) is quite capable of transferring momentum. This phenomenon can be explained not only by the fact that the flow is sustained by propagating internal waves, which effectively transfer momentum under strong stratification conditions, but also by the fact that turbulence permanently arises in the free atmosphere and in the deep ocean at Ri g ? 1. 相似文献
15.
A modified three-parameter model of turbulence for a thermally stratified atmospheric boundary layer (ABL) is presented. The model is based on tensor-invariant parametrizations for the pressure-strain and pressure-temperature correlations that are more complete than the parametrizations used in the Mellor-Yamada model of level 3.0. The turbulent momentum and heat fluxes are calculated with explicit algebraic models obtained with the aid of symbol algebra from the transport equations for momentum and heat fluxes in the approximation of weakly equilibrium turbulence. The turbulent transport of heat and momentum fluxes is assumed to be negligibly small in this approximation. The three-parameter E ? ε ? 2> model of thermally stratified turbulence is employed to obtain closed-form algebraic expressions for the fluxes. A computational test of a 24-h ABL evolution is implemented for an idealized two-dimensional region. Comparison of the computed results with the available observational data and other numerical models shows that the proposed model is able to reproduce both the most important structural features of the turbulence in an urban canopy layer near the urbanized ABL surface and the effect of urban roughness on a global structure of the fields of wind and temperature over a city. The results of the computational test for the new model indicate that the motion of air in the urban canopy layer is strongly influenced by mechanical factors (buildings) and thermal stratification. 相似文献
16.
It is well known that the sudden intrusion of Kuroshio warm water into the Bungo Channel (kyucho) is regulated by spring–neap tidal forcing. In order to clarify the physical background behind this regulation, numerical experiments are carried out using a high-resolution non-hydrostatic three-dimensional model. We first reproduce the strong mixing region off the east coast of the Bungo Channel resulting from tidal flow interaction with complicated land configurations during spring tides; behind islands and headlands, small-scale eddies satisfying an approximate cyclostrophic balance are generated. As a result, averaged over the whole model domain, the tidal-mean energy dissipation rate reaches ≈1.6?×?10?6?W?kg?1. The model predicted energy dissipation rates at the location and times of direct microstructure measurements in the Bungo Channel are comparable to the observed values. We next examine whether or not strong tidal mixing thus reproduced can inhibit the northward intrusion of Kuroshio warm water in the Bungo Channel. It is shown that the Kuroshio warm water can (or cannot) pass through the tidal mixing regions off the east coast of the Bungo Channel during periods of weakened (or enhanced) tidal mixing at neap (or spring) tides. This indicates that taking into account the realistic spring–neap modulation of tidal mixing intensity is indispensable to further increase the ability of the existing forecast system for kyucho in the Bungo Channel. 相似文献
17.
Vegetation canopies control mean and turbulent flow structure as well as surface wave processes in coastal regions. A non-hydrostatic RANS model based on NHWAVE (Ma et al., 2012) is developed to study turbulent mixing, surface wave attenuation and nearshore circulation induced by vegetation. A nonlinear k − ϵ model accounting for vegetation-induced turbulence production is implemented to study turbulent flow within the vegetation field. The model is calibrated and validated using experimental data from vegetated open channel flow, as well as nonbreaking and breaking random wave propagation in vegetation fields. It is found that the drag-related coefficients in the k − ϵ model Cfk and Cfϵ can greatly affect turbulent flow structure, but seldom change the wave attenuation rate. The bulk drag coefficient CD is the major parameter controlling surface wave damping by vegetation canopies. Using the empirical formula of Mendez and Losada (2004), the present model provides accurate predictions of vegetation-induced wave energy dissipation. Wave propagation through a finite patch of vegetation in the surf zone is investigated as well. It is found that the presence of a finite patch of vegetation may generate strong pressure-driven nearshore currents, with an onshore mean flow in the unvegetated zone and an offshore return flow in the vegetated zone. 相似文献
18.
An unstructured-grid, finite-volume, nonhydrostatic, parallel coastal ocean simulator 总被引：4，自引：4，他引：4
A finite-volume formulation is presented that solves the three-dimensional, nonhydrostatic Navier–Stokes equations with the Boussinesq approximation on an unstructured, staggered, z-level grid, with the goal of simulating nonhydrostatic processes in the coastal ocean with grid resolutions of tens of meters. In particular, the code has been developed to simulate the nonlinear, nonhydrostatic internal wave field in the littoral ocean. The method is based on the formulation developed by Casulli, in that the free-surface and vertical diffusion are semi-implicit, thereby removing stability limitations associated with the surface gravity wave and vertical diffusion terms. The remaining terms in the momentum equations are discretized explicitly with the second-order Adams–Bashforth method, while the pressure-correction method is employed for the nonhydrostatic pressure in order to achieve overall second-order temporal accuracy. Advection of momentum is accomplished with an Eulerian discretization which conserves momentum in cells that do not contain the free surface, and scalar advection is discretized in a way that ensures consistency with continuity, thereby ensuring local and global mass conservation using a velocity field that conserves volume on a local and global basis. The nonhydrostatic pressure field is solved efficiently using a block-Jacobi preconditioner, and while stability is limited by the internal gravity wave speed and vertical advection of momentum, applications requiring relatively small time steps due to accuracy or stability constraints are run efficiently on parallel computers, since the present formulation is written entirely with the message-passing interface (MPI). The ParMETIS libraries are employed in order to achieve a load-balanced parallel partitioning that minimizes interprocessor communication, and the grid is reordered to optimize per-processor performance by limiting cache misses while accessing arrays in memory. Test cases demonstrate the ability of the code to efficiently and accurately compute the nonhydrostatic lock exchange and internal waves in idealized as well as real domains, and we evaluate the parallel efficiency of the code using up to 32 processors. 相似文献
19.
20.
Efficient computation of surf zone waves using the nonlinear shallow water equations with non-hydrostatic pressure 总被引：2，自引：0，他引：2
A numerical method for non-hydrostatic, free-surface, irrotational flow governed by the nonlinear shallow water equations including the effects of vertical acceleration is presented at the aim of studying surf zone phenomena. A vertical boundary-fitted grid is used with the water depth divided into a number of layers. A compact finite difference scheme is employed for accurate computation of frequency dispersion requiring a limited vertical resolution and hence, capable of predicting the onset of wave breaking. A novel wet–dry algorithm is applied for a proper handling of moving shoreline. Mass and momentum are strictly conserved at discrete level while the method only dissipates energy in the case of wave breaking. The numerical results are verified with a number of tests and show that the proposed model using two layers without ad-hoc assumptions enables to resolve propagating nonlinear shoaling, breaking waves and wave run-up within the surf and swash zones in an efficient manner. 相似文献