Boolean-Based Surface Procedure for the External Heat Transfer Analysis of Dams during Construction

The external heat transfer of dams during construction is complex because such transfer is location specific and time varying. An external thermal model is developed in this paper. Five types of external heat flux are considered in the mathematical model: airside convection, electromagnetic radiation, absorbed solar input, water-side convection, and surface insulation effect. A method for extracting and classifying the external surfaces of dams on the basis of Boolean operations is proposed. Heat transfer conditions can be automatically set up for each step according to the proposed method, and the method can be used as a preprocessing facility for finite element analysis. A 285m high arch dam in Southwest China is examined as a study case. The model is implemented and found to correctly identify different types of external surfaces. Simulation result agrees well with the monitored temperatures.


Introduction
Concrete is the most widely used building material in the world.Rough statistics [1] show that over two billion ton quantities of concrete are produced all over the world per year.Temperature plays an important role in the construction of concrete structures [2].Considerable heat is released during the hydration reaction of young-age concrete.This heat results in temperature gradients near the structure surface because of external heat transfer under varying atmospheric conditions.A large volume of concrete increases the gradients and increases the risk of thermal-stress-induced cracks [3], thereby decreasing the long-term strength and durability of concrete.Thus, developing an appropriate thermal analysis method for concrete, particularly for mega structures (e.g., concrete dams), during the construction phase is essential and necessary.
A dam comprises an extremely large amount of concrete and is referred to as a typical case of massive concrete.To complete this colossal project, the structure must be divided into numerous blocks cast at different times.The entire construction usually requires a number of years to complete, which means that the climate condition changes with seasonal variation during the construction phase.A good example would be the construction of the famous Three Gorges Dam, which took a total of 14 years to complete 27.2 million cubic meters of concrete [4].
The boundary conditions should be consistent with the practical situation to ensure the accuracy of simulating the heat transfer problem.However, the setup work poses difficulty in simulating large and complex structures such as concrete dams.On one hand, the dam is bounded by a group of surfaces, which are subjected to various external thermal conditions.On the other hand, the dam constantly grows with the gradually changing boundary surfaces.Therefore, external thermal conditions are location specific and time varying, which makes the selection and setup of the boundary conditions complex.The surfaces of different steps are manually chosen in conventional thermal simulations, which is a tedious and time-consuming task that does not guarantee absolute accuracy.
A number of experimental and numerical studies on the external heat transfer problem of concrete were conducted.For instance, Demir [5] investigated the natural Mathematical Problems in Engineering convection temperature/velocity field around a horizontal concrete cylinder that is heated with a heat source.Lee et al. [6] conducted an experiment to investigate the convective heat transfer coefficient during the hydration process.G. F. Jones and R. W. Jones [7] studied a heat transfer problem for an insulated and reinforced concrete wall, in which the surface boundary is solved by using Newton's cooling law.However, these studies mainly focused on stable laboratory test conditions where the external heat transfer boundary condition is relatively simple.
Much effort was given to study the external heat transfer of dams.Léger and Seydou [8] and Léger and Leclerc [9] investigated the 1D/2D temperature field of some dams during the operating phase.Sheibany and Ghaemian [10] intensively studied the environmental effect of heat transfer in arch dam engineering but failed to consider the surface insulation effect.In another study, Jin et al. [11] investigated the solar effect on exposed dam surfaces during the operating phase.However, the objects of these studies were dams that had already been built; that is, the dams were already under the postconstruction or operation phase, where the situation of the external surface boundary regions barely changes.
The finite element (FE) method is widely used to simulate the thermal problem of dams during the construction phase.Numerous relevant studies address the FE method [2,[12][13][14][15][16][17][18][19][20][21].Cervera et al. [12] translated a 3D thermal problem into a 1D/2D FE analysis problem by simplifying the longitudinal direction.Wu and Luna [2] and Luna and Wu [15] studied the numerical implementation of the thermal stress problem in roller-compacted concrete (RCC) dams without considering solar and electromagnetic radiation.Saetta et al. [16] proposed the constituting factors of the external thermal flux.However, they did not discuss exact estimation methods.Noorzaei and Bayagoob [18,19] simulated RCC dams that are under construction to study wind and solar effects.However, the solar input in their study was simplified to a constant.
The numerical thermal simulation of dams had been studied, but the setup problem of the external boundary conditions at each step was not technically discussed.Limited studies focus on the preliminary procedure of external surfaces for the external heat transfer of dams.Therefore, this study presents a method for the automatic setup of external heat transfer boundaries during the analysis of the entire construction phase.The external heat transfer of dams is classified, and the basic mathematical model is established.Boolean-based algorithms are proposed to extract different types of external surfaces and to set appropriate heat transfer conditions.A super-high arch dam is examined as a study case by using the developed program to test the automatic method in practical engineering applications.

Mathematical Model
The governing equation of 3D unsteady heat transfer can be derived from the energy conservation principle and the Fourier law of heat conduction, which is written as follows: where  is the temperature;  is the thermal conductivity, which is always considered isotropic in concrete material; that is,  =   =   =   ;  is the internal heat source produced per volume, homogeneously distributed over the concrete. generally consists of two main parts: hydration heat generation and the cooling effect of embedded coolpipes, which was investigated in our previous work [14,17,20].
Massive concrete has complex boundary conditions during the construction phase.The external heat transfer of the boundary surfaces can generally be expressed as where  is the normal direction of the external surfaces and  is the external heat flux, which is a combination of five components.These components are discussed in detail in this section.Various combinations are applied in different surface areas to simulate a dam during construction.The rule of the combination setup is presented in detail in Section 4.

Air-Side Convection.
Heat convection occurs upon exposed concrete surfaces because of the bulk motion of air fluids that carry and disperse the heat [22].The airside convective flux   is proportional to the temperature difference between the concrete and the surrounding air.Heat convection can be given by using Newton's law of cooling or by using the Robin boundary condition: where   denotes the ambient temperature and ℎ  denotes the convective coefficient, which is considered constant in some previous studies on dams [2,11,14].However, the convective coefficient is actually affected by numerous factors, among which wind speed was highlighted in recent studies as the major factor [23,24].The thermal simulation of dams should consider this wind effect because dams are commonly built in canyons where wind is naturally strong.Nusselt-Jürges correlation [23], a well-established form of convective coefficient that considers wind speed, is used in this paper and can be defined as where   is the wind speed; , , and  are the parameters to be decided based on the applied condition;  []   is notably the degree Kelvin of ambient temperature.The following equation gives the conversion formula based on degree Celsius: Air temperature always varies and hardly remains constant even throughout a single day.In this work, the   of each step is used as the mean temperature over the increment of time.Constantly recording the ambient temperature is impractical; thus, we use the daily mean temperature to represent the entire day for convenience.We can employ a sinusoid function [25] curve with a period of 24 h for typical temperature history (see Figure 1).The mean daily temperature can be fit with two measured temperatures  1 and  2 at times  1 and  1 by using the following equation: where  Mean is the middle time between the peak and valley.When the record time is at the middle time such as  1 =  mean , the temperature measured at that time can be directly considered the mean temperature; that is,   =  1 .
Unfortunately, recording or accessing detailed local weather data in some engineering projects is difficult and cannot be guaranteed even though only one temperature measurement per day is taken.That implies that daily temperatures cannot be determined from (6).We can apply the following sinusoidal curve to approximate   over the entire year: where  denotes the day of the year (e.g., 1 for January 1st and 32 for February 1st).

Exposed Radiation to Environment.
Electromagnetic energy is constantly emitted by concrete through exposed surfaces and can be calculated by using the Stefan-Boltzmann law: where   is the emissivity of the surface and  = 5.670373 × 10 −8 W m −2 K −4 , which is referred to as the Stefan-Boltzmann constant [26].The temperatures applied in (8) should notably be used in degree Kelvin.

Absorbed
Radiation from the Sun.Solar radiation was not considered in numerous thermal studies on dams [2,14,15].A recent work [11] showed that irradiation has a great thermal effect on the exposed surfaces of a dam.Proper calculation of the absorbed solar radiation of the structure surface is complex because of considerations such as the sun's position, environmental factors, and surface conditions.The revolution orbit and rotation circle of the earth are not on the same plane because the equatorial plane tilts at an angle of 23.45 ∘ to the orbital plane.Solar declination , which is defined as the angle between the earth-sun line and the equatorial plane, consequently varies from 0 ∘ to 23.45 ∘ during one period of revolution (see Figure 2).The daily declination is obtained through approximation throughout a year as where  is the day of year.The equation shows that if 81 <  < 263.5 then  > 0, which is the time when the earth-sun line crosses the Northern Hemisphere, namely, the summer time for the Northern Hemisphere.Figures 2 and 3 show the position map of the sun and the structure surface of the earth.The relation among these angles through a 3D geometry deduction can be written as sin  = cos  cos  cos  + sin  sin , cos  = cos  cos  sin  − sin  cos  cos  , cos  = cos  cos  sin  + sin  cos , where  is the hour angle (the detailed calculation of this variable is tallied in the appendix);  is the latitude value (in degrees: positive for the northern site and negative for the southern site);  is the inclination angle of the structure surface;  is the angle of incidence, which is the angle between the earth-sun line and the line normal to the surface;  is the surface azimuth, which is defined as the angular displacement from the south direction to the normal direction of the structure surface on the horizontal plane;  is the solar altitude angle, which represents the angle between the horizontal plane and the emanating line from the sun;  is the surface-solar azimuth angle;  is the solar azimuth angle.
We can obtain the exact values of , , , , and  for a specific time on a specific location on earth.Then,  can be calculated by (10),  can be found by (11) and (12), and  is determined by (14).So far, all the solar angles involved in the irradiation computation are obtained.When  = 0,  notably reaches the daily maximum angle of /2 − | − |, which is the moment called the solar noon.
When the sun's position and solar angles are determined as described, we can calculate the total solar energy input that strikes the structure surface.This energy input can be measured by using three components: where  Beam ,  Diff , and  Refl represent the beam (direct), diffuse, and ground reflector components of solar radiation from the sky dome and  Beam and  Diff represent the irradiation coefficient, considering the surface angle.These components can be obtained through mathematical and geometry relations with the extraterrestrial solar radiation from a clear-sky model proposed by ASHRAE [27]: where  0 is the extraterrestrial solar radiant flux,  is the relative air mass,  is the parameter that represents optical depths, and  is the ground reflectance of the foreground.These four parameters can be investigated as described in the following paragraphs.
The extraterrestrial solar radiant flux  0 varies throughout the year owing to the slightly elliptical orbit path, which can be given by where  sc = 1367 W/m 2 and is referred to as the solar constant [28].
The relative air mass  represents the ratio of the actual air mass to the air mass with the sun shining directly overhead. can be considered a function of the solar altitude angle  [29] as follows: The ground reflectance  depends on the reflectance capability of the foreground near the site, ranging from 0 to 1. Several reflectance parameters of typical landforms were listed in a previous study [30], where water and dry ground were 0.07 and 0.20, respectively.
The parameters that represent the optical depths  Beam and  Diff are location specific, and the data from 11,124 stations worldwide can be found in the ASHRAE handbook [27].
In terms of reaching solar energy  Sun , the structure surface would absorb a part of this energy while reflecting the rest.The absorbed energy can then yield a correlation expression as where the scale factor  is the absorptivity of the exposed surface.Solar radiation input occurs only during daytime with a clear sky.  is small during night and rainy days and can be neglected.In engineering applications, the daily energy can be summed up on the basis of the hourly  Sun value: (24) where  Rise and  Set denote the sunrise and sunset hours, respectively, and  Sun denotes the time average value of  Sun during daytime.

Water-Side Convection.
Several dams begin to fill the reservoir during the construction phase to improve the performance of power production, which results in the concrete exchanging heat with the reservoir through the waterside surface below the reservoir level.This phenomenon is considered a convection boundary with liquid fluid, which can also be described in Newton's law of cooling: where   denotes the water temperature and ℎ  denotes the convection coefficient between the concrete and water.Coefficient ℎ  is relatively large for water-side convection unlike that for air-side convection.The water-side convection coefficient is 20 times to 1000 times that of the air-side convection [31,32].In this situation, the boundary can be directly dealt as a Dirichlet boundary condition (first-type boundary condition): In engineering applications,   can be obtained by monitoring the water temperature or through the empirical prediction method [33,34].

Surface Protection.
The concrete surface may be covered with insulation protections during construction to prevent thermal cracks induced by atmospheric temperature shock.Materials with high thermal resistances, such as cotton quilts, polystyrene foam plastic (PFP) [35] boards, and sprayed polyurethane, are ideal options for protection in engineering applications.Figure 4 shows two typical thermal protections used in dams.Tiling the cotton quilt on horizontal surfaces, such as the top faces of new-cast lifts, is convenient.And polystyrene foam boards are suitable for pasting clung to the nonhorizontal surfaces.We can consider the boundary condition as a composite wall comprising miscellaneous materials in a series.Similar to the air-side convection, we can apply the following Robin boundary condition: where ℎ  denotes the equivalent convection coefficient of the protection series and is given as the reciprocal of the sum of thermal resistances: where ℎ  denotes the nonprotective convective coefficient in a nonwind condition (  = 0 in (4)),   denotes the width of each protection, and   denotes the thermal conduction coefficient of each insulation layer.

Finite Element Methods
Solving (1) through analytical methods with the aforementioned geometrically complex dam structures and various boundary conditions is impractical.The FE method, which is the most widely applied numerical tool in engineering fields, is chosen as the numerical approach to solve the heat transfer problem of the dam under construction in this study.The 3D space  can be subdivided into numerous elements.Let   be the temperature of an arbitrary location within an element, which can be interpolated as a function of the nodal temperatures on this element mesh: where {}  denotes the nodal temperature arrays and [] denotes the interpolation function matrix.We can also transpose (29) into (1).The classic FE formulation can be formed as where [] and [] denote the heat capacity and thermal conductivity matrices, respectively and {} denotes the heat flux vector: where Given the boundary conditions mentioned in (2), {} can be rewritten as where Γ represents the external surface regions, which have various types.The method of determining the different types of surface regions at different times is proposed in Section 4.
Equation (30) yields the following expression with the help of the backward difference scheme through the discretization of the time dimension: where Δ denotes the increment time and the subscript inc is the time increment.The initial condition and thermal boundary conditions are known.The temperature field of any time can be deduced increment by increment.

Boolean-Based Surface Procedure
The external heat transfer conditions during the construction phase of a large dam are complex.For a proper FE analysis of a dam under construction, the casting of a new lift means activating the lift elements, which changes the external surface regions.In some large projects, the dam is composed of thousands of lifts, which makes the manual selection of different types of surfaces at each step impractical.Therefore, an automatic surface extracting and taxonomy method based on Boolean operations is presented in this study.

Boolean Operations.
Boolean operations are a group of symbolic computation rules used to operate polygons.Three basic operations are always conducted as illustrated in Figure 5: intersection, union, and difference.In the programming of a Boolean-based operation between two surface sets (e.g.,  and ), two labels  and  are attached on each surface element.If the element belongs to the set ,  = 1; otherwise,  = 0. Similarly,  = 1 refers to the elements in the internal area of .Therefore, the rules of the Boolean-based operation can be defined as follows.
(i) Intersection: it is denoted as  ∧ , which generates the set of elements located in the area where  and  overlap.This rule means that the sets are the elements with  +  = 2.
(ii) Union: it is denoted as  ∨ , which generates the set of overall elements, which is a member of , , or both.This rule means that the sets are the elements with  +  ≥ 1.
(iii) Difference: it is denoted as  − , which generates the set of elements located where  remains when (1) Read surface sets already been defined,  1 ,  2 ,  3 and  4 , according to Section 4.

(2)
Detect outlet hole surfaces in the dam body, and store them as  5 , according to Algorithm 3. (3) Detect joint surfaces between adjacent monoliths, and store them as  6 , according to Algorithm 3. (4) Detect foundation surfaces between dam and rock, and store them as S rock, according to Algorithm 3. ( 5) Kill all the dam elements.// the elements whose lift ID > 0. ( 6) for  = 1, 2, . . .,  do (7) Activate the lift cast at this step.// the elements whose lift ID = .(8) Recognize external surfaces of activated elements, store them as S all according to Algorithm 2. (9) Look for the element surfaces with lift ID =  from S all, and note them as S lift.(

mono ID Lift dependent
The variable represents the monolith to which the lift belongs, and mono ID [0] = 0.

Monolith dependent
The variable represents the highest elevation of each monolith during construction, where all variables are initially set as 0.
is substracted.This rule means that the sets are the elements with  −  = 1.Similarly,  −  = 1 yields the result of the  −  operation.

Preliminary Defining Surfaces.
Four simple preliminary surface groups should be identified to obtain precisely different types of external surfaces in the simulation.Figure 6 shows the preliminary surface groups of a typical arch dam layout.
(i) Model-side surfaces: the in-rock boundary surfaces of the model should be identified and noted as  1 .
(ii) Dam-top surfaces: the top surfaces of the dam body should be identified and noted as  2 , and usually those surfaces are on the same elevation.
(iii) Upstream surfaces: the upstream surfaces of the dam should be identified and noted as  3 .
(iv) Downstream surfaces: the downstream surfaces of the dam should be identified and noted as  4 .
Death/rebirth technology is used to simulate the construction process in FE analysis.The entire analysis is divided into numerous steps on the basis of the casting-time plan.The lift elements are activated at the beginning of each step.A   variable is created for each element to store the lift casting sequence.Two additional arrays are used to represent the space relations, as shown in Table 1.

Surface Procedure and FEM Implementation.
We are now ready to present the surface extraction and taxonomy procedures based on Boolean operations.First we set a nomenclature rule of the surface group identification.
(1) Postfix including integer numbers are the surfaces remain the same during the analysis, or we can name them as static surfaces, like  1 ,  2 , and so forth.(2) Postfix including characters are the surfaces change during the analysis, or we can name them as dynamic surfaces, like  up,  down, and so forth.The pseudo-code description of the procedure can be summarized by the following algorithm.
Various types of external heat transfer occur on different surfaces during construction.More than one type of heat transfer may co-exist on one surface.The surfaces stored in Algorithm 1 are an efficient taxonomy method for each step.We can set the boundary condition as summarized in Table 2.A few remarks should be made on the setup rule in Table 2: (i) The in-rock vertical surfaces generated by the modelside cutting are adiabatic because the horizontal heat transfer in the infinite rock is symmetrical to these surfaces.(ii) The ground temperature beneath the upper 6 m of the surface of the Earth maintains a nearly constant temperature between 10 ∘ C and 16 ∘ C [36], depending on the site latitude.Let   denote the upper level of the constant temperature layer, and the thermal boundaries below   on  1 can be expressed as the first-type boundary condition.(iii)  up and  down are used to represent the upstream water level and downstream water level, and the impounding plan should be drawn up when making prediction.(iv) A protection scheme should be established based on the design plan when a prediction is conducted The two faces have the same node list, so they are noted as the internal faces 9, 10, 3 7, 12, 11 6, 11, 9 3, 7, 6 or based on the actual condition when simulating the past.() is the variable that denotes whether protection is covered at time .In large engineering projects, () is also relevant to the surface locations.Some projects may be dependent on the lift and others on the elevation, in which case (,  ) or (, ) is used to represent the protection status.
(v) ℎ ℎ represents the height of the hole.When concrete is cast over the top of the holes, the sun barely shines into the spillway holes constructed inside the dam body; thus, solar irradiation is not considered in these areas.
(vi)   differs from all the surfaces because the inclination of every element face is not the same.
Algorithm 1 mentions that the external surfaces should be initially recognized.Each internal face for the FE meshes must have the same node list with only one face of another element.One way to extract the external surfaces is to check all the internal faces by comparing all the faces of the model.However, this method costs too much time for large models.In this study, a minimum node number is used as a label to accelerate the extracting process, the details of which are described in Algorithm 2. Figure 7 shows an example where two elements are adjacent on one internal face.
The detection algorithm of  5 ,  6 , and  rock mentioned in Algorithm 1 can be executed in Algorithm 3.
A remark should be noticed that the algorithms can also be applied as a preprocessing facility for a third-party FEM software package.An improvement which should be made on the origin Algorithm 1 is that the surface sets extracted at each step should be separately stored and without executing the FEM solver.

Software Development.
In this study, a user-friendly interface for data input is developed based on the EXCEL software.The embedded VBA programming language can be utilized to provide a convenient user interface and to handle the data automatically.Several forms are generated for users to collect the information required in the simulation, including climate condition, construction procedure, material properties, and predefined surfaces.The form style is concise, explicit, and easy-to-use.
A C-language-based program is compiled to process the data input from the spreadsheet and the FE meshes.Algorithm 1 to Algorithm 3 are implemented in the program.
Several commercial or free source FEM packages exist.Our developed program can be exported as requested provided that the input format is known.In this study, MSC.Marc [37] is used as the third-party FEM solver because of its quickness and stable benchmark.
Common users found the postprocessing module of MSC.Marc to be inconvenient and awkward.A program is developed to decode the .t16result output file for fast, automatic, and easy-to-use postprocessing.A customized form, which can be found in the EXCEL interface, offers two options: desired nodes and desired contour time.The developed program can then extract and plot the data, including the time-temperature curves and temperature contours.
Figure 8 summarizes the entire flow process.

Case Study
5.1.Engineering Description.The site of the Xiluodu hydropower project is located in the lower reach of the Jinsha River, Yunnan Province in Southwest China.The geographic coordinate is 28 ∘ 13  02  N, 103 ∘ 33  46  E. The project is designed as a dome dam that is double curved in both horizontal and vertical planes and with a maximum height of 285.5 m and a crest length of 700 m.The footing sits upon the riverbed foundation on a 324.5 m elevation.A total of 31 monoliths are planted across the river.Seven 12.5 m × 13.5 m surface spillways and eight 6 m × 6.7 m deep outlet holes are arranged in the arch dam body.The transversal section of the dam is curve shaped and is 14 m wide at the crest and 60 m wide at the footing.The upstream cofferdam was demolished in March 2012, and the reservoir was impounded.The upstream water level was measured and recorded (see Figure 9).

Climate Conditions.
The site is located in a V-shaped canyon of the Jinsha River and is subjected to a subtropical climate.A nearby hydrologic station records the weather data daily.The ambient temperatures and wind speeds from January 2009 to May 2013 are shown in Figure 10, in which the temperatures over the years can be considered repeated sinusoidal curves on the basis of (7).
A meteorological station is located near the dam site in Zhaotong, and the fundamental solar measurement data were used as the reference for the simulation.Table 3 shows the data for the solar irradiation calculation, which acts as the basic design parameters in the simulation.Table 3 shows the calculation result of the solar radiation, where  * Beam and  * Diff are the radiation values at solar noon,  Beam and  Diff are the average radiation value's during daytime, and   Sun and   Sun are, respectively, the average total solar input for the vertical and horizontal structure surfaces.
Beam and  Diff are the station measurements, based on which the values of  * Beam and  * Diff are calculated by ( 16) and ( 18) by setting  = 0.  Rise and  Set can then be obtained through a multiple trial computation when  reaches 0 ∘ .Averaging  Beam ,  Diff , and  Sun from  Rise to  Set yields the result of  Beam ,  Diff , and  Sun .Recognize all external surfaces of elements with mono ID > 0, and noted them as S base.// Searching for only the base elements based on Algorithm 2, that means the foundation faces the dam sits upon will also be included.(6) S rock =   −  1 .
Algorithm 3: Initial detection of the hole surfaces  5 , transverse joint surfaces  6 , and rock-exposed surfaces  rock.11 shows the instrument layout of the upstream surface of the dam.

Material Properties.
A double exponent hydration model [38] of concrete is applied in this engineering, which can be expressed as where  is the adiabatic temperature rise and  0 , ,  1 , and  2 are the fitting parameters.The technical reports of the material experiments reveal that three main types of concrete are used in the dam construction.The material properties are listed in Table 4.

Simulation Model.
The 3D dam-foundation system is discretized into hexahedral meshes.A total of 154,008 elements represent the dam body, and 28,026 elements represent the foundation.The total number of nodes is 216,837.The FE mesh layout of the dam body is shown in Figure 12.

Results and Discussions
6.1.Surface Taxonomy.The surface procedure mentioned in this study is implemented as a preprocessing facility.The surface groups classified at each step are stored.Figure 13 shows the results of the surface taxonomy map at the 500th and 700th steps, in which different colors represent different surface groups.The FE model is split along the central river to show clearly the upstream and downstream faces.The  (i) External thermal flux has great impact on the area underneath the boundary surfaces, for the temperature varies following the seasonal change.(ii) The arch dam seals the joint layer by layer and starts cooling the sealed regions before the sealing time.We can see that the weather-affected border is obviously narrow in some lifts because of the strong internal cooling source, like (d)∼(f) in the figure.(iii) The after-sealed internal concrete, located a particular height below the top elevation, is not sensitive to the external weather.The temperature does not change much throughout a whole year.(iv) The surface temperature gradient is more obvious in summer than in winter, for the dam's sealing temperature is 13 ∘ C, much lower than the average summer air temperature and closer to the average winter air temperature.(v) The new-cast lift can reach a relatively high temperature peak due to the hydration reaction and induce a temperature gradient with the old concrete and the external surfaces.(vi) The ground temperature has impact on the footing region where the temperature is higher than the sealing temperature.(vii) The upstream face regions which contact with water are mainly affected by the water temperature after March 2012, when the reservoir starts impounding, like (j)∼(i) in the figure.(viii) The blank area is the outlet hole, like (h)∼(i) in the figure.The hollow hole makes the hole wall contact with the air, which generates a temperature gradient near the surface.

Conclusions
The external heat transfer of dams during the construction phase is complex because such transfer is location specific and time varying.A sound model for the external heat transfer of dams during construction is developed in this study.First, a mathematical model of five external thermal flux factors, namely, air-side convection, electromagnetic radiation, absorbed solar input, water-side convection, and surface insulation effect, is proposed.A surface procedure for FE analysis is also established based on Boolean operations.The proposed procedure can extract external surfaces and can create a taxonomy of surfaces for proper boundary conditions.An arch dam that is under construction in China with a height of 285 m is used to test the proposed method.
The method can precisely extract the external surface and can correctly classify these surfaces in each step.The numerical result is compared with the monitoring temperatures.The result shows good agreement, thereby verifying the validity of the proposed method.

Figure 1 :
Figure 1: Ambient temperature sinusoid curve for one day (or one year).

Figure 2 :Figure 3 :
Figure 2: Positions of the Sun and the Earth (this figure shows a winter time for the north because the Earth-Sun line crosses the Southern Hemisphere).

Figure 4 :
Figure 4: Surface protection and its equivalent model.

Figure 7 :
Figure 7: Example of recognizing the inside faces.The red bold number is the minimum node number of the face.

Figure 8 :
Figure 8: Flowchart of the software development.
Time when 1000th lift cast

Figure 13 :
Figure 13: Sketch of surface taxonomy map.

( 1 )
for  = 1, 2, . . .,   do // Searching all the elements (2) Request the element type and the count of faces in one element is stored as .(3) for  = 1, 2, . . .,  do // Searching all the faces (4) Obtain the minimum node number of the face nodes, labeled as  , .(5) end loop  (6) end loop  (7) Search for the faces with the same value of  , .(8) Compare the topology of the faces.If two faces have the same nodes list, note them as internal faces.(9)All the faces that have not been noted can be stored as the external surfaces.Algorithm 2: Recognition of the external surfaces.

Figure 15 :
Figure 15: Temperature contour of the section in monolith 15 throughout a year from June 2011 to May 2012.
The label of the casting sequence and lift ID = 0 represent the nonconcrete elements (i.e., rock base).

Table 2 :
External surface taxonomies for boundary conditions setup.

Table 3 :
Calculation of solar irradiation.
mometers are arranged at the same elevation in the same concrete lift: one at the upstream surface and the other at the downstream surface.Monoliths number 6, number 9, number 16, number 22, and number 27 are selected as monitoring sections.Figure

Table 4 :
Material properties of concrete.June 2011 to May 2012.Conclusions can be drawn from these contour maps as follows. from