2.1 Eulerian-Eulerian approach

The Eulerian-Eulerian approach is a powerful method for modeling multiphase flows. It treats each phase as a separate continuum, allowing for complex interactions between phases while considering their spatial and temporal evolution.

This approach uses conservation equations, constitutive relations, and averaging techniques to describe multiphase systems. It's widely applied in various fields, from chemical engineering to environmental science, despite computational challenges and modeling complexities.

Eulerian-Eulerian approach fundamentals

• The Eulerian-Eulerian approach is a framework for modeling multiphase flows where each phase is treated as a separate continuum and governed by its own set of conservation equations
• This approach allows for the modeling of complex interactions between phases, such as mass, momentum, and energy transfer, while considering the spatial and temporal evolution of each phase
• The Eulerian-Eulerian approach is particularly useful for systems with high volume fractions of dispersed phases or when the interactions between phases are significant

Continuum mechanics foundation

Top images from around the web for Continuum mechanics foundation

• 

  Time-integration of elastodynamics equation — Numerical tours of continuum mechanics using ... View original

  Is this image relevant?

• 

  Continuum mechanics - Wikipedia View original

  Is this image relevant?

• 

  Multiphase flow modeling of asphaltene precipitation and deposition | Oil & Gas Science and ... View original

  Is this image relevant?

• 

  Time-integration of elastodynamics equation — Numerical tours of continuum mechanics using ... View original

  Is this image relevant?

• 

  Continuum mechanics - Wikipedia View original

  Is this image relevant?

1 of 3

Top images from around the web for Continuum mechanics foundation

• 

  Time-integration of elastodynamics equation — Numerical tours of continuum mechanics using ... View original

  Is this image relevant?

• 

  Continuum mechanics - Wikipedia View original

  Is this image relevant?

• 

  Multiphase flow modeling of asphaltene precipitation and deposition | Oil & Gas Science and ... View original

  Is this image relevant?

• 

  Time-integration of elastodynamics equation — Numerical tours of continuum mechanics using ... View original

  Is this image relevant?

• 

  Continuum mechanics - Wikipedia View original

  Is this image relevant?

1 of 3

• Continuum mechanics provides the mathematical foundation for the Eulerian-Eulerian approach by treating each phase as a continuous medium with its own properties and field variables
• The principles of mass, momentum, and energy conservation are applied to each phase separately, while accounting for the interactions between phases
• Constitutive relations are used to describe the material behavior and close the system of equations

Interpenetrating continua concept

• In the Eulerian-Eulerian approach, phases are treated as interpenetrating continua, meaning that each phase can occupy the same physical space simultaneously
• This concept allows for the modeling of phase interactions and the exchange of mass, momentum, and energy between phases
• The interpenetrating continua assumption is valid when the length scales of the dispersed phases are much smaller than the characteristic length scales of the flow

Local instantaneous variables

• The Eulerian-Eulerian approach uses local instantaneous variables to describe the state of each phase at a given point in space and time
• These variables include density, velocity, temperature, and phase volume fraction, among others
• Local instantaneous variables are governed by the conservation equations and constitutive relations specific to each phase

Phasic volume fractions

• Phasic volume fractions represent the proportion of each phase present in a given control volume
• The sum of all phasic volume fractions must equal unity at any point in space and time
• Volume fractions are used to account for the presence of multiple phases and their influence on the flow behavior and interphase interactions

Conservation equations in Eulerian-Eulerian

• The Eulerian-Eulerian approach is based on the conservation of mass, momentum, and energy for each phase separately
• The conservation equations are derived by applying the principles of continuum mechanics to each phase and accounting for the interactions between phases
• The conservation equations form a coupled system of partial differential equations that describe the spatial and temporal evolution of each phase

Continuity equation for phases

• The continuity equation for each phase ensures the conservation of mass within the system
• It accounts for the temporal and spatial changes in the phasic density and velocity, as well as any mass transfer between phases
• The continuity equation for phase $k$ is given by: $\frac{\partial}{\partial t}(\alpha_k \rho_k) + \nabla \cdot (\alpha_k \rho_k \mathbf{v}_k) = \Gamma_k$

