Mathematical Problems in Engineering
Volume 2021 (2021), Article ID 5566754, 10 pages
https://doi.org/10.1155/2021/5566754

Research Article

Test Method for the Solubility Model of Physical Blowing Agent of Self-Expanding Polymer

Heyang Jia, Xiaolong Li, Meimei Hao, Yang Li, Yanhui Zhong, and Bei Zhang

School of Water Conservancy Engineering, Zhengzhou University, 100 Science Avenue, Zhengzhou City 450001, China

Correspondence should be addressed to Xiaolong Li; wennuandeshang@hotmail.com

Received 12 January 2021; Revised 24 February 2021; Accepted 10 March 2021; Published 26 April 2021

Academic Editor: Isabella Torcicollo

Copyright © 2021 Heyang Jia et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper aims to present a solubility model of physical blowing agent (PBA) for a kind of commonly used self-expanding polymer on engineering. The self-expanding polymer contains Component A (isocyanate) and Component B (polyhydric alcohol, PBA, water, and catalyst). Component B grout of the polymer, which contains PBA, was heated to measure its temperature and volume variations. Based on the principle of mass conservation and Clapeyron equation, the solubility curve of PBA with respect to temperature was calculated. The solubility curve was then applied to simulate the foaming process. A two-component polymer grout foaming experiment was then carried out to verify the applicability of the measured solubility curve. The simulated changes of temperature and density with respect to time of polymer grout were analyzed and compared with experimental results. The error of both sets of curves is within 5%, which shows a good agreement among them and proves the feasibility of the solubility model. This study provides a relatively complete test and verification process for the solubility model of PBA, which lays a theoretical foundation for both the polymer grouting diffusion mechanism and engineering application.

1. Introduction

Grouting is an important method for antiseepage and reinforcement of engineering rock mass. It has been widely used in water conservancy, transportation, mining, and other fields. In recent years, a self-expanding nonaqueous reactive polyurethane grouting material (namely, polymer) and its high-pressure injection technique have developed rapidly [1–5]. Polymer material is a popular process for highway maintenance, tunnel disengaging repair, dam reinforcement, and other engineering fields due to its early strength, high expansibility, and good durability advantages. Polymer grouting has become one of the most characteristic development directions in geotechnical engineering field [6–9].

The research on the diffusion mechanism of polymer grout is significant for grouting design and construction. Compared with the constant density grout such as cement and sodium silicate, which are driven completely by static pressure, the grouting mechanism of self-expanding polymer is more complex. After the polymer grout is injected, chemical reactions occur rapidly. As the reaction goes on, a large amount of heat releases. The PBA gradually vaporizes with the increase of grout temperature, foaming a large number of tiny closed-cell bubbles suspended in the grout liquid. This results in the continuous expansion of grout. The grout was driven to flow as well. Numerical simulation on the chemical reaction process of polymer grout is an important way to study its expansion and diffusion mechanism. Obtaining the PBA solubility is one of the prerequisites for numerical solution of polymer grout expansion process.

Some scholars have conducted researches on the solubility model of PBA. Qian et al. [10] proposed a Flory-Huggins solubility model. Tesser et al. [11–14] determined the corresponding parameters of PBA such as CF3CH2F, CF3CHF2, CF3CH3, and CH2F2. Gupta and Khakhar [15, 16] established a solubility empirical formula for n-pentane PBA according to a solubility model given by V. K. Gupta. A. J. Rojas et al. presented a temperature-related solubility model. Rojas et al. [17, 18] presented relevant parameters of this solubility model for CFC-11a and HCFC-141b. These papers do not provide the specific test and verification process of the solubility model of PBA; instead, they just give the solubility model and its parameters directly. Therefore, it is necessary to conduct complete test and verification process in combination with different grout types.

In this paper, a two-component polyurethane engineering grouting material was taken as the study object, which involved the following steps:

