Investigating Castillo-Grone's Mimetic Difference Operators in Development of Geophysical Fluid Dynamics Models Implemented on GPGPUs

of the Dissertation Investigating Castillo-Grone’s Mimetic Difference Operators in Development of Geophysical Fluid Dynamics Models Implemented on GPGPUs by Mohammad Abouali San Diego State University and Claremont Graduate University This research investigates the performance of Castillo-Grone’s Mimetic (CGM) difference operators in 2D and 3D fully curvilinear grids. A software package is provided that assists the user in generating Laplacian operator in fully curvilinear grids using higher order CGM operators. Furthermore, a stability analysis is performed using von Neumann’s analysis to understand the behavior of these operators in solving advection equation. A shallow water equation solver is provided that harnesses the power of the general purpose graphics processing units (GPGPUs) and some timing are provided.

[1]  J. E. Castillo,et al.  A Matrix Analysis Approach to Higher-Order Approximations for Divergence and Gradients Satisfying a Global Conservation Law , 2003, SIAM J. Matrix Anal. Appl..

[2]  T. Coakley,et al.  Turbulence Modeling Validation, Testing, and Development , 1997 .

[3]  D. Wilcox.  Reassessment of the scale-determining equation for advanced turbulence models , 1988 .

[4]  Louis J. Wicker,et al.  Time-Splitting Methods for Elastic Models Using Forward Time Schemes , 2002 .

[5]  Patrick M. Knupp,et al.  Fundamentals of Grid Generation , 2020 .

[6]  Greg Holloway,et al.  Representing topographic stress for large-scale ocean models , 1992 .

[7]  Alexander Kurganov,et al.  Semidiscrete Central-Upwind Schemes for Hyperbolic Conservation Laws and Hamilton-Jacobi Equations , 2001, SIAM J. Sci. Comput..

[8]  Charles Meneveau,et al.  A scale-dependent Lagrangian dynamic model for large eddy simulation of complex turbulent flows , 2005 .

[9]  G. Mellor,et al.  Development of a turbulence closure model for geophysical fluid problems , 1982 .

[10]  J. Berntsen.  Internal Pressure Errors in Sigma-Coordinate Ocean Models , 2002 .

[11]  Malcolm J. Bowman,et al.  Rim current and coastal eddy mechanisms in an eddy-resolving Black Sea general circulation model , 2001 .

[12]  A. Thurnherr,et al.  Boundary Mixing and Topographic Blocking on the Mid-Atlantic Ridge in the South Atlantic* , 2003 .

[13]  James C. McWilliams.  Modeling the Oceanic General Circulation , 1996 .

[14]  Mikhail Shashkov,et al.  The sensitivity and accuracy of fourth order finite-difference schemes on nonuniform grids in one dimension , 1995 .

[15]  M. Merrifield,et al.  Model Estimates of M2 Internal Tide Generation over Mid-Atlantic Ridge Topography , 2009 .

[16]  J. M. Toole,et al.  Spatial Variability of Turbulent Mixing in the Abyssal Ocean , 1997, Science.

[17]  R. Inghilesi,et al.  Large Eddy Simulation of Turbulent Mixing in an Estuary Region , 2010 .

[18]  F. Menter.  Two-equation eddy-viscosity turbulence models for engineering applications , 1994 .

[19]  J. Castillo,et al.  Steady state diffusion problems on non-trivial domains: support operator method integrated with direct optimized grid generation , 2002 .

[20]  A. Mestas-Nuñez,et al.  A high resolution numerical study of Gulf of Mexico fronts and eddies , 1997 .

[21]  R. Winther,et al.  Numerical methods for incompressible viscous flow , 2002 .

[22]  Mustafa S. Altinakar,et al.  Efficient shallow water simulations on GPUs: Implementation, visualization, verification, and validation , 2012 .

[23]  Y. Notay.  An aggregation-based algebraic multigrid method , 2010 .

[24]  T. H. Bell,et al.  Lee waves in stratified flows with simple harmonic time dependence , 1975, Journal of Fluid Mechanics.

[25]  F. Menter.  Zonal Two Equation Kappa-Omega Turbulence Models for Aerodynamic Flows , 1993 .

[26]  Akio Arakawa,et al.  Computational Design of the Basic Dynamical Processes of the UCLA General Circulation Model , 1977 .

[27]  A. Blumberg,et al.  Diagnostic and prognostic numerical circulation studies of the South Atlantic Bight , 1983 .

[28]  James C. McWilliams,et al.  The Regional Ocean Modeling System : A Split-Explicit , Free-Surface , Topography-Following-Coordinate Ocean Model , 2022 .

[29]  Germán Larrazábal,et al.  Large sparse linear systems arising from mimetic discretization , 2007, Comput. Math. Appl..

