Multiphase simulations

Simulations of multiphase flow are, in general, very poor, with a few exceptions. Basically, there are three different kinds of multiphase models Euler-Lagrange, Euler-Euler, and volume of fluid (VOF) or level-set methods. The Euler-Lagrange and Euler-Euler models require that the particles (solid or fluid) are smaller than the computational grid and a finer resolution below that limit will not give a 

Multiscale modeling is an approach to minimize system-dependent empirical correlations for drag, particle-particle, and particle-fluid interactions [19]. This approach is visualized in Eigure 15.6. A detailed model is developed on the smallest scale. Direct numerical simulation (DNS) is done on a system containing a few hundred particles. This system is sufficient for developing models for particle-particle and particle-fluid interactions. Here, the grid is much smaller 

VOF or level-set models are used for stratified flows where the phases are separated and one objective is to calculate the location of the interface. In these models, the momentum equations are solved for the separated phases and only at the interface are additional models used. Additional variables, such as the volume fraction of each phase, are used to identify the phases. The simplest model uses a weight average of the viscosity and density in the computational cells that are shared between the phases. Very fine resolution is, however, required for systems when surface tension is important, since an accurate estimation of the curvature of the interface is required to calculate the normal force arising from the surface tension. Usually, VOF models simulate the surface position accurately, but the space resolution is not sufficient to simulate mass transfer in liquids. 

For fluid particles that continuously coalesce and breakup and where the bubble size distributions have local variations, there is still no generally accepted model available and the existing models are contradictory [20]. A population density model is required to describe the changing bubble and drop size. Usually, it is sufficient to simulate a handful of sizes or use some quadrature model, for example, direct quadrature method of moments (DQMOM) to decrease the number of variables. 

Owing to the high computational load, it is tempting to assume rotational symmetry to reduce to 2D simulations. However, the symmetrical axis is a wall in the simulations that allows slip but no transport across it. The flow in bubble columns or bubbling fluidized beds is never steady, but instead oscillates everywhere, including across the center of the reactor. Consequently, a 2D rotational symmetry representation is never accurate for these reactors. A second problem with axis symmetry is that the bubbles formed in a bubbling fluidized bed are simulated as toroids and the mass balance for the bubble will be problematic when the bubble moves in a radial direction. It is also problematic to calculate the void fraction with these models. 

---

Hence the missing baryons could be in the WHIM, which would be correspondingly enriched. Unfortunately such strong winds are not supported by hydrodynamical simulations. These use supernovae as the energy source that drives the wind. However the current multiphase simulations lack sufficient fine-scale resolution, as discussed below. 

Ranade and Van den Akker [74], for example, used the snapshot method for simulating gas-liquid flows in baffled stirred tanks using a time after volume averaged two-fluid model for incompressible flows (as described in sect 3.3). These multiphase simulations also predicted the near-impeller flows with fair accuracy. Most important, the cavities due to the accumulation of gas in the low-pressure region behind the impeller blades were detected. 

Tacite multiphase simulator for complex transient flow phenomena. 

Theoretical and experimental studies demonstrated that segmented gas/liquid flow enhances the heat transfer considerably compared to a single flow. Multiphase flow simulations revealed mainly two mechanisms explaining the increase of the Nusselt number, namely the circulation within the liquid slugs and the disturbing and renewing of the fluid layer near the wall by the gas bubbles. Detailed multiphase simulations for cylindrical chaimels by Lakehal et al. [19] lead to the... 

Where new models are being developed, a qualitatively reasonable result is the first requirement, and convergence criteria are typically looser. Luo et al., who reported an early sliding mesh simulation in a short note, used the flow pattern became cyclically repeatable as their convergence criterion [13]. Gosman et al., who report early multiphase simulations, simply required that residuals in the equations solved become smaller than a prescribed tolerance. [9]... 

Previous chapters have discussed both polymer degradation and adsorption in some detail but not in the context of oil recovery. Earlier in this chapter, it was indicated how terms describing these effects are incorporated into the multiphase simulation equations for polymer flooding. Here, some calculations which illustrate and quantify the effects of both polymer adsorption and degradation in both cross-sectional and areal reservoir models will be presented. In connection with polymer degradation, the effects of local cooling as described in Section 8.5.5 and shown in Figure 8.13 will also be discussed. 

