A finite element solver for ice sheet dynamics to be integrated with MPAS
|
|
- Felix Kelley
- 5 years ago
- Views:
Transcription
1 A finite element solver for ice sheet dynamics to be integrated with MPAS Mauro Perego in collaboration with FSU, ORNL, LANL, Sandia February 6, CESM LIWG Meeting, Boulder (CO), 202
2 Outline Introduction of finite element ice sheets component Integration of finite element ice sheet component with MPAS Results on ISMIP HOM and Greenland/Antarctica geometries Non linear solvers: features and issues
3 Ice Sheet Modeling Main components of an ice model: Ice flow equations (momentum and mass balance) r ¾ = ½g and r u = 0; with ¾ = pi = 2¹(") _ "_ pi, where ¹ viscosity, "_ shear rate Model for the evolution of the boundaries (thickness evolution = Hf lux Z u dz z ½cu rt + 2"¾ _
4 Ice Sheet Modeling Reference model: FULL STOKES Approximations based on their accuracy with respect to the ice sheet aspect ratio ± O(± 2 ) FO, Blatter Pattyn first order model2 (3D PDE, in horizontal velocities) O(±) Zeroth order, depth integrated models: SIA, Shallow Ice Approximation (slow sliding regimes), SSA Shallow Shelf Approximation (2D PDE) (fast sliding regimes) ' O(± 2 ) High order, depth integrated (2D) models: LL23, (LL)... Gagliardini and Zwinger, The Cryosphere. 2 Dukowicz, Price and Lipscomb, 200. J. Glaciol. 3 Schoof and Hindmarsh, 200. Q. J. Mech. Appl. Math.
5 Implementation Overview ice sheets FE dynamics component LifeV: Parallel, object oriented, C++ Finite Element Library: linear and quadratic finite elements assembling of finite element matrices handling of boundary conditions
6 Implementation Overview Trilinos: Parallel Data Structures (EPETRA) Parallel Linear Solvers (GMRES, CG...) Preconditioners (Multilevel, Multigrid, Incomplete LU) Nonlinear Solvers (NOX package: Newton, JFNK methods) ice sheets FE dynamics component LifeV: Parallel, object oriented, C++ Finite Element Library: linear and quadratic finite elements assembling of finite element matrices handling of boundary conditions Work in Progress MPAS (land ice component): Voronoi unstructured grids Evolution equation solvers (temperature and thickness equation)...
7 Interface MPAS LIFEV MPAS LIFEV 2D CVT mesh (Stereographic projection) thickness/elevation/layers bedrock sliding coefficient Solver options: model (FO, LL2, SSA, SIA) nonlinear solver (Newton, Picard, JFNK) Boundary condition (free-slip, no-slip, robin, coulomb) velocity heat dissipation viscosity ice-sheets component Land ice component temperature/ice flow factor
8 Interface Grids MPAS CVT Mesh (OK for Finite Volumes)
9 Interface Grids MPAS CVT Mesh (OK for Finite Volumes) velocity velocity Temperature Lifev Dual triangular Mesh (OK for Finite Elements)
10 Interface Grids MPAS CVT Mesh (OK for Finite Volumes) velocity velocity Temperature Based on 2D grid and thickness and layers build vertically structured 3D grid. Build prisms with triangular base and split them in tetrahedra. Lifev Dual triangular Mesh (OK for Finite Elements)
11 ISMIP HOM: Test C Surface velocity L = 5 km L = 40 km L = 0 km L = 80 km L = 20 km L = 60 km Perego, Gunzburger, Burkardt, Journal of Glaciology, 202.
12 Towards realistic simulations: Greenland, 5km resolution, sliding case Comparisons between different models SIA SSA LL2 FO
13 Towards realistic simulations: Greenland, 2km resolution, no slip SIA FO
14 Towards realistic simulations: Antarctica, 5km resolution, no slip SIA FO
15 Convergence of the nonlinear method (test C ISMIP HOM) Solve F (u) = 0: Newton method: J k (±u)k+ = F (uk ): (NOX library) The exact Jacobian matrix must be assembled at each nonlinear iteration. In order to increase robustness of Newton method, the Newton step is halved to achieve monotonic convergence. FO LL2 CPU time: 93s CPU time: 2s CPU time: 4s CPU time: 0.98s
16 Greenland, FO: convergence of the nonlinear method
17 First Order Equations FO is a nonlinear system of elliptic equations in the horizontal velocities: > > < r (2¹ "_ ) = > > : r (2¹ "_2 ) = n ( n ) ¹ = A "_e 2 q "_e = "_2xx + "_2yy + "_xx "_yy + "_2xy + "_2xz + "_2yz where s is the ice surface and, "_i;j = (@j ui uj ) ; 2 i; j 2 fx; y; zg; 2 6 "_ = 6 4 2"_xx + "_yy "_xy "_xz 3 Goldsby D. L. Kohlsted, JGR, 200. "_yx ; "_2 = 6 "_xx + 2"_yy "_yz Remark The nonlinear viscosity ¹ is singular when "_e = 0, however, ¹ "_ is not singular and the PDE is well de ned. ³p ( n ) ( ) n Viscosity regularization: "_e ¼ "_2e + ± 2 2
18 Greenland, FO: convergence of the nonlinear method Simpli ed Problem (similar kind of nonlinearity): x jxj k+ Picard method: x ( n ) solution: = CjCjn n ( ): = Cjx j k+ Newton method: x = C: k n ): ( = ( n)x + ncjx j k k µ x 2 Picard convergence: lim = = when n = 3 k k! x n 3 µ k+ x Newton convergence: lim = k 2 k! (x ) 2 n k+ regularization: x jxj ( n ³p ( n ) )¼x x2 + ± 2
19 Greenland, FO: convergence of the nonlinear method (comparison using different regularizations) Err: m/a Err: 7.2e 4 m/a Err: 0.34 m/a Viscosity regularization: ( n ) "_e ¼ ³p "_2e + ± 2 ( n )
20 Greenland, FO: convergence of the nonlinear method (comparison using different regularizations) ± variable : At each newton iteration ± is decreased from e-4 to e-9
21 Greenland, FO: convergence of the nonlinear method (LOCA continuation method) ±=e-9 ±=e-9 ±=e-4 ±=e-9 ±=.3e-7 ±=e-5 The parameter ± is decreased by LOCA from e-4 to e-9
22 Acknowledgment J. Burkardt, M. Gunzburger (FSU) B. Lipscomb, S. Price, X. Asay Davis, T. Ringler (LANL) K. Evans, J. Nichols (ORNL) A. Salinger (Sandia)
THE SCIENCE OF ICE SHEETS: THE MATHEMATICAL MODELING AND COMPUTATIONAL SIMULATION OF ICE FLOWS
THE SCIENCE OF ICE SHEETS: THE MATHEMATICAL MODELING AND COMPUTATIONAL SIMULATION OF ICE FLOWS Max Gunzburger Department of Scientific Computing, Florida State University Collaborators: Mauro Perego (Florida
More informationA posteriori estimator for the accuracy of the shallow shelf approximation
Math Geosci manuscript No. will be inserted by the editor A posteriori estimator for the accuracy of the shallow shelf approximation Guillaume Jouvet Marco Picasso Received: date / Accepted: date Abstract
More informationIce Sheet Modeling and Sea Level Rise. William Lipscomb, NCAR CESM Sea Level Session 10 January 2018
Ice Sheet Modeling and Sea Level Rise William Lipscomb, NCAR CESM Sea Level Session 10 January 2018 Ice sheets in IPCC AR4 The IPCC Fourth Assessment Report (AR4) projected 0.18 to 0.59 m of sea level
More informationElmer/Ice Chinese-Finnish seminar, October 30 th 2012
Elmer/Ice Chinese-Finnish seminar, October 30 th 2012 NEW GENERATION ICE SHEET MODELS Shallow ice approximation is so 80 s! Thomas Zwinger 1) and... 1) CSC IT Center for Science Ltd. Espoo, Finland ...
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 informationModeling land ice in the Community Earth System Model
Modeling land ice in the Community Earth System Model William Lipscomb Los Alamos National Laboratory CESM Tutorial 5 August 2011 Outline Motivation for simulating land ice in Earth-system models Introduction
More informationIterative Methods and Multigrid
Iterative Methods and Multigrid Part 3: Preconditioning 2 Eric de Sturler Preconditioning The general idea behind preconditioning is that convergence of some method for the linear system Ax = b can be
More informationThe Effect of a Basal-friction Parameterization on. Grounding-line Dynamics in Ice-sheet Models.
The Effect of a Basal-friction Parameterization on Grounding-line Dynamics in Ice-sheet Models. by Gunter Leguy Submitted in Partial Fulfillment of the Requirements for the Degree of Ph.D of Applied &
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
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 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 informationIMPLEMENTATION OF A VERTICALLY INTEGRATED ICE SHEET MOMENTUM BALANCE MODEL
University of Montana ScholarWorks at University of Montana Graduate Student Theses, Dissertations, & Professional Papers Graduate School 013 IMPLEMENTATION OF A VERTICALLY INTEGRATED ICE SHEET MOMENTUM
More informationUT-LANL DOE/ASCR ISICLES Project PARADICE: Uncertainty Quantification for Large-Scale Ice Sheet Modeling and Simulation.
UT-LANL DOE/ASCR ISICLES Project PARADICE: Uncertainty Quantification for Large-Scale Ice Sheet Modeling and Simulation Omar Ghattas Representing the rest of the UT-Austin/LANL team: Don Blankenship Tan
More informationDiscretization of PDEs and Tools for the Parallel Solution of the Resulting Systems
Discretization of PDEs and Tools for the Parallel Solution of the Resulting Systems Stan Tomov Innovative Computing Laboratory Computer Science Department The University of Tennessee Wednesday April 4,
More informationModelling of surface to basal hydrology across the Russell Glacier Catchment
Modelling of surface to basal hydrology across the Russell Glacier Catchment Sam GAP Modelling Workshop, Toronto November 2010 Collaborators Alun Hubbard Centre for Glaciology Institute of Geography and
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationEfficient FEM-multigrid solver for granular material
Efficient FEM-multigrid solver for granular material S. Mandal, A. Ouazzi, S. Turek Chair for Applied Mathematics and Numerics (LSIII), TU Dortmund STW user committee meeting Enschede, 25th September,
More informationOUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative methods ffl Krylov subspace methods ffl Preconditioning techniques: Iterative methods ILU
Preconditioning Techniques for Solving Large Sparse Linear Systems Arnold Reusken Institut für Geometrie und Praktische Mathematik RWTH-Aachen OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative
More informationModelling ice shelf basal melt with Glimmer-CISM coupled to a meltwater plume model
Modelling ice shelf basal melt with Glimmer-CISM coupled to a meltwater plume model Carl Gladish NYU CIMS February 17, 2010 Carl Gladish (NYU CIMS) Glimmer-CISM + Plume February 17, 2010 1 / 24 Acknowledgements
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationVerifying PISM, once known as COMMVNISM
Verifying PISM, once known as COMMVNISM A C++ Object-oriented Multi-Modal, Verifiable Numerical Ice Sheet Model Jed Brown 1 Ed Bueler 1 Craig Lingle 2 1 Department of Mathematics University of Alaska Fairbanks
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 informationPostprint.
http://www.diva-portal.org Postprint This is the accepted version of a paper published in Journal of Computational Physics. This paper has been peer-reviewed but does not include the final publisher proof-corrections
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 informationverification, validation, and basal strength in models for the present state of ice sheets
verification, validation, and basal strength in models for the present state of ice sheets Ed Bueler 1 Constantine Khroulev 2 Andy Aschwanden 3 1 Dept of Mathematics and Statistics, University of Alaska
More informationSolving Large Nonlinear Sparse Systems
Solving Large Nonlinear Sparse Systems Fred W. Wubs and Jonas Thies Computational Mechanics & Numerical Mathematics University of Groningen, the Netherlands f.w.wubs@rug.nl Centre for Interdisciplinary
More informationAn overview of the. package for non-conformal mesh tying or simple contact problems
An overview of the package for non-conformal mesh tying or simple contact problems Glen Hansen 2012 Trilinos User Group Meeting Oct. 30, 2012 SAND 2012-9202P (UUR) Sandia National Laboratories is a multi-program
More informationSplitting Iteration Methods for Positive Definite Linear Systems
Splitting Iteration Methods for Positive Definite Linear Systems Zhong-Zhi Bai a State Key Lab. of Sci./Engrg. Computing Inst. of Comput. Math. & Sci./Engrg. Computing Academy of Mathematics and System
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 informationMECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso
MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Preliminary Math Concept of Stress Stress Components Equilibrium
More informationIndefinite and physics-based preconditioning
Indefinite and physics-based preconditioning Jed Brown VAW, ETH Zürich 2009-01-29 Newton iteration Standard form of a nonlinear system F (u) 0 Iteration Solve: Update: J(ũ)u F (ũ) ũ + ũ + u Example (p-bratu)
More informationWhat s New in Isorropia?
What s New in Isorropia? Erik Boman Cedric Chevalier, Lee Ann Riesen Sandia National Laboratories, NM, USA Trilinos User s Group, Nov 4-5, 2009. SAND 2009-7611P. Sandia is a multiprogram laboratory operated
More informationAUTOMATIC DIFFERENTIATION FOR NONLINEAR, MULTIPHYSICS SIMULATION
AUTOMATIC DIFFERENTIATION FOR NONLINEAR, MULTIPHYSICS SIMULATION Derek Gaston October 24 2016 Multiphysics Object Oriented Simulation Environment MOOSE is a finite-element, multiphysics framework that
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 informationMultilevel low-rank approximation preconditioners Yousef Saad Department of Computer Science and Engineering University of Minnesota
Multilevel low-rank approximation preconditioners Yousef Saad Department of Computer Science and Engineering University of Minnesota SIAM CSE Boston - March 1, 2013 First: Joint work with Ruipeng Li Work
More informationCoupling of Multi fidelity Models Applications to PNP cdft and local nonlocal Poisson equations
Coupling of Multi fidelity Models Applications to PNP cdft and local nonlocal Poisson equations P. Bochev, J. Cheung, M. D Elia, A. Frishknecht, K. Kim, M. Parks, M. Perego CM4 summer school, Stanford,
More informationANALYSIS AND NUMERICAL METHODS FOR SOME CRACK PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING, SERIES B Volume 2, Number 2-3, Pages 155 166 c 2011 Institute for Scientific Computing and Information ANALYSIS AND NUMERICAL METHODS FOR SOME
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationUQ Benchmark Problems for Multiphysics Modeling
SIAMUQ UQ Challenge Benchmarks UQ Benchmark Problems for Multiphysics Modeling Maarten Arnst March 31, 2014 SIAMUQ UQ Challenge Benchmarks 1 / 25 Motivation Previous presentation at USNCCM2013 UQ Challenge
More informationIce in the climate system. Summary so far. Today. The Cryosphere. 1. Climate history of the Earth. 2. Paleo observations (1)
Ice in the climate system 1. Climate history of the Earth 2. Paleo observations (1) 3. Paleo observations (2) 4. Ice ages 5. Climate sensitivity 6. Ice in the climate system Summary so far Radiation (Milankovitch
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
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 informationThe solution of the discretized incompressible Navier-Stokes equations with iterative methods
The solution of the discretized incompressible Navier-Stokes equations with iterative methods Report 93-54 C. Vuik Technische Universiteit Delft Delft University of Technology Faculteit der Technische
More informationJae Heon Yun and Yu Du Han
Bull. Korean Math. Soc. 39 (2002), No. 3, pp. 495 509 MODIFIED INCOMPLETE CHOLESKY FACTORIZATION PRECONDITIONERS FOR A SYMMETRIC POSITIVE DEFINITE MATRIX Jae Heon Yun and Yu Du Han Abstract. We propose
More informationLecture 8: Fast Linear Solvers (Part 7)
Lecture 8: Fast Linear Solvers (Part 7) 1 Modified Gram-Schmidt Process with Reorthogonalization Test Reorthogonalization If Av k 2 + δ v k+1 2 = Av k 2 to working precision. δ = 10 3 2 Householder Arnoldi
More informationn i,j+1/2 q i,j * qi+1,j * S i+1/2,j
Helsinki University of Technology CFD-group/ The Laboratory of Applied Thermodynamics MEMO No CFD/TERMO-5-97 DATE: December 9,997 TITLE A comparison of complete vs. simplied viscous terms in boundary layer
More informationEfficient BDF time discretization of the Navier Stokes equations with VMS LES modeling in a High Performance Computing framework
Introduction VS LES/BDF Implementation Results Conclusions Swiss Numerical Analysis Day 2015, Geneva, Switzerland, 17 April 2015. Efficient BDF time discretization of the Navier Stokes equations with VS
More informationarxiv: v1 [physics.comp-ph] 26 Mar 2015
arxiv:1503.07743v1 [physics.comp-ph] 26 Mar 2015 Dynamically coupling the non-linear Stokes equations with the Shallow Ice Approximation in glaciology: Description and first applications of the ISCAL method
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 2014 Quiz #2 April 14, 2014 Name: SOLUTION ID#: PS1.:
More informationParameterization of lateral drag in flowline models of glacier dynamics
Journal of Glaciology, Vol. 58, No. 212, 2012 doi: 10.3189/2012JoG12J018 1119 Parameterization of lateral drag in flowline models of glacier dynamics Surendra ADHIKARI, Shawn J. MARSHALL Department of
More informationMultigrid Methods and their application in CFD
Multigrid Methods and their application in CFD Michael Wurst TU München 16.06.2009 1 Multigrid Methods Definition Multigrid (MG) methods in numerical analysis are a group of algorithms for solving differential
More informationLecture 17: Iterative Methods and Sparse Linear Algebra
Lecture 17: Iterative Methods and Sparse Linear Algebra David Bindel 25 Mar 2014 Logistics HW 3 extended to Wednesday after break HW 4 should come out Monday after break Still need project description
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 informationCitation for the original published paper (version of record):
http://www.diva-portal.org Postprint This is the accepted version of a paper published in Journal of Computational Physics. This paper has been peer-reviewed but does not include the final publisher proof-corrections
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 informationA DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any
Y Y Y X X «/ YY Y Y ««Y x ) & \ & & } # Y \#$& / Y Y X» \\ / X X X x & Y Y X «q «z \x» = q Y # % \ & [ & Z \ & { + % ) / / «q zy» / & / / / & x x X / % % ) Y x X Y $ Z % Y Y x x } / % «] «] # z» & Y X»
More informationAE/ME 339. Computational Fluid Dynamics (CFD) K. M. Isaac. Momentum equation. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.
AE/ME 339 Computational Fluid Dynamics (CFD) 9//005 Topic7_NS_ F0 1 Momentum equation 9//005 Topic7_NS_ F0 1 Consider the moving fluid element model shown in Figure.b Basis is Newton s nd Law which says
More informationIntegration of Sequential Quadratic Programming and Domain Decomposition Methods for Nonlinear Optimal Control Problems
Integration of Sequential Quadratic Programming and Domain Decomposition Methods for Nonlinear Optimal Control Problems Matthias Heinkenschloss 1 and Denis Ridzal 2 1 Department of Computational and Applied
More informationDynamics of Glaciers
Dynamics of Glaciers McCarthy Summer School 01 Andy Aschwanden Arctic Region Supercomputing Center University of Alaska Fairbanks, USA June 01 Note: This script is largely based on the Physics of Glaciers
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
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 informationGlacier Hydrology II: Theory and Modeling
Glacier Hydrology II: Theory and Modeling McCarthy Summer School 2018 Matt Hoffman Gwenn Flowers, Simon Fraser Operated by Los Alamos National Security, LLC for the U.S. Department of Energy's NNSA Observations
More informationIntegration of PETSc for Nonlinear Solves
Integration of PETSc for Nonlinear Solves Ben Jamroz, Travis Austin, Srinath Vadlamani, Scott Kruger Tech-X Corporation jamroz@txcorp.com http://www.txcorp.com NIMROD Meeting: Aug 10, 2010 Boulder, CO
More informationConvergence and accuracy of sea ice dynamics solvers. Martin Losch (AWI) Annika Fuchs (AWI), Jean-François Lemieux (EC)
Convergence and accuracy of sea ice dynamics solvers Martin Losch (AWI) Annika Fuchs (AWI), Jean-François Lemieux (EC) Sea ice thickness with MITgcm (JPL) Review of 2D momentum equations m Du Dt = mf k
More informationNumerical Simulation of Rhone s Glacier and G. Jouvet 1 M. Picasso 1 J. Rappaz 1 H. Blatter 2 M. Funk 3 M. Huss 3. Paris, ENPC, jan.
between 1874 and 2100 G. Jouvet 1 M. Picasso 1 J. Rappaz 1 H. Blatter 2 M. Funk 3 M. Huss 3 1 Institut d Analyse et Calcul Scientifique EPF Lausanne, Switzerland 2 Institute for Atmospheric and Climate
More informationMulti-Modal Flow in a Thermocoupled Model of the Antarctic Ice Sheet, with Verification
Multi-Modal Flow in a Thermocoupled Model of the Antarctic Ice Sheet, with Verification Craig Lingle 1 Jed Brown 2 Ed Bueler 2 1 Geophysical Institute University of Alaska Fairbanks, USA 2 Department of
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 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 informationNUMERICAL METHOD FOR THREE DIMENSIONAL STEADY-STATE TWO-PHASE FLOW CALCULATIONS
' ( '- /A NUMERCAL METHOD FOR THREE DMENSONAL STEADY-STATE TWO-PHASE FLOW CALCULATONS P. Raymond,. Toumi (CEA) CE-Saclay DMT/SERMA 91191 Gif sur Yvette, FRANCE Tel: (331) 69 08 26 21 / Fax: 69 08 23 81
More informationAE/ME 339. K. M. Isaac Professor of Aerospace Engineering. 12/21/01 topic7_ns_equations 1
AE/ME 339 Professor of Aerospace Engineering 12/21/01 topic7_ns_equations 1 Continuity equation Governing equation summary Non-conservation form D Dt. V 0.(2.29) Conservation form ( V ) 0...(2.33) t 12/21/01
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 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 informationAn Empirical Comparison of Graph Laplacian Solvers
An Empirical Comparison of Graph Laplacian Solvers Kevin Deweese 1 Erik Boman 2 John Gilbert 1 1 Department of Computer Science University of California, Santa Barbara 2 Scalable Algorithms Department
More informationPoisson Equation in 2D
A Parallel Strategy Department of Mathematics and Statistics McMaster University March 31, 2010 Outline Introduction 1 Introduction Motivation Discretization Iterative Methods 2 Additive Schwarz Method
More informationUsing AmgX to accelerate a PETSc-based immersed-boundary method code
29th International Conference on Parallel Computational Fluid Dynamics May 15-17, 2017; Glasgow, Scotland Using AmgX to accelerate a PETSc-based immersed-boundary method code Olivier Mesnard, Pi-Yueh Chuang,
More informationFast 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 informationFinite Element Method (FEM)
Finite Element Method (FEM) The finite element method (FEM) is the oldest numerical technique applied to engineering problems. FEM itself is not rigorous, but when combined with integral equation techniques
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 informationStructure preserving preconditioner for the incompressible Navier-Stokes equations
Structure preserving preconditioner for the incompressible Navier-Stokes equations Fred W. Wubs and Jonas Thies Computational Mechanics & Numerical Mathematics University of Groningen, the Netherlands
More informationThe flexible incomplete LU preconditioner for large nonsymmetric linear systems. Takatoshi Nakamura Takashi Nodera
Research Report KSTS/RR-15/006 The flexible incomplete LU preconditioner for large nonsymmetric linear systems by Takatoshi Nakamura Takashi Nodera Takatoshi Nakamura School of Fundamental Science and
More informationNewton-Krylov-Schwarz Method for a Spherical Shallow Water Model
Newton-Krylov-Schwarz Method for a Spherical Shallow Water Model Chao Yang 1 and Xiao-Chuan Cai 2 1 Institute of Software, Chinese Academy of Sciences, Beijing 100190, P. R. China, yang@mail.rdcps.ac.cn
More informationReactive transport in porous media: formulations, non-linear solvers and preconditioners
Reactive transport in porous media: formulations, non-linear solvers and preconditioners Michel Kern with L. Amir, B. Gueslin, A. Taakili Institut National de Recherche en Informatique et Automatique High
More informationChapter 5. Methods for Solving Elliptic Equations
Chapter 5. Methods for Solving Elliptic Equations References: Tannehill et al Section 4.3. Fulton et al (1986 MWR). Recommended reading: Chapter 7, Numerical Methods for Engineering Application. J. H.
More informationLesson Rigid Body Dynamics
Lesson 8 Rigid Body Dynamics Lesson 8 Outline Problem definition and motivations Dynamics of rigid bodies The equation of unconstrained motion (ODE) User and time control Demos / tools / libs Rigid Body
More informationAltered Jacobian Newton Iterative Method for Nonlinear Elliptic Problems
Altered Jacobian Newton Iterative Method for Nonlinear Elliptic Problems Sanjay K Khattri Abstract We present an Altered Jacobian Newton Iterative Method for solving nonlinear elliptic problems Effectiveness
More informationAn advanced ILU preconditioner for the incompressible Navier-Stokes equations
An advanced ILU preconditioner for the incompressible Navier-Stokes equations M. ur Rehman C. Vuik A. Segal Delft Institute of Applied Mathematics, TU delft The Netherlands Computational Methods with Applications,
More informationFAS and Solver Performance
FAS and Solver Performance Matthew Knepley Mathematics and Computer Science Division Argonne National Laboratory Fall AMS Central Section Meeting Chicago, IL Oct 05 06, 2007 M. Knepley (ANL) FAS AMS 07
More informationMultilevel Low-Rank Preconditioners Yousef Saad Department of Computer Science and Engineering University of Minnesota. Modelling 2014 June 2,2014
Multilevel Low-Rank Preconditioners Yousef Saad Department of Computer Science and Engineering University of Minnesota Modelling 24 June 2,24 Dedicated to Owe Axelsson at the occasion of his 8th birthday
More informationTrust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization
Trust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization Denis Ridzal Department of Computational and Applied Mathematics Rice University, Houston, Texas dridzal@caam.rice.edu
More informationCalculation of Sound Fields in Flowing Media Using CAPA and Diffpack
Calculation of Sound Fields in Flowing Media Using CAPA and Diffpack H. Landes 1, M. Kaltenbacher 2, W. Rathmann 3, F. Vogel 3 1 WisSoft, 2 Univ. Erlangen 3 inutech GmbH Outline Introduction Sound in Flowing
More informationHigh Performance Nonlinear Solvers
What is a nonlinear system? High Performance Nonlinear Solvers Michael McCourt Division Argonne National Laboratory IIT Meshfree Seminar September 19, 2011 Every nonlinear system of equations can be described
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationt x 0.25
Journal of ELECTRICAL ENGINEERING, VOL. 52, NO. /s, 2, 48{52 COMPARISON OF BROYDEN AND NEWTON METHODS FOR SOLVING NONLINEAR PARABOLIC EQUATIONS Ivan Cimrak If the time discretization of a nonlinear parabolic
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More informationX. He and M. Neytcheva. Preconditioning the incompressible Navier-Stokes equations with variable viscosity. J. Comput. Math., in press, 2012.
List of papers This thesis is based on the following papers, which are referred to in the text by their Roman numerals. I II III IV V VI M. Neytcheva, M. Do-Quang and X. He. Element-by-element Schur complement
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationProgress in Parallel Implicit Methods For Tokamak Edge Plasma Modeling
Progress in Parallel Implicit Methods For Tokamak Edge Plasma Modeling Michael McCourt 1,2,Lois Curfman McInnes 1 Hong Zhang 1,Ben Dudson 3,Sean Farley 1,4 Tom Rognlien 5, Maxim Umansky 5 Argonne National
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationDepartment of Computer Science, University of Illinois at Urbana-Champaign
Department of Computer Science, University of Illinois at Urbana-Champaign Probing for Schur Complements and Preconditioning Generalized Saddle-Point Problems Eric de Sturler, sturler@cs.uiuc.edu, http://www-faculty.cs.uiuc.edu/~sturler
More informationA robust multilevel approximate inverse preconditioner for symmetric positive definite matrices
DICEA DEPARTMENT OF CIVIL, ENVIRONMENTAL AND ARCHITECTURAL ENGINEERING PhD SCHOOL CIVIL AND ENVIRONMENTAL ENGINEERING SCIENCES XXX CYCLE A robust multilevel approximate inverse preconditioner for symmetric
More information