Multi-scale modeling: nested-grid and unstructured-mesh approaches

Abstract

In 1969, the Journal of Computational Physics published a seminal article by K. Bryan presenting the first ocean general circulation model. Since then, many numerical studies of the World Ocean, as well as regional or coastal flows, used models directly or indirectly inspired by the work of Bryan and his colleagues. A number of these models have evolved into highly modular and versatile computational systems, including multiple physical modules and options as well as varied biogeochemical, ecosystem and acoustics modeling capabilities. Several modeling systems are now well-documented tools, which are widely used in research institutions and various organizations around the world. The list of such modeling systems is large and too long to be summarized in this editorial. Over the last three decades, significant progress has been made in the parameterization of subgrid-scale processes, in data assimilation methodologies and in boundary condition schemes, as well as in the efficient implementation of algorithms on fast vector and subsequently parallel computers, allowing higher and higher resolution in space and time. However, many of today’s popular modeling systems can still be regarded as members of the first generation of ocean models: at their core, rather similar geophysical fluid dynamics equations are solved numerically using a conservative finite-difference method on a structured grid. Today, several aspects of structured-grid models could benefit from significant upgrades, learning from major advances in computational fluid dynamics. In particular, the use of a structured grid limits the flexibility in the spatial resolution and does not allow one to take full advantage of numerical algorithms such as finite volumes and finite elements, which can achieve their best performance when implemented on unstructured meshes. Even though many of today’s complex marine modeling and data assimilation systems have evolved significantly since Bryan’s prototype, it would be challenging to modify them step-by-step from a structured-grid approach to an unstructured-grid one. Therefore, novel marine model design research is underway, paving the way for the second generation of ocean modeling systems. It is difficult to predict today if this new generation of ocean models will achieve its chief objective: widening the range of resolved scales of motion with increased efficiencies and accuracies, possibly allowing multi-resolution, multi-scale, and multidynamics numerical simulations of marine flows, all occurring seamlessly within distributed computing environments. In fact, hybrid approaches merging the advantages of structured and unstructured-grid modeling may be the way forward. Whether or not unstructured mesh approaches will prevail is all the more difficult to predict now that structured mesh modelers have developed powerful solutions for increasing the resolution when and where needed. For instance, grid embedding is still a popular and useful method for enhancing model resolution. It can involve multiply nested domains and allows the relatively Ocean Dynamics (2008) 58:335–336 DOI 10.1007/s10236-008-0170-5