A Vortex Method for Bi-phasic Fluids Interacting with Rigid Bodies

We present an accurate Lagrangian method based on vortex particles, level-sets, and immersed boundary methods, for animating the interplay between two fluids and rigid solids. We show that a vortex method is a good choice for simulating bi-phase flow, such as liquid and gas, with a good level of realism. Vortex particles are localized at the interfaces between the two fluids and within the regions of high turbulence. We gain local precision and efficiency from the stable advection permitted by the vorticity formulation. Moreover, our numerical method straightforwardly solves the two-way coupling problem between the fluids and animated rigid solids. This new approach is validated through numerical comparisons with reference experiments from the computational fluid community. We also show that the visually appealing results obtained in the CG community can be reproduced with increased efficiency and an easier implementation

A vortex level set method for the two-way coupling of an incompressible fluid with colliding rigid bodies

International audienceWe present a vortex method for the simulation of the interaction of an incompressible flow with rigid bodies. The method is based on a penalization technique where the system is considered as a single flow, subject to the Navier-Stokes equation with a penalization term that enforces continuity at the solid-fluid interface and rigid motion inside the solid. Level set functions are used to capture interfaces, compute rigid motions inside the solid bodies and model collisions between bodies. A vortex in cell algorithm is built on this method. Numerical comparisons with existing 3D methods on problems of sedimentation and collision of spheres are provided to illustrate the capabilities of the method

Constitutive Relationships and Models in Continuum Theories of Multiphase Flows

In April, 1989, a workshop on constitutive relationships and models in continuum theories of multiphase flows was held at NASA's Marshall Space Flight Center. Topics of constitutive relationships for the partial or per phase stresses, including the concept of solid phase pressure are discussed. Models used for the exchange of mass, momentum, and energy between the phases in a multiphase flow are also discussed. The program, abstracts, and texts of the presentations from the workshop are included

Towards large eddy simulation of dispersed gas -liquid two-phase turbulent flows

This study presents a detailed investigation of all essential components of computational and modeling issues necessary for a successful large-eddy simulation (LES) of dispersed two-phase turbulent flows. In particular, a two-layer concept is proposed to enable the LES capability in two-phase flows involving dispersed bubbles that are relatively large compared to the mesh size. The work comprises three major parts.;Part I focuses on the development and verification of a transient, three-dimensional, finite-volume-method (FVM) based accurate Navier-Stokes solver, named DREAM II (second generation of the DREAM code). Several high-order schemes are implemented for both the spatial and temporal discretization. Solution of the coupled partial differential equations is attacked with a fractional step (projection) method. The developed solver is verified against various benchmarks including Taylor\u27s vortex, free-shear layer, backward-facing step flow and square cavity. A second-order overall accuracy is achieved in both space and time.;Part II concerns the modeling and LES of single-phase turbulent flows. A review of the LES theory and subgrid-scale (SGS) models is presented. Three SGS models, namely, Smagorinsky model, dynamic model and implicit model, are implemented and investigated. Then turbulent channel flow, plane mixing layer, and flow past a square cylinder are simulated, and comparisons of the first-, second-order statistics, and characteristic flow structures are made with direct numerical simulation (DNS) and/or benchmark experiments. The test results show superior quality of the present LES.;Part III delves into the theory, modeling and simulation of dispersed two-phase flow systems. A conceptual review of the characteristics and description of such system is made, considering both Eulerian-Eulerian (E-E) and Eulerian-Lagrangian (E-L) approaches, but with an emphasis on the latter. Various hydrodynamic forces acting on particles or bubbles are summarized and interpreted. Formulations regarding interphase coupling is discussed in depth. Typical computational treatments of modeled two-way couplings in an E-L DNS/LES are reviewed. Issues related to the interpolation are addressed. A general Lagrangian particle-tracking (LPT) program, named PART, is developed and verified using analytical solutions. (Abstract shortened by UMI.)

A generalised immersed boundary method for flows of dense suspension of solid particles