[30]  J. Crank,et al.  A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type , 1947 .

[31]  D. Wright,et al.  Internal Tide Generation over Topography: Experiments with a Free-Surface z-Level Ocean Model , 2001 .

[32]  K. Hoffmann,et al.  Computational Fluid Dynamics for Engineers , 1989 .

[33]  Mikhail Shashkov,et al.  Fourth- and sixth-order conservative finite difference approximations of the divergence and gradient , 2001 .

[34]  J. Neumann,et al.  Numerical Integration of the Barotropic Vorticity Equation , 1950 .

[35]  R. Hallberg.  Some aspects of the circulation in ocean basins with isopycnals intersecting the sloping boundaries , 1995 .

[36]  R. Stull.  An Introduction to Boundary Layer Meteorology , 1988 .

[37]  José E. Castillo,et al.  Unified Curvilinear Ocean Atmosphere Model (UCOAM): A vertical velocity case study , 2013, Math. Comput. Model..

[38]  Patrick J. Hogan,et al.  The Gulf Stream Pathway and the Impacts of the Eddy-Driven Abyssal Circulation and the Deep Western Boundary Current , 2008 .

[39]  U. Schumann,et al.  Numerical simulation of turbulent convection over wavy terrain , 1992, Journal of Fluid Mechanics.

[40]  G. Petrova,et al.  A SECOND-ORDER WELL-BALANCED POSITIVITY PRESERVING CENTRAL-UPWIND SCHEME FOR THE SAINT-VENANT SYSTEM ∗ , 2007 .

[41]  B. Geurts.  Elements of direct and large-eddy simulation , 2003 .

[42]  P. Moin,et al.  Application of a Fractional-Step Method to Incompressible Navier-Stokes Equations , 1984 .

[43]  Alexander F. Shchepetkin,et al.  Model evaluation experiments in the North Atlantic Basin : simulations in nonlinear terrain-following coordinates , 2000 .

[44]  K. Kuwahara,et al.  Computation of high Reynolds number flow around a circular cylinder with surface roughness , 1984 .

[45]  Alistair G.L. Borthwick,et al.  River and reservoir flow modelling using the transformed shallow water equations , 1992 .

[46]  Robert L. Haney,et al.  On the Pressure Gradient Force over Steep Topography in Sigma Coordinate Ocean Models , 1991 .

[47]  Alistair Adcroft,et al.  A New Treatment of the Coriolis Terms in C-Grid Models at Both High and Low Resolutions , 1999 .

[48]  A. Thurnherr,et al.  Hydrography and high-temperature heat flux of the Rainbow hydrothermal site (36°14′N, Mid-Atlantic Ridge) , 2001 .

[49]  Jean-Christophe Golaz,et al.  One-Way Nested Large-Eddy Simulation over the Askervein Hill , 2009 .

[50]  C. Wunsch,et al.  How well does a 1/4° global circulation model simulate large-scale oceanic observations? , 1996 .

[51]  A. Chorin.  A Numerical Method for Solving Incompressible Viscous Flow Problems , 1997 .

[52]  R. Pielke,et al.  The forward-in-time upstream advection scheme:extension to higher orders , 1987 .

[53]  José E. Castillo,et al.  A high performance GPU implementation of Surface Energy Balance System (SEBS) based on CUDA-C , 2013, Environ. Model. Softw..

[54]  G. Powers,et al.  A Description of the Advanced Research WRF Version 3 , 2008 .

[55]  R. Salmon,et al.  Geophysical Fluid Dynamics , 2019, Classical Mechanics in Geophysical Fluid Dynamics.

[56]  J. Marshall,et al.  Atmosphere, Ocean and Climate Dynamics: An Introductory Text , 1961 .

[57]  B. Cushman-Roisin.  Introduction to Geophysical Fluid Dynamics , 1994 .

[58]  P. Roache.  Fundamentals of computational fluid dynamics , 1998 .

[59]  José E. Castillo,et al.  Linear Systems Arising for Second-Order Mimetic Divergence and Gradient Discretizations , 2005, J. Math. Model. Algorithms.

[60]  D. Birchall,et al.  Computational Fluid Dynamics , 2020, Radial Flow Turbocompressors.

[61]  H. Hasumi,et al.  Developments in ocean climate modelling , 2000 .

[62]  A. Arakawa,et al.  A Potential Enstrophy and Energy Conserving Scheme for the Shallow Water Equations , 1981 .

[63]  B. Launder,et al.  Application of the energy-dissipation model of turbulence to the calculation of flow near a spinning disc , 1974 .

[64]  J. Castillo,et al.  MIMETIC SCHEMES ON NON-UNIFORM STRUCTURED MESHES , 2009 .