where $\alpha_k$ is the volume fraction, $\rho_k$ is the density, $\mathbf{v}_k$ is the velocity, and $\Gamma_k$ is the mass transfer term for phase $k$

Momentum equations for phases

• The momentum equations for each phase describe the conservation of momentum and the forces acting on the phase
• They account for the temporal and spatial changes in the phasic velocity, pressure gradient, viscous stresses, and interphase momentum transfer
• The momentum equation for phase $k$ is given by: $\frac{\partial}{\partial t}(\alpha_k \rho_k \mathbf{v}_k) + \nabla \cdot (\alpha_k \rho_k \mathbf{v}_k \mathbf{v}_k) = -\alpha_k \nabla p + \nabla \cdot \boldsymbol{\tau}_k + \alpha_k \rho_k \mathbf{g} + \mathbf{M}_k$

where $p$ is the pressure, $\boldsymbol{\tau}_k$ is the viscous stress tensor, $\mathbf{g}$ is the gravitational acceleration, and $\mathbf{M}_k$ is the interphase momentum transfer term for phase $k$

Energy equations for phases

• The energy equations for each phase ensure the conservation of energy within the system
• They account for the temporal and spatial changes in the phasic internal energy, heat flux, and interphase energy transfer
• The energy equation for phase $k$ is given by: $\frac{\partial}{\partial t}(\alpha_k \rho_k e_k) + \nabla \cdot (\alpha_k \rho_k e_k \mathbf{v}_k) = -\alpha_k \nabla \cdot \mathbf{q}_k + \alpha_k \rho_k \mathbf{g} \cdot \mathbf{v}_k + Q_k$

where $e_k$ is the specific internal energy, $\mathbf{q}_k$ is the heat flux vector, and $Q_k$ is the interphase energy transfer term for phase $k$

Interphase transfer terms

• Interphase transfer terms appear in the conservation equations to account for the exchange of mass, momentum, and energy between phases
• Mass transfer terms ($\Gamma_k$) represent the rate of mass exchange between phases due to processes such as evaporation, condensation, or chemical reactions
• Momentum transfer terms ($\mathbf{M}_k$) represent the forces acting between phases, such as drag, lift, and virtual mass forces
• Energy transfer terms ($Q_k$) represent the heat exchange between phases due to processes such as convection, conduction, or radiation

Constitutive relations

• Constitutive relations are essential in the Eulerian-Eulerian approach to close the system of conservation equations and describe the material behavior and interphase interactions
• These relations provide additional equations that relate the unknown variables to the known quantities and depend on the specific properties of the phases and the nature of their interactions
• Constitutive relations are often derived from empirical correlations, theoretical models, or experimental data

Interphase momentum transfer

• Interphase momentum transfer describes the forces acting between phases, which are crucial for accurately modeling the behavior of multiphase flows
• Drag force is the most significant interphase force and represents the resistance experienced by the dispersed phase as it moves through the continuous phase
• Other interphase forces include lift force (due to velocity gradients), virtual mass force (due to relative acceleration), and turbulent dispersion force (due to turbulent fluctuations)

Interphase heat and mass transfer

• Interphase heat and mass transfer are important in multiphase flows involving phase change, chemical reactions, or significant temperature gradients
• Heat transfer between phases can occur through convection (due to fluid motion), conduction (due to temperature gradients), or radiation (due to electromagnetic waves)
• Mass transfer between phases can occur through processes such as evaporation, condensation, dissolution, or chemical reactions

Granular flows and kinetic theory

• Granular flows are multiphase systems consisting of solid particles suspended in a fluid or gas
• The kinetic theory of granular flows is used to derive constitutive relations for the solid phase, accounting for particle collisions, friction, and energy dissipation
• The kinetic theory introduces additional variables, such as the granular temperature (a measure of the particle velocity fluctuations), and provides closure models for the solid phase stress tensor and interphase interactions

Turbulence modeling in multiphase

