A generalized method for advective-diffusive computations in engineering
|
|
- Bennett McDaniel
- 5 years ago
- Views:
Transcription
1 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 Mathematical Methods, University of A Coruña, Spain Abstract Besides the computational cost of solving advective-diffusive problems in convection dominated situations, the standard statement for this phenomenon leads to the result that mass can propagate at an infinite speed. This paradoxical result occurs as a consequence of using Fick s law [1] and it is related to the appearance of boundary layers on outflow borders when convection dominates diffusion. For these reasons we propose to use Cattaneo s law [2] instead of Fick s law as the constitutive equation of the problem. This approach leads to a totally hyperbolic system of partial differential equations. A finite diffusive velocity can be defined by using this approach. Several problems have been solved to show that the proposed formulation leads to stable results in convection dominated situations. Keywords: advection-diffusion, Taylor Galerkin, Cattaneo s equation. 1 Introduction Transport problems involving convective and diffusive processes in fluid media have a great applicability in engineering. This kind of phenomena can be modeled by using the so-called advective-diffusive equation. Unfortunately, it is very complicated to obtain an accurate and stable solution for this equation when the convective term becomes important. To overcome this problem, many interesting numerical methods have been introduced, but completely satisfactory results have not been achieved (for a detailed presentation of most of these methods see reference [3]). In this paper we review the formulation of the advective-diffusive equations applied to the spillage of a pollutant into a fluid medium. In particular, we notice that Fick s law leads to the result that mass can propagate at an infinite speed. This fact is related to the appearance of spurious oscillations in the numerical
2 564 Fluid Structure Interaction and Moving Boundary Problems solution of this equation [4, 5]. To overcome the infinite speed paradox we use a generalized constitutive equation proposed by Cattaneo in [2]. Cattaneo s law has been used in practical applications for pure-diffusive problems instead of Fick s law, but, up to the authors knowledge, it has not been used for problems with convective term. In this study the advective-diffusive problem has been formulated by using Cattaneo s law as the constitutive equation. This approach leads to a totally hyperbolic system of conservation laws. Therefore, a finite diffusive velocity can be defined. Finally, several practical examples have been solved. 2 Standard formulation of the advective-diffusive problem 2.1 Problem statement The classic description for the advective-diffusive problem under the assumption of incompressibility is based on the equations t + a x (u)+ x (q) =0 (1) q = K x (u) (2) where Eqn (1) is the pollutant mass conservation law and Eqn (2) is the constitutive equation known as Fick s law. In the above, u is the pollutant concentration, a is the velocity vector which satisfies the hydrodynamic equations of an incompressible fluid, q is the diffusive flux per unit fluid density and K is the diffusivity tensor which is assumed to be positive definite. Clearly, the system of Eqns (1) (2) is fully decoupled as we can introduce Eqn (2) into Eqn (1) and solve t + a x (u) x (K x (u)) = 0 (3) which is a scalar equation. It is well known that the Eqn (3) is a parabolic one. Therefore, boundary conditions must be imposed everywhere on the boundary of the domain. 2.2 A pure-diffusive example to analyze the infinite speed paradox In what follows we will show that the above formulation leads to mass propagation at an infinite speed. Let us consider an (incompressible) homogeneous, isotropic (hence, if I is the identity tensor, K = ki for a certain k > 0) and one-dimensional medium. We consider a pure-diffusive situation. We suppose the domain to be long enough to be approximated as infinitely long. Finally, we assume that the pollutant is added to the medium as a rapid pulse. In this case we should solve the following problem: find u: R [0, ) R such that t k 2 u =0 x2 x R t>0 (4) u(x, 0) = δ(x) x R (5)
3 Fluid Structure Interaction and Moving Boundary Problems 565 where δ is the Dirac distribution. The solution of Eqns (4) (5) is u(x, t) = 1 4πkt e x2 4kt, x R, t > 0. (6) which is the Gauss distribution function at each time t. So, it is straightforward (see [5] for further details) that the mean velocity of the pollutant is not bounded. 3 Proposed formulation of the advective-diffusive problem 3.1 Problem statement We will derive this formulation by substituting Eqn (2) by Cattaneo s law which involves a tensorial function τ. This mapping transforms each point (x,t) of the domain into that point relaxation tensor. The coordinates of the relaxation tensor are specific diffusion process times. Up to the authors knowledge, Cattaneo s equation has been only used in non-advective thermal problems. Thus, we had to find Cattaneo s equation with convective term [5]. Hence, basic equations for the transport problem described by using Cattaneo s law are t + a x (u)+ x (q) =0 (7) ( ) q q + τ t + x (q) a = K x (u) (8) where Eqn (8) is Cattaneo s law with convective term. It should be noted that we are using a generalized constitutive equation as we recover the classic formulation when τ is the zero tensor. 3.2 A pure-diffusive example with a finite propagation velocity In order to compare the solution of the classic formulation with the solution of the generalized formulation we now solve the Cattaneo counterpart of Eqns (4) (5). In this simple case the system of Eqns (7) (8) can be reduced to a second order partial differential equation (for further details, see [5]). Now we need two initial conditions because this problem involves second order derivatives with respect to the time. Then, we consider an (incompressible) homogeneous, isotropic, onedimensional and non-convective medium. With the above assumptions we can state this problem as [5]: find u: R [0, ) R such that τ 2 u t 2 + t k 2 u =0 x2 x R t>0 (9) u(x, 0) = δ(x) x R (10) (x, 0) = 0 t x R (11)
4 566 Fluid Structure Interaction and Moving Boundary Problems If we call c = k/τ mass wave celerity and we denote s = c 2 t 2 x 2,the solution of Eqns (9) (11) is [ 1 2 e c 2 2k t δ( x ct)+ c ( ) ( )] c 2k I 0 2k s + c2 c 2ks ti 1 2k s, x ct u(x, t) = 0, x >ct (12) where I 0 and I 1 are the modified Bessel functions of the first kind of order 0 and 1. We compare in Figure 1 solutions of Eqns (4) (5) and Eqns (9) (11) at t =4and t =10respectively. Clearly, if we use Cattaneo s equation a wave front exists which advances with a celerity c. 3.3 Study of the proposed model as a system of conservation laws The system of Eqns (7) (8) can not be reduced to a second order partial differential equation in multidimensional problems with a non-constant velocity field. In this case we must solve a coupled system of first order partial differential equations. If we suppose the medium to be (incompressible) homogeneous and isotropic, the system of Eqns (7) (8) can be written as a system of conservation laws. As we consider two-dimensional problems, we will use the notation q =(q 1,q 2 ) T, a = (a 1,a 2 ) T. Then, the system of Eqns (7) (8) can be written as U t + x (F )=S (13) where u ua 1 + q 1 ua 2 + q 2 0 U = τq 1 ; F = τq 1 a 1 + ku τq 1 a 2 ; S = q 1 (14) τq 2 τq 2 a 1 τq 2 a 2 + ku q 2 It can be shown that the system of Eqns (13) is totally hyperbolic because its associated characteristic equation [6] has 3 different real solutions. At this point, it is very useful to introduce the dimensionless number H = a c (15) which plays a similar role to Mach number in gas dynamics or Froude number in shallow water problems. Then, we define supercritical flow as one characterized by a local H number such that H > 1. In the same way, we define subcritical flow as one characterized by H<1 and critical flow as the flow which satisfies H =1. It should be noted that in supercritical flow it is not possible a pollutant transport towards upstream as the circular (diffusive) mode of propagation travels with a velocity lesser than the convective one.
5 Fluid Structure Interaction and Moving Boundary Problems Cattaneo s law Fick s law Cattaneo s law Fick s law Concentration Concentration Spatial coordinate Spatial coordinate Figure 1: Comparison at t =4(left) and at t =10(right) between the solution of Eqns (4) (5) (dashed line) and the solution of Eqns (9) (11) (solid line). Parameters k and τ have a value of one. The question of boundary conditions which can be prescribed for the system of Eqns (13) is not trivial. A very frequent situation is to have to impose a boundary condition at a fixed solid wall. In this situation only the normal component of flux q must be prescribed. Then, boundary condition in this case will be q n =0, being n the unit outward normal vector to that boundary. However, computational domains are usually limited by inflow and outflow boundaries, as well. We call inflow boundary the part of the boundary in which a n < 0 and outflow boundary the part in which a n > 0. In these types of boundaries the issue of boundary conditions is more complicated. However, by means of a Riemann analysis [6] in the direction of the outward normal we conclude the following: on a supercritical inflow boundary all components of U must be prescribed while no components should be imposed on a supercritical outflow boundary. In contrast, on a subcritical outflow boundary only one component of U will be prescribed while 2 components of U must be imposed on a subcritical inflow boundary. Independently of boundary conditions, an initial condition must be imposed to Eqns (13). 4 Numerical examples This section is devoted to the presentation of the results obtained by using the second order Taylor Galerkin method proposed by Donea in [7]. We will present two practical examples which consist of a pollutant transport in a rectangular channel. The unstructured meshes were generated by using the code GEN4U based on the algorithm by Sarrate and Huerta [8]. 4.1 Subcritical flow in a rectangular channel We consider a rectangular channel over the domain [0, 10] [0, 2]. A typical mesh of quadrilaterals has been used. Grid steps size are x = y =0.25 m.
6 568 Fluid Structure Interaction and Moving Boundary Problems The values of the parameters are k = 1 and τ = 1 which implies c = 1. The velocity field employed in the calculations is a(x 1,x 2 )=(0.8x 2 (2 x 2 ), 0) T (16) which is a divergence free field. By using Eqn (15) and Eqn (16) it is straightforward that the considered flow is subcritical at each point of the domain. As initial conditions we use u(x 1,x 2 )=e (x2 1 +x2 2 ) ; q(x 1,x 2 )=0. (17) As boundary conditions, we prescribe concentration and q 2 on the inflow boundary; q 1 on the outflow boundary and q 2 = 0 at the solid wall (in all cases we impose the values which result of initial conditions). A time step of t =0.1s has been used in the calculations. We step in time until the steady state has been reached. The solution is considered to be steady when Residual = U n+1 U n U n (18) We show in Figure 2 concentration solutions at t =3s, t =6s, t =9sandt =12 s Further, we plot the steady state solution. 4.2 Transcritical flow in a rectangular channel We consider a rectangular channel over the domain [0, 10] [0, 2]. The computational mesh (1375 elements) is plotted in Figure 3. The values of the parameters are k =10 2 and τ =1. Therefore, the mass wave celerity is c =0.1. We will solve this problem by using the velocity field a(x 1,x 2 )=(2x 2 (2 x 2 ), 0) T (19) which is a solenoidal one. By using Eqn (15) we find the flow to be supercritical except near the walls of the channel. In fact, this is a high velocity problem as it satisfies H = 20 at the straight line x 2 = 1. We use Eqn (17) as initial conditions. As boundary conditions, we prescribe all unknowns on the supercritical inflow boundary (we impose the values which result of the initial conditions). At the solid wall we impose q 2 =0which is also consistent with initial conditions. A time step of t = s has been used in the calculations. We march in time until the Eqn (18) is satisfied. We show in Figure 4 the initial condition involving concentration and concentration solutions at t =1.5s, t =3s, t =4.5s, t =6s and t =12sFurther, we plot the steady state solution. 5 Conclusions In this paper we propose to use Cattaneo s law as the constitutive equation of the advective-diffusive problem. This approach leads to a wave-like solution and it
7 Fluid Structure Interaction and Moving Boundary Problems 569 Figure 2: Subcritical channel flow. Initial condition and concentration solutions at t =3s, t =6s, t =9sandt =12s in downward direction. We plot the steady state solution in the lower graphic. Figure 3: Transcritical channel flow. Computational mesh (1375 elements). avoids the infinite speed paradox which is reached by using the standard formulation. The proposed formulation constitutes a generalized approach for advectivediffusive phenomena as the standard formulation can be considered as a subcase
8 570 Fluid Structure Interaction and Moving Boundary Problems Figure 4: Transcritical channel flow. Initial condition and concentration solutions at t =1.5 s, t =3s, t =4.5 s, t =6sandt =12sindownward direction. We plot the steady state solution in the lower graphic. of the proposed one. When Cattaneo s law is used two kinds of flow can occur: subcritical flow (pollutant can propagate towards any direction) and supercritical flow (propagation towards upstream is not possible). From a numerical point of view, two problems have been solved to show that the Cattaneo-type model leads, on the whole, to stable solutions even if we use the standard Galerkin method for the spatial discretization.
9 Acknowledgements This work has been partially supported by Grant # PGDIT01PXI11802PR of the SXID of the Xunta de Galicia, Grant # DPI of the SGPICT of the Ministerio de Ciencia y Tecnología ofthe Spanish Government, research fellowships of the Universidad de A Coruña and research fellowships of the Fundación de la Ingeniería Civil de Galicia. References Fluid Structure Interaction and Moving Boundary Problems 571 [1] Fick, A., Uber diffusion, Poggendorff s Annalen der Physik und Chemie, 94, pp , [2] Cattaneo, M.C., Sur une forme de l equation de la chaleur éliminant le paradoxe d une propagation instantaneé, Comptes Rendus de L Academie des Sciences: Series I-Mathematics, 247, pp , [3] Donea, J. and Huerta, A., Finite element methods for flow problems, John Wiley & Sons: Chichester, [4] Gómez, H., Colominas, I., Navarrina, F. and Casteleiro, M., On the intrinsic instability of the advection diffusion equation, Proc. of the 4 th European Congress on Computational Methods in Applied Sciences and Engineering (CDROM), eds. P. Neittaanmäki,T.Rossi,S.Korotov,E.Oñate,J.Périaux y D. Knörzer, Jyväskylä, Finland, [5] Gómez, H., A new formulation for the advective-diffusive transport problem, Technical Report (in Spanish), University of A Coruña, [6] Courant R. and Hilbert D. Methods of mathematical physics. Vol II., John Wiley & Sons: New York, [7] Donea, J., A Taylor Galerkin method for convective transport problems, International Journal for Numerical Methods in Engineering, 20, pp , [8] Sarrate, J. and Huerta, A., Efficient unstructured quadrilateral mesh generation, International Journal for Numerical Methods in Engineering, 49, pp , 2000.
RECENT ACHIEVEMENTS ON THE USE OF PURE HYPERBOLIC CATTANEO-TYPE CONVECTION-DIFFUSION MODELS IN CFD
8th. World Congress on Computational Mechanics (WCCM8) 5th European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS 2008) June 30 July 5, 2008 Venice, Italy RECENT ACHIEVEMENTS
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 information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Partial differential equations arise in a number of physical problems, such as fluid flow, heat transfer, solid mechanics and biological processes. These
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 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 informationFlood Routing by the Non-Linear Muskingum Model: Conservation of Mass and Momentum
Archives of Hydro-Engineering and Environmental Mechanics Vol. 56 (29), No. 3 4, pp. 121 137 IBW PAN, ISSN 1231 3726 Flood Routing by the Non-Linear Muskingum Model: Conservation of Mass and Momentum Dariusz
More informationNUMERICAL 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 informationDiffusion / Parabolic Equations. PHY 688: Numerical Methods for (Astro)Physics
Diffusion / Parabolic Equations Summary of PDEs (so far...) Hyperbolic Think: advection Real, finite speed(s) at which information propagates carries changes in the solution Second-order explicit methods
More information1- Constitutive Relations of Heat Flux
1- Constitutive Relations of Heat Flux By the second law of thermodynamics, there exists a physical quantity Q that is, at a given time instant, associated with each surface in a non-isothermal body. This
More informationCharacteristics for IBVP. Notes: Notes: Periodic boundary conditions. Boundary conditions. Notes: In x t plane for the case u > 0: Solution:
AMath 574 January 3, 20 Today: Boundary conditions Multi-dimensional Wednesday and Friday: More multi-dimensional Reading: Chapters 8, 9, 20 R.J. LeVeque, University of Washington AMath 574, January 3,
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 informationThe Finite Difference Method
Chapter 5. The Finite Difference Method This chapter derives the finite difference equations that are used in the conduction analyses in the next chapter and the techniques that are used to overcome computational
More informationThe Euler Equation of Gas-Dynamics
The Euler Equation of Gas-Dynamics A. Mignone October 24, 217 In this lecture we study some properties of the Euler equations of gasdynamics, + (u) = ( ) u + u u + p = a p + u p + γp u = where, p and u
More informationScientific Computing I
Scientific Computing I Module 10: Case Study Computational Fluid Dynamics Michael Bader Winter 2012/2013 Module 10: Case Study Computational Fluid Dynamics, Winter 2012/2013 1 Fluid mechanics as a Discipline
More informationTHREE-DIMENSIONAL FINITE DIFFERENCE MODEL FOR TRANSPORT OF CONSERVATIVE POLLUTANTS
Pergamon Ocean Engng, Vol. 25, No. 6, pp. 425 442, 1998 1998 Elsevier Science Ltd. All rights reserved Printed in Great Britain 0029 8018/98 $19.00 + 0.00 PII: S0029 8018(97)00008 5 THREE-DIMENSIONAL FINITE
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
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 informationFINITE 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 informationNumerical simulation of earthing grids
Proceedings of the 6th WSEAS International Conference on Power Systems, Lisbon, Portugal, September 22-24, 2006 78 Numerical simulation of earthing grids IGNASI COLOMINAS, FERMÍN NAVARRINA & MANUEL CASTELEIRO
More informationCavitation modeling using compressible Navier Stokes and Korteweg equations
Computational Methods in Multiphase Flow VIII 425 Cavitation modeling using compressible Navier Stokes and Korteweg equations T. Ait Ali 1, S. Khelladi 1, L. Ramirez 2 & X. Nogueira 2 1 DynFluid Laboratory,
More informationA posteriori error estimates applied to flow in a channel with corners
Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,
More informationDiffusion Processes. Lectures INF2320 p. 1/72
Diffusion Processes Lectures INF2320 p. 1/72 Lectures INF2320 p. 2/72 Diffusion processes Examples of diffusion processes Heat conduction Heat moves from hot to cold places Diffusive (molecular) transport
More informationMulti-D MHD and B = 0
CapSel DivB - 01 Multi-D MHD and B = 0 keppens@rijnh.nl multi-d MHD and MHD wave anisotropies dimensionality > 1 non-trivial B = 0 constraint even if satisfied exactly t = 0: can numerically generate B
More information2. Conservation of Mass
2 Conservation of Mass The equation of mass conservation expresses a budget for the addition and removal of mass from a defined region of fluid Consider a fixed, non-deforming volume of fluid, V, called
More informationLattice Boltzmann Method for Fluid Simulations
Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas April 14, 2011 1 Introduction In the last two decades, the Lattice Boltzmann method (LBM) has emerged as a promising tool
More informationNumerical Methods for Problems with Moving Fronts Orthogonal Collocation on Finite Elements
Electronic Text Provided with the Book Numerical Methods for Problems with Moving Fronts by Bruce A. Finlayson Ravenna Park Publishing, Inc., 635 22nd Ave. N. E., Seattle, WA 985-699 26-524-3375; ravenna@halcyon.com;www.halcyon.com/ravenna
More informationGauge finite element method for incompressible flows
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 34: 701 710 Gauge finite element method for incompressible flows Weinan E a, *,1 and Jian-Guo Liu b,2 a Courant Institute
More informationIntroduction to Partial Differential Equation - I. Quick overview
Introduction to Partial Differential Equation - I. Quick overview To help explain the correspondence between a PDE and a real world phenomenon, we will use t to denote time and (x, y, z) to denote the
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 informationA Taylor Galerkin approach for modelling a spherically symmetric advective dispersive transport problem
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2008; 24:49 63 Published online 20 November 2006 in Wiley InterScience (www.interscience.wiley.com)..955 A Taylor Galerkin
More informationCONVECTIVE HEAT TRANSFER
CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 4 HEAT TRANSFER IN CHANNEL FLOW BASIC CONCEPTS BASIC CONCEPTS Laminar
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 informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationDraft Notes ME 608 Numerical Methods in Heat, Mass, and Momentum Transfer
Draft Notes ME 608 Numerical Methods in Heat, Mass, and Momentum Transfer Instructor: Jayathi Y. Murthy School of Mechanical Engineering Purdue University Spring 00 c 1998 J.Y. Murthy and S.R. Mathur.
More informationNumerical Solutions for the Lévêque Problem of Boundary Layer Mass or Heat Flux
Excerpt from the Proceedings of the COMSOL Conference 2008 Hannover Numerical Solutions for the Lévêque Problem of Boundary Layer Mass or Heat Flux Ekkehard Holzbecher Weierstrass Institute for Applied
More informationSimulation Study on the Generation and Distortion Process of the Geomagnetic Field in Earth-like Conditions
Chapter 1 Earth Science Simulation Study on the Generation and Distortion Process of the Geomagnetic Field in Earth-like Conditions Project Representative Yozo Hamano Authors Ataru Sakuraba Yusuke Oishi
More informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
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 informationCFD Analysis of Forced Convection Flow and Heat Transfer in Semi-Circular Cross-Sectioned Micro-Channel
CFD Analysis of Forced Convection Flow and Heat Transfer in Semi-Circular Cross-Sectioned Micro-Channel *1 Hüseyin Kaya, 2 Kamil Arslan 1 Bartın University, Mechanical Engineering Department, Bartın, Turkey
More informationHomework 4 in 5C1212; Part A: Incompressible Navier- Stokes, Finite Volume Methods
Homework 4 in 5C11; Part A: Incompressible Navier- Stokes, Finite Volume Methods Consider the incompressible Navier Stokes in two dimensions u x + v y = 0 u t + (u ) x + (uv) y + p x = 1 Re u + f (1) v
More informationLATTICE BOLTZMANN METHOD AND DIFFUSION IN MATERIALS WITH LARGE DIFFUSIVITY RATIOS
THERMAL SCIENCE: Year 27, Vol. 2, No. 3, pp. 73-82 73 LATTICE BOLTZMANN METHOD AND DIFFUSION IN MATERIALS WITH LARGE DIFFUSIVITY RATIOS by Edouard WALTHER *, Rachid BENNACER, and Caroline DE SA LMT, ENS
More informationAbstract. Introduction
Stability Criterion for Explicit Schemes (Finite-Difference Method) on the Solution of the Advection-Diffusion Equation L.F. León, P.M. Austria Mexican Institute of Water Technology, Paseo Cuauhnáhuac
More information1. INTRODUCTION TO CFD SPRING 2018
1. INTRODUCTION TO CFD SPRING 018 1.1 What is computational fluid dynamics? 1. Basic principles of CFD 1.3 Stages in a CFD simulation 1.4 Fluid-flow equations 1.5 The main discretisation methods Appendices
More informationOn limiting for higher order discontinuous Galerkin method for 2D Euler equations
On limiting for higher order discontinuous Galerkin method for 2D Euler equations Juan Pablo Gallego-Valencia, Christian Klingenberg, Praveen Chandrashekar October 6, 205 Abstract We present an implementation
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 information(1:1) 1. The gauge formulation of the Navier-Stokes equation We start with the incompressible Navier-Stokes equation 8 >< >: u t +(u r)u + rp = 1 Re 4
Gauge Finite Element Method for Incompressible Flows Weinan E 1 Courant Institute of Mathematical Sciences New York, NY 10012 Jian-Guo Liu 2 Temple University Philadelphia, PA 19122 Abstract: We present
More informationChapter 17. Finite Volume Method The partial differential equation
Chapter 7 Finite Volume Method. This chapter focusses on introducing finite volume method for the solution of partial differential equations. These methods have gained wide-spread acceptance in recent
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
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 informationDimensionality influence on energy, enstrophy and passive scalar transport.
Dimensionality influence on energy, enstrophy and passive scalar transport. M. Iovieno, L. Ducasse, S. Di Savino, L. Gallana, D. Tordella 1 The advection of a passive substance by a turbulent flow is important
More informationNUMERICAL SIMULATION OF HEAT TRANSFER MECHANISMS DURING FEMTOSECOND LASER HEATING OF NANO-FILMS USING 3-D DUAL PHASE LAG MODEL
Proceedings of HTFED04 004 ASME Heat Transfer/ Fluids Engineering Summer Conference Charlotte, North Carolina, USA, July 11-15, 004 HT-FED004-5683 NUMERICAL SIMULATION OF HEAT TRANSFER MECHANISMS DURING
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 informationThird-Order Active-Flux Scheme for Advection Diffusion: Hyperbolic Diffusion, Boundary Condition, and Newton Solver
Preprint accepted for publication in Computers and Fluids, 16 http://d.doi.org/1.116/j.compfluid.15.1. Third-Order Active-Flu Scheme for Advection Diffusion: Hyperbolic Diffusion, Boundary Condition, and
More informationFINITE-VOLUME SOLUTION OF DIFFUSION EQUATION AND APPLICATION TO MODEL PROBLEMS
IJRET: International Journal of Research in Engineering and Technology eissn: 39-63 pissn: 3-738 FINITE-VOLUME SOLUTION OF DIFFUSION EQUATION AND APPLICATION TO MODEL PROBLEMS Asish Mitra Reviewer: Heat
More informationSolving Pure Torsion Problem and Modelling Radionuclide Migration Using Radial Basis Functions
International Workshop on MeshFree Methods 3 1 Solving Pure Torsion Problem and Modelling Radionuclide Migration Using Radial Basis Functions Leopold Vrankar (1), Goran Turk () and Franc Runovc (3) Abstract:
More informationModule 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 1: Finite Difference Method
file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture1/1_1.htm 1 of 1 6/19/2012 4:29 PM The Lecture deals with: Classification of Partial Differential Equations Boundary and Initial Conditions Finite Differences
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 information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More 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 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 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 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 informationA Finite Element Analysis on MHD Free Convection Flow in Open Square Cavity Containing Heated Circular Cylinder
American Journal of Computational Mathematics, 2015, 5, 41-54 Published Online March 2015 in SciRes. http://www.scirp.org/journal/ajcm http://dx.doi.org/10.4236/ajcm.2015.51003 A Finite Element Analysis
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 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 informationFINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES
FINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES SHRUTI JAIN B.Tech III Year, Electronics and Communication IIT Roorkee Tutors: Professor G. Biswas Professor S. Chakraborty ACKNOWLEDGMENTS I would like
More information3.4. Monotonicity of Advection Schemes
3.4. Monotonicity of Advection Schemes 3.4.1. Concept of Monotonicity When numerical schemes are used to advect a monotonic function, e.g., a monotonically decreasing function of x, the numerical solutions
More informationSelf-similar shock wave solutions for heat conductions in solids and for non-linear Maxwell equations
Self-similar shock wave solutions for heat conductions in solids and for non-linear Maxwell equations Imre Ferenc Barna 1,2 3 & Robert Kersner 1) Wigner Research Center of the Hungarian Academy of Sciences
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 informationNon-fourier Heat Conduction in a Long Cylindrical Media with Insulated Boundaries and Arbitrary Initial Conditions
Australian Journal of Basic and Applied Sciences, 3(): 65-663, 009 ISSN 1991-8178 Non-fourier Heat Conduction in a Long Cylindrical Media with Insulated Boundaries and Arbitrary Initial Conditions Mehdi
More informationBasic Theorems in Dynamic Elasticity
Basic Theorems in Dynamic Elasticity 1. Stress-Strain relationships 2. Equation of motion 3. Uniqueness and reciprocity theorems 4. Elastodynamic Green s function 5. Representation theorems Víctor M. CRUZ-ATIENZA
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 informationEngineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical
Outline Geurts Book Department of Fluid Mechanics, Budapest University of Technology and Economics Spring 2013 Outline Outline Geurts Book 1 Geurts Book Origin This lecture is strongly based on the book:
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationFluid Thermal Interaction of High Speed Compressible Viscous Flow Past Uncooled and Cooled Structures by Adaptive Mesh
The 20th Conference of echanical Engineering Network of Thailand 18-20 October 2006, Nakhon Ratchasima, Thailand luid Thermal nteraction of High peed Compressible Viscous low Past Uncooled and Cooled tructures
More informationNumerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders
Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders A. Jugal M. Panchal, B. A M Lakdawala 2 A. M. Tech student, Mechanical Engineering Department, Institute
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 informationCRITERIA FOR THE CHOICE OF FLOOD ROUTING METHODS IN
Criteria for the choice of flood routing methods in natural... CRITERIA FOR THE CHOICE OF FLOOD ROUTING METHODS IN NATURAL CHANNELS WITH OVERBANK FLOWS Roger Moussa 1 Abstract: The classification of river
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 informationTo study the motion of a perfect gas, the conservation equations of continuity
Chapter 1 Ideal Gas Flow The Navier-Stokes equations To study the motion of a perfect gas, the conservation equations of continuity ρ + (ρ v = 0, (1.1 momentum ρ D v Dt = p+ τ +ρ f m, (1.2 and energy ρ
More informationAPPLIED PARTIM DIFFERENTIAL EQUATIONS with Fourier Series and Boundary Value Problems
APPLIED PARTIM DIFFERENTIAL EQUATIONS with Fourier Series and Boundary Value Problems Fourth Edition Richard Haberman Department of Mathematics Southern Methodist University PEARSON Prentice Hall PEARSON
More information5. Advection and Diffusion of an Instantaneous, Point Source
1 5. Advection and Diffusion of an Instantaneous, Point Source In this chapter consider the combined transport by advection and diffusion for an instantaneous point release. We neglect source and sink
More informationOPEN CHANNEL FLOW. Computer Applications. Numerical Methods and. Roland Jeppson. CRC Press UNIVERSITATSB'BUOTHEK TECHNISCHE. INFORMATlONSBiBUOTHEK
OPEN CHANNEL FLOW Numerical Methods and Computer Applications Roland Jeppson TECHNISCHE INFORMATlONSBiBUOTHEK UNIVERSITATSB'BUOTHEK HANNOVER Si. i. CRC Press Taylor &.Francis Group Boca Raton London New
More informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
More informationBlock 3 Open channel flow
Numerical Hydraulics Block 3 Open channel flow Markus Holzner Contents of the course Block 1 The equations Block Computation of pressure surges Block 3 Open channel flow (flow in rivers) Block 4 Numerical
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 informationLinear Hyperbolic Systems
Linear Hyperbolic Systems 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, 2014 1 / 56 We study some basic
More information0.2. CONSERVATION LAW FOR FLUID 9
0.2. CONSERVATION LAW FOR FLUID 9 Consider x-component of Eq. (26), we have D(ρu) + ρu( v) dv t = ρg x dv t S pi ds, (27) where ρg x is the x-component of the bodily force, and the surface integral is
More informationAccuracy of Muskingum-Cunge flood routing
Alternative Hydraulics Paper 3, 02.03.2011 Accuracy of Musingum-Cunge flood routing Institute of Hydraulic and Water Resources Engineering Vienna University of Technology Karlsplatz 13/222, 1040 Vienna,
More informationDirect numerical simulation database for supercritical carbon dioxide
Direct numerical simulation database for supercritical carbon dioxide S. Pandey 1, X. Chu 2, E. Laurien 3 Emails: sandeep.pandey@ike.uni-stuttgart.de 1 xu.chu@itlr.uni-stuttgart.de 2 laurien@ike.unistuttgart.de
More informationCFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE
GANIT J. Bangladesh Math. Soc. (ISSN 1606-3694) 31 (2011) 33-41 CFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE Sharmina Hussain Department of Mathematics and Natural Science BRAC University,
More informationMATLAB Solution of Flow and Heat Transfer through a Porous Cooling Channel and the Conjugate Heat Transfer in the Surrounding Wall
MATLAB Solution of Flow and Heat Transfer through a Porous Cooling Channel and the Conjugate Heat Transfer in the Surrounding Wall James Cherry, Mehmet Sözen Grand Valley State University, cherryj1@gmail.com,
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 informationECE539 - Advanced Theory of Semiconductors and Semiconductor Devices. Numerical Methods and Simulation / Umberto Ravaioli
ECE539 - Advanced Theory of Semiconductors and Semiconductor Devices 1 General concepts Numerical Methods and Simulation / Umberto Ravaioli Introduction to the Numerical Solution of Partial Differential
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 7
Numerical Fluid Mechanics Fall 2011 Lecture 7 REVIEW of Lecture 6 Material covered in class: Differential forms of conservation laws Material Derivative (substantial/total derivative) Conservation of Mass
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 9 (214 ) 599 64 1th International Conference on Mechanical Engineering, ICME 213 Validation criteria for DNS of turbulent heat
More informationA note on stability in three-phase-lag heat conduction
Universität Konstanz A note on stability in three-phase-lag heat conduction Ramón Quintanilla Reinhard Racke Konstanzer Schriften in Mathematik und Informatik Nr. 8, März 007 ISSN 1430-3558 Fachbereich
More information