Thermal Simulation of Rolled Concrete Dams: Inﬂuence of the Hydration Model and the Environmental Actions on the Thermal Field

: Mathematical models for the simulation of the thermal evolution of roller-compacted concrete (RCC) dams during construction constitute an important tool for preventing excessive temperature rise, which may lead to cracking and losses of functionality. Here, we present a framework for the simulation of the thermal process. We deﬁne the boundary conditions of the problem using a careful description that incorporates the main heat exchange mechanisms. We adopt both a non-adiabatic and an adiabatic heat generation model for the simulation of the cement hydration. Our numerical framework lets us study the effect of the adopted heat generation model on the thermal ﬁeld. Moreover, we study the inﬂuence of the weather conditions on the evolution of the hydration, and on the starting date of construction. Our simulations have shown that the hydration model has an important inﬂuence over the temperature ﬁeld during the construction and the heat generation rate. Moreover, the hydration process and the temperature evolution are driven by the weather conditions. Once the next lift is cast, its thermal insulation effect makes the hydration take place under quasi-adiabatic conditions. As expected, dams built in cold months are prone to dissipate more heat than those built in warm seasons.


Introduction
Roller-compacted concrete (RCC) dams are especially prone to cracking while relatively new. The RCC technology is economically profitable thanks to the high production rates, which are achieved by the use of high capacity placing and compaction equipment for construction. Moreover, longitudinal joints and pipe cooling techniques are avoided, and transverse joints are kept to a minimum [1]. The high construction velocity, together with the massive nature of the structure, makes the heat released during concrete hydration a key aspect in the design of these structures [2].
The high construction rates of RCC dams may lead to significant temperature rises. Such an increase occurs over the first days after setting when the stiffness of concrete is low and viscous effects are important. As a result, the deformations of concrete associated with the temperature rise produce small thermal stresses. Several months after the beginning of the hydration, stiffness is significantly increased. However, the temperature of the concrete is still relatively high in the core of the dam because the massive nature of the structure does not stimulate the dissipation of the hydration heat. The inner temperature may increase up to the adiabatic temperature rise of the concrete in those points where the reaction occurs under adiabatic conditions. Instead, the hydration heat of the The NAH model proposed by Cervera et al. [20] and Cervera and Garcia-Soriano [21] has been used to simulate the thermo-mechanical behavior of the Rialb RCC gravity dam (Spain) during the construction and the first years of operation [22], and to simulate the Urugua-í RCC gravity dam located in Argentina [2,20]. This model has also been used to conduct thermal probabilistic simulations of RCC dams under construction [23]. Azenha [24] adopted the formulation of Reinhardt et al. [18] to compute the thermo-mechanical behavior at the early stages of RCC dams' lifetimes [24]; that of slabs [25] and the foundations of wind turbines [26]; and that of laboratory-based cubic samples of concrete [27]. The NAH approach incorporates the concrete temperature as an input variable, among others. The weather conditions are then considered when hardening occurs under no-adiabatic conditions. The heat generation rate is not given by an unique evolution, resulting in a widely variety of curves. The concrete temperature rise may considerably differ from the evolution provided by the AH approach, since the exposure of concrete to weather conditions varies during the construction: the last lift built is exposed to weather until the cast of the next lift.
Boundary conditions play a vital role in the simulations. Since the dissipation of the hydration heat depends on their magnitude, reliable methodologies for computing them are of vital importance. Concrete temperature varies due to the heat released during the hydration process and the exchange between the structure and the environment. The exchange is caused by the convection phenomenon, the long wave radiation exchange, the solar energy, and the evaporative cooling. Furthermore, the night radiative cooling should be considered when computing the long wave radiation exchange. Santillan et al. [28] proposed a methodology for computing the heat fluxes in concrete dams, taking into account all the previous phenomena. The methodology was validated by the comparison between the outputs of the models and the recorded data in an arch dam using both numerical [28,29] and analytical approaches [30].
Two types of boundary conditions were imposed in previous studies. In some studies, the temperature of the faces of the dam was set equal to the ambient temperature [12]. In most cases, convection between the faces of the dam and the environment was established with a constant value of the convection coefficient. The ambient temperature was taken as the environmental one [2,11,13,20,31,32], and in some cases, it was increased by some degrees to take into account the solar radiation [6,14,22]. Only a few studies considered all the reported heat exchange mechanisms. Salete and Lancha [33] incorporated all of them, except the night radiative cooling, yet Azenha [24] accounted for all them. Nevertheless, the effects of these heat exchange mechanisms on the hydration reaction have not been assessed so far.
Here, we study the temperature evolution of concrete during the construction of RCC gravity dams. We conduct numerical simulations of the construction process with high temporal and spatial resolution. We adopt both a non-adiabatic hydration formulation and an adiabatic model for concrete hydration. Moreover, our methodology for computing the boundary conditions considers the convection phenomenon, the long wave radiation exchange, the solar energy, and the evaporative cooling. The simulations allow us to characterize and quantify the influence of the hydration model on the temperature fields. Furthermore, the effect of the environmental actions on the hydration process is also assessed. The research queries are, then, how the hydration model and the environmental actions can alter the evolution of the temperature field, and whether such a detailed methodology is necessary or not.
The paper is organized as follows. Section 2 provides the governing equations of our simulation framework. First, we introduce the concrete hydration models. Afterward, we present the methodology for computing the heat fluxes and the initial temperature of the foundation. In the last part of this section, we describe the numerical implementation of the governing equations. The case study where we apply the methodology is described in Section 3. The results of the numerical simulations and their analysis are reported in Section 4. Finally, some conclusions are drawn in Section 5.

