共查询到20条相似文献，搜索用时 62 毫秒
1.
Within the framework of the Godunov-type cell-centered finite volume (CCFV) scheme, this paper proposes a 2D well-balanced shallow water model for unstructured grids. In this model, the face-based van Albada limiting scheme is employed in conjunction with a directional correction to reconstruct second order spatial values at the midpoint of the considered face. The Harten, Lax and van Leer approximate Riemann solver with the Contact wave restored (HLLC) is applied to compute the fluxes of mass and momentum, while the splitting implicit method is utilized to solve the friction source terms. The novel aspects of the model include the new limited directional correction with which the new local extrema caused by the unlimited correction are prevented efficiently, the simplified non-negative water depth reconstruction used to get rid of numerical instabilities and in turn to preserve mass conservation at wet–dry interfaces and the novel slope source term treatment which suits complex unstructured grids well by transforming the slope source of a cell into fluxes at its faces. This model is able to preserve the C-property and mass conservation, to achieve good convergence to steady state, to capture discontinuous flows and to handle complex flows involving wetting and drying over uneven beds on unstructured grids with poor connectivity in an accurate, efficient and robust way. These capabilities are verified against analytical solutions, numerical results of alternative models and experimental and field data. 相似文献
2.
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. 相似文献
3.
A two-dimensional (2D) numerical model has been developed to solve shallow water equations for simulation of dam-break flows. The spatial derivatives are discretized using a well-balanced explicit central upwind conservative scheme. The scheme is Riemann solver free and guarantees the positivity of the flow depth over complex topography if the Courant number is kept less than 0.25. The time integration is performed by Euler’s scheme. The model is verified against analytical results for water surface elevation and discharge for three benchmark test cases. A good agreement between analytical solutions and computed results is observed. The property of well-balancing in still water over an uneven bottom is also confirmed. The model is then validated by simulating a laboratory experiment in which a dam break flow propagates over a triangular obstacle. The model performance was found to be satisfactory. A dam break laboratory experimental test case on a frictionless horizontal bottom is also simulated for 2D validation of the model, and good agreement between simulation and the experimental data is observed. The suitability of the proposed model for real life applications is demonstrated by simulating the Malpasset dam-break event, which occurred in 1959 in France. The computed arrival time of the flood wave front and the maximum flow depths at various observation points matched well with the measurements on a 1/400 scale physical model. The overall performance indicates that this model can be applied for simulation of dam-break waves in real life cases. 相似文献
4.
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. 相似文献
5.
An important part in the numerical simulation of tsunami and storm surge events is the accurate modeling of flooding and the appearance of dry areas when the water recedes. This paper proposes a new algorithm to model inundation events with piecewise linear Runge–Kutta discontinuous Galerkin approximations applied to the shallow water equations. This study is restricted to the one-dimensional case and shows a detailed analysis and the corresponding numerical treatment of the inundation problem.The main feature is a velocity based “limiting” of the momentum distribution in each cell, which prevents instabilities in case of wetting or drying situations. Additional limiting of the fluid depth ensures its positivity while preserving local mass conservation. A special flux modification in cells located at the wet/dry interface leads to a well-balanced method, which maintains the steady state at rest. The discontinuous Galerkin scheme is formulated in a nodal form using a Lagrange basis. The proposed wetting and drying treatment is verified with several numerical simulations. These test cases demonstrate the well-balancing property of the method and its stability in case of rapid transition of the wet/dry interface. We also verify the conservation of mass and investigate the convergence characteristics of the scheme. 相似文献
6.
The shallow water equations are used to model flows in rivers and coastal areas, and have wide applications in ocean, hydraulic engineering, and atmospheric modeling. These equations have still water steady state solutions in which the flux gradients are balanced by the source term. It is desirable to develop numerical methods which preserve exactly these steady state solutions. Another main difficulty usually arising from the simulation of dam breaks and flood waves flows is the appearance of dry areas where no water is present. If no special attention is paid, standard numerical methods may fail near dry/wet front and produce non-physical negative water height. A high-order accurate finite volume weighted essentially non-oscillatory (WENO) scheme is proposed in this paper to address these difficulties and to provide an efficient and robust method for solving the shallow water equations. A simple, easy-to-implement positivity-preserving limiter is introduced. One- and two-dimensional numerical examples are provided to verify the positivity-preserving property, well-balanced property, high-order accuracy, and good resolution for smooth and discontinuous solutions. 相似文献
7.
Dam-break flows usually propagate along rivers and floodplains, where the processes of fluid flow, sediment transport and bed evolution are closely linked. However, the majority of existing two-dimensional (2D) models used to simulate dam-break flows are only applicable to fixed beds. Details are given in this paper of the development of a 2D morphodynamic model for predicting dam-break flows over mobile beds. In this model, the common 2D shallow water equations are modified, so that the effects of sediment concentrations and bed evolution on the flood wave propagation can be considered. These equations are used together with the non-equilibrium transport equations for graded sediments and the equation of bed evolution. The governing equations are solved using a matrix method, thus the hydrodynamic, sediment transport and morphological processes can be jointly solved. The model employs an unstructured finite volume algorithm, with an approximate Riemann solver, based on the Roe-MUSCL scheme. A predictor–corrector scheme is used in time stepping, leading to a second-order accurate solution in both time and space. In addition, the model considers the adjustment process of bed material composition during the morphological evolution process. The model was first verified against results from existing numerical models and laboratory experiments. It was then used to simulate dam-break flows over a fixed bed and a mobile bed to examine the differences in the predicted flood wave speed and depth. The effects of bed material size distributions on the flood flow and bed evolution were also investigated. The results indicate that there is a great difference between the dam-break flow predictions made over a fixed bed and a mobile bed. At the initial stage of a dam-break flow, the rate of bed evolution could be comparable to that of water depth change. Therefore, it is often necessary to employ the turbid water governing equations using a coupled approach for simulating dam-break flows. 相似文献
8.
Numerical resolution of well-balanced shallow water equations with complex source terms 总被引：3，自引：0，他引：3
This paper presents a well-balanced numerical scheme for simulating frictional shallow flows over complex domains involving wetting and drying. The proposed scheme solves, in a finite volume Godunov-type framework, a set of pre-balanced shallow water equations derived by considering pressure balancing. Non-negative reconstruction of Riemann states and compatible discretization of slope source term produce stable and well-balanced solutions to shallow flow hydrodynamics over complex topography. The friction source term is discretized using a splitting implicit scheme. Limiting value of the friction force is derived to ensure stability. This new numerical scheme is validated against four theoretical benchmark tests and then applied to reproduce a laboratory dam break over a domain with irregular bed profile. 相似文献
9.
A shallow flow generally features complex hydrodynamics induced by complicated domain topography and geometry. A numerical scheme with well-balanced flux and source term gradients is therefore essential before a shallow flow model can be applied to simulate real-world problems. The issue of source term balancing has been exhaustively investigated in grid-based numerical approaches, e.g. discontinuous Galerkin finite element methods and finite volume Godunov-type methods. In recent years, a relatively new computational method, smooth particle hydrodynamics (SPH), has started to gain popularity in solving the shallow water equations (SWEs). However, the well-balanced problem has not been fully investigated and resolved in the context of SPH. This work aims to discuss the well-balanced problem caused by a standard SPH discretization to the SWEs with slope source terms and derive a corrected SPH algorithm that is able to preserve the solution of lake at rest. In order to enhance the shock capturing capability of the resulting SPH model, the Monotone Upwind-centered Scheme for Conservation Laws (MUSCL) is also explored and applied to enable Riemann solver based artificial viscosity. The new SPH model is validated against several idealized benchmark tests and a real-world dam-break case and promising results are obtained. 相似文献
10.
W. Fellin 《Pure and Applied Geophysics》2002,159(7-8):1737-1748
— Godunov's method, a numerical method for solving conservation laws, is applied to nonlinear and inelastic wave propagation in soil. The solution is restricted to the one-dimensional case. An approximate Riemann solver for Godunov's method is presented. The capability of the numerical method is shown by a comparison with the analytical solution of a linear inelastic wave propagation. Finally the behaviour of the nonlinear inelastic soil is described by a hypoplastic constitutive law. 相似文献
11.
We propose a novel computational method for the efficient simulation of two-phase flow in fractured porous media. Instead of refining the grid to capture the flow along the faults or fractures, we represent the latter as immersed interfaces, using a reduced model for the flow and suitable coupling conditions. We allow for non matching grids between the porous matrix and the fractures to increase the flexibility of the method in realistic cases. We employ the extended finite element method for the Darcy problem and a finite volume method that is able to handle cut cells and matrix-fracture interactions for the saturation equation. Moreover, we address through numerical experiments the problem of the choice of a suitable numerical flux in the case of a discontinuous flux function at the interface between the fracture and the porous matrix. A wrong approximate solution of the Riemann problem can yield unphysical solutions even in simple cases. 相似文献
12.
Seokkoo KangAnne Lightbody Craig HillFotis Sotiropoulos 《Advances in water resources》2011,34(1):98-113
We develop an efficient and versatile numerical model for carrying out high-resolution simulations of turbulent flows in natural meandering streams with arbitrarily complex bathymetry. The numerical model solves the 3D, unsteady, incompressible Navier-Stokes and continuity equations in generalized curvilinear coordinates. The method can handle the arbitrary geometrical complexity of natural streams using the sharp-interface curvilinear immersed boundary (CURVIB) method of Ge and Sotiropoulos (2007) [1]. The governing equations are discretized with three-point, central, second-order accurate finite-difference formulas and integrated in time using an efficient, second-order accurate fractional step method. To enable efficient simulations on grids with tens of millions of grid nodes in long and shallow domains typical of natural streams, the algebraic multigrid (AMG) method is used to solve the Poisson equation for the pressure coupled with a matrix-free Krylov solver for the momentum equations. Depending on the desired level of resolution and available computational resources, the numerical model can either simulate, via direct numerical simulation (DNS), large-eddy simulation (LES), or unsteady Reynolds-averaged Navier-Stokes (URANS) modeling. The potential of the model as a powerful tool for simulating energetic coherent structures in turbulent flows in natural river reaches is demonstrated by applying it to carry out LES and URANS in a 50-m long natural meandering stream at resolution sufficiently fine to capture vortex shedding from centimeter-scale roughness elements on the bed. The accuracy of the simulations is demonstrated by comparisons with experimental data and the relative performance of the LES and URANS models is also discussed. 相似文献
13.
The sensitivity of a model output (called a variable) to a parameter can be defined as the partial derivative of the variable with respect to the parameter. When the governing equations are not differentiable with respect to this parameter, problems arise in the numerical solution of the sensitivity equations, such as locally infinite values or instability. An approximate Riemann solver is thus proposed for direct sensitivity calculation for hyperbolic systems of conservation laws in the presence of discontinuous solutions. The proposed approach uses an extra source term in the form of a Dirac function to restore sensitivity balance across the shocks. It is valid for systems such as the Euler equations for gas dynamics or the shallow water equations for free surface flow. The method is first detailed and its application to the shallow water equations is proposed, with some test cases such as dike- or dam-break problems with or without source terms. An application to a two-dimensional flow problem illustrates the superiority of direct sensitivity calculation over the classical empirical approach. 相似文献
14.
１ＩＮＴＲＯＤＵＣＴＩＯＮＲｉｖｅｒｓｉｎＴａｉｗａｎａｒｅｒｅｌａｔｉｖｅｌｙｓｔｅｅｐｅｒｃｏｍｐａｒｅｄｔｏｔｈｏｓｅｉｎｏｔｈｅｒｃｏｎｔｉｎｅｎｔ．Ｌｏｃａｌｙｏｃｕｒｅｄｓｕｐｅｒｃｒｉｔｉｃａｌｆｌｏｗａｒｅｆａｉｒｌｙｃｏｍｍｏｎｉｎ... 相似文献
15.
Despite the intensive research over the past decades in the field of stochastic subsurface hydrology, our ability to analyze and model heterogeneous groundwater systems remains limited. Most existing theories are either too restrictive to handle practical complexity or too expensive to be applied to realistic problem sizes. In this paper we present approximate, closed-form equations that allow modeling 2D nonstationary flows in statistically inhomogeneous aquifers, including composite aquifers containing multiple zones characterized by different statistical models. The composite representation has the effect of decreasing the variance of deviations from the mean, relaxing the limitation of the small-perturbation assumption. The simple formulas are illustrated with a number of examples and compared with a corresponding first-order nonstationary numerical analysis and Monte Carlo simulation. The results show that, despite the gross simplifications, the closed-form equations are robust and able to capture complex variance dynamics, reproducing surprisingly well the first-order numerical solutions and the Monte Carlo simulation even in highly nonstationary, variable situations. 相似文献
16.
Cheng-Chin Wu 《地球物理与天体物理流体动力学》2013,107(1):37-61
We present a high order accurate weighted essentially non-oscillatory (WENO) finite difference scheme for solving the equations of incompressible fluid dynamics and magnetohydrodynamics (MHD). This scheme is a direct extension of a WENO scheme that has been successfully applied to compressible fluids, with or without magnetic fields. A fractional time-step method is used to enforce the incompressibility condition. Two basic elements of the WENO scheme, upwinding and wave decomposition, are shown to be important in solving the incompressible systems. Numerical results demonstrate that the scheme performs well for one-dimensional Riemann problems, a two-dimensional double-shear flow problem, and the two-dimensional Orszag–Tang MHD vortex system. They establish that the WENO code is numerical stable even when there are no explicit dissipation terms. It can handle discontinuous data and attain converged results with a high order of accuracy. 相似文献
17.
During the last four decades, several numerical formulations and specialized software have been developed in response to studies about dam break (DB) wave propagation and its hydraulic and environmental impacts on downstream hydraulic structures and valleys. These methods cannot, however, be used to predict wave propagation within partially covered channels or reservoirs located upstream of hydraulic structures. In fact, such problems require the modelling of the complex transition from a free surface flow into a pressurized one. Because rivers or channels partially covered with ice sheets are typical examples commonly met in winter in northern climates, it is vitally important to assess ice-cover effects on the DB wave propagation and develop appropriate tools to predict resulting hydrodynamic loads on hydraulic structures downstream. This paper proposes an original numerical formulation to model wave propagation and hydrodynamic pressure in partially covered channels. The proposed formulation uses one-dimensional St. Venant equations to simulate open-water flow and water hammer equations to simulate pressure flow within the partially covered channel. To illustrate the use of the hydrodynamic pressures obtained, a case study is presented where a channel cover and a dam located downstream are modelled using finite elements to investigate their dynamic structural response. 相似文献
18.
Infiltration systems are widely used as an effective urban stormwater control measure. Most design methods and models roughly approximate the complex physical flow processes in these systems using empirical equations and fixed infiltration rates to calculate emptying times from full. Sophisticated variably saturated flow models are available, but rarely applied owing to their complexity. This paper describes the development and testing of an integrated one‐dimensional model of flow through the porous storage of a typical infiltration system and surrounding soils. The model accounts for the depth in the storage, surrounding soil moisture conditions and the interaction between the storage and surrounding soil. It is a front‐tracking model that innovatively combines a soil‐moisture‐based solution of Richard's equation for unsaturated flow with piston flow through a saturated zone as well as a reservoir equation for flow through a porous storage. This allows the use of a simple non‐iterative numerical solution that can handle ponded infiltration into dry soils. The model is more rigorous than approximate stormwater infiltration system models and could therefore be valuable in everyday practice. A range of test cases commonly used to test soil water flow models for infiltration in unsaturated conditions, drainage from saturation and infiltration under ponded conditions were used to test the model along with an experiment with variable depth in a porous storage over saturated conditions. Results show that the model produces a good fit to the observed data, analytical solutions and Hydrus. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
19.
Shallow sediment transport flow computation using time-varying sediment adaptation length 总被引：1，自引：0，他引：1
Based on the common approach,the adaptation length in sediment transport is normally estimated astemporally independent.However,this approach might not be theoretically justified as the process of reaching the sediment transport equilibrium stage is affected by the flow conditions in time,especially for fast moving flows,such as scour-hole developing flows.In this study,the two-dimensional（2D） shallow water formulation together with a sediment continuity-concentration（SCC） model were applied to flow with mobile sediment boundary.A timevarying approach was proposed to determine the sediment transport adaptation length to simulate the sediment erosion-deposition rate.The proposed computational model was based on the Finite Volume（FV） method.The Monotone Upwind Scheme of Conservative Laws（MUSCL）-Hancock scheme was used with the Harten Lax van Leer-contact（HLLC） approximate Riemann solver to discretize the FV model.In the flow applications of this paper,a highly discontinuous dam-break,fast sediment transport flow was used to calibrate the proposed timevarying sediment adaptation length model.Then the calibrated model was further applied to two separate experimental sediment transport flow applications documented in the literature,i.e.a highly concentrated sediment transport flow in a wide alluvial channel and a sediment aggradation flow.Good agreement with the experimental data were obtained with the proposed model simulations.The tests prove that the proposed model,which was calibrated by the discontinuous dam-break bed scouring flow,also performed well to represent rapid bed change and steady sediment mobility conditions. 相似文献
20.
A finite difference implicit scheme is presented in this paper for solution of the shallow water equations in one dimensional (1D) form. The present model has many advantages like, handling of discontinuous and complex bed topography, satisfying C-property (preservation of motionless water surface over a wet or dry bed) and capability of handling large value of temporal step etc. Another very important feature of the present model is that, no special treatment of the source vector of the governing equations is required here to deal with very less water depth. To investigate the performance of the present model in diverse situations, it is used to replicate four different problems of known analytical solution, and the model is found to be quite capable for varied situations. 相似文献