An alternative smooth particle hydrodynamics formulation to simulate chemotaxis in porous media.

Chemotaxis, the microorganisms autonomous motility along or against the concentration gradients of a chemical species, is an important, yet often neglected factor controlling the transport of bacteria through saturated porous media. For example, chemotactic bacteria could enhance bioremediation by directing their own motion to residual contaminants trapped in low hydraulic conductive zones of contaminated aquifers. The aim of the present work is to develop an accurate numerical scheme to model chemotaxis in saturated porous media and other advective dominating flow systems. We propose to model chemotaxis by using a new class of meshless Lagrangian particle methods we recently developed for applications in fluid mechanics. The method is based on the Smooth Particle Hydrodynamics (SPH) formulation of (Ben Moussa et al., Int Ser Numer Math, 13(1):29-62, 2006), combined with a new Weighted Essentially Non-Oscillatory (WENO) reconstruction technique on moving point clouds in multiple space dimensions. The purpose of this new numerical scheme is to fully exploit the advantages of SPH among traditional mesh-based and mesh-free schemes and to overcome drawbacks related to the use of standard SPH for modeling chemotaxis in porous media. First, we test the new scheme against analytical reference solutions. Then, under the assumption of complete mixing at the Darcy scale, we perform two-dimensional conservative solute transport simulations under steady-state flow conditions, to show the capability of the proposed new scheme to model chemotaxis.

The role of 3D-hydraulics in habitat modelling of hydropeaking events.

One way to study ecological implications induced by hydropeaking represents the coupling of hydrodynamic models with habitat suitability models, in which hydrodynamic parameters are typically used to describe the physical habitat of indicator species. This article discusses the differences in habitat suitability assessment between 2D and 3D CFD modelling as input for the habitat simulation tool CASiMiR. In the first part of the article, the accuracy of the hydraulic model is evaluated by comparing the model results with laboratory (model of a laboratory channel with erodible bed) and field measurements (Valsura River, Bolzano, Italy). In the second part, the habitat suitability for the Valsura River case study (affected by hydropeaking), is analyzed comparing different approaches for the reconstruction of the velocity field (depth-averaged velocities from 2D modelling, bottom velocity field reconstruction with log-law approach from 2D modelling and bottom velocity field from 3D modelling). The results show that the habitat suitability index (HSI) using 2D or 3D hydrodynamic models can be significantly different. These differences can be ascribed to a higher capability to depict the features of the flow field with highly variable and heterogeneous boundary conditions and to the possibility to simulate the near bed hydrodynamic parameters, which are relevant for certain target species. In particular, the HSI-values using 3D hydraulics lead to larger areas of highly suitable habitats compared to 2D simulations. Moreover, considering the entire flow range of hydropeaking events, the habitat simulations with bottom flow velocities from 3D modelling provide suitable habitats over the entire flow range representing the availability of stable suitable habitats, while the habitat availability of 2D modelled flow velocity is continuously decreasing with increasing flow rates.

Efficient conservative ADER schemes based on WENO reconstruction and space-time predictor in primitive variables.

We present a new version of conservative ADER-WENO finite volume schemes, in which both the high order spatial reconstruction as well as the time evolution of the reconstruction polynomials in the local space-time predictor stage are performed in primitive variables, rather than in conserved ones. To obtain a conservative method, the underlying finite volume scheme is still written in terms of the cell averages of the conserved quantities. Therefore, our new approach performs the spatial WENO reconstruction twice: the first WENO reconstruction is carried out on the known cell averages of the conservative variables. The WENO polynomials are then used at the cell centers to compute point values of the conserved variables, which are subsequently converted into point values of the primitive variables. This is the only place where the conversion from conservative to primitive variables is needed in the new scheme. Then, a second WENO reconstruction is performed on the point values of the primitive variables to obtain piecewise high order reconstruction polynomials of the primitive variables. The reconstruction polynomials are subsequently evolved in time with a novel space-time finite element predictor that is directly applied to the governing PDE written in primitive form. The resulting space-time polynomials of the primitive variables can then be directly used as input for the numerical fluxes at the cell boundaries in the underlying conservative finite volume scheme. Hence, the number of necessary conversions from the conserved to the primitive variables is reduced to just one single conversion at each cell center. We have verified the validity of the new approach over a wide range of hyperbolic systems, including the classical Euler equations of gas dynamics, the special relativistic hydrodynamics (RHD) and ideal magnetohydrodynamics (RMHD) equations, as well as the Baer-Nunziato model for compressible two-phase flows. In all cases we have noticed that the new ADER schemes provide less oscillatory solutions when compared to ADER finite volume schemes based on the reconstruction in conserved variables, especially for the RMHD and the Baer-Nunziato equations. For the RHD and RMHD equations, the overall accuracy is improved and the CPU time is reduced by about 25 %. Because of its increased accuracy and due to the reduced computational cost, we recommend to use this version of ADER as the standard one in the relativistic framework. At the end of the paper, the new approach has also been extended to ADER-DG schemes on space-time adaptive grids (AMR).

An efficient semi-implicit method for three-dimensional non-hydrostatic flows in compliant arterial vessels.

Blood flow in arterial systems can be described by the three-dimensional Navier-Stokes equations within a time-dependent spatial domain that accounts for the elasticity of the arterial walls. In this article, blood is treated as an incompressible Newtonian fluid that flows through compliant vessels of general cross section. A three-dimensional semi-implicit finite difference and finite volume model is derived so that numerical stability is obtained at a low computational cost on a staggered grid. The key idea of the method consists in a splitting of the pressure into a hydrostatic and a non-hydrostatic part, where first a small quasi-one-dimensional nonlinear system is solved for the hydrostatic pressure and only in a second step the fully three-dimensional non-hydrostatic pressure is computed from a three-dimensional nonlinear system as a correction to the hydrostatic one. The resulting algorithm is robust, efficient, locally and globally mass conservative, and applies to hydrostatic and non-hydrostatic flows in one, two and three space dimensions. These features are illustrated on nontrivial test cases for flows in tubes with circular or elliptical cross section where the exact analytical solution is known. Test cases of steady and pulsatile flows in uniformly curved rigid and elastic tubes are presented. Wherever possible, axial velocity development and secondary flows are shown and compared with previously published results.

Semi-implicit numerical modeling of axially symmetric flows in compliant arterial systems.

Blood flow in arterial systems is described by the three-dimensional Navier-Stokes equations within a time-dependent spatial domain that accounts for the viscoelasticity of the arterial walls. These equations are simplified by assuming cylindrical geometry, axially symmetric flow, and hydrostatic equilibrium in the radial direction. In this paper, an efficient semi-implicit method is formulated in such a fashion that numerical stability is obtained at a minimal computational cost. The resulting computer model is relatively simple, robust, accurate, and extremely efficient. These features are illustrated on nontrivial test cases where the exact analytical solution is known and by an example of a realistic flow through a complex arterial system.