A conservative and consistent implicit Cartesian cut-cell method for moving geometries with reduced spurious pressure oscillations

A conservative and consistent three-dimensional Cartesian cut-cell method is presented for reducing the spurious pressure oscillations often observed in moving body simulations in sharp-interface Cartesian grid methods. By analysing the potential sources of the oscillation in the cut-cell framework, an improved moving body algorithm is proposed for the cut-cell method for the temporal discontinuity of the solid volume change. Strict conservation of mass and momentum for both fluid and cut cells is enforced through pressure-velocity coupling to reduce local mass conservation errors. A consistent mass and momentum flux computation is employed in the finite volume method. In contrary to the commonly cut-cell methods, an implicit time integration scheme is employed in the present method, which prevents numerical instability without any additional small cut-cell treatment. The effectiveness of the present cut-cell method for reducing spurious pressure oscillations is demonstrated by simulating various two- and three-dimensional benchmark cases (in-line and transversely oscillating cylinder, oscillating and free-falling sphere), with good agreement with previous experimental measurements and other numerical methods available in the literature

Response of Flow and Saltating Particle Characteristics to Bed Roughness and Particle Spatial Density

With the goal to explore the effects of natural bed roughness on bedload transport, numerical simulations of flow and particle saltation are carried out with varying bed roughness and particle spatial density. A combination of Eulerian and Lagrangian point-particle methods is applied to solve the equations of motion of the fluid and the particles within the large-eddy simulation framework. Flows over smooth and rough beds with four particle densities are considered. As the bed roughness increases, there is a leftward shift of the double-averaged streamwise velocity profiles against dimensionless vertical distance, which is scaled with the bed roughness height, an upward shift of the peak values of double-averaged Reynolds stresses, and a fragmentation and disappearance of coherent structures in the form of high-speed and low-speed near-bed streaks. These observations are consistent with those of previous studies. As the bed roughness increases, the mean resting time of saltating particles increases, however the particles' saltation length, velocity, and angular velocity decrease, while their saltation height remains almost unchanged. Saltation height, saltation length, particle angular velocity, and resting time exhibit linear, gamma, normal, and exponential distributions, respectively. Further, as the bed roughness increases, the kurtosis and skewness of some particle parameters vary, and the particle velocity shifts from a symmetrical normal to an asymmetrical gamma distribution

A Cartesian cut-cell based multiphase flow model for large-eddy simulation of three-dimensional wave-structure interaction

A multiphase flow numerical approach for performing large-eddy simulations of three-dimensional (3D) wave-structure interaction is presented in this study. The approach combines a volume-of-fluid method to capture the air-water interface and a Cartesian cut-cell method to deal with complex geometries. The filtered Navier–Stokes equations are discretised by the finite volume method with the PISO algorithm for velocity-pressure coupling and the dynamic Smagorinsky subgrid-scale model is used to compute the unresolved (subgrid) scales of turbulence. The versatility and robustness of the presented numerical approach are illustrated by applying it to solve various three-dimensional wave-structure interaction problems featuring complex geometries, such as a 3D travelling wave in a closed channel, a 3D solitary wave interacting with a vertical circular cylinder, a 3D solitary wave interacting with a horizontal thin plate, and a 3D focusing wave impacting on an FPSO-like structure. For all cases, convincing agreement between the numerical predictions and the corresponding experimental data and/or analytical or numerical solutions is obtained. In addition, for all cases, water surface profiles and turbulent vortical structures are presented and discussed

A Cartesian cut-cell based multiphase flow model for large-eddy simulation of three-dimensional wave-structure interaction

A multiphase flow numerical approach for performing large-eddy simulations of three-dimensional (3D) wave-structure interaction is presented in this study. The approach combines a volume-of-fluid method to capture the air-water interface and a Cartesian cut-cell method to deal with complex geometries. The filtered Navier–Stokes equations are discretised by the finite volume method with the PISO algorithm for velocity-pressure coupling and the dynamic Smagorinsky subgrid-scale model is used to compute the unresolved (subgrid) scales of turbulence. The versatility and robustness of the presented numerical approach are illustrated by applying it to solve various three-dimensional wave-structure interaction problems featuring complex geometries, such as a 3D travelling wave in a closed channel, a 3D solitary wave interacting with a vertical circular cylinder, a 3D solitary wave interacting with a horizontal thin plate, and a 3D focusing wave impacting on an FPSO-like structure. For all cases, convincing agreement between the numerical predictions and the corresponding experimental data and/or analytical or numerical solutions is obtained. In addition, for all cases, water surface profiles and turbulent vortical structures are presented and discussed