[65]  J. Smagorinsky,et al.  GENERAL CIRCULATION EXPERIMENTS WITH THE PRIMITIVE EQUATIONS , 1963 .

[66]  Chi-Wang Shu,et al.  Strong Stability-Preserving High-Order Time Discretization Methods , 2001, SIAM Rev..

[67]  Artem Napov,et al.  An Algebraic Multigrid Method with Guaranteed Convergence Rate , 2012, SIAM J. Sci. Comput..

[68]  Jean-Marc Molines,et al.  Circulation characteristics in three eddy-permitting models of the North Atlantic , 2001 .

[69]  Yvan Notay,et al.  Aggregation-Based Algebraic Multigrid for Convection-Diffusion Equations , 2012, SIAM J. Sci. Comput..

[70]  Hernan G. Arango,et al.  Developments in terrain-following ocean models: intercomparisons of numerical aspects , 2002 .

[71]  A. Arakawa.  Computational design for long-term numerical integration of the equations of fluid motion: two-dimen , 1997 .

[72]  D. Sandwell,et al.  Global correlation of mesoscale ocean variability with seafloor roughness from satellite altimetry , 2000 .

[73]  L. Perelman,et al.  A finite-volume, incompressible Navier Stokes model for studies of the ocean on parallel computers , 1997 .

[74]  Louis J. Wicker,et al.  A Time-Splitting Scheme for the Elastic Equations Incorporating Second-Order Runge–Kutta Time Differencing , 1998 .

[75]  Andrew F. Loughe,et al.  A Reduced-Gravity Isopycnal Ocean Model: Hindcasts of El Niño , 1995 .

[76]  Robert Hallberg,et al.  Stable Split Time Stepping Schemes for Large-Scale Ocean Modeling , 1997 .

[77]  Justin Brent Runyan,et al.  A novel higher order Finite Difference Time Domain method based on the Castillo-Grone mimetic curl operator with applications concerning the time-dependent Maxwell Equations , 2011 .

[78]  John D. McCalpin.  A comparison of second‐order and fourth‐order pressure gradient algorithms in a σ‐co‐ordinate ocean model , 1994 .

[79]  J. Castillo.  Mathematical Aspects of Numerical Grid Generation , 1991, Frontiers in Applied Mathematics.

[80]  Elias Balaras,et al.  The inner–outer layer interface in large-eddy simulations with wall-layer models , 2003 .

[81]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[82]  Frank J. Millero,et al.  International one-atmosphere equation of state of seawater , 1981 .

[83]  George G. O'Brien,et al.  A Study of the Numerical Solution of Partial Differential Equations , 1950 .

[84]  A. Arakawa,et al.  Numerical methods used in atmospheric models , 1976 .

[85]  G. F. Lane-Serff,et al.  Flow and mixing in the rift valley of the Mid-Atlantic Ridge , 2002 .

[86]  N. P. Fofonoff,et al.  Algorithms for Computation of Fundamental Properties of Seawater. Endorsed by Unesco/SCOR/ICES/IAPSO Joint Panel on Oceanographic Tables and Standards and SCOR Working Group 51. Unesco Technical Papers in Marine Science, No. 44. , 1983 .

[87]  J. Klemp,et al.  The Simulation of Three-Dimensional Convective Storm Dynamics , 1978 .

[88]  Rainer Feistel,et al.  Accurate and Computationally Efficient Algorithms for Potential Temperature and Density of Seawater , 2003 .

[89]  Frank J. Millero,et al.  A new high pressure equation of state for seawater , 1982 .

[90]  Malcolm J. Bowman,et al.  Mediterranean Overflow Water (MOW) simulation using a coupled multiple-grid Mediterranean Sea/North Atlantic Ocean model , 2008 .

[91]  George L. Mellor,et al.  A generalized coordinate ocean model and a comparison of the bottom boundary layer dynamics in terrain-following and in z-level grids , 2004 .

[92]  Giacomo Bertoldi,et al.  Accounting for atmospheric boundary layer variability on flux estimation from RS observations , 2008 .

[93]  Guriĭ Ivanovich Marchuk,et al.  Numerical Methods in Weather Prediction , 1974 .

[94]  M. G. Marietta,et al.  AN OCEAN MODELLING SYSTEM WITH TURBULENT BOUNDARY LAYERS AND TOPOGRAPHY: NUMERICAL DESCRIPTION , 1987 .

[95]  W. Large,et al.  Large-eddy simulation of the equatorial ocean boundary layer: Diurnal cycling, eddy viscosity, and horizontal rotation , 1996 .

[96]  Guillermo Miranda,et al.  Mimetic Discretization Methods , 2013 .

[97]  Richard J. Greatbatch,et al.  CANDIE: A New Version of the DieCAST Ocean Circulation Model , 1998 .