A Nodal High-Order Discontinuous Galerkin Dynamical Core for Climate Simulations
|
|
- Clementine Edwards
- 5 years ago
- Views:
Transcription
1 A Nodal High-Order Discontinuous Galerkin Dynamical Core for Climate Simulations Institute for Mathematics Applied to the Geosciences (IMAGe) National Center for Atmospheric Research (NCAR) Boulder CO 80305, USA. International Conference On Spectral And High Order Methods June 18-22, 2007 Chinese Academy of Sciences, Beijing, China.
2 HOMME: High-Order Method Modeling Environment Supported by the DoE SciDAC and CCPP programs and NSF sponsorship of NCAR HOMME Developers Henry Tufo (PI, CU Boulder & NCAR ) John Dennis, Amik St-Cyr, Joe Tribbia & Ram Nair (Co-Is, NCAR) Jim Edward (IBM/NCAR), Hae-Won Choi (Postdoc, CU) Mike Levy and Theron Voran (Grad Students), Mark Taylor (SA Collaborator, Sandia National Laboratories) Project URL:
3 HOMME Motivation The DG Baroclinic Model (HOMME) 1 Vertical aspects (Lagrangian Dynamics, Remapping) 2 Horizontal Aspects (DGM, Discretization) Validation 1 2D 2 3D Parallel Performance Some Computational Issues Summary
4 Motivation HOMME Why do we need a new climate model? Because, the existing models have serious limitations to satisfy all of the following properties: 1 Local and global conservation 2 High-order accuracy 3 High parallel efficiency 4 Geometric flexibility ( Local method) 5 Monotonic (non-oscillatory) advection Discontinuous Galerkin Method (DGM) based model has the potential to address all of the above issues Recently, the Spectral Element (SE) model in HOMME shown to efficiently scale O (10, 000) processors on IBM BG/L.
5 HOMME (High-Order Method Modeling Environment) The Discontinuous Galerkin (DG) model is the conservative option in the HOMME framework HOMME Grid: The sphere is decomposed into 6 identical regions, using the equiangular projection (Sadourny, 1972) Local coordinate systems are free of singularities Creates a non-orthogonal curvilinear coordinate system Cubed Sphere Geometry: Logical cube-face orientation F 5 (Top) z F4 F 1 F 2 F 3 F 1 y F 6 (Bottom) x
6 Grid System HOMME Metric Tensor G ij, [Cubed-Sphere Sphere] Transform Central angles x 1,x 2 [ π/4,π/4] are the independent variables. G ij = R 2 ρ 4 cos 2 x 1 cos 2 x 2 [ 1 + tan 2 x 1 tanx 1 tanx 2 tanx 1 tanx tan 2 x 2 where ρ 2 = 1 + tan 2 x 1 + tan 2 x 2, i,j {1,2} ] Metric tensor in terms of longitude-latitude (λ,θ): [ G ij = A T R cos θ λ/ x 1 R cos θ λ/ x A; A = 2 R θ/ x 1 R θ/ x 2 ] The matrix A is used for transforming spherical velocity (u,v) to the covariant (u 1,u 2 ) and contravariant (u 1,u 2 )
7 3D Prognostic Equations for the DG Model Momentum Equations: No explicit vertical advection terms u 1 t + c E 1 u 2 t + c E 2 = Gu 2 (f + ζ) R T x 1(lnp) = Gu 1 (f + ζ) R T x2(ln p) c ( ) x 1, x 2, E 1 = (E,0), E 2 = (0,E), E = Φ + 1 ( u1 u 1 + u 2 u 2) 2 G = det(gij ), E is the energy term, Φ is the geopotential, and ζ is the relative vorticity.
8 3D Prognostic Equations: Flux-Form Continuity Equations Temperature field is advected with the mass variable p t ( p) + c (U i p ) = 0 t (Θ p) + c (U i Θ p ) = 0 t (q p) + c (U i q p ) = 0 where U i = ( u 1, u 2), p = Gδp, δp is the pressure thickness, and Θ is the potential temperature. Vertical layers are coupled with the hydrostatic relations: Φ = C p Θ Π, Φ = RT lnp where (p/p 0 ) κ and T Denotes the layer mean temperature.
9 Vertical Lagrangian Coordinates Vertically moving Lagrangian Surfaces The velocity fields (u 1,u 2 ), and total energy (Γ E ) are remapped onto the reference coordinates using the 1-D conservative cell-integrated semi-lagrangian (CISL) method (Nair & Machenhauer, 2002) P= Pressure thickness Lagrangian Surface t + t t P L 2 P E2 L 1 E 1 Topography Terrain following Lagrangian control volume coordinates
10 Computational Grid Structure for DG Model k 1/2 ( φ,p) k 1/2 k ( u, v, θ, q, δp ) k k+1/2 ( φ, p) k+1/2 GLL Grid box The remapping frequency is O(10) t Potential temperature Θ is retrieved from the remapped total energy Γ E = c p T + δ(pφ) δp + K E
11 The Vertical Lagrangian Dynamics The hydrostatic pressure at Lagrangian surface, Lin (MWR, 2004) p l = p top + l δp k, l = 1,2,3,...,N k=1 where p top represents the pressure at the model top, p l denotes the pressure at each Lagrangian surface. There are total N + 1 Lagrangian surfaces span N layers. The geopotential height at Lagrangian surface: Φ l = Φ s + l Φ k, l = N,N 1,...,1 k=n where Φ s represents the surface geopotential height at the model bottom and Φ l denotes the geopotential height at each Lagrangian surface.
12 DG-3D Model Horizontal Aspects: Shallow Water Model Flux-form SW equations (Vector invariant form) Nair et al. (MWR, 2005) u 1 t + x 1E = G u 2 (f + ζ) u 2 t + x 2E = G u 1 (f + ζ) t ( G h) + x 1( G u 1 h) + x 2( G u 2 h) = 0 where G = det(g ij ), h is the height, f Coriolis term; energy term and vorticity are defined as E = Φ (u 1 u 1 + u 2 u 2 ), ζ = 1 [ u2 G x 1 u ] 1 x 2.
13 The DG & SE Methods in HOMME DG SE Boundary Discontinuity FV 0 C Continuous DGM is a hybrid approach (DG SE + FV) The domain D is partitioned into non-overlapping elements Ω ij such that the element boundaries are discontinuous. Based on conservation laws but exploits the spectral expansion of SE method and treats the element boundaries using FV tricks.
14 DG Spatial Discretization for an Element Ω SW Equations can be written as 2D Scalar conservation law U t + F(U) = S(U), in Ω (0,T); (x1,x 2 ) Ω where U = U(x 1,x 2,t), ( / x 1, / x 2 ), F = (F,G) is the flux function, and S is the source term. Weak Galerkin formulation: Multiplication of the basic equation by a test function ϕ h V h and integration over an element Ω. Z Z Z Z U h ϕ h dω F(U h ) ϕ h dω + F(U h ) n ϕ h dγ = S(U h ) ϕ h dω t Ω Ω Γ Ω where U h is an approximate solution in Ω.
15 DGM: The Flux Term HOMME Element (Left) Element (Right) Element (Left) Element (Right) _ U h U + h U h U h After Num. Flux Operation Along the boundaries (Γ) of an element the solution U h is discontinuous (U h and U + h are the left and right limits). Therefore, the analytic flux F(U h ) n must be replaced by a numerical flux such as the Lax-Friedrichs Flux: F(U h ) n = 1 i h(f(u h 2 ) + F(U+ h )) n α(u+ h U h ). For the SW system, α is the upper bound on the absolute value of eigenvalues of the flux Jacobian F (U); (Nair et al., 2005) α 1 = max u 1 + Φ G 11, α 2 = max u 2 + Φ G 22
16 DG-3D Model: Computational Domain ( 1, 1) GLL Grid (8x8) (1, 1) η ( 1, 1) Cubed-Sphere (N e = 5) with 8 8 GLL points ξ (1, 1) Flux is the only communicator at the element edges Each face of the cubed-sphere is partitioned into N e N e rectangular non-overlapping elements (i.e., total 6 N 2 e ). Each element is mapped onto the Gauss-Lobatto-Legendre (GLL) grid defined by 1 ξ,η 1, for integration.
17 DGM: Spatial Discretization High-order nodal basis set The nodal basis set is constructed using a tensor-product of Lagrange-Legendre polynomials (h i (ξ)) with roots at Gauss-Lobatto quadrature points. (ξ 2 1)P N h i (ξ) = (ξ) 1 N(N + 1)P N (ξ i )(ξ ξ i ) ; h i (ξ)h j (ξ) = w i δ ij. where P N (ξ) is the N th order Legendre polynomial, and w i weights associated with the Gauss quadrature. The approximate solution (U h ) and test function (ϕ h ) are represented in terms of nodal basis set. U ij (ξ, η) = NX i=0 j=0 1 NX U ij h i (ξ)h j (η) for 1 ξ, η 1, The nodal version was shown to be more computationally efficient than the modal in (Dennis et al., 2006).
18 DG-3D Model: Explicit Time Integration Final form for the nodal discretization d dt U 4 ij(t) = xi 1 [I Grad + I Flux + I Source ], x2 j w i w j Solve the ODE system d U = L(U) dt in (0,T) Ω Time integration: Explicit third-order Runge-Kutta (SSP) scheme (Gottlieb et al., 2001) Options for explicit diffusion ( 2 or 4 ). Boyd-Vandeevan spatial Filter
19 2D results Horizontal Advection: Moving Vortices on the Sphere Initial field and DG solution after 12 days. Max error is O(10 5 ) Deformational Flow Test: Nair & Jablonowski (MWR, 2007) The vortices are located at diametrically opposite sides of the sphere, the vortices deform as they move along a prescribed trajectory. Analytical solution is known and the trajectory is chosen to be a great circle along the NE direction (α = π/4).
20 2D results SW Test-2: Geostrophic Flow (Nair et al. 2005) High-order accuracy and spectral convergence Steady state geostrophic flow (α = π/4). Max height error is O(10 6 ) m.
21 2D results SW Test-5: Flow over a Mountain (Dennis et al. 2006) No spectral ringing for the height fields Flow over a mountain ( 0.5 o ). Initial height field (left) initial and after 15 days of integration (right)
22 DG-3D: Baroclinic Instability Test (JW-Test) Jablonowski & Williamson (QJRMS, 2006) To assess the evolution of an idealized baroclinic wave in the Northern Hemisphere. The initial conditions are quasi-realistic and defined by analytic expressions. Analytic solutions do not exist. Initial Conditions
23 JW-Test: Evolution of Surface Pressure over the NH Baroclinic waves are triggered by perturbing the velocity field at (20 E, 40 N) This test case recommends up to 30 days of model simulation Ne = Nv = 8 (approx. 1.6 ) with 26 vertical levels and t = 30 Sec.
24 DG-3D Model Vs. NCAR Spectral Model The HOMME-DG dynamical core successfully simulates baroclinic instability. Simulated temperature (K) and surface pressure (hpa) at day 8 for a baroclinic instability test with the HOMME-DG model and the NCAR global spectral model (right). The horizontal resolution is approximately 1.4. Note that the DG solution is free of spectral ringing.
25 Jablonowski-Williamson Baroclinic Test (Convergence)
26 DG Model Vs. NCAR Climate Models Simulated surface pressure at day 11 for a baroclinic instability test with DG model, and NCAR global spectral model and a FV model. The models use 26 vertical levels and with approximate horizontal resolution of 0.7.
27 Parallel Performance (3D) - Frost [IBM BG/L] DG-3D parallel performance: Sustained Mflops on IBM BG/L (1024 DP nodes, 700 MHz PPC 440s): Approx. 9% peak Held-Suarez (preliminary) test: 800 days idealized climate simulation (1 resolution, 26 vertical levels, t = 10 Sec) 300 Sustaind FLOP Per Processor MFLOP elements: 1 task/node (CO) elements: 2 task/node (VN) 7776 elements: 1 tasks/node (CO) 7776 elements: 2 tasks/node (VN) Processor Parallel performance (strong scaling) results for JW-Test Held-Suarez test (800 days).
28 Some computational issues with limiting (Levy et al. 2007) Monotonic limiters do exist for relatively low-order DGM (N 4). The WENO based limiters (Qiu & Shu, 2005) are not only computationally expensive but results in a parallel communication bottleneck. Extending to high-order (N > 4) is impractical. In HOMME N > DG advection of a non-smooth Doswell vortex combined with a third-order WENO limiter (40 40) elements, limiting an element requires a 5 5 stencil!
29 Summary HOMME The DG-3D model successfully simulates the Baroclinic instability test and the results are comparable with that of the NCAR global spectral model. The preliminary scaling results are impressive and comparable to the SE version in HOMME. The R-K time integration scheme is robust for the DG-3D model, but very time-step restrictive. More efficient time integration schemes are required for practical climate simulations (semi-implicit). Efficient monotonic limiter for advection? Future Work: Time integration schemes. Coupling of the CAM/CCSM physics for the real climate simulations. Targeting for large-scale parallelism with O(100K) processors.
Cubed-Sphere Grids and Galerkin Approaches
2008 Global Atmospheric Core Workshop National Center for Atmospheric Research (NCAR) September 23-24, Boulder CO, USA. Institute for Mathematics Applied to the Geosciences (IMAGe) Cubed-Sphere: Central
More informationDevelopment of a Scalable Global Discontinuous Galerkin Atmospheric Model
Development of a Scalable Global Discontinuous Galerkin Atmospheric Model J.M. Dennis, R.D. Nair and H.M. Tufo Scientific Computing Division, National Center for Atmospheric Research 1850 Table Mesa Dr.,
More informationA Discontinuous Galerkin Global Shallow Water Model
876 M O N T H L Y W E A T H E R R E V I E W VOLUME 133 A Discontinuous Galerkin Global Shallow Water Model RAMACHANDRAN D. NAIR, STEPHEN J. THOMAS, AND RICHARD D. LOFT Scientific Computing Division, National
More informationA spectral finite volume transport scheme on the cubed-sphere
Applied Numerical Mathematics 57 (2007) 1021 1032 www.elsevier.com/locate/apnum A spectral finite volume transport scheme on the cubed-sphere Vani Cheruvu a,, Ramachandran D. Nair b, Henry M. Tufo a,b
More informationA High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere
A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere Michael Levy University of Colorado at Boulder Department of Applied Mathematics August 10, 2007 Outline 1 Background
More informationA Global Atmospheric Model. Joe Tribbia NCAR Turbulence Summer School July 2008
A Global Atmospheric Model Joe Tribbia NCAR Turbulence Summer School July 2008 Outline Broad overview of what is in a global climate/weather model of the atmosphere Spectral dynamical core Some results-climate
More informationNCAR Global Atmospheric Core Workshop, Boulder, June 2008
NCAR Global Atmospheric Core Workshop, Boulder, June 2008 D. Majewski based on Christiane Jablonowski (email: cjablono@umich.edu) University of Michigan Goals of the Test Suite NASA/GFDL Test cases should
More informationEmerging Numerical Methods for Atmospheric Modeling
Chapter 9 Emerging Numerical Methods for Atmospheric Modeling Ramachandran D. Nair, Michael N. Levy and Peter H. Lauritzen Abstract This chapter discusses the development of discontinuous Galerkin (DG)
More informationAMERICAN METEOROLOGICAL SOCIETY
AMERICAN METEOROLOGICAL SOCIETY Monthly Weather Review EARLY ONLINE RELEASE This is a preliminary PDF of the author-produced manuscript that has been peer-reviewed and accepted for publication. Since it
More informationA High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations
Motivation Numerical methods Numerical tests Conclusions A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations Xiaofeng Cai Department of Mathematics
More informationOverview of the Numerics of the ECMWF. Atmospheric Forecast Model
Overview of the Numerics of the Atmospheric Forecast Model M. Hortal Seminar 6 Sept 2004 Slide 1 Characteristics of the model Hydrostatic shallow-atmosphere approimation Pressure-based hybrid vertical
More informationThe Shallow Water Equations
If you have not already done so, you are strongly encouraged to read the companion file on the non-divergent barotropic vorticity equation, before proceeding to this shallow water case. We do not repeat
More informationAtmospheric Dynamics with Polyharmonic Spline RBFs
Photos placed in horizontal position with even amount of white space between photos and header Atmospheric Dynamics with Polyharmonic Spline RBFs Greg Barnett Sandia National Laboratories is a multimission
More informationA Strategy for the Development of Coupled Ocean- Atmosphere Discontinuous Galerkin Models
A Strategy for the Development of Coupled Ocean- Atmosphere Discontinuous Galerkin Models Frank Giraldo Department of Applied Math, Naval Postgraduate School, Monterey CA 93943 Collaborators: Jim Kelly
More informationAtmospheric Dynamics
Atmospheric Dynamics the CAM-FV dynamical core Peter Hjort Lauritzen Atmospheric Modeling and Predictability Section Climate and Global Dynamics Division National Center for Atmospheric Research Community
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 informationA semi-implicit non-hydrostatic covariant dynamical kernel using spectral representation in the horizontal and a height based vertical coordinate
A semi-implicit non-hydrostatic covariant dynamical kernel using spectral representation in the horizontal and a height based vertical coordinate Juan Simarro and Mariano Hortal AEMET Agencia Estatal de
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Dept. of
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
More informationA Class of Deformational Flow Test Cases for Linear Transport Problems on the Sphere
A Class of Deformational Flow Test Cases for Linear Transport Problems on the Sphere Ramachandran D. Nair,a, Peter H. Lauritzen b a Institute for Mathematics Applied to Geosciences (IMAGe), National Center
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationConservation Laws & Applications
Rocky Mountain Mathematics Consortium Summer School Conservation Laws & Applications Lecture V: Discontinuous Galerkin Methods James A. Rossmanith Department of Mathematics University of Wisconsin Madison
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationFormulation and performance of the Variable-Cubic Atmospheric Model
Formulation and performance of the Variable-Cubic Atmospheric Model John McGregor CSIRO Marine and Atmospheric Research Aspendale, Melbourne Southern Hemisphere PDEs on the Sphere NCAR 11 April 2014 CSIRO
More informationICON-IAP: A non-hydrostatic global. Almut Gassmann IAP Kühlungsborn, Germany. model designed for energetic consistency
ICON-IAP: A non-hydrostatic global Almut Gassmann IAP Kühlungsborn, Germany model designed for energetic consistency ICON-IAP is a model on a Voronoi mesh: The divergence lives on Voronoi cells (mostly
More informationCURRICULUM VITAE. Ramachandran D. Nair
CURRICULUM VITAE Ramachandran D. Nair Contact Information Mailing Address: Institute for Mathematics Applied to the Geosciences National Center for Atmospheric Research (NCAR) 1850 Table Mesa Drive, Boulder,
More informationMixed Mimetic Spectral Elements for Geophysical Fluid Dynamics
for Geophysical Fluid Dynamics Dave Lee Los Alamos National Laboratory Outline Connection of finite volumes to differential forms Key ideas of differential forms Differential forms for discrete data Construction
More informationPositivity-preserving high order schemes for convection dominated equations
Positivity-preserving high order schemes for convection dominated equations Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Xiangxiong Zhang; Yinhua Xia; Yulong Xing; Cheng
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationDevelopment of Yin-Yang Grid Global Model Using a New Dynamical Core ASUCA.
Development of Yin-Yang Grid Global Model Using a New Dynamical Core ASUCA. M. Sakamoto, J. Ishida, K. Kawano, K. Matsubayashi, K. Aranami, T. Hara, H. Kusabiraki, C. Muroi, Y. Kitamura Japan Meteorological
More informationChallenges and Opportunities in Modeling of the Global Atmosphere
Challenges and Opportunities in Modeling of the Global Atmosphere Zavisa Janjic, Vladimir Djurdjevic, Ratko Vasic, and Tom Black NOAA/NCEP/EMC, College Park, MD 20740, U.S.A. (Email: zavisa.janjic@noaa.gov)
More informationHYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS
1 / 36 HYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS Jesús Garicano Mena, E. Valero Sánchez, G. Rubio Calzado, E. Ferrer Vaccarezza Universidad
More informationSpectral transforms. Contents. ARPEGE-Climat Version 5.1. September Introduction 2
Spectral transforms ARPEGE-Climat Version 5.1 September 2008 Contents 1 Introduction 2 2 Spectral representation 2 2.1 Spherical harmonics....................... 2 2.2 Collocation grid.........................
More informationA high-order fully explicit incremental-remap -based semi-lagrangian shallow-water model
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 0000; 00:1 33 Published online in Wiley InterScience (www.interscience.wiley.com). A high-order fully explicit incremental-remap
More informationA primal-dual mixed finite element method. for accurate and efficient atmospheric. modelling on massively parallel computers
A primal-dual mixed finite element method for accurate and efficient atmospheric modelling on massively parallel computers John Thuburn (University of Exeter, UK) Colin Cotter (Imperial College, UK) AMMW03,
More informationHybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation
Numerical Analysis and Scientific Computing Preprint Seria Hybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation W. Guo J. Qiu Preprint #21 Department of Mathematics University
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationDynamical Core Developments at KIAPS
Dynamical Core Developments at KIAPS 1/ 48 Dynamical Core Developments at KIAPS Preliminary results T.-J. Oh, S.-J. Choi, T.-H. Yi, S.-H. Kang, and Y.-J. Kim Korea Institute of Atmospheric Prediction Systems
More informationUNCORRECTED PROOF. A Fully Implicit Compressible Euler Solver for 2 Atmospheric Flows * 3. 1 Introduction 9. Chao Yang 1,2 and Xiao-Chuan Cai 2 4
1 A Fully Implicit Compressible Euler Solver for 2 Atmospheric Flows * 3 Chao Yang 1,2 and Xiao-Chuan Cai 2 4 1 Institute of Software, Chinese Academy of Sciences, Beijing 100190, P. R. China, 5 yangchao@iscas.ac.cn
More informationAcceleration of Time Integration
Acceleration of Time Integration edited version, with extra images removed Rick Archibald, John Drake, Kate Evans, Doug Kothe, Trey White, Pat Worley Research sponsored by the Laboratory Directed Research
More informationCirculation and Vorticity
Circulation and Vorticity Example: Rotation in the atmosphere water vapor satellite animation Circulation a macroscopic measure of rotation for a finite area of a fluid Vorticity a microscopic measure
More informationBound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu
Bound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu Division of Applied Mathematics Brown University Outline Introduction Maximum-principle-preserving for scalar conservation
More informationJupiter s Red Spot. Luke McGuire Program in Applied Mathematics University of Arizona. May 21, 2009
Jupiter s Red Spot Luke McGuire Program in Applied Mathematics University of Arizona May 21, 2009 Abstract The Great Red Spot, the largest of Jupiter s long lived, approximately two dimensional, coherent
More informationWeighted ENO Schemes
Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University, The State University of New York February 7, 014 1 3 Mapped WENO-Z Scheme 1D Scalar Hyperbolic
More informationCESM2 release of CAM-SE (& CAM-SE-CSLAM)
CESM2 release of CAM-SE (& CAM-SE-CSLAM) Peter Hjort Lauritzen National Center for Atmospheric Research Boulder, Colorado, USA Collaborators: S. Goldhaber (NCAR), J. Bacmeister (NCAR), R.D.Nair (NCAR),
More informationEarth System Modeling Domain decomposition
Earth System Modeling Domain decomposition Graziano Giuliani International Centre for Theorethical Physics Earth System Physics Section Advanced School on Regional Climate Modeling over South America February
More informationFast Numerical Methods for Stochastic Computations
Fast AreviewbyDongbinXiu May 16 th,2013 Outline Motivation 1 Motivation 2 3 4 5 Example: Burgers Equation Let us consider the Burger s equation: u t + uu x = νu xx, x [ 1, 1] u( 1) =1 u(1) = 1 Example:
More informationPetascale atmospheric models for the Community Climate System Model: new developments and evaluation of scalable dynamical cores
Journal of Physics: Conference Series Petascale atmospheric models for the Community Climate System Model: new developments and evaluation of scalable dynamical cores To cite this article: M A Taylor et
More informationA parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows.
A parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows Tao Xiong Jing-ei Qiu Zhengfu Xu 3 Abstract In Xu [] a class of
More informationSpectral element schemes for the. Korteweg-de Vries and Saint-Venant equations
Spectral element schemes for the Korteweg-de Vries and Saint-Venant equations R. PASQUETTI a a. Université Côte d Azur, CNRS, Inria, LJAD, France, richard.pasquetti@unice.fr... Résumé : Les sytèmes hyperboliques
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
More informationYINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationOn simplifying incremental remap -based transport schemes
On simplifying incremental remap -based transport schemes Peter H. Lauritzen a,, Christoph Erath b, Rashmi Mittal c a Climate and Global Dynamics Division, National Center for Atmospheric Research, P.O.
More informationA second-order asymptotic-preserving and positive-preserving discontinuous. Galerkin scheme for the Kerr-Debye model. Abstract
A second-order asymptotic-preserving and positive-preserving discontinuous Galerkin scheme for the Kerr-Debye model Juntao Huang 1 and Chi-Wang Shu Abstract In this paper, we develop a second-order asymptotic-preserving
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationQuasi-3D Multiscale Modeling Framework as a Physics Option in CAM-SE
Quasi-3D Multiscale Modeling Framework as a Physics Option in CAM-SE Joon-Hee Jung, Celal S. Konor, Don Dazlich, David A. Randall, Department of Atmospheric Sciences, Colorado State University, USA Peter
More informationA class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes
Science in China Series A: Mathematics Aug., 008, Vol. 51, No. 8, 1549 1560 www.scichina.com math.scichina.com www.springerlink.com A class of the fourth order finite volume Hermite weighted essentially
More information1/3/2011. This course discusses the physical laws that govern atmosphere/ocean motions.
Lecture 1: Introduction and Review Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study of motion without regard for the cause. Dynamics: On the other hand, dynamics
More informationA Proposed Test Suite for Atmospheric Model Dynamical Cores
A Proposed Test Suite for Atmospheric Model Dynamical Cores Christiane Jablonowski (cjablono@umich.edu) University of Michigan, Ann Arbor PDEs on the Sphere Workshop Monterey, CA, June/26-29/2006 Motivation
More informationImproving Dynamical Core Scalability, Accuracy, and Limi:ng Flexibility with the ADER- DT Time Discre:za:on
Improving Dynamical Core Scalability, Accuracy, and Limi:ng Flexibility with the ADER- DT Time Discre:za:on Matthew R. Norman Scientific Computing Group National Center for Computational Sciences Oak Ridge
More informationA positivity-preserving high order discontinuous Galerkin scheme for convection-diffusion equations
A positivity-preserving high order discontinuous Galerkin scheme for convection-diffusion equations Sashank Srinivasan a, Jonathan Poggie a, Xiangxiong Zhang b, a School of Aeronautics and Astronautics,
More informationGeneralised Summation-by-Parts Operators and Variable Coefficients
Institute Computational Mathematics Generalised Summation-by-Parts Operators and Variable Coefficients arxiv:1705.10541v [math.na] 16 Feb 018 Hendrik Ranocha 14th November 017 High-order methods for conservation
More informationA Robust CFL Condition for the Discontinuous Galerkin Method on Triangular Meshes
A Robust CFL Condition for the Discontinuous Galerkin Method on Triangular Meshes N. Chalmers and L. Krivodonova Abstract When the discontinuous Galerkin (DG) method is applied to hyperbolic problems in
More informationSIMULATION OF GEOPHYSICAL TURBULENCE: DYNAMIC GRID DEFORMATION
SIMULATION OF GEOPHYSICAL TURBULENCE: DYNAMIC GRID DEFORMATION Piotr K Smolarkiewicz, and Joseph M Prusa, National Center for Atmospheric Research, Boulder, Colorado, U.S.A. Iowa State University, Ames,
More informationRecent developments of the global semi-lagrangian atmospheric model
Recent developments of the global semi-lagrangian atmospheric model Mikhail Tolstykh, Vladimir Shashkin (Inst. of Numerical Math. RAS, Hydrometcentre of Russia), Alla Yurova (Hydrometcentre of Russia)
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationThe Advanced Research WRF (ARW) Dynamics Solver
Dynamics: Introduction The Advanced Research WRF (ARW) Dynamics Solver 1. What is a dynamics solver? 2. Variables and coordinates 3. Equations 4. Time integration scheme 5. Grid staggering 6. Advection
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationA stable treatment of conservative thermodynamic variables for semi-implicit semi-lagrangian dynamical cores
A stable treatment of conservative thermodynamic variables for semi-implicit semi-lagrangian dynamical cores Kevin Viner Naval Research Laboratory, Monterey, CA September 26, 2012 Kevin Viner (NRL) PDE
More informationNew variables in spherical geometry. David G. Dritschel. Mathematical Institute University of St Andrews.
New variables in spherical geometry David G Dritschel Mathematical Institute University of St Andrews http://www-vortexmcsst-andacuk Collaborators: Ali Mohebalhojeh (Tehran St Andrews) Jemma Shipton &
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationScalable Non-Linear Compact Schemes
Scalable Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory International Conference on Spectral and High Order Methods
More informationIMPACTS OF SIGMA COORDINATES ON THE EULER AND NAVIER-STOKES EQUATIONS USING CONTINUOUS/DISCONTINUOUS GALERKIN METHODS
Approved for public release; distribution is unlimited IMPACTS OF SIGMA COORDINATES ON THE EULER AND NAVIER-STOKES EQUATIONS USING CONTINUOUS/DISCONTINUOUS GALERKIN METHODS Sean L. Gibbons Captain, United
More informationThe Lattice Boltzmann method for hyperbolic systems. Benjamin Graille. October 19, 2016
The Lattice Boltzmann method for hyperbolic systems Benjamin Graille October 19, 2016 Framework The Lattice Boltzmann method 1 Description of the lattice Boltzmann method Link with the kinetic theory Classical
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationA Stable Spectral Difference Method for Triangles
A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna
More informationHigh order finite difference Hermite WENO schemes for the Hamilton-Jacobi equations on unstructured meshes
High order finite difference Hermite WENO schemes for the Hamilton-Jacobi equations on unstructured meshes Feng Zheng, Chi-Wang Shu and Jianxian Qiu 3 Abstract In this paper, a new type of high order Hermite
More informationLecture 2: Reconstruction and decomposition of vector fields on the sphere with applications
2013 Dolomites Research Week on Approximation : Reconstruction and decomposition of vector fields on the sphere with applications Grady B. Wright Boise State University What's the problem with vector fields
More informationPOD Model Reduction of Large Scale Geophysical Models. Ionel M. Navon
POD Model Reduction of Large Scale Geophysical Models Ionel M. Navon School of Computational Science Florida State University Tallahassee, Florida Thanks to Prof. D. N. Daescu Dept. of Mathematics and
More informationCENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationBackground Error Covariance Modelling
Background Error Covariance Modelling Mike Fisher Slide 1 Outline Diagnosing the Statistics of Background Error using Ensembles of Analyses Modelling the Statistics in Spectral Space - Relaxing constraints
More informationMathematical Concepts & Notation
Mathematical Concepts & Notation Appendix A: Notation x, δx: a small change in x t : the partial derivative with respect to t holding the other variables fixed d : the time derivative of a quantity that
More informationCase Study: Numerical Simulations of the Process by which Planets generate their Magnetic Fields. Stephan Stellmach UCSC
Case Study: Numerical Simulations of the Process by which Planets generate their Magnetic Fields Stephan Stellmach UCSC 1 The Scientific Problem 2 The Magnetic Field of the Earth http://en.wikipedia.org/wiki/image:magnetosphere_rendition.jpg
More informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationOutline. 1 Boundary Value Problems. 2 Numerical Methods for BVPs. Boundary Value Problems Numerical Methods for BVPs
Boundary Value Problems Numerical Methods for BVPs Outline Boundary Value Problems 2 Numerical Methods for BVPs Michael T. Heath Scientific Computing 2 / 45 Boundary Value Problems Numerical Methods for
More informationarxiv: v2 [math.na] 9 Apr 2018
oname manuscript o. will be inserted by the editor) The BR1 Scheme is Stable for the Compressible avier-stokes Equations Gregor J. Gassner Andrew R. Winters Florian J. Hindenlang David A. Kopriva arxiv:1704.03646v
More informationLocally Linearized Euler Equations in Discontinuous Galerkin with Legendre Polynomials
Locally Linearized Euler Equations in Discontinuous Galerkin with Legendre Polynomials Harald Klimach, Michael Gaida, Sabine Roller harald.klimach@uni-siegen.de 26th WSSP 2017 Motivation Fluid-Dynamic
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationLinear model for investigation of nonlinear NWP model accuracy. Marko Zirk, University of Tartu
Linear model for investigation of nonlinear NWP model accuracy Marko Zirk, University of Tartu Introduction A method for finding numerical solution of non-hydrostatic linear equations of atmospheric dynamics
More informationPAPER 333 FLUID DYNAMICS OF CLIMATE
MATHEMATICAL TRIPOS Part III Wednesday, 1 June, 2016 1:30 pm to 4:30 pm Draft 21 June, 2016 PAPER 333 FLUID DYNAMICS OF CLIMATE Attempt no more than THREE questions. There are FOUR questions in total.
More informationCAM-SE: Lecture I. Peter Hjort Lauritzen
CAM-SE: Lecture I Peter Hjort Lauritzen Atmospheric Modeling and Predictability Section Climate and Global Dynamics Laboratory National Center for Atmospheric Research 2nd WCRP Summer School on Climate
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationThe Spectral-Element Method: Introduction
The Spectral-Element Method: Introduction Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Computational Seismology 1 / 59 Outline 1 Introduction 2 Lagrange
More informationExperience with DNS of particulate flow using a variant of the immersed boundary method
Experience with DNS of particulate flow using a variant of the immersed boundary method Markus Uhlmann Numerical Simulation and Modeling Unit CIEMAT Madrid, Spain ECCOMAS CFD 2006 Motivation wide range
More informationAICS TECHNICAL REPORT
AICS TECHNICAL REPORT N0. 2017-001 Modification of a baroclinic test case to a height-based coordinate system without topography implementation By Ryuji Yoshida 1,2, Hisashi Yashiro 1, Masaki Satoh 3,
More informationDust devils, water spouts, tornados
Balanced flow Things we know Primitive equations are very comprehensive, but there may be a number of vast simplifications that may be relevant (e.g., geostrophic balance). Seems that there are things
More informationAtmospheric modeling in the Climate System. Joe Tribbia NCAR.ESSL.CGD.AMP
Atmospheric modeling in the Climate System Joe Tribbia NCAR.ESSL.CGD.AMP The climate represents a coupled system consisting of an atmosphere, hydrosphere, biosphere, and cryosphere What is CCSM? Bio Geochemistry
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More information