• Turbulence modeling is crucial in multiphase flows due to the complex interactions between phases and the presence of turbulent fluctuations
• Turbulence can be modeled using various approaches, such as the $k-\epsilon$ model, the $k-\omega$ model, or the Reynolds stress model, which are extended to account for the presence of multiple phases
• Turbulent dispersion force is an additional interphase force that arises due to the turbulent fluctuations and can be modeled using gradient diffusion or more advanced approaches

Averaging approaches

• Averaging approaches are used in the Eulerian-Eulerian framework to derive the conservation equations and constitutive relations for multiphase flows
• Averaging is necessary to handle the complex and chaotic nature of multiphase flows and to obtain tractable equations that describe the macroscopic behavior of the system
• Different averaging techniques can be employed depending on the specific characteristics of the flow and the desired level of detail in the model

Time averaging vs volume averaging

• Time averaging is used to derive the conservation equations for stationary or slowly evolving multiphase flows
• It involves averaging the instantaneous equations over a time interval that is much larger than the characteristic time scales of the flow (e.g., turbulence or bubble dynamics)
• Volume averaging is used to derive the conservation equations for spatially inhomogeneous multiphase flows
• It involves averaging the instantaneous equations over a representative volume that is much larger than the characteristic length scales of the dispersed phases (e.g., bubble or particle size)

Ensemble averaging and its applications

• Ensemble averaging is a statistical approach used to derive the conservation equations for multiphase flows with random or chaotic behavior
• It involves averaging the instantaneous equations over a large number of realizations or experiments with similar initial and boundary conditions
• Ensemble averaging is particularly useful for flows with significant turbulence or for systems with a large number of dispersed particles or bubbles

Favre averaging for compressible flows

• Favre averaging is a density-weighted averaging technique used for compressible multiphase flows
• It is defined as the mass-weighted average of a variable, which simplifies the conservation equations by eliminating the density fluctuation terms
• Favre averaging is particularly useful for high-speed multiphase flows, such as those encountered in aerospace applications or supersonic combustion

Numerical methods for Eulerian-Eulerian

• Numerical methods are essential for solving the complex system of partial differential equations that arise from the Eulerian-Eulerian approach
• The choice of numerical method depends on the specific characteristics of the flow, the desired accuracy, and the computational resources available
• Commonly used numerical methods for Eulerian-Eulerian simulations include finite volume, finite element, and finite difference methods

Finite volume discretization

• Finite volume method is a popular discretization technique for Eulerian-Eulerian simulations due to its conservation properties and ability to handle complex geometries
• The computational domain is divided into a number of control volumes, and the conservation equations are integrated over each control volume
• The fluxes across the control volume faces are approximated using interpolation schemes, such as upwind, central, or high-resolution schemes

Pressure-velocity coupling algorithms

• Pressure-velocity coupling algorithms are used to solve the coupled system of equations for pressure and velocity in incompressible or low Mach number flows
• Common pressure-velocity coupling algorithms include SIMPLE (Semi-Implicit Method for Pressure-Linked Equations), PISO (Pressure-Implicit with Splitting of Operators), and Fractional Step Method
• These algorithms involve an iterative procedure to obtain a converged solution for pressure and velocity fields

Implicit vs explicit time integration

• Time integration schemes are used to advance the solution in time and can be either implicit or explicit
• Explicit schemes (e.g., Forward Euler) are simple to implement but require small time steps to maintain stability, which can be computationally expensive
• Implicit schemes (e.g., Backward Euler, Crank-Nicolson) are more complex but allow for larger time steps and better stability properties
• The choice of time integration scheme depends on the specific problem, the required accuracy, and the computational resources available

Convergence and stability considerations

• Convergence and stability are critical aspects of numerical simulations and must be carefully considered when setting up and running Eulerian-Eulerian simulations
• Convergence refers to the property of the numerical solution approaching the exact solution as the mesh is refined or the time step is reduced
• Stability refers to the ability of the numerical scheme to produce bounded solutions without numerical oscillations or divergence
• Convergence and stability can be improved by using higher-order discretization schemes, appropriate time step sizes, and robust solution algorithms

Applications of Eulerian-Eulerian approach

