Abstract. 1. Introduction
|
|
- Virginia Phelps
- 6 years ago
- Views:
Transcription
1 Dispersion of pollutants in coastal waters near a river mouth: a simplified numerical model M. Di Natale Dipartimento dilngegneria Civile, Seconda Universita dinapoli, Aversa (CE), Italy Abstract A mathematical model is presented for studying the hydrodynamic dispersion of pollutants discharged in coastal areas near a river mouth during flooding. The hydrodynamic process is represented in a 2D depth integrated system and the solute is assumed non-reactive. The model is numerically solved by finite differences characterized by an appropriate algorithm which is explicit to describe the field of motion and implicit to integrate the transport equation. 1. Introduction Within the framework of pollution problems in coastal waters, areas near river mouths assume particular technical and scientific interest. Due to the possible presence of many pollution sources in the hydrographic basin (eg. domestic sewage, toxic substances from industrial discharges, surface runoff of pesticides and fertilisers from agricultural land etc.), water-courses become carriers of pollutants which, without appropriate treatment plants, are transported to the mouth and then out to sea. The study of such problems is of considerable importance so as to define the extent of the polluted stretch of water and to ascertain the environmental impact upon bordering coastal zones. In some cases, dispersion phenomena in the marine ambient may be enhanced by wave motion or coastal currents. Thus, the most serious conditions may occur when the sea is calm. Moreover, it should be noted that the transport of pollutants is strictly correlated to variation of river flows in time, particularly during flooding. As is well known, during the initial phase of a flood, pollutant concentrations increase due to the washing effect of the first rains on sealed areas in urban zones, on vegetation and agricultural land. The
2 156 Environmental Problems in Coastal Regions subsequent phases, during the rising limb and the recession limb of the hydrograph, are characterized by a reduction in the polluting load. In this paper, a mathematical model is presented to describe the above phenomenon under the hypothesis of a calm sea into which flows a flood wave transporting a polluting load with a time-variable concentration. Hydrodynamic analysis of the problem is characterized by two phases: in thefirstphase, we obtain the equations of the field of motion in the coastal water stretch deriving from the in-flow of the flood Q(t); subsequently, using the equation of advective-diffusive transport, we examine hydrodynamic dispersion of time-variable concentration C(t) in the above-mentioned field of motion. The proposed model is not constrained either by the morphology of the mouth and the sea-bottom, or by Q(t) and C(t) functions which may be completely arbitrary. The introduced approximations concern the hypothesis of shallow water conditions, in which the field of motion and the concentrations refer to a 2D depth integrated scheme. Moreover, buoyancy effects deriving from small density differences betweenfluvialand marine water were ignored. The pollutants were assumed non-reactive in the sense that they do not modify the field of motion, although their chemical and biological transformation in time is made possible by introducing an appropriate decay term in the transport equation. The mathematical model is solved numerically and an application is reported to show its validity in a real case. 2. Hydrodynamic formulation With reference to the scheme in fig. 1, in the most general case of a 3- dimensional field of motion, the equations which describe the hydrodynamic process are [1]: which expresses the continuity condition, and : du dv dw + +%- = o (i) dx dy dz du _gu _gu gu Idp d( du} d du d( du ^^= (2-a) dv dv _5v 0v 1 dp d ( dv} d dv d ( dv (2-b) dv/ dv/ dv/ dv/ 1 dp d ( dv/} d ( dv/\ d ( dvf m} -- i-u +v +w = - + c.. H--g. + e at dx dy dz g dz M * dxj dy{ * dy ) dz{ * fr) dz (2-c)
3 Environmental Problems in Coastal Regions 157 which indicate the dynamic equilibrium. The terms u, v, w and p represent, respectively, the time means of the velocity components and the pressure. By using the Boussinesq's approximation to describe the turbulent stresses, E^ and e^ are horizontal and vertical eddy viscosity coefficients. Supposing L%, Ly and L% are, respectively, the characteristic scales of the lengths in directions x, y and z, in the case in question we obtain [2]: O[LJ»O[L,} O[L,]>O[LJ (3) from which it follows that w (x,y,z,t) = 0 and hence eqn. (2c) becomes: ar-«* (4) which characterizes a hydrostatic distribution of the pressures. z A Figure 1: Hydrodynamic scheme By integrating eqns. (2a-b) along the vertical and introducing the mean velocity U = J udz and V = j vdz, being h = hg -ho -ho (e.g. with reference to component U):, we also obtain
4 158 Environmental Problems in Coastal Regions 3x1 3x 0[LJ o -k dy( * dy Thus, from eqn. (3), we have the following: He.,^-»0 OL (5-a) O a f au a ( au ax (5-b) As regards the terms of resistance, it should be noted that, integrating along the vertical, it is: ^ r r -HQ a av p p TS and ty being the shear stresses on the surface and on the bottom. In practical applications, it is commonly considered that: and TS and %,, can be expressed as [3]: (6) (7-a) (7-b) in which p is the water density and f the non-dimensional friction coefficient which is assumed [4] (f = 1(H -* 10-^)
5 Environmental Problems in Coastal Regions 159 Finally, taking into consideration (4), (5a-b), (6) and (7a-b), the equations of motion and continuity in the 2D depth integrated model are as follows [5], [6]: I bx t/ \-/ 6V gn T_ g^y c^ <9Uh -^+ + = 0 gt gx 9y (io) 3. Mathematical model for turbulent advective-diffusive pollution transport Let C(x,y,t) be the instantaneous concentration (defined as the pollutant mass per unit of volume of the fluid). The solute mass conservation equation is [7]: dc Eqn. (11) obviously refers to the above 2D-depth integrated model and thus the C values (x,y,t) represent the means along the depth. In (11) the molecular diffusion effect was ignored since it is very slight compared with the effect of turbulence [8]. The terms Kx and Ky represent the turbulent dispersion coefficients while 9 (dimension time'*) is a coefficient which characterizes the time decay of the pollutant. Kx and Ky values, which represent one of the most difficult unknowns to estimate, are commonly obtained through the distribution assigned to the velocity profiles u(x,y,z,t) and v(x,y,z,t) and through the eddy diffusivity coefficients D% and Dy which, in accordance with Boussinesq's approximation, describe the diffusive terms in the most complete 3D form of eqn. 01). It is also worth noting that the differential eqn. (11) is a mixed type in the sense that the terms of the left-hand member (hyperbolic part) represent advective transport while the terms of the right-hand member (parabolic part) describe the diffusion effects. The relative importance between advective and diffusion terms is highlighted by the Pec let number: Pe = U;L/K in which U«,, L and K represent, respectively, the characteristic velocity of the fluid, the characteristic length of the flow domain and the dispersion coefficient. When Pe > 1, the advection effect prevails over diffusion and vice versa.
6 160 Environmental Problems in Coastal Regions 4. Numerical solution The solution of the PDE (9 a-b) and (10), which describe the field of motion and the solution of the transport differential equation (11) is obtained numerically by using finite difference methods that are particularly straightforward and efficient for the case in question [9], [10], [11]. In particular for eqns. (9 a-b) and (10), an explicit form is used. The flow domain is characterized by an orthogonal staggered grid for space-time discretizations with mesh sides Ax, Ay (fig. 2) and time step At. The calculation points are fixed by the i, j, m indices: the first and second for spatial coordinates x, y, the third for time t. The finite difference forms, according to the explicit scheme presented, are [12]: gat ( -4 "41 2f (12-a) vr'=vr- ^ (12-b) (12-c) where:
7 Environmental Problems in Coastal Regions 161 Ay i. J+1 Ax 4- Figure 2: Orthogonal staggered grid for spatial discretization The boundary conditions are fixed as follows: - the velocities at the mouth are known V" = VQ ; - the velocity components orthogonal to the shoreline are zero VJ" = 0; -in the open sea U^,V^,^ = U?,V^ along x-direction and U"u, V y,t y = Uf,V»,1$ alongy^lirection. Due to the explicit form of the integration, the space steps Ax, Ay and the time steps At are fixed so as to verify the stability condition CFL [13]: w< jj being water depth in the generic mesh i, j.
8 162 Environmental Problems in Coastal Regions As regards the transport eqn. (11), it should be pointed out that, as is well known, explicit numerical finite difference methods generate significant errors due to the effect of numerical diffusion. The proposed calculation method is therefore implicit with centered differences and is always stable for each Ax, Ay and At value. In such applications, with a view to reducing the dimensions of the arrays, simplified algorithms are used such as the ADI method. This method alternately involves at time step m an implicit scheme along x and an explicit scheme along y, while at time step m+1 an implicit scheme along y and an explicit scheme along x. The solution proposed here is fully implicit and, with reference to fig. 2, the inductive relation is [12]: 'ym+l pm+1 _ ym+\ pm+1 Mj+1 ^ij+1 Mj-1 ^ij-1 At 2Ax 2Ay K x Ay' 03) The system of linear algebraic equations supplied by (13) is characterized by a large, sparse, unsymmetric array. To reduce calculation times, use is made of a particular simplifying algorithm [14] which involves the transformation of the initial arrays into two packed arrays, one to memorize the non-zero elements and the other to fix the real positions of their indices. Boundary conditions are fixed as follows: - concentration values are known at the mouth: C = CQ ; - the concentration flux in the direction orthogonal to the shoreline is null: 8C ^ * - in the open sea the diffusive transport component is null, that is: g^c <^C r = 0 along y-direction and = 0 along x-direction. ox dy 5. Applicative example The results obtained from the application of the proposed numerical model in an applicative example are reported below. The geometric scheme in question is that infig.2, where it is assumed that a = 50 m, b = 600 m and c=1000 m. The water depth at therivermouth is 5 m and the sea-bottom has a uniform 1% slope increasing in the offshore direction. The flood hydrograph and the corresponding time variation of concentrations CQ at the mouth are reported in the histograms in fig. 3 a-b. The following assumptions are made: Ax = Ay = 25 m; moreover, it is considered that K = Kx = Ky and two different values of the dispersion coefficients are chosen, K = 1 nf/s and K = 10 nf/s, in
9 Environmental Problems in Coastal Regions 163 order to highlight the effect of the diffusive component for two different Peclet numbers. L t 3- s 2 - i _ 0 - k ^ t(h) Figure 3-a: Time variations of VQ i? I 1^ '" ^ d, n«0 ' 0,6-0, fe» , J 0 Figure 3-b: Time variations of GO The results are shown in the graphs in fig. 4, 5 and 6 where the spatial patterns of the concentration at fixed values of the time t are reported by means of isoconcentration curves. It should be noted that, at increasing times, the zones with high concentrations move towards the open sea. Such an effect is more evident for the lower K value. Moreover, as is to be expected, for K = 10 m^/s the greater effect of diffusive transport reduces the extent of the areas with high concentrations.
10 164 Environmental Problems in Coastal Regions 6. Conclusions Dispersion of pollutants transported by a river flood to a coastal zone near a mouth may be studied with good approximation by using the proposed mathematical model which is based on a 2D depth integrated representation of the field of motion. The proposedfinitedifference numerical model, which is both reliable and not excessively time-consuming, requires the use of an explicit scheme for defining the flow domain and a fully implicit scheme for estimating the distribution of pollutant concentrations. References 1. Ligget, J.A., Fluid Mechanics, McGraw-Hill, Pedlosky, J., Geophysical Fluid Dynamics, Springer-Verlag, Csanady, G.T., Circulation in the Coastal Ocean, Reidel Pub.Company, Soulsby, R.L., Tidal current boundary layers, Ocean Eng. Science vol.9-a, Abbott, M.B.,Price, W.A.,Coastal Estuarial and Harbour Engineers, Reference Book,FN SPON, Di Natale, M., A 2D stochastic modelfor the hydrodynamic dispersion of pollutants on the sea surface, Second International Conference on Computer Modelling of Seas and Coastal Region II, Cancun-Mexico, Holly, P.M. Jr. Dispersion in river and coastal waters-physical principles, Developments in Hydraulic Engineering 3, Elsevier Pub. New York, Fisher, H.B.,List, E.J., Koh, R.C.Y. Imberger, J.,Brooks, N.H., Mixing in Inland and Coastal Waters, Academic Press.Inc., Sauvaget, P., Dispersion in river and coastal -waters-numerical simulation, Developments in Hydraulic Engineering 3, Elsevier Pub. New York, Kowalik, Z.,Murty, T.S., Numerical Modelling of Ocean Dynamics,Advanced Series on Ocean Engineering, World Scientific Pub., Brebbia, C.A., Computational Hydrodynamics, Butterworths, Koutitas, C.G. Mathematical Models in Coastal Engineeringfeatech Press, Smith, G., Numerical Solution ofpartial Differencial Equations, Oxford, Gupta, S.K., Tanji K.K., Computer program for solution of large, sparse, unsymmetric linear systems of equations, International Journal for Numerical Methods in Engineering, vol. 11, 1987
11 Environmental Problems in Coastal Regions 165 Figure 4 - Visualization of the isoconcentration lines at t = 1 h a) values of C for k = 10 m2/s b) values of C for k = 1 m2/s
12 166 Environmental Problems in Coastal Regions Figure 5 - Visualization of the isoconcentration lines at t = 3 h a) values of C for k = 10 m2/s b) values of C for k = 1 m2/s
13 Environmental Problems in Coastal Regions 167 Figure 6 - Visualization of the isoconcentration lines at t = 6 h a) values of C for k = 10 m2/s b) values of C for k= 1 m2/s
3.3 Classification Diagrams Estuarine Zone Coastal Lagoons References Physical Properties and Experiments in
Contents 1 Introduction to Estuary Studies... 1 1.1 Why to Study Estuaries?.... 1 1.2 Origin and Geological Age... 4 1.3 Definition and Terminology... 7 1.4 Policy and Actions to Estuary Preservation....
More informationIntroduction to BASEMENT Basic Simulation Environment for Computation of Environmental Flow and Natural Hazard Simulation
Introduction to BASEMENT Basic Simulation Environment for Computation of Environmental Flow and Natural Hazard Simulation Numerical Hydraulics Autumn semester 2016 Prof. Dr. Markus Holzner Author: Pascal
More informationModels of ocean circulation are all based on the equations of motion.
Equations of motion Models of ocean circulation are all based on the equations of motion. Only in simple cases the equations of motion can be solved analytically, usually they must be solved numerically.
More informationAbstract. Introduction
Stability Criterion for Explicit Schemes (Finite-Difference Method) on the Solution of the Advection-Diffusion Equation L.F. León, P.M. Austria Mexican Institute of Water Technology, Paseo Cuauhnáhuac
More information1. INTRODUCTION TO CFD SPRING 2018
1. INTRODUCTION TO CFD SPRING 018 1.1 What is computational fluid dynamics? 1. Basic principles of CFD 1.3 Stages in a CFD simulation 1.4 Fluid-flow equations 1.5 The main discretisation methods Appendices
More informationGeomorphological Modelling in Coastal Waters
Abstract Geomorphological Modelling in Coastal Waters Morteza Kolahdoozan 1, Roger A. Falconer 2 (Fellow), Yiping Chen 3 Details are given herein of the development and application of a three dimensional
More informationEffect of turbulence on tidal suspended transport. B.A. O'Connor Department of Civil Engineering, University of Liverpool, UK
Effect of turbulence on tidal suspended transport B.A. O'Connor Department of Civil Engineering, University of Liverpool, UK Abstract The effect of enhanced turbulence upon tidal suspended transport has
More informationCalculation and Analysis of Momentum Coefficients in the Changjiang River Estuary
Third Chinese-German Joint Symposium on Coastal and Ocean Engineering National Cheng Kung University, Tainan November 8-16, 2006 Calculation and Analysis of Momentum Coefficients in the Changjiang River
More information1. INTRODUCTION TO CFD SPRING 2019
1. INTRODUCTION TO CFD SPRING 2019 1.1 What is computational fluid dynamics? 1.2 Basic principles of CFD 1.3 Stages in a CFD simulation 1.4 Fluid-flow equations 1.5 The main discretisation methods Appendices
More informationSecuring Manoeuverability of a Deep Draft Ship in a Sediment loaded Tidal River Berth
Securing Manoeuverability of a Deep Draft Ship in a Sediment loaded Tidal River Berth O. Stoschek 1, A. Matheja 1 & C. Zimmermann 1 1 Franzius-Institute for Hydraulic, Waterways and Coastal Engineering,
More informationSection 11.1 What is a Differential Equation?
1 Section 11.1 What is a Differential Equation? Example 1 Suppose a ball is dropped from the top of a building of height 50 meters. Let h(t) denote the height of the ball after t seconds, then it is known
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationTHREE-DIMENSIONAL FINITE DIFFERENCE MODEL FOR TRANSPORT OF CONSERVATIVE POLLUTANTS
Pergamon Ocean Engng, Vol. 25, No. 6, pp. 425 442, 1998 1998 Elsevier Science Ltd. All rights reserved Printed in Great Britain 0029 8018/98 $19.00 + 0.00 PII: S0029 8018(97)00008 5 THREE-DIMENSIONAL FINITE
More informationOCEAN HYDRODYNAMIC MODEL
Jurnal Teknologi Pengelolaan Limbah (Journal of Waste Management Technology), ISSN 1410-9565 Volume 10 Nomor 1 Juli 2007 (Volume 10, Number 1, July, 2007) Pusat Teknologi Limbah Radioaktif (Radioactive
More informationWQMAP (Water Quality Mapping and Analysis Program) is a proprietary. modeling system developed by Applied Science Associates, Inc.
Appendix A. ASA s WQMAP WQMAP (Water Quality Mapping and Analysis Program) is a proprietary modeling system developed by Applied Science Associates, Inc. and the University of Rhode Island for water quality
More informationAbstract. 1 Introduction
Mass Conservation in shallow water storm surge computations Jeffrey P. Laible Department of Civil and Environmental Engineering University of Vermont Burlington, Vermont 05405 E-mail: laible@emba. uvm.
More informationNumerical Hydraulics
ETH Zurich, Fall 2017 Numerical Hydraulics Assignment 2 Numerical solution of shallow water wave propagation (www.surfertoday.com) 1 Introduction 1.1 Equations Understanding the propagation of shallow
More informationTransactions on Engineering Sciences vol 18, 1998 WIT Press, ISSN
Simulation of natural convection in a reservoir P. Jelmek*, V. Havlik\ R. Cerny\ P. Pfikryl" * Czech Technical University, Faculty of Civil Engineering, Department of Physics, Thdkurova 7, 166 29 Prague
More informationTransactions on Modelling and Simulation vol 10, 1995 WIT Press, ISSN X
A new Crank-Nicholson algorithm for solving the diffusive wave flood routing equation along a complex channel network R. Moussa," C. BouquilW* " Institut National de la Recherche Agronomique, 34060 Montpellier,
More informationThe Shallow Water Equations
The Shallow Water Equations Clint Dawson and Christopher M. Mirabito Institute for Computational Engineering and Sciences University of Texas at Austin clint@ices.utexas.edu September 29, 2008 The Shallow
More informationQUESTIONS FOR THE TESTS-ECI 146
QUESTIONS FOR THE TESTS-ECI 46 IMPORTANT: Please note that the questions in the exams WILL BE DIFFERENT than these. (These questions are valid for the midterm and for the final exam) Remember: A) both
More informationComputational hydraulic techniques for the Saint Venant Equations in arbitrarily shaped geometry
Università degli Studi di Trento Dipartimento di Matematica Dottorato di Ricerca in Matematica Ph.D. Thesis Computational hydraulic techniques for the Saint Venant Equations in arbitrarily shaped geometry
More informationLecture 9: Tidal Rectification, Stratification and Mixing
Lecture 9: Tidal Rectification, Stratification and Mixing Chris Garrett 1 Additional Notes on Tidal Rectification This lecture continues the discussion of long-wavelength tidal flow over comparatively
More informationESCI 485 Air/Sea Interaction Lesson 1 Stresses and Fluxes Dr. DeCaria
ESCI 485 Air/Sea Interaction Lesson 1 Stresses and Fluxes Dr DeCaria References: An Introduction to Dynamic Meteorology, Holton MOMENTUM EQUATIONS The momentum equations governing the ocean or atmosphere
More informationNumerical Simulations of a Stratified Oceanic Bottom Boundary Layer. John R. Taylor - MIT Advisor: Sutanu Sarkar - UCSD
Numerical Simulations of a Stratified Oceanic Bottom Boundary Layer John R. Taylor - MIT Advisor: Sutanu Sarkar - UCSD Motivation Objective I: Assess and improve parameterizations of the bottom boundary
More informationSediment transport and river bed evolution
1 Chapter 1 Sediment transport and river bed evolution 1.1 What is the sediment transport? What is the river bed evolution? System of the interaction between flow and river beds Rivers transport a variety
More informationFluid Mechanics Answer Key of Objective & Conventional Questions
019 MPROVEMENT Mechanical Engineering Fluid Mechanics Answer Key of Objective & Conventional Questions 1 Fluid Properties 1. (c). (b) 3. (c) 4. (576) 5. (3.61)(3.50 to 3.75) 6. (0.058)(0.05 to 0.06) 7.
More informationThe shallow water equations Lecture 8. (photo due to Clark Little /SWNS)
The shallow water equations Lecture 8 (photo due to Clark Little /SWNS) The shallow water equations This lecture: 1) Derive the shallow water equations 2) Their mathematical structure 3) Some consequences
More informationDynamics of the Ems Estuary
Dynamics of the Ems Estuary Physics of coastal systems Jerker Menninga 0439738 Utrecht University Institute for Marine and Atmospheric research Utrecht Lecturer: Prof. dr. H.E. de Swart Abstract During
More informationHydrodynamic Characteristics of Gradually Expanded Channel Flow
International Journal of Hydraulic Engineering 2013, 2(2): 21-27 DOI: 10.5923/j.ijhe.20130202.01 Hydrodynamic Characteristics of Gradually Expanded Channel Flow Edward Ching-Rue y LUO Department of Civil
More informationModeling Complex Flows! Direct Numerical Simulations! Computational Fluid Dynamics!
http://www.nd.edu/~gtryggva/cfd-course/! Modeling Complex Flows! Grétar Tryggvason! Spring 2011! Direct Numerical Simulations! In direct numerical simulations the full unsteady Navier-Stokes equations
More informationContinuum mechanism: Stress and strain
Continuum mechanics deals with the relation between forces (stress, σ) and deformation (strain, ε), or deformation rate (strain rate, ε). Solid materials, rigid, usually deform elastically, that is the
More informationSAMPLE CHAPTERS UNESCO EOLSS WAVES IN THE OCEANS. Wolfgang Fennel Institut für Ostseeforschung Warnemünde (IOW) an der Universität Rostock,Germany
WAVES IN THE OCEANS Wolfgang Fennel Institut für Ostseeforschung Warnemünde (IOW) an der Universität Rostock,Germany Keywords: Wind waves, dispersion, internal waves, inertial oscillations, inertial waves,
More informationDynamics Rotating Tank
Institute for Atmospheric and Climate Science - IACETH Atmospheric Physics Lab Work Dynamics Rotating Tank Large scale flows on different latitudes of the rotating Earth Abstract The large scale atmospheric
More informationMath 266: Ordinary Differential Equations
Math 266: Ordinary Differential Equations Long Jin Purdue University, Spring 2018 Basic information Lectures: MWF 8:30-9:20(111)/9:30-10:20(121), UNIV 103 Instructor: Long Jin (long249@purdue.edu) Office
More informationFluid Mechanics. Jim Emery 9/3/ Bernoulli s Law 2. 2 Torricelli s Law 3. 3 Time to Empty a Tank 3. 4 Viscosity 4. 5 The Acceleration 5
Fluid Mechanics Jim Emery 9/3/2010 Contents 1 Bernoulli s Law 2 2 Torricelli s Law 3 3 Time to Empty a Tank 3 4 Viscosity 4 5 The Acceleration 5 6 The Equation of Motion Due to Cauchy 7 7 The Deformation
More informationAssessment of the performance of a turbulence closure model: along the tidally-influenced Kaipara River to the estuary, NZ
Assessment of the performance of a turbulence closure model: along the tidally-influenced Kaipara River to the estuary, NZ Berengere S. Dejeans 1, Julia C. Mullarney 2, Iain T. MacDonald 3 and Glen M.
More informationExtreme Mixing Events in Rivers
PUBLS. INST. GEOPHYS. POL. ACAD. SC., E-10 (406), 2008 Extreme Mixing Events in Rivers Russell MANSON 1 and Steve WALLIS 2 1 The Richard Stockton College of New Jersey, School of Natural Sciences and Mathematics,
More informationA two-fluid model of turbulent two-phase flow for simulating turbulent stratified flows
Ocean Engineering 30 (2003) 153 161 www.elsevier.com/locate/oceaneng A two-fluid model of turbulent two-phase flow for simulating turbulent stratified flows Y.M. Shen a,, C.-O. Ng b, A.T. Chwang b a State
More informationCOMPARISON BETWEEN EXPERIMENTAL AND NUMERICAL RESULTS OF FLOOD FLOW IN A COMPOUND MEANDERING CHANNEL
COMPARISON BETWEEN EXPERIMENTAL AND NUMERICAL RESULTS OF FLOOD FLOW IN A COMPOUND MEANDERING CHANNEL Alex George Mutasingwa, Graduate student, Hiroshima University Shoji Fukuoka, Professor, Department
More informationCranfield ^91. College of Aeronautics Report No.9007 March The Dry-Bed Problem in Shallow-Water Flows. E F Toro
Cranfield ^91 College of Aeronautics Report No.9007 March 1990 The Dry-Bed Problem in Shallow-Water Flows E F Toro College of Aeronautics Cranfield Institute of Technology Cranfield. Bedford MK43 OAL.
More informationTurbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing.
Turbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing. Thus, it is very important to form both a conceptual understanding and a quantitative
More informationModelling the Channel and the Bay of Biscay using the MARS model. Contributions: Dyneco Physed / Pelagos / Benthos EMH
Modelling the Channel and the Bay of Biscay using the MARS model Contributions: Dyneco Physed / Pelagos / Benthos EMH Overview 1. The MANGA Configuration 2. Model validation: Sea Surface Temperature 3.
More informationPHYSFLU - Physics of Fluids
Coordinating unit: 230 - ETSETB - Barcelona School of Telecommunications Engineering Teaching unit: 748 - FIS - Department of Physics Academic year: Degree: 2018 BACHELOR'S DEGREE IN ENGINEERING PHYSICS
More informationSand Ripple Dynamics on the Inner Shelf
Sand Ripple Dynamics on the Inner Shelf Donald N. Slinn Department of Civil and Coastal Engineering, University of Florida Gainesville, FL 32611-6590, Phone: (352) 392-9537 x 1431 Fax: (352) 392-3466 E-mail:
More informationTransactions on Ecology and the Environment vol 2, 1993 WIT Press, ISSN
Laboratory experiments related to the injection of buoyant fluid layers in stratified flows M. Priven," J. Atkinson/ G.A. Bemporad,' H. Rubin" " CAMERI - Coastal and Marine Engineering Research Institute,
More informationIntroduction to Marine Hydrodynamics
1896 190 1987 006 Introduction to Marine Hydrodynamics (NA35) Department of Naval Architecture and Ocean Engineering School of Naval Architecture, Ocean & Civil Engineering Shanghai Jiao Tong University
More informationESTIMATION OF SUSPENDED SEDIMENT CONCENTRATION IN ESTUARY
1 ESTIMATION OF SUSPENDED SEDIMENT CONCENTRATION IN ESTUARY ZHIYAO SONG,JUN KONG, WEISHENG ZHANG and YUN XING College of Traffic and Ocean Engineering and Eco-environmental Modeling center, Hohai University,
More informationAbstract. 1. Introduction
Flow Characteristics and Flushing Processes in Marinas and Coastal Embayments J.Yin Cranfield University, Cranfield, Bedfordshire MK43 OAL, UK E Mail: J. Yin@Cranfield. ac. ukl R A Falconer Cardiff School
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More information7. Basics of Turbulent Flow Figure 1.
1 7. Basics of Turbulent Flow Whether a flow is laminar or turbulent depends of the relative importance of fluid friction (viscosity) and flow inertia. The ratio of inertial to viscous forces is the Reynolds
More informationValidation of the Hydrodynamic Module of a 3D Ecological Model for the Ria de Aveiro Lagoon (Portugal)
International Congress on Environmental Modelling and Software Brigham Young University BYU ScholarsArchive th International Congress on Environmental Modelling and Software - Barcelona, Catalonia, Spain
More informations and FE X. A. Flow measurement B. properties C. statics D. impulse, and momentum equations E. Pipe and other internal flow 7% of FE Morning Session I
Fundamentals of Engineering (FE) Exam General Section Steven Burian Civil & Environmental Engineering October 26, 2010 s and FE X. A. Flow measurement B. properties C. statics D. impulse, and momentum
More information1 of 7 2/8/2016 4:17 PM
1 of 7 2/8/2016 4:17 PM Continuous issue-4 January - March 2015 Abstract The similarity solution of Diffusion equation using a technique of infinitesimal transformations of groups The similarity solution
More informationOn the validation study devoted to stratified atmospheric flow over an isolated hill
On the validation study devoted to stratified atmospheric flow over an isolated hill Sládek I. 2/, Kozel K. 1/, Jaňour Z. 2/ 1/ U1211, Faculty of Mechanical Engineering, Czech Technical University in Prague.
More informationParts Manual. EPIC II Critical Care Bed REF 2031
EPIC II Critical Care Bed REF 2031 Parts Manual For parts or technical assistance call: USA: 1-800-327-0770 2013/05 B.0 2031-109-006 REV B www.stryker.com Table of Contents English Product Labels... 4
More information1D Heat equation and a finite-difference solver
1D Heat equation and a finite-difference solver Guillaume Riflet MARETEC IST 1 The advection-diffusion equation The original concept applied to a property within a control volume from which is derived
More informationA Study on Residual Flow in the Gulf of Tongking
Journal of Oceanography, Vol. 56, pp. 59 to 68. 2000 A Study on Residual Flow in the Gulf of Tongking DINH-VAN MANH 1 and TETSUO YANAGI 2 1 Department of Civil and Environmental Engineering, Ehime University,
More information17th International Conference on Harmonisation within Atmospheric Dispersion Modelling for Regulatory Purposes 9-12 May 2016, Budapest, Hungary
17th International Conference on Harmonisation within Atmospheric Dispersion Modelling for Regulatory Purposes 9-12 May 2016, Budapest, Hungary SIMULATIONS OF TRAFFIC RELATED POLLUTANTS IN A MAIN STREET
More informationDynamique des rivières. res
Dynamique des rivières res 1 Rivers are enormously diverse size: varies by many orders of magnitude geometry: highly variable substrate: bedrock or sediment sediment type: sediment size ranges from mud
More informationThe estimation of vertical eddy diffusivity in estuaries
The estimation of vertical eddy diffusivity in estuaries A. ~temad-~hahidi' & J. 1rnberger2 l Dept of Civil Eng, Iran University of Science and Technology, Iran. 2 Dept of Env. Eng., University of Western
More informationNUMERICAL SIMULATION OF OPEN CHANNEL FLOW BETWEEN BRIDGE PIERS
TASK QUARTERLY 15 No 3 4, 271 282 NUMERICAL SIMULATION OF OPEN CHANNEL FLOW BETWEEN BRIDGE PIERS MICHAŁ SZYDŁOWSKI Faculty of Civil and Environmental Engineering, Gdansk University of Technology, Narutowicza
More informationSolving Pure Torsion Problem and Modelling Radionuclide Migration Using Radial Basis Functions
International Workshop on MeshFree Methods 3 1 Solving Pure Torsion Problem and Modelling Radionuclide Migration Using Radial Basis Functions Leopold Vrankar (1), Goran Turk () and Franc Runovc (3) Abstract:
More informationNew computation method for flood flows and bed variations in a low-lying river with complex river systems
River Flow 2014 Schleiss et al. (Eds) 2014 Taylor & Francis Group, London, ISBN 978-1-138-02674-2 New computation method for flood flows and bed variations in a low-lying river with complex river systems
More informationESTIMATION OF MORPHOLOGICAL IMPACT OF GROYNE LENGTHENING I. RÁTKY, ÉVA RÁTKY
ESTIMATION OF MORPHOLOGICAL IMPACT OF GROYNE LENGTHENING I. RÁTKY, ÉVA RÁTKY Abstract. Hydraulic-morphological calculations in open channel flows still cause problems for modellers, partially because of
More informationVector Calculus. A primer
Vector Calculus A primer Functions of Several Variables A single function of several variables: f: R $ R, f x (, x ),, x $ = y. Partial derivative vector, or gradient, is a vector: f = y,, y x ( x $ Multi-Valued
More informationThe shallow water approximation applied to the aluminium electrolysis process
The shallow water approximation applied to the aluminium electrolysis process L Leboucher, V. Bojarevics, K Pericleous Centre for Numerical Modelling and Process Analysis, University of Greenwich, Wellington
More informationRules of Differentiation
Rules of Differentiation The process of finding the derivative of a function is called Differentiation. 1 In the previous chapter, the required derivative of a function is worked out by taking the limit
More informationEffects of possible land reclamation projects on siltation in the Rotterdam harbour area. A model study.
Effects of possible land reclamation projects on siltation in the Rotterdam harbour area. A model study. J.M. de Kok
More informationNumerical Experiment on the Fortnight Variation of the Residual Current in the Ariake Sea
Coastal Environmental and Ecosystem Issues of the East China Sea, Eds., A. Ishimatsu and H.-J. Lie, pp. 41 48. by TERRAPUB and Nagasaki University, 2010. Numerical Experiment on the Fortnight Variation
More informationB-1. Attachment B-1. Evaluation of AdH Model Simplifications in Conowingo Reservoir Sediment Transport Modeling
Attachment B-1 Evaluation of AdH Model Simplifications in Conowingo Reservoir Sediment Transport Modeling 1 October 2012 Lower Susquehanna River Watershed Assessment Evaluation of AdH Model Simplifications
More informationWIND EFFECTS ON CHEMICAL SPILL IN ST ANDREW BAY SYSTEM
WIND EFFECTS ON CHEMICAL SPILL IN ST ANDREW BAY SYSTEM PETER C. CHU, PATRICE PAULY Naval Postgraduate School, Monterey, CA93943 STEVEN D. HAEGER Naval Oceanographic Office, Stennis Space Center MATHEW
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More informationOptimal control problems with PDE constraints
Optimal control problems with PDE constraints Maya Neytcheva CIM, October 2017 General framework Unconstrained optimization problems min f (q) q x R n (real vector) and f : R n R is a smooth function.
More informationSOLVING DIFFERENTIAL EQUATIONS
SOLVING DIFFERENTIAL EQUATIONS ARICK SHAO 1. Solving Ordinary Differential Equations Consider first basic algebraic equations with an unknown, e.g., 3(x + 2) = 18, 3x 2 = cos x. Such equations model some
More informationRational derivation of the Boussinesq approximation
Rational derivation of the Boussinesq approximation Kiyoshi Maruyama Department of Earth and Ocean Sciences, National Defense Academy, Yokosuka, Kanagawa 239-8686, Japan February 22, 2019 Abstract This
More information3.1 Derivative Formulas for Powers and Polynomials
3.1 Derivative Formulas for Powers and Polynomials First, recall that a derivative is a function. We worked very hard in 2.2 to interpret the derivative of a function visually. We made the link, in Ex.
More informationApplication of the Tidal Random Walk Model to Calculate Water Dispersion Coefficients in the Ems Estuary
Estuarine, Coastal and Shelf Science (997) 45, 23 33 Application of the Tidal Random Walk Model to Calculate Water Dispersion Coefficients in the Ems Estuary H. E. de Swart a, V. N. De Jonge b and M. Vosbeek
More informationNumerical Modeling of Inclined Negatively Buoyant Jets
Numerical Modeling of Inclined Negatively Buoyant Jets Presentation by: Hossein Kheirkhah Graduate Student in Civil Engineering Dep. of Civil Engineering University of Ottawa CANADA ICDEMOS April2014 Outline
More informationThe most common methods to identify velocity of flow are pathlines, streaklines and streamlines.
4 FLUID FLOW 4.1 Introduction Many civil engineering problems in fluid mechanics are concerned with fluids in motion. The distribution of potable water, the collection of domestic sewage and storm water,
More informationTwo-dimensional urban flood simulation: Fukuoka flood disaster in 1999
Flood Recovery, Innovation and Response I 59 Two-dimensional urban flood simulation: Fukuoka flood disaster in 999 H. Hashimoto & K. Park Department of Civil Engineering, Kyushu University, Japan Abstract
More informationIntroduction to Marine Hydrodynamics
1896 1920 1987 2006 Introduction to Marine Hydrodynamics (NA235) Department of Naval Architecture and Ocean Engineering School of Naval Architecture, Ocean & Civil Engineering Shanghai Jiao Tong University
More informationThe Fundamental Theorem of Calculus: Suppose f continuous on [a, b]. 1.) If G(x) = x. f(t)dt = F (b) F (a) where F is any antiderivative
1 Calulus pre-requisites you must know. Derivative = slope of tangent line = rate. Integral = area between curve and x-axis (where area can be negative). The Fundamental Theorem of Calculus: Suppose f
More informationUNIFORM FLOW CRITICAL FLOW GRADUALLY VARIED FLOW
UNIFORM FLOW CRITICAL FLOW GRADUALLY VARIED FLOW Derivation of uniform flow equation Dimensional analysis Computation of normal depth UNIFORM FLOW 1. Uniform flow is the flow condition obtained from a
More informationNumerical Methods for Problems with Moving Fronts Orthogonal Collocation on Finite Elements
Electronic Text Provided with the Book Numerical Methods for Problems with Moving Fronts by Bruce A. Finlayson Ravenna Park Publishing, Inc., 635 22nd Ave. N. E., Seattle, WA 985-699 26-524-3375; ravenna@halcyon.com;www.halcyon.com/ravenna
More informationconservation of linear momentum 1+8Fr = 1+ Sufficiently short that energy loss due to channel friction is negligible h L = 0 Bernoulli s equation.
174 Review Flow through a contraction Critical and choked flows The hydraulic jump conservation of linear momentum y y 1 = 1+ 1+8Fr 1 8.1 Rapidly Varied Flows Weirs 8.1.1 Broad-Crested Weir Consider the
More informationChapter 5. Shallow Water Equations. 5.1 Derivation of shallow water equations
Chapter 5 Shallow Water Equations So far we have concentrated on the dynamics of small-scale disturbances in the atmosphere and ocean with relatively simple background flows. In these analyses we have
More informationCHAPTER TWO HUNDRED FOUR
CHAPTER TWO HUNDRED FOUR Lateral Distributions of Water, Salt and Sediment Transport in a Partly Mixed Estuary R.J. Uncles, R.C.A. Elliott and S.A. Weston The transverse structure of the residual transport
More informationEXAMPLES (SEDIMENT TRANSPORT) AUTUMN 2018
EXAMPLES (SEDIMENT TRANSPORT) AUTUMN 2018 Q1. Using Cheng s formula estimate the settling velocity of a sand particle of diameter 1 mm in: (a) air; (b) water. Q2. Find the critical Shields parameter diameter
More informationHorizontal buoyancy-driven flow along a differentially cooled underlying surface
Horizontal buoyancy-driven flow along a differentially cooled underlying surface By Alan Shapiro and Evgeni Fedorovich School of Meteorology, University of Oklahoma, Norman, OK, USA 6th Baltic Heat Transfer
More informationB O S Z. - Boussinesq Ocean & Surf Zone model - International Research Institute of Disaster Science (IRIDeS), Tohoku University, JAPAN
B O S Z - Boussinesq Ocean & Surf Zone model - Volker Roeber 1 Troy W. Heitmann 2, Kwok Fai Cheung 2, Gabriel C. David 3, Jeremy D. Bricker 1 1 International Research Institute of Disaster Science (IRIDeS),
More informationVARIED FLOW IN OPEN CHANNELS
Chapter 15 Open Channels vs. Closed Conduits VARIED FLOW IN OPEN CHANNELS Fluid Mechanics, Spring Term 2011 In a closed conduit there can be a pressure gradient that drives the flow. An open channel has
More informationNumerical Hydraulics
ETHZ, Fall 017 Numerical Hydraulics Assignment 4 Numerical solution of 1D solute transport using Matlab http://www.bafg.de/ http://warholian.com Numerical Hydraulics Assignment 4 ETH 017 1 Introduction
More informationA NOTE ON THE CONTRIBUTION OF DISPERSIVE FLUXES TO MOMENTUM TRANSFER WITHIN CANOPIES. Research Note
A NOTE ON THE CONTRIBUTION OF DISPERSIVE FLUXES TO MOMENTUM TRANSFER WITHIN CANOPIES Research Note D. POGGI Dipartimento di Idraulica, Trasporti ed Infrastrutture Civili, Politecnico di Torino, Torino,
More informationLecture two. January 17, 2019
Lecture two January 17, 2019 We will learn how to solve rst-order linear equations in this lecture. Example 1. 1) Find all solutions satisfy the equation u x (x, y) = 0. 2) Find the solution if we know
More informationForced Convectional Heat Transfer in Solid Oxide Fuel Cells: An Analytical Treatment.
Ionics 9 (2003) 83 Forced Convectional Heat Transfer in Solid Oxide Fuel Cells: An Analytical Treatment. F.A. Coutelieris 1, A.K. Demin 2, S.L. Douvartzides 1 and P.E. Tsiakaras 1 1 University of Thessalia,
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationBeyond Ocean Modelling: Multi-Scale/Physics Numerical Simulation of the Hydrosphere I. Numerics and Preliminary Applications
1 Mathematics in Waterland COSSE Workshop Delft University of Technology, Delft, The Netherlands February 7-10, 2011 Beyond Ocean Modelling: Multi-Scale/Physics Numerical Simulation of the Hydrosphere
More information2. Governing Equations
1. Introduction Submarine pipeline, unlike any other hydraulic structures that are vertically erected, are laid horizontally on the bed of oceans and rivers. Hence, the design of submarine pipelines associated
More information39.1 Gradually Varied Unsteady Flow
39.1 Gradually Varied Unsteady Flow Gradually varied unsteady low occurs when the low variables such as the low depth and velocity do not change rapidly in time and space. Such lows are very common in
More information