Search

Coastal lagoons are particularly vulnerable to climate change, in particular, Sea Level Rise (SLR) due to their shallowness. Lake Burullus provides a variety of socio-economic services as the second largest coastal lagoon in Egypt. Recently, it has experienced significant ecological deterioration. Thus, its ecosystem is fragile in the face of anthropogenic induced changes. The main objective of the current study is to investigate the climate change impacts on characteristics of Lake Burullus. A depth averaged hydro-ecological modeling system, MIKE21, was applied to develop an eco - hydrodynamic model for the lake. The developed model was calibrated and verified for two successive years: July 2011–June 2012 and July 2012–June 2013. The model simulations exhibited good agreement with the measurements during the calibration and verification processes. Six different Regional Climate Models (RCMs) were compared, using six different statistical metrics, to determine the most accurate one for the study area. The required meteorological input, including surface air temperature, precipitation, and evaporation were derived from the selected RCM. The meteorological input was extracted for two different years in the 21 st century considering one Representative Concentration Pathways (RCPs) scenario, based on the Intergovernmental Panel on Climate Change (IPCC) 5th Report. Regional SLR projections for the Mediterranean Sea for the selected RCP scenario and the two studied years were obtained. These future climate change estimates were used to modify the validated model of the lake. A sensitivity analysis was applied to assess effect of future climatic conditions and SLR, separately. The results revealed that the lake water depths will increase and it will be warmer and more saline. Significant spatial variability of the studied parameters under climate change forcing is expected. Consequently, climate change is going to restrict the lake's ability to preserve the present-day species. An urgent management plan involving adaptation works, should be implemented to reduce such potential species losses in Egyptian lagoons.

An analytical model was developed to predict the boundary shear stress and the depth-averaged streamwise velocity in the breach. A more general expression of the depth-integrated form of the streamwise Reynolds averaged Navier-Stokes (RANS) equation was rewritten by considering the secondary flow effects. It accounts for secondary flow effects via the dimensionless secondary flow parameters k1 and k2. The averaged lateral shear stress, the bottom shear stress and the depth-averaged streamwise velocity of the breach was obtained. According to sensitivity analysis of parameters, it is obtained that the secondary flow parameter k2 has a great influence on the averaged lateral shear stress, but few influences on the bottom shear stress and the depth-averaged streamwise velocity. The dimensionless streamwise velocity variation parameter c1 has great influences on the averaged lateral shear stress, bottom shear stress and the depth-averaged streamwise velocity. Comparison of the analytical results with experimental data shows that the proposed model predicts the averaged lateral shear stress and the bottom shear stress well.

Estimation of streamwise depth-averaged velocity and boundary shear stress are important requisites for modelling of river flows associated with flood events. Recently, many methods have emerged to predict these flow variables with great accuracy, but they provide unsatisfactory results in shear layer regions. This paper presents an improved methodology to predict depth-averaged velocity and bed shear stress for a straight compound channel flow. An analytical solution to the depth-integrated turbulent form of the Navier–Stokes equation is obtained. The transverse shear stress in the mixing region is modelled using an effective eddy viscosity concept that contains horizontal coherent structures and three-dimensional bottom turbulence. The secondary flow term is modelled by considering the log-law profile for streamwise velocity and half cosine curve for the transverse velocity component. The analytical solution is successfully applied to a wide range of experimental compound channels and field cases. The efficacy of the present solution has been successfully tested by comparing with observed values.

The distribution of depth-averaged velocity is crucial to many engineering problems. This paper presents an analytical solution for the lateral depth-averaged velocity distribution in a 120° curved channel with asymmetric trapezoidal cross-sections in which the curved flow is not fully developed. The proposed analytical model was derived from the streamwise depth-integrated Reynolds-averaged Navier–Stokes equation, which considers the bed stress, Reynolds stress, secondary flow and local flow acceleration. In addition, practical boundary conditions are proposed to consider the trapezoidal cross-sections with linearly varying side beds and the flow separation at the inner bank near the bend entrance. The analytical solutions using the proposed model were found to be in good agreement with measured data. The main contributions of this study are the development of a conventional analytical model for sophisticated circumstances and the proposal of practical boundary conditions that consider flow separation.

