Search results

Dams are constructed for many purposes such as for power generation, irrigation, water supply and flood control. However, dams can also impose risks to the public, and the situation could be disastrous if dam failure occurred. The study area, Bertam Valley, is located downstream of hydroelectric dam known as Sultan Abu Bakar Dam, Cameron Highlands. The key objectives of the study are to determine the potential risk area at downstream and to assess the flooding impact on damage to buildings and infrastructures due to dam break event. ArcGIS application and output from two-dimensional flood modelling have been used as an integrated approach to analyse the impact due to dam break flood, by creating flood severity grid analysis. The result obtained shows that the estimated inundated area is about 0.28 km², and almost 197 buildings are potentially affected. Results from this study show that in the event of dam break, the huge volume of impounding water will pound to the downstream areas, threatening the populations, and environment along its path. The finding is useful to assist the local authorities and emergency responders in formulating an emergency procedure to save the people during an emergency.

This study aims to predict the consequences associated with the propagation of the flood wave that may occur after the failure of the Taksebt dam and suggest an efficient emergency action plan for mitigation purposes.

To achieve the objectives of this study, the hydrodynamic model HEC-RAS 2D was used for the flood routing of the dam-break wave, which gave an estimate of the hydraulic characteristics downstream the Taksebt dam.

Geospatial analysis of the simulation results conducted in a geographic information system (GIS) environment showed that many residential areas are considered to be in danger in case of the Taksebt dam-break event. Based on the obtained results, an emergency actions plan was suggested to moderate the causalities in the downstream area at risk.

Overall, this study showed that the integration of 2D hydraulic modeling and GIS provides great capabilities in providing realistic view of the dam-break wave propagation that enhances assessing the associated risks and proposing appropriate mitigation measures.

The paper aims to apply a modified version of the MPS method to a double-dam-breaking test in which high dispersion zones and high natural clusterization zones are present, such as when the water column collapses into two sides and the two solitary waves collide, respectively.

The work takes advantage of the mixed source term from the cheaper computational version of the moving particle semi-implicit (MPS) method to reduce one step from the MPS classical algorithm. The proposed test can be successfully simulated by applying modifications to the variance parameter in the Laplacian operator and gradient model.

The results show stable behavior in dispersion and clusterization zones. Also, the collision and merging produced by solitary waves was successfully simulated.

The main limitation in this work was the development of a comparison between the obtained results and the simulations with the original cheaper computational version of the MPS, this limitation is due to the overestimation of inter particle repulsive forces from its gradient model.

The application of solitary waves is of paramount importance in a number of applications, and this stems from the fact that the interaction of solitary waves with ships and other floating structures could generate highly deformed and complex free surface flows.

For future work, the modified version of the MPS method can be applied in flow over sill base simulations, in close and open channels, and in simulating breaking waves to determine impact pressures by using solitary wave propagation.

The simulation of interaction of large groups of particles as in the case when two solitary waves collide could cause severe instability problems in pressure, causing the computer analysis to stop. MPS classical algorithm takes into account an explicit step that, in this case, may provoke the problem. For this reason, the cheaper version of MPS method is used to correctly simulate solitary wave interactions.

The purpose of this paper is to prove the validity of the front‐tracking variant of the lattice Boltzmann method (LBM) to simulate free surface hydraulic flows (i.e. dam break flows).

In this paper, an algorithm for free surface simulations with the LBM method is presented. The method is chosen for its computational efficiency and ability to deal with complex geometries. The LBM is combined to a surface‐tracking technique applied to a fixed Eulerian mesh in order to simulate free surface flows.

The numerical method is then validated against two typical cases of environmental‐hydraulic interest (i.e. dam break) by comparing LBM results with experimental data available in literature. The results show that the model is able to reproduce the observed water levels and the wave fronts with reasonable accuracy in the whole period of the transient simulations, thus highlighting that the present method may be a promising tool for practical dam break analyses.

Even if the main philosophy of the proposed method is equal to the volume of fluid technique usually coupled to Navier‐Stokes models, no additional differential equation is needed to determine the relative volume fraction of the two phases, or phase fraction, in each computational cell, as the free‐surface tracking is automatically performed. This results in a method very simple to be coded with high computational efficiency. The results presented in this paper are the first, to the best of the authors' knowledge, in the field of hydraulic engineering.

Flood Emergency Response Plan (ERP) is a plan that guides responsibilities for proper operation of Sultan Abu Bakar (SAB) dam in response to emergency incidents affecting the dam by high water storage capacity. Based on this study, four major responsibilities are needed for SAB dam owing to protect any probable risk for downstream which can be Incident Commander, Deputy Incident Commander, On-scene Commander and Civil Engineer. Having organisation charts based on ERP exercise can be helpful for decreasing the probable risks in any projects such as Abu Bakar Dam and it is a way to identify and suspected and actual dam safety emergencies. A dam safety emergency is an event, which could potentially lead to dam break and need to be taken care of a massive plan. To mitigate the hydro hazard due to dam break, UNITEN has developed a new application software known as INSPiRE (Interactive Dam Safety Decision Support System). INSPiRE, as an intelligent dam safety software, is developed to address emergencies, which demand fast, decision making and effective multi-agency collaboration due to SAB dam break event. INSPiRE will contribute towards the sustainability of SAB dam’s owner as corporate reputations can be ruined through dam structural failures that can affect the economy of the nation and enhance the quality of life of the people.

The purpose of this paper is to find the best solver for parallelizing particle methods based on solving Pressure Poisson Equation (PPE) by taking Moving Particle Semi-Implicit (MPS) method as an example because the solution for PPE is usually the most time-consuming part difficult to parallelize.

