^{1}

^{2}

^{2}

^{*}

From results obtained in the rheological characterization of a 4% dispersion of flamboyant gum with the HaakeRT20 viscometer, for different conditions of pH = 3.0 and 9.0, temperature 5°C, 25°C and 45°C, the data of the rheological behavior of the gum dispersions were fitted to the power law model. To understand and predict the behavior of this gum, a model of Lattice Boltzmann D2Q9 was developed for the behavior, in addition to simulations for the conditions handled in the experiments performed with the HaakeRT20 viscometer.

Gums or hydrocolloids are polysaccharides derived from plants, among which are guar, carob, gum Arabic from Acacia trees, the pectin of the citrus peel and the apple pulp. In the case of carrageenan and alginates which are extracted from marine algae, xanthan gum is produced by microbial fermentation [

Polysaccharides are used in many industrial processes such as thickeners, gelling agents, emulsifiers, adherents, etc. Its high molecular weight and the interactions between the polymer chains allow altering rheological properties such as viscosity when dissolved or dispersed [

Rheology relates the behavior of materials when subjected to forces and deformations. These concepts are of great importance in rheological evaluations. The application areas are also extremely wide and diverse and require an important contribution of Engineers with extensive experience, although the Chemical and Process Engineers, by their role in the handling and processing of complex materials (such as foams, suspensions, emulsions, molten polymers and solutions, etc.), have a huge interest in this subject [

There are factors that can influence the rheological properties such as, variety, composition, temperature, instrumental techniques among others.

The rheometer or viscometer is used to measure the rheological properties of fluids under the principle of resistance to flow by the action of a known force or stress produced by a known amount of flow. Between the viscometers used to measurement the properties, are the capillaries and the descending ball, in the case of rheometers, the rotary and oscillatory rheometer is the more used. Commonly in the rotational rheometers, the controlled rate and controlled stress approach are used. In the first approach, the material under study is placed between two plates, one of them is rotated at a fixed speed and the other plate measures the torsional force. In the second approach, the situation is reversed [

There are investigations to know the rheological properties of solutions of guar gum, xanthan and alginate with the addition of salep for different concentrations, using the power law model [

Other investigations to evaluate the rheological properties for food hydrocolloids (carrageenan, pectin, gelatin, starch and xanthan) using a rotational viscometer for different concentrations and temperatures. Finding that mainly the model that was most adjusted to most of the experiments was the power model, and to describe the effect of temperature the Arrhenius model was used. To evaluate the effect of the concentration on the apparent viscosity we used the power, exponential and polynomial model [

The effect of calcium salt and pH on the rheological properties of the solutions of mixtures of xantham-carboxymethylcellulose has also been investigated. The gum solutions in all proportions showed a behavior to the Herschel-Bulkley model. In their research they found that the addition of salt and the decrease in pH cause an apparent reduction in viscosity [

Other studies focus on the rheological properties of the exudate of bitter almond gum (BAG) at different shear rates, concentrations, temperatures and pH in the presence of several salts. The rheological data was adjusted with the Power law model. It was found that the apparent viscosity of the BAG solutions increased with the increase in gum concentrations and decreased with the increasing shear rate at a specific temperature. The salts caused a reduction in viscosity. All the treatments had significant effects on the rheological parameters [

Our investigations were based on Flamboyant seeds (Delonix regia), which were harvested from trees in the city of Cuernavaca, Morelos, Mexico. The flamboyant’s gum (FG) was extracted according to Corzo-Rios, et al., 2014 [

Subsequently, a rheological characterization of a 4% dispersion of flamboyant gum was carried out with the HaakeRT20 viscometer equipment, for different conditions of pH = 3.0 and 9.0 and temperatures 5.0˚C, 25˚C and 45˚C. Obtaining that the behavior found corresponds to a non-Newtonian fluid that was adjusted to that of the Power law.

The Lattice Boltzmann Method has been presented in many investigations as an alternative to address problems concerning computational fluid dynamics, which has proven to be an effective tool in the past two decades.

In works concerning the application of rheology, the droplet generation in Newtonian and non-Newtonian fluid in T-junction and cross-junction microchannels was simulated by the Lattice Boltzmann method. The results revealed that it is necessary to take into account non-Newtonian rheology instead of the simple assumption of Newtonian fluid in numerical simulations [

They have also been investigated the effects of capillary number, viscosity ratio, gravitational acceleration, wettability, geometry, and non-Newtonian rheology property on finger pattern using the Lattice Boltzmann Method. The results provide an understanding of the viscous fingering phenomenon and confirm that the Lattice Boltzmann Method is an effective tool to investigate the immiscible displacement in complex geometries, as well as the non-Newtonian multiphase fluid flow behavior problems [

The LBM model with color-gradient has been used to simulate immiscible two-phase flows with power-law rheology. The deformation of a Newtonian droplet in a power law matrix subject to a shear flow, and the influence of power-law index, number capillary, wall confinement, and viscosity ratio. In addition to investigating the effect of non-Newtonian rheology on the droplet breakup under different confinements, and the results of droplet breakup in Newtonian systems are also shown for comparison with available in the literature data [

The velocity profiles of Newtonian and non-Newtonian fluids in microchannels have been investigated by the micro-particle image velocimetry technique (micro-PIV). Polyacrylamide aqueous solutions with different concentrations were used as non-Newtonian fluids. The comparison of the experimentally measured velocity profiles with the use of micro-PIV with the Lattice Boltzmann simulation and an approximate solution, showed acceptable results, which could be an alternative to the use of rheometers [

In this work, an LBM D2Q9 was developed, in the case of a power fluid, in addition to performing computational experiments for the pH and temperature conditions handled in the experiments done with the Haake RT20 viscometer, with the purpose of predicting the behavior of a flamboyant gum with the LBM.

Two of the most important equations in fluid mechanics are continuity and momentum. The first is the basic principle of the movement of fluids and requires that the mass be conserved, which is given by

∂ ρ ∂ t + ∇ ⋅ ρ u = 0 (1)

In Equation (1), ρ is the density and u is the velocity of the fluid. The second refers to the movement of the fluid, as well as the forces that produce these movements

ρ ( ∂ u ∂ t + u ⋅ ∇ u ) = − ∇ ⋅ τ ¯ ¯ − ∇ P + ρ g (2)

In the case of a Newtonian fluid, the ratio of the shear stress τ α β and the strain rate tensor e ˙ α β is given by the following relationship

τ α β = 2 μ e ˙ α β ; e ˙ α β = 1 2 ( ∂ u α ∂ x β + ∂ u β ∂ x α ) (3)

In this relation, the viscosity μ does not depend on the e ˙ α β , since μ remains constant. Incorporating this relation in Equation (2), we obtain the Navier Stokes equation

ρ ( ∂ u ∂ t + u ⋅ ∇ u ) = μ ∇ 2 u − ∇ P + ρ g (4)

For a non-Newtonian fluid, the power or model Ostwald of Waele is used, in which, the ratio of τ α β and e ˙ α β is given by

τ α β = 2 μ e f f e ˙ α β ; μ e f f = k e ˙ n − 1 ; e ˙ = e ˙ α β e ˙ α β (5)

In Equation (5), µ_{eff} is the effective viscosity, k and n are the rheological parameters corresponding to the consistency index and the fluid behavior index respectively. For the case of n = 1.0, it is a Newtonian fluid. If n < 1.0, it corresponds to a pseudoplastic liquid, in which, its effective viscosity decreases with e ˙ α β . If n > 1.0, corresponds to a dilatant liquid, in these, the effective viscosity increases with e ˙ α β [

This method evaluates the evolution of the distribution function of particle f i = ( x , t ) by means of the Boltzmann equation discretized in a lattice. That is the probability of finding a particle with velocity e i at position x and time t. The method can be summarized in two stages, the first is the advance of the particles to the next lattice site along the directions of motion for each time step Δt. The second stage is to simulate the collisions of the particles. The two mentioned stages can be described through Boltzmann’s discrete equation as

f i ( x + e i , t + 1 ) − f i ( x , t ) = − Ω i (6)

The collision term Ω_{i} that appears on the right side of Equation (6), can be replaced by

f i ( x + e i , t + 1 ) − f i ( x , t ) = − 1 τ [ f i ( x , t ) − f i e q ( x , t ) ] (7)

This simplification is known as the BGK approach (Bhatnagar-Gross-Krook). The rate of change referred to the equilibrium is 1/τ, which is the inverse of the relaxation time.

The macroscopic variables can be calculated in a direct way from the values of the distribution function as

ρ ( x , t ) = ∑ i = 1 n f i ( x , t ) u ( x , t ) = 1 ρ ( x , t ) ∑ i = 1 n f i ( x , t ) e i (8)

In the case of a Newtonian fluid, the relation of the kinematic viscosity and the relaxation time is given by υ = 1 / 3 ( τ − 1 / 2 ) .

One of the contributions of the method presented in this work involves a modification of the LBGK presented by Aharonov and Rothman [

τ = 1 2 + 3 k ( 2 e ˙ ) n − 1 (9)

In this work two types of boundaries were used, the “bounce-back” in the walls to ensure that the speeds in these are zero, and the periodic boundary conditions, in which the nodes located in the border will have their neighboring nodes in the opposite border, with a lattice size of 80 × 80. The values of “n” used were 0.2, 0.3, 0.4, 0.5, 0.6 and 0.7, with “k” values of 6.5 × 10^{−5}, 9.0 × 10^{−5}, 9.0 × 10^{−5}, 9.0 × 10^{−5}, 6.5 × 10^{−3} and 6.5 × 10^{−3} respectively. The steady state was reached after 90,000 steps of time.

The validation of the LBM was performed using the analytical solution of the plane Poiseuille flow, in dimensionless form for a power fluid

u 〈 u 〉 = ( A ) [ 1 − ( B − y C ) 1 n + 1 ] A + 1 ( C ) 1 n + 1 ( n 2 n + 1 ) [ − ( B ) 1 n + 2 − ( C ) 1 n + 2 ] (10)

In Equation (10), 〈 u 〉 is the average speed, A, B and C are constants related to the flow dimensions between two plates.

The comparison of the analytical velocity profiles Equation (10) and those performed with the LBM showed a relative error rate between both velocity profiles of less than 0.8%.

Simulation of Pseudoplastic Fluids with the Lattice Boltzmann Method with Experimental Data of Flamboyant GumsIn this work, the results obtained from the rheological characterization of flamboyant gum were used. The equipment used for these characterizations was the HaakeRT20 Viscometer, obtaining the conditions, as well as the results presented in

In

Flamboyant gum, 4.0%, pH = 9.0 @ 5.0˚C | Flamboyant gum, 4.0%, pH = 3.0 @ 5.0˚C | ||||||
---|---|---|---|---|---|---|---|

Sample A | Sample B | Sample A | Sample B | ||||

n = 0.5 | k = 16.06 | n = 0.52 | k = 16.32 | n = 0.26 | k = 201.59 | n = 0.22 | k = 286.51 |

Flamboyant gum, 4.0%, pH = 9.0 @ 25˚C | Flamboyant gum, 4.0%, pH = 3.0 @ 25˚C | ||||||

Sample A | Sample B | Sample A | Sample B | ||||

n = 0.4708 | k = 13.8158 | n = 0.4961 | k = 12.0731 | n = 0.2653 | k = 203.1150 | n = 0.2345 | k = 287.1152 |

Flamboyant gum, 4.0%, pH = 9.0 @ 45˚C | Flamboyant gum, 4.0%, pH = 3.0 @ 45˚C | ||||||

Sample A | Sample B | Sample A | Sample B | ||||

n = 0.5 | k = 7.84 | n = 0.46 | k = 13.47 | n = 0.24 | k = 264.18 | n = 0.25 | k = 212.78 |

In

In ^{−1}. In addition, it can be seen that the increase in shear stress is very pronounced which produces higher values of consistency index in the samples with pH = 3.0 than at pH = 9.0, which makes it difficult to determine the rheological behavior at higher shear rates at 100 s^{−1}.

An important aspect in the simulation of real systems is the correct transformation of units; that is, conversion of physical units to lattice units and vice versa. To do this, conversions were made of the experimental values of k_{phys} that have physical units with respect to the grid units k_{LB}, using the velocity of physical light and Lattice Boltzmann of 340.3 m/s and 1 / 3 respectively [

In

In

The shear stress and the shear rate were obtained through the analytical solution, with the simulations of the LBM and compared with the experimental values obtained with the HaakeRT20 viscometer.

Sample | n | k_{phys} (Pa∙s^{n}) | k_{LB} (u.r) | |
---|---|---|---|---|

Flamboyant gum, 4.0%, pH = 9.0 @ 5.0˚C | A | 0.50 | 16.06 | 8.1486 × 10^{−5} |

B | 0.52 | 16.32 | 1.0186 × 10^{−5} | |

Flamboyant gum, 4.0%, pH = 9.0 @ 25˚C | A | 0.47 | 13.82 | 5.1807 × 10^{−5} |

B | 0.49 | 12.07 | 5.8830 × 10^{−5} | |

Flamboyant gum, 4.0%, pH = 9.0 @ 45˚C | A | 0.50 | 7.84 | 3.9780 × 10^{−5} |

B | 0.46 | 13.47 | 4.5165 × 10^{−5} | |

Flamboyant gum, 4.0%, pH = 3.0 @ 5˚C | A | 0.26 | 201.59 | 8.5197 × 10^{−05} |

B | 0.22 | 286.51 | 8.0020 × 10^{−05} | |

Flamboyant gum, 4.0%, pH = 3.0 @ 25˚C | A | 0.27 | 203.12 | 9.0685 × 10^{−05} |

B | 0.23 | 287.12 | 9.3181 × 10^{−05} | |

Flamboyant gum, 4.0%, pH = 3.0 @ 45˚C | A | 0.24 | 264.18 | 9.0763 × 10^{−05} |

B | 0.25 | 212.78 | 8.1080 × 10^{−05} |

gum. The approximations are acceptable for the case of Figures 4(a)-(f), with the conditions of pH = 9.0, for the temperatures of 5.0˚C, 25˚C and 45˚C. However, for

A rheological characterization of the flamboyant gum was performed with the HaakeRT20 viscometer, for different conditions of pH = 3.0 and 9.0, temperature 5.0˚C, 25˚C and 45˚C, and it was found that this gum exhibited a behavior of power law of shear thinning type. The experimental data were compared with analytical solutions and with simulations carried out with the Lattice Boltzmann Method, which showed acceptable results, except for values of pH = 3.0 and temperature of 5.0˚C where differences can be observed between the 2 samples.

The Lattice Boltzmann Method, turns out to be an alternative to the conventional methods of computational fluid mechanics, in addition to presenting facility in terms of implementation and programming. The technique promises to be a predictive tool in the area of rheology.

The present work was developed under the sponsorship of CONACYT, IPN and FESC-UNAM.

The authors declare no conflicts of interest regarding the publication of this paper.

Corzo-Ríos, L.J., Ortiz, M.E.R. andVelázquez-Ortega, J.L. (2019) Simulation of Flamboyant Gums (Power Fluid) with the Lattice Boltzmann Method. Open Journal of Fluid Dynamics, 9, 193-205. https://doi.org/10.4236/ojfd.2019.92013