This study proposes a numerical model based on the depth-integrated non-hydrostatic shallow water equations with an improvement of wave breaking dissipation. Firstly, studies of parameter sensitivity were carried out using the proposed numerical model for simulation of wave breaking to understand the effects of the parameters of the breaking model on wave height distribution. The simulated results of wave height near the breaking point were very sensitive to the time duration parameter of wave breaking. The best value of the onset breaking parameter is around 0.3 for the non-hydrostatic shallow water model in the simulation of wave breaking. The numerical results agreed well with the published experimental data, which confirmed the applicability of the present model to the simulation of waves in near-shore areas.

Along sandy coastlines, submerged, shore-parallel sandbars play an essential role in shoreline morphology by dissipating wave energy through depth-induced wave breaking. While wave breaking and sediment transport around sandbars are complex three-dimensional (3D) processes, shoreline morphology is typically simulated with depth-averaged models that feature lower computational demand than do 3D models. In this context, this study examines the implications of depth-averaging the flow field and approximating the breaking process in nonhydrostatic models (e.g., XBeach nonhydrostatic) for the hydro- and morphodynamic processes around sandbars. The implications are drawn based on reproducing large-scale experiments of a barred beach profile using the single-layer (XBNH) and the reduced two-layer (XBNH+) modes of XBeach. While hydrodynamic processes were predicted with high accuracy on the sandbar's seaward side, wave heights were overpredicted on the bar's landward side. The overestimation was due to the simplified reproduction of the complex breaking process near the sandbar's peak, particularly in terms of the generated turbulence in the water column. Moreover, the velocity profile with a strong undertow could only be represented in a simplified way even using the two-layer mode XBNH+, thus resulting in inaccurate predictions of sediment loads around the sandbar. A parametric study is performed, and it revealed which model parameters control the simulation of the wave-breaking process. Thus, wave height predictions could be improved by tuning the energy-dissipation parameters. However, flow velocities and morphodynamic predictions could not be improved accordingly. Thus, this study identifies possible hydrodynamic model improvements, such as incorporating a roller dissipation model. Moreover, it improves understanding of key drivers and processes that should be included in nonhydrostatic depth-averaged models to simulate morphological changes around sandbars more efficiently.

Summary Two popular and computationally inexpensive class of methods for approximating the propagation of surface waves over two-dimensional variable bathymetry are ‘step approximations’ and ‘depth-averaged models’. In the former, the bathymetry is discretised into short sections of constant depth connected by vertical steps. Scattering across the bathymetry is calculated from the product of $2 \times 2$ transfer matrices whose entries encode scattering properties at each vertical step taken in isolation from all others. In the latter, a separable depth dependence is assumed in the underlying velocity field and a vertical averaging process is implemented leading to a second-order ordinary differential equation (ODE). In this article, the step approximation is revisited and shown to be equivalent to an ODE describing a depth-averaged model in the limit of zero-step length. The ODE depends on how the solution to the canonical vertical step problem is approximated. If a shallow water approximation is used, then the well-known linear shallow water equation results. If a plane-wave variational approximation is used, then a new variant of the mild-slope equations is recovered.

The published version of this article, unfortunately, contained error. Figure 3 of “GPU-hosted workstation for parallel computing” was lost. Then the sequences of the Figures 4-10 were wrong. Given in this article are the correct figures.

Several numerical models have been developed to simulate the dispersion of radionuclides in the aqueous environments. In this research, the 2-dimensional depth averaged hydrodynamic model is applied to obtain the current components in a part of Persian Gulf taking into account the tidal effects. The hydrodynamic equations solved by a finite differential approach. In this model the time step is 5 s (totally for 3 days) and spatial step is 10 Km. From the mainstream in the study area, it is also found that the velocity of water in the central region and near Iran's shore lines is higher than the other places. Using the average value of the speed from hydrodynamic equation, the concentration of some radioactive elements is also obtained by CROM tool.

