ABSTRACT We propose a flux vector splitting (FVS) for the solution of film flows radially spreading on a flat surface created by an impinging jet usin
|
|
- Rosaline Hunter
- 5 years ago
- Views:
Transcription
1 Analysis of Thin Film Flows Using A Flux Vector Splitting Λ J. Rafael Pacheco y Departamento de Ingenier a Mecánica Instituto Tecnológico de Monterrey Monterrey, N.L. 89, México Arturo Pacheco-Vega z Department of Aerospace and Mechanical Engineering University of Notre Dame Notre Dame, IN 55 Printed September 19, 00 Λ Address correspondence to ASU P.O. Box 10, Tempe AZ y rpacheco@asu.edu. Assistant Professor, ASME member z apacheco@nd.edu. Assistant Professor Visiting. ASME member
2 ABSTRACT We propose a flux vector splitting (FVS) for the solution of film flows radially spreading on a flat surface created by an impinging jet using the shallow-water approximation. The governing euations along with the boundary conditions are transformed from the physical to the computational domain and solved in a rectangular grid. A first-order upwind finite difference scheme is used at the point of the shock while a second-order upwind differentiation is applied elsewhere. Higher-order spatial accuracy is achieved by introducing a MUSCL approach. Three thin film flow problems 1) one-dimensional dam break problem, ) radial flow without jump, and ) radial flow with jump, are investigated with emphasis in the prediction of hydraulic jumps. Results demonstrate that the method is useful and accurate in solving the shallow water euations for several flow conditions. NOMENCLATURE ^A; ^B ~^A; ~^B ^A + ; ^B + ^A ; ^B a C Jacobian matrices modified Jacobian matrices Jacobians corresponding to positive eigenvalues Jacobians corresponding to negative eigenvalues speed of propagation [m/s] Chezy constant E; F; Q; S vectors in the physical domain ^E; ^F; ^Q; ^S vectors ^E + ; ^F + ^E ; ^F ^E ; ^F fl e z Fr in the computational domain Split flux vector corresponding to positive eigenvalues Split flux vector corresponding to negative eigenvalues matrices of corresponding eigenvectors unit vector in the vertical direction Froude number g gravitational acceleration [m/s ] ii
3 h J N film height [m] Jacobian Number of cells p pressure [N/m ] Q i Q w r i r o components of the depth averaged variables [m/s] volumetric flow rate [m /s] radial location of initial measurements [m] outer radius of the disk [m] t; fi time [s] U; V depth averaged velocities [m/s] U c ;V c u characteristic velocities [m/s] Cartesian velocity vector [m/s] u; v; w velocity components in the x, y, z directions [m/s] x; y; z Cartesian coordinates [m] z b z i reference height [m] initial film height in experiments [m] Greek symbols fl; ffi; ; r ffl ffi ^Λ; ^ ^Λ + ; ^ + ^Λ ; ^ ~^Λ; ~^ ν eigenvalues of the Jacobian [m/s] finite difference operators error nonlinear limiting parameter diagonal matrices of eigenvalues diagonal matrices of positive eigenvalues diagonal matrices of negative eigenvalues diagonal matrices of modified eigenvalues kinematic viscosity [m /s] iii
4 ρ density of the fluid [kg/m ] 1 ; weight function fi bx,fi by bottom shear stresses [N/m ] ο; ; curvilinear coordinates [m] Subscripts and superscripts i; j; k index for the Cartesian coordinates x; y; z max n x maximum index of time step x direction first derivative in the x direction y y direction first derivative in the y direction z z direction + positive eigenvalues negative eigenvalues Abbreviations CFL cfs FVS FDS MUSCL Courant-Friedrich-Levy number cubic feet per second flux vector splitting flux difference splitting monotone upstream-centered scheme for conservation laws iv
5 INTRODUCTION A phenomenon of interest in free-surface film flows under certain conditions is the formation of a sudden discontinuity in the depth of a flowing liuid; i.e. a hydraulic jump, as the flow in the radial direction makes a transition from supercritical to a subcritical regime. Significant back-flow may be present at the jump location with a corresponding loss of mechanical energy due to the abrupt change in depth. In general, the region of increased film thickness corresponds to a zone of substantially amplified turbulence. From a thermal standpoint, the location of the hydraulic jump is critical because of the degraded transport characteristics that exist there. This region is unsteady, and depends heavily on the drainage configuration downstream. Thin film flows can be characterized as shallow-water flows in which the vertical dimension is much smaller than any typical horizontal length scale. Important applications of the shallow-water euations include the prediction of weather, simulation of tidal flows, storm surges, river flows, and dam-break waves. For engineering applications numerical solution of the governing euations is a useful tool for aiding in the prediction of local flow properties of the fluid. From a numerical standpoint, a major difficulty in analyzing the hydraulic jump stems from the discontinuity of the fluid depth across the jump, which must be accurately captured by the numerical method. Due to the similarity between the depth-averaged euations and the two-dimensional compressible Navier-Stokes euations, the methods developed to solve the compressible flow may also be used to study the shallow water euations and therefore to analyze hydraulic jumps in thin flows. Several finite-difference schemes have been developed for the simulation of gas dynamics euations and are available in the literature. The main difference among them is the way they address the problem of shock waves formation, which is of fundamental importance in the overall accuracy of the calculations. These techniues are based either on centraldifference formulations [1, ], or on upwind-difference discretizations. Though the earlier 1
6 schemes give fairly accurate results in the case of smooth and weak-shocks, the precision given by the latter methods, even for problems containing strong shock waves, is superior []. Representatives of upwind schemes are those by Roe [] and Yee and Harten [5], based on flux difference splitting (FDS) techniues, and those by Steger and Warming [], van Leer [], and Walters et al. [8] structured on the flux vector splitting (FVS) approach. Since the shallow-water euations have the property of being non-homogeneous functions of the primitive variables, the application of FVS schemes to the solution of these euations has been limited, mainly due to the difficulty in finding the proper split of the flux vectors. However, a number of numerical studies for the shallow-water euations using schemes based on the FDS method have been carried out. Examples include the study of flow inachannel of infinite width [9], the simulation of open channel flows coupled with flux limiters [10], and the analysis of two-dimensional free-surface flow euations employing a high-order Godunovtype scheme based on MUSCL variable extrapolation and slope limiters [11], among others. In the current work we are interested in the analysis of free-surface thin film flows, for which ahydraulic jump may exist, using the shallowwater euations as mathematical model. To this end, we use a finite-difference scheme based on a new split of the flux-vector to solve the governing euations. A first-order upwind differentiation is enforced at the location of the shock whereas a higher order upwind differentiation is employed elsewhere. A straight forward formulation of the governing euations is presented followed by a discussion of the split of the flux vector that accounts for the inhomogeneity of the shallow-water model. The effects of the FVS on the accuracy and consistency of the spatial discretization are studied next. Finally, three thin film problems 1) one-dimensional dam break problem, ) radial flow without jump, and ) radial flow with jump, are investigated with emphasis in the capability of the techniue to capture discontinuities in the flow field. MATHEMATICAL MODEL GOVERNING EQUATIONS
7 We present the governing Navier-Stokes and continuity euations for incompressible fluid with constant r u = 0 + u ru = rp + νr u + g e z () where e z is the unit vector in the vertical direction; u = (u; v; w) represents the Cartesian velocity vector components; p is the pressure; ν is the kinematic viscosity; ρ is the density and g is the gravitational acceleration. In the following, the physical space will be denoted by (x; y; z), and the computational space by (ο; ; ). Assuming that vertical acceleration, wind stresses and Coriolis effects are negligible, the bottom slope is small, and considering that the friction losses can be computed using uniform steady flow formulas [1], the set of euations (1) () can be vertically averaged to obtain the two-dimensional unsteady shallow flow euations. The depth-averaged euations of motion can be written in conservation law form as where Q = h hu hv 5 ; E = hu hu ; F = = S; () hv huv hv + gh 5 ; S = 0 fi bx ρ fi by ρ 5 : () Here h is the water depth measured vertically, U and V are the Cartesian depth-averaged velocity components, z b is the channel bottom elevation, t is time, g is the gravitational acceleration, and fi bx and fi by are the bottom shear stress components. The bottom stress is an unknown uantity in the shallow water euations and has to be expressed in terms of the other variables in order to close" the system. In the current formulation the effective stress components are neglected and the shear stress on the bottom is approximated using the Chezy formulas fi bx = g C U(U + V ) 1= ; fi by = g C V (U + V ) 1= ; (5)
8 where C is the Chezy constant [1]. The advantage of this approach (assuming the pressure distribution to be hydrostatic), is that leads to a considerable simplification in both the mathematical formulation and the numerical solution without degrading the accuracy of the overall solution, for most engineering applications [1]. E. () can be written in generalized coordinates = ^S: () where ^Q represents the standard form of the conserved variables. Because of the particular structure of these euations, it is possible to = = ; () where ^A and ^B are the flux Jacobian matrices, i.e. ^Aij ^Ei =@ ^Qj and ^Bij ^Fi =@ ^Qj. The elements of Es. ()-() are shown in Appendix A. The system of euations () is hyperbolic and therefore the direction of propagation of the information is determined by the sign of the eigenvalues of the Jacobian matrices, ^A and ^B, which are all real. For ^A and ^B a similarity transformation exists such that ^A = ^E ^Λ^E 1 and ^B = ^F ^ ^F 1 fl fl, where ^Λ and ^ are the diagonal matrices of eigenvalues, and ^E and ^F fl the associated eigenvectors. If we define U c = ο x U + ο y V, V c = x U + y V and a = p gh, the characteristic directions in the ο- and -coordinates respectively can be written as: 1 = U c ; ; = U c ± a οx + οy (8) fl 1 = V c ; fl ; = V c ± a x + y (9) where U c and V c are the characteristic velocities. BOUNDARY CONDITIONS The system of euations () cannot be solved without assigning appropriate initial and boundary conditions. Forahyperbolic system the reuired number of boundary conditions is related to the behavior of the characteristics. At any point on the boundary one must specify
9 as many conditions as there are characteristic planes entering the region. The characteristic directions of propagation are given in euations (8) (9). Considering the radial velocity as positive if the flow is outward, we mayhave the situations shown in Table 1. Cases 1 and, in the table, are concerned with supercritical flow, Fr > 1 where Fr p =U= gh is the Froude number. Supercritical flows do not occur very often. For the more freuent cases, and (Fr < 1), we must specify either one or two boundary conditions depending on whether the flow is directed outward or inward. For supercritical flow (Fr > 1), the boundary conditions used were: (a) Upstream boundary: the longitudinal velocity U, transverse velocity V and the depth h are specified and do not change during computations, and (b) Downstream boundary: all the variables, i.e. U, V, and h are extrapolated from the interior points. For subcritical flow (Fr < 1), the corresponding boundary conditions are: (a) Upstream boundary: the depth h is specified and the longitudinal velocity U is extrapolated from the interior points, and (b) Downstream boundary: the velocities U and V are specified and the depth h is extrapolated from the interior points. SPLITTING SCHEME The splitting of the flux vectors ^E and ^F, in E. (), can be accomplished in a number of ways. The most common choice is in accord to the sign of the eigenvalues of the system, where the flux vectors can be split into sub-vectors. Each sub-vector is associated with the positive or negative eigenvalues of the flux matrix Jacobian. Thus, ^A can be split as ^A = ^A + + ^A where ^A + = ^E ^Λ+ ^E 1 and ^A = ^E ^Λ ^E 1, where the elements of ^Λ + and ^Λ are the positive and negative eigenvalues respectively. Similarly ^B = ^B + + ^B, where ^B + = ^F fl ^ + ^F 1 fl and ^B = ^F fl ^ ^F 1 fl. In the case of the Euler euations, the fluxes associated with these Jacobians are obtained from the property of homogeneity of the flux vector ^E( ^Q). If the flux vector ^E is an homogeneous function of degree one in ^Q, one has ^E = ^A ^Q. In a similar fashion, ^F 5
10 can be written as ^F = ^B ^Q. The property of homogeneity was used by Steger and Warming [] to split the flux vector into forward and backward contributions. The eigenvalues of the Jacobian matrix of the full flux were split into non-negative and non-positive groups. An alternate flux splitting with continuously differentiable flux contributions that lead to smoother solutions at sonic points was introduced in van Leer []. An important remark to be made is that the concept of flux splitting as previously defined, is totally dependent on the fact that the fluxes are homogeneous functions of degree one in ^Q. Although the system of euations () is hyperbolic and has some similarities with the Euler euations for compressible flows, there are some differences: the Jacobian matrices ^A and ^B of the shallow-water euations are not homogeneous functions of degree one in ^Q. Thence, the flux vectors of E. () cannot be expressed as ^E = ^A ^Q, and ^F = ^B ^Q. It is possible to use the FVS method to solve the Shallow-Water euations by splitting the flux vectors ^E and ^F in E. () according to ^E = ^E + + ^E and ^F = ^F + + ^F, where ^E + = ~^A+ ^Q and ^E = ~^A ^Q; (10) and ^F + = ~^B+ ^Q and ^F = ~^B ^Q: (11) To this end, we introduce the modified Jacobian matrices ~^A and ~^B, where ~^A = ^E ~^Λ^E 1 (1) with the diagonal matrix of eigenvalues ~^Λ = ο xq +ο yq ο x Q +ο y Q +1= p g(ο x +ο y ) ο x Q +ο y Q 1= p g(ο x +ο y ) 5 ; (1) and the corresponding expression for the modified Jacobian ~^B ~^B = ^F fl ~^ ^F 1 fl ; (1)
11 with the diagonal matrix of eigenvalues ~^ = xq + yq x Q + y Q +1= p g( x + y ) x Q + y Q 1= p g( x + y ) 5 : (15) It is now possible to obtain a well-behaved solution by splitting the spatial derivatives in two parts associated with the characteristic speeds. The flux vectors can be split into sub-vectors, such that each sub-vector is associated with positive or negative eigenvalues of the modified Jacobian matrix. Therefore, ~^A can be split according to ~^A = ~^A+ + ~^A where ~^A + = ^E ~^Λ+ ^E 1, and ~^A = ^E ~^Λ ^E 1. The elements of the diagonal matrices ~^Λ+ and ~^Λ are the positive and negative eigenvalues of the modified Jacobian matrix. Similarly, we write for the -coordinate ~^B = ~^B+ + ~^B where ~^B+ = ^Ffl ~^ + ^F 1 fl and ~^B = ^Ffl ~^ ^F 1 fl. ~^Λ and ~^ were modified from the original eigenvalues of the Jacobians ^A and ^B (see Appendix A) to introduce a new splitting of the flux vectors ^E and ^F. In a later section we will discuss the effect of modifying the characteristic speed on the consistency and accuracy of the scheme. DISCRETIZATION SCHEMES The governing euations () are solved by marching in time from some known set of initial conditions. The general form of the time-difference approximation is given by ^Q n+1 + ( 1 fi ffi 1+ ffiο fi 1+ ψ ffi ffiο h ( ^A+ + ^A ) n ^Q n+1i + ffi h^e+ + ^E i + ffi ffi ffi h^f+ + ^F i! n h ( ^B+ + ^B ) n ^Q n+1i) = (1 + 1) fi 1+ ^Sn 1 fi 1+ ^Sn ^Qn ^Q n 1 (1) where the difference approximations for ffi(^e + + ^E )=ffiο are shown in the next section. For the numerical simulations presented here, an explicit scheme is chosen ( 1 =0and =0) to solve euation (1).
12 SPATIAL DIFFERENCING SCHEME The upwind differencing of the convective terms is based on the flux splitting techniue, where the spatial derivatives of ^E + and ^E are approximated with backward and forward difference operators respectively. Combining the standard first- and second-order backward/forward differences the spatial derivatives can be expressed as ψ! ffi ^E+ (Q) ffiο i = ^E+ i ^E + i 1 ο + ffi i ^E+ i ^E + i 1 ο ffi i 1 ^E+ i 1 ^E + i ο (1) and ψ! ffi ^E (Q) ffiο i = ^E i+1 ^E i ο Euation (1) can also be written as with ψ! ffi ^E+ ffiο ^E + i+1= = ^E + i + ffi+ i i = ^E i+1 ^E i ο ^E + i+1= ^E + i 1= ο ffi+ i+1 ^E i+ ^E i+1 : (18) ο (19) ffi i (^E + i ^E + i 1); (0) which is euivalent to obtaining the value of ^E at the center of the cell followed by a forward extrapolation. The parameter ffi appearing in euations (1)-(18) limits the higher-order terms by acting as a nonlinear switch. ffi enhances shock capturing by dropping from a second-order to a first-order scheme near discontinuities so that oscillations in large gradient regions such as in the vicinity of a shock are suppressed. The spatial differencing is secondorder for ffi = 1 and first-order for ffi =0. The value of ffi is determined locally from ffi = " # n k οq i r ο Q i k + ffl (1) ( οq i ) +(r ο Q i ) + ffl where ffl is a small number O(10 ) [8], and n is a positive number (0 < n» 1) used to capture the shock profile. The operators ο and r ο are the standard forward and backward difference operators respectively. A spatial-discretization alternative is first, to extrapolate nodal point values of ^Q toward the face of the cell, and then to obtain ^E ± at the interfaces. This is implemented as a flux 8
13 balance across the cell and leads to the following approximation ψ! ffi ^E( ^Q) ffiο i = h^e+ ( ^Q i+1= ) ^E + ( ^Q i 1= )+^E ( ^Q + i+1= ) ^E ( ^Q + i 1= )i ο () where ^E ± ( ^Q ) denotes ^E ± evaluated at ^Q, and ^Q i+1=;j = ^Q i;j + ffi i [(1»ffi i )r ο +(1+»ffi i ) ο] ^Q i;j () ^Q + i+1=;j = ^Q i+1;j ffi+ i+1 [(1 +»ffi+ i+1) ο +(1»ffi + i+1)r ο ] ^Q i+1;j : () The aforementioned method, developed by van Leer [], is often referred as the MUSCL-type approach (Monotone Upstream-centered Schemes for Conservation Laws). In Es. ()- (),» determines the spatial accuracy of the difference approximation. For fully upwind differencing, the choice» = 1 is used whereas» = 1= corresponds to the upwind-biased third-order scheme [15, 1]. Analogous expressions to (1)-() may be obtained for the gradient ofthefluxes ffi(f + + ^F )=ffi. A comparison of numerical results for the dam-break problem, obtained with the discretized forms of Es. (1)-(18), and E. () will be presented in the section of numerical experiments. ACCURACY AND CONSISTENCY ANALYSES We now analyze the impact on the order of accuracy and consistency of the scheme due to the modification of the characteristic speeds (eigenvalues) of the modified Jacobian matrices ~^A and ~^B of our flux-splitting. Consider for simplicity the one-dimensional shallow-water ^E discretized in space and time with the following first-order scheme as ^Q n+1 i ^Q n i = fi h (^E+ i ^E + i 1)+(^E i+1 ^E i ) i n ο = fi ο =0 (5) h (^E+ i + ^E i+1) (^E + i 1 + ^E i ) i n : () 9
14 On expressing the absolute value of ^E as j^ej = ^E + ^E = j ~^Aj ^Q, where j ~^Aj = ^E j ~^Λj^E 1, and using ^E = ~^A+ ^Q + ~^A ^Q with ~^A+ ^Q = ^E+ and ~^A ^Q = ^E, we can recast E. () as ^Q n+1 i ^Q n i = t ^Ei+1 ^E n fi n i 1 + j^eji+1 j^ej i + j^ej i 1 ; () ο ο to show an artificial viscosity term proportional to the second derivative of the absolute flux value j^ej. It can be noticed from E. () that the consistency of the scheme is maintained. After a Taylor expansion of (j^ej i+1 j^ej i ) and (j^ej i j^ej i 1), we can transform E. () in a form that explicits the absolute value of the modified Jacobian j ~^Aj; e.g., j^ej i+1 j^ej i = j ~^Aji+ 1 ^Qi+1 ^Q i 1 + O( ο) ^Qi+1 ^Q i 1 ; (8) j^ej i j^ej i 1 = j ~^Aji 1 ^Qi+1 ^Q i 1 + O( ο) ^Qi+1 ^Q i 1 : (9) Using Es. (8) (9), E. () becomes, to the first order accuracy ^Q n+1 i ^Q n i = t» ^Ei+1 i 1 ^E n fi + j ~^Aji+ 1 ^Qi+1 ^Q i ο ο j ~^Aji ^Qi 1 ^Q n i 1 + O( ^Q ): (0) Euations () and () are euivalent formulations of the first-order conservativeupwind scheme (5) for smooth variations of ^Q. One can write E. () in terms of the numerical flux, as with ^E i+ 1 ^Q n+1 i ^Q n i = fi ^Ei+ 1 ^E n ο i 1 (1) = ^E + i + ^E i+1, or up to second order accuracy as ^E i+ 1 = 1 ^Ei + ^E 1 i+1 ~^A i+ 1 ^Qi+1 ^Q i + O( ^Q ): () It can be observed that the modification of the characteristic speed does not affect the order of accuracy of the scheme, as it is part of an artificial dissipation term [1]. NUMERICAL EXPERIMENTS 10
15 Three problems are solved in order to examine the proposed FVS. The first is the analysis of the well-known one-dimensional dam-breaking, the second corresponds to a radial flow within a supercritical region in the domain for which no hydraulic jump exists, and the third is a radial flow in the presence of supercritical and subcritical regions which induce a hydraulic jump. ONE-DIMENSIONAL DAM-BREAKING PROBLEM This problem consists of a dam located at the origin (x = 0) in a rectangular channel, that extends to ±1, in which fluid is initially at rest on both sides. The x-coordinate is in the downstream direction of the channel. The initial height of the fluid on the left (x <0) is 1 m while the height on the right (x > 0) is 0.8 m. For this problem, [1] developed an exact solution based on the theory of characteristics to analyze the behavior of the fluid after the dam breaks. The numerical simulations were carried out using the proposed FVS along with four different schemes: a first-order upwind, a second-order upwind, and two different schemes that use slope limiters. Figures 1(a) 1(d) illustrate comparisons of the results obtained from the numerical procedure with the exact solution at time t = 0:10 seconds using 0 mesh points. It can be observed from the figures that, in general, the discontinuities are captured by all the schemes, but some differences in accuracy are evident. For instance, Figure 1(a) shows that the first-order upwind scheme has good properties of monotonicity, however it suffers from severe inaccuracies because of the high level of truncation error. On the other hand, the second-order upwind scheme, shown in Figure 1(b), yields low levels of numerical dissipation, but steep gradients produce an oscillating behavior due to its numerical dispersion. In Figure 1(c) the slope limiter techniue based on extrapolation of the flux vector ^E shows a substantial reduction in the amplitude of the oscillations at the discontinuities but cannot be completely suppressed. Finally, the results obtained from coupling FVS with the MUSCL techniue indicate that the method is able to capture the discontinuities in a sharp way, avoiding any unrealistic oscillation as illustrated in Figure 1(d). 11
16 The maximum relative error in h, ffl max, at the points of discontinuity and time of 0:10 s is plotted in Figure as a function of mesh refinement, for the MUSCL scheme. Similar results were obtained for the velocity U. A uniform grid was used in the computation. The time step was refined proportional to the grid spacing. Figure shows that the present method is first order accurate at the shock locations. The maximum relative error in h at time 0:10 s away from the shock is presented in Table as function of mesh refinement. The experimentally established orders or accuracy are also shown in this table as function of mesh refinement. In our computation, we used» = 1= which corresponds to the upwind-biased third-order scheme [15]. RADIAL FLOW WITHOUT JUMP We simulate the radial flow with the proposed flux vector splitting method coupled with the MUSCL techniue and compare the results with the experiments of Sadler and Higgins [18]. Alayout of their experimental setup is shown in Figure ; the reported geometrical parameters and flow rates are presented in Table as case 1. It consists of a : m (8 ft) suare sheet of Plexiglas mounted horizontally on a 15: cm ( in) diameter vertical pipe. Both Plexiglas and pipe were fixed to a structure. Water was discharged from a nozzle, from below, onto a 0: cm (8 in) diameter disk located at a height of5mm(=1 in) above the sheet of Plexiglas. The water was then collected through a series of drains arranged in a circular fashion located at a distance of 8: cm ( in) from the center. In the experiment, the depth of the flow was measured with a point gauge mounted on an arm under steady state conditions. For simulation purposes, we use a flow rate of Q w =0:0098 m /s (0:5 cfs), as reported by Sadler and Higgins [18], and a value of the Chezy constant C = 80 [19] which is in the range commonly used for practical engineering calculations. It must be noted that a simplification of the problem can be obtained since the flow is axisymmetric, so that a one-dimensional version of the shallow-water euations in cylindrical coordinates may beused. Nonetheless, we consider a two-dimensional discretization of the domain in order to illustrate potential 1
17 applications of the method for non-axisymmetric flows [0]. The upstream values of h, U, and V are assigned as initial conditions for the numerical runs at all the grid points. The upstream boundary condition is of the Dirichlet type (case 1 of Table 1). Since the flow is supercritical within the domain, the boundary condition downstream corresponds to the Neumann type. Therefore the values of the variables are extrapolated using the values at the interior points. In the azimuthal direction periodic boundary conditions were implemented. The numerical scheme was tested until grid independent solutions were obtained by successively refining the mesh. Figure shows the film height for two different sets of mesh points. The FVS-based model was run until steady state was reached. A comparison between the numerical calculations obtained from the proposed FVS method and the experimental results of Sadler and Higgins [18] is presented in Figure. The water surface profiles determined from the model and the surface profile found experimentally for the current analysis are shown in this figure. The numerical results display a very good agreement with experiments. RADIAL FLOW WITH JUMP The radially spreading flow of a thin liuid film is simulated in the presence of supercritical and subcritical regions. The shock-capturing capabilities of our FVS method are tested and the results are compared with the laboratory test data obtained by [1]. A schematic of the experimental layout is shown in Figure 5. It consisted of a:1 m(ft)diameter by :5 cm (1 in) thickness circular plywood disc placed in a circular steel tank of diameter :9 m (:5 ft). A ring of 1:8 m ( ft) inside diameter, :1 m ( ft) outside diameter, and :5 cm (1 in) thickness was placed on the plywood disk to control the downstream depth and conseuently the location of the jump. A 10: cm ( in) diameter pipe was used to convey the water onto the disc. The distance between the nozzle and the disc was kept at cm (1=1 in). For the numerical simulation we use the flow conditions reported by Ahmad [1] which are presented as case in Table. Since the flow is supercritical in the upstream region, the 1
18 variables U, V, and h, once specified, remain unchanged during the iteration process. At the downstream boundary, U and V are defined whereas h is extrapolated from the interior points. Periodic boundary conditions are prescribed in the azimuthal direction. For the present case, the value of the Chezy constant was chosen as C =90[19], and the time step fi was restricted to a CFL=0.5. Results from this run are shown in Figure. The values of the speed of propagation (eigenvalues and fl ) are shown in Figure as a function of the radial position. From Figure, it is clear that the change of sign of the eigenvalue coincides with the location of the jump, as expected. It is also important to note that the location of the jump coincides at a location where the characteristics intersect. In our test, we start from an initial condition and stop the computation when the change in water depth reaches a prescribed tolerance of ffl =1: m. The computed results from the proposed scheme and the measurements of Ahmad [1] are shown in Figure for comparison. From the figure it can be observed first that the film heights from the simulations using a mesh of points, are almost identical when increasing the number of grid points in the radial-azimuthal directions to 1 11, for the same initial conditions indicating grid independence. Moreover, when comparing the numerical simulations with the experimental measurements becomes clear that the method simulates the supercritical flow accurately, and that it is able to determine the location of the jump. The sharp transition from supercritical to subcritical flow is clearly shown in Figure. In the region of transition from supercritical to subcritical flow, the computed water-depth magnitude is very close to the measured height reported in the experiment. In the subcritical region where the flow is turbulent, the difference in heights between simulations and experiments becomes larger. These differences may be due on one side to experimental uncertainties which are increased in the turbulent region, and on the other, to the energy dissipated at the location of the jump, which is not accounted for in the shallow-water formulation. CONCLUDING REMARKS 1
19 The hydraulic jump is an intriguing complex phenomenon commonly found in nature. The complexity stems from the change in the character of the flow at the jump location with the corresponding production of distorted free surfaces and flow separation. Because of its important role in engineering applications, e.g. the design of water-courses, energy dissipators, etc., it is crucial to determine accurately the characteristics of the jump such as position and thickness. The shallowwater euations are mathematically similar to the euations for compressible fluid flow. They are useful in predicting the location of shocks, at which the height of the free surface and the mean fluid velocity are discontinuous. In contrast to the compressible fluid flow euations, the shallow water euations do not posses the property of homogeneity and hence common flux vector splittings (FVS) cannot be used to solve them. In this work we have introduced a new FVS that allows the analysis of thin film flows, in particular, the formation of hydraulic jumps within the flow domain using the shallow-water theory. The proposed splitting of the flux vectors was obtained by modified Jacobian matrices that account for the inhomogeneity of the system. Details of the methodology and theoretical analyses have been included to illustrate its implementation, consistency and accuracy. The usefulness of the proposed FVS is tested by coupling the method to several finite difference schemes, two of them with slope limiters, that yield high-order accuracy away from the location of the shock. The numerical predictions have shown favorable comparisons with the theoretical solution of the dam-break problem. When the proposed FVS is used in conjunction with MUSCL, the resulting scheme is more efficient and produces accurate results in the location of the discontinuity and the height of the flow, avoiding unrealistic oscillations. An excellent agreement with laboratory experiments and analytical solutions confirmed the efficiency of the method. The techniue proposed here has been successfully applied to analyze axi-symmetric and non-axi-symmetric density currents [0]. 15
20 ACKNOWLEDGMENTS The authors would like to thank Professor Kyle Suires for helpful discussions on this work and to Professor Miguel Angel Arredondo of Universidad Iberoamericana, León, for access to computer facilities. We are grateful for the helpful comments of the anonymous referees which improved the uality ofthemanuscript. REFERENCES [1] MacCormack R. W. The effect of viscosity in hypervelocity impact cratering. AIAA Paper, pages 9 5, 199. [] Beam R. M. and Warming R. F. An implicit finite-difference algorithm for hyperbolic systems in conservation law form. J. Comput. Phys., :8 110, 19. [] Yee H. C. A class of high-resolution explicit and implicit shock-capturing methods. Technical report, von Karman Institute for Fluid Dynamics, [] Roe P. L. Approximate Riemann solvers, parameter vectors, and difference schemes. J. Comput. Phys., :5, [5] Yee H. C. and Harten A. Implicit TVD schemes for hyperbolic conservation-laws in curvilinear coordinates. AIAA J., 5:, 198. [] Steger J. L. and Warming R. F. Flux vector splitting of the inviscid gasdynamic euation with application to finite difference methods. Technical report, NASA TM-805, 199. [] Van Leer B. Flux-vector splitting for the Euler euations. In Krause E., editor, Lecture Notes in Physics., volume 10, pages Springer, Berlin, 198. [8] Walters R. W., Thomas J. L., and Van Leer B. An implicit flux-split algorithm for the compressible Euler and Navier-Stokes euations. In Zhuang F.G. and Zhu Y.L, editors, Lecture Notes in Physics., volume, pages 8 5. Springer, Berlin, 198. [9] Glaister P. Approximate Riemann solutions of the shallow water euations. J. of Hydraulic Research, :9 0, [10] Glaister P. Flux difference splitting for open-channel flows. Int. J. for Num. Meth. Fluids, 1:9 5, 199. [11] Alcrudo F. and Garcia-Navarro P. A high-resolution Godunov-type scheme in finite volumes for the D shallow-water euations. Int. J. for Num. Meth. Fluids, 1:89 505, 199. [1] Chapman R. S. and Kuo C. Y. Application of the two-euation» ffl turbulence model to a twodimensional, steady, free surface flow problem with separation. Int. J. for Num. Meth. Fluids, 5:5 8, [1] Vreugdenhill C. B. Numerical Methods for Shallow-Water Flows. Kluwer Academic Publishers, Dordrecht, Netherlands, 199. [1] Cunge J. Rapidly varying flow in power and pumping canals. In Mahmood K. and Yevjevich V., editors, Unsteady flow in open channels, pages 59 58, Littleton, CO., 195. Water Resources Publications. [15] Thomas J.L. and Walters R.W. Upwind relaxation algorigthms for the Navier-Stokes euations. AIAA Paper, pages , [1] Hirsch C. Numerical Computation of Internal and External Flows. John Wiley & Sons, Inc., New York, U.S.A., [1] Stoker J. J. Water Waves. Interscience, New York, U.S.A., 195. [18] Sadler C. D. and Higgins M. S. Radial free surfaces flows. Master's thesis, Massachussetts Institute of Technology, Cambridge, MA, 19. 1
21 [19] Younus M. and Chaudhry M. H. A depth-averaged» ffl turbulence model for the computation of free surface flow. J. of Hydraulic Research, ():15, 199. [0] Hunt J. C. R., Pacheco J. R., Mahalov A., and Fernando H. J. S. Effects of rotation and sloping terrain on fronts of density currents. J. Fluid Mech. (submitted), 00. [1] Ahmad D. Circular hydraulic jump. Master's thesis, Colorado State University, Fort Collings, CO, 19. 1
22 Case Condition Number of Number of Number of outward inward boundary characteristics characteristics conditions 1 U c < a οx + οy 0 a οx + οy» U c < 0 1 0» U c <a οx + οy 1 1 U c > a οx + οy 0 0 Table 1: Boundary conditions for the shallow-water euations. Case r i [m] r o [m] z i [m] Q w [m /s] Experiment 1 0:10 0: :0098 Sadler and Higgins [18] 0:08 0: :010 Ahmad [1] Table : Conditions from experiments for the solution of the shallow-water euations. 18
23 N ffl max [m] Order Table : Convergence test results for the Dam-break problem away from the discontinuity region. 19
24 h(x,t) x (a) Exact solution; Π First order scheme. 0
25 h(x,t) x (b) Exact solution; ffi Second order scheme. 1
26 h(x,t) x (c) Exact solution; 5 scheme based on flux extrapolation ^E with limiter.
27 h(x,t) x (d) Exact solution; scheme based on variable extrapolation ^Q (MUSCL) with limiter. Figure 1: Water depth of the dam-break problem in a rectangular channel at time t =0:10 second and 0 grid points.
28 ε max N Figure : Maximum relative error of h as a function of mesh refinement for the the dam-break problem on a rectangular channel at time t =0:10 s at the points of discontinuity.
29 z r i z i r o r Figure : Schematic of the experiment layout of Sadler and Higgins [18]. 5
30 [m] h r [m] Figure : Comparison of numerical solutions and experiments for Q w =0:0098 m /s (0:5 cfs); 5 (11 11) grid points; ffi (1 11) grid points; Π Sadler and Higgins [18].
31 z z i r i r o r Figure 5: Schematic of the experiment layout of Ahmad [1].
32 [m] h [m/s] λ,γ r Figure : Evolution of the eigenvalues and appearance of the hydraulic jump; ffi ; / fl ; height h. FVS with slope limiter solution; (1 11) grid points. [m] -50 8
33 [m] h r Figure : Comparison of water surface profile of a hydraulic jump for Q w = 0:010 m /s (0:0 cfs) between experiments and numerical solutions for different number of grid points: (11 11) grid mesh; ffi (1 11) grid mesh; Π observed by Ahmad [1]. [m] 9
34 APPENDIX A: EXPANDED VECTORS We present next the expanded form of the expressions often encountered in solving the shallow-water euations. The flux vectors are given by: ^Q = 1 J Q Q 5 = 1 J h hu hv ο x Q + ο y Q ^E = 1 ο x Q + ο y Q Q J ο x Q Q + ο y Q x Q + y Q ^F = 1 x Q + y Q Q J x Q Q + y Q ^S = 1 J 0 g Q C Q 1 g Q C Q 1 5 () + ο x g + ο y g + x g + y g 5 5 Q + Q Q + Q () (5) 5 : () The Jacobians of flux vectors are: ^A ^Q = 0 ο x ο y οxq +οy Q Q gq 1 Q 1 οx Q Q +ο y Q gq 1 Q 1 ο x Q +ο y Q ο y Q ο x Q ο x Q + ο y Q 5 : () The diagonal matrix of eigenvalues is, ^Λ = ο xq +ο yq ο x Q +ο y Q + p g(ο x +ο y ) ο x Q +ο y Q p g(ο x +ο y ) 5 (8) and the corresponding set of eigenvectors can be written as, ^E = p p Q ο y + οx (ο x +ο)gq y 1 Q (ο x +ο y) οx p p Q ο x + οy (ο x +ο)gq y 1 Q οy (ο x +ο y) (ο x +ο)gq y 1 (ο x +ο y) (ο x +ο)gq y 1 (ο x +ο y) 5 : (9) 0
35 The inverse of the left eigenvectors is expressed as: ^E 1 = Q ο y ο x Q (ο x +ο y) οy Q (ο x +ο y)g +ο x Q (ο x +ο y)g ο y Q + (ο x +ο y)g +ο x Q (ο x +ο y)g ο y (ο x +ο y) ο x (ο x +ο y)g (ο ο x x +ο y)g οx (ο x +ο y) ο y (ο x +ο y)g ο y (ο x +ο y)g 5 : (0) The Jacobians of flux vectors are: ^B ^Q = 0 x y xq + y Q Q gq 1 Q 1 x Q Q + y Q gq 1 Q 1 x Q + y Q y Q x Q x Q + y Q 5 : (1) The diagonal matrix form of the Jacobians of flux vectors can be written as ^ = xq + yq 0 0 p x Q + y Q + Q1 g( 0 x + ) y x Q + y Q p Q1 g( x + y ) 5 : () The set of eigenvectors corresponding to the eigenvalues is: ^F fl = p p Q y + x ( x + )gq y 1 Q ( x + y) x p p Q x + y ( x + )gq y 1 Q y ( x + y) ( x + )gq y 1 ( x + y) ( x + )gq y 1 ( x + y) 5 : () The inverse of the set of eigenvectors is: ^F 1 fl = y Q Q y x Q ( x + y) ( x + y)g + x Q ( x + y)g y Q + ( x + y)g + x Q y ( x + y) x ( x + y)g x ( x + y)g ( x + y)g x ( x + y) y ( x + y)g y ( x + y)g 5 () In order to show that the shallow-water euations do not have the property of homo- 1
36 geneity we write: (^E ^Λ^E 1 ) ^Q = 1 J (^F fl ^ ^F 1 fl ) ^Q = 1 J ο x Q + ο y Q ο x Q + ο y Q Q + ο x gq 1 ο x Q Q + ο y Q + ο y gq 1 x Q + y Q x Q + y Q Q + x gq 1 x Q Q + y Q + y gq (5) () which are different from () and (5).
Fluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationGODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS
GODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS ARTURO S. LEON Dept. of Civil and Envir. Engng., Univ. of Illinois at Urbana-Champaign, 2519 Hydrosystems Lab., MC-250. 205 North Mathews Av., Urbana,
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 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 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 information^Laboratorio Nacional de Computaqao Cientifica, Rua Lauro Muller 455, Rio de Janeiro, RJ , Brasil
A space-timefiniteelement formulation for the shallow water equations F.L.B. Ribeiro,' A.C. Galeao," L. Landau* "Programa de Engenharia Civil, COPPE / Universidade Federal do Rio de Janeiro, Caixa Postal
More informationTwo-Dimensional Simulation of Subcritical Flow at a Combining Junction: Luxury or Necessity?
Two-Dimensional Simulation of Subcritical Flow at a Combining Junction: Luxury or Necessity? Rabih Ghostine 1 ; Robert Mose ; José Vazquez 3 ; Abdellah Ghenaim 4 ; and Caroline Grégoire 5 Abstract: Classically,
More informationComparison of cell-centered and node-centered formulations of a high-resolution well-balanced finite volume scheme: application to shallow water flows
Comparison of cell-centered and node-centered formulations of a high-resolution well-balanced finite volume scheme: application to shallow water flows Dr Argiris I. Delis Dr Ioannis K. Nikolos (TUC) Maria
More informationForce analysis of underwater object with supercavitation evolution
Indian Journal of Geo-Marine Sciences Vol. 42(8), December 2013, pp. 957-963 Force analysis of underwater object with supercavitation evolution B C Khoo 1,2,3* & J G Zheng 1,3 1 Department of Mechanical
More informationA Multiphysics Strategy for Free Surface Flows
A Multiphysics Strategy for Free Surface Flows Edie Miglio, Simona Perotto, and Fausto Saleri MOX, Modeling and Scientific Computing, Department of Mathematics, Politecnico of Milano, via Bonardi 9, I-133
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 informationNPTEL Quiz Hydraulics
Introduction NPTEL Quiz Hydraulics 1. An ideal fluid is a. One which obeys Newton s law of viscosity b. Frictionless and incompressible c. Very viscous d. Frictionless and compressible 2. The unit of kinematic
More informationA note on critical flow section in collector channels
Sādhan ā, Vol. 26, Part 5, October 2001, pp. 439 445. Printed in India A note on critical flow section in collector channels 1. Introduction SUBHASISH DEY Department of Civil Engineering, Indian Institute
More informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More informationEigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations
NASA/CR-1998-208741 ICASE Report No. 98-51 Eigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations David L. Darmofal Massachusetts Institute of Technology, Cambridge,
More informationProjection Dynamics in Godunov-Type Schemes
JOURNAL OF COMPUTATIONAL PHYSICS 142, 412 427 (1998) ARTICLE NO. CP985923 Projection Dynamics in Godunov-Type Schemes Kun Xu and Jishan Hu Department of Mathematics, Hong Kong University of Science and
More informationAn evaluation of a conservative fourth order DNS code in turbulent channel flow
Center for Turbulence Research Annual Research Briefs 2 2 An evaluation of a conservative fourth order DNS code in turbulent channel flow By Jessica Gullbrand. Motivation and objectives Direct numerical
More informationThe Computations of Jet Interaction on a Generic Supersonic Missile
The Computations of Jet Interaction on a Generic Supersonic Missile *Jinbum Huh 1) and Seungsoo Lee 2) 1), 2) Department of Aerospace Engineering, Inha Univ., Incheon, Korea 2) slee@inha.ac.kr ABSTRACT
More informationHydraulic Modelling for Drilling Automation
Hydraulic Modelling for Drilling Automation CASA Day Harshit Bansal April 19, 2017 Where innovation starts Team: Supervisors at TU/e : W.H.A. Schilders, N. van de Wouw, B. Koren, L. Iapichino Collaborators:
More informationApplication of the relaxat ion met hod to model hydraulic jumps
Application of the relaxat ion met hod to model hydraulic jumps P. J. Montgomery Mathematics and Computer Science Program, University of Northern British Columbia, Prince George, Canada. Abstract A finite
More informationNAPC Numerical investigation of axisymmetric underexpanded supersonic jets. Pratikkumar Raje. Bijaylakshmi Saikia. Krishnendu Sinha 1
Proceedings of the 1 st National Aerospace Propulsion Conference NAPC-2017 March 15-17, 2017, IIT Kanpur, Kanpur NAPC-2017-139 Numerical investigation of axisymmetric underexpanded supersonic jets Pratikkumar
More informationVISCOUS FLUX LIMITERS
VISCOUS FLUX LIMITERS E. F. Toro Department of Aerospace Science College of Aeronautics Cranfield Institute of Technology Cranfield, Beds MK43 OAL England. Abstract We present Numerical Viscosity Functions,
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
More informationIMPLICIT NUMERICAL SCHEME FOR REGULATING UNSTEADY FLOW IN OPEN CHANNEL Mohamed. T. Shamaa 1, and Hmida M. Karkuri 2
IMPLICIT NUMERICAL SCHEME FOR REGULATING UNSTEADY FLOW IN OPEN CHANNEL Mohamed. T. Shamaa 1, and Hmida M. Karkuri 2 1 Associated Professor, Irrigation and Hydraulic Department, College of Technical Engineering,
More informationScienceDirect. A pipe network simulation model with dynamic transition between free surface and pressurized flow
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 70 ( 2014 ) 641 650 12th International Conference on Computing and Control for the Water Industry, CCWI2013 A pipe network simulation
More informationResolving the dependence on free-stream values for the k-omega turbulence model
Resolving the dependence on free-stream values for the k-omega turbulence model J.C. Kok Resolving the dependence on free-stream values for the k-omega turbulence model J.C. Kok This report is based on
More informationApplication of a Laser Induced Fluorescence Model to the Numerical Simulation of Detonation Waves in Hydrogen-Oxygen-Diluent Mixtures
Supplemental material for paper published in the International J of Hydrogen Energy, Vol. 30, 6044-6060, 2014. http://dx.doi.org/10.1016/j.ijhydene.2014.01.182 Application of a Laser Induced Fluorescence
More informationA Solution Method for the Reynolds-Averaged Navier-Stokes Equation
A Solution Method for the Reynolds-Averaged Navier-Stokes Equation T.-W. Lee Mechanical and Aerospace Engineering, Arizona State University, Tempe, AZ, 8587 Abstract- Navier-Stokes equations are difficult
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationCOURSE NUMBER: ME 321 Fluid Mechanics I 3 credit hour. Basic Equations in fluid Dynamics
COURSE NUMBER: ME 321 Fluid Mechanics I 3 credit hour Basic Equations in fluid Dynamics Course teacher Dr. M. Mahbubur Razzaque Professor Department of Mechanical Engineering BUET 1 Description of Fluid
More informationEFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL
EFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL Vahid SHOKRI 1*,Kazem ESMAEILI 2 1,2 Department of Mechanical Engineering, Sari Branch, Islamic Azad
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationWave propagation methods for hyperbolic problems on mapped grids
Wave propagation methods for hyperbolic problems on mapped grids A France-Taiwan Orchid Project Progress Report 2008-2009 Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan ISCM
More informationBasics on Numerical Methods for Hyperbolic Equations
Basics on Numerical Methods for Hyperbolic Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8,
More information3D Numerical Simulation of Supercritical Flow in Bends of Channel
3D Numerical Simulation of Supercritical Flow in Bends of Channel Masoud. Montazeri-Namin, Reyhaneh-Sadat. Ghazanfari-Hashemi, and Mahnaz. Ghaeini- Hessaroeyeh Abstract An attempt has been made to simulate
More informationClosed duct flows are full of fluid, have no free surface within, and are driven by a pressure gradient along the duct axis.
OPEN CHANNEL FLOW Open channel flow is a flow of liquid, basically water in a conduit with a free surface. The open channel flows are driven by gravity alone, and the pressure gradient at the atmospheric
More informationMULTIGRID CALCULATIONS FOB. CASCADES. Antony Jameson and Feng Liu Princeton University, Princeton, NJ 08544
MULTIGRID CALCULATIONS FOB. CASCADES Antony Jameson and Feng Liu Princeton University, Princeton, NJ 0544 1. Introduction Development of numerical methods for internal flows such as the flow in gas turbines
More informationTHE HYDRAULIC PERFORMANCE OF ORIENTED SPUR DIKE IMPLEMENTATION IN OPEN CHANNEL
Tenth International Water Technology Conference, IWTC10 2006, Alexandria, Egypt 281 THE HYDRAULIC PERFORMANCE OF ORIENTED SPUR DIKE IMPLEMENTATION IN OPEN CHANNEL Karima Attia 1 and Gamal El Saied 2 1
More informationIntensely swirling turbulent pipe flow downstream of an orifice: the influence of an outlet contraction
13 th Int. Symp. on Appl. Laser Techniques to Fluid Mechanics, Lisbon, Portugal, June 26-29, 26 Intensely swirling turbulent pipe flow downstream of an orifice: the influence of an outlet contraction Marcel
More informationInternational Engineering Research Journal
Special Edition PGCON-MECH-7 Development of high resolution methods for solving D Euler equation Ms.Dipti A. Bendale, Dr.Prof. Jayant H. Bhangale and Dr.Prof. Milind P. Ray ϯ Mechanical Department, SavitribaiPhule
More informationNUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD
NUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD V. G. Guedes a, G. C. R. Bodstein b, and M. H. Hirata c a Centro de Pesquisas de Energia Elétrica Departamento de Tecnologias
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 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 informationFEDSM COMPUTATIONAL AEROACOUSTIC ANALYSIS OF OVEREXPANDED SUPERSONIC JET IMPINGEMENT ON A FLAT PLATE WITH/WITHOUT HOLE
Proceedings of FEDSM2007: 5 th Joint ASME/JSME Fluids Engineering Conference July 30-August 2, 2007, San Diego, CA, USA FEDSM2007-37563 COMPUTATIONAL AEROACOUSTIC ANALYSIS OF OVEREXPANDED SUPERSONIC JET
More informationAn efficient numerical method for hydraulic transient computations M. Ciccotelli," S. Sello," P. Molmaro& " CISE Innovative Technology, Segrate, Italy
An efficient numerical method for hydraulic transient computations M. Ciccotelli," S. Sello," P. Molmaro& " CISE Innovative Technology, Segrate, Italy Abstract The aim of this paper is to present a new
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 informationFormation Of Hydraulic Jumps On Corrugated Beds
International Journal of Civil & Environmental Engineering IJCEE-IJENS Vol:10 No:01 37 Formation Of Hydraulic Jumps On Corrugated Beds Ibrahim H. Elsebaie 1 and Shazy Shabayek Abstract A study of the effect
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 informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationClosed duct flows are full of fluid, have no free surface within, and are driven by a pressure gradient along the duct axis.
OPEN CHANNEL FLOW Open channel flow is a flow of liquid, basically water in a conduit with a free surface. The open channel flows are driven by gravity alone, and the pressure gradient at the atmospheric
More informationFORMATION OF HYDRAULIC JUMPS ON CORRUGATED BEDS
International Journal of Civil & Environmental Engineering IJCEE-IJENS Vol: 10 No: 01 40 FORMATION OF HYDRAULIC JUMPS ON CORRUGATED BEDS Ibrahim H. Elsebaie 1 and Shazy Shabayek Abstract A study of the
More information1.060 Engineering Mechanics II Spring Problem Set 8
1.060 Engineering Mechanics II Spring 2006 Due on Monday, May 1st Problem Set 8 Important note: Please start a new sheet of paper for each problem in the problem set. Write the names of the group members
More information10.52 Mechanics of Fluids Spring 2006 Problem Set 3
10.52 Mechanics of Fluids Spring 2006 Problem Set 3 Problem 1 Mass transfer studies involving the transport of a solute from a gas to a liquid often involve the use of a laminar jet of liquid. The situation
More informationThe deposition efficiency and spatial thickness distribution of films created by Directed
Chapter 8 Vapor Transport Model Development The deposition efficiency and spatial thickness distribution of films created by Directed Vapor Deposition synthesis have been shown to be sensitive functions
More informationChapter 3 Bernoulli Equation
1 Bernoulli Equation 3.1 Flow Patterns: Streamlines, Pathlines, Streaklines 1) A streamline, is a line that is everywhere tangent to the velocity vector at a given instant. Examples of streamlines around
More informationComputation of Incompressible Flows: SIMPLE and related Algorithms
Computation of Incompressible Flows: SIMPLE and related Algorithms Milovan Perić CoMeT Continuum Mechanics Technologies GmbH milovan@continuummechanicstechnologies.de SIMPLE-Algorithm I - - - Consider
More informationAnalysis of Head Loss in Pipe System Components: Quasi-one-Dimensional Nozzle Flow
Analysis of Head Loss in Pipe System Components: Quasi-one-Dimensional Nozzle Flow Charles Ndambuki Muli School of Education, The Presbyterian University of East Africa, P.O. Box 387-00902, Kikuyu, Kenya.
More information28.2 Classification of Jumps
28.2 Classification of Jumps As mentioned earlier, the supercritical flow Froude number influences the characteristics of the hydraulic jump. Bradley and Peterka, after extensive experimental investigations,
More informationAerodynamic force analysis in high Reynolds number flows by Lamb vector integration
Aerodynamic force analysis in high Reynolds number flows by Lamb vector integration Claudio Marongiu, Renato Tognaccini 2 CIRA, Italian Center for Aerospace Research, Capua (CE), Italy E-mail: c.marongiu@cira.it
More informationPhysical Diffusion Cures the Carbuncle Phenomenon
Physical Diffusion Cures the Carbuncle Phenomenon J. M. Powers 1, J. Bruns 1, A. Jemcov 1 1 Department of Aerospace and Mechanical Engineering University of Notre Dame, USA Fifty-Third AIAA Aerospace Sciences
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 information6.1 Momentum Equation for Frictionless Flow: Euler s Equation The equations of motion for frictionless flow, called Euler s
Chapter 6 INCOMPRESSIBLE INVISCID FLOW All real fluids possess viscosity. However in many flow cases it is reasonable to neglect the effects of viscosity. It is useful to investigate the dynamics of an
More informationWell-balanced central finite volume methods for the Ripa system
Well-balanced central finite volume methods for the Ripa system R. Touma a C. Klingenberg b a Lebanese American University, Beirut, Lebanon b Würzburg University, Würzburg, Germany This paper is dedicated
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 informationGODUNOV-TYPE SOLUTIONS FOR TRANSIENT FLOWS IN SEWERS
GODUNOV-TYPE SOLUTIONS FOR TRANSIENT FLOWS IN SEWERS Arturo S. León 1, Mohamed S. Ghidaoui, M. ASCE 2, Arthur R. Schmidt, M. ASCE 3, Marcelo H. García, M. ASCE 4 Abstract This work is part of a long term
More information1.060 Engineering Mechanics II Spring Problem Set 4
1.060 Engineering Mechanics II Spring 2006 Due on Monday, March 20th Problem Set 4 Important note: Please start a new sheet of paper for each problem in the problem set. Write the names of the group members
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 informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More informationPresented by: Civil Engineering Academy
Presented by: Civil Engineering Academy Open-Channel Flow Uniform Flow (See CERM Ch. 19) Characterized by constant depth volume, and cross section. It can be steady or unsteady Non-uniform Flow *Not on
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 informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Mathematics, Vol.4, No.3, 6, 39 5. ANTI-DIFFUSIVE FINITE DIFFERENCE WENO METHODS FOR SHALLOW WATER WITH TRANSPORT OF POLLUTANT ) Zhengfu Xu (Department of Mathematics, Pennsylvania
More informationSelf-Excited Vibration in Hydraulic Ball Check Valve
Self-Excited Vibration in Hydraulic Ball Check Valve L. Grinis, V. Haslavsky, U. Tzadka Abstract This paper describes an experimental, theoretical model and numerical study of concentrated vortex flow
More informationPRESSURE AND VELOCITY AMPLITUDES OF THE INCOMPRESSIBLE FLUID IN CONCENTRIC ANNULAR PASSAGE WITH OSCILLATORY BOUNDARY: TURBULENT FLOW
Journal of Engineering Science and Technology Vol. 9, No. 2 (2014) 220-232 School of Engineering, Taylor s University PRESSURE AND VELOCITY AMPLITUDES OF THE INCOMPRESSIBLE FLUID IN CONCENTRIC ANNULAR
More informationCapSel Roe Roe solver.
CapSel Roe - 01 Roe solver keppens@rijnh.nl modern high resolution, shock-capturing schemes for Euler capitalize on known solution of the Riemann problem originally developed by Godunov always use conservative
More informationA note on the carbuncle in shallow water simulations
A note on the carbuncle in shallow water simulations Friedemann Kemm Institute for Applied Mathematics and Scientific Computing, Brandenburgische Technische Universität Cottbus, Platz der Deutschen Einheit
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 informationA Central Rankine Hugoniot Solver for Hyperbolic Conservation Laws
A Central Rankine Hugoniot Solver for Hyperbolic Conservation Laws S. Jaisankar, S. V. Raghurama Rao CFD Centre, Department of Aerospace Engineering, Indian Institute of Science, Bangalore, India-56 Abstract
More informationAnalysis of the Cooling Design in Electrical Transformer
Analysis of the Cooling Design in Electrical Transformer Joel de Almeida Mendes E-mail: joeldealmeidamendes@hotmail.com Abstract This work presents the application of a CFD code Fluent to simulate the
More informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
More informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
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 informationWell-balanced shock-capturing hybrid finite volume-finite difference schemes for Boussinesq-type models
NUMAN 2010 Well-balanced shock-capturing hybrid finite volume-finite difference schemes for Boussinesq-type models Maria Kazolea 1 Argiris I. Delis 2 1 Environmental Engineering Department, TUC, Greece
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 informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More informationWaves in a Shock Tube
Waves in a Shock Tube Ivan Christov c February 5, 005 Abstract. This paper discusses linear-wave solutions and simple-wave solutions to the Navier Stokes equations for an inviscid and compressible fluid
More informationA Non-Intrusive Polynomial Chaos Method For Uncertainty Propagation in CFD Simulations
An Extended Abstract submitted for the 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada January 26 Preferred Session Topic: Uncertainty quantification and stochastic methods for CFD A Non-Intrusive
More informationOpen Channel Flow Part 2. Ch 10 Young, notes, handouts
Open Channel Flow Part 2 Ch 10 Young, notes, handouts Uniform Channel Flow Many situations have a good approximation d(v,y,q)/dx=0 Uniform flow Look at extended Bernoulli equation Friction slope exactly
More informationColloquium FLUID DYNAMICS 2012 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 2012 p.
Colloquium FLUID DYNAMICS 212 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 212 p. ON A COMPARISON OF NUMERICAL SIMULATIONS OF ATMOSPHERIC FLOW OVER COMPLEX TERRAIN T. Bodnár, L. Beneš
More informationSIMULATION OF PRECESSION IN AXISYMMETRIC SUDDEN EXPANSION FLOWS
Second International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 6-8 December 1999 SIMULATION OF PRECESSION IN AXISYMMETRIC SUDDEN EXPANSION FLOWS Baoyu GUO, Tim
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 informationSignature: (Note that unsigned exams will be given a score of zero.)
Neatly print your name: Signature: (Note that unsigned exams will be given a score of zero.) Circle your lecture section (-1 point if not circled, or circled incorrectly): Prof. Dabiri Prof. Wassgren Prof.
More informationDEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING Urban Drainage: Hydraulics. Solutions to problem sheet 2: Flows in open channels
DEPRTMENT OF CIVIL ND ENVIRONMENTL ENGINEERING Urban Drainage: Hydraulics Solutions to problem sheet 2: Flows in open channels 1. rectangular channel of 1 m width carries water at a rate 0.1 m 3 /s. Plot
More informationA wave propagation method for compressible multicomponent problems on quadrilateral grids
A wave propagation method for compressible multicomponent problems on quadrilateral grids K.-M. Shyue Department of Mathematics, National Taiwan University, Taipei 16, Taiwan Abstract We describe a simple
More informationIntroduction to the Finite Volumes Method. Application to the Shallow Water Equations. Jaime Miguel Fe Marqués
Introduction to the Finite Volumes Method. Application to the Shallow Water Equations. Jaime Miguel Fe Marqués Contents Preliminary considerations 3. Study of the movement of a fluid................ 3.
More informationLECTURE 9: Open channel flow: Uniform flow, best hydraulic sections, energy principles, Froude number
LECTURE 9: Open channel flow: Uniform flow, best hydraulic sections, energy principles, Froude number Assist. Prof. Neslihan SEMERCİ Marmara University Department of Environmental Engineering Open channel
More informationNumerical Simulation of Rapidly Varied Water Flow in the Wild River Type Water Slide
Archives of Hydro-Engineering and Environmental Mechanics Vol. 50 (2003), No. 1, pp. 3 23 Numerical Simulation of Rapidly Varied Water Flow in the Wild River Type Water Slide Kazimierz Burzyński, Michał
More informationFE Fluids Review March 23, 2012 Steve Burian (Civil & Environmental Engineering)
Topic: Fluid Properties 1. If 6 m 3 of oil weighs 47 kn, calculate its specific weight, density, and specific gravity. 2. 10.0 L of an incompressible liquid exert a force of 20 N at the earth s surface.
More informationNumerical Simulations of Turbulent Flow in Volcanic Eruption Clouds
Numerical Simulations of Turbulent Flow in Volcanic Eruption Clouds Project Representative Takehiro Koyaguchi Authors Yujiro Suzuki Takehiro Koyaguchi Earthquake Research Institute, University of Tokyo
More informationPage 1. Neatly print your name: Signature: (Note that unsigned exams will be given a score of zero.)
Page 1 Neatly print your name: Signature: (Note that unsigned exams will be given a score of zero.) Circle your lecture section (-1 point if not circled, or circled incorrectly): Prof. Vlachos Prof. Ardekani
More informationPredictionof discharge coefficient of Venturimeter at low Reynolds numbers by analytical and CFD Method
International Journal of Engineering and Technical Research (IJETR) ISSN: 2321-0869, Volume-3, Issue-5, May 2015 Predictionof discharge coefficient of Venturimeter at low Reynolds numbers by analytical
More information