(1)Component B grout was heated, and the solubility of the PBA of polymer was tested according to the temperature and volume of Component B grout(2)The solubility model was applied to the simulation of chemical reaction process of polymer; the simulated temperature and density curves with time are obtained(3)Carry out the foaming experiment of self-expanding polymer grout to obtain its temperature and density changes with respect to time(4)The applicability of the measured solubility curve was verified by comparing the simulated value in (2) with the experimental one in (3)

Figure 1 shows the general sketch of this research.

Figure 1: General sketch of this research.

2. Chemistry Reaction Mechanism

The self-expanding polymer material is mainly constituted of isocyanate, polyhydric alcohol, blowing agent, water, and catalyst. These materials are combined with two components, called Components A and B. Component A consists of isocyanate; Component B consists of polyhydric alcohol, PBA, water, and catalyst.

The blowing agent can be divided into physical blowing agent and chemical blowing agent according to the blowing principle. Water is the common chemical blowing agent of polymer, which reacts with isocyanate to produce carbon dioxide and promote grout polymerization as well. The commonly used physical blowing agents are CFC-11, CFC-13, HCFC-141b, C5H10, and so forth. With the increase of the grout temperature, the PBA dissolved in the polymer grout gradually vaporizes. A large number of tiny bubbles are generated and dispersed in the grout. This causes continuous grout expansion. When studying the polymer diffusion mechanism, it is necessary to calculate the number of bubbles during the foaming process in order to accurately simulate the grout expansion process [19, 20]. Therefore, establishing an accurate solubility model of PBA in polymer grout is needful.

The chemical reaction of self-expanding polymer grout is complex, mainly including two primary reactions: gelling reaction and blowing reaction. Gelling reaction refers to the reaction between isocyanate and polyhydric alcohol to produce polyurethane, which is as follows:

Ignoring the effect of molecular diffusion, its chemical reaction kinetic equation is

Blowing reaction exists when water is used for the chemical blowing agent of polymer grout. Water reacts with isocyanate to form carbon dioxide, and the reaction equation is

The corresponding chemical reaction kinetic equation is

where AOH and AW are the preexponential factor, EOH and EW are the activation energy, ci,0 is the initial concentration of each component, ci is the current concentration of each component, and Xi is their fractional conversion: , , and . R is the universal gas constant, T is the grout temperature, ri is the mass fraction of each component, and ρi is their density [21, 22].

Both gelling reaction and blowing reaction are exothermic reactions. With the progress of chemical reaction, a large amount of heat energy was produced, which causes continuous temperature increase of grout system.

The polymer satisfies the following energy balance during the reaction process:

where Ci is the heat capacity of each component, (−ΔH)i represents the heat of reaction, and λ is the heat of evaporation.

On the macro level, the polymer grout can be regarded as a homogeneous fluid composed of three liquids and two gases. The liquid part consists of water, polyurethane, and liquid PBA. The gas part consists of gaseous PBA and carbon dioxide. The density of the polymer grout at any time is given by

where ρBL is the density of the PBA in the liquid phase, ρW is the density of water, and ρP is the density of polyurethane. rBL,0 and rW,0 are the initial mass fractions of liquid PBA and water, respectively, rBL is the current mass fraction of liquid PBA (namely, solubility of liquid PBA), rBG is the current mass fraction of gaseous PBA, and the relationship between rBG and rBL is

where rW and rCO2 are the mass fractions of water and carbon dioxide, respectively. They can be calculated from concentration and fractional conversion, which are as follows:

where Mi is the molecular weight and rCO2,D is the initial mass fraction of carbon dioxide.

It can be seen from equation (6) that the density equation is closely related to rBG and rBL. This further proved the significance of obtaining an accurate PBA solubility model for solving the polymer expansion process.

3. Test Method

Assuming that the PBA solubility in Component B is the same as that in the mixture of Components A and B, in order to eliminate the interference of carbon dioxide caused by blowing reaction, the solubility test was carried out in Component B grout.

The changes of grout volume and temperature with respect to time were obtained by heating Component B grout. The volume change, as well as the volume of PBA, can be calculated as follows:

where ΔV(t) is the volume of the PBA at time t, V₀ is the initial volume of Component B, and V(t) is the measured volume of Component B at time t.