Methodology
Heat transport through a solid domain Ω ⊂ δ with spatial dimensions δ ∈ {2, 3} is governed by the heat diffusion equation or Fourier equation as follows [5]: where ρ is the density of the solid, c is the specific heat, θ is the temperature, t is time, k is the thermal conductivity, and Q are the sources or sinks of heat. The boundary conditions simulate the interactions of the solid with the surrounding environment. We impose heat fluxes at dam faces, and temperature at the bottom boundary of the foundation.
In this section, we describe the main ingredients of our numerical framework. We begin with the hydration process, responsible for the term ∂Q/∂t =Q = q in the Fourier equation. We continue with the description of the methodology to compute the heat fluxes between the dam faces and the surrounding environment. We also introduce the approach to compute the temperature at the bottom part of the foundation. Finally, we present the numerical technique used to solve the Fourier equation.

Concrete Hydration Models
AH models describe the hydration reaction as a function of the age of concrete. Among the available formulations to simulate the AH temperature rise, the exponential model is often used to simulate the thermal evolution of RCC dams under construction. The adiabatic temperature rise of the concrete θ ad in K reads [7,8] where θ ad,M is the maximum temperature rise under adiabatic conditions in K, β is a parameter that accounts for the heat generation rate in day −1 , and τ is the age of the concrete in days. The cumulative heat generated Q in J/m 3 due to the hydration reaction is provided by and the heat generation rate ∂Q/∂t = q in W/m 3 reads AH models consider that concrete hydrates under adiabatic conditions, and consequently, the reaction is only driven by time. The velocity of the reaction is controlled by the parameter β, whose value depends on the initial temperature of concrete. In contrast, NAH models simulate the hydration reaction as a function of the temperature of the concrete and the degree of hydration α.
The formulation proposed by Reinhardt et al. [18] states that the heat generation rate may be described by the product of two independent functions. The first function f (α), named normalized heat generation rate, reads whereQ M is the maximum heat generation rate in W/m 3 . f (α) accounts for the degree of heat development. The second function g(θ) reproduces the influence of the temperature of the concrete on the hydration process. It provides the maximum heat generation rate at temperature θ. g(θ) is formulated with an Arrhenius-type law as follows where A T is a rate constant in W/m 3 , E a is the apparent activation energy in J/mol, and R is the ideal gas constant in J/(molK). The heat generation rate is then provided by The degree of hydration, α, is the ratio between the amount of cement that has reacted at time t and the initial amount of cement [34]. Since it is not straightforward to measure α with experiments, it is approximated by the degree of heat development r [35]. r is the ratio between the heat released until a time t and the total heat released in the hydration reaction Q ∞ : This formulation of α may be differentiated in regard to time to rendeṙ Therefore, combining Equations (7) and (9), it may be stated thaṫ

