scholarly journals A numerical geotechnical model for computing soil slides at banks of water reservoirs

2020 ◽  
Vol 11 (1) ◽  
Author(s):  
Nils Reidar B. Olsen ◽  
Stefan Haun

AbstractSoil slides can occur when the water level in a lake or a reservoir is lowered. This may take place in situations when a reservoir is flushed to remove sediments. The current study describes a three-dimensional numerical model used for the simulation of reservoir flushing that includes the slide movements. The geotechnical failure algorithms start with modelling the groundwater levels at the banks of the reservoir. A limit equilibrium approach is further used to find the location of the slides. The actual movement of the sediments is computed by assuming the soil to be a viscous liquid and by solving the Navier–Stokes equations. The resulting bed elevation changes from the slides are computed in adaptive grids that change as a function of water level, bed erosion and slide movements. The numerical model is tested on the Bodendorf reservoir in Austria, where field measurements are available of the bank elevations before and after a flushing operation. The results from the numerical simulations are compared with these observations. A parameter test shows that the results are very sensitive to the cohesion and less sensitive to the E and G modules of the soil.

2000 ◽  
Vol 31 (1) ◽  
pp. 57-72 ◽  
Author(s):  
N. R. B. Olsen ◽  
D. K. Lysne

A three-dimensional numerical model was used to model water circulation and spatial variation of temperature in Lake Sperillen in Norway. A winter situation was simulated, with thermal stratification and ice cover. The numerical model solved the Navier-Stokes equations on a 3D unstructured non-orthogonal grid with hexahedral cells. The SIMPLE method was used for the pressure coupling and the k-ε model was used to model turbulence, with a modification for density stratification due to the vertical temperature profile. The results were compared with field measurements of the temperature in the lake, indicating the location of the water current. Reasonably good agreement was found.


Energies ◽  
2021 ◽  
Vol 14 (3) ◽  
pp. 639
Author(s):  
Dennis Powalla ◽  
Stefan Hoerner ◽  
Olivier Cleynen ◽  
Nadine Müller ◽  
Jürgen Stamm ◽  
...  

The objective of the present paper is to develop a validated numerical model of a water vortex power plant that serves as a digital twin for further studies such as assessments of the ethohydraulic characteristics or the performance of such devices. The reference for the validation process is a large-scale hydraulic installation equipped with a full-scale water vortex power plant prototype installed in Dresden (Germany), where flow field measurements were carried out using three-dimensional Acoustic Doppler Velocimetry. The numerical model was implemented within the software package Star-CCM+. The unsteady, two-phase flow was solved with the Reynolds-Averaged Navier–Stokes equations in a Eulerian Multiphase approach, deploying a Volume of Fluid method to describe the free-surface flow. Water level and flow velocities were systematically compared in key areas of the device, demonstrating that the simulation is in good agreement with experimental observations. Relative differences are limited to at most 4% regarding water height in the system, and even the much more challenging velocity fields are reproduced with typical relative errors of roughly 10%. This validates the ability of the model to model the challenging flow conditions found in a water vortex power plant, enabling subsequent studies of the characteristics of this power plant concerning fish migration.


2010 ◽  
Vol 41 (2) ◽  
pp. 92-103 ◽  
Author(s):  
Peggy Zinke ◽  
Nils Reidar Bøe Olsen ◽  
Jim Bogen ◽  
Nils Rüther

A 3D numerical model was used to compute the discharge distribution in the channel branches of Lake Øyeren's delta in Norway. The model solved the Navier–Stokes equations with the k–ɛ turbulence model on a 3D unstructured grid. The bathymetry dataset for the modelling had to be combined from different data sources. The results for three different flow situations in 1996 and 1997 showed a relative accuracy of the computed discharges within the range of 0 to±20% compared with field measurements taken by an ADCP at 13 cross sections of the distributary channels. The factors introducing the most error in the computed results are believed to be uncertainties concerning the bathymetry. A comparison between the computational results of the older morphology data from 1985–1990 and the model morphology from 1995–2004 indicated that morphological changes in this period had already had consequences for the flow distribution in some channels. Other important error sources were the inevitable use of averaged water level gradients because of unavailable water level measurements within the delta.


2013 ◽  
Vol 353-356 ◽  
pp. 2496-2501
Author(s):  
Biao Lv

A three dimensional non-hydrostatic numerical model is presented based on the incompressible Navier-Stokes equations and mass transport equations. An unstructured finite-volume technique is used to discretized the governing equations with good adaptable to complicated boundary. A conservative scalar transport algorithm is also applied in this model. An integral method of the top- layer pressure is applied to reduce the number of vertical layers. Three classical examples including periodic waves propagating over a submerged bar and non-hydrostatic lock exchange are used to demonstrate the capability and efficiency of the model. The simulation results are in good agreement with the analytical solution and experimental data.