Two different mathematical models of fluid mechanics are now being investigated at the WSL Institute for Snow and Avalanche Research in Davos to model powder-snow avalanches. The first approach is to solve the full three-dimensional multiphase (ice-dust, air) incompressible Navier-Stokes equations; the second approach is to apply depth-averaged models to simulate both the formation and independent propagation of the powder cloud. The final goal of both models is to predict the dynamics of powder avalanches in three-dimensional terrain and specifically cloud impact pressures. Both models are driven by the same set of terrain dependent mass and momentum exchanges defined by the flow state (speed, density, height) of the avalanche core. The great advantage of the depth-average approach is computational speed, allowing the investigation of different hazard scenarios involving variable release locations, snow temperature and entrainment depths. This fact has allowed the widespread application of the depth-average model to many historical case studies, including the avalanches measured at the Vallée de la Sionne (VdlS) test site. However, a central modelling problem needs to be resolved: both air-entrainment (cloud height and density) and drag (cloud speed) are intimately linked to the turbulence created during the cloud formation phase.In this presentation, we present a depth-averaged turbulence model proposed by V. M. Teshukov [1] and extended by Richard and Gavrilyuk [2] and Gavrilyuk et al. [3], Ivanova et al. [5, 6]. The mathematical model is a 2D hyperbolic non-conservative system of equations that is mathematically equivalent to the Reynolds-averaged model of barotropic turbulent flows. The system is non-conservative, extending the classical shallow water equations to contain three independent components of the symmetric Reynolds stress tensor. We simulate the measured powder cloud heights of two VdlS avalanches using both the incompressible Navier-Stokes and turbulent shallow-water models, capturing the unsteady formation of billow height and width measured by ground based photogrammetry [4]. This can only be achieved by making air-entrainment dependent on the vorticity predicted by the turbulence model. We conclude by summarizing why we believe shallow-water type models can be applied for practical hazard engineering problems.References:[1] V. M. Teshukov in "Gas-dynamics analogy for vortex free-boundary flows.", J. Appl. Mech. Tech. Phys., 2007.[2] G. L. Richard, S. L. Gavrilyuk in "A new model of roll waves: comparison with Brock’s experiments", Journal of Fluid Mechanics, 2012.[3] S.L. Gavrilyuk, K.A. Ivanova, N. Favrie in "Multi-dimensional shear shallow water flows : problems and solutions", Journal of Computational Physics, 2018.[4] Dreier, L., Bühler, Y., Ginzler, C., and Bartelt, P.: Comparison of simulated powder snow avalanches with photogrammetric measurements, Annals of Glaciology, 57, 371 - 381, 10.3189/2016AoG71A532, 2016.][5] K.A. Ivanova, S.L. Gavrilyuk, ”Structure of the hydraulic jump in convergent radial flows”,Journal of Fluid Mechanics, Volume 860, 10 February2019 , pp. 441-464.[6] K.A. Ivanova, S.L. Gavrilyuk, B. Nkonga, G.L. Richard, ”Formation and coarsening of roll-waves in shear shallow water flows down an inclinedrectangular channel”,Computers& Fluids, 159, pp 189203, 2017

Debris flows are landslide phenomena which occur worldwide, posing a major threat to mountain settlements. They consist of flowing fine and coarse sediment saturated with water, which propagate mainly in channelized paths. Because of their high velocity and unpredictability, the evacuation of local populations is often impossible. Losses of human lives and economical damages can be avoided if a correct risk mitigation procedure is adopted. Hence, mitigation structures, such as filter barriers or flexible barriers are often installed in high-risk areas. The primary goal of these structures is to reduce the flow energy and to retain the coarsest boulders. Their design process, which is still frequently based only on empirical or simplified models, would greatly benefit from the support of a reliable numerical model. In this framework, continuum-based Depth-Averaged Models (DAMs) have been the dominant numerical tool since the 90s. DAMs can simulate events propagating over a wide area while keeping the computational time low, even on complex topographies (Pirulli, 2010). Nevertheless, the averaging process applied to velocity and pressure causes a loss of information, which is critical when the flow impact against structures is evaluated. A full 3D model would allow for a more accurate resolution of fluid-structure interaction (Leonardi et al., 2016). However, debris flows may propagate up to kilometres, and a complete 3D analysis would therefore require exceedingly long computational times. To bypass the shortcomings mentioned above, this work aims to couple DAMs to a 3D model based on the Lattice Boltzmann Method (LBM). Thus, the domain is split into two parts. First, DAMs describes the flow evolution from its initialization to the transport phase. In this portion of the domain, no structures are present. When the flow approaches a structure, DAMs is coupled to a 3D model. To verify the coupling procedure accuracy, the model is benchmarked on the laboratory tests conducted by Moriguchi et al. (2009). These laboratory tests targeted the flow of dry sand on a steep chute, evaluating the flow impact on a barrier. Preliminary results suggest that the coupled model reproduces the laboratory results reasonably well. Keywords: debris flow, coupled numerical modelling, depth-averaged method, 3D Lattice-Boltzmann Method REFERENCESLeonardi, A., Wittel, F. K., Mendoza, M., Vetter, R., & Herrmann, H. J. (2016). Particle-Fluid-Structure Interaction for Debris Flow Impact on Flexible Barriers. Computer-Aided Civil and Infrastructure Engineering, 31(5), 323–333.Moriguchi, S., Borja, R. I., Yashima, A., & Sawada, K. (2009). Estimating the impact force generated by granular flow on a rigid obstruction. Acta Geotechnica, 4(1), 57–71.Pirulli, M. (2010). On the use of the calibration-based approach for debris-flow forward-analyses. Natural Hazards and Earth System Science, 10(5), 1009–1019.