Heat Fluxes
The dam faces in contact with air exchange heat with the surrounding environment by four mechanisms: convection, long wave radiation, solar radiation, and evaporation of water. The total heat flux q T across a unit area and time in W/m 2 reads where q c is the convective heat transfer, q r is the long wave radiation exchange, q s is the heat flux generated by the solar radiation, and q ev is the evaporative cooling. The convective heat transfer between the surface of the concrete and the air in contact with it is modeled by the Newton cooling law: where H c is the heat transfer coefficient in W/(m 2 K), θ a is the ambient temperature in K, and θ is the temperature of the concrete's surface. We determine H c with the Kehlbeck equations [36], which read if the convection is forced, and H c = 3.50 (14) if natural. In the above equations V w is the wind velocity in m/s. The long wave radiation exchange, q r , between the surface of the concrete and the air in contact is governed by Stefan-Boltzmann law. A simplified expression for computing q r is [37]: θ sk being the sky temperature in K, and H r the radiation exchange coefficient in W/(m 2 K).
where e is the emissivity of the dam surface. The sky temperature can be approximated to the ambient temperature during the day [24]. However, at night θ sk is lower than the ambient one. This phenomenon is called night cooling. We compute θ sk with [38] θ sk = 0.736 + 0.00577θ dp where θ dp is the dew point temperature, which is assessed by the Clausius-Clapeyron equation as function of the relative humidity and the ambient temperature [39]. The heat fluxes due to the solar radiation and the evaporative cooling are computed following the methodology of [28]. The solar radiation generates a heat flux given by where a is the solar absorptivity of the concrete surface (dimensionless) and I T is the incident solar radiation on the surface in W/m 2 . I T is determined from the hourly diffuse and beam components of the solar radiation on a horizontal surface through the Reindl model [40]. The beam component is not considered when the surface is shaded. The evaporative cooling is heat loss due to the evaporation of the water stored on the surface of the concrete. q ev is computed through where q m is the moisture evaporative flux in kg/(m 2 s) and h w is the latent heat of evaporative water in J/kg. q m is given by E s being the moisture emission coefficient in kg/(m 2 sPa), e s the saturation vapor pressure in Pa, and h r the relative humidity. The Lewis' relationship provides E s as function of the convective heat transfer coefficient, the specific heat capacity of air c s (1007 J/(kg K)), and the total air pressure P in Pa. It reads [41] H c E s ≈ 29c s P 18 .
We can use an equivalent air temperature, θ eq , which is the temperature that with only the convection mechanism provides the same heat flux as the four described mechanisms [33]. Introducing Equations (12) and (15) into Equation (11) and expressing q T by the Newton law of cooling using θ eq , the following expression is achieved.
The equivalent ambient temperature reads and the equivalent convection coefficient H eq is then given by

Foundation Temperature
Kusuda and Achenbach [42] analyzed underground temperature sets from 63 recording stations spread across the United Stated. They found that the earth temperature can be estimated through where θ g (z, t) is the temperature of the ground in K at depth z in m and time t in days, θ g,m is the annual average ground temperature at z = 0 m, A g,sm is the annual amplitude of the ground surface temperature at z = 0 m, χ g is the thermal diffusivity of the ground in m 2 /day, and t 0 is the time lag in days. We infer from Equation (25) that the amplitude of θ g exponentially decreases with depth. At z = 20 m, the amplitude is almost zero, which means that θ g is approximately equal to θ g,m . Moreover, θ g,m can be approximated by the annual average air temperature [42]. Equation (25) also allows us to approximate the initial temperature of the foundation in our simulations.

Numerical Implementation
The governing Equation (1) is a partial differential equation. We adopt a finite volume scheme on a fixed regular quadrilateral mesh for the spatial discretization and a backward Euler scheme for the time integration. The non-adiabatic hydration model constitutes a source of non-linearity in our problem. We handle the non-linearity with a staggered strategy for updating the involved fields.
We adopt the staggered scheme proposed by Miehe et al. [43], initially used for solving highly nonlinear fracture mechanics problems. The framework has been validated with experimental observations [44]. This strategy has also been successfully applied for solving highly-nonlinear, fluid-driven fracture propagation simulations in elastic [45] and poroelastic solids [46], and results have been compared with analytical solutions [47].
The staggered scheme in the time interval [t n , t n+1 ] consists of three steps. We assume that all fields are known at time t n .

1.
First, α is updated with the values of α n and θ n by using Equation (9).

2.
Then,Q is computed with the updated values of α and θ n by using Equation (7).
This strategy requires small time steps to properly capture the non-linearity of the problem.

