Liquid-gas-solid flows with lattice Boltzmann: Simulation of floating bodies

This paper presents a model for the simulation of liquid-gas-solid flows by means of the lattice Boltzmann method. The approach is built upon previous works for the simulation of liquid-solid particle suspensions on the one hand, and on a liquid-gas free surface model on the other. We show how the two approaches can be unified by a novel set of dynamic cell conversion rules. For evaluation, we concentrate on the rotational stability of non-spherical rigid bodies floating on a plane water surface - a classical hydrostatic problem known from naval architecture. We show the consistency of our method in this kind of flows and obtain convergence towards the ideal solution for the measured heeling stability of a floating box.Comment: 22 pages, Preprint submitted to Computers and Mathematics with Applications Special Issue ICMMES 2011, Proceedings of the Eighth International Conference for Mesoscopic Methods in Engineering and Scienc

GPU-accelerated large-eddy simulation of ship-ice interactions

This paper reports on the applicability of the Lattice Boltzmann based free surface flow solver elbe to the simulation of complex ship-ice interactions in marine engineering. In order to model the dynamics of these colliding rigid multi-body systems, elbe is coupled to the ODE physics engine. First, basic validations of the ODE collision and friction models are presented, particularly focusing on interacting triangle meshes that later will serve to describe the ice floes. Then, the basic methodology and initial validation of the fluid-structure coupling of elbe and ODE is presented. Finally, performance is addressed: As elbe uses graphics processing units (GPUs) to accelerate the numerical calculations, the coupled numerical tool allows for investigations of ship-ice interactions in very competitive computational time and on off-the-shelf desktop hardware

Real-time physical engine for floating objects with two-way fluid-structure coupling

A method to simulate graphic animations of objects floating in a water surface in real time is presented. The fluid is simulated by means of the Lattice Boltzmann method for the shallow-waters equations, and the movement of the floating objects is calculated with a Newtonian physical engine suitable for the mechanics of rigid bodies. A two-way interaction between the fluid surface and the object structures is achieved by providing inputs to the Newtonian engine representing buoyancy, drag and lift forces calculated from the solution of the Lattice Boltzmann scheme, which in turn is perturbed by displacement forces acting at the objects boundaries. The method is tested in animation scenes of boats and different adrift objects, showing excellent rendering rates in desktop computers.Fil: Lazo, Marcos Gonzalo. Universidad Nacional del Centro de la Provincia de Buenos Aires; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Tandil; ArgentinaFil: Garcia Bauza, Cristian Dario. Universidad Nacional del Centro de la Provincia de Buenos Aires; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Tandil; ArgentinaFil: Boroni, Gustavo Adolfo. Universidad Nacional del Centro de la Provincia de Buenos Aires; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Tandil; ArgentinaFil: Clausse, Alejandro. Comision Nacional de Energia Atomica. Gerencia Quimica. CAC; Argentina. Universidad Nacional del Centro de la Provincia de Buenos Aires; Argentin

Direct simulation of liquid-gas-solid flow with a free surface lattice Boltzmann method

Direct numerical simulation of liquid-gas-solid flows is uncommon due to the considerable computational cost. As the grid spacing is determined by the smallest involved length scale, large grid sizes become necessary -- in particular if the bubble-particle aspect ratio is on the order of 10 or larger. Hence, it arises the question of both feasibility and reasonability. In this paper, we present a fully parallel, scalable method for direct numerical simulation of bubble-particle interaction at a size ratio of 1-2 orders of magnitude that makes simulations feasible on currently available super-computing resources. With the presented approach, simulations of bubbles in suspension columns consisting of more than $100\,000$ fully resolved particles become possible. Furthermore, we demonstrate the significance of particle-resolved simulations by comparison to previous unresolved solutions. The results indicate that fully-resolved direct numerical simulation is indeed necessary to predict the flow structure of bubble-particle interaction problems correctly.Comment: submitted to International Journal of Computational Fluid Dynamic

Link-wise Artificial Compressibility Method

The Artificial Compressibility Method (ACM) for the incompressible Navier-Stokes equations is (link-wise) reformulated (referred to as LW-ACM) by a finite set of discrete directions (links) on a regular Cartesian mesh, in analogy with the Lattice Boltzmann Method (LBM). The main advantage is the possibility of exploiting well established technologies originally developed for LBM and classical computational fluid dynamics, with special emphasis on finite differences (at least in the present paper), at the cost of minor changes. For instance, wall boundaries not aligned with the background Cartesian mesh can be taken into account by tracing the intersections of each link with the wall (analogously to LBM technology). LW-ACM requires no high-order moments beyond hydrodynamics (often referred to as ghost moments) and no kinetic expansion. Like finite difference schemes, only standard Taylor expansion is needed for analyzing consistency. Preliminary efforts towards optimal implementations have shown that LW-ACM is capable of similar computational speed as optimized (BGK-) LBM. In addition, the memory demand is significantly smaller than (BGK-) LBM. Importantly, with an efficient implementation, this algorithm may be one of the few which is compute-bound and not memory-bound. Two- and three-dimensional benchmarks are investigated, and an extensive comparative study between the present approach and state of the art methods from the literature is carried out. Numerical evidences suggest that LW-ACM represents an excellent alternative in terms of simplicity, stability and accuracy.Comment: 62 pages, 20 figure