To find the best solver, the authors compare six Krylov solvers, namely, Conjugate Gradient method (CG), Scaled Conjugate Gradient method (SCG), Bi-Conjugate Gradient Stabilized (BiCGStab) method, Conjugate Gradient Squared (CGS) method with Symmetric Lanczos Algorithm (SLA) method and Incomplete Cholesky Conjugate Gradient method (ICCG) in terms of convergence, time consumption, parallel efficiency and memory consumption for the semi-implicit particle method. The MPS method is parallelized by the hybrid Open Multi-Processing (OpenMP)/Message Passing Interface (MPI) model. The dam-break flow and channel flow simulations are used to evaluate the performance of different solvers.

It is found that CG converges stably, runs fastest in the serial way, uses the least memory and has highest OpenMP parallel efficiency, but its MPI parallel efficiency is lower than SLA because SLA requires less synchronization than CG.

With all these criteria considered and weighed, the recommended parallel solver for the MPS method is CG.

The purpose of this paper is to present a new numerical model for shallow water flows over heterogeneous sedimentary layers. It is already several years since the single-layered models have been used to model shallow water flows over erodible beds. Although such models present a real opportunity for shallow water flows over movable beds, this paper is the first to propose a multilayered solver for this class of flow problems.

Multilayered beds formed with different erodible soils are considered in this study. The governing equations consist of the well-established shallow water equations for the flow, a transport equation for the suspended sediments, an Exner-type equation for the bed load and a set of empirical equations for erosion and deposition terms. For the numerical solution of the coupled system, the authors consider a non-homogeneous Riemann solver equipped with interface-tracking tools to resolve discontinuous soil properties in the multilayered bed. The solver consists of a predictor stage for the discretization of gradient terms and a corrector stage for the treatment of source terms.

This paper reveals that modeling shallow water flows over multilayered sedimentary topography can be achieved by using a coupled system of partial differential equations governing sediment transport. The obtained results demonstrate that the proposed numerical model preserves the conservation property, and it provides accurate results, avoiding numerical oscillations and numerical dissipation in the approximated solutions.

A novel implementation of sediment handling is presented where both averaged and separate values for sediment species are used to ensure speed and precision in the simulations.

This paper aims to develop a robust set of advanced numerical tools to simulate multiphase flows under the superimposition of external uniform magnetic fields.

The flow has been simulated in a fully Eulerian framework by a {\it variational multi-scale} method, which allows to take into account the small-scale turbulence without explicitly model it. The multi-fluid problem has been solved through the convectively re-initialized level-set method to robustly deal with high density and viscosity ratio between the phases and the surface tension has been modelled implicitly in the level-set framework. The interaction with the magnetic field has been modelled through the classic induction equation for 2D problems and the time step computation is based on the electromagnetic interaction to guarantee convergence of the method. Anisotropic mesh adaptation is then used to adapt the mesh to the main problem’s variables and to reach good accuracy with a small number of degrees of freedom. Finally, the variational multiscale method leads to a natural stabilization of the finite elements algorithm, preventing numerical spurious oscillations in the solution of Navier–Stokes equations (fluid mechanics) and the transport equation (level-set convection).

The methodology has been validated, and it is shown to produce accurate results also with a low number of degrees of freedom. The physical effect of the external magnetic field on the multiphase flow has been analysed.

The dam-break benchmark case has been extended to include magnetically constrained flows.

A numerical method is developed to capture the interaction of solid object with two-phase flow with high density ratios. The current computational tool would be the first step of accurate modeling of wave energy converters in which the immense energy of the ocean can be extracted at low cost.

The full two-dimensional Navier–Stokes equations are discretized on a regular structured grid, and the two-step projection method along with multi-processing (OpenMP) is used to efficiently solve the flow equations. The level set and the immersed boundary methods are used to capture the free surface of a fluid and a solid object, respectively. The full two-dimensional Navier–Stokes equations are solved on a regular structured grid to resolve the flow field. Level set and immersed boundary methods are used to capture the free surface of liquid and solid object, respectively. A proper contact angle between the solid object and the fluid is used to enhance the accuracy of the advection of the mass and momentum of the fluids in three-phase cells.

The computational tool is verified based on numerical and experimental data with two scenarios: a cylinder falling into a rectangular domain due to gravity and a dam breaking in the presence of a fixed obstacle. In the former validation simulation, the accuracy of the immersed boundary method is verified. However, the accuracy of the level set method while the computational tool can model the high-density ratio is confirmed in the dam-breaking simulation. The results obtained from the current method are in good agreement with experimental data and other numerical studies.

The computational tool is capable of being parallelized to reduce the computational cost; therefore, an OpenMP is used to solve the flow equations. Its application is seen in the following: wind energy conversion, interaction of solid object such as wind turbine with water waves, etc.

A high efficient CFD approach method is introduced to capture the interaction of solid object with a two-phase flow where they have high-density ratio. The current method has the ability to efficiently be parallelized.

An area of great interest in current computational fluid dynamics research is that of free-surface modelling (FSM). Semi-implicit pressure-based FSM flow solvers typically involve the solution of a pressure correction equation. The latter being computationally intensive, the purpose of this paper is to involve the implementation and enhancement of an algebraic multigrid (AMG) method for its solution.

All AMG components were implemented via object-oriented C++ in a manner which ensures linear computational scalability and matrix-free storage. The developed technology was evaluated in two- and three-dimensions via application to a dam-break test case.

AMG performance was assessed via comparison of CPU cost to that of several other competitive sparse solvers. The standard AMG implementation proved inferior to other methods in three-dimensions, while the developed Freeze version achieved significant speed-ups and proved to be superior throughout.

A so-called Freeze method was developed to address the computational overhead resulting from the dynamically changing coefficient matrix. The latter involves periodic AMG setup steps in a manner that results in a robust and efficient black-box solver.