^Laboratorio Nacional de Computaqao Cientifica, Rua Lauro Muller 455, Rio de Janeiro, RJ , Brasil
|
|
- Gloria Wood
- 6 years ago
- Views:
Transcription
1 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 68506, Rio de Janeiro, RJ ^Laboratorio Nacional de Computaqao Cientifica, Rua Lauro Muller 455, Rio de Janeiro, RJ , Brasil Abstract This paper presents a space-time finite element formulation for problems governed by the shallow water equations. A constant time-discontinuous approximation is adopted, while linear three node triangles are used for the spatial discretization. The streamline upwind Petrov-Galerkin (SUPG) method is applied in its equivalent variational form to fit the time discretization. Also, the correspondent semi-discrete SUPG version is established, and some numerical results are presented in order to compare the performance of these methods. 1 Introduction As it is well known, the use of the classical Galerkin method to approximate convection-dominated phenomena leads to numerical solutions contaminated by spurious oscillations that are spread over the entire computational domain. In the context of weighted residual methods, a remarkable improvement on the numerical solution of such problems was provided by the consistent variational SUPG method proposed in by Brooks and Hughesfl]. Since then, many SUPG based methods have been used for multi variable systems of equations.. In Sharkib[2] and Almeida and Galeao[3], space-time Petrov-Galerkin (STPG) formulations were derived for the compressible Euler and Navier-Stokes equations, showing their inherent control of derivatives along the characteristics. This fact turns out to be an improvement over the classical semidiscrete SUPG formulation, where only streamline derivatives are controlled. The discontinuity capturing approach proposed in those references lead to
2 404 Computer Techniques in Environmental Studies stable and accurate methods to solve all details of sharp layers and/or shock discontinuities. If we realize that the mathematical structure of the above mentioned equations is identical to those that govern the shallow water problem, it can be immediately concluded that for such problems, P-G weighted residual methods will also perform well. This was done in Bova and Carey[4] and Saleri[5], where semi-discrete P-G (SDPG) finite elements were employed, and in Carbonel, Galeao and Loula[6], where a space-time P-G formulation was derived, using linear interpolation for both spatial and temporal discretization. In this paper, the continuous linear interpolation is retained for the spatial discretization, but constant time-discontinuous interpolation is adopted. For this choice of interpolation, the weighting function gives no contribution to the terms involving time derivatives. Even so, the resulting discrete space-time Galerkin (STG) equations coincide with the semi-discrete Galerkin approximation with an Euler backward difference scheme to approximate the first time derivatives. Nevertheless, this will not occur with the P-G weighting residual methods. The numerical examples that will be presented later will show that, under these circumstances, the STPG formulation will give more accurate results than the SDPG model. 2 Problem statement Let (x, y) G Q e W* define a set of points on an horizontal plane and let z e [- A, 77] denote the vertical direction, where h(x> y) represents the water depth and rj(x,y,z) is the water surface elevation, both measured from the undisturbed water surface. We start from the 3-D incompressible Navier-Stokes equations, after turbulent time-averaging, integrating these equations along the z direction using depth-averaged horizontal velocities. Under the simplifying assumption of a hydrostatic pressure distribution (negligible vertical acceleration), we arrive at the shallow water equations: 1 * - l KA + ^(qw vt' + ^g u fo+ //(*/,*,+f^) (la) -^Ww In these equations, H = A» 77 is the total water depth and u is the depthaveraged velocity, with components w and v in x and y directions respectively. The gravitational acceleration is given by g and/is the Coriolis parameter. The wind velocity is w, with components v/ and w/ a and C are, respectively, the surface and Chezy friction coefficients, and // is the eddy viscosity.
3 Computer Techniques in Environmental Studies 405 Multiplying the third equation by g and observing that, - cclc} O\ - K^;>, W where c = (gh), and considering similar expressions for (g//),* and (gff),y, we obtain the shallow water equations in the velocity-celerity variables (see Saleri [5] ) which, in matrix form, can be written as: U,,+A.VU = F ; A.VU = A*VU (3 a) (3b) u 0 c 0 w 0 c 0, /i c o" c (3c) C (3d) Once an initial state Ua(x, y) is specified at / = 0 and appropriate boundary conditions are prescribed, the system of equations above can be solved to give the unknown column vector U. To obtain the space-time description of (3a-d) we introduce the variable change: s = (I)/, where (1) has units of velocity. Then if we define: where I is the (3x3) identity matrix, we will say that the space-time solution of the original problem (3a-d) is the (3x1) column vector U that satisfies the transformed equation: (4a) in (4b)
4 406 Computer Techniques in Environmental Studies 3 Petrov-Galerkin finite element model In order to construct the space-time finite element subspace, let us consider partitions 0 = ^ < t\ <...( < * +, of 5R* and denote by / = ((, ; +,) the n* time interval. For each n the space-time integration domain is the "slab" S» = Q x /, with boundary F = DC/,,. If we define S* as the e'** element in, e = 1,... n, where (TVy is the total number of elements in 5%, then for n = 0, 1,2,... (/) the space-timefiniteelement partition H^ ^ is such that: ^,=^^';,%=0,,x/,,; Q=UQ,;Q,nQ.=0 g=l e=l ybrz#y (5a) (ii) the space-time finite element subspace consists of continuous piecewise polynomials on the slab $, and may be discontinuous in time across the time levels t^ that is: tf* EC"(S.); U" ^ep*(5;); f/" -=O (5b) where P^ is the set of polynomials of degree less than or equal to k. According to the above definitions, the variational STPG formulation for the problem (4a-b) reads: Find t/* e»* such that for n = 0, 1,2,... W. X J(T ^. V L/" ). A*d&dt + (6a) where, = A.VU* - F" = Uj+AW*).VU'' - F* (6c)
5 Remarks Computer Techniques in Environmental Studies 407 (1) If in (6a) the integrals are taken over Q and fig instead of & and 5^, respectively; the P-G weighting function (rz.vf)'') is replaced by (ra.vu**)', and in (5b) the finite element subspace is defined for all / (making the jumping term in (6a) to disappear); and finally, if we approximate the partial time-derivative by a time-differencing operator, we generate the semi-discrete SUPG method. (2) If in (6a) we do not consider the added SUPG contribution represented by the second term under the summation symbol, we reproduce the timediscontinuous Galerkin method. If in this case, for instance, we use constant time-discontinuous interpolation, the resulting space-time finite element method will be identical to the backward Euler semi-discretefiniteelement procedure. Since constant time interpolation is used, C/,J' = 0, and therefore it is clear that the jumping condition, represented by the third integral in (6a), is the term responsible for this equivalence. Nevertheless, the STPG method and the correspondent SDPG method will be different. (3) The definition of the matrix T will be also different in these two formulations. This point is focused in the next section. 4 Purely hyperbolic problems To simplify our analysis let us assume F = 0. Using the T matrix definition found in Sharkib[2], we have: for the space-time formulation, _ s ^ -1-1/2 tfvaf 1 (7a) where x<, = t : *i = x ; X2 = y ; & (k = 0,1,2) are the local coordinates of the parent element S^ ; and A\, AI are the Jacobian flux matrices defined in (3c). for the semi-discrete formulation, - -1/2 (7b)
6 408 Computer Techniques in Environmental Studies where / (/ = 1,2) are the local coordinates of the parent element i% If, once again, for the sake of simplicity, we consider unidimensional problems and 11*'^ partitions of equal (hexat) S* elements, (7a-b) simplify to: (8a) (8b) where, 4 4 (9a) (9b) are the eigenvalues and correspondent eigenvetors of (z With these definitions we have, )"^ and /4,^ respectively. *. 1 7^-^ C «-c u-c + - u + c H-Cl-lU U 4- C \ + \ U - C c (lob) Now let us introduce the non-dimensional factor a = -~-, or what is the same, CFL = =>a>\> where CFL denotes the well known Courant-Friedrichsa Levy number. Using this factor, the definition of T*t becomes:
7 Computer Techniques in Environmental Studies 409 A 28 (lla) (lib) =av (lie)»remarks (i) Notice the intrinsic dependence of r,, on the used time step A/, which is not accounted for in %</. In the limit, as At -» 0 ; T,, -» 0, in a consistent way. This does not occur with the semi-discrete formulation because r^ is independent on At. (ii) Although not realizable from the practical point of view, %,= r^ if and only if a = 0. Remind that a must be greater than one, in order to the CFL condition be attained. We restate this comment saying that r,, approaches %</ as the time step At becomes larger, or, in other words, when accuracy decreases. (in) If we assume that u (( c, then r,, and r«/ become almost diagonal matrices, and can be replaced by, I * ^ I i; r,u=- i Because a > 1, the above ratio between brackets is always less than one. Even in the most unfavorable situation, CFL = 1 = a, r>, introduces less dissipation than r,/. In order to get a deeper insight about the performance of the STPG and SDPG methods, some numerical experiments will be performed in the next section. For these examples, constant discontinuous time interpolation and piecewise linear continuous spatial interpolation will be adopted for the STPG method. For the correspondent SDPG formulation this same spatial discretization will be used, combined, with an implicit backward finite difference scheme for time discretization.
8 410 Computer Techniques in Environmental Studies 5 Numerical results Our first example is the well known dam break problem, which consists of a wall separating two undisturbed water levels that is suddenly removed (Figure la). Friction effects are neglected and the spatial discretization is given by a 2x100 triangular elements mesh, as illustrated in Figure Ib. Figures 2-3 show the results for / = 2.50, respectively, comparing the solutions obtained with the Galerkin, the space-time and the semi-discrete formulations. For a time step Af = 0.10 (Figure 2), Galerkin solutions exhibit oscillations in the entire computational domain. This does not occur for both, the semi-discrete and the space-time P-G solutions, which accurately approximate the high gradients between the three horizontal water levels. For this time step, the semi-discrete solution is sharper than that obtained with the space-time formulation. The effect of reducing the time-step is shown in Figure 3, where the results corresponding to Af = 0.05 are plotted for / = For the Galerkin solutions, the oscillations grow up. For the semi-discrete solutions some localized oscillations appear near the sharpest layer, while a sharper solution without oscillations is obtained with the space-time method. g = 10 L= 100 (a) h= 1 L/2 4 L/2 L/ 100 I l^~~ JM L/ 100 Figure 1: Dam break problem.
9 Computer Techniques in Environmental Studies = 2.50 Galerkin Space-time Semi-discrete Figure 2: Solution for time / = 2.50, A/= = 2.50 Galerkin Space-time Semi-discrete Figure 3: Solution for time / = 2.50, A/ = 0.05.
10 412 Computer Techniques in Environmental Studies The second example, illustrated by Figure 4a, is the problem of a reflecting wave in a frictionless horizontal channel of length L = 5000, discretized with 2x10 elements, as shown in Figure 4b. The channel is open at the inflow boundary and closed at the opposed boundary. The system is subjected to a boundary condition at point A, raising the water level suddenly from the initial state of rest (H = 10) to H = 10.1, within one time step. The results can be seen in Figures 5-6. In these figures, the time-history responses for the water surface elevation at point B are depicted. For A/ = 10 (Figure 5), the curves of both, Galerkin and semi-discrete solutions, almost coincide. The STPG method presents an overdiffusive behavior for this time step, leading to a solution that progressively damps along time. However, a completely different behavior occurs when A/ = 1. For this time step, the STPG solution reaches the rectangular pulse form, while the Galerkin and SDPG solutions present some oscillations, as can be observed in Figure 6. (a) h= 10 L = 5000 (b) ~>-J zz^l JZz^J L^: Lc=rc2 Lf L/25 Figure 4: Reflecting wave in a frictionless channel.
11 Computer Techniques in Environmental Studies 413 At = 10 Galerkin -Space-time Semi-discrete J Figure 5: Solution at point B, with A/ = 10. Galerkin Space-time Semi-discrete o.oo Figure 6: Solution at point B, with A/ = I.
12 414 Computer Techniques in Environmental Studies 6 Conclusions In this paper, a STPGfiniteelement model was derived for problems governed by the shallow water equations. A piecewise linear continuous interpolation (in the space variables) was used, and piecewise constant discontinuous functions were adopted for time interpolation. In addition, the correspondent semidiscrete P-G model was also presented. The numerical results showed in this work indicate that the STPG method performs better than the SDPG formulation. Finally we should comment about the necessity of using a consistent time-space shock-capturing PG weighting function such that, in the limit, when Af-> 0, the approximated numerical characteristics approaches the true characteristics. This can not be actually provided by the SUPG method. The CAU generalized method proposed in Almeida and Galeao [3] for compressible flows, fulfill these requirements and, can therefore be successfully applied to shallow water problems. References 1. Brooks, A. N., Hughes, T. J., Streamline Upwind Petrov-Galerkin Formulation for Convection-Dominated Flows with Particular Emphasis on the Incompressible Navier-Stokes Equations, Compiit. Meth. Appl. Mech. Engrg, Vol. 32, pp , Shakib, F, Finite Element Analysis of the Compressible Euler and Navier- Stokes Equations, Ph.D. Thesis, Stanford University, Almeida, R. C, Galeao, A. C, An Adaptive Petrov-Galerkin Formulation for the Compressible Euler and Navier-Stokes Equations. Comput. Meth. Appl Mech. Engrg, Vol. 129, pp , Bova, S. W., Carey, G. F., An entropy Variable Formulation and Petrov- Galerkin Methods for the Shallow Water Equations, in: Finite Element Modeling of Environmental Problems-Surface and Subsurface Flow and Transport, ed. G Carey, John Wiley, London, England, Saleri, F., Some Stabilization Techniques in Computational Fluid Dynamics, Proceedings of the 9"' International Conference on Finite Elements in Fluids, Venezia, Carbonel, C, Galeao, A. C, Loula, A. D., A Two-dimensional Finite Element Model for Shallow Water Waves, Proceedings of the 73"' Brazilian Congress and 2"** Iberian American Congress of Mechanical Engineering, Belo Horizonte, 1995.
FINITE ELEMENT SUPG PARAMETERS COMPUTED FROM LOCAL DOF-MATRICES FOR COMPRESSIBLE FLOWS
FINITE ELEMENT SUPG PARAMETERS COMPUTED FROM LOCAL DOF-MATRICES FOR COMPRESSIBLE FLOWS Lucia Catabriga Department of Computer Science, Federal University of Espírito Santo (UFES) Av. Fernando Ferrari,
More informationOn discontinuity capturing methods for convection diffusion equations
On discontinuity capturing methods for convection diffusion equations Volker John 1 and Petr Knobloch 2 1 Universität des Saarlandes, Fachbereich 6.1 Mathematik, Postfach 15 11 50, 66041 Saarbrücken, Germany,
More informationIn Proc. of the V European Conf. on Computational Fluid Dynamics (ECFD), Preprint
V European Conference on Computational Fluid Dynamics ECCOMAS CFD 2010 J. C. F. Pereira and A. Sequeira (Eds) Lisbon, Portugal, 14 17 June 2010 THE HIGH ORDER FINITE ELEMENT METHOD FOR STEADY CONVECTION-DIFFUSION-REACTION
More informationNUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE
Proceedings of the Czech Japanese Seminar in Applied Mathematics 2005 Kuju Training Center, Oita, Japan, September 15-18, 2005 pp. 69 76 NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 219 On finite element methods for 3D time dependent convection diffusion reaction equations
More informationOn finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion
On finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion Volker John and Ellen Schmeyer FR 6.1 Mathematik, Universität des Saarlandes, Postfach 15 11
More informationStabilized Formulations and Smagorinsky Turbulence Model for Incompressible Flows
WCCM V Fifth World Congress on Computational Mechanics July 7 2, 22, Vienna, Austria Eds.: H.A. Mang, F.G. Rammerstorfer, J. Eberhardsteiner Stabilized Formulations and Smagorinsky Turbulence Model for
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 informationADAPTIVE DETERMINATION OF THE FINITE ELEMENT STABILIZATION PARAMETERS
ECCOMAS Computational Fluid Dynamics Conference 2001 Swansea, Wales, UK, 4-7 September 2001 c ECCOMAS ADAPTIVE DETERMINATION OF THE FINITE ELEMENT STABILIZATION PARAMETERS Tayfun E. Tezduyar Team for Advanced
More informationJ. Liou Tulsa Research Center Amoco Production Company Tulsa, OK 74102, USA. Received 23 August 1990 Revised manuscript received 24 October 1990
Computer Methods in Applied Mechanics and Engineering, 94 (1992) 339 351 1 A NEW STRATEGY FOR FINITE ELEMENT COMPUTATIONS INVOLVING MOVING BOUNDARIES AND INTERFACES THE DEFORMING-SPATIAL-DOMAIN/SPACE-TIME
More informationRening the submesh strategy in the two-level nite element method: application to the advection diusion equation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2002; 39:161 187 (DOI: 10.1002/d.219) Rening the submesh strategy in the two-level nite element method: application to
More informationIX CONGRESSO BRASILEIRO DE ENGENHARIA E CIÊNCIAS TÉRMICAS. 9th BRAZILIAN CONGRESS OF THERMAL ENGINEERING AND SCIENCES
IX CONGRESSO BRASILEIRO DE ENGENHARIA E CIÊNCIAS TÉRMICAS 9th BRAZILIAN CONGRESS OF THERMAL ENGINEERING AND SCIENCES Paper CIT02-0184 A NUMERICAL METHODOLOGY TO SOLVE THERMAL POLLUTION PRO- BLEMS R. S.
More informationSUPG STABILIZATION PARAMETERS CALCULATED FROM THE QUADRATURE-POINT COMPONENTS OF THE ELEMENT-LEVEL MATRICES
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 24 P. Neittaanmäki, T. Rossi, K. Majava, and O. Pironneau (eds.) W. Rodi and P. Le Quéré (assoc. eds.) Jyväskylä,
More informationSingularly Perturbed Partial Differential Equations
WDS'9 Proceedings of Contributed Papers, Part I, 54 59, 29. ISN 978-8-7378--9 MTFYZPRESS Singularly Perturbed Partial Differential Equations J. Lamač Charles University, Faculty of Mathematics and Physics,
More informationElement diameter free stability parameters. for stabilized methods applied to uids
Element diameter free stability parameters for stabilized methods applied to uids by Leopoldo P. Franca Laboratorio Nacional de Computac~ao Cientica (LNCC/CNPq) Rua Lauro Muller 455 22290 Rio de Janeiro,
More informationJURJEN DUINTJER TEBBENS
Proceedings of ALGORITMY 005 pp. 40 49 GMRES ACCELERATION ANALYSIS FOR A CONVECTION DIFFUSION MODEL PROBLEM JURJEN DUINTJER TEBBENS Abstract. When we apply the GMRES method [4] to linear systems arising
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 informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationA Review of Petrov-Galerkin Stabilization Approaches and an Extension to Meshfree Methods
ScientifiComputing A Review of Petrov-Galerkin Stabilization Approaches and an Extension to Meshfree Methods Thomas-Peter Fries, Hermann G. Matthies Institute of Scientific Computing Technical University
More information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More information3. FORMS OF GOVERNING EQUATIONS IN CFD
3. FORMS OF GOVERNING EQUATIONS IN CFD 3.1. Governing and model equations in CFD Fluid flows are governed by the Navier-Stokes equations (N-S), which simpler, inviscid, form is the Euler equations. For
More informationFluid 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 informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
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 informationThe two-dimensional streamline upwind scheme for the convection reaction equation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2001; 35: 575 591 The two-dimensional streamline upwind scheme for the convection reaction equation Tony W. H. Sheu*,1
More informationIMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1
Computational Methods in Applied Mathematics Vol. 1, No. 1(2001) 1 8 c Institute of Mathematics IMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1 P.B. BOCHEV E-mail: bochev@uta.edu
More informationInfluence of the Heat Transfer on the Pressure Field in Radial Diffusers Flows
Purdue University Purdue e-pubs International Compressor Engineering Conference School of Mechanical Engineering 2010 Influence of the Heat Transfer on the Pressure Field in Radial Diffusers Flows Jose
More informationA unified finite element formulation for compressible and incompressible flows using augmented conservation variables
. 2% mm d BB ELSEVIER Comput. Methods Appl. Mech. Engrg. 161 (1998) 229-243 Computer methods in applied mechanics and englneerlng A unified finite element formulation for compressible and incompressible
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationEvaluating the LCD algorithm for solving linear systems of equations arising from implicit SUPG formulation of compressible
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2003; 00: 22 [Version: 2002/09/8 v2.02] Evaluating the LCD algorithm for solving linear systems of equations arising
More informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
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 informationSlip flow boundary conditions in discontinuous Galerkin discretizations of the Euler equations of gas dynamics
Slip flow boundary conditions in discontinuous Galerkin discretizations of the Euler equations of gas dynamics J.J.W. van der Vegt and H. van der Ven Nationaal Lucht- en Ruimtevaartlaboratorium National
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 informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
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 hydraulic techniques for the Saint Venant Equations in arbitrarily shaped geometry
Università degli Studi di Trento Dipartimento di Matematica Dottorato di Ricerca in Matematica Ph.D. Thesis Computational hydraulic techniques for the Saint Venant Equations in arbitrarily shaped geometry
More 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 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 informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationA Finite-Element based Navier-Stokes Solver for LES
A Finite-Element based Navier-Stokes Solver for LES W. Wienken a, J. Stiller b and U. Fladrich c. a Technische Universität Dresden, Institute of Fluid Mechanics (ISM) b Technische Universität Dresden,
More informationStabilization and shock-capturing parameters in SUPG formulation of compressible flows
Comput. Methods Appl. Mech. Engrg. 195 (006) 161 163 www.elsevier.com/locate/cma Stabilization and shock-capturing parameters in SUPG formulation of compressible flows Tayfun E. Tezduyar *, Masayoshi Senga
More informationImproved discontinuity-capturing finite element techniques for reaction effects in turbulence computation
Comput. Mech. (26) 38: 356 364 DOI.7/s466-6-45- ORIGINAL PAPER A. Corsini F. Rispoli A. Santoriello T.E. Tezduyar Improved discontinuity-capturing finite element techniques for reaction effects in turbulence
More informationABSTRACT 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
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
More informationANALYSIS OF FLOW IN A CONCENTRIC ANNULUS USING FINITE ELEMENT METHOD
Nigerian Journal of Technology (NIJOTECH) Vol 35, No 2, April 2016, pp 344 348 Copyright Faculty of Engineering, University of Nigeria, Nsukka, Print ISSN: 0331-8443, Electronic ISSN: 2467-8821 wwwnijotechcom
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 informationBack Matter Index The McGraw Hill Companies, 2004
INDEX A Absolute viscosity, 294 Active zone, 468 Adjoint, 452 Admissible functions, 132 Air, 294 ALGOR, 12 Amplitude, 389, 391 Amplitude ratio, 396 ANSYS, 12 Applications fluid mechanics, 293 326. See
More informationFINITE ELEMENT ANALYSIS OF MIXED CONVECTION HEAT TRANSFER ENHANCEMENT OF A HEATED SQUARE HOLLOW CYLINDER IN A LID-DRIVEN RECTANGULAR ENCLOSURE
Proceedings of the International Conference on Mechanical Engineering 2011 (ICME2011) 18-20 December 2011, Dhaka, Bangladesh ICME11-TH-014 FINITE ELEMENT ANALYSIS OF MIXED CONVECTION HEAT TRANSFER ENHANCEMENT
More informationSpace-time finite element computation of compressible flows involving moving boundaries and interfaces*
Computer Methods in Applied Mechanics and Engineering 107 (1993) 209-223 North-Holland CMA 409 Space-time finite element computation of compressible flows involving moving boundaries and interfaces* S.K.
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 informationA note on the numerical treatment of the k-epsilon turbulence model Λ
A note on the numerical treatment of the k-epsilon turbulence model Λ Adrián J. Lew y, Gustavo C. Buscaglia z and Pablo M. Carrica Centro Atómico Bariloche and Instituto Balseiro, 8400 Bariloche, Argentina.
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 informationRECENT DEVELOPMENTS IN COMPUTATIONAL REACTOR ANALYSIS
RECENT DEVELOPMENTS IN COMPUTATIONAL REACTOR ANALYSIS Dean Wang April 30, 2015 24.505 Nuclear Reactor Physics Outline 2 Introduction and Background Coupled T-H/Neutronics Safety Analysis Numerical schemes
More informationThe Discontinuous Galerkin Method for Hyperbolic Problems
Chapter 2 The Discontinuous Galerkin Method for Hyperbolic Problems In this chapter we shall specify the types of problems we consider, introduce most of our notation, and recall some theory on the DG
More informationAnalysis of the flow and heat transfer characteristics for MHD free convection in an enclosure with a heated obstacle
Nonlinear Analysis: Modelling and Control, 2011, Vol. 16, No. 1, 89 99 89 Analysis of the flow and heat transfer characteristics for MHD free convection in an enclosure with a heated obstacle S. Parvin,
More informationIntroduction to Aerodynamics. Dr. Guven Aerospace Engineer (P.hD)
Introduction to Aerodynamics Dr. Guven Aerospace Engineer (P.hD) Aerodynamic Forces All aerodynamic forces are generated wither through pressure distribution or a shear stress distribution on a body. The
More informationFinite element stabilization parameters computed from element matrices and vectors
Comput. Methods Appl. Mech. Engrg. 190 (2000) 411±430 www.elsevier.com/locate/cma Finite element stabilization parameters computed from element matrices and vectors Tayfun E. Tezduyar *, Yasuo Osawa Mechanical
More informationcentrifugal acceleration, whose magnitude is r cos, is zero at the poles and maximum at the equator. This distribution of the centrifugal acceleration
Lecture 10. Equations of Motion Centripetal Acceleration, Gravitation and Gravity The centripetal acceleration of a body located on the Earth's surface at a distance from the center is the force (per unit
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
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 informationFluid flows through unsaturated porous media: An alternative simulation procedure
Engineering Conferences International ECI Digital Archives 5th International Conference on Porous Media and Their Applications in Science, Engineering and Industry Refereed Proceedings Summer 6-24-2014
More informationON THE BENEFIT OF THE SUMMATION-BY-PARTS PROPERTY ON INTERIOR NODAL SETS
6th European Conference on Computational Mechanics (ECCM 6 7th European Conference on Computational Fluid Dynamics (ECFD 7 11 15 June 018, Glasgow, UK ON THE BENEFIT OF THE SUMMATION-BY-PARTS PROPERTY
More informationNUMERICAL STUDIES OF VARIATIONAL-TYPE TIME-DISCRETIZATION TECHNIQUES FOR TRANSIENT OSEEN PROBLEM
Proceedings of ALGORITMY 212 pp. 44 415 NUMERICAL STUDIES OF VARIATIONAL-TYPE TIME-DISCRETIZATION TECHNIQUES FOR TRANSIENT OSEEN PROBLEM NAVEED AHMED AND GUNAR MATTHIES Abstract. In this paper, we combine
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 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 informationREDUCING ORDER METHODS APPLIED TO RESERVOIR SIMULATION
REDUCING ORDER METHODS APPLIED TO RESERVOIR SIMULATION Lindaura Maria Steffens Dara Liandra Lanznaster lindaura.steffens@udesc.br daraliandra@gmail.com Santa Catarina State University CESFI, Av. Central,
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationOn the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations
On the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations Lutz Tobiska Institut für Analysis und Numerik Otto-von-Guericke-Universität
More informationHigh-Order Finite-Volume Methods! Phillip Colella! Computational Research Division! Lawrence Berkeley National Laboratory!
High-Order Finite-Volume Methods! Phillip Colella! Computational Research Division! Lawrence Berkeley National Laboratory! Why Higher Order?! Locally-refined grids, mapped-multiblock grids smooth except
More informationModeling, Simulating and Rendering Fluids. Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan
Modeling, Simulating and Rendering Fluids Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan Applications Mostly Hollywood Shrek Antz Terminator 3 Many others Games Engineering Animating Fluids is
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 informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationOn the solution of incompressible two-phase ow by a p-version discontinuous Galerkin method
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2006; 22:741 751 Published online 13 December 2005 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/cnm.846
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 informationEntropy 2011, 13, ; doi: /e OPEN ACCESS. Entropy Generation at Natural Convection in an Inclined Rectangular Cavity
Entropy 011, 13, 100-1033; doi:10.3390/e1305100 OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Article Entropy Generation at Natural Convection in an Inclined Rectangular Cavity Mounir
More informationPressure corrected SPH for fluid animation
Pressure corrected SPH for fluid animation Kai Bao, Hui Zhang, Lili Zheng and Enhua Wu Analyzed by Po-Ram Kim 2 March 2010 Abstract We present pressure scheme for the SPH for fluid animation In conventional
More information( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #
Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue
More informationProceedings of ICFDP8 8th ASME International Congress of Fluid Dynamics and Propulsion December 14-17, 2006, Sharm El-Sheikh, Egypt
Proceedings of ICFDP8 8th ASME International Congress of Fluid Dynamics and Propulsion December 14-17, 6, Sharm El-Sheikh, Egypt ICFDP8-EG-148 ITERATIVE SOLVERS AND INFLOW BOUNDARY CONDITIONS FOR PLANE
More informationAntony Jameson. AIAA 18 th Computational Fluid Dynamics Conference
Energy Estimates for Nonlinear Conservation Laws with Applications to Solutions of the Burgers Equation and One-Dimensional Viscous Flow in a Shock Tube by Central Difference Schemes Antony Jameson AIAA
More informationNumerical Mathematics
Alfio Quarteroni Riccardo Sacco Fausto Saleri Numerical Mathematics Second Edition With 135 Figures and 45 Tables 421 Springer Contents Part I Getting Started 1 Foundations of Matrix Analysis 3 1.1 Vector
More informationUne méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé
Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé CMCS/IACS Ecole Polytechnique Federale de Lausanne Erik.Burman@epfl.ch Méthodes Numériques
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationSeveral forms of the equations of motion
Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationA STABILIZED FINITE ELEMENT METHOD FOR INCOMPRESSIBLE VISCOUS FLOWS USING A FINITE INCREMENT CALCULUS FORMULATION
A STABILIZED FINITE ELEMENT METHOD FOR INCOMPRESSIBLE VISCOUS FLOWS USING A FINITE INCREMENT CALCULUS FORMULATION Eugenio Oñate International Centre for Numerical Methods in Engineering Universidad Politécnica
More informationFINITE ELEMENT FORMULATION FOR THE VORTICITY -STREAM FUNCTION FORM OF THE INCOMPRESSIBLE EULER EQUATIONS ON MULTIPLY-CONNECTED DOMAINS* T.E.
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING 73 (1989) 331-339 NORTH-HOLLAND FINITE ELEMENT FORMULATION FOR THE VORTICITY -STREAM FUNCTION FORM OF THE INCOMPRESSIBLE EULER EQUATIONS ON MULTIPLY-CONNECTED
More informationLevel Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows
Level Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows Abstract Maged Ismail Claremont Graduate University Level Set and Phase Field methods are well-known interface-capturing
More informationContinuous adjoint based error estimation and r-refinement for the active-flux method
Continuous adjoint based error estimation and r-refinement for the active-flux method Kaihua Ding, Krzysztof J. Fidkowski and Philip L. Roe Department of Aerospace Engineering, University of Michigan,
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 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 informationReduced Order Modeling of Steady and Unsteady Flow over a Sphere
ILASS-Americas 22nd Annual Conference on Liquid Atomization and Spray Systems, Cincinnati, OH, May 21 Reduced Order Modeling of Steady and Unsteady Flow over a Sphere B. T. Helenbrook and D. R. Witman
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 informationFinite volumes for complex applications In this paper, we study finite-volume methods for balance laws. In particular, we focus on Godunov-type centra
Semi-discrete central schemes for balance laws. Application to the Broadwell model. Alexander Kurganov * *Department of Mathematics, Tulane University, 683 St. Charles Ave., New Orleans, LA 708, USA kurganov@math.tulane.edu
More informationIncompressible Flow Solver by Means of Pseudo-Compressibility Method
Incompressible Flow Solver by Means of Pseudo-Compressibility Method Árpád Veress Department of Aircraft and Ships Budapest University of Technology and Economics H Budapest, Sztoczek u. 6. J. ép., Hungary
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 informationSome notes about PDEs. -Bill Green Nov. 2015
Some notes about PDEs -Bill Green Nov. 2015 Partial differential equations (PDEs) are all BVPs, with the same issues about specifying boundary conditions etc. Because they are multi-dimensional, they can
More informationBasic Features of the Fluid Dynamics Simulation Software FrontFlow/Blue
11 Basic Features of the Fluid Dynamics Simulation Software FrontFlow/Blue Yang GUO*, Chisachi KATO** and Yoshinobu YAMADE*** 1 FrontFlow/Blue 1) is a general-purpose finite element program that calculates
More informationScalable Non-Linear Compact Schemes
Scalable Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory International Conference on Spectral and High Order Methods
More informationA finite-volume algorithm for all speed flows
A finite-volume algorithm for all speed flows F. Moukalled and M. Darwish American University of Beirut, Faculty of Engineering & Architecture, Mechanical Engineering Department, P.O.Box 11-0236, Beirut,
More informationExcerpt from the Proceedings of the COMSOL Users Conference 2006 Boston
Using Comsol Multiphysics to Model Viscoelastic Fluid Flow Bruce A. Finlayson, Professor Emeritus Department of Chemical Engineering University of Washington, Seattle, WA 98195-1750 finlayson@cheme.washington.edu
More information