A spatially-averaged two-fluid model for dense large-scale gas-solid flows

AIChE Journal ◽  
2017 ◽  
Vol 63 (8) ◽  
pp. 3544-3562 ◽  
Author(s):  
Simon Schneiderbauer
Author(s):  
Aparicio Henriques ◽  
Pierre Coste ◽  
Sylvain Pigny ◽  
Jacques Magnaudet

An attempt to represent multiphase multiscale flow, filling the gap between Direct Numerical Simulation (DNS) and averaged approaches, is the purpose of this paper. We present a kind of Large Interface (LI) simulation formalism obtained after a filtering process on local instantaneous conservation equations of the two-fluid model which distinguishes between small scales and large scales contributions. LI surface tension force is also taken into account. Small scale dynamics call for modelisation and large scale for simulation. Joined to this formalism, a criterion to recognize LI’s is developped. It is used in an interface recognition algorithm which is qualified on a sloshing case and a bubble oscillation under zero-gravity. This method is applied to a rising bubble in a pool that collapses at a free surface and to a square-base basin experiment where splashing and sloshing at the free surface are the main break-up phenomena.


Author(s):  
F. J. Moraga ◽  
M. Lopez de Bertodano ◽  
D. A. Drew ◽  
R. T. Lahey

Two-fluid model simulations of a bubbly vertical boundary layer with point injection are presented. A new bubble turbulence dispersion model, designed to be used with RANS type turbulence models, was formulated and compared with recent data of [1] and [2]. These data showed that bubble migration toward the wall is controlled by the coherent large scale liquid structures within the boundary layer. The model is based on the application of a kinetic transport equation, similar to Boltzmann’s equation, and the idea that by selectively removing bubbles from the liquid eddies within the boundary layer, bubble capture at the wall introduces a preferential direction of migration and/or nonhomogeneous, anisotropic dispersion. This is the first model capable of predicting all the types of void fraction profiles observed experimentally for point injection. It is shown that without this new model, two-fluid model simulations fail to predict the experimental data. In addition, a new physical interpretation of the data of [1] and [2] is presented, which strongly suggests that the quantity controlling bubble migration toward the wall and bubble dispersion, is the boundary layer drift parameter (i.e., the ratio of the bubble’s terminal velocity to the free-stream liquid velocity).


Author(s):  
Yifei Duan ◽  
Zhi-Gang Feng

There are two primary approaches in modeling fluid-solid flows based on the method of treating particles suspended in flows. The first approach is the Eulerian-Lagrangian or Discrete Element Method (DEM) approach that tracks individual particles by solving the equations of motion of these particles. The second approach is the Eulerian-Eulerian approach or two-fluid model (TFM) that considers particles as another continuum phase or fluid. The TFM is preferred in modeling and predicting gas-solid flow behaviors in many engineering applications because of its efficiency in handling large-scale complex systems with large number of particles. However, one of the challenges in TFM is the uncertainty related to the selection of closure laws and transport properties of solid phases. In this study we employ the MFIX code, a general-purpose TFM computer code developed at the National Energy Technology Laboratory, to investigate the effect of different drag models and heat transfer models on the simulation results on the flow hydrodynamics and heat transfer of gas-solid fluidized beds. Three drag models (Gidspow model, Syamlal-O’Brien model, and Koch-Hill model) and two heat transfer model (Gunn model and a recently developed model by Sun et al., 2015) are tested. Simulation results from these models are compared with experimental measurements. The accuracy and applicability of these models are assessed and discussed.


2021 ◽  
Vol 502 (2) ◽  
pp. 2733-2749
Author(s):  
Siddhartha Gupta ◽  
Prateek Sharma ◽  
Andrea Mignone

ABSTRACT Cosmic rays (CRs) are frequently modelled as an additional fluid in hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations of astrophysical flows. The standard CR two-fluid model is described in terms of three conservation laws (expressing conservation of mass, momentum, and total energy) and one additional equation (for the CR pressure) that cannot be cast in a satisfactory conservative form. The presence of non-conservative terms with spatial derivatives in the model equations prevents a unique weak solution behind a shock. We investigate a number of methods for the numerical solution of the two-fluid equations and find that, in the presence of shock waves, the results generally depend on the numerical details (spatial reconstruction, time stepping, the CFL number, and the adopted discretization). All methods converge to a unique result if the energy partition between the thermal and non-thermal fluids at the shock is prescribed using a subgrid prescription. This highlights the non-uniqueness problem of the two-fluid equations at shocks. From our numerical investigations, we report a robust method for which the solutions are insensitive to the numerical details even in absence of a subgrid prescription, although we recommend a subgrid closure at shocks using results from kinetic theory. The subgrid closure is crucial for a reliable post-shock solution and also its impact on large-scale flows because the shock microphysics that determines CR acceleration is not accurately captured in a fluid approximation. Critical test problems, limitations of fluid modelling, and future directions are discussed.


2021 ◽  
Vol 33 (3) ◽  
pp. 033324
Author(s):  
Alejandro Clausse ◽  
Martín López de Bertodano

2021 ◽  
Vol 33 (3) ◽  
pp. 037116
Author(s):  
Victor L. Mironov

2021 ◽  
Vol 3 (4) ◽  
Author(s):  
R. Ponalagusamy ◽  
Ramakrishna Manchi

AbstractThe present communication presents a theoretical study of blood flow through a stenotic artery with a porous wall comprising Brinkman and Darcy layers. The governing equations describing the flow subjected to the boundary conditions have been solved analytically under the low Reynolds number and mild stenosis assumptions. Some special cases of the problem are also presented mathematically. The significant effects of the rheology of blood and porous wall of the artery on physiological flow quantities have been investigated. The results reveal that the wall shear stress at the stenotic throat increases dramatically for the thinner porous wall (i.e. smaller values of the Brinkman and Darcy regions) and the rate of increase is found to be 18.46% while it decreases for the thicker porous wall (i.e. higher values of the Brinkman and Darcy regions) and the rate of decrease is found to be 10.21%. Further, the streamline pattern in the stenotic region has been plotted and discussed.


Sign in / Sign up

Export Citation Format

Share Document