Case Study
We apply our numerical framework to the construction of a theoretical case study inspired in Enciso dam, a RCC gravity dam. The dam was built on Cidacos river in Spain. The dam axis direction, from upstream to downstream, is north 66 • east. The latitude of the dam site is 42.15 • . The foundation is composed of limonites and sandstones, and a grout curtain provides a supplemental barrier to water flow through the dam foundation. The closest weather station to the dam site is Fitero, about 80 km far from the dam, and provides daily climatic recorded data. These are: air temperature (mean, maximum and time, minimum and time), relative humidity (mean, maximum and time, minimum and time), rainfall, mean wind speed, and global solar insolation. Data sets are available from 1 January 2014, to 31 December 2016. We interpolate data using the piece-wise cubic Hermite interpolating polynomial, which enjoys no overshoots and a shape-preserving aspect [48].
The height of the dam from the foundation to the crest is 98.4 m; the crest is 3.45 m in width; the length of the dam at the base is 81.21 m; and the height and width of the steps are 1.20 and 0.96 m respectively. Other dimensions of the dam are included in Figure 1a. We simulate the construction by placing lifts of 30 cm thick every 12 h, at 8 a.m. and 8 p.m. The construction is stopped when the ambient temperature is lower than 2 • C or higher than 36 • C. Concrete is cooled by replacing part of mixing water by ice. The placing temperature, θ p , is equal to the ambient temperature, θ a , when the latter is between 2 • C and 18 • C; otherwise, θ p = 18 • C when θ a is between 18 • C and 27 • C, and θ p = θ a − 9 when θ a is between 27 • C and 36 • C. The concrete curing lasts 14 days, and we simulate it by considering the horizontal surfaces to be humid.
Our domain compromises the dam and a portion of the foundation, Figure 1a. The domain is meshed with squared volumes of size 3 cm. The boundary conditions are: (1) convection at the surfaces in contact with the air at the corresponding equivalent ambient temperature and convection coefficient; (2) no heat flow at the vertical boundaries of the foundation; and (3) imposed temperature at the bottom of the foundation equal to the annual average air temperature. The initial temperature of the foundation is computed with Equation (25). We adopt a time step of 5 min. The properties of the foundation are: (ρc) f = 2120 kJ/(m 3 K) and k f = 2.97 W/(mK). The parameters of the hydration reaction of the concrete are: (ρc) c = 2300 kJ/(m 3 K), k c = 3.00 W/(mK), A = 7.22 × 10 10 W/m 3 , Q ∞ = 6.79 × 10 7 J/m 3 , and E a = 45.51 kJ/mol. The normalized heat generation function, f (α), against the degree of hydration, α, is represented in Figure 1b. These magnitudes were obtained from laboratory experiments conducted over samples of an RCC dam located in Portugal [24].
We simulated the adiabatic temperature rise with the previous parameters of the NAH for a sample whose initial temperatures were 16 • C and 25 • C, and we fit the AH exponential model to the simulated data. The best fits were found for β = 7.69 × 10 −6 s −1 and β = 1.21 × 10 −5 s −1 , respectively.

Results and Discussion
Our methodology is the starting point for the simulation of the thermal evolution in our case study during the construction. The simulations allow us to study the influence of the selection hydration model on the concrete temperature field and on the heat generation rate field. Moreover, the simulation of the hydration reaction with a NAH model lets us study the influence of the environmental actions and the starting date of the construction on the hydration reaction and the concrete temperature evolution.