The Ideal Gas Law was put forward by Clapeyron in 1834, and the equation is

where p is atmospheric pressure, R is the universal gas constant with 8.31 J/(mol·k), and T is the temperature.

The mole number of the gaseous PBA can be calculated from equation (10) with T(t) and ΔV(t), which were obtained previously.

The mass of the gaseous PBA at time t is obtained as follows:

where n(t) is the mole number and MB is the molecular weight of PBA.

Since the total mass of the PBA remains constant, the mass of the liquid PBA dissolved in Component B at time t can be obtained with

Given the total mass of Component B grout with m, the solubility of the PBA at time t can be calculated as

where mB is the total mass of the PBA and mBL is the mass of the PBA in the liquid phase.

By repeating the above steps, the solubility of the PBA with corresponding temperature at each moment can be calculated, so as to obtain the solubility curve of the PBA with the change of temperature [23].

4. Test Process

In order to test the solubility curve of the PBA, the heating and testing device of PBA is designed, as shown in Figure 2. The testing device consists of graduated metal measuring cylinder, the temperature sensor, the induction cooker, and so on. Component B grout was heated by the electromagnetic heating. The grout temperature increased after heating. The PBA vaporized and formed tiny bubbles distributed in the grout, which caused the progressively grout expansion.

Figure 2: Heating and testing device.

The test of the solubility model of the PBA is mainly divided into the following steps:

(i)Inject Component B grout into a graduated metal measuring cylinder.(ii)Turn on the induction cooker; the metal measuring cylinder generates heat energy with the electromagnetic effect, which causes the expansion of Component B grout.(iii)Use the temperature sensor to read and record the temperature value of Component B grout during the heating process.(iv)Read the height of the grout level of Component B on the measuring cylinder at the same time; the volume value of Component B grout can be calculated according to the height value.(v)Calculate the solubility of PBA corresponding to the temperature according to the recorded temperature and volume value; thus, the solubility model is obtained.

The conditions of the used polymer grout are as follows.

At the initial moment, the mass of the PBA in Component B grout is 7.2 g, the molecular weight of the PBA is 116.9 g/mol, the mass of Component B grout is 372 g, the initial temperature is 299.65 k, the volume of Component B grout is 150 cm³, and the standard atmospheric pressure is 1.01325 × 10⁵ Pa.

Figure 3 shows the volume and temperature variations of Component B grout at different times. It can be seen that, with the continuous temperature decrease, Component B grout expands constantly. This indicates that, during the heating process, the vaporized PBA inside Component B leads to the bubble formation and the grout expansion.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 3: The volume and temperature changes of Component B grout at different times.

The grout is estimated to be composed of a cylindrical base and half of an oblate ellipsoid of revolution. Based on this characteristic, Ireka et al. [24] proposed a calculation method of grout volume, as shown in the following equation:

where R is the radius, h(t) is the height from the bottom to the highest point of the liquid level, and h₁(t) is the height of cylindrical base.

Figure 4 shows the schematic diagram of how to read the grout height.

Figure 4: The schematic diagram of grout volume reading.

The temperature change of Component B grout is displayed in Figure 5; it is obvious that the grout temperature increases as time goes on.

Figure 5: The temperature change curve with respect to time of the Component B grout.

Figure 6 shows the variation of volume with respect to time of Component B grout. The volume of grout increases slowly from 0 s to 150 s. After 150 s, the volume increases rapidly and gradually slows down.

Figure 6: The volume change curve with respect to time of the Component B grout.

The curve of Component B grout volume change with respect to the change of temperature was obtained according to the method of part two, which is plotted in Figure 7. It shows that the grout volume expands slowly at the beginning and tends to be fast when the temperature reaches about 360 k.

Figure 7: The volume change curve of Component B grout with respect to the temperature change.

Figure 8 is the calculation curve of the PBA with respect to the temperature change. The points and the curves represent the real experimental value and the fitting value, respectively. At the initial moment, the solubility of the PBA decreases slowly with the increase of temperature. When the temperature is greater than about 360 k, the solubility starts to decrease rapidly.