A two-dimensional diffusive wave framework is proposed for modeling interaction between flow through and outside porous media. The depth-integrated parabolic model was developed with an ass...

An efficient time-adaptive multigrid algorithm is used to solve a range of normal and oblique droplet impacts on dry surfaces and liquid films using the Depth-Averaged Form (DAF) method of the governing unsteady Navier–Stokes equations. The dynamics of a moving three-phase contact line on dry surfaces is predicted by a precursor film model. The method is validated against a variety of experimental results for droplet impacts, looking at factors such as crown height and diameter, spreading diameter and splashing for a range of Weber, Reynolds and Froude numbers along with liquid film thicknesses and impact angles. It is found that, while being a computationally inexpensive methodology, the DAF method produces accurate predictions of the crown and spreading diameters as well as conditions for splash, however, underpredicts the crown height as the vertical inertia is not included in the model.

Flow-like landslides commonly happen in mountainous areas and may cause economic and human life losses in the impacted areas. Computer modelling has become an effective tool for landslide risk assessment and reduction. Models based on discrete element method (DEM) have been widely used for landslide prediction; however, this method is computationally too demanding for large-scale applications. Depth-averaged models (DAMs) have been widely reported for simulating run-out and deposition of flow-like landslides over large spatial domains due to its relatively higher computational efficiency. To combine the advantages of both types of modelling approaches, this work introduces a novel landslide model developed by coupling a DEM model with DAM for simulation of flow-like landslides, in which the DEM is employed in the landslide initiation area to better simulate the failure mechanism of slope, and the DAM is adopted in the landslide runout and deposition phase, where the landslide has developed into flow-like landslide with fluid-like behaviour. Finally, the new coupled model is validated against an experimental test case. Satisfactory results have been obtained, demonstrating that the coupled model is able to accurately capture the detailed dynamics of flow-like landslides.

In this paper, a novel combination method based on the AdaBoost algorithm and the general regression neural network (GRNN) is proposed to predict the depth-averaged current (DAC) for the next profile using the underwater gliders (UGs). Firstly, considering changes of the seawater density and the hull deformation, the velocity calculation model of UGs and the calculation method of DAC are given. Secondly, two DAC forecasting models are established with historical data using the support vector machine (SVM) and the back-propagation neural network (BPNN) methods, respectively. Then, the forecasting models are combined nonlinearly by adopting the GRNN. Finally, outputs from the GRNN are combined using the AdaBoost algorithm, which further taken for the DAC forecasting. A sea trial was conducted using Petrel-II glider in the South China Sea, which verified the accuracy of the three forecasting models. Results demonstrate that the proposed combination method has a better performance in DAC forecasting than the other two single forecasting models.

This paper analyses an important aspect of the continuum numerical modelling of rapid landslides as debris flows: “By using the same rheological parameter values, are the results obtained with codes that implement the same constitutive equations, but a different numerical solver, equal?” To answer this question, the two numerical codes RASH3D and GeoFlow_SPH are used here to back-analyse the debris flow event that occurred in the Nora stream (northwestern Italian Alps) in October 2000. Comparison of results evidenced that the RASH3D best-fit rheological values for the Nora event back-analysis overestimated both the final depositional heights and the simulated flow velocities if used in GeoFlow_SPH. To obtain thickness values comparable with those measured in situ, it was necessary to re-calibrate GeoFlow_SPH rheological parameter values. This way, with the exception of a larger lateral spreading of the sliding mass given by RASH3D, both thickness and velocity values were similar for the two numerical codes.