Influence of the Hydration Model over the Temperature Field
We study the influence of the hydration model on the temperature field by simulating the construction of an RCC dam using three models. The construction starts on 1 April 2014, and lasts 165 days. In the first simulation we adopt a NAH model for the hydration of concrete (first column of Figure 2). In the second simulation we employ an AH model whose β-parameter is calibrated for an initial concrete temperature equal to the mean placing temperature, 16 • C (second column of Figure 2), and in the third we employ one an AH model with β-parameter calibrated for a higher temperature, 25 • C. We plot a vertical temperature profile of the dam computed with the three models in the third column of Figure 2. The temperature fields at mid-height construction time (2000 h) present significant differences between the three simulations. The simulation with the NAH model, Figure 2a, provides higher temperatures in the core of the dam than the ones with the AH model. The NAH model accounts for the effect of the weather conditions on the hydration process, the heat generation rate being higher as the concrete temperature increases. Because the placing temperature of the concrete is relatively low, the generation rate after setting is also low. Then, the hydration reaction slowly evolves for the first 12 h as concrete remains cold because most of the heat is dissipated into the atmosphere. After those 12 h, a new lift is placed and works as a thermal insulation lift; then the temperature of the concrete rises and the heat generation rate also does. As a result, the hydration reaction feeds back itself: as the concrete temperature is higher, the heat generation rate also does and makes the temperature increase.
The simulation with the AH model whose β-parameter is calibrated for the mean placing temperature, Figure 2b, provides lower temperatures in the core of the dam than the simulation with the NAH model, Figure 2c. The formulation of the AH model, Equation (4), has the maximum heat generation rate occur at the placing time, i.e., for τ = 0 s, and decrease exponentially with time. Moreover, during the 12 h after pouring concrete, the top surface of the lift is subjected to the environmental actions, and consequently, a significant part of the hydration heat is dissipated into the atmosphere. Furthermore, the relatively low ambient temperature in spring stimulates the dissipation of the heat into the atmosphere.
The simulation with the AH model with the highest β-parameter results in lower concrete temperatures in the core of the dam (Figure 2c). The β-parameter controls the heat generation rate, and q is higher as the parameter increases. Consequently, more heat is generated for the 12 h after placing, and more heat is dissipated into the environment, resulting in a lower concrete temperature rise.
Once construction has finished (time 4000 h), the core of the dam is hotter than the faces in all the simulations (Figure 2d,e). Nevertheless, the simulation with the NAH model results in higher temperatures than the simulations with the AH model in most parts of the dam; see Figure 2f. The difference accentuates for the AH model with the higher β-parameter.
The distribution of the heat generation rate at time 2000 h is highly influenced by the adopted hydration model. The NAH model- Figure 2g-captures the effect of the concrete temperature on the heat generation rate, whereas the AH model does not- Figure 2h-and provides a uniform rate within a given lift. The concrete temperature is also greatly influenced by the environmental actions. At time 2000 h, the last lift was cast at 8 a.m. on June 22, 2014, and the upstream vertical faces were exposed to sunbeams. Therefore, the concrete temperature rises, and the heat generation rate also does, as can be appreciated in the snapshot of Figure 2. The horizontal faces are also exposed to sunbeams, but the evaporation of the curing water cools the concrete and counteracts the effect of the solar radiation.

Influence of the Weather Conditions on the Evolution of the Hydration Reaction and Concrete Temperature
We illustrate the effects of the weather conditions on the evolution of the hydration reaction and the concrete temperature by analyzing the evolution of each variable at several points and under different environmental conditions; see Figure 3. We run two simulations with our case study: in the first one the construction starts on 1 April 2014, and takes place in spring; and in the second simulation it starts on 1 October 2014, and extend s into autumn. We plot the evolutions of several variables at four representative points. Two points are located at the top and the bottom of the lift  (Figure 3d). The NAH model provides a rich-variety of evolutions of q, in contrast to the AH model, which gives only one evolution. Moreover, q also drives the evolution of α, and in turn, the evolution of the mechanical properties.
The evolution of the heat generation rate q- Figure 3-is driven by the environmental conditions and the concrete temperature, θ, during the first 12 h. During that period, the top surface of the last lift is directly exposed to the weather conditions. Afterward, a new lift is placed and the concrete hydrates under quasi-adiabatic conditions. This situation is not given in those points closed to the dam faces, since concrete hydrates under quasi-isothermal conditions. Several evolution patterns are given, depending on the weather conditions-spring or autumn; casting time of concrete; day or night; and position within the lift, top or bottom.
The heat generation rate for the points of the dam built in spring, cast in the morning, and located in the top of the lift- Figure 3a-evolves following the ambient temperature, θ a . Since q is low during the first five hours, the concrete temperature follows θ a . Afterward, the concrete temperature exceeds the ambient one, as the heat released by the concrete is higher than the heat exchanged with the environment. Nevertheless, the drop in the ambient temperature makes the concrete cool down because the heat released into the environment is more than the generated heat. After 12 h, the next lift is cast and the concrete hydrates under quasi-adiabatic conditions. Although q decreases with time, the concrete temperature increases because the heat does not dissipate into the atmosphere. The evolution of q for the points located in the bottom of the lift is not so influenced by the environmental conditions. The concrete temperature gradually increases with time, since the thermal isolation effect of the mass over the point and the sound weather conditions do not allow the concrete to release all the generated heat.  The evolution of q for the points cast at night- Figure 3b-present a similar evolution to the previous ones. The point in the top follows the ambient temperature: initially θ decreases due to the low generated heat rate and the decrease in θ a , but afterward the concrete warms thanks to the increase in both q and θ a . This effect is more pronounced for the point in the bottom of the lift due to the thermal isolation of the concrete over it.
The heat generated in the dam built in autumn is also greatly influenced by the weather conditions. The concrete temperature of the points in the top of the lift evolves following the ambient tendency, it increases during the day- Figure 3c-and decreases at night-Figure 3d. Nevertheless, the temperature at the bottom always rises because of the thermal isolation of the concrete. In contrast to the dam built in spring, the heat generation rate is very low during the first 12 h, as the lift is directly exposed to the low ambient temperature that slows down the hydration reaction. Afterward, thanks to the thermal isolation of the new lift, the concrete warms and the reaction accelerates.

