Theoretical advances. To illustrate our approach, consider the scalar ODE model,
|
|
- Anabel Easter Fields
- 5 years ago
- Views:
Transcription
1 Final Report GR/R69/ : Analysis of Numerical Methods for Incompressible Fluid Dynamics Personnel supported: Professor Philip Gresho (nominated Visiting Fellow; visited UK: 7/6/ 5//). Dr David Kay (visited UMIST: /7/ 6/7/; /9/ /9/; /6/ /6/ ) Introduction. Simulation of the motion of an incompressible fluid remains an important but very challenging problem. The resources required for accurate three-dimensional simulation of practical flows test even the most advanced of supercomputer hardware. The necessity for reliable and efficient solvers is widely recognised. The aim of this project was to study the performance of adaptive time-stepping methods and state-of-the-art iterative solvers in the context of numerical solution methods for the Navier-Stokes equations. In both research threads we have made important advances commensurate with those envisaged on the original proposal these are detailed below. The research effort was initiated and stimulated by Phil Gresho s six month visit to the UK. This visit was split roughly evenly between Dundee and Manchester, and discussions between Silvester, Griffiths and Gresho have continued electronically since Gresho s return to California. The main outcome of these technical exchanges has been a detailed analytic study of a particular adaptive time-stepping scheme that is a candidate best method in the context of low order finite element approximation in space. This work is currently being written up (in []), and a sample set of results is reviewed below. In parallel with this theoretical work, Silvester and Kay took the opportunity to discuss a variety of issues with Gresho associated with practical implementation of the adaptive time-stepping methodology. As a result, the approach has been implemented into a research code that is being developed by Kay at the University of Sussex. Whilst testing of this code is ongoing, the effectiveness of the methodology, in combination with our fast iterative solver technology, is already apparent. A discussion of practical issues is given in []. Theoretical advances. To illustrate our approach, consider the scalar ODE model,. y = λy y() = y ; (a) (b) with the simple solution y(t) = y e λt, () where λ is generally complex: λ = τ + i ω. (3) Here τ corresponds to a decay time constant and is meant to model simple diffusion, whereas the frequency, ω, is a simple model for advection. The combination of the two corresponds to advection-diffusion. The basic ODE method for our smart integration is the well known (second-order accurate, implicit) trapezoid rule (TR) chosen for several reasons: It is unconditionally stable (A-stable) stable when the ODE is. It is non-dissipative important when modelling pure advection (τ = ), or even advectiondominated problems ( τ ω). It is self starting.
2 There are no extraneous/spurious roots. Also relevant are the following facts: There are no explicit linear multi-step methods that are A-stable. The highest order A-stable method is two and the method is necessarily implicit. Of all second-order, A-stable methods, TR is the most accurate it has the smallest local truncation error coefficient. Applied to any ODE, the Trapezoid rule is simply y n+ = y n + t n ( ẏn + ẏ n+) where the ẏ n+ term on the right-hand-side makes the integrator implicit. Substituting the ODE for ẏ from (a) gives yn+ = λ t n + λ t n (4) y n. (5) At this point, if t n = t = constant, (5) merely represents the fixed- t version ( dumb version) of the integrator. Our real goal is a smart integrator that automatically changes t such that the local error is bounded for all time by a user-specified tolerance. To make progress towards this goal, we start from the truncation error term d n := y n+ y (t n+ ) = t3 n... y n +O ( t 4 ) n where y(t n+ ) is the exact solution at t = t n+, and match-up the Trapezoid rule with an appropriate explicit method that is the same order in terms of local truncation error. One possibility here is the second order Adams-Bashforth method, AB: ( y p n+ = y n + t n ẏ n + t n ẏn ẏ ) n, (7) t n where the superscript p signifies that it is a predictor it predicts a good estimate (for nonlinear equations) for y n+, the final (corrected) Trapezoid rule result. The basic idea here is that y p n+ leads to a technique that hopefully gives a good estimate of the otherwise unknown... y n+ in (6), from which local error control is easily developed. The first problem that we looked at is the heat equation (in R ): the aim is to compute u(x, t) satisfying the parabolic PDE u t = ν u x (8) on the unit interval, with zero initial data u(x, ), and with boundary data u(, t) =, u(, t) =. The motivation for studying (8) is that the PDE solution u can be analytically computed thus enabling us to compare numerical solutions with true physics. Taking a uniform linear finite element discretization in space (so that h = /N, with N = 64 (black) and N = 8 (red) representing the number of elements), the behaviour of the TR AB integrator with error tolerance tol = 5 is illustrated in Fig.. Notice the three distinct phases of the time-step evolution: The first phase is a rapid growth from the intial time-step t = 9, to a level fixed by the dominant eigenmode of the PDE problem. Specifically, t h exp(4t/3h )( tol) /3 (indicated with a solid line). (6)
3 Fig.. Heat equation ν = : log t vs. log t A second phase where the time-step follows the PDE dynamics; specifically t grows like t / independent of h. This is consistent with the notion of parabolic smoothing. A final phase (also independent of h) is where the time-step t (/π ) exp(π t/3)( tol) /3 (indicated in green), associated with the smallest eigenmode of the PDE problem. Note that if error control is tightened (by reducing tol) then the same qualitative behaviour is observed except that the timestep is reduced by a factor tol /3. Having established that the TR AB integrator is able to follow the physics of a diffusion problem, we moved on to the tougher case of advection-diffusion (again in R ): find u(x, t) satisfying u t = ν u x u x, (9) together with zero boundary data u(, t) = and u(, t) =. Although the exact form of the solution u is not available in this case, insight is nevertheless possible using classical asymptotic techniques (shown to us by a leading expert, Professor John Dold). A good example of what is possible is illustrated in Fig., where numerical and asymptotic solutions are compared in the case of a Gaussian wave profile being diffused and advected from left to right, prior to hitting a wall at the right hand boundary. We use a standard linear finite element discretization defined on a so-called Shishkin grid. (This is made up of two uniform grids of 64 elements, with a grid spacing chosen so that approximately half of the nodes lie in the boundary layer.) The agreement between numerics (dotted) and asymptotics (solid line) is impressive. The behaviour of our smart integrator with an error tolerance tol = 4 is illustrated in Fig. 3. There are four distinct phases of the time-step evolution in this case. Rapid growth from the intial time-step, t =, to a level fixed by the error tolerance and the discretization size h. A second phase where the time-step follows the (small-time) diffusion time-scale dynamics. As in Fig., t grows like t / independent of h (the straight line has slope /). A third phase (also independent of h) where the time-step is governed by the convection timescale dynamics. Here the profile is being convected from left to right, and the time-step stays essentially constant until the profile is wiped out at the right-hand boundary. 3
4 Fig.. Splat solution. ν = 4 : tol = 4 : t =.77. t = x BL solution x Fig. 3. Advection-diffusion equation ν = 4 : log t vs. log t. 3 4 t A final phase where the time-step rapidly grows (where the solution is essentially zero). Once again it is evident that the TR AB integrator is following the physics. An important issue that is considered in [] is the interplay between the temporal and the spatial discretization. Specifically, spectral properties of the discrete operators on stretched grids (favoured by Engineers) and Bakhvalov-Shishkin grids (favoured by mathematicians) have proved to be highly thought provoking both from an approximation viewpoint and a linear algebra perspective (the matrices are highly non-normal). We anticipate that the results in [] will be a fertile ground for both of these research communities in the future. Practical advances. Efficient implementation of the TR AB integrator into a code for solving the incompressible Navier-Stokes equation is the final aim of the project. Conventional codes typically use semi-implicit time integration leading to a Poisson or Stokes-type system at every time-step, but with a stability restriction on the time-step. In contrast, there is no such time-step restriction in our case. The price that must be paid for this improved robustness is the need to solve a so-called Oseen system (representing the linearised Navier-Stokes equations) at every time-step, see []. Initiated by David Kay s PDRA research (EPSRC grants GR/K96 and GR/L567), optimally efficient solvers for such Oseen problems have become a reality in the the last five years. (This is ongoing joint research involving Kay, Silvester and Howard Elman 4 t
5 from the University of Maryland, together with Daniel Loghin and Andy Wathen from Oxford University Computing Laboratory, see [3].) Specifically the preconditioning framework that we have developed offers the possibility of uniformly fast convergence independent of the problem parameters (namely; the mesh size, the time step and the Reynolds number). In contrast, conventional multigrid approaches to this problem tend to work well if the time-step is sufficiently small, but efficiency deteriorates rapidly if either the time-step or the Reynolds number is increased. Fig. 4. Snapshot of vortex shedding Re 4: tol = 3. Fig. 5. Effect of error tolerance on the flow solution. 3 Time step tol =.5 Time step tol = 3.5 Time step tol = 4 X Velocity Y Velocity Pressure
6 Guided by Phil Gresho s experience, we have tested Kay s code on a range of flow problems. One example is that of a flow being forced around a cylinder and ultimately reaching a periodic state of vortex shedding. (A snap-shot of the numerical solution in the shedding state is illustrated in Fig. 4.) Specific issues that we have addressed include alternative linearization approaches; the stability of the spatial discretization (adding streamline diffusion improves solver performance but can impinge on solution accuracy); mesh quality (big jumps in mesh-size need to be avoided in critical areas near the cylinder); and solver robustness with respect to the mesh-size and Reynolds number. Another point to note is that there is a delicate balance between the stopping criteria for the linear equation solver and the error control tolerance. The really fundamental issue is the need to balance spatial and temporal errors. This is illustrated by the results in Fig. 5, showing a point value of the velocity and pressure (on a fixed spatial mesh) with three choices of adaptive time-stepping tolerance. (Ideally we would like to automatically compute a tolerance so as to ensure that spatial error dominates.) If the error tolerance is not small enough ( 4 ) then the shedding frequency is wrong (cf. Fig. 4)! Dissemination aspects. Invited Talks given by David Silvester: Fast and robust solvers for time-discretised incompressible Navier-Stokes equations, Symposium on Finite Element Methods, Technical University Chemnitz, Germany, September ; Fast black-box preconditioners for self-adjoint PDE problems, Minisymposium on Numerical Linear Algebra, BAMC Meeting, University of Warwick, April, 5th Householder Symposium, Peebles, Scotland, June. Research Seminars given by David Silvester: Incompressible flow modelling is a dodgy business, University of Sussex, November ; A posteriori error estimation for elliptic PDEs, University of Manchester, February. Research Seminars given by Philip Gresho: Thermal convection in an 8: sidewall-heated enclosure, Scottish Computational Mathematics Symposium, Glasgow, September. On the stability of a particular steady incompressible flow at Reynolds number 8, University of Manchester, October ; A -D time periodic thermal convection problem, University of Manchester, October ; Some interesting aspects of a common mathematically ill-posed problem: impulsive starts of an incompressible fluid, UMIST, October ; How do to better CFD status and potential future prospects/directions for solving the incompressible Navier-Stokes equations, UMIST, November. (These seminars were widely advertised and well-attended. For the abstracts, see Publications: GR/R69/. Gresho, P., Griffiths, D. and Silvester, D. Adaptive time-stepping methods for solving incompressible CFD problems; part I: theoretical motivation, in preparation.. Kay, D. and Silvester, D. Adaptive time-stepping methods for solving incompressible CFD problems; part II: practical performance, in preparation. 3. Wathen, A., Loghin, D, Kay, D., Elman, H., and Silvester, D., A preconditioner for the 3D Oseen equations, Oxford University Computing Laboratory Report # /4,. 6
Fast solvers for steady incompressible flow
ICFD 25 p.1/21 Fast solvers for steady incompressible flow Andy Wathen Oxford University wathen@comlab.ox.ac.uk http://web.comlab.ox.ac.uk/~wathen/ Joint work with: Howard Elman (University of Maryland,
More informationFast Solvers for Unsteady Incompressible Flow
ICFD25 Workshop p. 1/39 Fast Solvers for Unsteady Incompressible Flow David Silvester University of Manchester http://www.maths.manchester.ac.uk/~djs/ ICFD25 Workshop p. 2/39 Outline Semi-Discrete Navier-Stokes
More informationA Review of Solution Techniques for Unsteady Incompressible Flow
Zeist 2009 p. 1/57 A Review of Solution Techniques for Unsteady Incompressible Flow David Silvester School of Mathematics University of Manchester Zeist 2009 p. 2/57 PDEs Review : 1966 1999 Update : 2000
More informationA Review of Preconditioning Techniques for Steady Incompressible Flow
Zeist 2009 p. 1/43 A Review of Preconditioning Techniques for Steady Incompressible Flow David Silvester School of Mathematics University of Manchester Zeist 2009 p. 2/43 PDEs Review : 1984 2005 Update
More informationFast Iterative Solution of Saddle Point Problems
Michele Benzi Department of Mathematics and Computer Science Emory University Atlanta, GA Acknowledgments NSF (Computational Mathematics) Maxim Olshanskii (Mech-Math, Moscow State U.) Zhen Wang (PhD student,
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationA Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers
Applied and Computational Mathematics 2017; 6(4): 202-207 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20170604.18 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Robust Preconditioned
More informationNavier-Stokes equations
1 Navier-Stokes equations Introduction to spectral methods for the CSC Lunchbytes Seminar Series. Incompressible, hydrodynamic turbulence is described completely by the Navier-Stokes equations where t
More informationUncertainty Quantification: Does it need efficient linear algebra?
USA trip October 2014 p. 1/62 Uncertainty Quantification: Does it need efficient linear algebra? David Silvester University of Manchester, UK Catherine Powell University of Manchester,UK Yes. USA trip
More informationPALADINS: Scalable Time-Adaptive Algebraic Splitting and Preconditioners for the Navier-Stokes Equations
2013 SIAM Conference On Computational Science and Engineering Boston, 27 th February 2013 PALADINS: Scalable Time-Adaptive Algebraic Splitting and Preconditioners for the Navier-Stokes Equations U. Villa,
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationON THE ROLE OF COMMUTATOR ARGUMENTS IN THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS FOR STOKES CONTROL PROBLEMS
ON THE ROLE OF COUTATOR ARGUENTS IN THE DEVELOPENT OF PARAETER-ROBUST PRECONDITIONERS FOR STOKES CONTROL PROBLES JOHN W. PEARSON Abstract. The development of preconditioners for PDE-constrained optimization
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationAdditive Manufacturing Module 8
Additive Manufacturing Module 8 Spring 2015 Wenchao Zhou zhouw@uark.edu (479) 575-7250 The Department of Mechanical Engineering University of Arkansas, Fayetteville 1 Evaluating design https://www.youtube.com/watch?v=p
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationSOE3213/4: CFD Lecture 3
CFD { SOE323/4: CFD Lecture 3 @u x @t @u y @t @u z @t r:u = 0 () + r:(uu x ) = + r:(uu y ) = + r:(uu z ) = @x @y @z + r 2 u x (2) + r 2 u y (3) + r 2 u z (4) Transport equation form, with source @x Two
More informationEfficient Solvers for the Navier Stokes Equations in Rotation Form
Efficient Solvers for the Navier Stokes Equations in Rotation Form Computer Research Institute Seminar Purdue University March 4, 2005 Michele Benzi Emory University Atlanta, GA Thanks to: NSF (MPS/Computational
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 11 Partial Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002.
More informationSimulating Solid Tumor Growth using Multigrid Algorithms
Simulating Solid Tumor Growth using Multigrid Applied Mathematics, Statistics, and Program Advisor: Doron Levy, PhD Department of Mathematics/CSCAMM University of Maryland, College Park September 30, 2014
More informationFinal Report: DE-FG02-95ER25239 Spectral Representations of Uncertainty: Algorithms and Applications
Final Report: DE-FG02-95ER25239 Spectral Representations of Uncertainty: Algorithms and Applications PI: George Em Karniadakis Division of Applied Mathematics, Brown University April 25, 2005 1 Objectives
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 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 informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationHessian-based model reduction for large-scale systems with initial condition inputs
Hessian-based model reduction for large-scale systems with initial condition inputs O. Bashir 1, K. Willcox 1,, O. Ghattas 2, B. van Bloemen Waanders 3, J. Hill 3 1 Massachusetts Institute of Technology,
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationSimulating Solid Tumor Growth Using Multigrid Algorithms
Simulating Solid Tumor Growth Using Multigrid Algorithms Asia Wyatt Applied Mathematics, Statistics, and Scientific Computation Program Advisor: Doron Levy Department of Mathematics/CSCAMM Abstract In
More informationChapter 6 - Ordinary Differential Equations
Chapter 6 - Ordinary Differential Equations 7.1 Solving Initial-Value Problems In this chapter, we will be interested in the solution of ordinary differential equations. Ordinary differential equations
More informationMultigrid and Iterative Strategies for Optimal Control Problems
Multigrid and Iterative Strategies for Optimal Control Problems John Pearson 1, Stefan Takacs 1 1 Mathematical Institute, 24 29 St. Giles, Oxford, OX1 3LB e-mail: john.pearson@worc.ox.ac.uk, takacs@maths.ox.ac.uk
More informationEfficient Solvers for Stochastic Finite Element Saddle Point Problems
Efficient Solvers for Stochastic Finite Element Saddle Point Problems Catherine E. Powell c.powell@manchester.ac.uk School of Mathematics University of Manchester, UK Efficient Solvers for Stochastic Finite
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
More informationMathematics and Computer Science
Technical Report TR-2007-002 Block preconditioning for saddle point systems with indefinite (1,1) block by Michele Benzi, Jia Liu Mathematics and Computer Science EMORY UNIVERSITY International Journal
More information5.1 2D example 59 Figure 5.1: Parabolic velocity field in a straight two-dimensional pipe. Figure 5.2: Concentration on the input boundary of the pipe. The vertical axis corresponds to r 2 -coordinate,
More informationc 2006 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 27, No. 5, pp. 1651 1668 c 2006 Society for Industrial and Applied Mathematics BLOCK PRECONDITIONERS BASED ON APPROXIMATE COMMUTATORS HOWARD ELMAN, VICTORIA E. HOWLE, JOHN SHADID,
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 informationFinding Rightmost Eigenvalues of Large, Sparse, Nonsymmetric Parameterized Eigenvalue Problems
Finding Rightmost Eigenvalues of Large, Sparse, Nonsymmetric Parameterized Eigenvalue Problems AMSC 663-664 Final Report Minghao Wu AMSC Program mwu@math.umd.edu Dr. Howard Elman Department of Computer
More informationComputational challenges in Numerical Weather Prediction
Computational challenges in Numerical Weather Prediction Mike Cullen Oxford 15 September 2008 Contents This presentation covers the following areas Historical background Current challenges Why does it
More informationCONVERGENCE BOUNDS FOR PRECONDITIONED GMRES USING ELEMENT-BY-ELEMENT ESTIMATES OF THE FIELD OF VALUES
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 CONVERGENCE BOUNDS FOR PRECONDITIONED GMRES USING ELEMENT-BY-ELEMENT
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
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 informationCS520: numerical ODEs (Ch.2)
.. CS520: numerical ODEs (Ch.2) Uri Ascher Department of Computer Science University of British Columbia ascher@cs.ubc.ca people.cs.ubc.ca/ ascher/520.html Uri Ascher (UBC) CPSC 520: ODEs (Ch. 2) Fall
More informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationModule 4: Numerical Methods for ODE. Michael Bader. Winter 2007/2008
Outlines Module 4: for ODE Part I: Basic Part II: Advanced Lehrstuhl Informatik V Winter 2007/2008 Part I: Basic 1 Direction Fields 2 Euler s Method Outlines Part I: Basic Part II: Advanced 3 Discretized
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationMinghao Wu Rostami Address: Webpage: Positions Tenure-Track Assistant Professor Postdoctoral Researcher Lecturer Graduate Assistant
Minghao Wu Rostami Address: Department of Mathematics, 215 Carnegie Building, Syracuse University, Syracuse, NY 13244 Email: mwrostam@syr.edu Webpage: http://mwrostam.mysite.syr.edu/ Positions Tenure-Track
More informationTwo-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5
Two-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5 Jingwei Zhu May 14, 2014 Instructor: Surya Pratap Vanka 1 Project Description The objective of
More informationTermination criteria for inexact fixed point methods
Termination criteria for inexact fixed point methods Philipp Birken 1 October 1, 2013 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany Department of Mathematics/Computer
More informationAdaptive Time Space Discretization for Combustion Problems
Konrad-Zuse-Zentrum fu r Informationstechnik Berlin Takustraße 7 D-14195 Berlin-Dahlem Germany JENS LANG 1,BODO ERDMANN,RAINER ROITZSCH Adaptive Time Space Discretization for Combustion Problems 1 Talk
More informationPreconditioning for Nonsymmetry and Time-dependence
Preconditioning for Nonsymmetry and Time-dependence Andy Wathen Oxford University, UK joint work with Jen Pestana and Elle McDonald Jeju, Korea, 2015 p.1/24 Iterative methods For self-adjoint problems/symmetric
More informationLinear Multistep Methods I: Adams and BDF Methods
Linear Multistep Methods I: Adams and BDF Methods Varun Shankar January 1, 016 1 Introduction In our review of 5610 material, we have discussed polynomial interpolation and its application to generating
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9
Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers
More informationIs My CFD Mesh Adequate? A Quantitative Answer
Is My CFD Mesh Adequate? A Quantitative Answer Krzysztof J. Fidkowski Gas Dynamics Research Colloqium Aerospace Engineering Department University of Michigan January 26, 2011 K.J. Fidkowski (UM) GDRC 2011
More informationUNIVERSITY OF MARYLAND NAVIER-STOKES EQUATIONS CS-TR #4073 / UMIACS TR #99-66 OCTOBER 1999
UNIVERSITY OF MARYLAND INSTITUTE FOR ADVANCED COMPUTER STUDIES DEPARTMENT OF COMPUTER SCIENCE EFFICIENT PRECONDITIONING OF THE LINEARIZED NAVIER-STOKES EQUATIONS CS-TR #4073 / UMIACS TR #99-66 DAVID SILVESTER
More informationROBUST DESIGN OF A MECHANICAL LATCH
LA-UR-15-21296 Page 1 of 20 ROBUST DESIGN OF A MECHANICAL LATCH François Hemez Technical Staff Member, XTD-Division, Los Alamos National Laboratory Adjunct Professor, University of California San Diego
More informationImplicitly Defined High-Order Operator Splittings for Parabolic and Hyperbolic Variable-Coefficient PDE Using Modified Moments
Implicitly Defined High-Order Operator Splittings for Parabolic and Hyperbolic Variable-Coefficient PDE Using Modified Moments James V. Lambers September 24, 2008 Abstract This paper presents a reformulation
More informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
More informationPRECONDITIONING TECHNIQUES FOR NEWTON S METHOD FOR THE INCOMPRESSIBLE NAVIER STOKES EQUATIONS
BIT Numerical Mathematics 43: 961 974, 2003. 2003 Kluwer Academic Publishers. Printed in the Netherlands. 961 PRECONDITIONING TECHNIQUES FOR NEWTON S METHOD FOR THE INCOMPRESSIBLE NAVIER STOKES EQUATIONS
More information7.4 The Saddle Point Stokes Problem
346 CHAPTER 7. APPLIED FOURIER ANALYSIS 7.4 The Saddle Point Stokes Problem So far the matrix C has been diagonal no trouble to invert. This section jumps to a fluid flow problem that is still linear (simpler
More informationTime-dependent Dirichlet Boundary Conditions in Finite Element Discretizations
Time-dependent Dirichlet Boundary Conditions in Finite Element Discretizations Peter Benner and Jan Heiland November 5, 2015 Seminar Talk at Uni Konstanz Introduction Motivation A controlled physical processes
More informationIncompressible Flow & Iterative Solver Software (IFISS) Installation & Software Guide 1
Incompressible Flow & Iterative Solver Software (IFISS) Installation & Software Guide 1 David J Silvester 2 Howard C Elman 3 Alison Ramage 4 Version 35 released 22 September 2016 1 This project was supported
More informationTurbulent Rankine Vortices
Turbulent Rankine Vortices Roger Kingdon April 2008 Turbulent Rankine Vortices Overview of key results in the theory of turbulence Motivation for a fresh perspective on turbulence The Rankine vortex CFD
More informationOn solving linear systems arising from Shishkin mesh discretizations
On solving linear systems arising from Shishkin mesh discretizations Petr Tichý Faculty of Mathematics and Physics, Charles University joint work with Carlos Echeverría, Jörg Liesen, and Daniel Szyld October
More informationDirect Numerical Simulation of fractal-generated turbulence
Direct Numerical Simulation of fractal-generated turbulence S. Laizet and J.C. Vassilicos Turbulence, Mixing and Flow Control Group, Department of Aeronautics and Institute for Mathematical Sciences, Imperial
More informationA STUDY OF MULTIGRID SMOOTHERS USED IN COMPRESSIBLE CFD BASED ON THE CONVECTION DIFFUSION EQUATION
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More informationStagnation-Point Pressure Distribution and Wall Shear Stress: Numerical Simulation and Similarity Solution
Stagnation-Point Pressure Distribution and Wall Shear Stress: Numerical Simulation and Similarity Solution Vai Kuong Sin and Tai Yin Tong Abstract Two-dimensional stagnation-point flow in three different
More informationElliptic Problems / Multigrid. PHY 604: Computational Methods for Physics and Astrophysics II
Elliptic Problems / Multigrid Summary of Hyperbolic PDEs We looked at a simple linear and a nonlinear scalar hyperbolic PDE There is a speed associated with the change of the solution Explicit methods
More informationPreconditioners for the incompressible Navier Stokes equations
Preconditioners for the incompressible Navier Stokes equations C. Vuik M. ur Rehman A. Segal Delft Institute of Applied Mathematics, TU Delft, The Netherlands SIAM Conference on Computational Science and
More informationA Space-Time Multigrid Solver Methodology for the Optimal Control of Time-Dependent Fluid Flow
A Space-Time Multigrid Solver Methodology for the Optimal Control of Time-Dependent Fluid Flow Michael Köster, Michael Hinze, Stefan Turek Michael Köster Institute for Applied Mathematics TU Dortmund Trier,
More informationLinear Solvers. Andrew Hazel
Linear Solvers Andrew Hazel Introduction Thus far we have talked about the formulation and discretisation of physical problems...... and stopped when we got to a discrete linear system of equations. Introduction
More informationProject 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions
Project 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions R. S. Sellers MAE 5440, Computational Fluid Dynamics Utah State University, Department of Mechanical and Aerospace Engineering
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 informationD-,+O K.1K. b9 V 378
CLSSIICATON ALD -A207 000 REPORT DOCUMENTATION PAGE ors Appove1 la. REPORT SECURITY CLASSIFICATION 1lb. RESTRICTIVE MARKINGS 2a. SECURITY CLASSIFICATION AUT _"ALECTE 3. DISTRIBUTION /AVAILABILITY OF REPOR
More informationFirst, Second, and Third Order Finite-Volume Schemes for Diffusion
First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.
More informationCode Verification of Multiphase Flow with MFIX
Code Verification of Multiphase Flow with MFIX A N I R U D D H A C H O U D H A R Y ( A N I R U D D @ V T. E D U ), C H R I S T O P H E R J. ROY ( C J R O Y @ V T. E D U ), ( V I R G I N I A T E C H ) 23
More informationUnsteady Incompressible Flow Simulation Using Galerkin Finite Elements with Spatial/Temporal Adaptation
Unsteady Incompressible Flow Simulation Using Galerkin Finite Elements with Spatial/Temporal Adaptation Mohamed S. Ebeida Carnegie Mellon University, Pittsburgh, PA 15213 Roger L. Davis and Roland W. Freund
More informationTime stepping methods
Time stepping methods ATHENS course: Introduction into Finite Elements Delft Institute of Applied Mathematics, TU Delft Matthias Möller (m.moller@tudelft.nl) 19 November 2014 M. Möller (DIAM@TUDelft) Time
More informationNumAn2014 Conference Proceedings
OpenAccess Proceedings of the 6th International Conference on Numerical Analysis, pp 198-03 Contents lists available at AMCL s Digital Library. NumAn014 Conference Proceedings Digital Library Triton :
More information20. A Dual-Primal FETI Method for solving Stokes/Navier-Stokes Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 23 DDM.org 2. A Dual-Primal FEI Method for solving Stokes/Navier-Stokes
More informationHigh-order ADI schemes for convection-diffusion equations with mixed derivative terms
High-order ADI schemes for convection-diffusion equations with mixed derivative terms B. Düring, M. Fournié and A. Rigal Abstract We consider new high-order Alternating Direction Implicit ADI) schemes
More informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationCHAPTER 5: Linear Multistep Methods
CHAPTER 5: Linear Multistep Methods Multistep: use information from many steps Higher order possible with fewer function evaluations than with RK. Convenient error estimates. Changing stepsize or order
More informationNewton s Method and Efficient, Robust Variants
Newton s Method and Efficient, Robust Variants Philipp Birken University of Kassel (SFB/TRR 30) Soon: University of Lund October 7th 2013 Efficient solution of large systems of non-linear PDEs in science
More informationJournal of Computational and Applied Mathematics. Multigrid method for solving convection-diffusion problems with dominant convection
Journal of Computational and Applied Mathematics 226 (2009) 77 83 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
More informationFREE BOUNDARY PROBLEMS IN FLUID MECHANICS
FREE BOUNDARY PROBLEMS IN FLUID MECHANICS ANA MARIA SOANE AND ROUBEN ROSTAMIAN We consider a class of free boundary problems governed by the incompressible Navier-Stokes equations. Our objective is to
More informationDan s Morris s Notes on Stable Fluids (Jos Stam, SIGGRAPH 1999)
Dan s Morris s Notes on Stable Fluids (Jos Stam, SIGGRAPH 1999) This is intended to be a detailed by fairly-low-math explanation of Stam s Stable Fluids, one of the key papers in a recent series of advances
More informationConvergence Acceleration for Flow Problems
Comprehensive Summaries of Uppsala Dissertations from the Faculty of Science and Technology 598 Convergence Acceleration for Flow Problems BY HENRIK BRANDÉN ACTA UNIVERSITATIS UPSALIENSIS UPPSALA 2001
More informationNonparametric density estimation for elliptic problems with random perturbations
Nonparametric density estimation for elliptic problems with random perturbations, DqF Workshop, Stockholm, Sweden, 28--2 p. /2 Nonparametric density estimation for elliptic problems with random perturbations
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More informationEfficient solution of stationary Euler flows with critical points and shocks
Efficient solution of stationary Euler flows with critical points and shocks Hans De Sterck Department of Applied Mathematics University of Waterloo 1. Introduction consider stationary solutions of hyperbolic
More informationApplication of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations
Application of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, CA, 94305
More informationStability of the Parareal Algorithm
Stability of the Parareal Algorithm Gunnar Andreas Staff and Einar M. Rønquist Norwegian University of Science and Technology Department of Mathematical Sciences Summary. We discuss the stability of the
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationNUMERICAL COMPUTATION OF A VISCOUS FLOW AROUND A CIRCULAR CYLINDER ON A CARTESIAN GRID
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2000 Barcelona, 11-14 September 2000 c ECCOMAS NUMERICAL COMPUTATION OF A VISCOUS FLOW AROUND A CIRCULAR CYLINDER
More informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationAssignment on iterative solution methods and preconditioning
Division of Scientific Computing, Department of Information Technology, Uppsala University Numerical Linear Algebra October-November, 2018 Assignment on iterative solution methods and preconditioning 1.
More information