Computer Simulation of Artificial Ground Freezing during the Building of Underground Tunnel
In Germany, during the construction of a new intersection between the regional carrier, Deutsche Bahn, and a city underground railway line, a new railway section with a total length of 1.3 kilometers was required. A specific feature of this section is that its path crosses a large complex of buildings protected by law and it was therefore decided to construct a tunnel at this section. However, this task was seriously complicated by high groundwater flow rates in one segment of this section. As a result of the high filtration rate during tunnel drilling in this segment, there is a distinct probability that the earth under the foundations will be subject to settling. This situation may compromise any structures built on top of this underlying area, which is of course, unacceptable.
As a solution to this nonstandard problem, an artificial ice shield over the drilling area was proposed. This measure allows the strengthening of the ground and securing the foundations while performing construction of the tunnel. The tunnel consists of two paths, each of which requiring separate drilling operation. The image below shows a vertical section of the ground in the tunnel drilling area.
The blue color shows the area to be frozen.
In order to justify the application of the proposed design solutions, it was necessary to perform computer simulation of artificial ground freezing in the tunnel construction area using specialized software.
When building a computer model, the following geometrical dimensions of the computational domain were set: • Length – 62 meters. • Width – 20.2 meters. • Height – 18.4 meters.
The magnitudes of the computational domain were chosen on the basis of the geometric dimensions of the tunnel and cooling devices and the geological and lithological structure of the ground. The influence of the boundaries on the thermal processes inside computational domain was also factored in.
Ground layers in the area under consideration have a complex structure. A geotechnical survey was carried out and the following table with borehole heights was obtained.
The plan of wells is shown in the following image:
Geostatistical tools were used to create a threedimensional reconstruction of the ground within the computational domain, as shown in the image below. This also shows the area of tunnel paths and freezing pipes to be laid.
23 freezing pipes with a total length of more than 1,200 meters were proposed for ground freezing and the creation of an ice shield. The diameter of each pipe is 3.5 inches. The cooling performance of the freezing aggregate is 465 kW, which allows the temperature of 40 ^{o}C. A 30% calcium solution was selected as the freezing agent.
Collocation of the freezing pipes is presented in the section diagram below.
The drilling of pipeholes is executed at a slight angle. A summary table of the initial and final depths of the freezing pipes to be laid is presented below.
Water filtration velocity in sand layer comprises 1.25 meters per day; the groundwater temperature is +12.9^{o}C.
Ground characteristics are shown in the table below.
Ρ_{dry} – ground density; k_{f} – filtration coefficient; λ_{2,1} – heat conductivity in the thawed and frozen ground; C_{2,1} — heat capacity of thawed and frozen ground; α, β – coefficients that describe the content of unfrozen moisture.
The Firsttype boundary condition is set on the surface of the computational domain as a dependence of temperature change over time, shown in the diagram below.
The Secondtype boundary condition in the form of heat flow equal to zero is set on the side boundary of the computational domain. The Firsttype boundary condition representing a constant temperature of +12.9 ^{o}C is set on the lower boundary.
As a result of computational domain discretization, we obtain a computational mesh consisting of 2,037,312 nodes.
The simulated period was 120 days (from February to June); the time spent to carry out the computation was 35 minutes
Cross section of the middle of the computational domain on the 30th day of operation of the cooling devices
The simulation results of the performance of the cooling devices showed that on the 120th day, the ice shield over both paths was fully formed. Moreover, the first path was almost completely frozen by this time. This can be clearly seen in the image representing the computation results below.
Simulation results on the 120th day of operation of the cooling devices
Isolines of the YZ plane section on the 120th day of computation
You will be also interested in: Computer Simulation of AGF during the Building of Tunnel (Article in PDF)


Examples of projects: