Implementation of on-site velocity boundary conditions for D3Q19 lattice Boltzmann

On-site boundary conditions are often desired for lattice Boltzmann simulations of fluid flow in complex geometries such as porous media or microfluidic devices. The possibility to specify the exact position of the boundary, independent of other simulation parameters, simplifies the analysis of the system. For practical applications it should allow to freely specify the direction of the flux, and it should be straight forward to implement in three dimensions. Furthermore, especially for parallelized solvers it is of great advantage if the boundary condition can be applied locally, involving only information available on the current lattice site. We meet this need by describing in detail how to transfer the approach suggested by Zou and He to a D3Q19 lattice. The boundary condition acts locally, is independent of the details of the relaxation process during collision and contains no artificial slip. In particular, the case of an on-site no-slip boundary condition is naturally included. We test the boundary condition in several setups and confirm that it is capable to accurately model the velocity field up to second order and does not contain any numerical slip.Comment: 13 pages, 4 figures, revised versio

Quantitative analysis of numerical estimates for the permeability of porous media from lattice-Boltzmann simulations

During the last decade, lattice-Boltzmann (LB) simulations have been improved to become an efficient tool for determining the permeability of porous media samples. However, well known improvements of the original algorithm are often not implemented. These include for example multirelaxation time schemes or improved boundary conditions, as well as different possibilities to impose a pressure gradient. This paper shows that a significant difference of the calculated permeabilities can be found unless one uses a carefully selected setup. We present a detailed discussion of possible simulation setups and quantitative studies of the influence of simulation parameters. We illustrate our results by applying the algorithm to a Fontainebleau sandstone and by comparing our benchmark studies to other numerical permeability measurements in the literature.Comment: 14 pages, 11 figure

Improved axisymmetric lattice Boltzmann scheme

This paper proposes an improved lattice Boltzmann scheme for incompressible axisymmetric flows. The scheme has the following features. First, it is still within the framework of the standard lattice Boltzmann method using the single-particle density distribution function and consistent with the philosophy of the lattice Boltzmann method. Second, the source term of the scheme is simple and contains no velocity gradient terms. Owing to this feature, the scheme is easy to implement. In addition, the singularity problem at the axis can be appropriately handled without affecting an important advantage of the lattice Boltzmann method: the easy treatment of boundary conditions. The scheme is tested by simulating Hagen-Poiseuille flow, three-dimensional Womersley flow, Wheeler benchmark problem in crystal growth, and lid-driven rotational flow in cylindrical cavities. It is found that the numerical results agree well with the analytical solutions and/or the results reported in previous studies.Comment: 31 pages, 7 figures

Seasonal diet pattern of non-native tubenose goby (Proterorhinus semilunaris) in a lowland reservoir (Mušov, Czech Republic)

The tubenose goby (Proterorhinus semilunaris) is a gobiid species currently extending its area of distribution in Central Europe. The objective of the study was to evaluate the annual pattern of its feeding habits in the newly colonised habitats of the Mušov reservoir on the Dyje River (the Danube basin, Czech Republic) with respect to natural food resources. In the reservoir, tubenose goby has established a numerous population, densely colonising stony rip-rap banks. Its diet was exclusively of animal origin with significant dominance of and preference for two food items – chironomid (Chironomidae) larvae and waterlouse (Asellus aquaticus), which contributed 40.2 and 27.6%, respectively, to the total food bulk ingested. The index of preponderance for the two items was also very high, amounting to 73.8 and 26.5, respectively. In the annual pattern, a remarkable preference for chironomid larvae was recorded in the summer period whilst waterlouse were consumed predominantly in winter months. The proportion of other food items was rather marginal – only corixids, copepods, ceratopogonids and cladocerans were of certain minor importance with proportions of 5.4, 4.3, 4.1 and 3.9%, respectively. Certain signs of cannibalism were also recorded, with 0.9 and 0.2% of the diet consisting of their own progeny and eggs, respectively

Seasonal diet pattern of non-native tubenose goby (

The tubenose goby (Proterorhinus semilunaris) is a gobiid species currently extending its area of distribution in Central Europe. The objective of the study was to evaluate the annual pattern of its feeding habits in the newly colonised habitats of the Mušov reservoir on the Dyje River (the Danube basin, Czech Republic) with respect to natural food resources. In the reservoir, tubenose goby has established a numerous population, densely colonising stony rip-rap banks. Its diet was exclusively of animal origin with significant dominance of and preference for two food items – chironomid (Chironomidae) larvae and waterlouse (Asellus aquaticus), which contributed 40.2 and 27.6%, respectively, to the total food bulk ingested. The index of preponderance for the two items was also very high, amounting to 73.8 and 26.5, respectively. In the annual pattern, a remarkable preference for chironomid larvae was recorded in the summer period whilst waterlouse were consumed predominantly in winter months. The proportion of other food items was rather marginal – only corixids, copepods, ceratopogonids and cladocerans were of certain minor importance with proportions of 5.4, 4.3, 4.1 and 3.9%, respectively. Certain signs of cannibalism were also recorded, with 0.9 and 0.2% of the diet consisting of their own progeny and eggs, respectively

Modeling the Impact of Restoration Efforts on Phosphorus Loading and Transport through Everglades National Park, FL, USA

Ecosystems of Florida Everglades are highly sensitive to phosphorus loading. Future restoration efforts, which focus on restoring Everglades water flows, may pose a threat to the health of these ecosystems. To determine the fate and transport of total phosphorus and evaluate proposed Everglades restoration, a water quality model has been developed using the hydrodynamic results from the M3ENP (Mike Marsh Model of Everglades National Park) — a physically-based hydrological numerical model which uses MIKE SHE/MIKE 11 software. Using advection–dispersion with reactive transport for the model, parameters were optimized and phosphorus loading in the overland water column was modeled with good accuracy (60%). The calibrated M3ENP-AD model was then modified to include future bridge construction and canal water level changes, which have shown to increase flows into ENP. These bridge additions increased total dissolved phosphorus (TP) load downstream in Shark Slough and decreased TP load in downstream Taylor Slough. However, there was a general decrease in TP concentration and TP mass per area over the entire model domain. The M3ENP-AD model has determined the mechanisms for TP transport and quantified the impacts of ENP restoration efforts on the spatial–temporal distribution of phosphorus transport. This tool can be used to guide future Everglades restoration decisions