FEM solution of the ψ-ω equations with explicit viscous diffusion 1
|
|
- Florence Dorsey
- 5 years ago
- Views:
Transcription
1 FEM solution of te ψ-ω equations wit explicit viscous diffusion J.-L. Guermond and L. Quartapelle 3 Abstract. Tis paper describes a variational formulation for solving te D time-dependent incompressible Navier Stokes equations expressed in terms of te stream function and vorticity. Te difference between te proposed approac and te standard one is tat te vorticity equation is interpreted as a dynamical equation governing te evolution of (te weak Laplacian of) te stream function wile te Poisson equation for te stream function is used as an expression to evaluate te distribution of te vorticity in te domain and on te boundary. A time discretization is adopted wit te viscous diffusion made explicit, wic leads to split te viscous effects from te incompressibility similarly to te fractional-step projection metods for te primitive variable equations. In some sense, te present metod generalizes to te variational framework a well-known idea tat is used in nite differences approximations and wic is based on a Taylor series expansion of te stream function on te boundary. Some error estimates and some numerical results are given. INTRODUCTION Te prototype nite element procedure for solving te D Stokes equations formulated in terms of te vorticity and stream function is te uncoupled solution metod for te biarmonic problem introduced by Glowinski and Pironneau []. Tis approac can compute te solution of te Stokes problem by an uncoupled direct metod. in Proc. Computational Fluid Dynamics 98, ECCOMAS96, (Septembre 998, Atne, Grce), K.D. Papailiou, D. Tsaalis, J. Priaux, M. Pandol (Eds.), Wiley, Vol. (998) Laboratoire d Informatique pour la Mécanique et les Sciences de l Ingénieur, CNRS, BP 33, 93, Orsay, France (guermond@limsi.fr). 3 Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci, 3, 33 Milano, Italy. In te case of te time-dependent equations, tis metod assumes necessarily an implicit treatment of te viscous term and determines te unknown boundary value of te vorticity by means of an operator associated wit tat part of te boundary were no-slip velocity conditions are prescribed. Wile te tecniques based on tis metod or on te related idea of an in uence matrix guarantee good stability properties, tey are not so easy to implement; ence one may be tempted to make a trade-off between stability and simplicity. In te present paper we investigate one possible alternative tecnique for solving te evolutionary Navier Stokes equations in two dimensions expressed in a weak variational form. In te context of nite differences, wic are known for teir simplicity, one classical approac is to assume an explicit treatment of te viscous term to derive vorticity boundary formulas. In tis case, te Neumann condition for te stream function can be used as te last piece of te time-stepping algoritm. More precisely, te derivative boundary condition can be interpreted as a relationsip specifying te boundary value of te new vorticity after te time advancement of te (internal distribution of) vorticity as been completed and after te new stream function as been determined, for details see e.g. Peyret and Taylor [7], E and Liu [] or Napolitano et al. [6]. Te aim of tis paper is to develop te variational counterpart of tis idea and to sow tat an interpretation alternative to te standard one can be given to te vorticity equation written in weak form wic leads to a new and very simple algoritm for te numerical solution of te unsteady ψ-ω equations by means of nite elements. VARIATIONAL PROBLEM Te proposed numerical metod stems from te following variational statement of te incompressible Navier Stokes equation for D ows in a simply connected, c 998 J.-L. Guermond and L. Quartapelle ECCOMAS 98. Publised in 998 by Jon Wiley & Sons, Ltd.
2 bounded domain Ω wit a smoot boundary Γ = Ω. For te initial solenoidal velocity u belonging to H (Ω) and te body force f L (, T ; L (Ω)), we ave te problem Find ψ L (, T ; H (Ω)) C(, T ; H (Ω)) and ω L (, T ; L (Ω)) C(, T ; H (Ω)), so tat φ H (Ω), (( ψ) t=, φ) = (u, φ ẑ), and tat, for all t >, ( ψ H ψ (Ω), t, ψ ) + ν ω, ψ + b(ω, ψ, ψ ) = (f, ψ ẑ), v L (Ω), (ω, v) + ( ψ, v) =. () Homogeneous velocity boundary conditions for te sake of simplicity. Te trilinear form b(ω, ψ, ψ ) associated wit te advection term is written in rotational form, namely, b(ω, ψ, ψ ) = (ω ψ, ψ ẑ), () to guarantee conservation of te (kinetic) energy, wen ν =, also in te spatially discrete case of H - conformal approximations. Note te unusual form of te evolutionary term wic involves mixed time and space derivatives of ψ. Actually, tis weak form of te dynamical equation is te most natural one witin te present variational setting since it stems from te momentum equation were 3 SPATIAL DISCRETIZATION Let W and Ψ, be two nite dimensional subspaces of H (Ω) and H (Ω), respectively. We assume tat Ψ, W, and W, Ψ, satisfy standard interpolation and inverse properties (described in detail in [5]). FIRST ORDER SCHEME In tis section we approximate te time derivative by means of te rst order Euler sceme. Let [, T ] be a nite time interval and N be an integer. We denote δt = T/N and t n = nδt for n N. For any function of time, ϕ(t), we denote ϕ n = ϕ(t n). Te fully discrete problem is formulated as follows. Te initialization step reads: Find ψ Ψ, so tat, φ Ψ,, (3) ( ψ, φ ) = (u, φ ẑ), Find ω W so tat, v W, (ω, v ) = ( ψ, v ). Ten for eac n, carry out te following two uncoupled steps: Find ψ n+ Ψ, so tat, φ Ψ,, ( (ψ n+ ψ), n φ ) δt + b(ω n, ψ n+, φ ) = ν ( ω n, φ ) + (f n+, φ ẑ), () (5) te velocity is replaced by ψ ẑ, te viscous term is written as te curl of ω, and te velocity test functions belonging to V = v H (Ω), v = } and Find ω n+ W so tat, v W, (ω n+, v ) = ( ψ n+, v ). (6) are expressed as ψ ẑ by virtue of te te wellknown isomorpim (see Girault and Raviart [3])... ẑ : H (Ω) V. At variance wit more usual ways of writing te vorticity transport equation, no integration by parts is to be performed; in oter words, te curl of te momentum equation as not been taken in a strong form. In tis way te momentum equation becomes an evolutionary equation for te stream function (actually for its Laplacian in weak form) wereas te Poisson equation is used as te de nition of vorticity. Te nonlinear term is accounted for in a semi-implicit form for te sake of simplicity. All tat is said afterwards olds wit minor modi cations if tis term is made explicit. From te teoretical point of view, te modi cations in question essentially amount to deriving sligtly sarper bounds for te nonlinear residuals. Observe tat, wen compared to te classical ω-ψ formulation, te present metod intercanges te role of te variables ψ and ω. Here, te dynamical equation for te transport of ω as turned into an equation governing te evolution of (te weak form of) ψ, wereas te Poisson equation for ψ as become an expression controlling ω explicitly. Explicit ψ-ω metod 59 J.-L. Guermond and L. Quartapelle
3 . Vorticity integral conditions Note tat te explicit evaluation of te new vorticity eld ω n+ troug te solution of te mass matrix problem (6) does enforce te integral conditions for te vorticity wic underlay te uncoupled metod due to Glowinski and Pironneau []. In fact, considering more general, i.e., nonomogeneous, boundary conditions ψ n+ Γ = a n+ and ψn+ = b n+, n Γ te vorticity problem would read Find ω n+ W so tat, v W, (ω n+, v ) = ( ψ n+, v ) b n+ v. Selecting te functions v in te subspace of te discrete armonic functions, namely, v = η W so tat ( η, v ) =, v W, te weak equation above gives ( (ω n+, η ) a n+ η n ) bn+ η, Γ since it can be sown tat Γ an+ η / n ( ψ n+, η ). One recovers te vorticity integral conditions for te transient problem at time t n+. Tus, te proposed metod, wit te viscous diffusion made explicit, allows te vorticity integral conditions to be ful- lled a posteriori, as already pointed out in [6]. In te present formulation te vorticity boundary value is determined in a way tat is very similar to te classical procedure used in te context of nite differences. In fact te vorticity boundary formula used in second-order accurate central differences is obtained by means of a Taylor series expansion as follows ψ ( x) = ψ() + x ψ() x [ ] ( x) ω () + ψ() + O(( x) 3 ). y Tis argument uses te Poisson equation for ψ on te boundary togeter wit te Diriclet and Neumann boundary data for ψ. In some sense, te Taylor expansion above mimics te weak vorticity equation for a weigting function v suc tat v Γ. 5 STABILITY ANALYSIS Te following convergence result is establised in [5]: Γ Teorem Under convenient regularity ypotesis on te solution (ψ, ω) of te continuous problem (), tere is c s(ω) > and c e(t, ν, ψ, Ω) > so tat if δt c s /ν, ten ψ ψ l (H (Ω))+ ω ω l (L (Ω)) c e(δt+ l ). 6 SECOND ORDER BDF SCHEME Te present tecnique is not restricted to rst order; it can be modi ed to obtain ig order accuracy in time. Tis can be done simply by approximating te time derivative by a ig order nite differencing (Crank Nicolson, tree-level backward differencing, etc...) and by extrapolating te terms tat involve ω, accordingly. To illustrate tis possibility we present in te following a second order sceme based on te tree-level backward differencing of te time derivative and using a semi-implicit evaluation of te nonlinear term by means of linear extrapolation in time of te te vorticity. Initialize te sceme by evaluating (ψ, ω) and (ψ, ω). ψ and ω are evaluated from te initial data troug (3) and (). ψ can be obtained by many means; for instance, it can be calculated by using a second order Runge Kutta tecnique; from ψ one evaluates ω easily. Ten, for eac n, carry out te following two steps: Find ψ n+ Ψ, so tat, φ Ψ,, ( (3ψ n+ ψ n + ψ n ), φ ) δt + b(ω n ω n, ψ n+, φ ) = ν ( (ω n ω n ), φ ) + (f n+, φ ẑ), (7) and Find ω n+ W so tat, v W, (ω n+, v ) = ( ψ n+, v ). (8) Tis sceme is second order accurate in time. Its error analysis follows te same ideas as tose tat ave been used to analyze te rst order sceme and a uniform stability result can be found in [5]. 7 NUMERICAT TESTS We ave implemented te tree-level BDF sceme above in two different manners: rst by evaluating te nonlinear term in a fully explicit manner by time extrapolation Explicit ψ-ω metod 6 J.-L. Guermond and L. Quartapelle
4 e 7. e 8 e 8 e 7 e e 7 e 7 e 7 e e 7.. e 7 e 7.7 e Figure. Streamlines of FEM solution of te driven cavity problem for Re = at t = 6.5: fully explicit sceme wit direct solution of te symmetric linear systems. e 8 e. e 7 e e 7. e 8 5e 5 e e 7 e e 7 e 8 e e 7 e 8 Figure. Streamlines of FEM solution of te driven cavity problem for Re = at t = 6.5: semi-implicit sceme wit iterative solution of te linear systems..7 e. e 8 5e 5 e 7 e e e 7. e of bot ψ n+ and ω n+, leading to te weak equation, φ Ψ,, (3 ψ n+ /δt, φ ) = b(ω n ω n, ψ n ψ n, φ ) + g n (φ ), second, by evaluating te nonlinear term by means of a semi-implicit approximation wic gives a nonsymmetric contribution to te elliptic problem for ψ n+, as follows, (3 ψ n+ /δt, φ ) + b(ω n ω n, ψ n+, φ ) = g n (φ ). In bot cases te source term g n (φ ) is de ned by g n (φ ) = ( (ψ n ψ n /)/δt ν (ω n ω n ), φ ). Te rst sceme allows using direct algoritms for solving te two symmetric and time-independent linear sys- tems (stiffness and mass matrix) for ψ n+ and ω n+, respectively. Beside te mass matrix problem, te second sceme as a nonsymmetric system of linear equations to be solved at eac time step. As a consequence, we used te GMRES tecnique wit preconditioning based on incomplete factorization of te (constant) stiffness matrix. Te two algoritms ave been tested and compared by solving te unsteady driven cavity problem using nonuniform meses of 8 linear triangles, rst for Re =. In gures and we report te streamlines at time t = 6.5 calculated by te direct and te iterative FEM sceme using δt =. and δt =.5, respectively. Figures 3 and contains te vorticity distribution at te same time, to be compared wit te spectral solution calculated by a Galerkin spectral metod of Glowinski Pironneau type [] wit δt =. and 5 Legendre polynomials in eac direction, sown in gure 5. Te two FE metods are equally accurate, te sligt differ- Explicit ψ-ω metod 6 J.-L. Guermond and L. Quartapelle
5 ences between tem being caused by te use of different meses. Te second example is te calculation of te same problem for a moderately ig Reynolds number Re =. Te streamlines at time t = 5. provided by te explicit direct FEM sceme wit δt =. are sown in gure 6. Tis solution is in fair agreement wit a nite difference solution obtained on a uniform 5 grid and using a centered -accurate approximation of te fully explicit Jacobian reported in gure Figure 3. Vorticity of FEM solution of te driven cavity problem for Re = at t = 6.5: fully explicit sceme wit direct solution of te symmetric linear systems. 8 CONCLUSIONS In tis paper we ave presented a nite element sceme for solving te time-dependent Navier Stokes equations formulated in terms of te stream function and te vorticity. Te calculation of te stream function and te vorticity are uncoupled owing to an explicit treatment of te viscous diffusion togeter wit a non-standard writing of te evolutionary term in te weak form of te momentum equation. Te explicit treatment of viscous diffusion Explicit ψ-ω metod 6 J.-L. Guermond and L. Quartapelle
6 Figure. Vorticity of FEM solution of te driven cavity problem for Re = at t = 6.5: semi-implicit sceme wit iterative solution of te linear systems..5 Figure 5. Vorticity of Galerkin Legendre spectral solution of te driven cavity problem for Re = at t = 6.5: Glowinski Pironneau uncoupled metod. implies a stability condition of te type: νδt/ c. Tis stability constraint is te price to be paid for te extreme algoritmic simplicity of te proposed uncoupled sceme, especially wen compared to Glowinski Pironneau metod and related tecniques. Wile te stability restriction may be severe for creeping ows, te matter improves for convection dominated ows since te combination of te cell Reynolds number condition for adequate spatial resolution wit te stability condition gives a condition δt c. Te uncoupling strategy proposed in te paper is not limited to time discretizations of low order; a second order accurate sceme based on te tree-level backward difference formula as been presented. REFERENCES [] F. AUTERI AND L. QUARTAPELLE, Galerkin spectral metod for te ζ-ψ equations, submitted for publication, 998. [] W. E AND J.-G. LIU, Vorticity boundary conditions and related issues for nite differences scemes, J. Comput. Pys.,, 996, [3] V. GIRAULT AND P.-A. RAVIART, Finite Element Metods for Navier Stokes Equations, Springer Series in Computational Matematics, 5, Springer-Verlag, 986. [] R. GLOWINSKI AND O. PIRONNEAU, Numerical metods for te rst biarmonic equation and for te twodimensional Stokes problem, SIAM Review,, 979, 67. [5] J.-L. GUERMOND AND L. QUARTAPELLE, Weak approximation of te ψ-ω equations wit explicit viscous diffusion, submitted for publication, 998. [6] M. NAPOLITANO, G. PASCAZIO AND L. QUAR- TAPELLE, A review of vorticity conditions in te solution of te ζ-ψ equations, Computers and Fluids, to appear, 998. [7] R. PEYRET AND T.D. TAYLOR, Computational Metods for Fluid Flow, Springer-Verlag, Berlin, 983. Explicit ψ-ω metod 63 J.-L. Guermond and L. Quartapelle
7 Mon May 8 :: Figure 6. Streamlines of FEM solution of te driven cavity problem for Re = at t = 5.. Figure 7. Streamlines of FDM solution of te driven cavity problem for Re = at t = 5.. Explicit ψ-ω metod 6 J.-L. Guermond and L. Quartapelle
1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationNumerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems
Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for
More informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More information1 Introduction. J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods
J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods 1 Introduction Achieving high order time-accuracy in the approximation of the incompressible Navier Stokes equations by means of fractional-step
More informationA method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps
A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.
More informationJian-Guo Liu 1 and Chi-Wang Shu 2
Journal of Computational Pysics 60, 577 596 (000) doi:0.006/jcp.000.6475, available online at ttp://www.idealibrary.com on Jian-Guo Liu and Ci-Wang Su Institute for Pysical Science and Tecnology and Department
More informationGlowinski Pironneau method for the 3D ω-ψ equations
280 GUERMOND AND QUARTAPELLE Glowinski Pironneau method for the 3D ω-ψ equations Jean-Luc Guermond and Luigi Quartapelle 1 LIMSI CNRS, Orsay, France, and Dipartimento di Fisica, Politecnico di Milano,
More informationarxiv: v1 [math.na] 12 Mar 2018
ON PRESSURE ESTIMATES FOR THE NAVIER-STOKES EQUATIONS J A FIORDILINO arxiv:180304366v1 [matna 12 Mar 2018 Abstract Tis paper presents a simple, general tecnique to prove finite element metod (FEM) pressure
More informationSecond-order finite-element projection method for 3D flows
Second-order finite-element projection method for 3D flows J.-L. Guermond I and L. Quartapelle 2 1 LIMSI-CNRS, BP 133, 91403 Orsay, France, (guermond@limsi.fr) 2 Dipartimento di Fisica del Politecnico
More informationBlanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS
Opuscula Matematica Vol. 26 No. 3 26 Blanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS Abstract. In tis work a new numerical metod is constructed for time-integrating
More informationarxiv: v3 [math.na] 15 Dec 2009
THE NAVIER-STOKES-VOIGHT MODEL FOR IMAGE INPAINTING M.A. EBRAHIMI, MICHAEL HOLST, AND EVELYN LUNASIN arxiv:91.4548v3 [mat.na] 15 Dec 9 ABSTRACT. In tis paper we investigate te use of te D Navier-Stokes-Voigt
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationMANY scientific and engineering problems can be
A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial
More informationCLEMSON U N I V E R S I T Y
A Fractional Step θ-metod for Convection-Diffusion Equations Jon Crispell December, 006 Advisors: Dr. Lea Jenkins and Dr. Vincent Ervin Fractional Step θ-metod Outline Crispell,Ervin,Jenkins Motivation
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationSOME IMPLEMENTATIONS OF PROJECTION METHODS FOR NAVIER STOKES EQUATIONS
SOME IMPLEMENTATIONS OF PROJECTION METHODS FOR NAVIER STOKES EQUATIONS JEAN-LUC GUERMOND Abstract. Tis paper is concerned wit te implementation of spatially discrete versions of Corin Temam s projection
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationAN EFFICIENT AND ROBUST METHOD FOR SIMULATING TWO-PHASE GEL DYNAMICS
AN EFFICIENT AND ROBUST METHOD FOR SIMULATING TWO-PHASE GEL DYNAMICS GRADY B. WRIGHT, ROBERT D. GUY, AND AARON L. FOGELSON Abstract. We develop a computational metod for simulating models of gel dynamics
More informationUnconditional long-time stability of a velocity-vorticity method for the 2D Navier-Stokes equations
Numerical Analysis and Scientific Computing Preprint Seria Unconditional long-time stability of a velocity-vorticity metod for te D Navier-Stokes equations T. Heister M.A. Olsanskii L.G. Rebolz Preprint
More informationarxiv: v1 [math.na] 17 Jul 2014
Div First-Order System LL* FOSLL* for Second-Order Elliptic Partial Differential Equations Ziqiang Cai Rob Falgout Sun Zang arxiv:1407.4558v1 [mat.na] 17 Jul 2014 February 13, 2018 Abstract. Te first-order
More informationFinite Difference Methods Assignments
Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling Assignment 1: A one-dimensional eat equation
More informationGrad-div stabilization for the evolutionary Oseen problem with inf-sup stable finite elements
Noname manuscript No. will be inserted by te editor Grad-div stabilization for te evolutionary Oseen problem wit inf-sup stable finite elements Javier de Frutos Bosco García-Arcilla Volker Jon Julia Novo
More informationNumerical Solution to Parabolic PDE Using Implicit Finite Difference Approach
Numerical Solution to arabolic DE Using Implicit Finite Difference Approac Jon Amoa-Mensa, Francis Oene Boateng, Kwame Bonsu Department of Matematics and Statistics, Sunyani Tecnical University, Sunyani,
More informationAN OVERVIEW OF PROJECTION METHODS FOR INCOMPRESSIBLE FLOWS
AN OVERVIEW OF PROJECTION METHODS FOR INCOMPRESSIBLE FLOWS J.L. GUERMOND, P. MINEV 2, AND JIE SHEN 3 Abstract. We discuss in tis paper a series of important numerical issues related to te analysis and
More informationDepartment of Mathematical Sciences University of South Carolina Aiken Aiken, SC 29801
RESEARCH SUMMARY AND PERSPECTIVES KOFFI B. FADIMBA Department of Matematical Sciences University of Sout Carolina Aiken Aiken, SC 29801 Email: KoffiF@usca.edu 1. Introduction My researc program as focused
More information(4.2) -Richardson Extrapolation
(.) -Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Suppose tat lim G 0 and lim F L. Te function F is said to converge to L as
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationQuantum Numbers and Rules
OpenStax-CNX module: m42614 1 Quantum Numbers and Rules OpenStax College Tis work is produced by OpenStax-CNX and licensed under te Creative Commons Attribution License 3.0 Abstract Dene quantum number.
More informationDifferent Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method
WDS'10 Proceedings of Contributed Papers, Part I, 151 156, 2010. ISBN 978-80-7378-139-2 MATFYZPRESS Different Approaces to a Posteriori Error Analysis of te Discontinuous Galerkin Metod I. Šebestová Carles
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationSimulations of the turbulent channel flow at Re τ = 180 with projection-based finite element variational multiscale methods
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Met. Fluids 7; 55:47 49 Publised online 4 Marc 7 in Wiley InterScience (www.interscience.wiley.com). DOI:./fld.46 Simulations of te
More informationarxiv: v1 [math.ap] 4 Aug 2017
New Two Step Laplace Adam-Basfort Metod for Integer an Non integer Order Partial Differential Equations arxiv:178.1417v1 [mat.ap] 4 Aug 17 Abstract Rodrigue Gnitcogna*, Abdon Atangana** *Department of
More informationCS522 - Partial Di erential Equations
CS5 - Partial Di erential Equations Tibor Jánosi April 5, 5 Numerical Di erentiation In principle, di erentiation is a simple operation. Indeed, given a function speci ed as a closed-form formula, its
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 34, pp. 14-19, 2008. Copyrigt 2008,. ISSN 1068-9613. ETNA A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES FOR HIGH INDEX DIFFERENTIAL-ALGEBRAIC EQUATIONS
More informationEntropy and the numerical integration of conservation laws
Pysics Procedia Pysics Procedia 00 2011) 1 28 Entropy and te numerical integration of conservation laws Gabriella Puppo Dipartimento di Matematica, Politecnico di Torino Italy) Matteo Semplice Dipartimento
More informationInf sup testing of upwind methods
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Met. Engng 000; 48:745 760 Inf sup testing of upwind metods Klaus-Jurgen Bate 1; ;, Dena Hendriana 1, Franco Brezzi and Giancarlo
More informationA Hybridizable Discontinuous Galerkin Method for the Compressible Euler and Navier-Stokes Equations
A Hybridizable Discontinuous Galerkin Metod for te Compressible Euler and Navier-Stokes Equations J. Peraire and N. C. Nguyen Massacusetts Institute of Tecnology, Cambridge, MA 02139, USA B. Cockburn University
More informationNew Fourth Order Quartic Spline Method for Solving Second Order Boundary Value Problems
MATEMATIKA, 2015, Volume 31, Number 2, 149 157 c UTM Centre for Industrial Applied Matematics New Fourt Order Quartic Spline Metod for Solving Second Order Boundary Value Problems 1 Osama Ala yed, 2 Te
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationFINITE ELEMENT DUAL SINGULAR FUNCTION METHODS FOR HELMHOLTZ AND HEAT EQUATIONS
J. KSIAM Vol.22, No.2, 101 113, 2018 ttp://dx.doi.org/10.12941/jksiam.2018.22.101 FINITE ELEMENT DUAL SINGULAR FUNCTION METHODS FOR HELMHOLTZ AND HEAT EQUATIONS DEOK-KYU JANG AND JAE-HONG PYO DEPARTMENT
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationarxiv: v1 [math.na] 28 Apr 2017
THE SCOTT-VOGELIUS FINITE ELEMENTS REVISITED JOHNNY GUZMÁN AND L RIDGWAY SCOTT arxiv:170500020v1 [matna] 28 Apr 2017 Abstract We prove tat te Scott-Vogelius finite elements are inf-sup stable on sape-regular
More informationImplicit-explicit variational integration of highly oscillatory problems
Implicit-explicit variational integration of igly oscillatory problems Ari Stern Structured Integrators Worksop April 9, 9 Stern, A., and E. Grinspun. Multiscale Model. Simul., to appear. arxiv:88.39 [mat.na].
More informationNumerical analysis of a free piston problem
MATHEMATICAL COMMUNICATIONS 573 Mat. Commun., Vol. 15, No. 2, pp. 573-585 (2010) Numerical analysis of a free piston problem Boris Mua 1 and Zvonimir Tutek 1, 1 Department of Matematics, University of
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More information7 Semiparametric Methods and Partially Linear Regression
7 Semiparametric Metods and Partially Linear Regression 7. Overview A model is called semiparametric if it is described by and were is nite-dimensional (e.g. parametric) and is in nite-dimensional (nonparametric).
More informationA First-Order System Approach for Diffusion Equation. I. Second-Order Residual-Distribution Schemes
A First-Order System Approac for Diffusion Equation. I. Second-Order Residual-Distribution Scemes Hiroaki Nisikawa W. M. Keck Foundation Laboratory for Computational Fluid Dynamics, Department of Aerospace
More informationCONVERGENCE ANALYSIS OF FINITE ELEMENT SOLUTION OF ONE-DIMENSIONAL SINGULARLY PERTURBED DIFFERENTIAL EQUATIONS ON EQUIDISTRIBUTING MESHES
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume, Number, Pages 57 74 c 5 Institute for Scientific Computing and Information CONVERGENCE ANALYSIS OF FINITE ELEMENT SOLUTION OF ONE-DIMENSIONAL
More informationFlavius Guiaş. X(t + h) = X(t) + F (X(s)) ds.
Numerical solvers for large systems of ordinary differential equations based on te stocastic direct simulation metod improved by te and Runge Kutta principles Flavius Guiaş Abstract We present a numerical
More informationVariational Localizations of the Dual Weighted Residual Estimator
Publised in Journal for Computational and Applied Matematics, pp. 192-208, 2015 Variational Localizations of te Dual Weigted Residual Estimator Tomas Ricter Tomas Wick Te dual weigted residual metod (DWR)
More informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationOn the accuracy of the rotation form in simulations of the Navier-Stokes equations
On te accuracy of te rotation form in simulations of te Navier-Stokes equations William Layton 1 Carolina C. Manica Monika Neda 3 Maxim Olsanskii Leo G. Rebolz 5 Abstract Te rotation form of te Navier-Stokes
More informationCrouzeix-Velte Decompositions and the Stokes Problem
Crouzeix-Velte Decompositions and te Stokes Problem PD Tesis Strauber Györgyi Eötvös Loránd University of Sciences, Insitute of Matematics, Matematical Doctoral Scool Director of te Doctoral Scool: Dr.
More informationSuperconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract
Superconvergence of energy-conserving discontinuous Galerkin metods for linear yperbolic equations Yong Liu, Ci-Wang Su and Mengping Zang 3 Abstract In tis paper, we study superconvergence properties of
More information1. Introduction. Consider a semilinear parabolic equation in the form
A POSTERIORI ERROR ESTIMATION FOR PARABOLIC PROBLEMS USING ELLIPTIC RECONSTRUCTIONS. I: BACKWARD-EULER AND CRANK-NICOLSON METHODS NATALIA KOPTEVA AND TORSTEN LINSS Abstract. A semilinear second-order parabolic
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationEfficient, unconditionally stable, and optimally accurate FE algorithms for approximate deconvolution models of fluid flow
Efficient, unconditionally stable, and optimally accurate FE algoritms for approximate deconvolution models of fluid flow Leo G. Rebolz Abstract Tis paper addresses an open question of ow to devise numerical
More informationCOMPUTATIONAL COMPARISON BETWEEN THE TAYLOR HOOD AND THE CONFORMING CROUZEIX RAVIART ELEMENT
Proceedings of ALGORITMY 5 pp. 369 379 COMPUTATIONAL COMPARISON BETWEEN E TAYLOR HOOD AND E CONFORMING OUZEIX RAVIART ELEMENT ROLF KRAHL AND EBERHARD BÄNSCH Abstract. Tis paper is concerned wit te computational
More informationOptimal Shape Design of a Two-dimensional Asymmetric Diffsuer in Turbulent Flow
THE 5 TH ASIAN COMPUTAITIONAL FLUID DYNAMICS BUSAN, KOREA, OCTOBER 7 ~ OCTOBER 30, 003 Optimal Sape Design of a Two-dimensional Asymmetric Diffsuer in Turbulent Flow Seokyun Lim and Haeceon Coi. Center
More informationSTABILITY OF DISCRETE STOKES OPERATORS IN FRACTIONAL SOBOLEV SPACES
STABILITY OF DISCRETE STOKES OPERATORS IN FRACTIONAL SOBOLEV SPACES JEAN-LUC GUERMOND,, JOSEPH E PASCIAK Abstract Using a general approximation setting aving te generic properties of finite-elements, we
More informationSMAI-JCM SMAI Journal of Computational Mathematics
SMAI-JCM SMAI Journal of Computational Matematics Compatible Maxwell solvers wit particles II: conforming and non-conforming 2D scemes wit a strong Faraday law Martin Campos Pinto & Eric Sonnendrücker
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Matematics 94 (6) 75 96 Contents lists available at ScienceDirect Journal of Computational and Applied Matematics journal omepage: www.elsevier.com/locate/cam Smootness-Increasing
More informationSome Error Estimates for the Finite Volume Element Method for a Parabolic Problem
Computational Metods in Applied Matematics Vol. 13 (213), No. 3, pp. 251 279 c 213 Institute of Matematics, NAS of Belarus Doi: 1.1515/cmam-212-6 Some Error Estimates for te Finite Volume Element Metod
More informationExplicit Hyperbolic Reconstructed Discontinuous Galerkin Methods for Time-Dependent Problems
AIAA AVIATION Forum June 25-29 218 Atlanta Georgia 218 Fluid Dynamics Conference 1.2514/6.218-427 Explicit Hyperbolic Reconstructed Discontinuous Galerkin Metods for Time-Dependent Problems Jialin Lou
More information4.2 - Richardson Extrapolation
. - Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Definition Let x n n converge to a number x. Suppose tat n n is a sequence
More informationOverlapping domain decomposition methods for elliptic quasi-variational inequalities related to impulse control problem with mixed boundary conditions
Proc. Indian Acad. Sci. (Mat. Sci.) Vol. 121, No. 4, November 2011, pp. 481 493. c Indian Academy of Sciences Overlapping domain decomposition metods for elliptic quasi-variational inequalities related
More informationAN ANALYSIS OF NEW FINITE ELEMENT SPACES FOR MAXWELL S EQUATIONS
Journal of Matematical Sciences: Advances and Applications Volume 5 8 Pages -9 Available at ttp://scientificadvances.co.in DOI: ttp://d.doi.org/.864/jmsaa_7975 AN ANALYSIS OF NEW FINITE ELEMENT SPACES
More informationSome Applications of Fractional Step Runge-Kutta Methods
Some Applications of Fractional Step Runge-Kutta Metods JORGE, J.C., PORTERO, L. Dpto. de Matematica e Informatica Universidad Publica de Navarra Campus Arrosadia, s/n 3006 Pamplona Navarra SPAIN Abstract:
More informationA Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Hybrid Mixed Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationStability properties of a family of chock capturing methods for hyperbolic conservation laws
Proceedings of te 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August 0-, 005 (pp48-5) Stability properties of a family of cock capturing metods for yperbolic conservation
More informationarxiv: v1 [physics.flu-dyn] 3 Jun 2015
A Convective-like Energy-Stable Open Boundary Condition for Simulations of Incompressible Flows arxiv:156.132v1 [pysics.flu-dyn] 3 Jun 215 S. Dong Center for Computational & Applied Matematics Department
More informationAdvancements In Finite Element Methods For Newtonian And Non-Newtonian Flows
Clemson University TigerPrints All Dissertations Dissertations 8-3 Advancements In Finite Element Metods For Newtonian And Non-Newtonian Flows Keit Galvin Clemson University, kjgalvi@clemson.edu Follow
More informationA Projection FEM for Variable Density Incompressible Flows
Journal of Computational Physics 65, 67 88 (000) doi:0.006/jcph.000.6609, available online at http://www.idealibrary.com on A Projection FEM for Variable Density Incompressible Flows J.-L. Guermond and
More informationH(div) conforming and DG methods for incompressible Euler s equations
H(div) conforming and DG metods for incompressible Euler s equations Jonny Guzmán Filánder A. Sequeira Ci-Wang Su Abstract H(div) conforming and discontinuous Galerkin (DG) metods are designed for incompressible
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationNotes on Multigrid Methods
Notes on Multigrid Metods Qingai Zang April, 17 Motivation of multigrids. Te convergence rates of classical iterative metod depend on te grid spacing, or problem size. In contrast, convergence rates of
More informationMIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp
MIXED DISCONTINUOUS GALERIN APPROXIMATION OF THE MAXWELL OPERATOR PAUL HOUSTON, ILARIA PERUGIA, AND DOMINI SCHÖTZAU SIAM J. Numer. Anal., Vol. 4 (004), pp. 434 459 Abstract. We introduce and analyze a
More information27. A Non-Overlapping Optimized Schwarz Method which Converges with Arbitrarily Weak Dependence on h
Fourteent International Conference on Domain Decomposition Metods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 003 DDM.org 7. A Non-Overlapping Optimized Scwarz Metod wic Converges
More informationarxiv: v1 [math.na] 25 Jul 2014
A second order in time, uniquely solvable, unconditionally stable numerical sceme for Can-Hilliard-Navier-Stokes equation Daozi Han, Xiaoming Wang November 5, 016 arxiv:1407.7048v1 [mat.na] 5 Jul 014 Abstract
More informationFourier Type Super Convergence Study on DDGIC and Symmetric DDG Methods
DOI 0.007/s095-07-048- Fourier Type Super Convergence Study on DDGIC and Symmetric DDG Metods Mengping Zang Jue Yan Received: 7 December 06 / Revised: 7 April 07 / Accepted: April 07 Springer Science+Business
More informationA SADDLE POINT LEAST SQUARES APPROACH TO MIXED METHODS
A SADDLE POINT LEAST SQUARES APPROACH TO MIXED METHODS CONSTANTIN BACUTA AND KLAJDI QIRKO Abstract. We investigate new PDE discretization approaces for solving variational formulations wit different types
More informationStokes Equation. Chapter Mathematical Formulations
Capter 5 Stokes Equation 5.1 Matematical Formulations We introduce some notations: Let v = (v 1,v 2 ) T or v = (v 1,v 2,v 3 ) T. Wen n = 2 we define curlv = v 2 x v 1 y, ) curlη =. ( η y, η x In fact,
More informationDIRECT NUMERICAL SIMULATION IN A LID-DRIVEN CAVITY AT HIGH REYNOLDS NUMBER
Conference on Turbulence and Interactions TI26, May 29 - June 2, 26, Porquerolles, France DIRECT NUMERICAL SIMULATION IN A LID-DRIVEN CAVITY AT HIGH REYNOLDS NUMBER E. Leriche, Laboratoire d Ingénierie
More informationA Mixed-Hybrid-Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Mixed-Hybrid-Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More informationControllability of a one-dimensional fractional heat equation: theoretical and numerical aspects
Controllability of a one-dimensional fractional eat equation: teoretical and numerical aspects Umberto Biccari, Víctor Hernández-Santamaría To cite tis version: Umberto Biccari, Víctor Hernández-Santamaría.
More informationCONVERGENCE ANALYSIS OF YEE SCHEMES FOR MAXWELL S EQUATIONS IN DEBYE AND LORENTZ DISPERSIVE MEDIA
INTRNATIONAL JOURNAL OF NUMRICAL ANALYSIS AND MODLING Volume Number 4 ages 657 687 c 04 Institute for Scientific Computing and Information CONVRGNC ANALYSIS OF Y SCHMS FOR MAXWLL S QUATIONS IN DBY AND
More informationCONVERGENCE ANALYSIS OF YEE SCHEMES FOR MAXWELL S EQUATIONS IN DEBYE AND LORENTZ DISPERSIVE MEDIA
INTRNATIONAL JOURNAL OF NUMRICAL ANALYSIS AND MODLING Volume XX Number 0 ages 45 c 03 Institute for Scientific Computing and Information CONVRGNC ANALYSIS OF Y SCHMS FOR MAXWLL S QUATIONS IN DBY AND LORNTZ
More informationDifferential equations. Differential equations
Differential equations A differential equation (DE) describes ow a quantity canges (as a function of time, position, ) d - A ball dropped from a building: t gt () dt d S qx - Uniformly loaded beam: wx
More information3 Parabolic Differential Equations
3 Parabolic Differential Equations 3.1 Classical solutions We consider existence and uniqueness results for initial-boundary value problems for te linear eat equation in Q := Ω (, T ), were Ω is a bounded
More informationLEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS
SIAM J. NUMER. ANAL. c 998 Society for Industrial Applied Matematics Vol. 35, No., pp. 393 405, February 998 00 LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS YANZHAO CAO
More informationA SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION
A SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION GABRIEL R. BARRENECHEA, LEOPOLDO P. FRANCA 1 2, AND FRÉDÉRIC VALENTIN Abstract. Tis work introduces and analyzes novel stable
More informationThese errors are made from replacing an infinite process by finite one.
Introduction :- Tis course examines problems tat can be solved by metods of approximation, tecniques we call numerical metods. We begin by considering some of te matematical and computational topics tat
More informationMass Lumping for Constant Density Acoustics
Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial
More information