- 1Center for Durability & Sustainability Studies of Shandong Province, Qingdao University of Technology, Qingdao, China
- 2State Key Laboratory of Hydroscience and Engineering, Tsinghua University, Beijing, China
The old mortar attached to recycled aggregate (RA) is the main reason for the difference in water movement between RA concrete (RAC) and natural aggregate concrete. In this study, considering the old and new interfacial transition zones, a five-phase composite model for describing the water transport and distribution in RAC is established at the mesoscale. The key parameters describing water unsaturated transport in two types of mortar, saturated hydraulic conductivity (Ks) and van Genuchten model parameters (α, n), are obtained through the constant-head permeability test and isothermal adsorption test. By using the finite element method, the numerical simulations of unsaturated moisture movement in the homogeneous mortar, natural aggregate concrete, and five-phase RAC are systematically carried out. The proposed water transport model in the matrix is validated by comparison with the available experimental findings from the literature. The results show that the model can well predict unsaturated water transport in cement-based materials, including RAC. A parameter sensitivity analysis is undertaken to ascertain the main influencing factors of water transport in RAC. It is concluded that the RA replacement rate (Rra), the thickness of the old mortar (dm), and the aggregate volume fraction (Fa) are the primary parameters affecting moisture movement in RAC.
Introduction
With the rapid progress of infrastructure construction, two contradictions have gradually emerged. One is the imbalance between excessive and unsustainable exploitation of natural resources and the ever-growing requirements of building materials. Another is that the increasing wastes from construction and demolition conflicts with inadequate waste disposal methods. This phenomenon poses a severe threat to the economic development and ecological equilibrium in many countries and regions all over the world (Oikonomou, 2005; Xiao et al., 2012; Fraj and Idir, 2017; Tam et al., 2018; Lam et al., 2019; Song et al., 2020). Thus, processing the construction waste into recycled aggregates (RAs) instead of the natural aggregates (NAs) is known as an effective way to resolve these two contradictions. However, compared with NA concrete (NAC), RA concrete (RAC) whose applications are limited to low-level civil construction works are rarely used in modern infrastructure construction due to the inferior performance of RA caused by its complicated composition.
As is known to all, there are at least two interface transition zones (ITZs) in RAC (Tam et al., 2005; Xiao et al., 2013): 1) the old ITZ between NA and its adhered old cement mortar, which already exist in RA; and 2) the new ITZ between the RA and new cement mortar. Compared to traditional concrete, this caused more defects in mechanical properties and more significant hidden dangers of durability problems (Bao et al., 2020a; Wang et al., 2020). Thus far, investigations mainly focused on the mechanical properties of RAC (Folino and Xargay, 2014; Kim and Yun, 2014; Velay-Lizancos et al., 2018; Wang and Xiao, 2018; Munir et al., 2020), yet still limited in durability problems.
Most of the durability degradation mechanisms of cement-based materials can almost be attributed to the distribution and migration of moisture which either causes physical damage to the concrete structure such as freeze-thaw cycles (Zhang et al., 2017b; Wang et al., 2019; Bao et al., 2020b) (Zhang et al., 2017b; Bao et al., 2020b) or carries deleterious agents, such as chloride and sulfate, into the interior of cement-based materials (Phillipson et al., 2007; Tian et al., 2020; Zhang et al., 2020). Thus far, there are few experimental methods to study water transport directly. Neutron radiography is an excellent non-destructive testing method (Zhang et al., 2017a). However, inhomogeneous moisture distribution in the thickness direction caused by the complex structure of RAs will result in the moisture distribution overlap during the imaging process and eventually lead to abnormal images. As for other advanced non-destructive testing techniques such as X-ray computed tomography and nuclear magnetic resonance, there are also certain limitations. X-ray computed tomography (Smyl et al., 2016; Xue et al., 2020) has higher resolution than neutron imaging and can realize spatial detection of structures relatively faster, yet limited by the distinction between matrix and moisture. Nuclear magnetic resonance (Rucker-Gramm and Beddoe, 2010; Wyrzykowski et al., 2017; McDonald et al., 2020) is sensitive to hydrogen nucleus, however, it is limited by its parameters setting (such as the determination of surface relaxation rate) and imaging resolution. Thus, numerical methods may be an effective method for moisture distribution research.
In recent years, some contributions have been made to improve the mass transfer theory of cement-based materials by using numerical simulation methods. For example, considering the coupling effect of moisture and heat, Burkan Isgor and Razaqpur (2004) developed a macroscopic spatial and temporal simulation of the carbonation front advancement in cracked concrete. Based on Voronoi tessellation, Wang et al. (2016) presented a 2D numerical lattice model of mass transport in concrete on the mesoscale level to study the migration of moisture and chloride ions in cracked concrete under the action of dry and wet cycles. According to the theory of computational mechanics, Caggiano et al. (2018) conducted mesoscopic numerical simulation research on water transport of concrete under capillary action by combining zero-thickness interface and pipe elements. On the basis of the unsaturated flow theory, Li et al. (2017, 2018) carried out a numerical simulation study on the moisture distribution characteristic of ITZ and the crack network. Zhao K. et al. (2020) carried out study on chloride transport in cement mortar during drying process which evaluates the influence of diffusion and convection on chloride transport by experiments and numerical simulation. Zhao H. et al. (2020) proposed a coupled hygro-thermo-chemical model, and accurately predicted the three stage of coupled humidity and temperature filed via numerical and experimental method. However, there are still few studies on numerical simulation of the mass transfer process in RAC, especially lacking the water unsaturated transport research.
For the purpose of revealing the water migration in a heterogeneous RAC of multiple phases, this work is aimed to find the water transport characteristic in RAC and calculate the distribution of water content at any elapsed time instead of developing a complicated model for water transport in real RAC. The model in this paper focused on the mesoscale, whereby RAC was idealized as a five-phase composite porous media made up of two kinds of ITZs, two types of mortars, and NCs. Nonlinear diffusion equations based on unsaturated flow theory, were used to describe the water movement, and solved by the finite element method (FEM). And the model parameters were determined by the experimental method in this paper.
Unsaturated Flow Theory
In principle, the moisture migration in cement-based materials is dominated by two mechanisms: permeability and absorption. If the cement-based materials are thoroughly saturated with water, Darcy’s law can be used to describe the process of moisture movement, i.e., permeability. Nevertheless, saturated flow is the exception, and most building constructions can not serve in a saturated state so that capillary absorption should be considered as the primary mechanism of water transport. According to the work of Swartzendruber (1969) in the field of soil, the extended Darcy’s law can be used as the governing equation for describing water unsaturated transport,
where u is the flow velocity (LT−1); K(Θ) is the unsaturated hydraulic conductivity which is also known as unsaturated permeability (LT−1);
where
In fluid mechanics, the law of conservation of mass is expressed as a Continuity Equation, which is as follows:
Combining Eq. 1 with Eq. 3 leads to the fundamental equation of unsaturated flow, namely the Richards’ equation
Usually, Richards’ equation can also be expressed in another two equivalent forms:
where
A serious limitation of the Θ-based Richards Equation is that the materials of interest must be homogeneous. However, RAC is a kind of multiple-phase composite material which can not be considered as a homogeneous material. In this way, across the interfaces between different phases, the distribution of water content becomes discontinuous. In other words, the Ψ-based form, which can balance the saturation and the capillary potential well, is more suitable for this study.
In Eq. 1, the unsaturated hydraulic conductivity K(Θ) can hardly be obtained experimentally for cementitious materials. According to soil mechanics theory, it can be well estimated by the most widely used model, Mualem equation (Mualem, 1976),
where
with α (L−1) and n as fitting constants and properties of the material. These parameters can adequately describe the water retention curve of cement-based materials, so as to express the unsaturated water transport characteristics of different phase under capillary suction.
Determination and Validation of Parameters
In order to solve the governing equation, the first issue is to determine the parameters (i.e., Ks, α, and n) of different phases. Since the data for cementitious materials are sparse, to ensure the accuracy and effectiveness of the simulation results, adsorption tests, and static pressure penetration tests were conducted in this work. The process is as shown in Figure 1.
Materials and Specimen Preparation
Ordinary Portland cement 52.5 type and river sand with a maximum grain size of 5 mm were used to prepare the mortar specimen. Cube specimens with the side length of 70.7 mm and cylinder specimens with the size of Φ100 mm × 50 mm were prepared separately. All these specimens were demolded after 24 h and then cured in water at an ambient temperature of 20 ± 2°C. After 7 days, specimens were cut into slices of 5 mm thin and then soaked in water for 90 days to minimize the negative impact of hydration on subsequent tests. The mix proportion was determined, according to the research of Zhang et al. (2017a), as shown in Table 1.
TABLE 1. Mix proportion of cement mortar, kg/m3 (The permission has been obtained from the copyright holders).
Isothermal Adsorption Test
The adsorption isotherm is the basis for studying the properties related to water transport of materials. In order to obtain the van Genuchten parameters α and n, it is necessary to achieve the time-dependent adsorption curves. The following experiments were therefore carried out.
In the beginning, different humidity environments need to be established. Five inorganic salts were selected to prepare supersaturated salt solutions namely MgCl2, NaBr, NaCl, KCl, and K2SO4, by whose theoretical humidity that they can regulate is 33, 58, 76, 86, and 98%, respectively. Correspondingly, five humidity chambers with different salt solutions were prepared in a constant temperature and humidity room [T = 20 ± 2°C and relative humidity (RH) = 50 ± 2%]. The RH (–) of each chamber was measured repeatedly by using humidity sensors until it reached a steady state (the fluctuation in 24 h did not exceed ±1.5%). The actual RH were shown in Table 2.
Secondly, the cured square sheet specimens of M1 and M2 were dried to a constant weight at 105 C (mass variation in 7 days is less than 0.01 g), and the dried mass was recorded as
the relationship between
The relationship between hydrostatic pressure and RH can be described by the Kelvin-Laplace equation,
where Pw denotes the hydrostatic pressure (L−1 MT−2),
The value of α and n were obtained by fitting Eq. 9, as shown in Figure 2A.
The isotherm adsorption curves fitted in Figure 2A can be transformed into water retention curves in Figure 2B by Eq. 6. The saturation can be easily converted into capillary potential according to these curves, and the initial condition value and boundary condition value of the matrix can be determined quantitatively. It can be manifested from the profiles that when the capillary potential reaches Ψ = −2.16 × 10−6 mm, the saturation of M1 and M2 is almost equal to zero (that is, the dry state).
The Static Pressure Penetration Test
In Richards’ equation, the value of saturated hydraulic conductivity Ks is also an important parameter to define the properties of the matrix material, which can be obtained by the static pressure penetration test. The pressure penetration apparatus shown in Figure 3A is manufactured according to the research of Ludirdja and Young (1989).
FIGURE 3. Static pressure test: (A) penetration apparatus, and (B) curves of cumulative volume of water permeated vs. time for M1 and M2.
Firstly, each mortar specimen was wrapped with the rubber ring and sealed at the junction by epoxy resin to prevent water leakage. And then, the upper chamber, the rubber ring, and the lower chamber were connected in sequence. The lower chamber was filled with water until the liquid level in the stainless-steel pipe reached the top. In the same way, the upper chamber and the pipette were filled with water carefully to 305 mm. The top of the stainless-steel pipe was kept horizontal with the upper surface of the specimen so that the water at the top will penetrate through the specimen into the lower chamber under the single action of water pressure and flow out from the stainless-steel pipe. The volume change of the water in the pipette was recorded every 24 h as qi (L3), and then water was filled again to 305 mm height. Due to the little change of water level within 24 h, the specimens can be considered to be under the load of a constant head at all times. The above procedure was repeated until qi reached a constant value. The relationship between
where L denotes the thickness of the specimen (L); A denotes the permeability area (L2), which is the cut surface area of each specimen; h denotes the water head (L), which is the distance from the liquid surface at the top of the pipette to each specimen.
Validation of the Models
It is difficult to measure the distribution of moisture in non-transparent cement-based materials by experimental methods, and the process of non-destructive testing is also limited. Therefore, in order to verify the accuracy of the theoretical model and parameters, the results of the neutron radiography test carried out by Zhang et al. (2017a) were chosen to compare with the simulation results in this paper. The specimens used in the analysis are mortar specimens with a water-cement ratio of 0.6, whose mix proportion is the same as that of M1 in Table 1. The prepared mortar specimen was cut into thin slices of 100 mm × 100 mm × 25 mm and sealed with aluminum foil tape around to prevent moisture intrusion from other boundaries instead of the bottom surface (Figure 4A). The sealed specimen was placed in a flume for the capillary water absorption experiment, and the neutron radiography test was administered in the meantime.
FIGURE 4. (A) Schematic and (B) model of water absorption test, and the comparison of experimental and numerical results: (C) the neutron images of water distribution by Zhang et al. (2017a) (The permission has been obtained from the copyright holders), (D) water distribution cloud maps by numerical simulation, and (E) one-dimensional water distribution profiles by Zhang et al. (2017a) (The permission has been obtained from the copyright holders).
According to the capillary absorption test, the cement mortar absorption model can be established, as shown in Figure 4B. The matrix is initially in a dry state and can only contact the water through the bottom boundary. The remaining boundaries do not exchange water with the outside environment. Therefore, the initial condition and boundary condition can be expressed by
where
The cloud maps of water distribution at different times can be obtained by employing numerical simulation, and the results are compared with the experimental results of neutron imaging. By comparing Figures 4C,D, it can be observed intuitively that the moisture in the specimen is absorbed upward from the lower surface of the specimen under the action of the capillary adsorption force. Due to the inhomogeneity of the mortar composite, the water front is not completely parallel to the water absorption surface, and the water distribution at different coordinates is also slightly different. However, these differences did not substantially affect the results. The distribution of water front illustrates that the water transport is almost one-dimensional parallel to the bottom of the boundary, and the movement of the water front slows down gradually with time elapsing. The simulated position of the water front at different times shows good agreement with the experimental results.
The marked central area shown in Figures 4C,D is selected for quantitative analysis to obtain numerical and experimental water distribution shown in Figure 4E. It can be seen quantitatively that the water distribution along the y-axis in the simulation results is roughly the same as that in the test results. The difference between the simulation results and the test results mainly comes from the following points: a) The mortar model is simplified to a homogeneous model structure whereby water is uniformly distributed in the horizontal direction. However, in the capillary absorption test, due to various factors such as uneven distribution of sand particles and insufficient hydration, the results will inevitably be different from the idealized model’s. Sand and insufficient hydration will cause the lag of unsaturated water transport. b) In the capillary water absorption test, the sample needs to be dried by oven. High temperature will damage the material, enhance the connectivity of the pores, and accelerate the process of water transport in the matrix. c) The inhomogeneity of the mortar specimens causes the saturation water content at each position to be different, so it is inevitable to underestimate or overestimate the saturation of certain areas after normalization by the formula. However, the influence of these errors is very weak and within an acceptable range, thus overall results show well agreement. This agreement suggests that the numerical model and substitute parameters correctly describe the capillary absorption process in the cement mortar matrix, which provides a guarantee for the accuracy of subsequent studies on water unsaturated transport in RAC.
Transport Analysis for Recycled Aggregate Concrete
Mesoscale Recycled Aggregate Concrete Model
Considering that the polygonal aggregate has a huge amount of calculation for the simulation and the limited influence on water transport (Li et al., 2016; Liu et al., 2018), the round aggregates are used in our study to simplify the calculation. As mentioned before, the essential difference between RAC and NAC is the coarse aggregate. Compared to NA, RAC can not be considered as traditional aggregate-matrix two-phase material structure nor aggregate-ITZ-matrix three-phase material structure because of the old mortar attached to old NA. The old ITZ initially existed between the old aggregate and the old mortar, and new ITZ gradually forms between the old mortar and the new mortar after setting. As a consequence, in order to describe the migration of water in RAC more accurately, the five-phase RAC model consisting of old aggregate, old mortar paste, old ITZ, new mortar paste, and new ITZ were established by simplifying the actual RA into a concentric circle model (Figure 5A).
FIGURE 5. Mesoscale structure of recycled aggregate concrete: (A) schematic diagrams and (B) finite element meshes diagrams.
By means of the Monte Carlo method, in other words, cyclically putting concentric circles of different sizes into a square area until it reaches the expected area fraction (i.e., coarse aggregate volume fraction), the two-dimensional RAC random aggregate model can be established. This process was achieved by randomly generating NC aggregate center coordinates, radius, and adhesive mortar thickness and running a set of MATLAB codes. The meshed RAC model is shown in Figure 5B.
Definition of Interface Transition Zones Phase
According to the description of the smooth plate crack flow method (Chen et al., 2013; Li et al., 2017), a variant of extended Darcy’s law was adopted to describe the water flow in ITZ. In this way, ITZs in the model will have no geometric width, and water transport in it might be simplified to flow along interior boundaries between the two different phases. The governing equation of water transport in the crack is as follows,
where
where
where l is a constant equal to −0.5 for cementitious materials (–);
where μ is the viscous coefficient of water. With these formulations above, the ITZ domain can be idealized as a porous and zero-model-thickness interface with higher hydraulic conductivity in RAC.
Definition of Natural Aggregate and Mortar Phases
The pore structure determines the mass transfer performance of cementitious material. Since the pore structure of RC is different from NC, the water transfer law in RAC is naturally different from NAC. The migration velocity of water in the matrix can be regarded as a function of the characteristics of the RAC’s different phases. In other words, correctly distinguishing and defining the properties of each phase to represent its pore structure characteristic in numerical simulation is an essential factor in ascertaining the water distribution in RAC. In this model, both the cement mortar phase and the aggregate phase are considered as the porous material so that the water transport in them follows the Richards’ Equation (i.e., Eq. 3). Consequently, combining the previously mentioned theory of water transport in the matrix, the governing equation of water transport in recycled concrete can be obtained, which is expressed as follows,
And the initial and boundary condition can be expressed as
In the Eq. 16, if the governing equation in domain
Compared with the new cement mortar matrix, due to various factors such as crushing and grinding, the micro-cracks of the old adhesive mortar and the old ITZ increase, and their saturated water content and permeability also increase accordingly. This means that the saturated water content
Validation of the Recycled Aggregate Concrete Models
To verify the accuracy of the model, we compared the simulation results with the capillary water absorption test results of Bao et al. (2020a) which is shown in Figures 6A,B. It is illustrated that, both the macroscopic water absorption and the microscopic water penetration depth, the simulation results agree well with experimental results. The differences are mainly due to the following factors: a) There is a slight difference in water cement ratio between experiment (W/C = 0.39) and simulation (W/C = 0.40) which might cause differences in results; b) There are some uncertainties in practice, such as the insufficient hydration and the heterogeneity of mortar matrix, which is difficult to achieve by numerical model; c) Although a complex RA model has been established, it is still a simplified form of practical RA. These factors may cause small changes in the results but will not have a substantial impact on the conclusions. Hence one can conclude that the RAC model used in this study is accurate and reliable.
FIGURE 6. Comparison of numerical results and the experimental results of Bao et al. (2020a) on (A) water penetration depth and (B) amount of water absorption (The permission has been obtained from the copyright holders).
Modeling Results of Water Uptake in Recycled Aggregate Concrete and Natural Aggregate Concrete
Equation 16 was implemented in the commercial finite element software, COMSOL Multiphysics, the multi-field coupled FEM simulation program, and the cloud map of water distribution of NAC and RAC at different times were obtained to describe the 2D water distribution in RAC and NAC at different times (Figures 7A,B).
FIGURE 7. Comparison of numerical simulation results between recycled aggregate concrete (RAC) and natural aggregate concrete (NAC): cloud maps of water distribution in (A) RAC and (B) NAC, the arrow diagram of flow flux in (C) RAC and (D) NAC, and the schematic of water movement characteristic between aggregates in (E) RAC and (F) NAC.
It can be seen that, as t = 30 min, the water front in RAC has wholly surpassed the marked aggregate, and the water front behind the aggregate penetrates deeper than it in the mortar matrix. While the moisture penetrating speed in NAC slows down, and the water front in the mortar penetrates is deeper. When t = 60 min, the difference in the penetration depth is more prominent. The peak of the water front has reached a depth of about 13 mm in RAC, while barely about 7 mm in NAC.
Despite considering the influence of ITZ of NAC, the hindrance effect of NA on water transport is much higher than the promoting effect of ITZ. By comparing the flow flux arrow diagram shown in Figures 7C,D, one can see in NAC that although the flow velocity of ITZ is higher than that of the ordinary mortar matrix, it is still lower than the moisture transmission in the narrow and long domain between NAs. However, in RAC, the flow velocity between RAs is relatively low, and the main flow velocity is concentrated in the old adhesive mortar and ITZ area. This phenomenon can be explained that water only moves in the matrix and single ITZ on condition that NA is assumed as impermeable, whereby the real water transport area is equivalent to narrowing. Thus, due to the limitation of water transport capacity, the vast majority of water still crowds in the mortar matrix, so that the convex surface of the water front is formed between the aggregates (Figure 7E). However, the old mortar and the double ITZs are all weak areas for RAC. Both of the porosity and hydraulic conductivity are relatively large. On the one hand, similar to NAC, the water absorption of RAC is accelerated due to the existence of ITZ. On the other hand, as moving in ITZ, water also transfers rapidly into the old mortar whose capillary pressure is higher than that of the new mortar, which further accelerates the penetration of water. Because of the higher penetration velocity in RA, the convex surface of the moisture front is formed behind RA (Figure 7F).
The relationship between the water penetration depth and the square root of the time were plotted in Figure 8, which is approximated to a positive linear correlation. The slope of the curve corresponds to the velocity of water uptake, which indicated that the water absorption velocity in RAC is much larger than the other two materials. As t = 600 min, the penetration depth of water in RAC even reaches 1.6 times that in NAC. By comparing NAC and cement mortar, it can be illustrated that the penetration depth in RAC always slightly lags behind in cement mortar, which confirms that the hindrance effect of NA on moisture is stronger than the promoting effect of single ITZ as mentioned before. However, the error bars of moisture penetration in NAC shows lower variability, which indicates that NA and single ITZ have limited effects on water uptake. Instead, this situation is quite different in RAC. “Convenient access” of RA occupies a dominant position in water transport, which leads to the increased influence of matrix heterogeneity, and then causes larger variability of penetration depth showed in profile. This feature can also be found intuitively in the cloud map in Figures 7A,B, as t = 480 min.
FIGURE 8. The profiles of water penetration depth in recycled aggregate concrete (RAC), natural aggregate concrete (NAC), and cement mortar at different time.
Parameter Sensitivity
The previous content illustrates the critical role RA plays in unsaturated absorption of RAC. Therefore, it is necessary to conduct a more detailed study on the impact of aggregate. For this purpose, parameter sensitivity analyses on RA replacement rate (
Effect of Recycled Aggregate Replacement Rate on Water Absorption
As shown in Figure 9A, it can be learned that with the
FIGURE 9. (A) 2D cloud diagram and (B) profiles of moisture distribution with different Rra (Fa = 50%).
Effect of Thickness of Old Mortar on Water Absorption
To investigate the influence of old adhesive mortar on water movement in RAC, without changing the size of the NC, the RAC model is generated. Figure 10A indicates that with the increase of
FIGURE 10. (A) 2D cloud diagram and (B) profiles of moisture distribution with different dm (Fa = 50%).
Effect of Recycled Aggregate Volume Fraction on Water Absorption
The cloud map in Figure 11A shows that the change in Fa causes the difference in moisture distribution. However, this change seems to be reflected only in the difference of water front tortuosity instead of the depth. The water penetration depth of different
Conclusions
In this paper, the influence of aggregates and ITZs on the water transport performance of concrete are studied by combining the FEM and unsaturated theory. The numerical simulation approach is reported to evaluate the unsaturated transport process in RAC. The following conclusion can be drawn:
• On the basis of the non-homogeneous characteristic of concrete, the five-phase RAC model consisting of old aggregate, old adhesive mortar, old ITZ, new mortar, and new ITZ is established on the mesoscopic level. The capillary absorption equation in unsaturated RAC is reasonably derived, and the model parameters required in the equation are determined by the experimental method. The accuracy of the model and parameters are validated by comparing experimental data.
• The double ITZs and old mortar in RAC provide a more convenient channel for moisture movement, which results in moisture preferentially transporting along the aggregate instead of the new mortar matrix and causing the forefront of the wet front always forming behind the RA. This transmission mode greatly accelerates the capillary absorption rate of water in the concrete, thereby seriously threatening the durability of this concrete.
• The parameter sensitivity analysis of recycled concrete is conducted. It is found that as the thickness of the old adhesive mortar and the replacement rate of RA increase, the water absorption velocity in RAC is significantly accelerated. The moisture uptake velocity increases with the growth of aggregate volume fraction, but the amplitude is relatively small, which shows that the hindrance effect of NA on moisture is offset by the promotion effect of adhesive mortar with dual ITZs.
Data Availability Statement
All datasets presented in this study are included in the article.
Author Contributions
The corresponding author is responsible for ensuring that the descriptions are accurate. ZL: writing original draft, implementation of the numerical simulation and data curation. PZ: Formulation of overarching research goals and aims. JB: development of methodology. YH: investigation, modifying and editing.
Funding
Financial support for on-going projects by the Natural Science Foundation of China (51922052, U1706222, 51778309), Natural Science Foundation of Shandong Province (ZR2018JL018, ZR2019PEE001) and Open Research Fund Program of State Key Laboratory of Hydroscience and Engineering (SKLHSE-2019-C-04) is gratefully acknowledged.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
Bao, J., Li, S., Zhang, P., Ding, X., Xue, S., Cui, Y., et al. (2020a). Influence of the incorporation of recycled coarse aggregate on water absorption and chloride penetration into concrete. Constr. Build. Mater. 239, 117845. doi:10.1016/j.conbuildmat.2019.117845
Bao, J., Xue, S., Zhang, P., Dai, Z., and Cui, Y. (2020b). Coupled effects of sustained compressive loading and freeze‐thaw cycles on water penetration into concrete. Struct. Concr. 2020, 1–11. doi:10.1002/suco.201900200
Burkan Isgor, O., and Razaqpur, A. G. (2004). Finite element modeling of coupled heat transfer, moisture transport and carbonation processes in concrete structures. Cem. Concr. Compos. 26 (1), 57–73. doi:10.1016/s0958-9465(02)00125-7
Caggiano, A., Said Schicchi, D., Mankel, C., Ukrainczyk, N., and Koenders, E. A. B. (2018). A mesoscale approach for modeling capillary water absorption and transport phenomena in cementitious materials. Comput. Struct. 200, 1–10. doi:10.1016/j.compstruc.2018.01.013
Chen, B., Song, E., and Cheng, X. (2013). “Plane-symmetrical simulation of flow and heat transport in fractured geological media: a discrete fracture model with comsol,” in Multiphysical testing of soils and shales (Berlin, Germany: Springer), 149–154.
Folino, P., and Xargay, H. (2014). Recycled aggregate concrete – mechanical behavior under uniaxial and triaxial compression. Constr. Build. Mater. 56, 21–31. doi:10.1016/j.conbuildmat.2014.01.073
Fraj, A. B., and Idir, R. (2017). Concrete based on recycled aggregates – recycling and environmental analysis: a case study of Paris’ region. Constr. Build. Mater. 157, 952–964. doi:10.1016/j.conbuildmat.2017.09.059
Kim, S.-W., and Yun, H.-D. (2014). Evaluation of the bond behavior of steel reinforcing bars in recycled fine aggregate concrete. Cem. Concr. Compos. 46, 8–18. doi:10.1016/j.cemconcomp.2013.10.013
Lam, P. T. I., Yu, A. T. W., Wu, Z., and Poon, C. S. (2019). Methodology for upstream estimation of construction waste for new building projects. J. Clean. Prod. 230, 1003–1012. doi:10.1016/j.jclepro.2019.04.183
Li, X., Chen, S., Xu, Q., and Xu, Y. (2017). Modeling the three-dimensional unsaturated water transport in concrete at the mesoscale. Comput. Struct. 190, 61–74. doi:10.1016/j.compstruc.2017.05.005
Li, X., Chen, S., Xu, Q., and Xu, Y. (2018). Modeling capillary water absorption in concrete with discrete crack network. J. Mater. Civ. Eng. 30 (1), 04017263. doi:10.1061/(asce)mt.1943-5533.0002122
Li, X., Xu, Y., and Chen, S. (2016). Computational homogenization of effective permeability in three-phase mesoscale concrete. Constr. Build. Mater. 121, 100–111. doi:10.1016/j.conbuildmat.2016.05.141
Liu, Q.-f., Feng, G.-l., Xia, J., Yang, J., and Li, L.-y. (2018). Ionic transport features in concrete composites containing various shaped aggregates: a numerical study. Compos. Struct. 183, 371–380. doi:10.1016/j.compstruct.2017.03.088
Lockington, D., Parlange, J.-Y., and Dux, P. (1999). Sorptivity and the estimation of water penetration into unsaturated concrete. Mater. Struct. 32 (5), 342. doi:10.1007/BF02479625
Ludirdja, D., Berger, R. L., and Young, J. F. (1989). Simple method for measuring water permeability of concrete. ACI Mater. J. 86 (5), 433–439. doi:10.14359/2000
McDonald, P. J., Istok, O., Janota, M., Gajewicz-Jaromin, A. M., and Faux, D. A. (2020). Sorption, anomalous water transport and dynamic porosity in cement paste: a spatially localised 1H NMR relaxation study and a proposed mechanism. Cem. Concr. Res. 133, 106045. doi:10.1016/j.cemconres.2020.106045
Mualem, Y. (1976). A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resour. Res. 12 (3), 513–522. doi:10.1029/WR012i003p00513
Munir, M. J., Kazmi, S. M. S., Wu, Y.-F., Patnaikuni, I., Zhou, Y., and Xing, F. (2020). Stress strain performance of steel spiral confined recycled aggregate concrete. Cem. Concr. Compos. 108, 103535. doi:10.1016/j.cemconcomp.2020.103535
Oikonomou, N. D. (2005). Recycled concrete aggregates. Cem. Concr. Compos. 27 (2), 315–318. doi:10.1016/j.cemconcomp.2004.02.020
Phillipson, M. C., Baker, P. H., Davies, M., Ye, Z., McNaughtan, A., Galbraith, G. H., et al. (2007). Moisture measurement in building materials: an overview of current methods and new approaches. Build. Serv. Eng. Res. Technol. 28 (4), 303–316. doi:10.1177/0143624407084184
Rucker-Gramm, P., and Beddoe, R. E. (2010). Effect of moisture content of concrete on water uptake. Cem. Concr. Res. 40 (1), 102–108. doi:10.1016/j.cemconres.2009.09.001
Smyl, D., Hallaji, M., Seppänen, A., and Pour-Ghaz, M. (2016). Quantitative electrical imaging of three-dimensional moisture flow in cement-based materials. Int. J. Heat Mass Tran. 103, 1348–1358. doi:10.1016/j.ijheatmasstransfer.2016.08.039
Song, Q., Zhao, H., Jia, J., Yang, L., Lv, W., Bao, J., et al. (2020). Pyrolysis of municipal solid waste with iron-based additives: a study on the kinetic, product distribution and catalytic mechanisms. J. Clean. Prod. 258, 120682. doi:10.1016/j.jclepro.2020.120682
Swartzendruber, D. (1969). “The flow of water in unsaturated soils,” in Flow through porous media. Editor R. J. M. De Weist (New York, NY: Academic Press), 215–292.
Tam, V. W. Y., Gao, X. F., and Tam, C. M. (2005). Microstructural analysis of recycled aggregate concrete produced from two-stage mixing approach. Cem. Concr. Res. 35 (6), 1195–1203. doi:10.1016/j.cemconres.2004.10.025
Tam, V. W. Y., Soomro, M., and Evangelista, A. C. J. (2018). A review of recycled aggregate in concrete applications (2000–2017). Constr. Build. Mater. 172, 272–292. doi:10.1016/j.conbuildmat.2018.03.240
Tian, Y., Zhang, P., Zhao, K., Du, Z., and Zhao, T. (2020). Application of Ag/AgCl sensor for chloride monitoring of mortar under dry-wet cycles. Sensors 20 (5), 1394. doi:10.3390/s20051394
van Genuchten, M. T. (1980). A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 44 (5), 892–898. doi:10.2136/sssaj1980.03615995004400050002x
Velay-Lizancos, M., Vazquez-Burgo, P., Restrepo, D., and Martinez-Lage, I. (2018). Effect of fine and coarse recycled concrete aggregate on the mechanical behavior of precast reinforced beams: comparison of FE simulations, theoretical, and experimental results on real scale beams. Constr. Build. Mater. 191, 1109–1119. doi:10.1016/j.conbuildmat.2018.10.075
Wang, C., and Xiao, J. (2018). Evaluation of the stress-strain behavior of confined recycled aggregate concrete under monotonic dynamic loadings. Cem. Concr. Compos. 87, 149–163. doi:10.1016/j.cemconcomp.2017.12.012
Wang, L., Bao, J., and Ueda, T. (2016). Prediction of mass transport in cracked-unsaturated concrete by mesoscale lattice model. Ocean Eng. 127, 144–157. doi:10.1016/j.oceaneng.2016.09.044
Wang, R., Yu, N., and Li, Y. (2020). Methods for improving the microstructure of recycled concrete aggregate: a review. Constr. Build. Mater. 242, 118164. doi:10.1016/j.conbuildmat.2020.118164
Wang, Y., Cao, Y., Zhang, P., Ma, Y., Zhao, T., Wang, H., et al. (2019). Water absorption and chloride diffusivity of concrete under the coupling effect of uniaxial compressive load and freeze-thaw cycles. Constr. Build. Mater. 209, 566–576. doi:10.1016/j.conbuildmat.2019.03.091
Witherspoon, P. A., Wang, J. S. Y., Iwai, K., and Gale, J. E. (1980). Validity of Cubic Law for fluid flow in a deformable rock fracture. Water Resour. Res. 16 (6), 1016–1024. doi:10.1029/WR016i006p01016
Wyrzykowski, M., McDonald, P. J., Scrivener, K. L., and Lura, P. (2017). Water redistribution within the microstructure of cementitious materials due to temperature changes studied with 1H NMR. J. Phys. Chem. C 121 (50), 27950–27962. doi:10.1021/acs.jpcc.7b08141
Xiao, J., Li, W., Corr, D. J., and Shah, S. P. (2013). Effects of interfacial transition zones on the stress–strain behavior of modeled recycled aggregate concrete. Cem. Concr. Res. 52, 82–99. doi:10.1016/j.cemconres.2013.05.004
Xiao, J., Li, W., Fan, Y., and Huang, X. (2012). An overview of study on recycled aggregate concrete in China (1996–2011). Constr. Build. Mater. 31, 364–383. doi:10.1016/j.conbuildmat.2011.12.074
Xue, S., Zhang, P., Bao, J., He, L., Hu, Y., and Yang, S. (2020). Comparison of Mercury Intrusion Porosimetry and multi-scale X-ray CT on characterizing the microstructure of heat-treated cement mortar. Mater. Charact. 160, 110085. doi:10.1016/j.matchar.2019.110085
Zhang, P., Dai, Y., Wang, W., Yang, J., Mo, L., Guo, W., et al. (2020). Effects of magnesia expansive agents on the self-healing performance of microcracks in strain-hardening cement-based composites (SHCC). Mater. Today Commun. 25, 101421. doi:10.1016/j.mtcomm.2020.101421
Zhang, P., Hou, D., Liu, Q., Liu, Z., and Yu, J. (2017a). Water and chloride ions migration in porous cementitious materials: an experimental and molecular dynamics investigation. Cem. Concr. Res. 102, 161–174. doi:10.1016/j.cemconres.2017.09.010
Zhang, P., Wittmann, F. H., Vogel, M., Müller, H. S., and Zhao, T. (2017b). Influence of freeze-thaw cycles on capillary absorption and chloride penetration into concrete. Cem. Concr. Res. 100, 60–67. doi:10.1016/j.cemconres.2017.05.018
Zhao, H., Jiang, K., Yang, R., Tang, Y., and Liu, J. (2020). Experimental and theoretical analysis on coupled effect of hydration, temperature and humidity in early-age cement-based materials. Int. J. Heat Mass Transf. 146, 118784. doi:10.1016/j.ijheatmasstransfer.2019.118784
Keywords: mesoscale model, recycled aggregate concrete, Richards’ equation, interfacial transition zone, water unsaturated transport
Citation: Liu Z, Zhang P, Bao J and Hu Y (2020) Numerical Simulation of Water Transport in Unsaturated Recycled Aggregate Concrete. Front. Mater. 7:560621. doi: 10.3389/fmats.2020.560621
Received: 09 May 2020; Accepted: 24 August 2020;
Published: 21 September 2020.
Edited by:
Hongyan Ma, Missouri University of Science and Technology, United StatesReviewed by:
Haitao Zhao, Hohai University, ChinaJinbo Yang, Shandong Agricultural University, China
Qiang Zeng, Zhejiang University, China
Copyright © 2020 Liu, Zhang, Bao and Hu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Peng Zhang, zhp0221@163.com