Precise prediction of bed shear stress (τ0) and velocity profile in open channels is critical in various applications. This paper deliberates on the depth-averaged von Karman coefficient in...

The study establishes different analytical model for velocity distribution, boundary shear stress, stage–discharge curves and their application to the trapezoidal compound channel. The analytical s...

A depth-averaged two-phase model is proposed for debris flows over fixed beds, explicitly incorporating interphase and particle-particle interactions, fluid and solid fluctuations and multi grain sizes. A first-order model based on the kinetic theory of granular flows is employed to determine the stresses due to solid fluctuations, while the turbulent kinetic energy - dissipation rate model is used to determine the stresses from fluid fluctuations. A well-balanced numerical algorithm is applied to solve the governing equations. The present model is benchmarked against USGS experimental debris flows over fixed beds. Incorporating the stresses due to fluid and solid fluctuations and properly estimating the bed shear stresses are shown to be crucial for reproducing the debris flows. Longitudinal particle segregation is resolved, demonstrating coarser sediments around the fronts and finer grains trailing the head. Based on extended modeling exercises, debris flow efficiency is shown to increase with initial volume, which is underpinned by observed datasets.

A tidal propagation characteristic of Hooghly estuary is presented using numerical (ADCIRC) and analytical models (Friedrichs and Aubrey, 1994) along with observations. The analytical model is based on Friedrichs and Aubrey (1994) that simplifies the governing hydrodynamic equations greatly by retaining only those terms that are significant without losing the overall understanding of the propagation process. The analytical model is compared with corresponding 2-D depth averaged numerical (ADCIRC) model that retains all non-linear terms. The assumptions for simplification are found to be reasonable in the light of close agreement among analytical, numerical models and observations. A plan-form geometrical characteristic as well as hydrodynamic variable of the Hooghly has been compared with that of Delaware estuary for corroborating similar tidal propagation process. The Hooghly estuary has flood dominant asymmetric tidal propagation and a positive amplitude growth factor (μ). The observed tidal celerity (phase speed) on an average is slightly more than frictionless celerity. Using the conventions of Toffolon et al. (2006), Hooghly can be classified into ‘strongly convergent – strongly dissipative’ estuary. From the results it can be construed that the estuary is yet to stabilise and reach its equilibrium morphology. It can be close to its equilibrium as very little amplification (0.1 m) is noticed in the predominant semi-diurnal constituent M2 over 78 km (barely 7%) in the estuary. The parameters of width variation (γ) and the ratio between friction and inertia (χ) have been used to define the marginal condition for amplification. The relative position of Hooghly in terms of marginal condition is consistent with similar set of estuaries elsewhere that have been grouped using the above parameters.

We use Lattice Boltzmann Method (LBM) MRT and Cumulant schemes to study the performance and accuracy of single-phase flow modeling for propped fractures. The simulations are run using both the two- and three-dimensional Stokes equations, and a 2.5D Stokes–Brinkman approximate model. The LBM results are validated against Finite Element Method (FEM) simulations and an analytical solution to the Stokes–Brinkman flow around an isolated circular obstacle. Both LBM and FEM 2.5D Stokes–Brinkman models are able to reproduce the analytical solution for an isolated circular obstacle. In the case of 2D Stokes and 2.5D Stokes–Brinkman models, the differences between the extrapolated fracture permeabilities obtained with LBM and FEM simulations for fractures with multiple obstacles are below 1%. The differences between the fracture permeabilities computed using 3D Stokes LBM and FEM simulations are below 8% . The differences between the 3D Stokes and 2.5 Stokes–Brinkman results are less than 7% for FEM study, and 8% for the LBM case. The velocity perturbations that are introduced around the obstacles are not fully captured by the parabolic velocity profile inherent to the 2.5D Stokes–Brinkman model.

