Sebastian Götschel

Lars Stietz, Sebastian Götschel, Peter Schleper et al.

Bathymetry reconstruction is an important problem in various fields, including oceanography and environmental monitoring. This paper presents a Bayesian inference approach to reconstructing bathymetries from point measurements of the water height. We test the method for parameterized and discretized bathymetries with synthetic data to evaluate its performance and limitations. Our results indicate that the Bayesian framework provides a robust approach to bathymetry reconstruction. Finally, we use the framework to reconstruct a real-world bathymetry in a wave flume from experimental measurements and compare its performance to an adjoint optimization method. The Bayesian approach improves the normalized root mean squared error (NRMSE) of the reconstruction and provides better qualitative features, while also quantifying uncertainty.

Niklas Göschel, Sebastian Götschel, Daniel Ruprecht

Machine-learning based methods like physics-informed neural networks and physics-informed neural operators are becoming increasingly adept at solving even complex systems of partial differential equations. Boundary conditions can be enforced either weakly by penalizing deviations in the loss function or strongly by training a solution structure that inherently matches the prescribed values and derivatives. The former approach is easy to implement but the latter can provide benefits with respect to accuracy and training times. However, previous approaches to strongly enforcing Neumann or Robin boundary conditions require a domain with a fully $C^1$ boundary and, as we demonstrate, can lead to instability if those boundary conditions are posed on a segment of the boundary that is piecewise $C^1$ but only $C^0$ globally. We introduce a generalization of the approach by Sukumar \& Srivastava (doi: 10.1016/j.cma.2021.114333), and a new approach based on orthogonal projections that overcome this limitation. The performance of these new techniques is compared against weakly and semi-weakly enforced boundary conditions for the scalar Darcy flow equation and the stationary Navier-Stokes equations.

Maximilian Witte, Fabricio Rodrigues Lapolli, Philip Freese et al.

Using the nonlinear shallow water equations as benchmark, we demonstrate that a simulation with the ICON-O ocean model with a 20km resolution that is frequently corrected by a U-net-type neural network can achieve discretization errors of a simulation with 10km resolution. The network, originally developed for image-based super-resolution in post-processing, is trained to compute the difference between solutions on both meshes and is used to correct the coarse mesh every 12h. Our setup is the Galewsky test case, modeling transition of a barotropic instability into turbulent flow. We show that the ML-corrected coarse resolution run correctly maintains a balance flow and captures the transition to turbulence in line with the higher resolution simulation. After 8 day of simulation, the $L_2$-error of the corrected run is similar to a simulation run on the finer mesh. While mass is conserved in the corrected runs, we observe some spurious generation of kinetic energy.