Figure 8: The solubility curve of the PBA with respect to the temperature change.

The fitting curve, which is shown in Figure 7, was fitted with the experimental data. The curve can be expressed by the following formula:

where a = −2.67781 × 10⁻⁵, b = 0.02725, c = −7.26908 × 10⁻⁶, d = 0.00478, and e = 0.75947.

5. Applicability Verification

The following steps were applied to verify the solubility model of the PBA:

(i)Carry out the foaming experiment of self-expanding polymer grout to obtain its temperature and volume changes with respect to time.(ii)Obtain the curve of the polymer density change with respect to time according to the volume change.(iii)Use the solubility curve of PBA to simulate the chemical reaction process of polymer. The calculated value of polymer temperature and density changes are gained according to simulation.(iv)Compare the calculated temperature and density variations with the corresponding experimental results, respectively. The measured solubility curve of the PBA can be judged to be applicable if the deviations are small.

The two-component polymer foaming experiment was carried out in a graduated measuring cylinder. Before starting the experiment, the measuring cylinder was placed on an electronic scale. The temperature sensor was placed inside the measuring cylinder. After a certain amount of polymer grout was injected into the measuring cylinder, the chemical reaction of the polymer grout occurred. The liquid height and the polymer temperature were recorded at different times. The grout volume was calculated according to the liquid height with equation (13). The corresponding density value could be obtained with the grout volume.

The foaming process was recorded in real time, as shown in Figure 9. It can be seen that the polymer begins to expand rapidly after a few seconds of injection. The grout temperature increases quickly as well. After about 60 s, both the rates of grout expansion and the temperature increase slow down and eventually become stable. The polymer mass remains almost constant during the reaction process, which is obvious in Figure 8. It indicates that the vaporized PBA is still located in the mixture and does not overflow.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 9: The temperature and volume changes of polymer grout at different times.

The reaction process of two-component polymer grout was then simulated numerically. Relevant parameters required for numerical simulation are listed in Table 1. The numerical code was written by using the programming language, Fortran.

Table 1: Relevant parameters required for numerical simulation.

Parameter
Value
Parameter
Value

cOH,0 (mol·m⁻³)
4079.6
CP (J/kg·k)
1800

cNCO,0 (mol·m⁻³)
4557.7
CCO2 (J/kg·k)
836.6

cW,0 (mol·m⁻³)
294.7
CBL (J/kg·k)
1159

T₀ (k)
300.75
CBG (J/kg·k)
1000

AOH (m³·s⁻¹·mol⁻¹)
2.7029
CW (J/kg·k)
4200

AW (m³·s⁻¹·mol⁻¹)
26.54
ρP (kg/m³)
1100

EOH (J·mol⁻¹·k⁻¹)
35195.707
ρBL (kg/m³)
1228

EW (J·mol⁻¹·k⁻¹)
42046.74
MB (g/mol)
116.90

(−ΔH)OH (J·mol⁻¹)
77050
MCO2 (g/mol)
44

(−ΔH)W (J·mol⁻¹)
86000
rBL,0
0.0196

The reaction temperature of the self-expanding polymer grout was calculated by considering the formulas of two chemical reaction kinetics equations with equations (2) and (4) and the energy balance equation with equation (5). The Runge-Kutta method was adopted to calculate the equation set. The formula of polymer density is a function of temperature, which is shown in equation (6). Therefore, the polymer density at a certain time could be obtained after the polymer reaction temperature was figured out.

Figure 10 shows the experimental and simulated results of polymer grout temperature during the reaction process. The polymer grout temperature increases sharply at the beginning of the reaction and then slows down. The simulated temperature is higher than that of experiment from about 25 s, but it gradually reaches agreement with the experimental result later. At around 45 s, the difference value reaches the maximum point; however, the error between the two is still within 2%. The simulated value is in good agreement with the experimental value, which proves the accuracy of the solubility model.

Figure 10: Comparison of the polymer reaction temperature.