• The Eulerian-Eulerian approach has been widely used to simulate various multiphase flow systems in different fields, such as chemical engineering, nuclear engineering, and environmental science
• The ability to model the complex interactions between phases and the spatial and temporal evolution of the flow makes the Eulerian-Eulerian approach suitable for a wide range of applications
• Some common applications of the Eulerian-Eulerian approach are discussed below

Gas-liquid flows and bubbly columns

• Gas-liquid flows are encountered in many industrial processes, such as bubble columns, airlift reactors, and aerated stirred tanks
• The Eulerian-Eulerian approach is used to model the distribution of gas bubbles in the liquid phase, the bubble coalescence and breakup, and the mass transfer between phases
• Bubbly columns are a specific type of gas-liquid flow system used for gas-liquid contacting, chemical reactions, and wastewater treatment

Fluidized beds and particle transport

• Fluidized beds are widely used in chemical and process industries for gas-solid contacting, heat transfer, and chemical reactions
• The Eulerian-Eulerian approach is used to model the fluidization behavior, the particle distribution, and the gas-particle interactions in fluidized beds
• Particle transport in pipelines and pneumatic conveying systems can also be modeled using the Eulerian-Eulerian approach, accounting for particle-particle and particle-wall interactions

Droplet and spray modeling

• Droplet and spray flows are important in various applications, such as fuel injection in combustion engines, spray drying, and agricultural spraying
• The Eulerian-Eulerian approach is used to model the droplet dispersion, evaporation, and coalescence in the carrier phase
• Advanced models can account for droplet breakup, collisions, and the effect of turbulence on droplet dynamics

Multiphase reacting flows

• Multiphase reacting flows involve chemical reactions between different phases and are encountered in various industrial processes, such as catalytic reactors, polymerization reactors, and combustion systems
• The Eulerian-Eulerian approach is used to model the species transport, heat transfer, and chemical reactions in multiphase reacting flows
• The interphase mass, momentum, and energy transfer, as well as the reaction kinetics, are accounted for in the model

Limitations and challenges

• Despite its wide applicability and success in modeling multiphase flows, the Eulerian-Eulerian approach has some limitations and challenges that need to be considered
• These limitations arise from the assumptions made in the model, the computational requirements, and the complexity of the physical phenomena involved
• Addressing these limitations and challenges is an active area of research in the multiphase flow community

Computational cost and scalability

• Eulerian-Eulerian simulations can be computationally expensive, especially for large-scale and three-dimensional problems
• The computational cost arises from the need to solve a large system of coupled partial differential equations for each phase and to account for the interphase interactions
• Parallel computing and high-performance computing techniques are often required to make Eulerian-Eulerian simulations feasible for practical applications

Modeling of complex physics

• The Eulerian-Eulerian approach relies on constitutive relations and closure models to describe the complex physical phenomena involved in multiphase flows
• These models are often based on empirical correlations or simplified theoretical assumptions and may not capture all the relevant physics accurately
• Modeling of turbulence, phase change, chemical reactions, and particle-particle interactions can be particularly challenging and may require advanced models and numerical techniques

Validation and verification strategies

• Validation and verification are essential steps in ensuring the accuracy and reliability of Eulerian-Eulerian simulations
• Validation involves comparing the simulation results with experimental data or analytical solutions to assess the model's ability to predict the real-world behavior
• Verification involves ensuring that the numerical implementation of the model is correct and that the discretization errors are within acceptable limits
• Comprehensive validation and verification studies are necessary to establish the credibility of Eulerian-Eulerian simulations for practical applications

Comparison with other multiphase approaches

• The Eulerian-Eulerian approach is one of several modeling approaches for multiphase flows, and its performance should be compared with other approaches to assess its strengths and weaknesses
• Other common approaches include the Eulerian-Lagrangian approach (where the dispersed phase is tracked as individual particles) and the Volume of Fluid (VOF) method (where the interface between phases is explicitly tracked)
• The choice of modeling approach depends on the specific characteristics of the flow, the desired level of detail, and the computational resources available
• A thorough understanding of the advantages and limitations of each approach is necessary to select the most appropriate method for a given multiphase flow problem