The traditional Alternating Direction Implicit (ADI) scheme has been proven to be incapable of modelling trans-critical flows. Its inherent lack of shock-capturing capability often results in spurious oscillations and computational instabilities. However, the ADI scheme is still widely adopted in flood modelling software, and various special treatments have been designed to stabilise the computation. Modification of the Boussinesq coefficient to adjust the amount of fluid inertia is a numerical treatment that allows the ADI scheme to be applicable to rapid flows. This study comprehensively examines the impact of this numerical treatment over a range of flow conditions. A shock-capturing TVD-MacCormack model is used to provide reference results. For unsteady flows over a frictionless bed, such as idealised dam-break floods, the results suggest that an increase in the value of the Boussinesq coefficient reduces the amplitude of the spurious oscillations. The opposite is observed for steady rapid flows over a frictional bed. Finally, a two-dimensional urban flooding phenomenon is presented, involving unsteady flow over a frictional bed. The results show that increasing the value of the Boussinesq coefficient can significantly reduce the numerical oscillations and reduce the predicted area of inundation. In order to stabilise the ADI computations, the Boussinesq coefficient could be judiciously raised or lowered depending on whether the rapid flow is steady or unsteady and whether the bed is frictional or frictionless. An increase in the Boussinesq coefficient generally leads to overprediction of the propagating speed of the flood wave over a frictionless bed, but the opposite is true when bed friction is significant.

Anatomical variations in the suprascapular nerve (SSN) and its depth in the suprascapular notch can make it difficult to target with ultrasonography (US). One alternative could be a proximal approach to the SSN, if US provides a reliable description of its origin (orSSN). The primary objective of this study was to demonstrate that US can reliably locate the orSSN. The secondary objective was to describe the features of the proximal SSN. Seventy brachial plexuses (BPs) from 30 healthy volunteers (60 BPs) and 5 cadavers (10 BPs) were included. There were two parts to this study: (1) description of the proximal SSN in healthy volunteers using US to determine the diameter, depth and location of the orSSN; (2) targeting of the orSSN with US in cadaver limbs to determine its distance from the needle, ink marking and locating the orSSN. In Part I, the diameter of the orSSN averaged 1.33 mm (1-9 mm) and its depth averaged 5.12 mm (2.7-10.6 mm). The orSSN was located in the upper trunk of the BP (53) or its posterior division (7). In Part II, the orSSN was successfully targeted in nine of the 10 specimens by US; the needle/orSSN distance averaged 3.8 mm (0-8 mm). The implanted needle was at the orSSN in two cases, proximal to it in seven and distal to it in one. US is a valid modality for describing and pinpointing the orSSN, irrespective of patient morphology. Clin. Anat. 30:747-752, 2017. © 2017Wiley Periodicals, Inc.

International audience; Avalanche models are increasingly employed for elaborating land-use maps and designing defense structures, but they rely on poorly known parameters. Careful uncertainty assessment is thus required but difficulty arises from the nature of the outputs of these models, which are commonly both functional and scalar. Hence, so far in the avalanche field, few sensitivity analyses have been performed. In this work, we propose to determine the most influential inputs of an avalanche model by estimating aggregated Sobol' indices. We propose a nonparametric estimation procedure based on the Nadaraya-Watson kernel smoother, which allows to estimate the aggregated Sobol' indices from a given random sample of small to moderate size. Due to the limited size of the sample, the kernel estimation is biased. Therefore , we propose a bootstrap based bias correction before selecting the bandwidth by cross-validation. After different test-cases showing the efficiency of our approach, it is applied to a real avalanche case. Results show that the friction parameters and the snow depth in the release zone are the most influential parameters determining the avalanche characteristics.

The complexity of transfer of momentum between the main channel and flood plain in a compound channel can be tackled through the analytical solution of Shiono and Knight Method. Solution of this approach counts on the calibration of three indispensable parameters i.e., secondary currents, eddy viscosity coefficient and bottom friction. Due to uniqueness of these parameters, the boundary friction factor changes unevenly throughout the lateral direction especially at the shear layer region. In this paper, the analytical solution eventuating in the Shiono Knight Method is solved considering a new panel of shear layer width at the junction. Reliable experimental measurements are employed for quantification of shear layer width and to procure the new calibrating coefficients for secondary flow and friction factor for the proposed panels. The model is found to provide satisfactory results when applied to experimental, FCF channels and natural river data sets.

Irregular, unauthorized and non-technical mining of sand and gravel from rivers play an important role in undesirable changes in morphology and environment. The present experimental study investigated sand mining pit modelling in the different scenarios. A two-dimensional flexible mesh, depth averaged model was set-up to simulate the experiments. The results exhibited an acceptable agreement between the numerical and experimental results, so that the maximum error was limited to less than 10%. After verifying the numerical model, a critical reach in Helleh River was selected as a case study. The effect of sand and gravel mining on a bridge at the site was studied. The results also indicated that for a flood event with a return period of 25 years, sand and gravel mining changed the bed profile up to 1.27 m at the bridge foundation. Moreover, initial signs of river meandering emerged in one scenario.

