Runge Kutta Methods Optimized For Advection Problems
|
|
- Holly Carter
- 5 years ago
- Views:
Transcription
1 Runge Kutta Methods Optimized For Advection Problems Brown University December 9, 2009
2 Introduction Classical Runge Kutta Methods Overview Standard Runge Kutta Schemes for Advection Equation Low Dissipation and Low Dispersion Runge Kutta() Schemes
3 Introduction Classical Runge Kutta Methods Introduction Computational acoustics is a recently emerging tool for acoustic problems Minimal dispersion and dissipation errors are desired and so high order spatial discretization schemes are used Need for accurate and efficient time advancing schemes
4 Introduction Classical Runge Kutta Methods 4th Order Runge Kutta Method In 4th order Runge Kutta method, we have for stability k t 2.83 Consider the Wave Equation u t + u x = 0 For k t = 1.6, method is stable but exhibit serious dispersion and dissipation errors dt=1.6; dx=
5 Introduction Classical Runge Kutta Methods Amplification Ũ n+1 k = Ũ n k (1 + p j=1 c j(ick t) j ) where c j are related to w i, β ij in U n+1 = u n + p i w i K i K i = tf (u n + i 1 j=1 β ijk j )i = 1, 2...p Amplification Factor r = Ũk n+1 Ũ k n = 1 + j=p j=1 c j( iσ) j where σ = ck t The exact amplification factor is r e = e ick t = e iσ
6 Introduction Alternating Schemes Low Storage Methods Low Dissipation and Dispersion Runge Kutta Methods Traditional choice of c j maximizes the possible order of accuracy Optimization is carried out instead on Γ p j=1 c j( iσ) j e iσ 2 dσ = MIN where Γ specifies the range of ck t in the optimization
7 Introduction Alternating Schemes Low Storage Methods Two Step Alternating Schemes Dissipation and dispersion errors can be further reduced and higher order of acccuracy can be obtained using alternate stepping schmes Amplification factors of the first and second steps are given by r 1 = 1 + p 1 j=1 a j( iσ) j r 2 = 1 + p 2 j=1 b j( iσ) j p 1, p 2 are the number of stages of the two steps. Combined Amplification factor = r 1 r 2 Optimization criteria Γ 0 (1 + p 1 j=1 a j( iσ) j )(1 + p 2 j=1 b j( iσ) j ) e 2iσ 2 dσ = MIN
8 Introduction Alternating Schemes Low Storage Methods Low Storage Methods 3N storage required for classical RK4 and covered so far. Consider du dt = F (t, U(t)); U(t 0 ) = U 0 For low storage schemes, Runge Kutta algorithm becomes w i = α i w i 1 + hf (t i 1, U i 1 ) U i = U i 1 + β i w i α 1 = 0; u 0 = u n 1 ; u n = u s ; t i = t n 1 + hc i
9 Example, Dissipation and Dispersion can be estimated by the amplification factor G(ω t) = ũn+1 (ω t) ũ n (ω t) = 1 + s j=1 γ j(iω t) j Dissipation = 1 G(ω t) Dispersion = ω t ω t /π
10 Example, Dissipation and Dispersion Limits Scheme Dissipation(10 4 ) Dispersion(10 4 ) RK RK46(S) RK26(B) RK46-NL
11 Numerical Example Example Linear Convective Equation with the methods that have been optimized is considered u t + u x2 x = 0 u(x, 0) = K exp( ) c 2 Analytic Solution u(x, t) = K exp( (x t)2 ) K = 0.5; c = 3; c 2 9 Point central difference of 8th order is used for spatial discretization u x (x j, t) 1 t [ 4 5 (u j 1 u j+1 ) (u j 2 u j+2 ) (u j 3 u j+3 ) (u j 4 u j+4 )] Domain x=-50 to x = 450 and dx = 1 Time interval = [0,400] and dt = 1.26
12 Absolute Errors Example x t LDD46(C) LDD25 RK4 1/ x x x / x x x / x x x / x x x 10 3
13 Example 4th Order Runge Kutta Method Classical RK4 scheme stable but not time accurate Time steps much smaller than stability limit required for long time integrations dt=1.26; dx=
14 5 Stage - Stanescu Example i α β c
15 LDD56-Stanescu Example Two step alternating scheme with five/six stages in the first/second step The scheme has fourth order accuracy in both steps
16 Further Research Advantages And Disadvantages Conclusion References Further Research Few more methods have been proposed to optimize Runge Kutta schemes for wave propogation phenomena Calvo et al - First maximize the stability range of the algorithm and then improve the range of well resolved Courant numbers Bogey and Bailly: Minimize the sum of the norms of dissipation and dispersion errors in a given range of frequencies Ramboer et al: Minimize the total dissipation and dispersion errors deriving from coupling with time integration
17 Further Research Advantages And Disadvantages Conclusion References Advantages/Disadvantages Advantages Minimize dispersion and dissipation Low storage requirements(only two memory slots per variable) Larger time steps can be used Disadvantages More computations/function evaluations as compared to Runge Kutta Fourth order accuracy ensured only for linear operators with large stability regions Weak stablility properties for non linear operators for schemes of order four. Does not strictly apply to problems where grid spacing is dictated by physical constraints.
18 Further Research Advantages And Disadvantages Conclusion References Analysis of the dissipation and dispersion properties of Runge Kutta methods. schemes are proposed for convection problems Efficient implementation of large scale wave propogation problems using 2N storage formulation.
19 Further Research Advantages And Disadvantages Conclusion References References Low Dissipation and Low Dispersion RUnge Kutta Schemes, Hu et al 2N Storage Low Dissipation and Dispersion Runge Kutta Schemes, Stanescu and Habashi A General Strategy for the Optimization of RUnge Kutta Schemes for Wave Propogation Phenomena, Bernardini et al Low Dissipation and Low Dispersion Fourth Order Runge Kutta Algorithm, Berland et al A New Minimum Storage RungeKutta Scheme For Computational Acoustics, Calvo et al
Turbulence et Génération de Bruit Equipe de recherche du Centre Acoustique LMFA, UMR CNRS 5509, Ecole Centrale de Lyon Simulation Numérique en Aéroacoustique Institut Henri Poincaré - 16 novembre 2006
More informationCS205b/CME306. Lecture 9
CS205b/CME306 Lecture 9 1 Convection Supplementary Reading: Osher and Fedkiw, Sections 3.3 and 3.5; Leveque, Sections 6.7, 8.3, 10.2, 10.4. For a reference on Newton polynomial interpolation via divided
More informationREVIEW OF NUMERICAL SCHEMES AND BOUNDARY CONDITIONS APPLIED TO WAVE PROPAGATION PROBLEMS
REVIEW OF NUMERICAL SCHEMES AND BOUNDARY CONDITIONS APPLIED TO WAVE PROPAGATION PROBLEMS O. de Almeida Universidade Federal de Uberlândia Departamento de Engenharia Mecânica FEMEC Campus Santa Mônica CP.
More informationMathematics for chemical engineers
Mathematics for chemical engineers Drahoslava Janovská Department of mathematics Winter semester 2013-2014 Numerical solution of ordinary differential equations Initial value problem Outline 1 Introduction
More informationOptimized Low Dispersion and Low Dissipation Runge- Kutta Algorithms in Computational Aeroacoustics
Appl. Math. Inf. Sci. 8, No. 1, 57-68 214 57 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/1.12785/amis/816 Optimized Low Dispersion Low Dissipation Runge- Kutta
More informationA-posteriori Diffusion Analysis of Numerical Schemes in Wavenumber Domain
2th Annual CFD Symposium, August 9-1, 218, Bangalore A-posteriori Diffusion Analysis of Numerical Schemes in Wavenumber Domain S. M. Joshi & A. Chatterjee Department of Aerospace Engineering Indian Institute
More informationNumerical Methods on the Image Processing Problems
Numerical Methods on the Image Processing Problems Department of Mathematics and Statistics Mississippi State University December 13, 2006 Objective Develop efficient PDE (partial differential equations)
More informationOn the simulation of wave propagation with a higher-order finite volume scheme based on Reproducing Kernel Methods
On the simulation of wave propagation with a higher-order finite volume scheme based on Reproducing Kernel Methods X. Nogueira a, I. Colominas,a, L. Cueto-Felgueroso b, S. Khelladi c a Group of Numerical
More informationParallel Adaptive Tsunami Modelling with Triangular Discontinuous Galerkin Schemes
Parallel Adaptive Tsunami Modelling with Triangular Discontinuous Galerkin Schemes Stefan Vater 1 Kaveh Rahnema 2 Jörn Behrens 1 Michael Bader 2 1 Universität Hamburg 2014 PDES Workshop 2 TU München Partial
More informationHigh Order Weighted Essentially Non-Oscillatory Schemes for Convection. Dominated Problems. Chi-Wang Shu 1
High Order Weighted Essentially Non-Oscillatory Schemes for Convection Dominated Problems Chi-Wang Shu Division of Applied Mathematics, Brown University, Providence, Rhode Island 09 ABSTRACT High order
More informationChE 400: Applied Chemical Engineering Calculations Tutorial 6: Numerical Solution of ODE Using Excel and Matlab
ChE 400: Applied Chemical Engineering Calculations Tutorial 6: Numerical Solution of ODE Using Excel and Matlab Tutorial 6: Numerical Solution of ODE Gerardine G. Botte This handout contains information
More informationChapter 6. Semi-Lagrangian Methods
Chapter 6. Semi-Lagrangian Methods References: Durran Chapter 6. Review article by Staniford and Cote (1991) MWR, 119, 2206-2223. 6.1. Introduction Semi-Lagrangian (S-L for short) methods, also called
More informationSemi-Conservative Schemes for Conservation Laws
Semi-Conservative Schemes for Conservation Laws Rosa Maria Pidatella 1 Gabriella Puppo, 2 Giovanni Russo, 1 Pietro Santagati 1 1 Università degli Studi di Catania, Catania, Italy 2 Università dell Insubria,
More informationFinite difference methods
Finite difference methods Siltanen/Railo/Kaarnioja Spring 8 Applications of matrix computations Applications of matrix computations Finite difference methods Spring 8 / Introduction Finite difference methods
More informationLevel Set Methods and Fast Marching Methods
Level Set Methods and Fast Marching Methods I.Lyulina Scientific Computing Group May, 2002 Overview Existing Techniques for Tracking Interfaces Basic Ideas of Level Set Method and Fast Marching Method
More informationMid-Year Report. Discontinuous Galerkin Euler Equation Solver. Friday, December 14, Andrey Andreyev. Advisor: Dr.
Mid-Year Report Discontinuous Galerkin Euler Equation Solver Friday, December 14, 2012 Andrey Andreyev Advisor: Dr. James Baeder Abstract: The focus of this effort is to produce a two dimensional inviscid,
More informationEuler s Methods (a family of Runge- Ku9a methods)
Euler s Methods (a family of Runge- Ku9a methods) ODE IVP An Ordinary Differential Equation (ODE) is an equation that contains a function having one independent variable: The equation is coupled with an
More informationODEs occur quite often in physics and astrophysics: Wave Equation in 1-D stellar structure equations hydrostatic equation in atmospheres orbits
Solving ODEs General Stuff ODEs occur quite often in physics and astrophysics: Wave Equation in 1-D stellar structure equations hydrostatic equation in atmospheres orbits need workhorse solvers to deal
More informationFault Tolerant Domain Decomposition for Parabolic Problems
Fault Tolerant Domain Decomposition for Parabolic Problems Marc Garbey and Hatem Ltaief Department of Computer Science, University of Houston, Houston, TX 77204 USA garbey@cs.uh.edu, ltaief@cs.uh.edu 1
More informationA Diagonal Split-cell Model for the High-order Symplectic FDTD Scheme
PIERS ONLINE, VOL. 2, NO. 6, 2006 715 A Diagonal Split-cell Model for the High-order Symplectic FDTD Scheme Wei Sha, Xianliang Wu, and Mingsheng Chen Key Laboratory of Intelligent Computing & Signal Processing
More informationCONSERVATIVE AND NON-CONSERVATIVE METHODS BASED ON HERMITE WEIGHTED ESSENTIALLY-NON-OSCILLATORY RECONSTRUCTION FOR VLASOV EQUATIONS
CONSERVATIVE AND NON-CONSERVATIVE METHODS BASED ON HERMITE WEIGHTED ESSENTIALLY-NON-OSCILLATORY RECONSTRUCTION FOR VLASOV EQUATIONS CHANG YANG AND FRANCIS FILBET Abstract. We develop weighted essentially
More informationJoint Optimization of the Spatial and the Temporal Discretization Scheme for Accurate Computation of Acoustic Problems
Commun. Comput. Phys. doi:.428/cicp.oa-27-9 Vol., No., pp. -27 2 Joint Optimization of the Spatial and the Temporal Discretization Scheme for Accurate Computation of Acoustic Problems Jitenjaya Pradhan,
More informationNew formulations of the semi-lagrangian method for Vlasov-type equations
New formulations of the semi-lagrangian method for Vlasov-type equations Eric Sonnendrücker IRMA Université Louis Pasteur, Strasbourg projet CALVI INRIA Nancy Grand Est 17 September 2008 In collaboration
More informationarxiv: v1 [math.na] 26 Jun 2014
for spectrally accurate wave propagation Vladimir Druskin, Alexander V. Mamonov and Mikhail Zaslavsky, Schlumberger arxiv:406.6923v [math.na] 26 Jun 204 SUMMARY We develop a method for numerical time-domain
More informationIntegral Equation Methods for Vortex Dominated Flows, a High-order Conservative Eulerian Approach
Integral Equation Methods for Vortex Dominated Flows, a High-order Conservative Eulerian Approach J. Bevan, UIUC ICERM/HKUST Fast Integral Equation Methods January 5, 2016 Vorticity and Circulation Γ =
More informationLong time integrations of a convective PDE on the sphere by RBF collocation
Long time integrations of a convective PDE on the sphere by RBF collocation Bengt Fornberg and Natasha Flyer University of Colorado NCAR Department of Applied Mathematics Institute for Mathematics Applied
More informationThe Jello Cube Assignment 1, CSCI 520. Jernej Barbic, USC
The Jello Cube Assignment 1, CSCI 520 Jernej Barbic, USC 1 The jello cube Undeformed cube Deformed cube The jello cube is elastic, Can be bent, stretched, squeezed,, Without external forces, it eventually
More informationWeno Scheme for Transport Equation on Unstructured Grids with a DDFV Approach
Weno Scheme for Transport Equation on Unstructured Grids with a DDFV Approach Florence Hubert and Rémi Tesson Abstract In this paper we develop a DDFV approach for WENO scheme on unstructred grids for
More informationThird Order WENO Scheme on Three Dimensional Tetrahedral Meshes
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 5, No. 2-4, pp. 86-848 Commun. Comput. Phys. February 29 Third Order WENO Scheme on Three Dimensional Tetrahedral Meshes Yong-Tao Zhang 1, and Chi-Wang Shu
More informationComputer Project 3. AA Computational Fluid Dyanmics University of Washington. Mishaal Aleem March 17, 2015
Computer Project 3 AA 543 - Computational Fluid Dyanmics University of Washington Mishaal Aleem March 17, 2015 Contents Introduction........................................... 1 3.1 Grid Generator.......................................
More informationOn the Construction, Comparison, and Local Characteristic Decomposition for High-Order Central WENO Schemes
Journal of Computational Physics 8, 87 09 (00) doi:0.006/jcph.00.79 On the Construction, Comparison, and Local Characteristic Decomposition for High-Order Central WENO Schemes Jianxian Qiu, and Chi-Wang
More informationThe WENO Method in the Context of Earlier Methods To approximate, in a physically correct way, [3] the solution to a conservation law of the form u t
An implicit WENO scheme for steady-state computation of scalar hyperbolic equations Sigal Gottlieb Mathematics Department University of Massachusetts at Dartmouth 85 Old Westport Road North Dartmouth,
More informationFast Multipole Method on the GPU
Fast Multipole Method on the GPU with application to the Adaptive Vortex Method University of Bristol, Bristol, United Kingdom. 1 Introduction Particle methods Highly parallel Computational intensive Numerical
More informationStability Analysis of the Muscl Method on General Unstructured Grids for Applications to Compressible Fluid Flow
Stability Analysis of the Muscl Method on General Unstructured Grids for Applications to Compressible Fluid Flow F. Haider 1, B. Courbet 1, J.P. Croisille 2 1 Département de Simulation Numérique des Ecoulements
More informationHorizontal Mixing in the WRF-ARW Model. Russ Schumacher AT April 2006
Horizontal Mixing in the WRF-ARW Model Russ Schumacher AT 730 5 April 2006 Overview Description of implicit and explicit horizontal mixing in WRF Evaluation of WRF s implicit mixing Model options for explicit
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
AIAA SciTech Forum 8 January 08, Kissimmee, Florida 08 AIAA Aerospace Sciences Meeting 0.54/6.08-0067 A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley and Dr.
More informationSemi-Lagrangian Advection. The Basic Idea. x 1.
3.2.6. Semi-Lagrangian Advection We have studied the Eulerian leapfrog scheme and found it to be conditionally stable. The criterion for stability was the CFL condition µ c t x 1. For high spatial resolution
More informationHarmonic Spline Series Representation of Scaling Functions
Harmonic Spline Series Representation of Scaling Functions Thierry Blu and Michael Unser Biomedical Imaging Group, STI/BIO-E, BM 4.34 Swiss Federal Institute of Technology, Lausanne CH-5 Lausanne-EPFL,
More informationOperator Upscaling and Adjoint State Method
Operator Upscaling and Adjoint State Method Tetyana Vdovina, William Symes The Rice Inversion Project Rice University vdovina@rice.edu February 0, 009 Motivation Ultimate Goal: Use 3D elastic upscaling
More informationThe Nagumo Equation with Comsol Multiphysics
The Nagumo Equation with Comsol Multiphysics Denny Otten 1 Christian Döding 2 Department of Mathematics Bielefeld University 33501 Bielefeld Germany Date: 25. April 2016 1. Traveling Front in the Nagumo
More informationDownloaded 10/27/17 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 37, No. 3, pp. C354 C383 c 2015 Society for Industrial and Applied Mathematics EFFICIENT IMPLEMENTATION OF NONLINEAR COMPACT SCHEMES ON MASSIVELY PARALLEL PLATFORMS DEBOJYOTI
More informationarxiv:math/ v1 [math.ap] 18 Jul 2006
arxiv:math/6747v [math.ap] 8 Jul 6 A linear dispersive mechanism for numerical error growth: spurious caustics Claire David, Pierre Sagaut, Tapan Sengupta Université Pierre et Marie Curie-Paris 6 Laboratoire
More informationCOMPARISON OF HIGH-ACCURACY FINITE-DIFFERENCE METHODS FOR LINEAR WAVE PROPAGATION
SIAM J. SCI. COMPUT. Vol. 22, No. 2, pp. 476 52 c 2 Society for Industrial and Applied Mathematics COMPARISON OF HIGH-ACCURACY FINITE-DIFFERENCE METHODS FOR LINEAR WAVE PROPAGATION DAVID W. ZINGG Abstract.
More informationObservations on the fifth-order WENO method with non-uniform meshes
Observations on the fifth-order WENO method with non-uniform meshes Rong Wang, a, Hui Feng, b, Raymond J. Spiteri a,, a Department of Computer Science, University of Saskatchewan, Saskatoon, SK, S7N 5C9,
More informationHigh Accuracy Schemes for DNS and Acoustics
Journal of Scientific Computing ( 6) DOI:.7/s95-5-498- High Accuracy Schemes for DNS and Acoustics T. K. Sengupta, S. K. Sircar, and A. Dipankar Received May, 4; accepted (in revised form) January 7, 5
More informationConservative high order semi-lagrangian finite difference WENO methods for advection in incompressible flow. Abstract
Conservative high order semi-lagrangian finite difference WENO methods for advection in incompressible flow Jing-Mei Qiu 1 and Chi-Wang Shu Abstract In this paper, we propose a semi-lagrangian finite difference
More informationReal time Ray-Casting of Algebraic Surfaces
Real time Ray-Casting of Algebraic Surfaces Martin Reimers Johan Seland Center of Mathematics for Applications University of Oslo Workshop on Computational Method for Algebraic Spline Surfaces Thursday
More informationAudio acoustic modeling using full-wave methods
Acoustics 8 Paris Audio acoustic modeling using full-wave methods Timo Lahivaara, Tomi Huttunen and Simo-Pekka Simonaho University of Kuopio, P.O.Box 1627, 7211 Kuopio, Finland simo-pekka.simonaho@uku.fi
More informationCUDA GPGPU. Ivo Ihrke Tobias Ritschel Mario Fritz
CUDA GPGPU Ivo Ihrke Tobias Ritschel Mario Fritz Today 3 Topics Intro to CUDA (NVIDIA slides) How is the GPU hardware organized? What is the programming model? Simple Kernels Hello World Gaussian Filtering
More informationComputing Nearly Singular Solutions Using Pseudo-Spectral Methods
Computing Nearly Singular Solutions Using Pseudo-Spectral Methods Thomas Y. Hou Ruo Li January 9, 2007 Abstract In this paper, we investigate the performance of pseudo-spectral methods in computing nearly
More informationSimulation in Computer Graphics. Particles. Matthias Teschner. Computer Science Department University of Freiburg
Simulation in Computer Graphics Particles Matthias Teschner Computer Science Department University of Freiburg Outline introduction particle motion finite differences system of first order ODEs second
More informationNew Very High-Order Upwind Multilayer Compact Schemes with Spectral-Like Resolution for Flow Simulations
New Very High-Order Upwind Multilayer Compact Schemes with Spectral-Lie Resolution for Flow Simulations Zeyu Bai and Xiaolin Zhong University of California, Los Angeles, CA, 995, USA Hypersonic boundary
More informationA Fast Phase Space Method for Computing Creeping Rays KTH. High Frequency Wave Propagation CSCAMM, September 19-22, 2005
A Fast Phase Space Method for Computing Creeping Rays Olof Runborg Mohammad Motamed KTH High Frequency Wave Propagation CSCAMM, September 19-22, 2005 Geometrical optics for high-frequency waves Consider
More informationNumerical Analysis of Shock Tube Problem by using TVD and ACM Schemes
Numerical Analysis of Shock Tube Problem by using TVD and Schemes Dr. Mukkarum Husain, Dr. M. Nauman Qureshi, Syed Zaid Hasany IST Karachi, Email: mrmukkarum@yahoo.com Abstract Computational Fluid Dynamics
More informationHigh Order Fixed-Point Sweeping WENO Methods for Steady State of Hyperbolic Conservation Laws and Its Convergence Study
Commun. Comput. Phys. doi:.48/cicp.375.6a Vol., No. 4, pp. 835-869 October 6 High Order Fixed-Point Sweeping WENO Methods for Steady State of Hyperbolic Conservation Laws and Its Convergence Study Liang
More informationVon Neumann Analysis for Higher Order Methods
1. Introduction Von Neumann Analysis for Higher Order Methods Von Neumann analysis is a widely used method to study how an initial wave is propagated with certain numerical schemes for a linear wave equation
More informationSkåne University Hospital Lund, Lund, Sweden 2 Deparment of Numerical Analysis, Centre for Mathematical Sciences, Lund University, Lund, Sweden
Volume Tracking: A New Method for Visualization of Intracardiac Blood Flow from Three-Dimensional, Time-Resolved, Three-Component Magnetic Resonance Velocity Mapping Appendix: Theory and Numerical Implementation
More informationA Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere
A Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere Ram Nair Computational and Information Systems Laboratory (CISL) National Center for Atmospheric Research
More informationPost Processing, Visualization, and Sample Output
Chapter 7 Post Processing, Visualization, and Sample Output Upon successful execution of an ADCIRC run, a number of output files will be created. Specifically which files are created depends upon how the
More informationOn the order of accuracy and numerical performance of two classes of finite volume WENO schemes
On the order of accuracy and numerical performance of two classes of finite volume WENO schemes Rui Zhang, Mengping Zhang and Chi-Wang Shu November 29, 29 Abstract In this paper we consider two commonly
More informationFrom Hyperbolic Diffusion Scheme to Gradient Method: Implicit Green-Gauss Gradients for Unstructured Grids
Preprint accepted in Journal of Computational Physics. https://doi.org/10.1016/j.jcp.2018.06.019 From Hyperbolic Diffusion Scheme to Gradient Method: Implicit Green-Gauss Gradients for Unstructured Grids
More informationBACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION
BACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION TODD F. DUPONT AND YINGJIE LIU Abstract. We propose a method that significantly
More informationSimulation of Turbulent Particulate Flow on HPC Systems
Simulation of Turbulent Particulate Flow on HPC Systems Konstantin Fröhlich, Lennart Schneiders, Matthias Meinke, Wolfgang Schröder k.froehlich@aia.rwth-aachen.de Institute of Aerodynamics RWTH Aachen
More informationEFFICIENT SOLUTION ALGORITHMS FOR HIGH-ACCURACY CENTRAL DIFFERENCE CFD SCHEMES
EFFICIENT SOLUTION ALGORITHMS FOR HIGH-ACCURACY CENTRAL DIFFERENCE CFD SCHEMES B. Treidler, J.A. Ekaterineris and R.E. Childs Nielsen Engineering & Research, Inc. Mountain View, CA, 94043 Abstract Preliminary
More informationAn Eulerian Approach for Computing the Finite Time Lyapunov Exponent (FTLE)
An Eulerian Approach for Computing the Finite Time Lyapunov Exponent (FTLE) Shingyu Leung Department of Mathematics, Hong Kong University of Science and Technology masyleung@ust.hk May, Shingyu Leung (HKUST)
More informationHigh-Order CENO Reconstruction Scheme For Three-Dimensional Unstructured Mesh
High-Order CENO Reconstruction Scheme For Three-Dimensional Unstructured Mesh by Al-Amin Aziz A thesis submitted in conformity with the requirements for the degree of Masters of Applied Science Graduate
More informationCircuit Placement: 2000-Caldwell,Kahng,Markov; 2002-Kennings,Markov; 2006-Kennings,Vorwerk
Circuit Placement: 2000-Caldwell,Kahng,Markov; 2002-Kennings,Markov; 2006-Kennings,Vorwerk Andrew A. Kennings, Univ. of Waterloo, Canada, http://gibbon.uwaterloo.ca/ akenning/ Igor L. Markov, Univ. of
More informationAn Orthogonal Curvilinear Terrain-Following Coordinate for Atmospheric Models!
Solution of Partial Differential Equations on the Sphere! An Orthogonal Curvilinear Terrain-Following Coordinate for Atmospheric Models! Yiyuan Li 1, Bin Wang 1, Donghai Wang 2! 1. Institute of Atmospheric
More informationA TSTT integrated FronTier code and its applications in computational fluid physics
Institute of Physics Publishing Journal of Physics: Conference Series 16 (2005) 471 475 doi:10.1088/1742-6596/16/1/064 SciDAC 2005 A TSTT integrated FronTier code and its applications in computational
More informationFinal Report. Discontinuous Galerkin Compressible Euler Equation Solver. May 14, Andrey Andreyev. Adviser: Dr. James Baeder
Final Report Discontinuous Galerkin Compressible Euler Equation Solver May 14, 2013 Andrey Andreyev Adviser: Dr. James Baeder Abstract: In this work a Discontinuous Galerkin Method is developed for compressible
More informationA gradient-augmented level set method with an optimally local, coherent advection scheme
A gradient-augmented level set method with an optimally local, coherent advection scheme The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters.
More informationThe jello cube. Undeformed cube. Deformed cube
The Jello Cube Assignment 1, CSCI 520 Jernej Barbic, USC Undeformed cube The jello cube Deformed cube The jello cube is elastic, Can be bent, stretched, squeezed,, Without external forces, it eventually
More informationCombos: Advection-Diffusion and operator splitting
96 Chapter 7 Combos: Advection-Diffusion and operator splitting Up to this point we have treated advection and diffusion as two separate processes and considered appropriate numerical techniques for solving
More informationPhysics Tutorial 2: Numerical Integration Methods
Physics Tutorial 2: Numerical Integration Methods Summary The concept of numerically solving differential equations is explained, and applied to real time gaming simulation. Some objects are moved in a
More informationROC-HJ: Reachability analysis and Optimal Control problems - Hamilton-Jacobi equations
ROC-HJ: Reachability analysis and Optimal Control problems - Hamilton-Jacobi equations Olivier Bokanowski, Anna Désilles, Hasnaa Zidani February 2013 1 Introduction The software ROC-HJ is a C++ library
More informationEdge Detection Free Postprocessing for Pseudospectral Approximations
Edge Detection Free Postprocessing for Pseudospectral Approximations Scott A. Sarra March 4, 29 Abstract Pseudospectral Methods based on global polynomial approximation yield exponential accuracy when
More informationA new class of central compact schemes with spectral-like resolution II: Hybrid weighted nonlinear schemes. Abstract
A new class of central compact schemes with spectral-like resolution II: Hybrid weighted nonlinear schemes Xuliang Liu 1, Shuhai Zhang, Hanxin Zhang 3 and Chi-Wang Shu 4 Abstract In this paper, we develop
More informationTimo Lähivaara, Tomi Huttunen, Simo-Pekka Simonaho University of Kuopio, Department of Physics P.O.Box 1627, FI-70211, Finland
Timo Lähivaara, Tomi Huttunen, Simo-Pekka Simonaho University of Kuopio, Department of Physics P.O.Box 627, FI-72, Finland timo.lahivaara@uku.fi INTRODUCTION The modeling of the acoustic wave fields often
More informationMedical Image Segmentation using Level Sets
Medical Image Segmentation using Level Sets Technical Report #CS-8-1 Tenn Francis Chen Abstract Segmentation is a vital aspect of medical imaging. It aids in the visualization of medical data and diagnostics
More informationHomework Set #2-3, Math 475B
Homework Set #2-3, Math 475B Part I: Matlab In the last semester you learned a number of essential features of MATLAB. 1. In this instance, you will learn to make 3D plots and contour plots of z = f(x,
More informationImplicit schemes for wave models
Implicit schemes for wave models Mathieu Dutour Sikirić Rudjer Bo sković Institute, Croatia and Universität Rostock April 17, 2013 I. Wave models Stochastic wave modelling Oceanic models are using grids
More informationCollocation and optimization initialization
Boundary Elements and Other Mesh Reduction Methods XXXVII 55 Collocation and optimization initialization E. J. Kansa 1 & L. Ling 2 1 Convergent Solutions, USA 2 Hong Kong Baptist University, Hong Kong
More informationA Semi-Lagrangian High-Order Method for Navier Stokes Equations
Journal of Computational Physics 172, 658 684 (2001) doi:10.1006/jcph.2001.6847, available online at http://www.idealibrary.com on A Semi-Lagrangian High-Order Method for Navier Stokes Equations Dongbin
More informationVector Visualization
Vector Visualization Vector Visulization Divergence and Vorticity Vector Glyphs Vector Color Coding Displacement Plots Stream Objects Texture-Based Vector Visualization Simplified Representation of Vector
More informationThis is an author-deposited version published in: Eprints ID: 4362
This is an author-deposited version published in: http://oatao.univ-toulouse.fr/ Eprints ID: 4362 To cite this document: CHIKHAOUI Oussama, GRESSIER Jérémie, GRONDIN Gilles. Assessment of the Spectral
More informationPartition Design and Optimization for High-Order Spectral Volume Schemes on Tetrahedral Grids
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition 4-7 January 200, Orlando, Florida AIAA 200-720 Partition Design and Optimization for High-Order Spectral Volume
More informationHigh-Order Finite Difference Schemes for computational MHD
High-Order Finite Difference Schemes for computational MHD A. Mignone 1, P. Tzeferacos 1 and G. Bodo 2 [1] Dipartimento di Fisica Generale, Turin University, ITALY [2] INAF Astronomic Observatory of Turin,,
More informationOptimally Blended Finite/Spectral Element Scheme for Wave Propagation
Modern Techniques for Numerical Solution of PDEs, ACMAC, Crete, September 2011 p. 1/54 Optimally Blended Finite/Spectral Element Scheme for Wave Propagation Mark Ainsworth (joint work with Hafiz Abdul
More informationLocal Selection for Heuristic Algorithms as a Factor in Accelerating Optimum Search
Local Selection for Heuristic Algorithms as a Factor in Accelerating Optimum Search Danuta Jama Institute of Mathematics Silesian University of Technology Kaszubska 23, 44-100 Gliwice, Poland Email: Danuta.Jama@polsl.pl
More information37 Self-Assessment. 38 Two-stage Runge-Kutta Methods. Chapter 10 Runge Kutta Methods
Chapter 10 Runge Kutta Methods In the previous lectures, we have concentrated on multi-step methods. However, another powerful set of methods are known as multi-stage methods. Perhaps the best known of
More informationRK-Time Integration and High-Order Spatial Discretization a New Dynamical Core for the LMK
RK-Time Integration and High-Order Spatial Discretization a New Dynamical Core for the LMK Jochen Förstner Deutscher Wetterdienst, Kaiserleistraße 35, 63067 Offenbach am Main jochen.foerstner@dwd.de COSMO
More informationComputational Methods. Sources of Errors
Computational Methods Sources of Errors Manfred Huber 2011 1 Numerical Analysis / Scientific Computing Many problems in Science and Engineering can not be solved analytically on a computer Numeric solutions
More informationAcoustic Simulation. COMP 768 Presentation Lakulish Antani April 9, 2009
Acoustic Simulation COMP 768 Presentation Lakulish Antani April 9, 2009 Acoustic Simulation Sound Synthesis Sound Propagation Sound Rendering 2 Goal Simulate the propagation of sound in an environment
More informationSpace Filling Curves and Hierarchical Basis. Klaus Speer
Space Filling Curves and Hierarchical Basis Klaus Speer Abstract Real world phenomena can be best described using differential equations. After linearisation we have to deal with huge linear systems of
More informationGeodetic Iterative Methods for Nonlinear Acoustic Sources Localization: Application to Cracks Eects Detected with Acoustic Emission
6th NDT in Progress 2011 International Workshop of NDT Experts, Prague, 10-12 Oct 2011 Geodetic Iterative Methods for Nonlinear Acoustic Sources Localization: Application to Cracks Eects Detected with
More informationSOLVING PARTIAL DIFFERENTIAL EQUATIONS ON POINT CLOUDS
SOLVING PARTIAL DIFFERENTIAL EQUATIONS ON POINT CLOUDS JIAN LIANG AND HONGKAI ZHAO Abstract. In this paper we present a general framework for solving partial differential equations on manifolds represented
More informationM.Sc. in Computational Science. Fundamentals of Atmospheric Modelling
M.Sc. in Computational Science Fundamentals of Atmospheric Modelling Peter Lynch, Met Éireann Mathematical Computation Laboratory (Opp. Room 30) Dept. of Maths. Physics, UCD, Belfield. January April, 2004.
More informationModule1: Numerical Solution of Ordinary Differential Equations. Lecture 6. Higher order Runge Kutta Methods
Module1: Numerical Solution of Ordinary Differential Equations Lecture 6 Higher order Runge Kutta Methods Keywords: higher order methods, functional evaluations, accuracy Higher order Runge Kutta Methods
More information21. Efficient and fast numerical methods to compute fluid flows in the geophysical β plane
12th International Conference on Domain Decomposition Methods Editors: Tony Chan, Takashi Kako, Hideo Kawarada, Olivier Pironneau, c 2001 DDM.org 21. Efficient and fast numerical methods to compute fluid
More informationHigh-Order Navier-Stokes Simulations using a Sparse Line-Based Discontinuous Galerkin Method
High-Order Navier-Stokes Simulations using a Sparse Line-Based Discontinuous Galerkin Method Per-Olof Persson University of California, Berkeley, Berkeley, CA 9472-384, U.S.A. We study some of the properties
More information