Global Existence and Regularity for the 3D Stochastic Primitive Equations of the Ocean and Atmosphere with Multiplicative White Noise

The Primitive Equations are a basic model in the study of large scale Oceanic and Atmospheric dynamics. These systems form the analytical core of the most advanced General Circulation Models. For this reason and due to their challenging nonlinear and anisotropic structure the Primitive Equations have recently received considerable attention from the mathematical community. In view of the complex multi-scale nature of the earth's climate system, many uncertainties appear that should be accounted for in the basic dynamical models of atmospheric and oceanic processes. In the climate community stochastic methods have come into extensive use in this connection. For this reason there has appeared a need to further develop the foundations of nonlinear stochastic partial differential equations in connection with the Primitive Equations and more generally. In this work we study a stochastic version of the Primitive Equations. We establish the global existence of strong, pathwise solutions for these equations in dimension 3 for the case of a nonlinear multiplicative noise. The proof makes use of anisotropic estimates, $L^{p}_{t}L^{q}_{x}$ estimates on the pressure and stopping time arguments.Comment: To appear in Nonlinearit

Spectral/hp element methods: recent developments, applications, and perspectives

The spectral/hp element method combines the geometric flexibility of the classical h-type finite element technique with the desirable numerical properties of spectral methods, employing high-degree piecewise polynomial basis functions on coarse finite element-type meshes. The spatial approximation is based upon orthogonal polynomials, such as Legendre or Chebychev polynomials, modified to accommodate C0-continuous expansions. Computationally and theoretically, by increasing the polynomial order p, high-precision solutions and fast convergence can be obtained and, in particular, under certain regularity assumptions an exponential reduction in approximation error between numerical and exact solutions can be achieved. This method has now been applied in many simulation studies of both fundamental and practical engineering flows. This paper briefly describes the formulation of the spectral/hp element method and provides an overview of its application to computational fluid dynamics. In particular, it focuses on the use the spectral/hp element method in transitional flows and ocean engineering. Finally, some of the major challenges to be overcome in order to use the spectral/hp element method in more complex science and engineering applications are discussed

Time Discrete Approximation of Weak Solutions for Stochastic Equations of Geophysical Fluid Dynamics and Applications

As a first step towards the numerical analysis of the stochastic primitive equations of the atmosphere and oceans, we study their time discretization by an implicit Euler scheme. From deterministic viewpoint the 3D Primitive Equations are studied with physically realistic boundary conditions. From probabilistic viewpoint we consider a wide class of nonlinear, state dependent, white noise forcings. The proof of convergence of the Euler scheme covers the equations for the oceans, atmosphere, coupled oceanic-atmospheric system and other geophysical equations. We obtain the existence of solutions weak in PDE and probabilistic sense, a result which is new by itself to the best of our knowledge

Multiscale Interactions in Geophysical Fluids

The dynamics of the atmosphere and ocean involves a broad range of spatial and temporal scales, many of which emerge through complex nonlinear mechanisms from forcings at very different scales. This poses major challenges for the numerical prediction of the weather, ocean state and climate: many processes have scales that are too small to be resolved yet they play an essential role in determining large-scale features. This workshop examined how modern mathematical methods – ranging from multiscale asymptotics to adaptive numerical methods and stochastic modelling – can be applied to represent the large-scale impact of these small-scale processes and improve both deterministic and probabilistic predictions

Stochastic representation of the Reynolds transport theorem: revisiting large-scale modeling

We explore the potential of a formulation of the Navier-Stokes equations incorporating a random description of the small-scale velocity component. This model, established from a version of the Reynolds transport theorem adapted to a stochastic representation of the flow, gives rise to a large-scale description of the flow dynamics in which emerges an anisotropic subgrid tensor, reminiscent to the Reynolds stress tensor, together with a drift correction due to an inhomogeneous turbulence. The corresponding subgrid model, which depends on the small scales velocity variance, generalizes the Boussinesq eddy viscosity assumption. However, it is not anymore obtained from an analogy with molecular dissipation but ensues rigorously from the random modeling of the flow. This principle allows us to propose several subgrid models defined directly on the resolved flow component. We assess and compare numerically those models on a standard Green-Taylor vortex flow at Reynolds 1600. The numerical simulations, carried out with an accurate divergence-free scheme, outperform classical large-eddies formulations and provides a simple demonstration of the pertinence of the proposed large-scale modeling