Effect of the Weather Conditions on the q-Peak and Time to Peak
The NAH model provides a rich-variety evolution of q, maximum heat generation rates q max , and timing of q max , t q max , in contrast to the AH model, which gives a unique evolution of q. The weather conditions have a direct impact on q max and t q max , as shown in Figure 4. The lifts cast on spring mornings- Figure 4a-reach higher q max than the lifts built in autumn. Moreover, most of the peaks in spring are given within the first 12 h after setting the concrete; i.e., before the next lift is built. In contrast, all the peaks in autumn occur after the first 12 h, once the thermal isolation effect of the next lift allows the concrete to warm up.  The lifts set at night present lower values of q max than those built in the mornings- Figure 4b. The night-time weather conditions make q max reach lower values. As in the previous case, the lifts built on spring mornings reach higher q max than the lifts built on autumn mornings. However, the timing of q max at the top of the lifts delays and peaks is reached after the next lift is cast. Regarding the bottom parts, t q max occurs both before and after casting the next lift in spring, whereas it is given after casting the next lift in autumn.
The maximum value of the heat generation rate within every lift for the dam built in spring and lifts placed in the mornings is higher at the bottom than the top- Figure 5a-due to the thermal isolation effect of the top concrete. In the same way, t q max is higher at the top than the bottom- Figure 5a-for those lifts where t q max at the top is lower than 12 h, but t q max is lower at the top than the bottom for those lifts where t q max at the top is higher than 12 h. This behavior comes from the weather conditions: cold days delay the timing of q max at the top and warm days bring it forward.
The lifts placed at night present a different behavior- Figure 5b. The lifts cast in autumn present in general higher values of q max and lower t q max at the bottom than at the top. In the case of the lifts disposed in spring, t q max is lower at the bottom than at the top, but no conclusion can be reached regarding q max .

Influence of the Boundary Conditions over the Concrete Temperature
We assess the influence of the boundary conditions over the concrete temperature evolution through a new set of simulations. We compute the boundary conditions with a simplified, straightforward methodology that is widely used. The heat exchange between the structure and the environment is boiled down to the convection between the surface of the dam and the air in contact with it. We compare the outputs of several thermal models that incorporate these boundary conditions with the outputs of other models whose boundary conditions are computed with our methodology.
We run two new simulations whose boundary conditions are assessed with the simplified approach: the first simulation corresponds to a dam built in spring, and the second one in autumn. We plot the evolution of the concrete temperature at four representative points. Two points are, The evolution of θ is highly influence by the boundary conditions, especially during the early ages of concrete. In the short term, less than 12 h after casting, both methodologies provide quite different results, whereas in the long term results are similar. Nevertheless, the simplified methodology always results in higher temperature rises. The differences arise from the heat exchange mechanisms: the simplified methodology does not account for the solar radiation, nor the long wave radiation exchange, nor the night cooling effect.
Our methodology provides higher concrete temperatures in the top points cast in the morning of May, Figure 6a, than the outputs provided by the model with the simplified approach. The higher temperature rise is explained by the effect of the solar radiation, which warms the concrete and speeds up the hydration reaction. This effect is not so noticeable in the top point cast in the morning of November, Figure 6c, since solar radiation in autumn is less intense than in spring. Nevertheless, after sunset our model provides lower concrete temperatures than the simplified methodology due to the night cooling effect, which cools down concrete and speeds down the hydration reaction.
The differences in temperatures at the bottom points of the lifts are also noticeable. Since our methodology incorporates the night cooling effect, concrete temperatures on the top surface of the dam are cooler than the outputs provided by the model with the simplified boundary conditions. Once the next lift is cast, the newly poured concrete releases heat to the bottom lift and the hydration reaction speeds down. Consequently, our methodology provides lower temperature rise for the bottom points.
The night cooling has also important effects on the lifts cast at night. The temperature of the top points computed with the model that incorporates our methodology drops at night, whereas the temperature rises in the models with the simplified methodology. This behavior is given in both the dam built in spring, Figure 6b, and the one built in autumn, Figure 6d.