Author(s):  
Nina Wolfrum ◽  
Giovanni Brignole ◽  
Karl Engel

A numerical model has been developed to reproduce the effects of complex casing treatments (CT) in steady RANS simulations of multistage compressors. While some CTs, such as circumferential grooves, can be described by a rotation surface and can thus easily be included in conventional steady simulations, the CFD analysis of other casing treatments like axial slot or recessed vanes, currently requires a time-resolving analysis of the interaction between such structures and rotating parts. At present unsteady simulations are still too time consuming to be used in the early phase of a compressor design. In the presented study a numerical model was developed for casing treatment applications, to introduce the unsteady effects caused by such casing treatments into steady CFD-simulations. With the help of the model, non-axisymmetric elements can be eliminated from the geometry allowing a steady simulation to be used. The flow acceleration and redirection caused by these geometrical elements is replaced with adequate source terms introduced into the three-dimensional Navier-Stokes equations. These source terms, derived from a consecutive time- and circumferential averaging of the three-dimensional unsteady Reynolds-averaged Navier-Stokes-equations, arise from the momentum and energy equations. Using these additional terms, the CT-model simulates both the pressure forces that the walls of the real casing treatment exert on the flow, and the effects of the mean blockage induced by the omitted geometry. Furthermore, the deterministic stresses, caused by a circumferentially inhomogeneous flow within the CT-structure, are modeled. The source terms consist of geometrical data that can be derived directly from the real geometry of the casing treatment as well as physical quantities of the time-averaged flow in the real casing treatment. The latter terms can be obtained from a reference unsteady simulation. In the presented case one unsteady simulation was sufficient to set up the model for a complete speed line. The model was implemented into the three-dimensional Navier-Stokes-code TRACE [5][12]. By using steady instead of unsteady CFD simulations, the time required for a speedline computation was reduced by a factor of 10. At the same time, the numerical results of the CT-model showed good alignment with the reference data. The model was evaluated for several different styles of compressors. In this paper various results are presented, including speedlines as well as radial inflow- and outflow-profiles.


Author(s):  
Roberto Biollo ◽  
Ernesto Benini

The progress of numerical methods and computing facilities has led to using Computational Fluid Dynamics (CFD) as a current tool for designing components of gas turbine engines. It is known, however, that a sophisticated numerical model is required to well reproduce the many complex flow phenomena which characterize compression systems, such as shock waves and their interactions with boundary layers and tip clearance flows. In this work, the flow field inside the NASA Rotor 37, a well known test case representative of complex three-dimensional viscous flow structures in transonic bladings, was simulated using a commercial CFD code based on the 3-D Reynolds-averaged Navier-Stokes equations. In order to improve the accuracy of predictions, different aspects of the numerical model were analyzed; in particular, an attempt was made to understand the influence of grid topology, number of nodes and their distribution, turbulence model, and discretization scheme of numerical solution on the accuracy of computed results. Existing experimental data were used to assess the quality of the solutions. The obtainment of a good agreement between computed and measured performance maps and downstream profiles was clearly shown. Also, detailed comparisons with experimental results indicated that the overall features of the three-dimensional shock structure, the shock-boundary layer interaction, and the wake development can be calculated very well in the numerical approach for all the operating conditions. The possibility for a numerical model to better understand the aerodynamic behaviour of existing transonic compressors and to help the design of new configurations was demonstrated. It was also pointed out that the development of an accurate model requires the knowledge of both the physical phenomena place within the flow field and the features of the code which model them.


2020 ◽  
Vol 14 (4) ◽  
pp. 7369-7378
Author(s):  
Ky-Quang Pham ◽  
Xuan-Truong Le ◽  
Cong-Truong Dinh

Splitter blades located between stator blades in a single-stage axial compressor were proposed and investigated in this work to find their effects on aerodynamic performance and operating stability. Aerodynamic performance of the compressor was evaluated using three-dimensional Reynolds-averaged Navier-Stokes equations using the k-e turbulence model with a scalable wall function. The numerical results for the typical performance parameters without stator splitter blades were validated in comparison with experimental data. The numerical results of a parametric study using four geometric parameters (chord length, coverage angle, height and position) of the stator splitter blades showed that the operational stability of the single-stage axial compressor enhances remarkably using the stator splitter blades. The splitters were effective in suppressing flow separation in the stator domain of the compressor at near-stall condition which affects considerably the aerodynamic performance of the compressor.


Sign in / Sign up

Export Citation Format

Share Document