Stochastic Transport in Upper Ocean Dynamics

This open access proceedings volume brings selected, peer-reviewed contributions presented at the Stochastic Transport in Upper Ocean Dynamics (STUOD) 2021 Workshop, held virtually and in person at the Imperial College London, UK, September 20–23, 2021. The STUOD project is supported by an ERC Synergy Grant, and led by Imperial College London, the National Institute for Research in Computer Science and Automatic Control (INRIA) and the French Research Institute for Exploitation of the Sea (IFREMER). The project aims to deliver new capabilities for assessing variability and uncertainty in upper ocean dynamics. It will provide decision makers a means of quantifying the effects of local patterns of sea level rise, heat uptake, carbon storage and change of oxygen content and pH in the ocean. Its multimodal monitoring will enhance the scientific understanding of marine debris transport, tracking of oil spills and accumulation of plastic in the sea. All topics of these proceedings are essential to the scientific foundations of oceanography which has a vital role in climate science. Studies convened in this volume focus on a range of fundamental areas, including: Observations at a high resolution of upper ocean properties such as temperature, salinity, topography, wind, waves and velocity; Large scale numerical simulations; Data-based stochastic equations for upper ocean dynamics that quantify simulation error; Stochastic data assimilation to reduce uncertainty. These fundamental subjects in modern science and technology are urgently required in order to meet the challenges of climate change faced today by human society. This proceedings volume represents a lasting legacy of crucial scientific expertise to help meet this ongoing challenge, for the benefit of academics and professionals in pure and applied mathematics, computational science, data analysis, data assimilation and oceanography

Stratified shallow flow modelling

Environmental hydraulics covers a very wide range of applications including free surface flows in rivers. estuaries and lakes. To find engineering solutions to environmental hydraulics problems. 3D numerical modelling is nowadays widely used. However. the computation of sharp spatial gradients (such as found in stratified estuaries and lakes. around plumes near outfalls along rivers and coasts or in exchange areas of high shear). and the modelling of these processes along steep bathymetric slopes (such as found at the edge of dredged channels or of the continental shelf) remains a challenge. In addition. crude assumptions (such as the hydrostatic assumption) are often made to the primary differential equations in order to simplify the problem and enable long term prediction of environmental hydraulic changes. In this thesis. a robust adaptive mesh displacement (AMD) method is implemented and validated against the lock exchange case in particular. The AMD method aims at vertically focusing nodes within each water column to capture sharp gradients. while reducing the number of nodes or requiring prior knowledge of the flow structure. Second. a direct computation of dynamic pressure is introduced based on the equation of vertical momentum and validated against the analytical potential flow theory solution of a source-sink pair. Dynamic pressure is necessary to model destratification recirculation devices. or flow over dredge channel. or solitary waves. for instance. This direct computation method makes the hydrostatic assumption redundant. Third. a new advection scheme is implemented. whose main advantage is simplicity averaging over Riemann problems without solving them. while excessive numerical viscosity is compensated for by using high-resolution MUSCL type reconstruction. Recommendations are made in this thesis to extend the advection scheme developed herein for tracer advection to the non-linear shallow water equations. to the diffusion terms and to turbulence closure laws within the same finite element framework

Stochastic Transport in Upper Ocean Dynamics

This open access proceedings volume brings selected, peer-reviewed contributions presented at the Stochastic Transport in Upper Ocean Dynamics (STUOD) 2021 Workshop, held virtually and in person at the Imperial College London, UK, September 20–23, 2021. The STUOD project is supported by an ERC Synergy Grant, and led by Imperial College London, the National Institute for Research in Computer Science and Automatic Control (INRIA) and the French Research Institute for Exploitation of the Sea (IFREMER). The project aims to deliver new capabilities for assessing variability and uncertainty in upper ocean dynamics. It will provide decision makers a means of quantifying the effects of local patterns of sea level rise, heat uptake, carbon storage and change of oxygen content and pH in the ocean. Its multimodal monitoring will enhance the scientific understanding of marine debris transport, tracking of oil spills and accumulation of plastic in the sea. All topics of these proceedings are essential to the scientific foundations of oceanography which has a vital role in climate science. Studies convened in this volume focus on a range of fundamental areas, including: Observations at a high resolution of upper ocean properties such as temperature, salinity, topography, wind, waves and velocity; Large scale numerical simulations; Data-based stochastic equations for upper ocean dynamics that quantify simulation error; Stochastic data assimilation to reduce uncertainty. These fundamental subjects in modern science and technology are urgently required in order to meet the challenges of climate change faced today by human society. This proceedings volume represents a lasting legacy of crucial scientific expertise to help meet this ongoing challenge, for the benefit of academics and professionals in pure and applied mathematics, computational science, data analysis, data assimilation and oceanography