Impact of Environmental Turbulence on the Performance and Loadings of a Tidal Stream Turbine

A large-eddy simulation (LES) of a laboratory-scale horizontal axis tidal stream turbine operating over an irregular bathymetry in the form of dunes is performed. The Reynolds number based on the approach velocity and the chord length of the turbine blades is approximately 60,000. The simulated turbine is a 1:30 scale model of a full-scale prototype and both turbines operate at very similar tip-speed ratio of λ ≈ 3. The simulations provide quantitative evidence of the effect of seabed-induced turbulence on the instantaneous performance and structural loadings of the turbine revealing how large-scale, energetic turbulence structures affect turbine performance and bending moments of the rotor blades. The data analysis shows that wake recovery is notably enhanced in comparison to the same turbine operating above a flat-bed and this is due to the higher turbulence levels generated by the dune. The results demonstrate the need for studying in detail the flow and turbulence characteristics at potential tidal turbine deployment sites and to incorporate observed large-scale velocity and pressure fluctuations into the structural design of the turbines

Turbulent Flow through Idealized Emergent Vegetation

This paper presents results of several large-eddy simulations (LES) of turbulent flow in an open channel through staggered arrays of rigid, emergent cylinders, which can be regarded as idealized vegetation. In this study, two cylinder Reynolds numbers, RD=1,340 and RD=500, and three vegetation densities are considered. The LES of the lowest density and at RD=1,340 corresponds to a recently completed laboratory experiment, the data of which is used to validate the simulations. Fairly good agreement between calculated and measured first- and second-order statistics along measurement profiles is found, confirming the accuracy of the simulations. The high resolution of the simulations enables an explicit calculation of drag forces, decomposed into pressure and friction drag, that are exerted on the cylinders. The effect of the cylinder Reynolds number and the cylinder density on the drag and hence on the flow resistance is quantified and in agreement with previous experimental studies. Turbulence structures are visualized through instantaneous pressure fluctuations, isosurfaces of the Q-criterion and contours of vertical vorticity in horizontal planes. Analysis of velocity time signals and distributions of drag and lift forces over time reveals that flow and turbulence are more influenced by the vegetation density than by the cylinder Reynolds number

Scalability of an Eulerian-Lagrangian large-eddy simulation solver with hybrid MPI/OpenMP parallelisation

Eulerian-Lagrangian approaches capable of accurately reproducing complex fluid flows are becoming more and more popular due to the increasing availability and capacity of High Performance Computing facilities. However, the parallelisation of the Lagrangian part of such methods is challenging when a large number of Lagrangian markers are employed. In this study, a hybrid MPI/OpenMP parallelisation strategy is presented and implemented in a finite difference based large-eddy simulation code featuring the immersed boundary method which generally employs a large number of Lagrangian markers. A master-scattering-gathering strategy is used to deal with the handling of the Lagrangian markers and OpenMP is employed to distribute their computational load across several CPU threads. A classical domain-decomposition-based MPI approach is used to carry out the Eulerian, fixed-mesh fluid calculations. The results demonstrate that by using an effective combination of MPI and OpenMP the code can outperform a pure MPI parallelisation approach by up to 20%. Outcomes from this paper are of interest to various Eulerian-Lagrangian applications including the immersed boundary method, discrete element method or Lagrangian particle tracking

The impact of sequencing depth on the inferred taxonomic composition and AMR gene content of metagenomic samples

Shotgun metagenomics is increasingly used to characterise microbial communities, particularly for the investigation of antimicrobial resistance (AMR) in different animal and environmental contexts. There are many different approaches for inferring the taxonomic composition and AMR gene content of complex community samples from shotgun metagenomic data, but there has been little work establishing the optimum sequencing depth, data processing and analysis methods for these samples. In this study we used shotgun metagenomics and sequencing of cultured isolates from the same samples to address these issues. We sampled three potential environmental AMR gene reservoirs (pig caeca, river sediment, effluent) and sequenced samples with shotgun metagenomics at high depth (~ 200 million reads per sample). Alongside this, we cultured single-colony isolates of Enterobacteriaceae from the same samples and used hybrid sequencing (short- and long-reads) to create high- quality assemblies for comparison to the metagenomic data. To automate data processing, we developed an open- source software pipeline, ‘ResPipe’