Parallel Adaptive Tsunami Modelling with Triangular Discontinuous Galerkin Schemes
|
|
- Annice Skinner
- 5 years ago
- Views:
Transcription
1 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 Differential Equations on the Sphere
2 Tsunami Modeling within ASCETE Advanced Simulation of Coupled Earthquake and Tsunami Events Geosciences Rupture Modeling Seismology ASCETE Mathema3cs Computer Science Tsunami Modeling Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
3 Tsunami Modeling within ASCETE Advanced Simulation of Coupled Earthquake and Tsunami Events Geosciences Rupture Modeling Seismology ASCETE Mathema3cs Computer Science Tsunami Modeling Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
4 Tsunami Modeling within ASCETE Advanced Simulation of Coupled Earthquake and Tsunami Events Geosciences Rupture Modeling Seismology ASCETE Mathema3cs Computer Science Tsunami Modeling Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
5 Numerical Simulation of Tsunami Events Some Requirements A numerical algorithm for tsunami simulation should... accurately represent wave propagation and indundation at the coast, be well-balanced (preserves the still water steady state), be mass conservative, be robust (non-oscillatory, positive,... ), high order of minor priority, be computationally efficient (small stencil,... ). Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
6 The Governing Equations Shallow Water Equations system of conservation laws with source term (balance laws): U t + F(U) = S(U) in Ω T R2 R + with U = (φ, φu) T, φ = gh, subject to initial and boundary conditions flux and source defined as ( ) φu F(U) = φu u φ2 I ( ) 0 S(U) = S b + S r + S f + S v z H(t, x) h(t,x) = H b S b = φg b S r = f(k φu) bathymetry Coriolis force b(t, x) x S f = γ 2 u u /φ 1/3 S v = ν (φu) bottom friction eddy viscosity Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
7 Space Discretization DG Formulation & Basis Function Representation Galerkin ansatz (residual orthogonal to some function space V ): ( ) U + F(U) S(U) ψ dx = 0 ψ V t Ω weak DG formulation (integration by parts): ( ) Uh F h S h ψ ω dx = ψ ω F h n dσ t T ω T ω Representation on reference element: nodal approach based on electrostatics / Fekete points local solution represented by Lagrange polynomials ψ i V h of degree p U(ξ, t) = r U(ξ i, t)ψ i (ξ) = i=1 r Ũ i (t)ψ i (ξ) i=1 following GIRALDO ET AL. [2002], GIRALDO [2006] Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
8 Time discretization Runge-Kutta Timestepping with Limiting SSP Multistage scheme: U (0) h U (i) h U n+1 h = Un h = ΛΠ h = U (k+1) h { i 1 l=0 α il U (l) h + β il t n H h (U (l) h ) } for i = 1... k + 1 Intermediate solution limited with ΛΠ h after each stage: TVB corrected minmod function [COCKBURN AND SHU, 1998]: m(u i,x, ū i+1 ū i, ū i ū i 1 ) limiting in conservative (h, hu) vs. hydrostatic variables (h + b, hu) positivity preservation [XING ET AL., 2010]: Ũ n (x, y) = θ ( ( ) U n (x, y) U n) + U n φ with θ = min 1, n φ n min x {φ n (x)} Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
9 2D Test Case Quasi Stationary Vortex Advected by Constant Background Flow using grid library amatos (conforming elements, refinement by bisection) linear elements, no limiting simple gradient based error indicator Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
10 1D Results Oscillatory Flow in a Parabolic Bowl domain [ 5000, 5000] with b(x) = 10(x/3000) cells, t = 1.0, CFL 0.2 limiting in (h, hu) analytical solution by Thacker (1981): h(x, t) + b(x) = h 0 B2 Bx 8h 0 (1 + cos(2ωt)) cos ωt, 4g 2a g ω = 2gh 0 /a, B = 5 initial momentum at t = 0 is set to zero Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
11 1D Results Oscillatory Flow in a Parabolic Bowl t = 0 t = 1000 t = 2000 t = 3000 Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
12 1D Results Tsunami Runup onto a Sloping Beach Benchmark problem 1 from The Third International Workshop on Long-Wave Runup Models (2004). uniformly sloping beach with b(x) = x prescribed initial surface elevation and momentum with (hu)(x, 0) 0 domain [ 500; 50000], 1010 cells t = 0.05, CFL 0.2 limiting in (h, hu) Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
13 1D Results Tsunami Runup onto a Sloping Beach t = 0 t = 160 t = 175 t = 220 Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
14 Wellbalancing with Wetting and Drying How to Discretize a Semi-dry Cell? Assuming a fixed spatial grid (no moving grid points), piecewise linear DG discretization:? Tsunami Runup at t = 50 limiting in (h, hu) Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
15 Wetting & Drying Solution Approach Well-balancing with wetting/drying: distinguish between flood type and dam-break type cells neglect gradient in surface elevation in flood-type cells inner cell redistribution of tendencies to wet nodes in flood-type cells limiting only in fully wet cells, adjustment of negative heights flood type dam-break type Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
16 1D Results Lake at Rest domain [0, 1] with periodic boundary conditions still water u 0 50 cells, t = 0.002, CFL 0.3, t max = 20 limiting in (h + b, hu) Test initialized with zero momentum field:... with random deviations in the momentum field of the order 10 8 : Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
17 Parallel AMR for Triangular Grids Challenges for Grid Refinement and Coarsening strong, dynamically adaptive refinement (to capture solution and geometrical details etc.) frequent re-meshing of large parts of the grid substantial change of problem size during simulation Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
18 Parallel AMR for Triangular Grids Challenges for Parallelisation trend to multi-/manycore; multiple layers of parallelism dynamic load balancing due to remeshing and varying computational load goal: retain locality properties (partitioning) Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
19 Structured Triangular Meshes Newest Vertex Bisection Refinement & Sierpinski Curves recursively structured triangular grids (newest vertex bisection) fully adaptive grid described by a corresponding refinement tree element orders (tree and grid cells) defined by Sierpinski space-filling curves minimum memory requirements triangle strips as data structure exploit locality properties for cache and parallel efficiency Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
20 Test Scenario Adaptive Tsunami Simulation with Finite Volumes basic shallow water model, Finite Volume discretisation augmented Riemann solver (George, 08) provided by GeoClaw dynamically adaptive grid with up to 100 Mio grid cells refinement/coarsening and load balancing after each time step Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
21 Test Scenario Adaptive Tsunami Simulation with Finite Volumes 1 Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
22 Test Scenario Adaptive Tsunami Simulation with Finite Volumes 2 Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
23 Test Scenario Adaptive Tsunami Simulation with Finite Volumes 3 Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
24 Test Scenario Adaptive Tsunami Simulation with Finite Volumes 4 Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
25 Test Scenario Adaptive Tsunami Simulation with Finite Volumes 5 Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
26 Test Scenario Adaptive Tsunami Simulation with Finite Volumes 6 Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
27 Test Scenario Adaptive Tsunami Simulation with Finite Volumes 7 Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
28 Test Scenario Adaptive Tsunami Simulation with Finite Volumes 8 Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
29 Test Scenario Adaptive Tsunami Simulation with Finite Volumes 9 Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
30 Test Scenario Adaptive Tsunami Simulation with Finite Volumes 10 Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
31 Summary Summary tsunami modeling on adaptive triangular grids discretization by discontinuous Galerkin method wetting & drying for DG (positivity vs. well-balancing) parallel adaptive triangular grids (Sierpinski SFC) locality properties for cache and parallel efficiency Outlook further development of wetting/drying treatment (positivity AND well-balancing) extension to 2D, realistic tsunami simulations implementation into parallel framework for HPC Vater et al., Parallel Adaptive Tsunami Modelling with DG, PDEs / 20
Efficiency of adaptive mesh algorithms
Efficiency of adaptive mesh algorithms 23.11.2012 Jörn Behrens KlimaCampus, Universität Hamburg http://www.katrina.noaa.gov/satellite/images/katrina-08-28-2005-1545z.jpg Model for adaptive efficiency 10
More informationKlima-Exzellenz in Hamburg
Klima-Exzellenz in Hamburg Adaptive triangular meshes for inundation modeling 19.10.2010, University of Maryland, College Park Jörn Behrens KlimaCampus, Universität Hamburg Acknowledging: Widodo Pranowo,
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 informationShallow Water Simulations on Graphics Hardware
Shallow Water Simulations on Graphics Hardware Ph.D. Thesis Presentation 2014-06-27 Martin Lilleeng Sætra Outline Introduction Parallel Computing and the GPU Simulating Shallow Water Flow Topics of Thesis
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 informationLecture 1: Finite Volume WENO Schemes Chi-Wang Shu
Lecture 1: Finite Volume WENO Schemes Chi-Wang Shu Division of Applied Mathematics Brown University Outline of the First Lecture General description of finite volume schemes for conservation laws The WENO
More informationA new multidimensional-type reconstruction and limiting procedure for unstructured (cell-centered) FVs solving hyperbolic conservation laws
HYP 2012, Padova A new multidimensional-type reconstruction and limiting procedure for unstructured (cell-centered) FVs solving hyperbolic conservation laws Argiris I. Delis & Ioannis K. Nikolos (TUC)
More informationSimulating the effect of Tsunamis within the Build Environment
Simulating the effect of Tsunamis within the Build Environment Stephen Roberts Department of Mathematics Australian National University Ole Nielsen Geospatial and Earth Monitoring Division GeoSciences
More informationCOMPARISON OF CELL-CENTERED AND NODE-CENTERED FORMULATIONS OF A HIGH-RESOLUTION WELL-BALANCED FINITE VOLUME SCHEME: APPLICATION TO SHALLOW WATER FLOWS
XVIII International Conference on Water Resources CMWR 010 J. Carrera (Ed) CIMNE, Barcelona 010 COMPARISON OF CELL-CENTERED AND NODE-CENTERED FORMULATIONS OF A HIGH-RESOLUTION WELL-BALANCED FINITE VOLUME
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 informationDependency of tsunami simulations on advection scheme, grid resolution, bottom friction and topography
IMUM-2010, MIT, 17-20 August Dependency of tsunami simulations on advection scheme, grid resolution, bottom friction and topography C. Wekerle, S. Harig, W. Pranowo, A. Androsov, A. Fuchs, N. Rakowsky,
More informationParallelizing Adaptive Triangular Grids with Refinement Trees and Space Filling Curves
Parallelizing Adaptive Triangular Grids with Refinement Trees and Space Filling Curves Daniel Butnaru butnaru@in.tum.de Advisor: Michael Bader bader@in.tum.de JASS 08 Computational Science and Engineering
More informationMemory Efficient Adaptive Mesh Generation and Implementation of Multigrid Algorithms Using Sierpinski Curves
Memory Efficient Adaptive Mesh Generation and Implementation of Multigrid Algorithms Using Sierpinski Curves Michael Bader TU München Stefanie Schraufstetter TU München Jörn Behrens AWI Bremerhaven Abstract
More informationDevelopment of an operational tsunami model for inclusion into the Indian Ocean Tsunami Early Warning System (GITEWS)
Development of an operational tsunami model for inclusion into the Indian Ocean Tsunami Early Warning System (GITEWS) Jörn Behrens for Polar and Marine Research Acknowledge contributions by Alexey Androsov,
More informationHigh Order Schemes for CFD: A Review. Juan Cheng 1. Institute of Applied Physics and Computational Mathematics, Beijing , China.
High Order Schemes for CFD: A Review Juan Cheng 1 Institute of Applied Physics and Computational Mathematics, Beijing 100088, China and Chi-Wang Shu 2 Division of Applied Mathematics, Brown University,
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 informationThis is a draft. The full paper can be found in Journal of Scientific Computing xx(x):xx xx:
EFFICIENT GPU-IMPLEMENTATION OF ADAPTIVE MESH REFINEMENT FOR THE SHALLOW-WATER EQUATIONS MARTIN L. SÆTRA 1,2, ANDRÉ R. BRODTKORB 3, AND KNUT-ANDREAS LIE 1,3 This is a draft. The full paper can be found
More informationRunge-Kutta discontinuous Galerkin method using a new type of WENO limiters on unstructured mesh 1
Runge-Kutta discontinuous Galerkin method using a new type of WENO limiters on unstructured mesh Jun Zhu, inghui Zhong 3, Chi-Wang Shu 4 and Jianxian Qiu 5 Abstract In this paper we generalize a new type
More informationContents. I The Basic Framework for Stationary Problems 1
page v Preface xiii I The Basic Framework for Stationary Problems 1 1 Some model PDEs 3 1.1 Laplace s equation; elliptic BVPs... 3 1.1.1 Physical experiments modeled by Laplace s equation... 5 1.2 Other
More informationEfficient Storage and Processing of Adaptive Triangular Grids using Sierpinski Curves
Efficient Storage and Processing of Adaptive Triangular Grids using Sierpinski Curves Csaba Attila Vigh, Dr. Michael Bader Department of Informatics, TU München JASS 2006, course 2: Numerical Simulation:
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 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 informationThe Shallow Water Equations and CUDA
The Shallow Water Equations and CUDA HPC - Algorithms and Applications Alexander Pöppl Technical University of Munich Chair of Scientific Computing January 11 th 2017 Last Tutorial Discretized Heat Equation
More informationHomework 4A Due November 7th IN CLASS
CS207, Fall 2014 Systems Development for Computational Science Cris Cecka, Ray Jones Homework 4A Due November 7th IN CLASS Previously, we ve developed a quite robust Graph class to let us use Node and
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 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 informationEvacuate Now? Faster-than-real-time Shallow Water Simulations on GPUs. NVIDIA GPU Technology Conference San Jose, California, 2010 André R.
Evacuate Now? Faster-than-real-time Shallow Water Simulations on GPUs NVIDIA GPU Technology Conference San Jose, California, 2010 André R. Brodtkorb Talk Outline Learn how to simulate a half an hour dam
More informationMichal A. Kopera a,, Francis X. Giraldo a
Analysis of Adaptive Mesh Refinement for IMEX Discontinuous Galerkin Solutions of the Compressible Euler Equations with Application to Atmospheric Simulations Michal A. Kopera a,, Francis X. Giraldo a
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 informationNonoscillatory Central Schemes on Unstructured Triangular Grids for Hyperbolic Systems of Conservation Laws
Nonoscillatory Central Schemes on Unstructured Triangular Grids for Hyperbolic Systems of Conservation Laws Ivan Christov 1,* Bojan Popov 1 Peter Popov 2 1 Department of Mathematics, 2 Institute for Scientific
More informationThe Shallow Water Equations and CUDA
The Shallow Water Equations and CUDA Alexander Pöppl December 9 th 2015 Tutorial: High Performance Computing - Algorithms and Applications, December 9 th 2015 1 Last Tutorial Discretized Heat Equation
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 informationEfficient Algorithmic Approaches for Flow Simulations on Cartesian Grids
Efficient Algorithmic Approaches for Flow Simulations on Cartesian Grids M. Bader, H.-J. Bungartz, B. Gatzhammer, M. Mehl, T. Neckel, T. Weinzierl TUM Department of Informatics Chair of Scientific Computing
More informationComputational Astrophysics 5 Higher-order and AMR schemes
Computational Astrophysics 5 Higher-order and AMR schemes Oscar Agertz Outline - The Godunov Method - Second-order scheme with MUSCL - Slope limiters and TVD schemes - Characteristics tracing and 2D slopes.
More informationLiterature Report. Daniël Pols. 23 May 2018
Literature Report Daniël Pols 23 May 2018 Applications Two-phase flow model The evolution of the momentum field in a two phase flow problem is given by the Navier-Stokes equations: u t + u u = 1 ρ p +
More informationRadial Basis Function-Generated Finite Differences (RBF-FD): New Opportunities for Applications in Scientific Computing
Radial Basis Function-Generated Finite Differences (RBF-FD): New Opportunities for Applications in Scientific Computing Natasha Flyer National Center for Atmospheric Research Boulder, CO Meshes vs. Mesh-free
More informationHigh-order, conservative, finite difference schemes for computational MHD
High-order, conservative, 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
More informationHigh-Order Numerical Algorithms for Steady and Unsteady Simulation of Viscous Compressible Flow with Shocks (Grant FA )
High-Order Numerical Algorithms for Steady and Unsteady Simulation of Viscous Compressible Flow with Shocks (Grant FA9550-07-0195) Sachin Premasuthan, Kui Ou, Patrice Castonguay, Lala Li, Yves Allaneau,
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 informationDebojyoti Ghosh. Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering
Debojyoti Ghosh Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering To study the Dynamic Stalling of rotor blade cross-sections Unsteady Aerodynamics: Time varying
More informationPartial Differential Equations
Simulation in Computer Graphics Partial Differential Equations Matthias Teschner Computer Science Department University of Freiburg Motivation various dynamic effects and physical processes are described
More informationThe Shallow Water Equations and CUDA
The Shallow Water Equations and CUDA Oliver Meister December 17 th 2014 Tutorial Parallel Programming and High Performance Computing, December 17 th 2014 1 Last Tutorial Discretized Heat Equation System
More informationNonoscillatory Central Schemes on Unstructured Triangulations for Hyperbolic Systems of Conservation Laws
Nonoscillatory Central Schemes on Unstructured Triangulations for Hyperbolic Systems of Conservation Laws Ivan Christov Bojan Popov Department of Mathematics, Texas A&M University, College Station, Texas
More informationFluent User Services Center
Solver Settings 5-1 Using the Solver Setting Solver Parameters Convergence Definition Monitoring Stability Accelerating Convergence Accuracy Grid Independence Adaption Appendix: Background Finite Volume
More information3D NURBS-ENHANCED FINITE ELEMENT METHOD
7th Workshop on Numerical Methods in Applied Science and Engineering (NMASE 8) Vall de Núria, 9 a 11 de enero de 28 c LaCàN, www.lacan-upc.es 3D NURBS-ENHANCED FINITE ELEMENT METHOD R. Sevilla, S. Fernández-Méndez
More informationShallow Water Equation simulation with Sparse Grid Combination Technique
GUIDED RESEARCH Shallow Water Equation simulation with Sparse Grid Combination Technique Author: Jeeta Ann Chacko Examiner: Prof. Dr. Hans-Joachim Bungartz Advisors: Ao Mo-Hellenbrand April 21, 2017 Fakultät
More informationTesting refinement criteria in adaptive Discontinuous Galerkin simulations of dry atmospheric convection
Testing refinement criteria in adaptive Discontinuous Galerkin simulations of dry atmospheric convection Andreas Müller a,, Jörn Behrens b, Francis X. Giraldo c, Volkmar Wirth a a Institute for Atmospheric
More informationComparison between adaptive and uniform discontinuous Galerkin simulations in dry 2D bubble experiments
Calhoun: The NPS Institutional Archive DSpace Repository Faculty and Researchers Faculty and Researchers Collection 2012-11-08 Comparison between adaptive and uniform discontinuous Galerkin simulations
More information1 Exercise: Heat equation in 2-D with FE
1 Exercise: Heat equation in 2-D with FE Reading Hughes (2000, sec. 2.3-2.6 Dabrowski et al. (2008, sec. 1-3, 4.1.1, 4.1.3, 4.2.1 This FE exercise and most of the following ones are based on the MILAMIN
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationDonna Calhoun. Boise State University
Wave propagation algorithms in two dimensions, adaptive refinement and wellbalancing Donna Calhoun Boise State University PASI 2013 Valparaiso, Chile January 2-13 Plan How do these Riemann solvers make
More informationA Novel Approach to High Speed Collision
A Novel Approach to High Speed Collision Avril Slone University of Greenwich Motivation High Speed Impact Currently a very active research area. Generic projectile- target collision 11 th September 2001.
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 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 informationINTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 3, 2012
INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 3, 2012 Copyright 2010 All rights reserved Integrated Publishing services Research article ISSN 0976 4399 Efficiency and performances
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 informationADAPTIVE MESH REFINEMENT TECHNIQUES FOR HIGH-ORDER FINITE-VOLUME WENO SCHEMES
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 informationDISCONTINUOUS GALERKIN SHALLOW WATER SOLVER ON CUDA ARCHITECTURES
9 th International Conference on Hydroinformatics HIC 2010, Tianjin, CHINA DISCONTINUOUS GALERKIN SHALLOW WATER SOLVER ON CUDA ARCHITECTURES D. SCHWANENBERG Operational Water Management, Deltares, Rotterdamseweg
More informationQUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD
QUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD Keisuke Yoshida, Tadaharu Ishikawa Dr. Eng., Tokyo Institute
More informationNIA CFD Seminar, October 4, 2011 Hyperbolic Seminar, NASA Langley, October 17, 2011
NIA CFD Seminar, October 4, 2011 Hyperbolic Seminar, NASA Langley, October 17, 2011 First-Order Hyperbolic System Method If you have a CFD book for hyperbolic problems, you have a CFD book for all problems.
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 informationmulti-moments moments or multi-moment moment constraints:
Isaac Newton Institute for Mathematical Sciences Program on Multiscale Numerics for the Atmosphere and Ocean Cambridge, Nov. 16, 2012 High order finite volume methods using multi-moments moments or multi-moment
More informationUniversity of Reading. Discretization On Non-uniform Meshes: Tests solving shallow-water equations
University of Reading School of Mathematics, Meteorology and Physics Discretization On Non-uniform Meshes: Tests solving shallow-water equations By FAWZI B ALBUSAIDI August 2008 This dissertation is submitted
More informationAn Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems
An Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems Long Chen University of California, Irvine chenlong@math.uci.edu Joint work with: Huayi Wei (Xiangtan University),
More informationTsunami coastal impact The use of VOF-URANS methods with examples
Tsunami coastal impact The use of VOF-URANS methods with examples Richard MARCER 2 Content Coastal impact : physics to simulate EOLE CFD : a 3D Navier-Stokes code Focus on the VOF free surface method Examples
More informationDual Interpolants for Finite Element Methods
Dual Interpolants for Finite Element Methods Andrew Gillette joint work with Chandrajit Bajaj and Alexander Rand Department of Mathematics Institute of Computational Engineering and Sciences University
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 informationLax-Wendroff and McCormack Schemes for Numerical Simulation of Unsteady Gradually and Rapidly Varied Open Channel Flow
Archives of Hydro-Engineering and Environmental Mechanics Vol. 60 (2013), No. 1 4, pp. 51 62 DOI: 10.2478/heem-2013-0008 IBW PAN, ISSN 1231 3726 Lax-Wendroff and McCormack Schemes for Numerical Simulation
More informationA New Trouble-Cell Indicator for Discontinuous Galerkin Methods for. Hyperbolic Conservation Laws ABSTRACT
A New Trouble-Cell Indicator for Discontinuous Galerkin Methods for Hyperbolic Conservation Laws Guosheng Fu and Chi-Wang Shu ABSTRACT We introduce a new troubled-cell indicator for the discontinuous Galerkin
More information2.7 Cloth Animation. Jacobs University Visualization and Computer Graphics Lab : Advanced Graphics - Chapter 2 123
2.7 Cloth Animation 320491: Advanced Graphics - Chapter 2 123 Example: Cloth draping Image Michael Kass 320491: Advanced Graphics - Chapter 2 124 Cloth using mass-spring model Network of masses and springs
More informationAlgebraic Flux Correction Schemes for High-Order B-Spline Based Finite Element Approximations
Algebraic Flux Correction Schemes for High-Order B-Spline Based Finite Element Approximations Numerical Analysis group Matthias Möller m.moller@tudelft.nl 7th DIAM onderwijs- en onderzoeksdag, 5th November
More informationDiscontinuous Fluctuation Distribution
Discontinuous Fluctuation Distribution Matthew Hubbard School of Computing, University of Leeds, Leeds, LS2 9JT, UK. Abstract This paper describes a new numerical scheme for the approximation of steady
More informationADH Rules of Thumb. US Army Corps of Engineers. Coastal and Hydraulics Laboratory Engineer Research & Development Center
ADH Rules of Thumb Refinement Levels Refinement Levels (ML mat# lev#) allows an element to be divided into a maximum of 2lev# of elements, you may not use all of them depending on the flow conditions and
More informationAN ADAPTIVE MESH REDISTRIBUTION ALGORITHM FOR CONVECTION-DOMINATED PROBLEMS. Zheng-Ru Zhang and Tao Tang
COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume, Number 3, september pp. 34 357 AN ADAPTIVE MESH REDISTRIBUTION ALGORITHM FOR CONVECTION-DOMINATED PROBLEMS Zheng-Ru Zhang
More informationRunge Kutta Methods Optimized For Advection Problems
Runge Kutta Methods Optimized For Advection Problems Brown University December 9, 2009 Introduction Classical Runge Kutta Methods Overview Standard Runge Kutta Schemes for Advection Equation Low Dissipation
More informationHierarchical Reconstruction for Spectral Volume Method on Unstructured Grids
Hierarchical Reconstruction for Spectral Volume Method on Unstructured Grids Zhiliang Xu, Yingjie Liu and Chi-Wang Shu April 4, 2009 Abstract The hierarchical reconstruction (HR) [, 24] is applied to a
More informationNUMERICAL SIMULATION OF THE SHALLOW WATER EQUATIONS USING A TIME-CENTERED SPLIT-IMPLICIT METHOD
18th Engineering Mechanics Division Conference (EMD007) NUMERICAL SIMULATION OF THE SHALLOW WATER EQUATIONS USING A TIME-CENTERED SPLIT-IMPLICIT METHOD Abstract S. Fu University of Texas at Austin, Austin,
More informationLevel set methods Formulation of Interface Propagation Boundary Value PDE Initial Value PDE Motion in an externally generated velocity field
Level Set Methods Overview Level set methods Formulation of Interface Propagation Boundary Value PDE Initial Value PDE Motion in an externally generated velocity field Convection Upwind ddifferencingi
More informationRunge Kutta discontinuous Galerkin method using WENO limiters II: unstructured meshes
Runge Kutta discontinuous Galerkin method using WENO limiters II: unstructured meshes Jun Zhu, Jianxian Qiu,Chi-WangShu 3 and Michael Dumbser 4 Abstract In [], Qiu and Shu investigated using weighted essentially
More informationTechnische Universität München Fakultät für Informatik. Computational Science and Engineering (Int. Master s Program)
Technische Universität München Fakultät für Informatik Computational Science and Engineering (Int. Master s Program) Parallel Refinement and Coarsening of recursively structured Adaptive triangular grids
More informationComparison of Central and Upwind Flux Averaging in Overlapping Finite Volume Methods for Simulation of Super-Critical Flow with Shock Waves
Proceedings of the 9th WSEAS International Conference on Applied Mathematics, Istanbul, Turkey, May 7, 6 (pp55665) Comparison of and Flux Averaging in Overlapping Finite Volume Methods for Simulation of
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 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-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 informationOpenFOAM in Wave and Offshore CFD
OpenFOAM in Wave and Offshore CFD Capabilities of the Naval Hydro Pack Hrvoje Jasak Wikki Ltd. United Kingdom Faculty of Mechanical Engineering and Naval Architecture, Uni Zagreb, Croatia University of
More informationForest-of-octrees AMR: algorithms and interfaces
Forest-of-octrees AMR: algorithms and interfaces Carsten Burstedde joint work with Omar Ghattas, Tobin Isaac, Georg Stadler, Lucas C. Wilcox Institut für Numerische Simulation (INS) Rheinische Friedrich-Wilhelms-Universität
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 informationDiscontinuous Galerkin Finite Element Methods for Shallow Water Flow: Developing a Computational Infrastructure for Mixed Element Meshes
Discontinuous Galerkin Finite Element Methods for Shallow Water Flow: Developing a Computational Infrastructure for Mixed Element Meshes A Thesis Presented in Partial Fulfillment of the Requirements for
More informationAdaptive-Mesh-Refinement Pattern
Adaptive-Mesh-Refinement Pattern I. Problem Data-parallelism is exposed on a geometric mesh structure (either irregular or regular), where each point iteratively communicates with nearby neighboring points
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 informationFAST ALGORITHMS FOR CALCULATIONS OF VISCOUS INCOMPRESSIBLE FLOWS USING THE ARTIFICIAL COMPRESSIBILITY METHOD
TASK QUARTERLY 12 No 3, 273 287 FAST ALGORITHMS FOR CALCULATIONS OF VISCOUS INCOMPRESSIBLE FLOWS USING THE ARTIFICIAL COMPRESSIBILITY METHOD ZBIGNIEW KOSMA Institute of Applied Mechanics, Technical University
More informationNUMERICAL INVESTIGATION OF ONE- DIMENSIONAL SHALLOW WATER FLOW AND SEDIMENT TRANSPORT MODEL USING DISCONTINOUOUS GALERKIN METHOD
Clemson University TigerPrints All Dissertations Dissertations 5-2014 NUMERICAL INVESTIGATION OF ONE- DIMENSIONAL SHALLOW WATER FLOW AND SEDIMENT TRANSPORT MODEL USING DISCONTINOUOUS GALERKIN METHOD Farzam
More informationEulerian Techniques for Fluid-Structure Interactions - Part II: Applications
Published in Lecture Notes in Computational Science and Engineering Vol. 103, Proceedings of ENUMATH 2013, pp. 755-762, Springer, 2014 Eulerian Techniques for Fluid-Structure Interactions - Part II: Applications
More informationDDFV Schemes for the Euler Equations
DDFV Schemes for the Euler Equations Christophe Berthon, Yves Coudière, Vivien Desveaux Journées du GDR Calcul, 5 juillet 2011 Introduction Hyperbolic system of conservation laws in 2D t W + x f(w)+ y
More informationNumerical Investigation of Transonic Shock Oscillations on Stationary Aerofoils
Numerical Investigation of Transonic Shock Oscillations on Stationary Aerofoils A. Soda, T. Knopp, K. Weinman German Aerospace Center DLR, Göttingen/Germany Symposium on Hybrid RANS-LES Methods Stockholm/Sweden,
More informationNON-OSCILLATORY HIERARCHICAL RECONSTRUCTION FOR CENTRAL AND FINITE VOLUME SCHEMES
NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION FOR CENTRAL AND FINITE VOLUME SCHEMES YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. This is the continuation of the paper central discontinuous
More informationAdaptive Discontinuous Galerkin Method for the Shallow Water Equations
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2003; 00:1 20 [Version: 2002/09/18 v1.01] Adaptive Discontinuous Galerkin Method for the Shallow Water Equations Jean-François
More informationSpectral(Finite) Volume Method for Conservation Laws on Unstructured Grids
Journal of Computational Physics 179, 665 697 (2002) doi:10.1006/jcph.2002.7082 Spectral(Finite) Volume Method for Conservation Laws on Unstructured Grids II. Extension to Two-Dimensional Scalar Equation
More informationLoad-balancing multi-gpu shallow water simulations on small clusters
Load-balancing multi-gpu shallow water simulations on small clusters Gorm Skevik master thesis autumn 2014 Load-balancing multi-gpu shallow water simulations on small clusters Gorm Skevik 1st August 2014
More informationAdaptive Mesh Astrophysical Fluid Simulations on GPU. San Jose 10/2/2009 Peng Wang, NVIDIA
Adaptive Mesh Astrophysical Fluid Simulations on GPU San Jose 10/2/2009 Peng Wang, NVIDIA Overview Astrophysical motivation & the Enzo code Finite volume method and adaptive mesh refinement (AMR) CUDA
More information