The comparison between the experimental value and the simulated value of polymer density is shown in Figure 11. Before about 20 s, the polymer density drops quickly, and it tends to be stable after that. The experimental and simulated curves are in a similar trend with respect to time. However, the calculated density is a little bit greater than the experimental value from 15 s to 25 s. The biggest error between the two is 5% at around 25 s. In general, the numerical simulation result is close to the experimental result, which further illustrates the applicability of the solubility model.

Figure 11: Comparison of the polymer density.

6. Summary and Conclusions

An accurate solubility model of PBA is the prerequisite for studying the diffusion mechanism of self-expanding polymer grout. The main work and conclusion are summarized as follows:

(1)This paper provides a fast and effective way to accurately obtain the PBA solubility model and its verification process of polymer, which is more complete than the previous researches.(2)The experiment of two-component polymer grout was carried out to test the temperature and density changes during the grout expansion process. The temperature and density values obtained by the simulation are compared with the experimental values. The error between the simulated and experimental values is within 5%, which verifies the applicability of the solubility model.(3)As can be seen from the comparison figures (Figures 10 and 11), some biases still exist. The main reason for this phenomenon comes from the error of heating experiment of Component B grout.(4)The solubility of different PBAs in the grout composed with different raw materials is different. In this paper, only one of the commonly used polymers has been studied. Further experiments of variety polymers should be conducted to improve this research.(5)The research lays a theoretical foundation for both the simulation analysis of diffusion polymer process and the grouting mechanism study.

Nomenclature

Ai: | Preexponential factor
Ei: | Activation energy
ci,0: | Initial concentration
ci: | Current concentration
Xi: | Fractional conversion
R: | Universal gas constant
T: | Grout temperature
ri: | Mass fraction
ρi: | Density
Ci: | Heat capacity
(−ΔH)i: | Heat of reaction
λ: | Heat of evaporation
Mi: | Molecular weight
p: | Atmospheric pressure
Vi: | Volume
n: | Mole number
mi: | Mass
hi: | Height

Indices

OH: | Hydroxyl group
W: | Water
CO₂: | Carbon dioxide
NCO: | Cyanate group
BL: | Liquid phase of physical blowing agent
BG: | Gas phase of physical blowing agent
B: | Physical blowing agent
P: | Polyurethane
F: | Polymer grout foam.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by the National Key Research and Development Plan (Grant nos. 2017YFC0504503 and 2018YFB1600200), the Program for Science and Technology Innovation Talents in Universities of Henan Province (Grant no. 19HASTIT041), the National Natural Science Foundation of China (Grant nos. 51908514, 51878624, 51878622, and U1404525), the Natural Science Foundation of Henan Province (Grant no. 202300410746), Henan Science Fund for Distinguished Young Scholars (Grant no. 202300410354), Postdoctoral Research Sponsorship in Henan Province (Grant no. 19030026), and Major Scientific and Technological Special Project in Henan (Grant no. 181100310400).

References

D. Klempner, Handbook of Polymeric Foams and Foam Technology, Hanser Publishers, Munich, Germany, 2nd edition, 2004.

G. Gustafson and H. Stille, “Stop criteria for cement grouting,” Felsbau, vol. 23, no. 3, pp. 62–68, 2005.

L. R. Shucai, “Protection against water or mud inrush in tunnels by grouting: a review,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 8, no. 5, pp. 753–766, 2016.

R. Pusch and Knutsson, “Study of cement-grout penetration into fractures under static and oscillatory conditions,” Tunnelling & Underground Space Technology, vol. 45, pp. 10–19, 2015.

W. H. Lee, S. W. Lee, T. J. Kang, K. Chung, and J. R. Youn, “Processing of polyurethane/polystyrene hybrid foam and numerical simulation,” Fibers and Polymers, vol. 3, no. 4, pp. 159–168, 2002.

H. Z. Saboorian-Jooybari, Z. M. Dejam, and Z. Chen, “Heavy oil polymer flooding from laboratory core floods to pilot tests and field applications: half-century studies,” Journal of Petroleum Science and Engineering, vol. 142, pp. 85–100, 2016.

