THREE-DIMENSIONAL FINITE DIFFERENCE MODEL FOR TRANSPORT OF CONSERVATIVE POLLUTANTS
|
|
- Jeffry McDaniel
- 5 years ago
- Views:
Transcription
1 Pergamon Ocean Engng, Vol. 25, No. 6, pp , Elsevier Science Ltd. All rights reserved Printed in Great Britain /98 $ PII: S (97) THREE-DIMENSIONAL FINITE DIFFERENCE MODEL FOR TRANSPORT OF CONSERVATIVE POLLUTANTS S. Sankaranarayanan, N. J. Shankar* and H. F. Cheong Department of Civil Engineering, National University of Singapore, 10 Kent Ridge Crescent, Singapore Abstract A three-dimensional finite difference transport model appropriate for the coastal environment is developed for the solution of the three-dimensional convection diffusion equation. A higher order upwind scheme is used for the convective terms of the convection diffusion equation, to minimise the numerical diffusion. The validity of the numerical model is verified through five test problems, whose exact solutions are known Elsevier Science Ltd. NOMENCLATURE C(x, y, z, t) pollutant concentration at location (x, y, z) at any time t C 1, C 2, C 3, constants C i,j,k concentration at grid location i, j, k G(x, y, z, t) point source or sink K scaling constant K x, K y, K z diffusion coefficients in the x, y and z directions L x, L y, L z length of the basin in x, y and z directions r(t), s(t), f(t) expression [see Equation (26)] T p period of rotation U, V, W flow velocities in the x, y and z directions x 0, y 0, z 0 coordinates of the centre of Gaussian pulse dimensionless parameter t time step x, y, z grid spacing in x, y and z directions free surface elevation 1. INTRODUCTION Mathematical modelling of the transport of salinity, pollutants and suspended matter in shallow waters involves the numerical solution of a convection diffusion equation. Many popular finite difference methods, such as the upwind scheme of Spalding (1972) and the flux-corrected scheme (Boris and Book, 1973) are available for the solution of the depthintegrated form of the convection diffusion equation. Another widely used approach is the split-operator approach (Sobey, 1983; Li and Chen, 1989), in which the convection and diffusion terms are solved by two different numerical methods. Noye and Tan (1988) used a weighted discretisation with the modified equivalent partial differential equation *Author to whom correspondence should be addressed. 425
2 426 S. Sankaranarayanan et al. approach for solving one-dimensional convection diffusion equations. Later, Noye and Tan (1989) extended this technique to solve two-dimensional convection diffusion equations. But the above mentioned techniques have difficulty in solving three-dimensional problems, because of extensive matrix inversions at each time step. Numerical studies show that the use of central differencing for the convective terms of the convection diffusion equation results in negative species concentration. Lam (1975) points out that the central difference approximation will overestimate the advective flux so much that it often causes a negative concentration to appear in the neighbouring cell. To circumvent such a shortcoming of central differencing, the upwind or donor cell method introduced by Gentry et al. (1966) is generally used. To overcome the shortcomings of numerical dispersion, Leonard (1979) introduced an upstream interpolation method, namely QUICK (Quadratic Upstream Interpolation Convective Kinematics) for one-dimensional unsteady flow. Later, Leonard (1988) gave an improved version of the QUICK scheme, eliminating the wiggles completely by introducing exponential integration into regions with sharp fronts. Chen and Falconer (1994) showed that the QUICK scheme is only second-order accurate in space and presented different forms of the third-order convection, second-order diffusion for the solution of the convection diffusion equation. A relatively recent phenomenon is the development of three-dimensional transport models. Lardner and Song (1991) used an algorithm splitting the horizontal convection and diffusion by an implicit finite element method, treating the horizontal convection diffusion explicitly and vertical convection diffusion by an implicit finite element method. It is to be noted that Lardner and Song (1991) used a first-order upwind scheme for the convection terms of the convection diffusion equations. Sommeijer and Kok (1995) made a detailed study on the use of various time-integration techniques for the numerical solution of threedimensional convection diffusion equations using finite differences. The numerical model was validated by comparing the results obtained with analytical solutions for the case of transport of a Gaussian pulse in unsteady and non-uniform flow. In the present study, a third-order upwind difference scheme as given in Kowalik and Murty (1993) has been used for the convection terms of the convection diffusion equation. Earlier, the authors Shankar et al. (1996) used a third-order upwind scheme for the convective terms of the shallow water momentum equations. 2. MATHEMATICAL MODELLING OF CONVECTION DIFFUSION PROCESSES 2.1. Governing equations The mathematical model describing the transport processes in three dimensions is given by the convection diffusion equation t + (CU) + (CV) + (CW) 2 C K x y z x x K 2 C 2 y y K 2 C 2 z z 2 = G(x, y, z, t) (1) subject to the initial conditions, C(x, y, z, 0), and the boundary conditions, Equations (2) (4). (x, y, t) at z = 0 on the top surface (2) z
3 3-D finite difference transport model 427 (x, y, t) at z = h on the bottom surface (3) z = 0 along the lateral boundaries (4) n where C is the unknown pollutant concentration at location (x, y, z), U, V and W are the flow velocities in the x, y and z directions, respectively, K x, K y and K z are the diffusion coefficients in the x, y and z directions, respectively, G represents any source or sink. x, y and z are the Cartesian coordinates with the xy plane horizontal and occupying the undisturbed position of the water surface, and the z-axis pointing vertically upwards, and n is the outward drawn normal in the xy plane. The position of the free surface is denoted by z = (x, y, t) as shown in Fig. 1 and that of the bottom by z = h(x, y). The values of, U, V and W are to be obtained from a hydrodynamic model such as that of Shankar et al. (1996). Equation (1) reduces to t + U x + V y + W z K 2 C x x K 2 C 2 y y K 2 C 2 z z = G(x, y, z, t) 2 (5) since U x + V y + W = 0 z 2.2. Finite difference formulation of the convection term using first-order upwinding The values of C, U, V and W are defined at all the grid points i, j and k, respectively. A typical convection term / x is evaluated using first-order differencing as x = C i,j,k C i 1,j,k x C i +1,j,k C i,j,k x + O( x) if U(i, j, k) 0 + O( x) if U(i, j, k) 0 (6) Fig. 1. Definition sketch.
4 428 S. Sankaranarayanan et al. The first-order upwind scheme as given by Equation (6) is very popular in the numerical solution of the convection diffusion equation (Lardner and Song, 1991) and shallow water hydrodynamic equations (Shankar et al., 1995; Horie, 1980). However, the use of a firstorder upwind scheme is strongly discouraged by numerical modellers, since it introduces numerical diffusion. In other words, when a first-order upwind scheme is used for the convection terms, it amounts to introducing one additional term in the governing equations (O Brien, 1986) Finite difference formulation of the convection term using third-order upwinding In the present study, a third-order upwind differencing has been used as given in Kowalik and Murty (1993). A typical convective term / x is discretised using a third-order differencing as x = (C i 2,j,k 6C i 1,j,k +3C i,j,k +2C i +1,j,k )/6( x) +o( x) 3 (7) for U i,j,k 0 and i = 3,nx 2. x = ( 2C i 1,j,k 3C i,j,k +6C i +1,j,k C i +2,j,k )/6( x) +o( x) 3 (8) for U i,j,k 0 and i = 3,nx 2. However, for the boundary points, a four-point upstream formula can be written such that either points to the left or to the right are considered in the finite difference approximation. Near the left boundary, i.e. for i = 2 x = ( 11C i,j,k +18C i +1,j,k 9C i +2,j,k +2C i +3,j,k )/6( x) +o( x) 3 (9) Near the right boundary, i.e. for i = nx 1 x = ( 2C i 3,j,k +9C i 2,j,k 18C i 1,j,k +11C i,j,k )/6( x) +o( x) 3 (10) where nx is the number of grid points in the x direction. The second derivatives occurring in diffusion terms are evaluated using a central difference. A typical second derivative 2 C/ x 2 is evaluated using 2 C x = C i +1,j,k 2C i,j,k + C i 1,j,k (11) 2 ( x) Stability criterion Hindmarsh et al. (1984) established the stability criterion for the multi-dimensional advection diffusion equation for an explicit scheme as t 2 K x ( x) 2 + K y ( y) K z ( z) 2 + U ( x) + V ( y) + W ( z) (12)
5 3-D finite difference transport model VALIDATION OF THE NUMERICAL MODEL It is necessary to compare the model results with analytical solutions, to validate the developed numerical model. The simplified idealised test cases, for which analytical solutions are available, present a valuable means to debug the developed computer program. Two such numerical examples used by Noye and Tan (1988, 1989) involving transport of a one-dimensional and two-dimensional Gaussian pulse, respectively, and two test cases used by Lardner and Song (1991), involving two-dimensional and three-dimensional transport of a point source pollutant are used to validate the model developed in the present study. The test cases involving the transport of point source pollutant concentration by convection and diffusion are taken, since they represent an important practical case of a sudden influx of a pollutant, such as due to oil spill. The four test cases used in the present study have analytical solutions, which makes it possible to compare the results obtained in the numerical solution. MKS units are used throughout the paper Test case involving one-dimensional convection and diffusion For one-dimensional convection diffusion, Equation (1) reduces to t + U x K x 2 C = G(x, y, z, t) (13) 2 x The analytical solution to the one-dimensional advection and diffusion of a Gaussian pulse of unit height, centred at x = 1 in a region bounded by 0 x 9 as given by Noye and Tan (1988) is C(x, t) = 1 4t +1 exp (x x 0 Ut) 2 K x (4t +1) (14) where U is the velocity in the x direction, x 0 is the centre of the initial Gaussian pulse, K x is the diffusion coefficient in the x direction and t is the time coordinate. The initial condition is given by C(x, t) = 1 4t +1 exp (x 1)2 K x (15) and the boundary condition at the two ends at any time t is obtained by substituting x = 0 and x = 9, respectively, in Equation (14). The values of the various parameters used are K x = m 2 /s, U = 0.8 m/s, and the Gaussian pulse of unit height centred at x = 1. The distribution of the Gaussian pulse at t = 5 s is computed using the analytical solution and compared with the concentration distribution obtained using the numerical solution as shown in Fig. 2. The space step and time step are taken to be m and s, respectively. The Courant number (U t/ x) works out to be 0.2 and the cell Reynolds number (U x/k x ) works out as Test case involving two-dimensional convection and diffusion of a Gaussian pulse For two-dimensional convection diffusion, Equation (1) reduces to t + U x + V y K 2 C x x K 2 y 2 C = G(x, y, z, t) (16) 2 y
6 430 S. Sankaranarayanan et al. Fig. 2. Comparison of analytical and numerical solutions for transport of one-dimensional Gaussian pulse for K x = m 2 /s, U = 0.8 m/s. An analytical solution for solving the two-dimensional convection diffusion equation in a rectangular region bounded by 0 x 2 and 0 y 2 for the transport of a Gaussian pulse of unit height (Fig. 3) centred at (0.5, 0.5) is given by C(x, y, t) = 1 4t +1 exp (x 0.5 Ut)2 (y 0.5 Vt)2 K x (4t +1) K y (4t +1) (17) Fig. 3. Three-dimensional perspective view of the initial Gaussian pulse.
7 3-D finite difference transport model 431 with the initial condition C(x, y,0)= exp (x 0.5)2 K x (y 0.5)2 K y (18) and the values of the concentrations along the boundaries at any time t are obtained by substituting the values of their coordinates in Equation (17). The parameters used in the simulation are U = V = 0.8 m/s, K x = K y = 0.01 m 2 /s. The space and time steps are taken to be m and s and hence the Courant number and the cell Reynolds number are calculated to be 0.4 and 2.0, respectively. According to the analytical solution (Fig. 4), the pulse has moved to x = 1.5 and y = 1.5 with a pulse height of at time t = 5 s, while the numerical solution predicts the pulse height as The concentration distributions computed using numerical and analytical solutions are given in Figs 4a and 4b, respectively. Fig. 4. Concentration distributions by analytical and numerical solutions after 5 s, x = m, t = s, K x = K y = 0.01 m 2 /s, U = V = 0.8 m/s: (a) analytical solution; (b) numerical solution.
8 432 S. Sankaranarayanan et al Test case involving two-dimensional convection and diffusion of a point source For an infinite region, the exact solution of the two-dimensional convection diffusion Equation (16), the concentration C(x, y, t) at a given time t for a point source pollutant placed at the centre of the domain (x 0, y 0 ) is given by C(x, y, t) = K 4 K x t 4 K y t exp (x x 0 Ut) 2 (y y 0 Vt) 2 4K x t 4K y t (19) where (x 0, y 0 ) is the location of the point source and K is a scaling factor, taken to be for this test case. The transport of an instantaneous point source pollutant, initially placed at the centre of a domain of size 10,000 10,000 m 2, by convection and diffusion is simulated using the model. Uniform horizontal velocities of U = V = 0.5 m/s, grid spacing x = y = 5000 m, time step t = 50 s and K x = K y = 10,000 m 2 /s are used. The solution for the bounded region is taken large enough so that the concentrations remain essentially zero Fig. 5. Concentration distributions for two-dimensional convection and diffusion by analytical and numerical solutions after 36,000 s, x = y = 5000 m, t = 50 s, K x = K y = 10,000 m 2 /s, U = V = 0.5 m/s: (a) analytical solution; (b) numerical solution.
9 3-D finite difference transport model 433 Fig. 6. Concentration distributions for three-dimensional convection and diffusion by analytical and numerical solutions, at z = 13 m from the surface, after 23,000 s, x = y = 2000 m, z = 6.5 m, t = 50 s, K x = K y = K y = 2000 m 2 /s, K z = 0.01 m 2 /s, U = V = 0.2 m/s, W = 0.0 m/s, t 0 = 5000 s: (a) analytical solution; (b) numerical solution. at the boundaries for the time interval of the computation (Lardner and Song, 1991). The concentration distribution obtained after t = 36,000 s using the model developed in the present study compares well with that given by the analytical solution, Equation (19), as shown in Fig. 5aFig. 5b, respectively. The maximum relative error between the analytical and numerical solutions is about 10% Test case involving three-dimensional convection and diffusion in uniform flow For an infinite region, the exact solution C(x, y, z, t) at a given time t for a point source placed at the centre of the domain is given by
10 434 S. Sankaranarayanan et al. Fig. 7. Concentration distributions for three-dimensional convection and diffusion by analytical and numerical solutions, at z = 32.5 m from the surface, after 23,000 s, x = y = 2000 m, z = 6.5 m, t = 50 s, K x = K y = 2000 m 2 /s, K z = 0.01 m 2 /s, U = V = 0.2 m/s, W = 0.0 m/s, t 0 = 5000 s: (a) analytical solution; (b) numerical solution. where C(x, y, z, t) = C x (x, t, U)C y (y, t, V)C z (z, t, W) (20) 1 C x (x, U, t) = 4 K x t exp (x x 0 Ut) 2 4K x t (21) C y (y, V, t) = C x (z, W, t) = 1 4 K y t exp (y y 0 Vt) 2 4K y t (22) 1 4 K z t exp (z z 0 Wt) 2 4K z t (23)
11 3-D finite difference transport model 435 Fig. 8. Concentration distributions for three-dimensional convection and diffusion by analytical and numerical solutions, at z = 52 m from the surface, after 23,000 s, x = y = 2000 m, z = 6.5 m, t = 50 s, K x = K y = 2000 m 2 /s, K z = 0.01 m 2 /s, U = V = 0.2 m/s, W = 0.0 m/s, t 0 = 5000 s. where (x 0, y 0 ) is the location of the point source and K is a scaling factor, taken to be for this test case. For applying the boundary conditions, Equations (2) and (3), the derivative of C(x, y, z, t) is evaluated as [C(x, y, z, t)] z = C x (x, t, U)C y (y, t, V)C z (z, t, W) [z 0 + w(t + t 0 ) z] 2[K z (t + t 0 )] Substituting the values of z = 0 and z = h in Equation (18), the boundary condition Equations (2) and (3) at the top and bottom surfaces are satisfied in the numerical solution. The initial condition is given by C(x, y, z, t 0 ), where if t 0 0, then the problem has a smooth initial condition. The transport of an instantaneous point source pollutant, initially placed at the centre of a domain of size 40,000 40, m 3, by convection and diffusion is simulated (24)
12 436 S. Sankaranarayanan et al. using the model. Uniform horizontal velocities of U = V = 0.2 m/s and W = 0 m/s, grid spacing x = y = 2000 m, time step t = 50 s and K x = K y = 2000 m 2 /s, K z = 0.01 m 2 /s are used. For a smooth initial condition, the value of t 0 is taken to be 5000 s. The solution for the bounded region is taken large enough so that the solution remains essentially zero at the boundaries for the time interval of the computation (Lardner and Song, 1991). The concentration distribution obtained after t = 18,000 s using the model developed in the present study compares well with that given by the analytical solution Equation (20) at various levels, showing good comparison as in Figs 6 8, respectively. The maximum relative error between the analytical and numerical solutions is about 13% Test case involving three-dimensional convection and diffusion in unsteady non-uniform flow The analytical solution for the transport of conservative pollutant in unsteady, nonuniform flow as given by Sommeijer and Kok (1995) is C(x, y, z, t) = exp z L z f(t) x L x r(t) 2 y L y s(t) 2 (25) To stress the local nature of the concentration, the value of is taken to be 30. r(t) = cos 2 t T p ; s(t) = sin 2 t T p and f(t) = 4t (T b + t) (26) The flow fields are given by U(x, y, z, t) = C 1 sin x L x + y L y sin z L z d(t) (27) Fig. 9. Concentration distribution C(x, y, z, t) for t = 0 at the surface (z = 0).
13 3-D finite difference transport model 437 Fig. 10. Concentration distribution C(x, y, z, t) for t = 0 at z = 50 m from the top. Fig. 11. Concentration distribution C(x, y, z, t) for t = 0 at z = 100 m from the top: (a) analytical solution; (b) numerical solution.
14 438 S. Sankaranarayanan et al. Fig. 12. Concentration distributions for three-dimensional convection and diffusion by analytical and numerical solutions, at the surface (z = 0) after 10,800 s, x = y = 2000 m, z = 10 m, t = 50 s, K x = K y = K z = 0.5 m 2 /s: (a) analytical solution; (b) numerical solution. where V(x, y, z, t) = C 1 cos x L x + y L y sin z L z d(t) (28) W(x, y, z, t) = C 1 cos x L x + y L y + C 2 sin x L x + y L y L z cos z L z d(t) d(t) = cos 2 t T p C 1 = 3, C 2 = 4, and = 0.05 The flow velocities given by Equations (27) (29) satisfy the continuity relation for the mass balance given by (29)
15 3-D finite difference transport model 439 Fig. 13. Concentration distributions for three-dimensional convection and diffusion by analytical and numerical solutions, at z = 50 m from the surface, after 10,800 s, x = y = 2000 m, z = 10 m, t = 50 s, K x = K y = K z = 0.5 m 2 /s: (a) analytical solution; (b) numerical solution. U x + V y + W z = 0 (30) The values of the various parameters are given by L x = L y = 20,000 m, L z = 100 m, K x = K y = K z = 0.5 m 2 /s, T p = 43,200 s, T b = 32,400 s. It is to be noted that the source term, G(x, y, z, t), is evaluated at each time step from Equation (1), as given in Appendix A. The initial concentration distribution at t = 0 is given by substituting t = 0 in Equation (25). The typical concentration distributions at the surface, 50 m from the top and at the
16 440 S. Sankaranarayanan et al. Fig. 14. Concentration distributions for three-dimensional convection and diffusion by analytical and numerical solutions, at the bottom (z = 100 m) after 10,800 s, x = y = 2000 m, z = 10 m, t = 50 s, K x = K y = K z = 0.5 m 2 /s. bottom of the basin, given as initial conditions are shown in Figs The space and time steps are taken to be 2000 m and 50 s. The concentration distribution at time t = 10,800 s obtained at depths of 0, 50 and 100 m from the top obtained in the present study compares well with the analytical solution as shown in Figs The maximum percentage error between analytical and numerical solutions is found to be 11%.
17 3-D finite difference transport model CONCLUSIONS A three-dimensional finite difference transport model appropriate for the coastal environment is developed for the numerical solution of the three-dimensional convection diffusion equation. A higher order (third-order) upwind scheme is used for the convective terms of the transport equation, to minimise the numerical diffusion. The validity of the numerical model is verified through four test problems, for the transport of point source pollutants in uniform flow, by comparing the results obtained in the present study with those of analytical solutions. The validity of the numerical solution is also tested against analytical solutions for the case of transport of a Gaussian pulse in unsteady and nonuniform flow. However, for problems involving transport of pollutants in coastal waters, the velocity field is to be obtained a priori from the solution of a compatible hydrodynamic model developed by the authors (Shankar et al., 1996). Work is currently in progress to couple the transport model developed in the present study with the multi-level three-dimensional hydrodynamic model. REFERENCES Boris, J. B. and Book, D. L. (1973) Flux corrected for transport algorithm that works. Journal of Computational Physics 11, Chen, Y. and Falconer, R. A. (1994) Modified forms of the third-order convection, second-order diffusion equation. Advances in Water Resources 17, Gentry, R. A., Martin, R. E. and Daly, B. J. (1966) An Eulerian differencing method for unsteady compressible flow problems. Journal of Computational Physics 8, Hindmarsh, A. C., Gresho, P. and Griffiths, D. F. (1984) The stability of explicit Euler integration for certain finite difference approximations of the multidimensional advection diffusion equation. International Journal for Numerical Methods in Fluids 4, Horie, T. (1980) Hydraulic investigation on sea flow and substance dispersion in estuarine and coastal region. Technical Note 360, Port and Harbour Research Institute, Japan. Kowalik, Z. and Murty, T. S. (1993) Numerical Modeling of Ocean Dynamics. World Scientific, Singapore. Lam, D. C. L. (1975) Computer modeling of pollutant transport in Lake Erie. Water Pollution 25, Lardner, R. W. and Song, Y. (1991) An algorithm for three-dimensional convection and diffusion with very different horizontal and vertical scales. International Journal for Numerical Methods in Engineering 32, Leonard, B. P. (1979) A stable and accurate convective modeling procedure based on upstream formulation. Computer Methods in Applied Mechanics and Engineering 19, Leonard, B. P. (1988) Simple high accuracy resolution program for convective modeling of discontinuities. International Journal for Numerical Methods in Fluids 8, Li, Y. S. and Chen, C. P. (1989) An efficient split operator scheme for 2D advection diffusion equation using finite elements and characteristics. Applied Mathematical Modeling 13, Noye, B. J. and Tan, H. H. (1988) A third-order semi-implicit finite difference method for solving the onedimensional convection diffusion equation. International Journal for Numerical Methods in Engineering 26, Noye, B. J. and Tan, H. H. (1989) Finite difference methods for the two-dimensional advection diffusion equation. International Journal for Numerical Methods in Fluids 9, O Brien, J. J. (1986) Advanced Physical Oceanographic Modeling, p Riedel, Boston. Shankar, N. J., Cheong, H. F. and Sankaranarayanan, S. (1995) Multilevel finite difference model for threedimensional hydrodynamic circulation. Internal Report, Department of Civil Engineering, National University of Singapore, p. 55. Shankar, N. J., Cheong, H. F. and Sankaranarayanan, S. (1996) Multilevel finite difference model for threedimensional hydrodynamic circulation. International Journal of Ocean Engineering, 24(9), Sobey, R. J. (1983) Fractional step algorithm for estuarine mass transport. International Journal for Numerical Methods in Fluids 3, Sommeijer, B. P. and Kok, J. (1995) Implementation and performance of the time integration of a 3D numerical transport model. International Journal for Numerical Methods in Fluids 21, Spalding, D. B. (1972) A novel finite difference formulation for differential expressions involving both first and second derivatives. International Journal for Numerical Methods in Engineering 4,
18 442 S. Sankaranarayanan et al. APPENDIX A From Equation (1) the value of G(x, y, z, t) is evaluated to be where G(x, y, z, t) = t 2 C K x x K 2 y + U(x, y, z, t) x 2 C y K 2 C 2 z z 2 + V(x, y, z, t) + W(x, y, z, t) y z (A1) C(x, y, z, t) = exp z L z f(t) x L x r(t) 2 y L y s(t) 2 (A2) r(t) = cos 2 t s(t) = 1 T p ; sin 2 t 4t T p and f(t) = (A3) (T b + t) t = C(x, y, z, t) 4t + T b + t x L x r(t) sin 2 t T p (A4) + T p y L y s(t) cos 2 t T p 4t (T b + t) + 2 T p x L x r(t) (A5) ( 2 ) = C(x, y, z, t) x L x ( 2 ) = C(x, y, z, t) y L y x L y s(t) (A6) z = C(x, y, z, t) 1 (A7) L z 2 C x = C(x, y, z, t) ( 2 ) (L x ) 2 (L x ) x 2 L x r(t) 2 (A8) 2 C y = C(x, y, z, t) ( 2 ) (L y ) 2 (L y ) y 2 L y s(t) 2 (A9) 2 C = C(x, y, z, t) (A10) 2 z (L z ) 2
Numerical Solution of One-dimensional Advection-diffusion Equation Using Simultaneously Temporal and Spatial Weighted Parameters
Australian Journal of Basic and Applied Sciences, 5(6): 536-543, 0 ISSN 99-878 Numerical Solution of One-dimensional Advection-diffusion Equation Using Simultaneously Temporal and Spatial Weighted Parameters
More informationChapter 3. Finite Difference Methods for Hyperbolic Equations Introduction Linear convection 1-D wave equation
Chapter 3. Finite Difference Methods for Hyperbolic Equations 3.1. Introduction Most hyperbolic problems involve the transport of fluid properties. In the equations of motion, the term describing the transport
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 informationSolution Methods. Steady convection-diffusion equation. Lecture 05
Solution Methods Steady convection-diffusion equation Lecture 05 1 Navier-Stokes equation Suggested reading: Gauss divergence theorem Integral form The key step of the finite volume method is to integrate
More information3.4. Monotonicity of Advection Schemes
3.4. Monotonicity of Advection Schemes 3.4.1. Concept of Monotonicity When numerical schemes are used to advect a monotonic function, e.g., a monotonically decreasing function of x, the numerical solutions
More informationBACHELOR OF TECHNOLOGY IN MECHANICAL ENGINEERING (COMPUTER INTEGRATED MANUFACTURING)
No. of Printed Pages : 6 BME-028 BACHELOR OF TECHNOLOGY IN MECHANICAL ENGINEERING (COMPUTER INTEGRATED MANUFACTURING) Term-End Examination December, 2011 00792 BME-028 : FLUID MECHANICS Time : 3 hours
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 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 informationA semi-lagrangian scheme for advection-diffusion equation
EPiC Series in Engineering Volume 3, 2018, Pages 162 172 HIC 2018. 13th International Conference on Hydroinformatics Engineering A semi-lagrangian scheme for advection-diffusion equation Ersin Bahar 1,
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 informationA Taylor series method for numerical fluid mechanics
Eighth Australasian Fluid Mechanics Conference University of Newcastle, N.S.W. 28 November 2 December 1983, Retyped in 2011. A Taylor series method for numerical fluid mechanics J.D. Fenton School Of Mathematics
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 information2. Conservation of Mass
2 Conservation of Mass The equation of mass conservation expresses a budget for the addition and removal of mass from a defined region of fluid Consider a fixed, non-deforming volume of fluid, V, called
More informationChapter 2 Finite-Difference Discretization of the Advection-Diffusion Equation
Chapter Finite-Difference Discretization of the Advection-Diffusion Equation. Introduction Finite-difference methods are numerical methods that find solutions to differential equations using approximate
More informationChapter 17. Finite Volume Method The partial differential equation
Chapter 7 Finite Volume Method. This chapter focusses on introducing finite volume method for the solution of partial differential equations. These methods have gained wide-spread acceptance in recent
More informationSolution Methods. Steady State Diffusion Equation. Lecture 04
Solution Methods Steady State Diffusion Equation Lecture 04 1 Solution methods Focus on finite volume method. Background of finite volume method. Discretization example. General solution method. Convergence.
More informationDiscretization of Convection Diffusion type equation
Discretization of Convection Diffusion type equation 10 th Indo German Winter Academy 2011 By, Rajesh Sridhar, Indian Institute of Technology Madras Guides: Prof. Vivek V. Buwa Prof. Suman Chakraborty
More informationA general technique for assessing the numerical accuracy of solute transport models
WATER RESOURCES RESEARCH, VOL. 35, NO. 12, PAGES 3961 3966, DECEMBER 1999 A general technique for assessing the numerical accuracy of solute transport models Feng Ruan Advanced Technology, Schlumberger,
More informationContents. Parti Fundamentals. 1. Introduction. 2. The Coriolis Force. Preface Preface of the First Edition
Foreword Preface Preface of the First Edition xiii xv xvii Parti Fundamentals 1. Introduction 1.1 Objective 3 1.2 Importance of Geophysical Fluid Dynamics 4 1.3 Distinguishing Attributes of Geophysical
More informationA generalized method for advective-diffusive computations in engineering
Fluid Structure Interaction and Moving Boundary Problems 563 A generalized method for advective-diffusive computations in engineering H. Gómez, I. Colominas, F. Navarrina & M. Casteleiro Department of
More informationA semi-implicit non-hydrostatic covariant dynamical kernel using spectral representation in the horizontal and a height based vertical coordinate
A semi-implicit non-hydrostatic covariant dynamical kernel using spectral representation in the horizontal and a height based vertical coordinate Juan Simarro and Mariano Hortal AEMET Agencia Estatal de
More informationNumerical Oscillations and how to avoid them
Numerical Oscillations and how to avoid them Willem Hundsdorfer Talk for CWI Scientific Meeting, based on work with Anna Mozartova (CWI, RBS) & Marc Spijker (Leiden Univ.) For details: see thesis of A.
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 information10.34 Numerical Methods Applied to Chemical Engineering. Quiz 2
10.34 Numerical Methods Applied to Chemical Engineering Quiz 2 This quiz consists of three problems worth 35, 35, and 30 points respectively. There are 4 pages in this quiz (including this cover page).
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More information3.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 informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 9 (214 ) 599 64 1th International Conference on Mechanical Engineering, ICME 213 Validation criteria for DNS of turbulent heat
More informationControl of saline water intrusion by vane-like barrier
Control of saline water intrusion by vane-like barrier M. Ifuku Department of Civil and Environmental Engineering, Ehime University, Japan. Abstract There have been many studies on density current in estuary
More informationdynamics of f luids in porous media
dynamics of f luids in porous media Jacob Bear Department of Civil Engineering Technion Israel Institute of Technology, Haifa DOVER PUBLICATIONS, INC. New York Contents Preface xvii CHAPTER 1 Introduction
More informationA new moving boundary shallow water wave equation numerical model
ANZIAM J. 48 (CTAC2006) pp.c605 C617, 2007 C605 A new moving boundary shallow water wave equation numerical model Joe Sampson 1 Alan Easton 2 Manmohan Singh 3 (Received 10 August 2006; revised 28 November
More informationLecture 12: Transcritical flow over an obstacle
Lecture 12: Transcritical flow over an obstacle Lecturer: Roger Grimshaw. Write-up: Erinna Chen June 22, 2009 1 Introduction The flow of a fluid over an obstacle is a classical and fundamental problem
More informationChapter 8 - pg. 1 CHAPTER 8 ESTUARIES. To paraphrase Pritchard, a pioneer in studies of estuarine circulation,
Chapter 8 - pg 1 CHAPTER 8 ESTUARIES Estuaries are semi-closed basins in which a rather complex interaction between river inputs, tidal currents and wind leads to the turbulent mixing of salt from the
More informationFinite difference method for solving Advection-Diffusion Problem in 1D
Finite difference method for solving Advection-Diffusion Problem in 1D Author : Osei K. Tweneboah MATH 5370: Final Project Outline 1 Advection-Diffusion Problem Stationary Advection-Diffusion Problem in
More informationSD Numerical Simulation Technique for Hydrodynamic Flow Gas-Solids Mixing
SD Numerical Simulation Technique for Hydrodynamic Flow Gas-Solids Mixing Presented by: Irla Mantilla N. irlamn@uni.edu.pe National University of Engineering, Lima - Perú Abstract We formulate a new mathematical
More informationUNIVERSITY OF EAST ANGLIA
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must
More informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
More informationFundamentals of Transport Processes Prof. Kumaran Indian Institute of Science, Bangalore Chemical Engineering
Fundamentals of Transport Processes Prof. Kumaran Indian Institute of Science, Bangalore Chemical Engineering Module No # 05 Lecture No # 25 Mass and Energy Conservation Cartesian Co-ordinates Welcome
More informationIMPORTANCE OF NUMERICAL EFFICIENCY FOR REAL TIME CONTROL OF TRANSIENT GRAVITY FLOWS IN SEWERS
1106 September 11~16, 2005, Seoul, Korea IMPORTANCE OF NUMERICAL EFFICIENCY FOR REAL TIME CONTROL OF TRANSIENT GRAVITY FLOWS IN SEWERS ARTURO S. LEON 1, MOHAMED S. GHIDAOUI 2, ARTHUR R. SCHMIDT 3 and MARCELO
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 information3.3. Phase and Amplitude Errors of 1-D Advection Equation
3.3. Phase and Amplitude Errors of 1-D Advection Equation Reading: Duran section 2.4.2. Tannehill et al section 4.1.2. The following example F.D. solutions of a 1D advection equation show errors in both
More informationFinite volume method for CFD
Finite volume method for CFD Indo-German Winter Academy-2007 Ankit Khandelwal B-tech III year, Civil Engineering IIT Roorkee Course #2 (Numerical methods and simulation of engineering Problems) Mentor:
More informationNodalization. The student should be able to develop, with justification, a node-link diagram given a thermalhydraulic system.
Nodalization 3-1 Chapter 3 Nodalization 3.1 Introduction 3.1.1 Chapter content This chapter focusses on establishing a rationale for, and the setting up of, the geometric representation of thermalhydraulic
More informationPAPER 333 FLUID DYNAMICS OF CLIMATE
MATHEMATICAL TRIPOS Part III Wednesday, 1 June, 2016 1:30 pm to 4:30 pm Draft 21 June, 2016 PAPER 333 FLUID DYNAMICS OF CLIMATE Attempt no more than THREE questions. There are FOUR questions in total.
More informationRichardson Extrapolated Numerical Methods for Treatment of One-Dimensional Advection Equations
Richardson Extrapolated Numerical Methods for Treatment of One-Dimensional Advection Equations Zahari Zlatev 1, Ivan Dimov 2, István Faragó 3, Krassimir Georgiev 2, Ágnes Havasi 4, and Tzvetan Ostromsky
More informationCFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate
International Journal Of Engineering Research And Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 13, Issue 9 (September 2017), PP.12-19 CFD Analysis for Thermal Behavior of Turbulent
More informationThermal Analysis Contents - 1
Thermal Analysis Contents - 1 TABLE OF CONTENTS 1 THERMAL ANALYSIS 1.1 Introduction... 1-1 1.2 Mathematical Model Description... 1-3 1.2.1 Conventions and Definitions... 1-3 1.2.2 Conduction... 1-4 1.2.2.1
More informationz x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.
Chapter 2 First order PDE 2.1 How and Why First order PDE appear? 2.1.1 Physical origins Conservation laws form one of the two fundamental parts of any mathematical model of Continuum Mechanics. These
More informationMAE 222 Mechanics of Fluids Final Exam with Answers January 13, Give succinct answers to the following word questions.
MAE 222 Mechanics of Fluids Final Exam with Answers January 13, 1994 Closed Book Only, three hours: 1:30PM to 4:30PM 1. Give succinct answers to the following word questions. (a) Why is dimensional analysis
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More informationFlood Routing by the Non-Linear Muskingum Model: Conservation of Mass and Momentum
Archives of Hydro-Engineering and Environmental Mechanics Vol. 56 (29), No. 3 4, pp. 121 137 IBW PAN, ISSN 1231 3726 Flood Routing by the Non-Linear Muskingum Model: Conservation of Mass and Momentum Dariusz
More informationNUMERICAL SOLUTION OF TWO-REGION ADVECTION-DISPERSION TRANSPORT AND COMPARISON WITH ANALYTICAL SOLUTION ON EXAMPLE PROBLEMS
Proceedings of ALGORITMY 2002 Conference on Scientific Computing, pp. 130 137 NUMERICAL SOLUTION OF TWO-REGION ADVECTION-DISPERSION TRANSPORT AND COMPARISON WITH ANALYTICAL SOLUTION ON EXAMPLE PROBLEMS
More informationUnit IV Derivatives 20 Hours Finish by Christmas
Unit IV Derivatives 20 Hours Finish by Christmas Calculus There two main streams of Calculus: Differentiation Integration Differentiation is used to find the rate of change of variables relative to one
More informationUnit IV Derivatives 20 Hours Finish by Christmas
Unit IV Derivatives 20 Hours Finish by Christmas Calculus There two main streams of Calculus: Differentiation Integration Differentiation is used to find the rate of change of variables relative to one
More informationComputational Astrophysics
16 th Chris Engelbrecht Summer School, January 2005 3: 1 Computational Astrophysics Lecture 3: Magnetic fields Paul Ricker University of Illinois at Urbana-Champaign National Center for Supercomputing
More informationAnnouncements. Topics: Homework: - sections 4.5 and * Read these sections and study solved examples in your textbook!
Announcements Topics: - sections 4.5 and 5.1-5.5 * Read these sections and study solved examples in your textbook! Homework: - review lecture notes thoroughly - work on practice problems from the textbook
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 informationMIKE 21 & MIKE 3 Flow Model. Advection-Dispersion Module. Scientific Documentation
MIKE 21 & MIKE 3 Flow Model Advection-Dispersion Module Scientific Documentation MIKE 2017 DHI headquarters Agern Allé 5 DK-2970 Hørsholm Denmark +45 4516 9200 Telephone +45 4516 9333 Support +45 4516
More information5. Advection and Diffusion of an Instantaneous, Point Source
1 5. Advection and Diffusion of an Instantaneous, Point Source In this chapter consider the combined transport by advection and diffusion for an instantaneous point release. We neglect source and sink
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationIntegrating tracer with remote sensing techniques for determining dispersion coefficients of the Dâmbovita River, Romania
Integrated Methods in Catchment Hydrology Tracer, Remote Sensing and New Hydrometric Techniques (Proceedings of IUGG 99 Symposium HS4, Birmingham, July 999). IAHS Publ. no. 258, 999. 75 Integrating tracer
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationSelf-similar solutions for the diffraction of weak shocks
Self-similar solutions for the diffraction of weak shocks Allen M. Tesdall John K. Hunter Abstract. We numerically solve a problem for the unsteady transonic small disturbance equations that describes
More informationLARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS
The 6th ASME-JSME Thermal Engineering Joint Conference March 6-, 3 TED-AJ3-3 LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS Akihiko Mitsuishi, Yosuke Hasegawa,
More informationKeywords: Large-eddy simulation, Turbulent coherent structure, Four quadrant analysis, Integral scale
The Eighth Asia-Pacific Conference on Wind Engineering, December 4, 3, Chennai, India NUMERICAL ANALYSIS OF THE MOMENTUM TRANSPORT AND TEMPORAL AND SPATIAL SCALES OF TURBULENT COHERENT STRUCTURES IN THE
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 informationCalculus I Midterm Exam. eftp Summer B, July 17, 2008
PRINT Name: Calculus I Midterm Exam eftp Summer B, 008 July 17, 008 General: This exam consists of two parts. A multiple choice section with 9 questions and a free response section with 7 questions. Directions:
More informationScientific Computing I
Scientific Computing I Module 10: Case Study Computational Fluid Dynamics Michael Bader Winter 2012/2013 Module 10: Case Study Computational Fluid Dynamics, Winter 2012/2013 1 Fluid mechanics as a Discipline
More informationThis chapter focuses on the study of the numerical approximation of threedimensional
6 CHAPTER 6: NUMERICAL OPTIMISATION OF CONJUGATE HEAT TRANSFER IN COOLING CHANNELS WITH DIFFERENT CROSS-SECTIONAL SHAPES 3, 4 6.1. INTRODUCTION This chapter focuses on the study of the numerical approximation
More informationTheoretical and Experimental Studies on Transient Heat Transfer for Forced Convection Flow of Helium Gas over a Horizontal Cylinder
326 Theoretical and Experimental Studies on Transient Heat Transfer for Forced Convection Flow of Helium Gas over a Horizontal Cylinder Qiusheng LIU, Katsuya FUKUDA and Zheng ZHANG Forced convection transient
More informationFundamentals of Fluid Dynamics: Waves in Fluids
Fundamentals of Fluid Dynamics: Waves in Fluids Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/ tzielins/ Institute
More informationDealing with Sedimental Transport Over Partly Non-Erodible Bottoms
Utah State University DigitalCommons@USU International Junior Researcher and Engineer Workshop on Hydraulic Structures Jun 17th, 12:00 AM - Jun 20th, 12:00 AM Dealing with Sedimental Transport Over Partly
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS NUMERICAL FLUID MECHANICS FALL 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139 2.29 NUMERICAL FLUID MECHANICS FALL 2011 QUIZ 2 The goals of this quiz 2 are to: (i) ask some general
More informationUNIT I FLUID PROPERTIES AND STATICS
SIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayanavanam Road 517583 QUESTION BANK (DESCRIPTIVE) Subject with Code : Fluid Mechanics (16CE106) Year & Sem: II-B.Tech & I-Sem Course & Branch:
More informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More informationOn the evaluation quadratic forces on stationary bodies
On the evaluation quadratic forces on stationary bodies Chang-Ho Lee AMIT Inc., Chestnut Hill MA, USA June 9, 006 Abstract. Conservation of momentum is applied to finite fluid volume surrounding a body
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationStable numerical methods for diffusive flows
WATERCOMP 89 The First Australian Conference on Technical Computing in the Water Industry, Melbourne 30 May-1 June 1989, pp 185-189, Retyped in 2011. Stable numerical methods for diffusive flows J.D. Fenton
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 informationSIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE CAVITY FLOWS
ICAS 2000 CONGRESS SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE CAVITY FLOWS H Yao, R K Cooper, and S Raghunathan School of Aeronautical Engineering The Queen s University of Belfast, Belfast BT7 1NN,
More informationH3: Transition to Steady State Tidal Circulation
May 7, Chapter 6 H3: Transition to Steady State Tidal Circulation Contents 6. Problem Specification............................. 6-6. Background................................... 6-6.3 Contra Costa Water
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 informationPollution. Elixir Pollution 97 (2016)
42253 Available online at www.elixirpublishers.com (Elixir International Journal) Pollution Elixir Pollution 97 (2016) 42253-42257 Analytical Solution of Temporally Dispersion of Solute through Semi- Infinite
More informationThree-dimensional Lagrangian Tracer Modelling in Wadden Sea Areas. Diploma thesis by Frank Wolk Carl von Ossietzky University Oldenburg
Three-dimensional Lagrangian Tracer Modelling in Wadden Sea Areas Diploma thesis by Frank Wolk Carl von Ossietzky University Oldenburg Hamburg, October 3, 23 i Contents List of Figures List of Tables 3
More informationA Two-Level Time-Stepping Method for Layered Ocean Circulation Models: Further Development and Testing
A Two-Level Time-Stepping Method for Layered Ocean Circulation Models: Further Development and Testing Robert L. Higdon Department of Mathematics Oregon State University Corvallis, Oregon 97331-4605 This
More informationChapter 8 STATISTICAL ASSESSMENT OF NUMERICAL DIFFUSION
84 hapter 8 STATISTIAL ASSESSMENT OF NUMERIAL DIFFUSION The assessment of the performance of high-order numerical methods regarding the numerical diffusion introduced in the solution of the scalar transport
More informationExperimental and Numerical Investigation of Two- Phase Flow through Enlarging Singularity
Purdue University Purdue e-pubs International Refrigeration and Air Conditioning Conference School of Mechanical Engineering 212 Experimental and Numerical Investigation of Two- Phase Flow through Enlarging
More informationUNIVERSITY of LIMERICK
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4607 SEMESTER: Autumn 2012-13 MODULE TITLE: Introduction to Fluids DURATION OF
More informationA note on benchmarking of numerical models for density dependent flow in porous media
Advances in Water Resources 29 (2006) 1918 1923 www.elsevier.com/locate/advwatres A note on benchmarking of numerical models for density dependent flow in porous media B. Ataie-Ashtiani *, M.M. Aghayi
More information1 Introduction. J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods
J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods 1 Introduction Achieving high order time-accuracy in the approximation of the incompressible Navier Stokes equations by means of fractional-step
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 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 informationNumerical Simulation of the Fluid Dynamics of 3D Free Surface Flows with an Upwind Convection Scheme
Numerical Simulation of the Fluid Dynamics of 3D Free Surface Flows with an Upwind Convection Scheme Laís Corrêa, Valdemir Garcia Ferreira Depto. de Matemática Aplicada e Estatística, ICMC, USP, 13566-590,
More informationMethod of Lines. Received April 20, 2009; accepted July 9, 2009
Method of Lines Samir Hamdi, William E. Schiesser and Graham W. Griffiths * Ecole Polytechnique, France; Lehigh University, USA; City University,UK. Received April 20, 2009; accepted July 9, 2009 The method
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 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 informationReservoir Oscillations with Through Flow
American Journal of Environmental Sciences 3 (): 37-42, 27 ISSN 553-345X 27 Science Publications Reservoir Oscillations with Through Flow A. A. Khan 28 Lowry Hall, epartment of Civil Engineering, Clemson
More informationSplit explicit methods
Split explicit methods Almut Gassmann Meteorological Institute of the University of Bonn Germany St.Petersburg Summer School 2006 on nonhydrostatic dynamics and fine scale data assimilation Two common
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 2014 15 FLUID DYNAMICS - THEORY AND COMPUTATION MTHA5002Y Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
More information