A FEniCS-based programming framework for modeling turbulent flow by the Reynolds-averaged Navier-Stokes equations

Finding an appropriate turbulence model for a given flow case usually calls for extensive experimentation with both models and numerical solution methods. This work presents the design and implementation of a flexible, programmable software framework for assisting with numerical experiments in computational turbulence. The framework targets Reynolds-averaged Navier-Stokes models, discretized by finite element methods. The novel implementation makes use of Python and the FEniCS package, the combination of which leads to compact and reusable code, where model- and solver-specific code resemble closely the mathematical formulation of equations and algorithms. The presented ideas and programming techniques are also applicable to other fields that involve systems of nonlinear partial differential equations. We demonstrate the framework in two applications and investigate the impact of various linearizations on the convergence properties of nonlinear solvers for a Reynolds-averaged Navier-Stokes model.     

Sedimentary and structural control on pockmark development—evidence from the Nyegga pockmark field, NW European margin

The Nyegga region, located at water depths of about 600–800 m on the NW European continental margin, contains more than 200 pockmarks. Recently collected TOPAS seismic profiles and EM1002 bathymetric records now provide high-resolution information on their seabed and shallow sub-seabed geological setting. The identified pockmarks are up to 15 m deep, between 30 m and 600 m across and reach a maximum area of ca. 315,000 m². The pockmarks are sediment-empty features. They do not have any preferred direction of orientation and show large variations in their shape. The pockmarks are restricted to a <16.2 cal. ka old sediment unit. This unit comprises sandy mud and is characterised by sedimentation rates of ca. 1 mm/year. The pockmarks are localised over a thick late Plio-Pleistocene prograding sediment package and a polygonal faulted Miocene-Oligocene ooze-rich unit. The late Plio-Plistocene deposits host bottom simulating reflectors, indicative of gas hydrate-bearing sediments. Inspection of the newly collected high-resolution dataset, combined with previously analysed sediment cores and 2D multichannel seismic profiles, reveals that the Nyegga pockmark field does not show any strong relationship between seabed features, sub-seabed structures and the sedimentary setting. This suggests a more complex evolution history of the Nyegga pockmark field then previously thought.