B. Zhang, L. Wang, and M. Hao, “An analytical solution for the radial flow of variable density grout in rock fractures,” Construction and Building Materials, vol. 206, no. 10, pp. 630–640, 2019.

J. M. Castro, “Studies of mold filling and curing in the reaction injection molding process,” AICHE Journal, vol. 28, no. 2, pp. 250–260, 1982.

L. Wang, M. Hao, and Y. Zhong, “A quasi-3D numerical model for grout injection in a parallel fracture based on finite volume method,” Complexity, vol. 2019, no. 3, Article ID 4139616, 18 pages, 2019.

C. Qian, S. J. Mumby, and B. E. Eichinger, “Phase diagrams of binary polymer solutions and blends,” Macromolecules, vol. 34, no. 12, pp. 2540–2545, 1991.

R. Tesser, M. Di Serio, A. Sclafani et al., “Modeling of polyurethane foam formation,” Journal of Applied Polymer Science, vol. 92, no. 3, pp. 1875–1886, 2004.

E. Santacesaria, M. D. Serio, and R. Gargiulo, “Vapour-liquid equilibrium measurements for binary mixtures of R32, R143a, R134a and R125 with a perfluoropolyether lubricant,” Journal of Fluorine Chemistry, vol. 121, no. 1, pp. 15–22, 2003.

E. Santacesaria, A. Passaretti, and M. Di Serio, “Phase equilibria in binary mixtures refrigerant + fluorinated lubricating oil: vapor−Liquid and Liquid−Liquid measurements,” Journal of Chemical & Engineering Data, vol. 49, no. 4, pp. 838–846, 2004.

E. Santacesaria, E. Musso, and M. D. Serio, “Description of the vapor-liquid equilibrium in binary refrigerant/lubricating oil systems by means of an extended Flory-Huggins model,” Journal of Fluorine Chemistry, vol. 99, no. 1, pp. 29–36, 1999.

V. K. Gupta and D. V. Khakhar, “Formation of integral skin polyurethane foams,” Polymer Engineering & Science, vol. 39, no. 1, pp. 164–176, 1999.

S. M. Geier, M. C. Winkler, and M. Piesche, “Numerical simulation of mold filling processes with polyurethane foams,” Chemical Engineering & Technology, vol. 32, no. 9, pp. 1438–1447, 2009.

A. J. Rojas, J. H. Marciano, and R. J. Williams, “Rigid polyurethane foams: a model of the foaming process,” Polymer Engineering and Science, vol. 22, no. 13, pp. 840–844, 1982.

D. J. R. Seo and J. R. Youn, “Numerical analysis on reaction injection molding of polyurethane foam by using a finite volume method,” Polymer, vol. 46, no. 17, pp. 6482–6493, 2005.

S. A. Baser and D. V. Khakhar, “Modeling of the dynamics of R-11 blown polyurethane foam formation,” Polymer Engineering & Science, vol. 34, no. 8, pp. 632–641, 1994.

S. A. Baser and D. V. Khakhar, “Modeling of the dynamics of water and R-11 blown polyurethane foam formation,” Polymer Engineering and Science, vol. 34, no. 8, pp. 642–649, 1994.

J. Kosek, M. Karimi, and D. L. Marchisio, “Multi-scale modelling of expanding polyurethane foams: coupling macro- and bubble-scales,” Chemical Engineering Science, vol. 148, pp. 55–64, 2016.

S. D. Lipshitz and C. W. Macosko, “Kinetics and energetics of a fast polyurethane cure,” Journal of Applied Polymer Science, vol. 21, no. 8, pp. 2029–2039, 1977.

M. D. L. Karimi and D. L. Marchisio, “A baseline model for the simulation of polyurethane foams via the population balance equation,” Macromolecular Theory and Simulations, vol. 24, no. 4, pp. 291–300, 2015.

I. E. Ireka, D. Niedziela, J. S. K. TröltzschSchäfer, F. Helbig, T. Chinyoka, and L. Kroll, “Computational modelling of the complex dynamics of chemically blown polyurethane foam,” Physics of Fluids, vol. 27, no. 11, pp. 113102–113192, 2015.