Effect of the Construction Date on the Thermal Field
The weather conditions during the construction influence the hydration process, and consequently, the temperature field. For this purpose, we simulate the construction of two RCC dams whose constructions start in October, 2014-autumn, and April, 2014-spring,respectively. The construction lasts 165 days for each. We adopt the NAH formulation.
The temperature field for the dam built in spring, Figure 7a,b, is higher or equal at every point compared to the dam built in autumn, Figure 7c In spring, weather becomes warmer with time. As a result, less heat is dissipated, as the placing time is later and concrete temperature after hydration is higher. For this reason, the hottest part of the core is in the upper part of the dam built in spring. Instead, during the fall the weather becomes colder as time advances. The dam is more prone to dissipate heat as the placing time is later. The hottest areas of the core are in the lower part, which was built with warmer weather conditions.

Conclusions
We have simulated the thermal evolutions of roller-compacted concrete (RCC) dams under construction. Two important ingredients of our numeral framework are: (1) the heat generation formulation for simulating the hydration reaction of the concrete, and (2) the boundary conditions of the thermal problem.
We have adopted both a non-adiabatic hydration (NAH) model and an adiabatic hydration (AH) model for simulating the hardening of the concrete. The model considers the temperature of the concrete, θ, and the hydration degree for computing the heat generation rate, q. In contrast, the AH model assumes that concrete hydrates under adiabatic conditions. The evolution of q is only driven by time. Moreover, we have incorporated to the computation of the boundary conditions the following heat flux mechanisms: convection, long wave radiation, solar radiation, and evaporation of water.
Our simulations have shown that the NAH model provides higher temperatures in the concrete than the AH formulation. The AH model assumes that the maximum heat generation rate takes place just after placing the concrete. A large amount of heat is dissipated until the next lift is placed. Then, the new lift works as an thermal insulation and the temperature of the concrete continues rising. In contrast, the simulations with the AH model have evidenced that the hydration reaction is slow after placing the concrete. The concrete progressively warms and q increases. Then, the peak of q does not occur just after placing. Consequently, the amount of heat dissipated before the next lift is placed is, in general, lower than in the simulations with the NAH model.
The weather conditions also affect the hydration reaction. The AH formulation leads to a rich variety of evolutions for q. Both the weather conditions and the insulation provided by the newly placed lifts explain the evolutions. In contrast, the AH model results in a unique evolution for q.
We have adopted a detailed methodology to compute the boundary conditions of the thermal problem. Our methodology captures the main mechanisms that control the heat exchange of the dam with the surroundings. We have compared the outputs of a model with boundary conditions assessed with our methodology with the results of a model with simplified boundary conditions that only account for the convection with the air. In the short term, less than 12 h after casting, our methodology provides higher concrete temperatures in the top points of lifts cast in the mornings due to the solar radiation. Moreover, in the case of top points cast at night, our methodology results in lower temperatures due to the night cooling effect. In the long term, both approaches provide similar temperatures, but our methodology results in slightly lower temperatures compared to the outputs provided by the simplified approach.
We have also confirmed that the starting date of the construction affects both q and θ. Those dams built in autumn result in lower and delayed peaks for q than those built in spring. Moreover, the temperate field at mid-construction is greater or equivalent at every point for the dams built in spring compared to the ones built in autumn. The position of the hottest area in the core of the dam is also influenced by the weather conditions. The hottest area of those dams built in autumn is located in the lower part of the core, whereas for those built in spring it is in the upper part.
Our simulations have shown that the choice of the hydration model is a key aspect. Results of the thermal simulations may be underestimates if the model is not properly selected. Moreover, our investigation has also illustrated how the weather conditions influence the evolution of the temperature field. Our approach also captures the key controlling mechanisms of the heat exchange between the dam and the surrounding environment and shows promise for simulating and understanding more complex processes in the evolution of the temperatures during the construction of concrete dams.