Chapter 3 An introduction to OpenFOAM

Chapter 3 is an introduction to OpenFOAM, as the most popular CFD tool in effluent discharge modelling. A decade ago, OpenFOAM was only an academic tool for studying complex fluid mechanics problems. However, it has grown quickly in recent years and has found its way into industry as well (e.g., many consulting firms have invested in creating their own solvers for the particular problems they often solve). This book covers the fundamentals of OpenFOAM related to effluent discharge modeling: the choice of available solvers and differences between them, mesh generation options and methodology in OpenFOAM and postprocessing the numerical results

Improving the Accuracy of Hydrodynamic Simulations in Data Scarce Environments Using Bayesian Model Averaging: A Case Study of the Inner Niger Delta, Mali, West Africa

In this paper, the study area was the Inner Niger Delta (IND) in Mali, West Africa. The IND is threatened by climate change, increasing irrigation, and dam operations. 2D hydrodynamic modelling was used to simulate water levels, discharge, and inundation extent in the IND. Three different digital elevation models (DEM) (SRTM, MERIT, and a DEM derived from satellite images were used as a source of elevation data. Six different models were created, with different sources of elevation data and different downstream boundary conditions. Given that the performance of the models varies according to the location in the IND, the variable under consideration and the performance criteria, Bayesian Model Averaging (BMA) was used to assess the relative performance of each of the six models. The BMA weights, along with deterministic performance measures, such as the Nash Sutcliffe coefficient (NS) and the Pearson’s correlation coefficient (r), provide quantitative evidence as to which model is the best when simulating a particular hydraulic variable at a particular location. After the models were combined with BMA, both discharge and water levels could be simulated with reasonable precision (NS > 0.8). The results of this work can contribute to the more efficient management of water resources in the IND

Conservative characteristic-based schemes for shallow flows

Tableau d'honneur de la Faculté des études supérieures et postdoctorales, 2005-2006Les équations en eaux peu profondes, encore appelées équations de Saint-Venant, sont utilisées dans de nombreux cas importants comme les fleuves, les lacs, les estuaires et les océans. La conservation de certaines quantités est une propriété importante qui est habituellement désirée pour assurer la précision des simulations à long terme et également pour le cas des écoulements complexes avec présence d'ondes de choc. Cette thèse examine tout d'abord la formulation de schémas semi-Lagrangiens, qui sont bien connus pour demeurer stables pour des nombres très élevés de CFL. Cependant, ces schémas perdent leur propriété de stabilité lorsque la conservation totale des quantités, qui est cruciale pour une simulation correcte les ondes de chocs, est imposée. Un schéma semi- Lagrangien entièrement conservatif est développé ici et ce dernier demeure stable pour des nombres élevés de CFL. L'approche proposée est ensuite étendue à la méthode des caractéristiques (MOC) et une version conservative du schéma MOC est développée. Contrairement au schéma MOC original, qui ne peut pas simuler correctement les ondes de choc à cause du manque de conservation, le schéma proposé les simule avec succès. De plus, le nouveau schéma présente des avantages sur le plan numérique, tant pour la diffusion et la dispersion que pour la stabilité. Le cas 2D est ensuite considéré, et la méthode de volume finie est utilisée à cause de son conservation inhérente. Le cas 2D est ensuite considéré, et la méthode de volumes finis est utilisée à cause de ses qualités inhérentes de conservation. La plupart des méthodes numériques disponibles sont sensibles au problème du déséquilibre entre les termes source et de flux, particulièrement en présence d'un maillage non structuré. D'autre part, la plupart des schémas numériques disponibles (par exemple les schémas HLL et ENO) induisent un niveau élevé de diffusion numérique en simulant des écoulements tourbillonnaires. Trois approches différentes, applicables sur des maillages non structurés sont développées ici. Elles peuvent simuler des conditions complexes d'écoulement comprenant les topographies variables, les écoulements tourbillonnaires, trans-critiques et discontinus. Finalement plusieurs méthodes de volumes finis upwind sont utilisées, via une analyse de type Fourier, pour évaluer le niveau d`amortissement des modes de Rossby. Contrairement aux bons résultats habituellement obtenus par les méthodes de volumes finis upwind dans iii le cas d'écoulements dominés par la convection, on remarque ici que les ondes de Rossby sont amorties de manière excessive.Shallow water equations arise in many important cases such as in rivers, lakes, estuaries and oceans. Conservation is an important property which is usually desired to ensure the accuracy of the long term simulations and also for the case of complex flows with shockwaves. This thesis begins with semi-Lagrangian schemes, which are well known to remain stable for very high CFL numbers. However, they lose their high stability property when the fully conservative property, which is crucial for a correct simulation of shock waves, is imposed. An inherently fully conservative semi-Lagrangian scheme is developed here which remains stable for high CFL numbers. The proposed approach is then extended to the method of characteristics (MOC) and a conservative extension of MOC is developed. Contrary to the original MOC, which is unable to simulate shockwaves due to the lack of conservation, the proposed scheme easily simulates them. Further, the new scheme presents favorable features in terms of numerical diffusion and dispersion. The 2D case is then considered, and the finite volume method is employed due to its inherent conservation properties. Most available numerical methods face the problem of imbalance between the source and flux terms, particularly when unstructured grids are used. On the other hand, most available numerical schemes (such as the HLL and the ENO schemes) induce a high level of numerical diffusion in simulating recirculating flows. Three different approaches using unstructured grids are successfully developed here. The new schemes can simulate complex flow conditions including recirculating, trans-critical and discontinuous flows over variable topographies. Finally, the performance of the upwind finite volume schemes, for Rossby waves, is studied using a Fourier analysis approach. Contrary to the usual good results obtained for those schemes in the case of convection dominated flows, it is observed here that they lead to an excessive damping of the Rossby modes

Buoyant Jets in Cross-Flows: Review, Developments, and Applications

Significant environmental effects from the use of marine outfall discharges have led to increased efforts by both regulatory bodies and research groups to minimize the negative impacts of discharges on the receiving water bodies. Understanding the characteristics of discharges under conditions representative of marine environments can enhance the management of discharges and mitigate the adverse impacts to marine biota. Thus, special attention should be given to ambient cross-flow effects on the mixing behaviors of jet discharges. A buoyant jet in cross-flow has different practical applications such as film cooling and dilution, and provide a higher mixing capability in comparison with free jets or discharges into stationary environments. The main reason for this is believed to be the existence of various complicated vortical structures including a counter-rotating vortex pair as the jet expands downstream. Although tremendous research efforts have been devoted to buoyant jets issuing into cross-flows over the past five decades, the mixing process of an effluent at the discharge point is not yet well understood because of the highly complex fluid interactions and dispersion patterns involved. Therefore, there is a need for a deeper understanding of buoyant jets in cross-flows in order to obtain better predictive methods and more accurate design guidelines. The main aims of this study were (i) to establish the background behind the subject of buoyant jets in cross-flows including the flow structures resulting from the interaction of jets and cross-flows and the impacts of current on mixing and transport behavior; (ii) to present a summary of relevant experimental and numerical research efforts; and finally, (iii) to identify and discuss research gaps and future research directions

Selenium Removal from Water and Wastewater by Different Technologies: A Systematic Review

Background: Selenium (Se) is an essential element playing a vital role in the metabolism of organisms. Se can generally be discharged in the potable water through natural and anthropogenic activities. Both excess and shortage of Se can cause significant adverse health effects in humans. Excess values of se may stimulate toxici-ty, leading to selenosis and alkali disease in humans and grazing animals, respectively. Methods: A review search was systematically carried out from the databases Embase, PubMed/MEDLINE, Scopus, PubMed Central (PMC), Google Scholar, as well as medRxiv by using the following keywords: “waste water”, “bioremediation”, “selenium removal”, “adsorption”, and “drinking water”. This study provides a review of the recent literature covering the period between 2011 and 2021. After screening the full text of the articles, 27 papers were enrolled. This study reviews the reported techniques for Se removal from water and wastewater, including adsorption, biological treatment, microbial reduction, bioreactors, fungal bioreactor, algal treatment, phytoremediation, and photocatalysis. Results: Biological and bioremediation techniques, such as microbial reduction, biotransformation, and fluid-ized bed reactor have removal efficiency about 100. The highest Se concentration of 15-7600 µg/L was achieved in ground waters in Ethiopia and the lowest level of 0.07 µg /L in Finland. Conclusion: The combination of biological treatment with chemical or physical technologies is envisaged to optimize se elimination and to ensure ecological protection and human health safety

Numerical Simulation of Flow in Parshall Flume Using Selected Nonlinear Turbulence Models

This study uses a computational fluid dynamics (CFD) approach to simulate flows in Parshall flumes, which are used to measure flowrates in channels. The numerical results are compared with the experimental data, which show that choosing the right turbulence model, e.g., v2−f and LC, is the key element in accurately simulating Parshall flumes. The Standard Error of Estimate (SEE) values were very low, i.e., 0.76% and 1.00%, respectively, for the two models mentioned above. The Parshall flume used for this experiment is a good example of a hydraulic structure for which the design can be more improved by implementing a CFD approach compared with a laboratory (physical) modeling approach, which is often costly and time-consuming

A Review of Numerical Simulations of Secondary Flows in River Bends

River bends are one of the common elements in most natural rivers, and secondary flow is one of the most important flow features in the bends. The secondary flow is perpendicular to the main flow and has a helical path moving towards the outer bank at the upper part of the river cross-section, and towards the inner bank at the lower part of the river cross-section. The secondary flow causes a redistribution in the main flow. Accordingly, this redistribution and sediment transport by the secondary flow may lead to the formation of a typical pattern of river bend profile. It is important to study and understand the flow pattern in order to predict the profile and the position of the bend in the river. However, there are a lack of comprehensive reviews on the advances in numerical modeling of bend secondary flow in the literature. Therefore, this study comprehensively reviews the fundamentals of secondary flow, the governing equations and boundary conditions for numerical simulations, and previous numerical studies on river bend flows. Most importantly, it reviews various numerical simulation strategies and performance of various turbulence models in simulating the flow in river bends and concludes that the main problem is finding the appropriate model for each case of turbulent flow. The present review summarizes the recent advances in numerical modeling of secondary flow and points out the key challenges, which can provide useful information for future studies

Stability analysis of unstructured finite volume methods for linear shallow water flows using pseudospectra and singular value decomposition

The discretization of the shallow water system on unstructured grids can lead to spurious modes which usually can affect accuracy and/or cause stability problems. This paper introduces a new approach for stability analysis of unstructured linear finite volume schemes for linear shallow water equations with the Coriolis Effect using spectra, pseudospectra, and singular value decomposition. The discrete operator of the scheme is the principal parameter used in the analysis. It is shown that unstructured grids have a large influence on operator normality. In some cases the eigenvectors of the operator can be far from orthogonal, which leads to amplification of solutions and/or stability problems. Large amplifications of the solution can be observed, even for discrete operators which respect the condition of asymptotic stability, and in some cases even for Lax-Richtmyer stable methods. The pseudospectra are shown to be efficient for the verification of stability of finite volume methods for linear shallow water equations. In some cases, the singular value decomposition is employed for further analysis in order to provide more information about the existence of unstable modes. The results of the analysis can be helpful in choosing the type of mesh, the appropriate placements of the variables of the system on the grid, and the suitable discretization method which is stable for a wide range of modes.The authors thank the anonymous reviewers for their valuable comments and suggestions to improve the quality of the paper. This publication was made possible by NPRP Grant 4-935-2-354 from the Qatar National Research Fund (a member of the Qatar Foundation). The statements made herein are solely the responsibility of the authors. Appendix AScopu