共查询到20条相似文献，搜索用时 15 毫秒
1.
Peter Mewis 《国际流体数值方法杂志》2013,72(8):864-882
Within the mixed FEM, the mini‐element that uses a bubble shape function for the solution of the shallow water wave equations on triangle meshes is simplified to a sparse element formulation. The new formulation has linear shape functions for water levels and constant shape functions for velocities inside each element. The suppression of decoupled spurious solutions is excellent with the new scheme. The linear dispersion relation of the new element has similar advantages as that of the wave equation scheme (generalised wave continuity scheme) proposed by Lynch and Gray. It is shown that the relation is monotonic over all wave numbers. In this paper, the time stepping scheme is included in the dispersion analysis. In case of a combined space–time staggering, the dispersion relation can be improved for the shortest waves. The sparse element is applied in the flow model Bubble that conserves mass exactly. At the same time, because of the limited number of degrees of freedom, the computational efficiency is high. The scheme is not restricted to orthogonal triangular meshes. Three test cases demonstrate the very good accuracy of the proposed scheme. The examples are the classical quarter annulus test case for the linearised shallow water equations, the hydraulic jump and the tide in the Elbe river mouth. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
2.
Matthew Horritt 《国际流体数值方法杂志》2004,44(11):1231-1255
Many environmental applications of shallow water flow modelling can be characterized as only slowly varying and everywhere sub‐critical. A simplified finite volume model is therefore developed that is capable of describing pertinent shallow water flow processes more efficiently than the usual Godunov/ Riemann characteristics approaches. The model is tested against a number of analytical and numerical solutions to the governing equations. The model reproduces accurately flow round a circular bend, flow over topography, flow up an initially dry beach and floodwave propagation down a meandering river reach, with mass conservative solutions. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
3.
4.
Michele La Rocca Claudia Adduce Valentina Lombardi Giampiero Sciortino Reinhard Hinkelmann 《国际流体数值方法杂志》2012,70(8):1048-1072
In this paper the dynamics of a two‐layered liquid, made of two immiscible shallow‐layers of different density, has been investigated within the framework of the lattice Boltzmann method (LBM). The LBM developed in this paper for the two‐layered, shallow‐water flow has been obtained considering two separate sets of LBM equations, one for each layer. The coupling terms between the two sets have been defined as external forces, acted on one layer by the other. Results obtained from the LBM developed in this paper are compared with numerical results obtained solving the two‐layered, shallow‐water equations, with experimental and other numerical results published in literature. The results are interesting. First, the numerical results obtained by the LBM and by the shallow‐water model can be considered as equivalent. Second, the LBM developed in this paper is able to simulate motion conditions on nonflat topography. Third, the agreement between the LBM (and also shallow‐water model) numerical results and the experimental results is good when the evolution of the flow does not depend on the viscosity, that is, during the initial phase of the flow, dominated by gravity and inertia forces. When the viscous forces dominate the evolution of the flow the agreement between numerical and experimental results depends strongly on the viscosity; it is good if the numerical LBM viscosity has the same order of magnitude of the liquid's kinematic viscosity. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
5.
Sébastien Erpicum Tobias Meile Benjamin J. Dewals Michel Pirotton Anton J. Schleiss 《国际流体数值方法杂志》2009,61(11):1227-1246
A 2D numerical flow model, developed at the Research unit of Hydrology, Applied Hydrodynamics and Hydraulic Constructions at ULg, has been applied to flows in a macro‐rough channel. The model solves the shallow water equations (SWE) with a two length scale, depth‐integrated k‐type approach for turbulence modeling. Data for the comparison have been provided by experiments conducted at the Laboratory of Hydraulic Constructions at EPFL. In the experiments with different non‐prismatic channel configurations, namely large‐scale cavities at the side walls, three different 2D flow characteristics could be observed in cavities. With the used numerical model features, especially regarding turbulence and friction modeling, a single set of bottom and side wall roughness could be found for a large range of discharges investigated in a prismatic channel. For the macro rough configurations, the numerical model gives an excellent agreement between experimental and numerical results regarding backwater curves and flow patterns if the side wall cavities have low aspect ratios. For configurations with high aspect ratios, the head loss generated by the preservation of important recirculation gyres in the cavities is slightly underestimated. The results of the computations reveal clearly that the separation of turbulence sources in the mathematical model is of great importance. Indeed, the turbulence related to 2D transverse shear effects and the 3D turbulence, generated by bed friction, can have very different amplitude. When separating these two effects in the numerical models, most of the flow features observed experimentally can be reproduced accurately. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
6.
We present a remeshed particle‐mesh method for the simulation of three‐dimensional compressible turbulent flow. The method is related to the meshfree smoothed particle hydrodynamics method, but the present method introduces a mesh for efficient calculation of the pressure gradient, and laminar and turbulent diffusion. In addition, the mesh is used to remesh (reorganise uniformly) the particles to ensure a regular particle distribution and convergence of the method. The accuracy of the presented methodology is tested for a number of benchmark problems involving two‐ and three‐dimensional Taylor‐Green flow, thin double shear layer, and three‐dimensional isotropic turbulence. Two models were implemented, direct numerical simulations, and Smagorinsky model. Taking advantage of the Lagrangian advection, and the finite difference efficiency, the method is capable of providing quality simulations while maintaining its robustness and versatility. 相似文献
7.
A mesh‐free particle method, based on the moving particle semi‐implicit (MPS) interaction model, has been developed for the simulation of two‐dimensional open‐boundary free‐surface flows. The incompressibility model in the original MPS has been replaced with a weakly incompressible model. The effect of this replacement on the efficiency and accuracy of the model has been investigated. The new inflow–outflow boundary conditions along with the particle recycling strategy proposed in this study extend the application of the model to open‐boundary problems. The final model is able to simulate open‐boundary free surface flow in cases of large deformation and fragmentation of free surface. The models and proposed algorithms have been validated and applied to sample problems. The results confirm the model's efficiency and accuracy. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
8.
In this paper, the flow/acoustics splitting method for predicting flow‐generated noise is further developed by introducing high‐order finite difference schemes. The splitting method consists of dividing the acoustic problem into a viscous incompressible flow part and an inviscid acoustic part. The incompressible flow equations are solved by a second‐order finite volume code EllipSys2D/3D. The acoustic field is obtained by solving a set of acoustic perturbation equations forced by flow quantities. The incompressible pressure and velocity form the input to the acoustic equations. The present work is an extension of our acoustics solver, with the introduction of high‐order schemes for spatial discretization and a Runge–Kutta scheme for time integration. To achieve low dissipation and dispersion errors, either Dispersion‐Relation‐Preserving (DRP) schemes or optimized compact finite difference schemes are used for the spatial discretizations. Applications and validations of the new acoustics solver are presented for benchmark aeroacoustic problems and for flow over an NACA 0012 airfoil. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
9.
When performing shallow flow simulations on adaptive grids, the C‐property (i.e. conservation property) and the mass conservation may not be simultaneously preserved, that is, either C‐property or mass conservation is likely to be violated following grid refining or coarsening. The cause of such a contradiction is analyzed in detail in this work, which essentially links to the reconstruction of bed and flow information in those newly created cells during grid adaptation. An effective approach is subsequently proposed to resolve the contradiction by locally modifying the bed elevation in certain problematic cells when reconstructing flow information using linear interpolation as part of the grid adaptation procedure. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
10.
This paper presents a first‐order HLLC (Harten‐Lax‐Van Leer with contact discontinuities) scheme to solve the Saint‐Venant shallow‐water equations, including morphological evolution of the bed by erosion and deposition of sediments. The Exner equation is used to model the morphological evolution of the bed, while a closure equation is needed to evaluate the rate of sediment transport. The system of Saint‐Venant–Exner equations is solved in a fully coupled way using a finite‐volume technique and a HLLC solver for the fluxes, with a novel wave‐speed estimator adapted to the Exner equation. Wave speeds are usually derived by computing the eigenvalues of the full system, which is highly time‐consuming when no analytical expression is available. In this paper, an eigenvalue analysis of the full system is conducted, leading to simple but still accurate wave‐speed estimators. The new numerical scheme is then tested in three different situations: (1) a circular dam‐break flow over movable bed, (2) an one‐dimensional bed aggradation problem simulated on a 2D unstructured mesh and (3) the case of a dam‐break flow in an erodible channel with a sudden enlargement, for which experimental measurements are available. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
11.
This paper presents the results of some studies on the development and application of a finite element method (FEM) with a closed-form solution technique for time discretization. The closed-form solution is based on the eigenvalues/vectors of a coefficient matrix. The method is first applied to the one-dimensional linearized shallow water equations and then extended to the two-dimensional shallow water equations. An attempt is made to improve its efficiency by incorporating time splitting and using the closed-form solution technique only for linear terms. Some case studies of a rectangular channel and harbour are presented to illustrate the satisfactory working of the method. © 1997 by John Wiley & Sons, Ltd. Int. j. numer. methods fluids 24: 953–963, 1997. 相似文献
12.
A new numerical scheme, namely space–time conservation element and solution element (CE/SE) method, has been used for the solution of the two‐dimensional (2D) dam‐break problem. Distinguishing from the well‐established traditional numerical methods (such as characteristics, finite difference, finite element, and finite‐volume methods), the CE/SE scheme has many non‐traditional features in both concept and methodology: space and time are treated in a unified way, which is the most important characteristic for the CE/SE method; the CEs and SEs are introduced, both local and global flux conservations in space and time rather than space only are enforced; an explicit scheme with a stagger grid is adopted. Furthermore, this scheme is robust and easy to implement. In this paper, an improved CE/SE scheme is extended to solve the 2D shallow water equations with the source terms, which usually plays a critical role in dam‐break flows. To demonstrate the accuracy, robustness and efficiency of the improved CE/SE method, both 1D and 2D dam‐break problems are simulated numerically, and the results are consistent with either the analytical solutions or experimental results. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
13.
The successful implementation of a finite element model for computing shallow‐water flow requires the identification and spatial discretization of a surface water region. Since no robust criterion or node spacing routine exists, which incorporates physical characteristics and subsequent responses into the mesh generation process, modelers are left to rely on crude gridding criteria as well as their knowledge of particular domains and their intuition. Two separate methods to generate a finite element mesh are compared for the Gulf of Mexico. A wavelength‐based criterion and an alternative approach, which employs a localized truncation error analysis (LTEA), are presented. Both meshes have roughly the same number of nodes, although the distribution of these nodes is very different. Two‐dimensional depth‐averaged simulations of flow using a linearized form of the generalized wave continuity equation and momentum equations are performed with the LTEA‐based mesh and the wavelength‐to‐gridsize ratio mesh. All simulations are forced with a single tidal constituent, M2. Use of the LTEA‐based procedure is shown to produce a superior (i.e., less error) two‐dimensional grid because the physics of shallow‐water flow, as represented by discrete equations, are incorporated into the mesh generation process. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
14.
Numerical modelling of shallow water flow in two dimensions is presented in this work with the results obtained in dam break tests. Free surface flow in channels can be described mathematically by the shallow‐water system of equations. These equations have been discretized using an approach based on unstructured Delaunay triangles and applied to the simulation of two‐dimensional dam break flows. A cell centred finite volume method based on Roe's approximate Riemann solver across the edges of the cells is presented and the results are compared for first‐ and second‐order accuracy. Special treatment of the friction term has been adopted and will be described. The scheme is capable of handling complex flow domains as shown in the simulation corresponding to the test cases proposed, i.e. that of a dam break wave propagating into a 45° bend channel (UCL) and in a channel with a constriction (LNEC‐IST). Comparisons of experimental and numerical results are shown. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
15.
16.
A novel efficient interface‐tracking method is developed to gain an insight into the interface in a multiphase or multifluid system, called the modified particle binary level set (MPBLS) method, in which the binary level set function is defined to distinguish the different phases or fluids and further modified by Lagrangian particles scattered along the interface for achieving higher accuracy. The validation of the MPBLS method is carried out first by simulating the free motion of a red blood cell (RBC) in the rotating, shear and Poiseuille flows, respectively. Subsequently, further validations are performed by comparing with the experimental and numerical results published previously. As one of important applications, the MPBLS method is employed to investigate the deformation behaviors of RBCs with different shapes in a capillary. The simulations show that the healthy RBC gradually changes the geometric shape from a biconcave to a steady parachute shape. It is thus guaranteed that the RBC successfully traverses through the smaller capillaries compared with undeformed RBC. However, the unhealthy RBC with the circular or elliptical shape has different deformation behaviors, in which the steady parachute shape is much less concave at the rear and more convex in the front. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
17.
Alberto Canestrelli Aukje Spruyt Bert Jagers Rudy Slingerland Mart Borsboom 《国际流体数值方法杂志》2016,81(3):151-177
In this work, an approach is proposed for solving the 3D shallow water equations with embedded boundaries that are not aligned with the underlying horizontal Cartesian grid. A hybrid cut‐cell/ghost‐cell method is used together with a direction‐splitting implicit solver: Ghost cells are used for the momentum equations in order to prescribe the correct boundary condition at the immersed boundary, while cut cells are used in the continuity equation in order to conserve mass. The resulting scheme is robust, does not suffer any time step limitation for small cut cells, and conserves fluid mass up to machine precision. Moreover, the solver displays a second‐order spatial accuracy, both globally and locally. Comparisons with analytical solutions and reference numerical solutions on curvilinear grids confirm the quality of the method. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
18.
The classical Saint‐Venant system is well suited for the modeling of dam breaks, hydraulic jumps, reservoir emptying, flooding etc. For many applications, the extension of the Saint‐Venant system to the case of non‐rectangular channels is necessary and this section‐averaged Saint‐Venant system exhibits additional source terms. The main difficulty of these equations consists of the discretization of these source terms. In this paper we propose a kinetic interpretation for the section averaged Saint‐Venant system and derive an associated numerical scheme. The numerical scheme—2nd order in space and time—preserves the positivity of the water height, and is well‐balanced. Numerical results including comparisons with analytic and experimental test problems illustrate the accuracy and the robustness of the numerical algorithm. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
19.
This study describes the influence of mass transfer on the steady two‐dimensional magnetohydrodynamic boundary layer flow of a Jeffery fluid bounded by a stretching sheet. A uniform magnetic field in the presence of chemical reaction is applied. The arising nonlinear partial differential equations are reduced to nonlinear ordinary differential equations by similarity variables. Similar solutions of velocity and concentration fields are derived by a homotopy analysis method. The values of surface mass transfer and gradient of mass transfer are also tabulated. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
20.
Vincent Guinot 《国际流体数值方法杂志》2003,41(11):1191-1219
Most existing algorithms for two‐dimensional shallow water simulations treat multi‐dimensional waves using wave splitting or time splitting. This often results in anisotropy of the computed flow. Both wave splitting and time splitting are based on a local decomposition of the multi‐dimensional problem into one‐dimensional, orthogonal problems. Therefore, these algorithms handle boundary conditions in a very similar way to classical one‐dimensional algorithms. This should be expected to trigger a dependence of the number of boundary conditions on the direction of the flow at the boundaries. However, most computational codes based on alternate directions do not exhibit such sensitivity, which seems to contradict the theory of existence and uniqueness of the solution. The present paper addresses these issues. A Riemann solver is presented that aims to convert two‐dimensional Riemann problems into a one‐dimensional equivalent Riemann problem (ERP) at the interfaces between the computational cells. The ERP is derived by applying the theory of bicharacteristics at each end of the interface and by performing a linear averaging along the interface. The proposed approach is tested against the traditional one‐dimensional approach on the classical circular dambreak problem. The results show that the proposed solver allows the isotropy of the solution to be better preserved. Use of the two‐dimensional solver with a first‐order scheme may give better results than use of a second‐order scheme with a one‐dimensional solver. The theory of bicharacteristics is also used to discuss the issue of boundary conditions. It is shown that, when the flow is subcritical, the number of boundary conditions affects the accuracy of the solution, but not its existence and uniqueness. When only one boundary condition is to be prescribed, it should not be the velocity in the direction parallel to the boundary. When two boundary conditions are to be prescribed, at least one of them should involve the component of the velocity in the direction parallel to the boundary. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献