Implementation of the LANS-alpha turbulence model in a primitive equation ocean model

This paper presents the first numerical implementation and tests of the Lagrangian-averaged Navier-Stokes-alpha (LANS-alpha) turbulence model in a primitive equation ocean model. The ocean model in which we work is the Los Alamos Parallel Ocean Program (POP); we refer to POP and our implementation of LANS-alpha as POP-alpha. Two versions of POP-alpha are presented: the full POP-alpha algorithm is derived from the LANS-alpha primitive equations, but requires a nested iteration that makes it too slow for practical simulations; a reduced POP-alpha algorithm is proposed, which lacks the nested iteration and is two to three times faster than the full algorithm. The reduced algorithm does not follow from a formal derivation of the LANS-alpha model equations. Despite this, simulations of the reduced algorithm are nearly identical to the full algorithm, as judged by globally averaged temperature and kinetic energy, and snapshots of temperature and velocity fields. Both POP-alpha algorithms can run stably with longer timesteps than standard POP. Comparison of implementations of full and reduced POP-alpha algorithms are made within an idealized test problem that captures some aspects of the Antarctic Circumpolar Current, a problem in which baroclinic instability is prominent. Both POP-alpha algorithms produce statistics that resemble higher-resolution simulations of standard POP. A linear stability analysis shows that both the full and reduced POP-alpha algorithms benefit from the way the LANS-alpha equations take into account the effects of the small scales on the large. Both algorithms (1) are stable; (2) make the Rossby Radius effectively larger; and (3) slow down Rossby and gravity waves.Comment: Submitted to J. Computational Physics March 21, 200

SpECTRE: A Task-based Discontinuous Galerkin Code for Relativistic Astrophysics

We introduce a new relativistic astrophysics code, SpECTRE, that combines a discontinuous Galerkin method with a task-based parallelism model. SpECTRE's goal is to achieve more accurate solutions for challenging relativistic astrophysics problems such as core-collapse supernovae and binary neutron star mergers. The robustness of the discontinuous Galerkin method allows for the use of high-resolution shock capturing methods in regions where (relativistic) shocks are found, while exploiting high-order accuracy in smooth regions. A task-based parallelism model allows efficient use of the largest supercomputers for problems with a heterogeneous workload over disparate spatial and temporal scales. We argue that the locality and algorithmic structure of discontinuous Galerkin methods will exhibit good scalability within a task-based parallelism framework. We demonstrate the code on a wide variety of challenging benchmark problems in (non)-relativistic (magneto)-hydrodynamics. We demonstrate the code's scalability including its strong scaling on the NCSA Blue Waters supercomputer up to the machine's full capacity of 22,380 nodes using 671,400 threads.Comment: 41 pages, 13 figures, and 7 tables. Ancillary data contains simulation input file

Distributed simulation of city inundation by coupled surface and subsurface porous flow for urban flood decision support system

We present a decision support system for flood early warning and disaster management. It includes the models for data-driven meteorological predictions, for simulation of atmospheric pressure, wind, long sea waves and seiches; a module for optimization of flood barrier gates operation; models for stability assessment of levees and embankments, for simulation of city inundation dynamics and citizens evacuation scenarios. The novelty of this paper is a coupled distributed simulation of surface and subsurface flows that can predict inundation of low-lying inland zones far from the submerged waterfront areas, as observed in St. Petersburg city during the floods. All the models are wrapped as software services in the CLAVIRE platform for urgent computing, which provides workflow management and resource orchestration.Comment: Pre-print submitted to the 2013 International Conference on Computational Scienc

Multiscale, Multiphysics Modelling of Coastal Ocean Processes: Paradigms and Approaches

This Special Issue includes papers on physical phenomena, such as wind-driven flows, coastal flooding, and turbidity currents, and modeling techniques, such as model comparison, model coupling, parallel computation, and domain decomposition. These papers illustrate the need for modeling coastal ocean flows with multiple physical processes at different scales. Additionally, these papers reflect the current status of such modeling of coastal ocean flows, and they present a roadmap with numerical methods, data collection, and artificial intelligence as future endeavors

some remarks about a community open source lagrangian pollutant transport and dispersion model

Nowadays fishes and mussels farming is very important, from an economical point of view, for the local social background of the Bay of Naples. Hence, the accurate forecast of marine pollution becomes crucial to have reliable evaluation of its adverse effects on coastal inhabitants' health. The use of connected smart devices for monitoring the sea water pollution is getting harder because of the saline environment, the network availability and the maintain and calibration costs2. To this purpose, we designed and implemented WaComM (Water Community Model), a community open source model for sea pollutants transport and dispersion. WaComM is a model component of a scientific workflow which allows to perform, on a dedicated computational infrastructure, numerical simulations providing spatial and temporal high-resolution predictions of weather and marine conditions of the Bay of Naples leveraging on the cloud based31FACE-IT workflow engine27. In this paper we present some remarks about the development of WaComM, using hierarchical parallelism which implies distributed memory, shared memory and GPGPUs. Some numerical details are also discussed. Peer-review under responsibility of the Conference Program Chairs

Modelling the Eddy Pattern in the harbour of Zeebrugge

Hydrodynamic