Experiment and simulation of turbulent flow in local scour around a spur dyke

The turbulent flow in the local scour hole around a single non-submerged spur dyke is investigated with both experimental and numerical methods. The experiments are conducted under clear-water scour regime with an impermeable spur dyke. The scour geometry and flow velocities are measured in details with a high-resolution laser displacement meter, electro-magnetic velocimetries and PIV （Particle image velocimetry）. A 3D non-linear k-ε model is developed to simulate the complex local flow field around the scour area. The numerical model is formulated using FVM （Finite volume method） on a collocated unstructured mesh, capable of resolving complex geometries and boundaries. It is found that the simulation results are reasonably consistent with those of the experimental measurements. Based on the study results, the nature of the flow structure around a spur dyke with local scour hole is analyzed.     

Well-balanced high-order centered schemes on unstructured meshes for shallow water equations with fixed and mobile bed

In this paper, we study the numerical approximation of the two-dimensional morphodynamic model governed by the shallow water equations and bed-load transport following a coupled solution strategy. The resulting system of governing equations contains non-conservative products and it is solved simultaneously within each time step. The numerical solution is obtained using a new high-order accurate centered scheme of the finite volume type on unstructured meshes, which is an extension of the one-dimensional PRICE-C scheme recently proposed in Canestrelli et al. (2009) [5]. The resulting first-order accurate centered method is then extended to high order of accuracy in space via a high order WENO reconstruction technique and in time via a local continuous space–time Galerkin predictor method. The scheme is applied to the shallow water equations and the well-balanced properties of the method are investigated. Finally, we apply the new scheme to different test cases with both fixed and movable bed. An attractive future of the proposed method is that it is particularly suitable for engineering applications since it allows practitioners to adopt the most suitable sediment transport formula which better fits the field data.     

ParBreZo: A parallel,unstructured grid,Godunov-type,shallow-water code for high-resolution flood inundation modeling at the regional scale

Topographic data are increasingly available at high resolutions (<10 m) over large spatial extents to support detailed flood inundation modeling and loss estimation analyses required for flood risk management. This paper describes ParBreZo, the parallel implementation of a two-dimensional, Godunov-type, shallow-water code, to address the computational demand of high-resolution flood modeling at the regional scale (10²–10⁴ km²). A systematic approach to unstructured grid partitioning (domain decomposition) is presented, and the Single Process Multiple Data (SPMD) paradigm of distributed-memory parallelism is implemented so the code can be executed on computer clusters with distributed memory, shared memory, or some combination of the two (now common with multi-core architectures). In a fully-wetted, load-balanced test problem, the code scales very well with a parallel efficiency of close to 100% on up to 512 processes (maximum tested). A weighted grid partitioning is used to partially address the load balancing challenge posed by partially wetted domains germane to flooding applications, where the flood extent varies over time, while the partitioning remains static. An urban dam-break flood test problem shows that weighted partitions achieve a parallel efficiency exceeding 70% using up to 48 processes. This corresponds to a 97% reduction in execution time so results are obtained in a matter of minutes, which is attractive for routine engineering analyses. A hurricane storm surge test problem shows that a 10 m resolution, 12 h inundation forecast for a 40 km length of coastline can be completed in under 2 h using 512 processors. Hence, if coupled to a hurricane forecast system capable of resolving storm surge, inundation forecasts could be made at 10 m resolution with at least a 10 h lead time.     

边坡稳定强度折减有限元分析中的若干问题讨论

强度折减有限元法在边坡稳定性分析中得到了广泛应用,但仍存在一些问题未能很好地解决。首先探讨了不同单元类型对边坡稳定安全系数计算精度的影响,指出采用三角形二阶单元,能够较好地模拟土坡的渐近破坏过程,且相对网格大小控制在0.1～0.15之间时能获得较精确的解;其次讨论了边坡失稳的4个判据问题,指出在合理网格大小区间内,有限元强度折减至土坡破坏时,先后出现塑性区贯通、等效塑性应变带贯通、位移陡增与数值计算不收敛等4种表征,以此分别作为边坡失稳判据,所得到的安全系数具有较好的一致性;最后讨论了边坡破坏时滑动面深浅、滑出点位置与坡体材料的黏聚力、内摩擦角及坡角间的关系,指出黏聚力越大、内摩擦角越小或坡度越小,滑动面越深,且滑出点越远离坡脚。     

