Optimization of the numerical algorithms of the ADREA-I mesoscale prognostic meteorological model for real-time applications

The real-time applicability of the ADREA-I prognostic mesoscale meteorological model was enhanced by applying the preconditioned BiCGSTAB method for the numerical solution of the pressure equation in combination with increasing the magnitude of the time steps up to the values allowed by the Courant number. The ILU, MILU, ILUT and ILUM preconditioning methods with different ordering strategies were used. The implementation was developed for arbitrarily complex geometries. The application of MILU(1) preconditioning and ILUT preconditioning with red-black ordering of the unknowns (RB+ILUT) has resulted in up to six times shorter overall computational time in comparison to the previously implemented line relaxation (LR) method. The feasibility of increasing the time steps has been proved by comparing the results of the 24-h ADREA-I forecasts with the observations during a real sea breeze case in Attiki, Greece: decreasing the time steps by a factor of 10 in comparison with the values allowed by the Courant number leads to decrease of the statistical error indicators by only 1-5%.

[1]  J. Ehrhardt,et al.  The RODOS System: Decision Support for Off-Site Emergency Management in Europe , 1997 .

[2]  J. Bartzis,et al.  Prediction of orographic precipitation using cartesian coordinates and a single prognostic equation for the water substance , 1991 .

[3]  G. Grell,et al.  A description of the fifth-generation Penn State/NCAR Mesoscale Model (MM5) , 1994 .

[4]  W. Raskob,et al.  The RODOS System : Decision Support For Nuclear Off-Site Emergency Management In Europe , 2000 .

[5]  G. L. Browning,et al.  The Accuracy and Efficiency of Semi-Implicit Time Stepping for Mesoscale Storm Dynamics , 2001 .

[6]  Nicolas Moussiopoulos,et al.  Direct solution of the Helmholtz equation using fourier analysis on the cyber 205 , 1988 .

[7]  Wolfgang Raskob,et al.  STATUS OF THE RODOS SYSTEM FOR OFF-SITE EMERGENCY MANAGEMENT AFTER NUCLEAR AND RADIOLOGICAL ACCIDENTS , 2006 .

[8]  E. F. F. Botta,et al.  Matrix Renumbering ILU: An Effective Algebraic Multilevel ILU Preconditioner for Sparse Matrices , 1999, SIAM J. Matrix Anal. Appl..

[9]  J. G. Bartzis,et al.  ADREA-I: A Three-Dimensional Transient Transport Code for Complex Terrain and Other Applications , 1991 .

[10]  James M. Ortega,et al.  Iterative solution of nonlinear equations in several variables , 2014, Computer science and applied mathematics.

[11]  Jun Zhang,et al.  A grid-based multilevel incomplete LU factorization preconditioning technique for general sparse matrices , 2001, Appl. Math. Comput..

[12]  C. Brand.  An incomplete-factorization preconditioning using repeated red-black ordering , 1992 .

[13]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[14]  S. Patankar.  Numerical Heat Transfer and Fluid Flow , 2018, Lecture Notes in Mechanical Engineering.

[15]  Bjarne Amstrup,et al.  Urban meteorological modelling for nuclear emergency preparedness. , 2006, Journal of environmental radioactivity.

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

[17]  C. F. Ratto,et al.  Mass-consistent models for wind fields over complex terrain: The state of the art , 1994 .

[18]  E. Bender.  Numerical heat transfer and fluid flow. Von S. V. Patankar. Hemisphere Publishing Corporation, Washington – New York – London. McGraw Hill Book Company, New York 1980. 1. Aufl., 197 S., 76 Abb., geb., DM 71,90 , 1981 .

[19]  Fred Wubs,et al.  A fully-implicit model of the global ocean circulation , 2003 .

[20]  O. Axelsson,et al.  A survey of some estimates of eigenvalues and condition numbers for certain preconditioned matrices , 1997 .

[21]  Roberto San José García,et al.  An operational real-time air quality modelling system for industrial plants , 2007, Environ. Model. Softw..

[22]  P. Wesseling.  Principles of Computational Fluid Dynamics , 2000 .

[23]  E. D'Azevedo,et al.  Towards a cost-effective ILU preconditioner with high level fill , 1992 .

[24]  John C. Statharas,et al.  Development and interaction of thermally driven flows over Attiki peninsula under northerly background wind—a case study , 1998 .

[25]  Piotr K. Smolarkiewicz,et al.  Spectral Preconditioners for Nonhydrostatic Atmospheric Models , 2003 .

[26]  J. Bartzis,et al.  The BOND project: Biogenic aerosols and air quality in Athens and Marseille greater areas , 2004 .

[27]  Yvan Notay,et al.  A multilevel block incomplete factorization preconditioning , 1999 .

[28]  Miguel Angel Medina,et al.  Enhanced TABS-MDS model for simulating large-scale free surface flows , 2006, Environ. Model. Softw..

[29]  Ivar Gustafsson,et al.  An incomplete factorization preconditioning method based on modification of element matrices , 1996 .

[30]  J. Bartzis,et al.  Simulation of the sea breeze under opposing synoptic conditions , 1993 .

[31]  Weeratunge Malalasekera,et al.  An introduction to computational fluid dynamics - the finite volume method , 2007 .