共查询到20条相似文献，搜索用时 31 毫秒
1.
The interaction between fluid and sediment particles is widely involved in hydraulic engineering problems. In the current study, an explicit incompressible mesh-free method in the framework of the Moving Particle Semi-implicit(MPS) method is proposed to simulate the interaction between the two phases in submerged conditions. The proposed method solves two sets of the continuity and momentum equations, respectively, for the fluid phase and the sediment phase according to the mixture theory. In th... 相似文献
2.
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. 相似文献
3.
Multiple-point geostatistical simulation of dykes: application at Sungun porphyry copper system,Iran 总被引：5，自引：3，他引：2
Hassan Rezaee Omid Asghari Mohammad Koneshloo Julián M. Ortiz 《Stochastic Environmental Research and Risk Assessment (SERRA)》2014,28(7):1913-1927
Variogram-based methods are not capable of capturing high (>2) order statistics since the variogram measures the relationship between two points at a time only. Multiple-point geostatistics (MPS) has brought new insights into many geological modeling problems. The application of MPS methods has been well documented in realizing complex geological patterns. These methods have often been used in reservoir characterization since their advent in recent decades. The frequent non-linear behaviors of geologic continuity are not limited to reservoirs, but mineral deposits bear complicated formations in many cases. Relying on the power of MPS methods and considering the complexity of geological scenarios in mineral deposits, we have applied MPS in the modeling of mineral deposits. A training image (TI) is produced using geological data from upper horizons of a porphyry copper ore deposit which have been mined out during the previous mining operations. In this study, the SNESIM algorithm has been used. A number of realizations are produced using this multiple-point geostatistical method. Extensive validation steps are performed considering the TI as the reference model. These validations first show that the TI is representative for the domain under study and also illustrates some degrees of similarity between the TI and the realizations. Despite simplifications made to the problem, the application of MPS in mineral deposit modeling still faces many challenges. 相似文献
4.
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. 相似文献
5.
《Advances in water resources》2003,26(6):635-647
A kinetic flux vector splitting (KFVS) scheme for shallow water flows based on the collisionless Boltzmann equation is formulated and applied. The scheme is explicit and first order in space and time with stability governed by the Courant condition. The consistency of the KFVS scheme with the shallow water equations is proven using the equivalent differential equations approach. The accuracy and efficiency of the KFVS scheme in modeling complex flow features are compared to those of the Boltzmann Bhatnagar–Gross–Krook (BGK) scheme as well as a Riemann-based scheme. In particular, all schemes are applied to (i) strong shock waves, (ii) extreme expansion waves, (iii) a combination of strong shock waves and extreme expansion waves, and (iv) a one-dimensional dam break problem. Additionally, the KFVS, BGK and Riemann schemes are applied to a one-dimensional dam break problem for which laboratory data is available. These test cases reveal that all three schemes provide solutions of comparable accuracy, but the KFVS model is 1.5–2 times faster to execute than the BGK scheme and 2–3 times faster than the Riemann-based scheme. The absence of the collision term from the Boltzmann equation not only makes the mathematical formulation of KFVS easy but also helps elucidate this approach to the novice. The accuracy, efficiency, and simplicity of the KFVS scheme indicate its potential in modeling an array of water resources problems. Due to the scalar nature of the Boltzmann equation, the extension of the KFVS scheme to 2-D surface water flows is straightforward. 相似文献
6.
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. 相似文献
7.
A two-phase numerical model using Smoothed Particle Hydrodynamics (SPH) is applied to two-phase liquid-sediments flows. The absence of a mesh in SPH is ideal for interfacial and highly non-linear flows with changing fragmentation of the interface, mixing and resuspension. The rheology of sediment induced under rapid flows undergoes several states which are only partially described by previous research in SPH. This paper attempts to bridge the gap between the geotechnics, non-Newtonian and Newtonian flows by proposing a model that combines the yielding, shear and suspension layer which are needed to predict accurately the global erosion phenomena, from a hydrodynamics prospective. The numerical SPH scheme is based on the explicit treatment of both phases using Newtonian and the non-Newtonian Bingham-type Herschel-Bulkley-Papanastasiou constitutive model. This is supplemented by the Drucker-Prager yield criterion to predict the onset of yielding of the sediment surface and a concentration suspension model. The multi-phase model has been compared with experimental and 2-D reference numerical models for scour following a dry-bed dam break yielding satisfactory results and improvements over well-known SPH multi-phase models. With 3-D simulations requiring a large number of particles, the code is accelerated with a graphics processing unit (GPU) in the open-source DualSPHysics code. The implementation and optimisation of the code achieved a speed up of x58 over an optimised single thread serial code. A 3-D dam break over a non-cohesive erodible bed simulation with over 4 million particles yields close agreement with experimental scour and water surface profiles. 相似文献
8.
M. ALIPARAST 《国际泥沙研究》2009,24(1):99-107
A numerical model based upon a second-order upwind ceil-center f＇mite volume method on unstructured triangular grids is developed for solving shallow water equations. The assumption of a small depth downstream instead of a dry bed situation changes the wave structure and the propagation speed of the front which leads to incorrect results. The use of Harten-Lax-vau Leer （HLL） allows handling of wet/dry treatment. By usage of the HLL approximate Riemann solver, also it make possible to handle discontinuous solutions. As the assumption of a very small depth downstream oftbe dam can change the nature of the dam break flow problem which leads to incorrect results, the HLL approximate Riemann solver is used for the computation of inviscid flux functions, which makes it possible to handle discontinuous solutions. A multidimensional slope-limiting technique is applied to achieve second-order spatial accuracy and to prevent spurious oscillations. To alleviate the problems associated with numerical instabilities due to small water depths near a wet/dry boundary, the friction source terms are treated in a fully implicit way. A third-order Runge-Kutta method is used for the time integration of semi-discrete equations. The developed numerical model has been applied to several test cases as well as to real flows. The tests are tested in two cases： oblique hydraulic jump and experimental dam break in converging-diverging flume. Numerical tests proved the robustness and accuracy of the model. The model has been applied for simulation of dam break analysis of Torogh in Iran. And finally the results have been used in preparing EAP （Emergency Action Plan）. 相似文献
9.
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. 相似文献
10.
Mesh-free discretization, flexibly distributing nodes without computationally expensive meshing process, is able to deal with staircase problem, oversampling and undersampling problems and saves plenty of nodes through distributing nodes suitably with respect to irregular boundaries and model parameters. However, the time-domain mesh-free discretization usually exhibits poorer stability than that in regular grid discretization. In order to reach unconditional stability and easy implementation in parallel computing, we develop the frequency-domain finite-difference method in a mesh-free discretization, incorporated with two perfectly matched layer boundary conditions. Furthermore, to maintain the flexibility of mesh-free discretization, the nodes are still irregularly distributed in the absorbing zone, which complicates the situation of artificial boundary reflections. In this paper, we implement frequency-domain acoustic wave modelling in a mesh-free system. First, we present the perfectly matched layer boundary condition to suppress spurious reflections. Moreover, we develop the complex frequency shifted–perfectly matched layer boundary condition to improve the attenuation of grazing waves. In addition, we employ the radial-basis-function-generated finite difference method in the mesh-free discretization to calculate spatial derivatives. The numerical experiment on a rectangle homogeneous model shows the effectiveness of the perfectly matched layer boundary condition and the complex frequency shifted–perfectly matched layer boundary condition, and the latter one is better than the former one when absorbing large angle incident waves. The experiment on the Marmousi model suggests that the complex frequency shifted–perfectly matched layer boundary condition works well for complicated models. 相似文献
11.
Desorption-diffusion model and lost gas quantity estimation of coalbed methane from coal core under drilling fluid medium 总被引：1，自引：0，他引：1
Geoff WANG 《中国科学:地球科学(英文版)》2010,(4)
The differences of coalbed methane(CBM) desorption-diffusion from coal drilling-core under various drilling fluid medium are not considered in the present calculating methods of lost CBM quantity,which leads possibly to the inaccuracy of CBM quantity in coal seam.Here we took the desorption of CBM from coal core under drilling fluid medium as a pressure-swing process,and based on the Langmuir equation and Fick-first law,established the desorption-diffusion model and numerical modeling method of lost gas(inc... 相似文献
12.
A time‐domain seismic SSI analysis method for inelastic bridge structures through the use of a frequency‐dependent lumped parameter model 
下载免费PDF全文
下载免费PDF全文 It is highlighted in the past that the soil–structure interaction phenomenon can produce a significant alteration on the response of a bridge structure. A variety of approaches has been developed in the past, which is capable of tackling the soil–structure interaction problem from different perspectives. The popular approach of a discretized truncated finite element model of the soil domain is not always a numerically viable solution, especially for computationally demanding simulations such as the probabilistic fragility analysis of a bridge structure or the real time hybrid simulation. This paper aims to develop a complete modeling procedure that is capable of coping with the soil–structure interaction problem of inelastic bridge structures through the use of a frequency dependent lumped parameter assembly. The proposed procedure encounters accuracy and global stability issues observed on past methods while maintaining the broad applicability of the method by any commercial FEM software. A case study of an overpass bridge structure under earthquake excitations is illustrated in order to verify the proposed method. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
13.
Large-scale shaking table model tests on seismically induced failure of Concrete-Faced Rockfill Dams
An increasing number of high Concrete-Faced Rockfill Dams (CFRDs) have been and are being built in highly seismic hazardous regions. Because there are few examples, the failure mechanisms for each damaged CFRD still remain unclear, which prevents a rational evaluation of aseismic performance of a CFRD. In this paper, a series of large-scale shaking table tests are conducted to determine the dynamic failure modes of CFRDs, especially the failure mechanism of the face slab. A type of model slab material is developed according to similitude rules. Particle Image Velocimetry (PIV) is used and improved to track the motion of each grain and subsequently measure the entire velocity field of the deforming cross-sections of dam models. The failure processes of the dams, the failure mechanisms of the slab fractures and the slab dislocations are discussed based on the experimental results. The experiments indicate that large deformations of the uppermost part of the dam caused the face slab to curve upward and crack. Another important fracture morphology of the face slab is the dislocation of its joint. The change in the granular fabric caused the granular slope to slide down and thrust outward, eventually dragging the upper half of the face slab down against its lower half; this motion is accompanied by the formation of a joint dislocation. Therefore, the stability and anti-deformation ability of both the downstream slope and the crest region of the upstream slope should be improved to withstand earthquake. The results of the shaking table model tests are consistent with the numerical simulation, and both results support the prototype dam׳s seismically induced failure characteristics. 相似文献
14.
Uplift mechanism for a shallow-buried structure in liquefiable sand subjected to seismic load: centrifuge model test and DEM modeling 总被引：2，自引：1，他引：1
Based on a centrifuge model test and distinct element method(DEM), this study provides new insights into the uplift response of a shallow-buried structure and the liquefaction mechanism for saturated sand around the structure under seismic action. In the centrifuge test, a high-speed microscopic camera was installed in the structure model, by which the movements of particles around the structure were monitored. Then, a two-dimensional digital image processing technology was used to analyze the microstructure of saturated sand during the shaking event. Herein, a numerical simulation of the centrifuge experiment was conducted using a two-phase(solid and fl uid) fully coupled distinct element code. This code incorporates a particle-fl uid coupling model by means of a "fi xed coarse-grid" fl uid scheme in PFC3D(Particle Flow Code in Three Dimensions), with the modeling parameters partially calibrated based on earlier studies. The physical and numerical models both indicate the uplifts of the shallow-buried structure and the sharp rise in excess pore pressure. The corresponding micro-scale responses and explanations are provided. Overall, the uplift response of an underground structure and the occurrence of liquefaction in saturated sand are predicted successfully by DEM modeling. However, the dynamic responses during the shaking cannot be modeled accurately due to the restricted computer power. 相似文献
15.
16.
Fatih Ünes 《洁净——土壤、空气、水》2010,38(3):296-308
Experimental findings and observations indicate that plunging flow is related to the formation of bed load deposition in dam reservoirs. The sediment delta begins to form in the plunging region where the inflow river water meets the ambient reservoir water. Correct estimation of dam reservoir flow, plunging point, and plunging depth is crucial for dam reservoir sedimentation and water quality issues. In this study, artificial neural network (ANN), multi‐linear regression (MLR), and two‐dimensional hydrodynamic model approaches are used for modeling the plunging point and depth. A multi layer perceptron (MLP) is used as the ANN structure. A two‐dimensional model is adapted to simulate density plunging flow through a reservoir with a sloping bottom. In the model, nonlinear and unsteady continuity, momentum, energy, and k–ε turbulence equations are formulated in the Cartesian coordinates. Density flow parameters such as velocity, plunging points, and plunging depths are determined from the simulation and model results, and these are compared with previous experimental and model works. The results show that the ANN model forecasts are much closer to the experimental data than the MLR and mathematical model forecasts. 相似文献
17.
陈江 《地震工程与工程振动》2010,30(2)
基于接触非线性有限元模型,以锦屏一级拱坝为例,库水分别采用附加质量模型、可压缩流体有限元模型、不可压缩流体有限元模型计算了正常蓄水位及运行低水位时坝体的动力响应,结果表明:库水模型对拱坝动力响应有较大影响,随库水深度的增大,各模型计算结果差异增大;相比于流体可压缩模型,采用不可压缩流体模型所得动力响应普遍偏大;运行低水位工况,由于静水压力减小导致拱效应减弱,从而降低了拱坝的整体性,因此运行低水位工况各缝开度普遍高于正常蓄水位工况,且其拉应力范围较大,因此,运行低水位工况将对抗震设计起控制作用。 相似文献
18.
1 INTRODUCTION 1. 1 Dam Removal The aging of the more than 75,000 dams in the U.S., coupled with the increasing awareness of their environmental costs, has made dam decommissioning and removal a topic of current interest to the scientific community, management agencies, and the general public. It is estimated that 85% of the dams in the U.S. will be near the end of their operational lives by the year 2020 (FEMA, 1999), necessitating thorough consideration of dam removal or repair for… 相似文献
19.
将时空守恒元/解元(CE/SE)方法推广到二维孔隙介质多相流问题的数值计算中,采用人工压缩法耦合速度和压力,同时结合杂交粒子水平集方法捕捉物质界面.提出一套完整的二维欧拉型孔隙介质非稳态多相不可压缩黏性流动计算方案.通过对溃坝和液滴在重力作用下的运动和变形问题的数值模拟,验证了方法的精度和有效性.在此基础上,提出了一个新的孔隙介质两相流物理模型——双层流体顶盖驱动方腔流. 相似文献
20.
A semi-analytical mesh-free series solution method is presented for modeling regional steady-state subsurface saturated–unsaturated flow in 2-D geometrically complex homogenous and stratified hill-slope cross sections. Continuous solutions for pressure in the saturated and unsaturated zone are determined iteratively, as is the location of the water table surface. Mass balance is satisfied exactly over the entire domain except along boundaries and interfaces between layers, where errors are in an acceptable range. The solutions are derived and demonstrated on multiple test cases. The errors for specific cases are assessed and discussed. 相似文献