-
Attached and separated rotating flow over a finite height ridge
Authors:
Stefan Frei,
Erik Burman,
Edward R Johnson
Abstract:
This paper discusses the effect of rotation on the boundary layer in high Reynolds number flow over a ridge using a numerical method based on stabilised finite elements that captures steady solutions up to Reynolds number of order $10^6$. The results are validated against boundary layer computations in shallow flows and for deep flows against experimental observations reported in Machicoane et al.…
▽ More
This paper discusses the effect of rotation on the boundary layer in high Reynolds number flow over a ridge using a numerical method based on stabilised finite elements that captures steady solutions up to Reynolds number of order $10^6$. The results are validated against boundary layer computations in shallow flows and for deep flows against experimental observations reported in Machicoane et al. (Phys. Rev. Fluids, 2018). In all cases considered the boundary layer remains attached, even at large Reynolds numbers, provided the Rossby number of the flow is sufficiently small. At any fixed Rossby number the flow detaches at sufficiently high Reynolds number to form a steady recirculating region in the lee of the ridge. At even higher Reynolds numbers no steady flow is found. This disappearance of steady solutions closely reproduces the transition to unsteadiness seen in the laboratory.
△ Less
Submitted 23 February, 2024;
originally announced February 2024.
-
Coupling finite and boundary element methods to solve the Poisson--Boltzmann equation for electrostatics in molecular solvation
Authors:
Michal Bosy,
Matthew W. Scroggs,
Timo Betcke,
Erik Burman,
Christopher D. Cooper
Abstract:
The Poisson--Boltzmann equation is widely used to model electrostatics in molecular systems. Available software packages solve it using finite difference, finite element, and boundary element methods, where the latter is attractive due to the accurate representation of the molecular surface and partial charges, and exact enforcement of the boundary conditions at infinity. However, the boundary ele…
▽ More
The Poisson--Boltzmann equation is widely used to model electrostatics in molecular systems. Available software packages solve it using finite difference, finite element, and boundary element methods, where the latter is attractive due to the accurate representation of the molecular surface and partial charges, and exact enforcement of the boundary conditions at infinity. However, the boundary element method is limited to linear equations and piecewise constant variations of the material properties. In this work, we present a scheme that couples finite and boundary elements for the Poisson--Boltzmann equation, where the finite element method is applied in a confined {\it solute} region, and the boundary element method in the external {\it solvent} region. As a proof-of-concept exercise, we use the simplest methods available: Johnson--Nédélec coupling with mass matrix and diagonal preconditioning, implemented using the Bempp-cl and FEniCSx libraries via their Python interfaces. We showcase our implementation by computing the polar component of the solvation free energy of a set of molecules using a constant and a Gaussian-varying permittivity. We validate our implementation against the finite difference code APBS (to 0.5\%), and show scaling from protein G B1 (955 atoms) up to immunoglobulin G (20\,148 atoms). For small problems, the coupled method was efficient, outperforming a purely boundary integral approach. For Gaussian-varying permittivities, which are beyond the applicability of boundary elements alone, we were able to run medium to large sized problems on a single workstation. Development of better preconditioning techniques and the use of distributed memory parallelism for larger systems remains an area for future work. We hope this work will serve as inspiration for future developments for molecular electrostatics with implicit solvent models.
△ Less
Submitted 10 May, 2023;
originally announced May 2023.
-
Unique continuation for the Lamé system using stabilized finite element methods
Authors:
Erik Burman,
Janosch Preuss
Abstract:
We introduce an arbitrary order, stabilized finite element method for solving a unique continuation problem subject to the time-harmonic elastic wave equation with variable coefficients. Based on conditional stability estimates we prove convergence rates for the proposed method which take into account the noise level and the polynomial degree. A series of numerical experiments corroborates our the…
▽ More
We introduce an arbitrary order, stabilized finite element method for solving a unique continuation problem subject to the time-harmonic elastic wave equation with variable coefficients. Based on conditional stability estimates we prove convergence rates for the proposed method which take into account the noise level and the polynomial degree. A series of numerical experiments corroborates our theoretical results and explores additional aspects, e.g. how the quality of the reconstruction depends on the geometry of the involved domains. We find that certain convexity properties are crucial to obtain a good recovery of the wave displacement outside the data domain and that higher polynomial orders can be more efficient but also more sensitive to the ill-conditioned nature of the problem.
△ Less
Submitted 24 April, 2023; v1 submitted 12 December, 2022;
originally announced December 2022.
-
Error estimates for the Smagorinsky turbulence model: enhanced stability through scale separation and numerical stabilization
Authors:
Erik Burman,
Peter Hansbo,
Mats G. Larson
Abstract:
In the present work we show some results on the effect of the Smagorinsky model on the stability of the associated perturbation equation. We show that in the presence of a spectral gap, such that the flow can be decomposed in a large scale with moderate gradient and a small amplitude fine scale with arbitratry gradient, the Smagorinsky model admits stability estimates for perturbations, with expon…
▽ More
In the present work we show some results on the effect of the Smagorinsky model on the stability of the associated perturbation equation. We show that in the presence of a spectral gap, such that the flow can be decomposed in a large scale with moderate gradient and a small amplitude fine scale with arbitratry gradient, the Smagorinsky model admits stability estimates for perturbations, with exponential growth depending only on the large scale gradient. We then show in the context of stabilized finite element methods that the same result carries over to the approximation and that in this context, for suitably chosen finite element spaces the Smagorinsky model acts as a stabilizer yielding close to optimal error estimates in the $L^2$-norm for smooth flows in the pre-asymptotic high Reynolds number regime.
△ Less
Submitted 29 January, 2021;
originally announced February 2021.