共查询到20条相似文献，搜索用时 31 毫秒
1.
This paper describes the application of the Smoothed Particle Hydrodynamics（SPH） method for modeling two dimensional waves caused by dam break over a movable bed in two dimensions.The two phase SPH method is developed to solve the Navier-Stokes equations.Both fluid and sediment phases are described by particles as weakly compressible fluids and the incompressibility is achieved by the equation of state.The sediment phase is modeled as a non-Newtonian fluid using three alternative approaches of artificial viscosity and Bingham Model.In this paper,the new formulations for two-phase flows are proposed.The numerical results obtained from the developed SPH model show acceptable accuracy with comparison to experimental data. 相似文献
2.
In this study an incompressible smoothed particle hydrodynamics (ISPH) approach coupled with the sediment erosion model is developed to investigate the sediment bed scour and grain movement under the dam break flows. Two-phase formulations are used in the ISPH numerical algorithms to examine the free surface and bed evolution profiles, in which the entrained sediments are treated as a different fluid component as compared with the water. The sediment bed erosion model is based on the concept of pick-up flow velocity and the sediment is initiated when the local flow velocity exceeds a critical value. The proposed model is used to reproduce the sediment erosion and follow-on entrainment process under an instantaneous dam break flow and the results are compared with those from the weakly compressible moving particle semi-implicit (WCMPS) method as well as the experimental data. It has been demonstrated that the two-phase ISPH model performed well with the experimental data. The study shows that the ISPH modelling approach can accurately predict the dynamic sediment scouring process without the need to use empirical sediment transport formulas. 相似文献
3.
Ocean Dynamics - In this paper, a two-phase model based on the SPH method is established using the Newtonian fluid and non-Newtonian fluid to describe water and sediment respectively. In this... 相似文献
4.
Ahmad ShakibaeiniaYee-Chung Jin 《Advances in water resources》2011,34(6):794-807
Mesh-free particle (Lagrangian) methods such as Moving Particle Semi-Implicit (MPS) and Smoothed Particle Hydrodynamics (SPH) are the latest generation of methods in the field of computational fluid dynamics that attracts lots of attention in modeling applications where large interfacial deformations and fragmentations exist. Due to their mesh-free nature, these methods are capable of simulating any kind of boundary/interface deformation and fragmentations. This study aims to develop a new mesh-free particle model based on the weakly compressible MPS (WC-MPS) formulation for modeling of a dam break over a mobile bed, which is a highly erosive and transient flow problem. A multiphase model, capable of handling the density and viscosity discontinuity and in which the solid (sediment) phase is treated as a non-Newtonian fluid, is introduced. The resulting model is first validated using a two-phase dam break problem and is then applied to the mobile-bed dam break problem with different conditions, comparing the results to those obtained from some experimental works. 相似文献
5.
《Advances in water resources》1998,21(5):351-362
Flow and displacement of non-Newtonian fluids in porous media occurs in many subsurface systems, related to underground natural resource recovery and storage projects, as well as environmental remediation schemes. A thorough understanding of non-Newtonian fluid flow through porous media is of fundamental importance in these engineering applications. Considerable progress has been made in our understanding of single-phase porous flow behavior of non-Newtonian fluids through many quantitative and experimental studies over the past few decades. However, very little research can be found in the literature regarding multi-phase non-Newtonian fluid flow or numerical modeling approaches for such analyses.For non-Newtonian fluid flow through porous media, the governing equations become nonlinear, even under single-phase flow conditions, because effective viscosity for the non-Newtonian fluid is a highly nonlinear function of the shear rate, or the pore velocity. The solution for such problems can in general only be obtained by numerical methods.We have developed a three-dimensional, fully implicit, integral finite difference simulator for single- and multi-phase flow of non-Newtonian fluids in porous/fractured media. The methodology, architecture and numerical scheme of the model are based on a general multi-phase, multi-component fluid and heat flow simulator — TOUGH2. Several rheological models for power-law and Bingham non-Newtonian fluids have been incorporated into the model. In addition, the model predictions on single- and multi-phase flow of the power-law and Bingham fluids have been verified against the analytical solutions available for these problems, and in all the cases the numerical simulations are in good agreement with the analytical solutions. In this presentation, we will discuss the numerical scheme used in the treatment of non-Newtonian properties, and several benchmark problems for model verification.In an effort to demonstrate the three-dimensional modeling capability of the model, a three-dimensional, two-phase flow example is also presented to examine the model results using laboratory and simulation results existing for the three-dimensional problem with Newtonian fluid flow. 相似文献
6.
将时空守恒元/解元(CE/SE)方法推广到二维孔隙介质多相流问题的数值计算中,采用人工压缩法耦合速度和压力,同时结合杂交粒子水平集方法捕捉物质界面.提出一套完整的二维欧拉型孔隙介质非稳态多相不可压缩黏性流动计算方案.通过对溃坝和液滴在重力作用下的运动和变形问题的数值模拟,验证了方法的精度和有效性.在此基础上,提出了一个新的孔隙介质两相流物理模型——双层流体顶盖驱动方腔流. 相似文献
7.
The first part of the research reported here consists of an experimental campaign to study the scouring of a granular bed(glass beads, sand) induced by a dam break in an open channel. Two configurations are considered: with and without cylinders. In the second part of this study, the volume of fluid method coupled with the shear stress transport turbulent model and the lagrangian particle tracking method is used to simulate the local scour processes. The four-way coupling is realized by consider... 相似文献
8.
The current study proposes a novel framework for the numerical model for estimating the temporal scour considering unsteady sediment inflow and the sediment sorting process. The framework was applied to local scour upstream of a slit weir. The scour model is based on an ordinary nonlinear differential equation derived from sediment continuity and scour rate equations. A one-dimensional(1-D)Boussinesq-type model coupled with nonequilibrium sediment transport was incorporated in the scour model to... 相似文献
9.
《国际泥沙研究》2021,36(6):696-710
Interactions between hydrokinetic turbines and near scale and far scale bed sediment particles are considered a critical area of assessment; however, a limited number of research studies have been published to address this issue. The current paper explores the influence of the blockage ratio on the morphology evolution downstream from a hydrokinetic turbine. A modeling framework is derived to predict the significant transport induced by a turbine installed on the erodible fluvial sandy bed surface, such as the Eulerian multi-phase model for the sediment and the Blade Element Momentum Theory (BEMT) for the turbine, using the open source platform OpenFOAM. Two configurations of different blockage ratios are considered. The configurations show differences in the bed morphology and the sediment transport. It has been shown for both configurations that the scour capabilities are enhanced below the turbine due to the acceleration of the flow and increasing local shear stress on the sediment. 相似文献
10.
Transport of suspended sediment under the dam‐break flow on an inclined plane bed of arbitrary slope
The problem of transport of suspended sediment after the break of a dam on an inclined bed is considered. To that end we use the shallow‐water approximation for arbitrary, constant slopes of the bottom, taking into consideration the effect of friction. The numerical technique and the frictional model are validated by comparison with available experimental data and asymptotic analytical solutions, with special attention to the numerical solution near the wetting front. The transport of suspended sediment down the inclined bed is obtained and discussed as a function of the slope of the bed for different values of the parameters characterizing the sediment and its transport properties. For sufficiently large times we always find the formation of roll waves near the water front, which affects the transport of sediments significantly. These strong oscillations are accurately computed with the numerical method used. The relative importance of the bed load (to the suspended load) sediment transport is also discussed as a function of the size of the sediment particles and the slope of the bed for different models on the initiation of sediment suspension from bed load. We also check the dilute sediment approach and characterize the conditions for its failure. Finally, the results of the present simplified model are intended to be used as tests of more complex numerical models. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
11.
Details are given of the development and application of a numerical model for predicting free-surface flows in estuarine and coastal basins using the finite volume method. Both second- and third-order accurate and oscillation free explicit numerical schemes have been used to solve the shallow water equations. The model deploys an unstructured triangular mesh and incorporates two types of mesh layouts, namely the ‘cell centred’ and ‘mesh vertex’ layouts, and provides a powerful mesh generator in which a user can adjust the mesh-size distribution interactively to create a desirable mesh. The quality of mesh has been shown to have a major impact on the overall performance of the numerical model.The model has been applied to simulate two-dimensional dam break flows for which transient water level distributions measured within a laboratory flume were available. In total 12 model runs were undertaken to test the model for various flow conditions. These conditions include: (1) different bed slopes (ranging from zero to 0.8%), (2) different upstream and downstream water level conditions, and (3) initially wet and dry bed conditions, downstream of the dam. Detailed comparisons have been made between model predicted and measured water levels and good agreement achieved between both sets of results. The model was then used to predict water level and velocity distributions in a real estuary, i.e. the Ribble Estuary, where the bed level varies rapidly at certain locations. In order to model the whole estuary, a 1-D numerical model has also been used to model the upper part of the estuary and this model was linked dynamically to the 2-D model. Findings from this application are given in detail. 相似文献
12.
After experiencing 8-day combined tidal current, circulation and wave actions, scour depth surrounding cylinder object freely resting on sandy seabed in the East China Sea （ECS） in January is numerically predicted using the DRAMBUIE model designed for scour burial, which has been widely used and verified by in-situ experiments. During the period of numerical integration, the value of time t is generally variable at every time step via the special time-stepped approach developed by this paper to eliminate the time error. The tidal current velocity, wave orbital velocity and the depth-averaged circulation in the ECS have been obtained by numerical simulations with Estuarine Coastal and Ocean Model （ECOM）, Simulating Waves Nearshore （SWAN） model and Regional Ocean Modeling System （ROMS） model respectively. The control experiment and several idealized test cases on influential factors in scour depth reveal that the dominant hydrodynamic factor is tidal current in the ECS under normal weather conditions, and the impacts of shelf circulation and wave motion on local scour almost can be ignored with an exception of the Kuroshio area where the high-speed mainstream of Kuroshio flows. It is also indicated that in sandy sediments, the distribution of scour depth nearly follows the pattern of tidal currents, while the secondary influencing factor on scour depth appears to be grain size of sandy sediment in the ECS. Numerical tests on sediment grain size further testify that much finer sand is more easily scoured, and an increasing trend for scour depth with reduction of grain size is displayed due to imposed resistance of larger sized particles. Three aspects explored by this paper, including the empirical equations in the Defense Research Agency Mine Burial Environment （DRAMBUIE） model, the accuracy of inputs and infill process can severely affect the prediction of scour depth surrounding cylinder objects freely resting on sandy seabed in the ECS. 相似文献
13.
《国际泥沙研究》2019,34(6):537-549
Dam removal can generate geomorphic disturbances, including channel bed and bank erosion and associated abrupt/pulsed release and downstream transfer of reservoir sediment, but the type and rate of geomorphic response often are hard to predict. The situation gets even more complex in systems which have been impacted by multiple dams and a long and complex engineering history. In previous studies one-dimensional (1-D) models were used to predict aspects of post-removal channel change. However, these models do not consider two-dimensional (2-D) effects of dam removal such as bank erosion processes and lateral migration. In the current study the impacts of multiple dams and their removal on channel evolution and sediment delivery were modeled by using a 2-D landscape evolution model (CAESAR-Lisflood) focusing on the following aspects: patterns, rates, and processes of geomorphic change and associated sediment delivery on annual to decadal timescales. The current modeling study revealed that geomorphic response to dam removal (i.e., channel evolution and associated rates of sediment delivery) in multiple dam settings is variable and complex in space and time. Complexity in geomorphic system response is related to differences in dam size, the proximity of upstream dams, related buffering effects and associated rates of upstream sediment supply, and emerging feedback processes as well as to the presence of channel stabilization measures. Modeled types and rates of geomorphic adjustment, using the 2-D landscape evolution model CAESAR-Lisflood, are similar to those reported in previous studies. Moreover, the use of a 2-D method showed some advantages compared to 1-D models, generating spatially varying patterns of erosion and deposition before and after dam removal that provide morphologies that are more readily comparable to field data as well as features like the lateral re-working of past reservoir deposits which further enables the maintenance of sediment delivery downstream. 相似文献
14.
Summary In the present paperTaylor's analysis of the dispersion of a soluble matter in Newtonian flow through a circular tube is extended in the case of non-Newtonian flows of Eyring and Reiner-Philippoff model fluids. It has been shown here that the results for the Newtonian fluid can be deduced from the corresponding results of both the two types non-Newtonian flows. Few specific cases of both the two types of fluids have been studied. Aris modification ofTaylor's analysis is also applicable to non-Newtonian flows discussed here. The results may be useful in connection with the study of the dispersion of soluble salts in blood vessels. It may also be useful to physicians who wish to study molecular diffusion coefficients. 相似文献
15.
This numerical investigation was carried out to advance mechanistic understanding of sediment transport under sheet flow conditions. An Euler–Euler coupled two-phase flow model was developed to simulate fluid–sediment oscillatory sheet flow. Since the concentration of sediment particles is high in such flows, the kinematics of the fluid and sediment phases are strongly coupled. This model includes interaction forces, intergranular stresses and turbulent stress closure. Each phase was modeled via the Reynolds-Averaged Navier–Stokes equations, with interphase momentum conservation accounting for the interaction between the phases. The generation and transformation of turbulence was modeled using the two-equation k–ε turbulence model. Concentration and sediment flux profiles were compared with experimental data for sheet flow conditions considering both symmetric and asymmetric oscillatory flows. Sediment and fluid velocity variations, concentration profiles, sediment flux and turbulence parameters of wave-generated sheet flow were studied numerically with a focus on sediment transport characteristics. In all applications, the model predictions compared well with the experimental data. Unlike previous investigations in which the flow is driven by a horizontal pressure gradient, the present model solves the Navier–Stokes equations under propagating waves. The model’s ability to predict sediment transport under oscillatory sheet flow conditions underscores its potential for understanding the evolution of beach morphology. 相似文献
16.
A two-dimensional lattice Boltzmann model (LBM) is presented for transient shallow water flows. The model is based on the shallow water equations coupled with the large eddy simulation model. In order to obtain accurate results efficiently, a multi-block lattice scheme is applied at the area where a local finer grid is needed for strong change in physical variables. The model is verified by applying to five cases with transient processes: (a) a tidal wave over steps; (b) a perturbation over a submerged hump; (c) partial dam break flow; (d) circular dam break flow; (e) interaction between a dam break surge and four square cylinders. The objectives of this study are to validate the two-dimensional LBM in transient flow simulation and provide the detailed transient processes in shallow water flows. 相似文献
17.
A two‐dimensional (2D) finite‐difference shallow water model based on a second‐order hybrid type of total variation diminishing (TVD) approximate solver with a MUSCL limiter function was developed to model flooding and inundation problems where the evolution of the drying and wetting interface is numerically challenging. Both a minimum positive depth (MPD) scheme and a non‐MPD scheme were employed to handle the advancement of drying and wetting fronts. We used several model problems to verify the model, including a dam break in a slope channel, a dam break flooding over a triangular obstacle, an idealized circular dam‐break, and a tide flow over a mound. Computed results agreed well with the experiment data and other numerical results available. The model was then applied to simulate the dam breaking and flooding of Hsindien Creek, Taiwan, with the detailed river basin topography. Computed flooding scenarios show reasonable flow characteristics. Though the average speed of flooding is 6–7 m s?1, which corresponds to the subcritical flow condition (Fr < 1), the local maximum speed of flooding is 14·12 m s?1, which corresponds to the supercritical flow condition (Fr ≈ 1·31). It is necessary to conduct some kind of comparison of the numerical results with measurements/experiments in further studies. Nevertheless, the model exhibits its capability to capture the essential features of dam‐break flows with drying and wetting fronts. It also exhibits the potential to provide the basis for computationally efficient flood routing and warning information. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
18.
Channel cross‐sectional changes since construction of Livingston Dam and Lake Livingston in 1968 were studied in the lower Trinity River, Texas, to test theoretical models of channel adjustment, and to determine controls on the spatial extent of channel response. High and average flows were not significantly modified by the dam, but sediment transport is greatly reduced. The study is treated as an opportunistic experiment to examine the effects of a reduction in sediment supply when discharge regime is unchanged. Channel scour is evident for about 60 km downstream, and the general phenomena of incision, widening, coarsening of channel sediment and a decrease in channel slope are successfully predicted, in a qualitative sense, by standard models of channel response. However, there is no consistent channel response within this reach, as various qualitatively different combinations of increases, decreases or no change in width, depth, slope and roughness occur. These multiple modes of adjustment are predicted by the unstable hydraulic geometry model. Between about 60 km and the Trinity delta 175 km downstream of the dam, no morphological response to the dam is observed. Rather than a diminution of the dam's effects on fluvial processes, this is due to a fundamental change in controls of the fluvial system. The downstream end of the scour zone corresponds to the upstream extent of channel response to Holocene sea level rise. Beyond 60 km downstream, the Trinity River is characterized by extensive sediment storage and reduced conveyance capacity, so that even after dam construction sediment supply still exceeds transport capacity. The channel bed of much of this reach is near or below sea level, so that sea level rise and backwater effects from the estuary are more important controls on the fluvial system than upstream inputs. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
19.
１ＩＮＴＲＯＤＵＣＴＩＯＮＲｉｖｅｒｓｉｎＴａｉｗａｎａｒｅｒｅｌａｔｉｖｅｌｙｓｔｅｅｐｅｒｃｏｍｐａｒｅｄｔｏｔｈｏｓｅｉｎｏｔｈｅｒｃｏｎｔｉｎｅｎｔ．Ｌｏｃａｌｙｏｃｕｒｅｄｓｕｐｅｒｃｒｉｔｉｃａｌｆｌｏｗａｒｅｆａｉｒｌｙｃｏｍｍｏｎｉｎ... 相似文献
20.
Numerical modeling of local scour of non-uniform graded sediment for two arrangements of pile groups
A three-dimensional(3D) non-hydrostatic numerical model is established to investigate local scour around four aligned circular piles in uniform and non-uniform sediment mixtures and to provide information for improving scour countermeasures design. In the current study, unsteady Reynolds averaged Navier-Stokes(URANS) equations along with a Re-normalization Group(RNG) k-ε model were applied to simulate the flow field. A non-uniform sediment transport model was applied to estimate the bedload tran... 相似文献