Immersed boundary method (IBM) provides computational advantages in approximating moving solid surfaces on fixed numerical meshes. It has been widely used for fully-resolved simulations of particulate flows. This thesis proposes a generalised formulation of IBM with improved applicability to flows with dense concentrations of particles and unstructured meshes. The new IBM formulation, which is based on the smooth-interface direct forcing approach, directly uses the algebraic discretised terms of the momentum equations in the evaluation of the forces on Lagrangian immersed boundary (IB) points, and evaluate the integral Lagrangian volumes based on these forces. Appropriate reconstructions of the boundary forces are adopted to ensure the compatibility with the momentum-weighted interpolation used for the finite-volume discretisation with a collocated mesh arrangement. A modified direct forcing formulation is also proposed, which results in an efficiency gain of a devised segregated flow-particle coupling scheme. The novel framework is applied to flows with stationary and moving IBs on both Cartesian and arbitrary triangular/tetrahedral meshes, and the results are similar or better than other related methods that are mostly developed for Cartesian meshes. Accurate and stable enforcement of the no-slip condition on the IB at every time-step is demonstrated, even for flows with strong transient behaviour and high velocity and pressure gradients. Local continuity in the vicinity of the IB is also preserved, ensuring local and global mass conservation alongside the local no-slip condition. Adaptations devised for unstructured meshes results in an accuracy close to that obtained on Cartesian meshes. The framework is successfully applied in the simulations of fluidisation of dense particle bed and a rising pack of light particles, showing robust stability. The issues related to the interfering regularised forces of different particle surfaces are not significant using the present formulation, hence eliminate unphysical flow patterns between aggregated particles.Open Acces

Proceedings of the 8th International Junior Researcher and Engineer Workshop on Hydraulic Structures

Full proceedings for the 8th International Junior Researcher and Engineer Workshop on Hydraulic Structures

Multiscale Modelling Of Platelet Aggregation

During clotting under flow, platelets bind and activate on collagen and release autocrinic factors such ADP and thromboxane, while tissue factor (TF) on the damaged wall leads to localized thrombin generation. Toward patient-specific simulation of thrombosis, a multiscale approach was developed to account for: platelet signaling (neural network trained by pairwise agonist scanning, PAS-NN), platelet positions (lattice kinetic Monte Carlo, LKMC), wall-generated thrombin and platelet-released ADP/thromboxane convection-diffusion (PDE), and flow over a growing clot (lattice Boltzmann). LKMC included shear-driven platelet aggregate restructuring. The PDEs for thrombin, ADP, and thromboxane were solved by finite element method using cell activation-driven adaptive triangular meshing. At all times, intracellular calcium was known for each platelet by PAS-NN in response to its unique exposure to local collagen, ADP, thromboxane, and thrombin. The model accurately predicted clot morphology and growth with time on collagen/TF surface as compared to microfluidic blood perfusion experiments. The model also predicted the complete occlusion of the blood channel under pressure relief settings. Prior to occlusion, intrathrombus concentrations reached 50 nM thrombin, ~1 μM thromboxane, and ~10 μM ADP, while the wall shear rate on the rough clot peaked at ~1000-2000 sec-1. Additionally, clotting on TF/collagen was accurately simulated for modulators of platelet cyclooxygenase-1, P2Y1, and IP-receptor. The model was then extended to a rectangular channel with symmetric Gaussian obstacles representative of a coronary artery with severe stenosis. The upgraded stenosis model was able to predict platelet deposition dynamics at the post-stenotic segment corresponding to development of artery thrombosis prior to severe myocardial infarction. The presence of stenosis conditions alters the hemodynamics of normal hemostasis, showing a different thrombus growth mechanism. The model was able to recreate the platelet aggregation process under the complex recirculating flow features and make reasonable prediction on the clot morphology with flow separation. The model also detected recirculating transport dynamics for diffusible species in response to vortex features, posing interesting questions on the interplay between biological signaling and prevailing hemodynamics. In future work, the model will be extended to clot growth with a patient cardio-vasculature under pulsatile flow conditions