A two-dimensional depth-averaged continuum framework is developed for simulating the interaction of flows outside and within porous media. A unique set of generalized equations applicable to both surface and subsurface systems is derived originating from conservation laws of microscopic variables with first volume-averaging and then integrating vertically. The model is free of ad-hoc coupling as interface conditions are satisfied implicitly owing to the generalized nature of the governing equations. The set of depth-integrated equations (unconditionally hyperbolic) is resolved by employing a Total-Variation-Diminishing MacCormack scheme. The developed model is validated with four one-dimensional experimental test-cases (containing both Darcy and non-Darcy regimes) and applied to a demonstrative two-dimensional scenario having various interface conditions. Results obtained from the proposed framework highlights the potential applicability of the generalized model for diverse time-scales and vertical/ sloping interface conditions.

Based on the Reynolds-averaged Navier–Stokes equations and the Shiono and Knight method, this paper describes the development of an analytical model to predict the lateral distribution of d...

This study builds upon established methods for using acoustic Doppler current profilers (ADCPs) to measure velocities to include another measurement technique—longitudinal measurements—whic...

Simultaneous effects of centrifugal force, floodplain interferences, and variation of the depth-averaged velocity are combined in an analytical model to derive a relationship for estimating...

Objective The present pilot study evaluated the impact of closed flap surgical technique on periodontal tissues and the esthetic outcome. Materials and methods Twenty-four consecutively treated patients with a mean age of 15 (range 12-30) underwent closed flap surgery for the exposure of impacted teeth. The study evaluated periodontal probing depths, bone level, zone of keratinized tissue, and esthetic appearance after the completion of orthodontic treatment. Results Tooth impaction was more frequent in females (No. 14) than males (No. 10) in this study. Palatal and buccal impactions were 77% and 23%, respectively. Post-orthodontic treatment, periodontal probing depth averaged 2.2 mm and zone of keratinized gingiva averaged 3.6 mm. No discernible radiographic bone loss or root resorption was noticed. Overall, gingival margins provided an esthetic appearance in all cases. Conclusion Closed flap surgical technique was associated with periodontal health and esthetic outcomes for an average of 10 years post-orthodontic treatment. Clinical significance The evolution of teeth and maxillary/mandibular bones has distinctively differed in origins and development contributing to tooth impaction and overcrowding. This study demonstrated that closed flap surgery is a reliable procedure for the treatment of impacted teeth that provides periodontal health and esthetic appearance post-orthodontic treatment.

천이류와 같은 급변류에 의한 하상변동을 예측하기 위한 이차정확도의 유한체적법 모형을 제시하였다. 부정류 조건하에서의 유사이송과 하상변동문제에 적용하기 위하여 유사이송모형을 천수방정식과 연계하였다. 지배방정식은 MUSCL 기반의 유한체적법을 이용하였고, 계산요소간 흐름률은 HLLC approximate Riemann solver를 이용하여 계산하였다. 일차원과 이차원 수로에서의 댐붕괴파에 의한 하상변동문제와 월류로 인한 하류부 댐사면의 침식문제에 적용한 결과, 적정한 매개변수를 이용하는 경우에 전반적으로 정확한 수치모의 결과를 얻을 수 있었다. 또한 전반적인 계산결과는 수치적으로 안정적이고 물리적으로 타당한 결과를 나타내었고, 이로부터 제시된 수치모의 기법이 상류와 사류조건하에서의 하상변동 문제에 적용이 가능할 것으로 판단된다. 【A stable second-order finite volume method was proposed to predict sediment transport under rapidly varied flow conditions such as transcritical flow. For the use under unsteady flow conditions, a sediment transport model was coupled with shallow water equations. HLLC approximate Riemann solver based on a monotone upstream-centered schemes for conservation laws (MUSCL) reconstruction was used for the computation of the flux terms. From the comparisons of dam break flow experiments on erodible beds in one- and two-dimensional channels, good agreements were obtained when proper parameters were provided. Lastly, dam surface erosion problem by overtopped water was simulated. Overall, the numerical solutions showed reasonable results, which demonstrated that the proposed numerical scheme could provide stable and physical results in the cases of subcritical and supercritical flow conditions.】