Computation of solution is commenced by iterating at 0.01 second time step size for 4000 time steps (for 1mm fabric thickness) and 2000 time steps (for 1.5 mm fabrie thiekness) with 3 iterations per each time step. The solution converged to the defined toleranees around 500-550 iterations for all multiphase simulation of RTM resin flow in textile fabries. A representative residual convergenee plot is as shown by Figure 13.5 [1]. Animated frame of resin volume fraction for the illustrated example at 1 seeond and 5 second are as shown in Figure 13.6 [1]. Figures 13.7 and 13.8 [1] provide the resin mass flow rate plots with respect to flow time at inlet and outlet, respeetively. 

Multiphase and nonequilibrium simulations are extremely difficult. These usually entail both a large amount of computing resources and a lot of technical expertise on the part of the researcher. Readers of this book are urged to refer such projects to specialists in this area. 

Knowledge of these types of reaetors is important beeause some industrial reaetors approaeh the idealized types or may be simulated by a number of ideal reaetors. In this ehapter, we will review the above reaetors and their applieations in the ehemieal proeess industries. Additionally, multiphase reaetors sueh as the fixed and fluidized beds are reviewed. In Chapter 5, the numerieal method of analysis will be used to model the eoneentration-time profiles of various reaetions in a bateh reaetor, and provide sizing of the bateh, semi-bateh, eontinuous flow stirred tank, and plug flow reaetors for both isothermal and adiabatie eonditions. 

J. Tiaden, B. Nestler, H. J. Diepers, I. Steinbach. Physica D 115 11, 1998 G. J. Schmitz, B. Nestler. Simulation of phase transitions in multiphase systems, peritectic solidification of YBaCuO-superconductors. Mater Sci Eng B 53 11, 1998. 

We have considered thermodynamic equilibrium in homogeneous systems. When two or more phases exist, it is necessary that the requirements for reaction equilibria (i.e., Equations (7.46)) be satisfied simultaneously with the requirements for phase equilibria (i.e., that the component fugacities be equal in each phase). We leave the treatment of chemical equilibria in multiphase systems to the specialized literature, but note that the method of false transients normally works quite well for multiphase systems. The simulation includes reaction—typically confined to one phase—and mass transfer between the phases. The governing equations are given in Chapter 11. 

Mathpati, C.S. and Joshi, J.B. (2007) Insight into theories of heat and mass transfer at the solid/fluid interface using direct numerical simulation and large eddy simulation. Joint 6th International Symposium on Catalysis in Multiphase Reactors/5th International Symposium on Multifunctional Reactors (CAMURE-6/ISMR-5-), 2007, Pune. 

Y. Inoue, Y. Chen, and H. Ohashi, A mesoscopic simulation model for immiscible multiphase fluids, J. Comput. Phys. 201, 191 (2004). 

Zhang K., Moridis G., et al. TOUGH+C02 A multiphase fluid-flow simulator for C02 geologic sequestration in saline aquifers. 2011 Computers Geosdnces 37 714— 723. 

Pruess K. TOUGH2 A genaral simulator for multiphase fluid and heat flow. 1991 Lawrence Berkeley National Laboratory Report, LBNL, California 29400. 

Helmig R., Class H., et al. Architecture of the modular program system MUFTE-UG for simulating multiphase flow and transport processes in heterogeneous prous media. 1998 MathematischeGeologie2 123-131. 

Kanai, A., and Mtyata, H. Numerical simulation of bubbles in a boundary layer by Maker-Density-Function . Proceedings of the 3rd International Conference on Multiphase Flow, Lion, France (1998). 

This review paper is restricted to stirred vessels operated in the turbulent-flow regime and exploited for various physical operations and chemical processes. The developments in the field of computational simulations of stirred vessels, however, are not separated from similar developments in the fields of, e.g., turbulent combustion, flames, jets and sprays, tubular reactors, and multiphase reactors and separators. Fortunately, there is a strong degree of synergy and mutual cross-fertilization between these various fields. This review paper focuses on aspects specific to stirred vessels (such as the revolving impeller, the resulting strong spatial variations in turbulence properties, and the macroinstabilities) and on the processes carried out in them. 

---