Air and Thermal Coupling
July 7, 2010
SoilVision systems Ltd. is pleased to announce the addition of full air and thermal coupling between the SVAirFlow® and SVHeat® Professional software packages. The addition extends the capabilities of the software to model flow of air in the pore space of soils. In particular, the flow of air as a result of thermal gradients can be accomplished. This addition allows benefits in modeling the following applications:
- Heap leaching: certain heap leaching operations consume oxygen during the reaction process. There is therefore the danger that the reactions will become oxygen starved and extraction will not reach capacity. Therefore air injection systems are sometimes introduced to enhance recovery. The performance of such systems can be modeled with the combined SVAirFlow® and SVHeat® software packages.
- Convective air flow in embankments: thermal convection in her embankments can now be modeled (Goering, 2003)
- Waste rock air movement: the reactions in many waste rock piles often produce heat which then affects the flow of air in such piles. The numerical modeling of such a process can aid in the long-term design of remediation scenarios for waste rock piles.
In order to examine this coupling further the following benchmark is presented. Further details of this model may be found in the SVAirFlow® Verification Manual. Goering (2003), Sun et al (2007), and Lebeau (2009) have conducted the numerical simulation of cooling an embankment using natural convection. The intent is to protect railway or roadway construction from thaw in permafrost area. This example presents the characteristics of air flow and temperature distribution as governed by natural convection within the embankment using the coupled software.
The model geometry of the right half of a symmetric roadway is shown in Figure 1, below. The embankment is constructed with ballast type of material that has a high air permeability. The pore-air in the embankment is allowed to approach equilibrium with the surrounding atmosphere. The climate temperature changes with time. The ground foundation is assumed to be uniform silt material with a high water content. On the natural ground surface, the air transfer on the boundary is assumed to impermeable, and the climate temperature varies with time. The air is assumed to be incompressible in natural convection.
Figure 1: Model Geometry
The model is simulated for 21 years to eliminate the effect of the initial condition on the temperature in ground foundation. The following briefly discusses and compares the results with other software.
The patterns of convective air flow in the embankment
During winter, the cool air at the embankment surface enters into the embankment. This happens because the air density increases with the decrease in temperature. The warm air then flows out from the lower region in the embankment to the atmosphere by the process of natural convection, as illustrated in Figures 2 to 6. The onset of natural convection occurs between October and November. After that time the extensity of convective airflow increases with the decrease in temperature until January. The convective flow reduces with the increase in temperature after the end of January. During the summer, because the temperature on the embankment surface is higher than the temperature at the lower region in the embankment, a small amount of convective airflow occurs in the embankment (see Figures 7 to 9).
Compare convective flow patterns with results from others
Figure 10 is the airflow simulated from SVAirflow/SVHeat, and Figure 11 is from Goering (2003). It can be seen that the patterns of convective airflow at different times are very similar. It should be noted that the figures from Goering example are presented only as a comparison of the flow pattern. Exact comparison to the flow vectors is not possible as Goring did not published enough specific of the modeling process. As illustrated in several of the above figures, multi-airflow loops are formed in wintertime. For the two neighboring convective airflow loops, if one flows clockwise, another must move anti-clockwise. This characteristic of natural convection was also discussed in other research (Chen et al., 2008). Figure 10 further illustrates the comparison of the average value of airflow velocity at various locations (see circle points in Figure 1) with the value obtained by Goering. It should be noted that the airflow velocity for Goering in Figure 11 is the average value of airflow velocity magnitude over the whole embankment. This may be why the airflow velocity simulated in this benchmark in the SVAirFlow software has greater variation than the Goering’s value. It should be noted that the purpose of this benchmark is to replicate the reasonable flow patterns published by Goering (2003). An exact duplication of flow vectors is not needed for this benchmark to be deemed successful.
Figure 2: Patterns of convective airflow during winter, mean air velocity = 3.2 m/s
Figure 3: Patterns of convective airflow during winter, mean air velocity = 6.26 m/s
Figure 4: Patterns of convective airflow during winter, mean air velocity = 12.69 m/s
Figure 5: Patterns of convective airflow during winter, mean air velocity = 17.78 m/s
Figure 6: Patterns of convective airflow during winter, mean air velocity = 6.57 m/s
Figure 7: Patterns of convective airflow during summer, mean air velocity = 4.97 m/s
Figure 8: Patterns of convective airflow during summer, mean air velocity = 6.29 m/s
Figure 9: Patterns of convective airflow during summer, mean air velocity = 5.36 m/s
Figure 10: Comparison of the patterns of convective airflow with the patterns from Goering, mean air velocity = 6.24 m/s
Figure 11: Patterns of convective airflow presented by Goering (2003)
SoilVision Systems Ltd. distributes a large database of numerical models. The above model can be found in our publicly distributed models under:
Other examples that illustrate this process can be seen in the following distribution models which can be found in the SVAirflow examples manual:
- Project: Geotherm Model: Embankment_Convective_Airflow
- Project: Geotherm Model: ConvectiveBoxHeatedBelow
Purpose: A simple box model is used to illustrate the convective airflow when a warm temperature (TH) is introduced to the bottom of the box and a cooling temperature (TL) is applied to the top. This example uses a SVAirFlow and SVHeat coupled model to simulate natural convection of air or thermo-buoyant motion under the different temperatures TH and TL.
Figure 12: Temperature contour of natural convection in the case of TH = 30 °C and TL = 0 °C
Figure 13: Convective air-flow in the case of TH = 30 °C and TL = 0 °C
Figure 14: Temperature contour of natural convection in the case of TH = 50 °C and TL = 0 °C
Figure 15: Convective air-flow in the case of TH = 50 °C and TL = 0 °C
- Project: Geotherm Model: ConvectiveBoxHeatedFromSide
Purpose: This example illustrates the thermo-buoyant motion when a box is heated from the right or left side.
Figure 16: Temperature contour heated from side (temperature at right side 15 °C and left 0 °C)
Figure 17: Anti-clockwise convective air-flow heated (temperature at right side 15 °C and left 0 °C)
Figure 18: Temperature contour heated side (temperature at right side 0 °C and left 15 °C)
Figure 19: Clockwise convective air-flow heated side (temperature at right side 0 °C and left 15 °C)
For more information on this feature, feel free to contact us directly, or visit our online store to purchase SVHeat® and SVAirFlow®.