自动网格法在大型滑坡模型试验位移测试中的应用

应用非接触式光学测量方法中自动网格法原理,研究了大型滑坡模型试验中的位移测试方法及其分布规律。以石榴树包滑坡模型试验为例,针对自动网格法用于大型滑坡模型试验位移测量时所存在的问题提出了解决措施,并对石榴树包滑坡在不同水位运行工况下的位移规律进行了研究分析。     

土石坝二维地震反应有限元计算的影响因素研究

用有限元法计算土石坝地震反应时,有限元模型的网格剖分、计算边界的处理是计算精度的主要影响因素,也是有限元法分析前处理的主要工作。基于河海大学岩土工程研究所开发的二维有限元土石坝计算程序,分析了网格剖分的形式、网格的合理尺寸、简单边界的计算取值范围对土石坝地震反应的影响。研究结果表明,在土石坝地震反应计算中,网格剖分的整体形式应该水平成层,网格尺寸以最大坝高的7 %～10 %为上限;在有覆盖层的土石坝工程中,简单边界的计算区域应包括覆盖层,并与覆盖层厚度有特定的关系。上述成果可使前处理过程更严谨、更易操作,可作为研究和实际工程计算的参考。     

A visualization-oriented 3D method for efficient computation of urban solar radiation based on 3D–2D surface mapping

The temporal and spatial distribution of solar energy in urban areas is highly variable because of the complex building structures present. Traditional GIS-based solar radiation models rely on two-dimensional (2D) digital elevation models to calculate insolation, without considering building facades and complicated three-dimensional (3D) shading effects. Inspired by the ‘texture baking’ technique used in computer graphics, we propose a full 3D method for computing and visualizing urban solar radiation based on image-space data representation. First, a surface mapping approach is employed to project each 3D triangular mesh onto a 2D raster surface whose cell size determines the calculation accuracy. Second, the positions and surface normal vectors of each 3D triangular mesh are rasterized onto the associated 2D raster using barycentric interpolation techniques. An efficient compute unified device architecture -accelerated shadow-casting algorithm is presented to accurately capture shading effects for large-scale 3D urban models. Solar radiation is calculated for each raster cell based on the input raster layers containing such information as slope, aspect, and shadow masks. Finally, a resulting insolation raster layer is produced for each triangular mesh and is represented as an RGB texture map using a color ramp. Because a virtual city can be composed of tens of thousands of triangular meshes and texture maps, a texture atlas technique is presented to merge thousands of small images into a single large image to batch draw calls and thereby efficiently render a large number of textured meshes on the graphics processing unit.     

A semi-implicit three-dimensional numerical model for non-hydrostatic pressure free-surface flows on an unstructured,sigma grid

A semi-implicit 3-D numerical formulation for solving non-hydrostatic pressure free-surface flows on an unstructured,sigma grid is proposed.Pressure-splitting and 9 semi-implicit methods are inherited and reformed from Casulli’s z-coordinate formulation.The non-orthogonal sigma-coordinate transformation leads to additional terms. The resulting linear system for the non-hydrostatic correction is diagonally dominant but unsymmetric,and it is solved by the BiCGstab method.In contrast with z-coordinate non-hydrostatic models,the new model fits vertical boundaries much better,which is important for the long-time simulation of sediment transport and riverbed deformation.A lock-exchange density flow is computed to determine whether the new scheme is able to simulate non-hydrostatic free-surface flows.The new model is further verified using the field data of a natural river bend of the lower Yangtze River.Good agreement between simulations and earlier research results,field data is obtained, indicating that the new model is applicable to hydraulic projects in real rivers.     

最小二乘配点法解一维非稳定流问题

最小二乘配点无网格法是一种新型高效的无网格法。该方法除节点外又在研究域内引入辅助点,近似函数仍然只通过节点构造,微分方程在所有节点和辅助点上满足。用最小二乘配点法计算河间承压非稳定流问题,算例表明,最小二乘配点法比有限差分法计算精度高,稳定性好。