Assessment of self-adapting local projection-based solvers for laminar and turbulent industrial flows

In this work, we study the performance of some local projection-based solvers in the Large Eddy Simulation (LES) of laminar and turbulent flows governed by the incompressible Navier\u2013Stokes Equations (NSE). On one side, we focus on a high-order term-by-term stabilization Finite Element (FE) method that has one level, in the sense that it is defined on a single mesh, and in which the projection-stabilized structure of standard Local Projection Stabilization (LPS) methods is replaced by an interpolation-stabilized structure. The interest of LPS methods is that they ensure a self-adapting high accuracy in laminar regions of turbulent flows, which turns to be of overall optimal high accuracy if the flow is fully laminar. On the other side, we propose a new Reduced Basis (RB) Variational Multi-Scale (VMS)-Smargorinsky turbulence model, based upon an empirical interpolation of the sub-grid eddy viscosity term. This method yields dramatical improvements of the computing time for benchmark flows. An overview about known results from the numerical analysis of the proposed methods is given, by highlighting the used mathematical tools. In the numerical study, we have considered two well known problems with applications in industry: the (3D) turbulent flow in a channel and the (2D/3D) recirculating flow in a lid-driven cavity

Das unstetige Galerkinverfahren für Strömungen mit freier Oberfläche und im Grundwasserbereich in geophysikalischen Anwendungen

Free surface flows and subsurface flows appear in a broad range of geophysical applications and in many environmental settings situations arise which even require the coupling of free surface and subsurface flows. Many of these application scenarios are characterized by large domain sizes and long simulation times. Hence, they need considerable amounts of computational work to achieve accurate solutions and the use of efficient algorithms and high performance computing resources to obtain results within a reasonable time frame is mandatory. Discontinuous Galerkin methods are a class of numerical methods for solving differential equations that share characteristics with methods from the finite volume and finite element frameworks. They feature high approximation orders, offer a large degree of flexibility, and are well-suited for parallel computing. This thesis consists of eight articles and an extended summary that describe the application of discontinuous Galerkin methods to mathematical models including free surface and subsurface flow scenarios with a strong focus on computational aspects. It covers discretization and implementation aspects, the parallelization of the method, and discrete stability analysis of the coupled model.Für viele geophysikalische Anwendungen spielen Strömungen mit freier Oberfläche und im Grundwasserbereich oder sogar die Kopplung dieser beiden eine zentrale Rolle. Oftmals charakteristisch für diese Anwendungsszenarien sind große Rechengebiete und lange Simulationszeiten. Folglich ist das Berechnen akkurater Lösungen mit beträchtlichem Rechenaufwand verbunden und der Einsatz effizienter Lösungsverfahren sowie von Techniken des Hochleistungsrechnens obligatorisch, um Ergebnisse innerhalb eines annehmbaren Zeitrahmens zu erhalten. Unstetige Galerkinverfahren stellen eine Gruppe numerischer Verfahren zum Lösen von Differentialgleichungen dar, und kombinieren Eigenschaften von Methoden der Finiten Volumen- und Finiten Elementeverfahren. Sie ermöglichen hohe Approximationsordnungen, bieten einen hohen Grad an Flexibilität und sind für paralleles Rechnen gut geeignet. Diese Dissertation besteht aus acht Artikeln und einer erweiterten Zusammenfassung, in diesen die Anwendung unstetiger Galerkinverfahren auf mathematische Modelle inklusive solcher für Strömungen mit freier Oberfläche und im Grundwasserbereich beschrieben wird. Die behandelten Themen umfassen Diskretisierungs- und Implementierungsaspekte, die Parallelisierung der Methode sowie eine diskrete Stabilitätsanalyse des gekoppelten Modells