Lateral distributions of depth-averaged velocity in open compound channels with submerged vegetated floodplains are analyzed, based on an analytical solution to the depth-integrated Reynolds-Averaged Navier-Stokes equation with a term included to account for the effects of vegetation. The cases of open channels are: rectangular channel with submerged vegetated corner, and compound channel with submerged vegetated floodplain. The present paper proposes a method for predicting lateral distribution of the depth-averaged velocity with submerged vegetated floodplains. The method is based on a two-layer approach where flow above and through the vegetation layer is described separately. An experiment in compound channel with submerged vegetated floodplain is carried out for the present research. The analytical solutions of the three cases are compared with experimental data. The corresponding analytical depth-averaged velocity distributions show good agreement with the experimental data.

Many geophysical problems and engineering applications involve fluid flow through evolving faults and fractures. These situations are typical in CO2 injection scenarios where supercritical CO2 or CO2–saturated fluids are able to dissolve the solid matrix. Dissolution or precipitation within a fracture could lead to self-reinforcing or self-limiting behavior due to feedbacks from changes in permeability (k). There is limited predictive understanding of this phenomenon. Here, we investigated the case of fracture dissolution using a two-dimensional depth-averaged reactive transport model. The model simulated an acidic fluid dissolving a three-dimensional rough fracture. We observed different dissolution patterns across a wide range of geometric properties and different Peclet and Damkohler numbers (960 cases total). Despite this variation, a persistent linear increase in k with time (t) roughly manifested once the dissolution front advanced through the fracture. The time rate of change of k (k/t) is roughly constant for all simulations. The modeled rate is consistent with a theoretical value based on the parallel plates model and with experimentally derived values. The results confirm the common presence of linear time-dependence of fracture k due to self-reinforcement.

Natural rivers have many branching junctions. The flow in branching junctions is complex, owing to significant changes associated with flow dynamics and sediment transport that result in erosion and deposition problems. A branching channel of the Tigris River in Missan, Iraq, was selected for investigation of the scouring and deposition zones. A two-dimensional (2D) numerical model was used to simulate the hydro-morphodynamics in the branching channel, where hypothetical vanes as control structures were included at the junction to control the scouring and deposition zones. The simulation results suggest the most effective location, dimension, and angle of the introduced vanes. For the studied junction, controlling morphological features was achieved by introducing a single vane with an inclination angle of 90&deg, on the flow direction of the Tigris River. The most effective location of the introduced vane was the location that caused considerable enhancement in the flow depth and velocity distribution.

This paper presents a practical method to predict lateral depth-averaged velocity distribution in trapezoidal meandering channels. Flow structure in meandering channels is more complex than straight channels due to three-dimentional nature of flow. Continuous variation of channel geometry along the flow path associated with secondary currents makes the depth-averaged velocity computation difficult. Design methods based on straight-wide channels incorporate large errors while estimating discharge in meandering channel. Hence, based on the present experimental results, a non-linear form of equation involving three parameters for estimating lateral depth-averaged velocity is formulated. The present experimental meandering channel is wide (aspect ratio = b/h > 5) and with high sinuosity of 2.04. The study serves for a better understanding of the flow and velocity patterns in trapezoidal meandering channel.

Smoothed Particle Hydrodynamics has been successfully used for various fluid-dynamics problems, such as breaking-waves, flooding etc., since it was originally proposed. While the Lagrangian approach is naturally suitable for free-surface flows, enforcing boundary conditions and poor approximations in the presence of discontinuities in the solution are major difficulties with the method. In this paper we present an enhanced conservative Godunov SPH based on the work of Inutsuka [S. Inutsuka, Reformulation of smoothed particle hydrodynamics with Riemann solver, Journal of Computational Physics 179 (2002) 238–267] that accurately resolves discontinuities without the need to use artificial viscosity, preserves partition of unity everywhere in the domain, correctly and flexibly enforces necessary essential and frictional slip boundary conditions to approximately solve free-surface granular flows. The development is motivated by the need to improve upon depth averaged grid based models of large scale debris flows and avalanches often characterized as granular flows. Simple validation of the results is obtained by comparison to table-top experiments.