Optimization of Geothermal Heat Pump Systems for Sustainable Urban Development in Southeast Asia

: This study examines the optimization of ground source heat pump (GSHP) systems and energy piles for sustainable urban development, focusing on Southeast Asia. GSHPs, which utilize geothermal energy for indoor HVAC needs, offer a sustainable alternative to traditional systems by utilizing consistent subsurface temperatures for heating and cooling. The study highlights the importance of understanding thermal movement within the soil, especially in soft marine clays prevalent in Southeast Asia, to improve GSHP system efficiency. Using a one-dimensional finite difference model, the study examines the effects of soil thermal conductivity and density on system performance. The results show that GSHP systems, especially when integrated with energy piles, significantly reduce electricity consumption and greenhouse gas emissions, underscoring their potential to mitigate the urban heat island effect in densely populated areas. Despite challenges posed by the region’s hot and humid climate, which could affect long-term effectiveness, the study highlights the need for further study, including field experiments and advanced modeling techniques, to optimize GSHP configurations and fully exploit geothermal energy in urban environments. The study’s insights into soil thermal dynamics and system design optimization contribute to advancing sustainable urban infrastructure development.


Introduction
In the construction of smart cities, implementing energy piles and assessing heat transfer within the ground is paramount for efficient and sustainable urban development.The use of shallow geothermal systems, also known as ground source heat pump (GSHP) systems, utilizes the geothermal energy available beneath the surface to provide either heating or cooling resources.This technology meets the thermal requirements for indoor Heating, Ventilation, and Air-Conditioning (HVAC) systems, as described by Sarbu and Sebarchievici [1].GSHP systems can be configured in several forms, including vertical GSHP systems integrated into the pile foundations, which are essential for energy piles in smart cities, horizontal GSHP systems, and GSHP systems that use a reservoir for cooling.For smaller residential buildings, a closed-loop system is often recommended due to its ease of maintenance and operation, as shown in Figure 1.Shallow GSHP systems, defined by Eugster and Sanner [2] as operating at a depth of 3 to 100 m, use the soil mass as a thermal reservoir.This approach provides a consistent temperature source throughout the year, in contrast to the variable air temperatures at the surface [3,4].Shallow geothermal systems are primarily used for heating and cooling by exploiting subsurface temperatures below 25 • C, making them suitable for small-scale and residential applications in various locations.Conversely, deep geothermal systems target temperatures above 25 • C and can reach up to 200 • C, serving different energy needs [2].various locations.Conversely, deep geothermal systems target temperatures above 25 °C and can reach up to 200 °C, serving different energy needs [2].This study explores the use of GSHP systems, with a focus on shallow geothermal systems, for smart city development.GSHPs use the constant temperature below the earth's surface for heating or cooling, efficiently meeting indoor HVAC requirements.The study emphasizes the integration of energy piles-structures that both support buildings and function as part of a GSHP system-and highlights their potential in sustainable urban construction.By evaluating thermal movement within the soil, particularly in soft marine clays common in Southeast Asia, the study advances the understanding of how GSHP systems can be optimized for tropical climates where the urban heat island effect is a significant challenge.A notable aspect of this study is the use of a one-dimensional finite difference model to predict the temperature distribution around energy piles over time, allowing for a simplified yet effective analysis.This approach facilitates the evaluation of the effect of soil thermal conductivity and density on GSHP system efficiency.The results indicate that soil density affects the thermal performance of the system.Denser soils improve thermal conductivity.By focusing on the practical application of GSHP systems in smart cities, especially in regions such as Bangkok with specific geotechnical challenges, this study provides valuable insights into the deployment of sustainable energy solutions.It highlights the need for further study to optimize GSHP configurations and operating parameters to realize the full potential of geothermal energy in urban development.The study's implications for smart city construction include the reduction of electricity consumption, mitigation of greenhouse gas emissions, and improvement of HVAC system efficiency, demonstrating a promising path towards sustainable urban living environments.
While existing studies focus primarily on the technical efficiency and economic aspects of GSHP systems in temperate climates, this study extends existing studies to tropical regions, addressing a significant knowledge gap.Unlike previous studies that rarely evaluate the interaction between soil thermal properties and GSHP performance in hot and humid climates, this study uses a one-dimensional finite difference model to analyze the effects of soil thermal conductivity and density on system efficiency.This approach is novel in its application to the soft marine clays prevalent in Southeast Asia, which differ significantly from the more commonly studied temperate zone soils.By integrating energy piles (which are typically used as part of a building's foundation) into the GSHP system, this study utilizes a dual-function infrastructure that enhances both structural This study explores the use of GSHP systems, with a focus on shallow geothermal systems, for smart city development.GSHPs use the constant temperature below the earth's surface for heating or cooling, efficiently meeting indoor HVAC requirements.The study emphasizes the integration of energy piles-structures that both support buildings and function as part of a GSHP system-and highlights their potential in sustainable urban construction.By evaluating thermal movement within the soil, particularly in soft marine clays common in Southeast Asia, the study advances the understanding of how GSHP systems can be optimized for tropical climates where the urban heat island effect is a significant challenge.A notable aspect of this study is the use of a one-dimensional finite difference model to predict the temperature distribution around energy piles over time, allowing for a simplified yet effective analysis.This approach facilitates the evaluation of the effect of soil thermal conductivity and density on GSHP system efficiency.The results indicate that soil density affects the thermal performance of the system.Denser soils improve thermal conductivity.By focusing on the practical application of GSHP systems in smart cities, especially in regions such as Bangkok with specific geotechnical challenges, this study provides valuable insights into the deployment of sustainable energy solutions.It highlights the need for further study to optimize GSHP configurations and operating parameters to realize the full potential of geothermal energy in urban development.The study's implications for smart city construction include the reduction of electricity consumption, mitigation of greenhouse gas emissions, and improvement of HVAC system efficiency, demonstrating a promising path towards sustainable urban living environments.
While existing studies focus primarily on the technical efficiency and economic aspects of GSHP systems in temperate climates, this study extends existing studies to tropical regions, addressing a significant knowledge gap.Unlike previous studies that rarely evaluate the interaction between soil thermal properties and GSHP performance in hot and humid climates, this study uses a one-dimensional finite difference model to analyze the effects of soil thermal conductivity and density on system efficiency.This approach is novel in its application to the soft marine clays prevalent in Southeast Asia, which differ significantly from the more commonly studied temperate zone soils.By integrating energy piles (which are typically used as part of a building's foundation) into the GSHP system, this study utilizes a dual-function infrastructure that enhances both structural support and energy efficiency.The results illustrate how variations in soil properties can affect system performance, which is a critical consideration for urban planners and engineers seeking to deploy these systems in similar urban settings.The originality of this study lies in its methodological innovation and its focus on a region where high ambient temperatures and urban heat island effects pose unique challenges.The integration of field experiments and advanced modeling provides new insights into the optimization of GSHP configurations, which is critical to realizing the full potential of geothermal energy in urban sustainability initiatives.This study not only fills an important research gap but also sets the stage for future innovations in geothermal energy use in tropical urban environments.

GSHP Systems and Thermal Movement Evaluation
Especially in the integration of energy piles, the evaluation of thermal movement in the ground becomes crucial.This evaluation is essential for effectively operating GSHP systems in cooling modes.The core architecture of the GSHP system consists of three essential elements: the heat source, the heat sink, and the heat exchanger.Specifically, the system's heat source is derived from the operation of the HVAC system that regulates the building's indoor environment.Conversely, the ground acts as a heat sink, absorbing excess heat.The heat exchanger transfers heat from the source to the sink.The heat transfer fluid, typically a water-antifreeze solution, is circulated by a pump, facilitating heat transfer by conduction and convection.Heat exchangers can be divided into two types: open-loop and closed-loop systems.An open-loop system uses groundwater drawn from wells to cool and returns the heated water to the aquifer.In contrast, a closed-loop system uses a network of High-Density Polyethylene (HDPE) pipes in series or parallel configuration, as shown in Figure 2 [5].The most basic design within this system is the U-shaped exchanger, which is vertically embedded in boreholes or integrated into the building's foundation piles, providing a path for the circulating medium.To improve heat dissipation, various configurations are used, such as the double U-pipe in parallel or cross arrangement or the spiral/coil structure, to maximize the contact surface area for efficient heat exchange.critical role of accurately assessing the potential of geothermal energy in urban planning and construction.

Regional Focus on Europe and Southeast Asia
The origin of GSHP technology can be traced to the European region, where GSHPs are mostly operated in heating mode to provide heat to the heating infrastructure.Sayegh et al. (2019) [10] discussed how the increasing popularity of district heating in Europe had a positive impact on the development and contributed to the advancement of heat pump technology, of which GSHP is a subset.The integration of heat pump technology into district heating networks led to an increase in efficiency and flexibility, resulting in a lower level of carbon emissions.The four basic heat pump integration scenarios are as follows: (i) heat pump integration into the existing network without major changes, (ii) heat pump integration into an extended network, (iii) deep renovation of the existing district heating system, and (iv) the design of a new district heating system supplied by a heat pump.In the current situation, heat pumps are considered a renewable energy technology in Europe, where they account for 5% to 20% of the European renewable energy targets.The success of this policy implementation will help Europe to achieve its sustainability goals faster, to the benefit of its inhabitants.
The viable heat source for GSHP systems used in residential and commercial buildings in Europe mainly comes from the shallow soil layer.Brys et al. (2019) [11] concluded

Energy piles in parallel configuration
The single energy pile with a U-shape heat exchanger Understanding and evaluating thermal movement in the ground is paramount in implementing energy piles.This understanding hinges on the concept of thermal conductivity, a material property that indicates how well heat can move through a substance.Clay soil, with thermal conductivity values ranging from approximately 1 to 2.5 W/(m•K), has a significantly higher capacity for heat transfer than air, which has a thermal conductivity of only 0.026 W/(m•K), demonstrating the superior efficiency of the ground in conducting heat by a factor of nearly 40 [6].This characteristic underscores the GSHP system's environmental benefits, which significantly reduces electricity consumption and greenhouse gas emissions compared to traditional air source heat pump systems.The versatility of GSHPs makes them an optimal choice for various settings, including residential, commercial, and industrial buildings, as highlighted by Tissen et al. [7].Specifically, GSHP systems are distinguished by their quieter operation, longer lifespan, and lower maintenance requirements than their air source counterparts.Although the initial investment for GSHP installations is generally higher than for air source systems, the long-term energy savings achieved typically offset the upfront costs.From 1995 to 2020, there was a remarkable increase in the adoption of ground source heat pumps, indicating a growing degree of recognition of their benefits and efficiency.According to the most recent data, the leading countries in terms of installed GSHP capacity include China, the United States, Sweden, Germany, and Finland, with capacities of 26,450, 20,230, 6680, 4400, and 2300 MWt, respectively [8,9].This global trend underscores the growing appreciation of GSHP systems as a sustainable and efficient solution in the development of smart cities and underscores the critical role of accurately assessing the potential of geothermal energy in urban planning and construction.

Regional Focus on Europe and Southeast Asia
The origin of GSHP technology can be traced to the European region, where GSHPs are mostly operated in heating mode to provide heat to the heating infrastructure.Sayegh et al. (2019) [10] discussed how the increasing popularity of district heating in Europe had a positive impact on the development and contributed to the advancement of heat pump technology, of which GSHP is a subset.The integration of heat pump technology into district heating networks led to an increase in efficiency and flexibility, resulting in a lower level of carbon emissions.The four basic heat pump integration scenarios are as follows: (i) heat pump integration into the existing network without major changes, (ii) heat pump integration into an extended network, (iii) deep renovation of the existing district heating system, and (iv) the design of a new district heating system supplied by a heat pump.In the current situation, heat pumps are considered a renewable energy technology in Europe, where they account for 5% to 20% of the European renewable energy targets.The success of this policy implementation will help Europe to achieve its sustainability goals faster, to the benefit of its inhabitants.
The viable heat source for GSHP systems used in residential and commercial buildings in Europe mainly comes from the shallow soil layer.Brys et al. (2019) [11] concluded the results of a thermal analysis of over 10 years (August 2007 to July 2017).They collected the measurements of net solar radiation (the difference between incoming solar radiation and outgoing reflected solar radiation) data of positive heat fluxes (G > 0) and negative heat fluxes (G < 0) in bare and grassy soil in Wroclaw, Poland.Heat fluxes (a measure of how much heat is transferred per unit area per unit time) varied seasonally and annually, with bare soil showing higher heat flux values than grass-covered soil.The ten-year average sums for negative soil heat fluxes, which are suitable for GSHP operation in heating mode, reached 91.8 kWh/m 2 for the bare soil and 47.4 kWh/m 2 for the turfed area.The negative heat flux reaches its minimum value in May and its maximum value in October.This confirms that the shallow soil layer has the potential to be an efficient and environmentally friendly heat source for GSHP operation throughout the year.
Focusing on the Southeast Asia region, an area that has recently begun to attract attention for exploring the implementation of GSHP systems, the effects of urbanization come to the fore.The region's rapid urban expansion, coupled with an influx of population into densely populated urban centers, contributes to the urban heat island effect.This phenomenon results in urban areas experiencing significantly higher temperatures than their rural counterparts due to the heat-trapping properties of buildings, roads, and other infrastructure.Such conditions not only reduce the efficiency of conventional air source HVAC systems, but also increase their energy consumption.
Due to the hot and humid climate, most buildings in Southeast Asia are equipped with air-conditioning systems.In urban areas, the density of these units is particularly high, leading to a concentration of noise.Poor maintenance of air conditioners can result in increased noise due to worn components, loose parts, or inefficient operation, combined with inadequate noise regulations and poorly enforced codes.This allows building owners and operators to neglect noise mitigation measures, exacerbating the problem.The outdoor installation of the air conditioner's heat exchanger, cooling tower, and compressor is located in confined spaces and close to residential areas, often on street-facing walls or in narrow alleys.This amplifies the noise and directs it into public spaces.This means that noise from AC units can easily reach residential areas, disturbing residents and affecting their quality of life.To measure the efficiency of the GSHP compared to the traditional air source heat pump, it would be necessary to use the index numbers.The Coefficient of Performance (COP) is the ratio of energy output (watts unit) to energy input (watts unit).A higher COP means the air conditioner can provide more cooling while using less energy at a given point in time under ideal conditions.A higher COP means the heat pump is producing more heat output for a given amount of electrical input, indicating higher efficiency, translating to lower electricity consumption, resulting in lower energy bills.On the other hand, the Seasonal Coefficient of Performance (SCOP) is a more comprehensive index that considers a heat pump's efficiency over an entire cooling season that lasts several months.It takes into account varying temperatures and weather conditions.A higher SCOP indicates better overall efficiency throughout the year and helps to reduce the greenhouse gas emissions and overall environmental impact.Another commonly used measurement index is the Energy Efficiency Ratio (EER).This metric is generally used to measure the cooling efficiency of a heat pump or GSHP at a given point in time.It is the ratio of cooling output (BTU unit) to energy input (watts unit).A higher EER indicates that the air conditioner can provide more cooling while using less energy.The Seasonal Energy Efficiency Ratio (SEER), similar to SCOP, is an index that evaluates the cooling efficiency of a heat pump over an entire cooling season, taking into account varying temperatures and conditions.A higher SEER indicates better cooling performance throughout the year.
COP/SCOP and EER/SEER are typically compared between ASHP and GSHP.The GSHP yields the higher score in all COP, SCOP, EER, and SEER indexes compared to traditional ASHP at 2.5 to 4.0 compared to 3.1-5.0 in the COP index, 2.5 to 3.5 compared to 3.5-4.5 in the SCOP index, 8.0 to 12 compared to 14.0-20.0 in the EER index, and 13.0 to 20 compared to 16.0-30.0in the SEER index.This is because the ground maintains a more constant temperature than the air, resulting in more efficient heat exchange.However, GSHPs typically have higher installation costs due to the need for ground loop systems but come with the benefit of lower operating costs compared to ASHPs, while ASHPs are usually cheaper to install but may have higher operating costs.

Challenges and Prospects in Southeast Asia
The current issue of energy piles and ground source heat pump (GSHP) systems in Southeast Asia is still in its early stages, with minimal research on both practical and simulated models.Given the region's high ambient temperatures, there are concerns about the efficiency and environmental impact of GSHP installations from both a technical and environmental perspective.As a result, available data is scarce.However, there is considerable room for development and these systems may become commercially viable.In the future, I expect the popularity of GSHPs to increase in Southeast Asia.Because of the high initial prices, they may be used in high-rise and large structures at first, but I think they will eventually become ubiquitous in residential areas as part of standard air-conditioning systems.
Faced with these challenges, GSHP systems emerge as a potential solution to mitigate the adverse effects of urban heat islands.However, the tropical, consistently hot, and humid climate of Southeast Asia presents a unique technical hurdle.The relatively low ground-to-air temperature differential characteristic of the region suggests that GSHPs would operate predominantly in cooling mode.As a result, the long-term effectiveness and efficiency of GSHP systems could be compromised as the ground temperature could eventually converge with the ambient air temperature, reducing the cooling potential of the system.Despite these challenges, the volume of field experiments and studies specific to GSHP performance in such climates remains limited.This gap underscores the need for a comprehensive study to evaluate the viability and optimization of GSHP systems under tropical conditions.

Case Studies and Soil Analysis in Southeast Asia
Figure 3 [4] shows the geographical distribution of studies on GSHP systems in the Southeast Asian region, showing a wide range of studies in Thailand [12][13][14][15][16], Indonesia [17], Malaysia [18], Singapore [19], and Vietnam.In particular, the study conducted in Kamphaeng Phet by Yasukawa et al. [16] is recognized as the first investigation of GSHP technology in the region, highlighting the potential benefits of such applications at the study site.These studies were conducted primarily at sites underlain by soft marine clay deposits, with soil samples collected at depths ranging from 5 to 20 m below the surface.Table 1 shows that the clay samples have a similar specific gravity (Gs), liquid limit (LL), plastic limit (PL), and sand-silt-clay percentage content to Bangkok clay.As shown in Table 1 [20][21][22][23][24][25][26][27][28][29][30][31][32], a commonality among these samples is their classification as marine deposit clays.Such clays, formed by sedimentation processes in marine contexts such as deltaic and coastal plains, include notable varieties such as Ariake clay from Japan, Singapore marine clay, Hong Kong marine clay, Kuantan clay from Indonesia, and Ho Chi Minh City marine clay.The latter has significant engineering similarities to Bangkok's soft clay, characterized by high plasticity, low permeability, and significant compressibility.These soil characteristics pose substantial challenges to construction efforts, particularly in urban areas where soil conditions play a critical role in ensuring the structural integrity of densely populated building projects.The inclusion of GSHP systems in such environments requires a thorough understanding of local geotechnical conditions to optimize design and functionality, underscoring the critical importance of evaluating thermal movement on the ground as part of smart city infrastructure development.The focus is on the use of energy piles and shallow GSHPs where soft marine clays are present.For example, one location that has started to gain attraction is around Bangkok, where the study of GSHP system implementation has begun.If the result of the study in Bangkok is obtained, it can be used in other places with similar soil.The background temperature of the soil in the Bangkok area was studied in 2011 [15].The temperature was obtained from 56 temperature depth profiles throughout the central region of Thailand; the measurements were taken at every 2 m interval with an accuracy of 0.01 • C. The temperature at shallow depths from 0 to 50 m below the ground surface is in the range of 27 • C to 30 • C, which is, on average, 5 • C lower than the monthly average maximum atmospheric temperature.
Smart Cities 2024, 7, FOR PEER REVIEW 7 U-tubes installed at the 50 m mark.Notably, this GSHP system was integrated into the air-conditioning infrastructure of a compact building on the Bangkok campus known as "Parot Racha".The results of this study showed that the GSHP system was more efficient than traditional air source heat pump systems, achieving a reduction in electricity consumption of approximately 30% until the heat pump inlet temperature reached 5 °C above ambient.This result underscores the significant potential of GSHP technology to reduce long-term energy consumption in buildings despite Bangkok's elevated daytime temperatures.Further progress in 2020 has been documented by Jotisankasa et al. [13], who conducted a comprehensive energy pile study in a residential building.The study used a 17.5 m-long, 152 mm diameter steel pipe pile equipped with a U-tube heat exchanger in a parallel configuration.The GSHP setup consisted of two different loops, labeled number 1 and number 2, with the former consisting of six interconnected piles and the latter consisting of seven.Preliminary evaluations of the GSHP efficiency, as measured by the thermal conductivity of the system, yielded values of 2.01 W/m-K for loop number 1 and 1.49 W/m-K for loop number 2. These evaluations were conducted with the piles positioned in the stiff clay layers of Bangkok, highlighting the significant impact of groundwater infiltration on the performance of the GSHP system.This study underscores the critical role of assessing thermal dynamics within the soil for the effective implementation of GSHP technologies in smart city developments, particularly in energy pile designs.For each clay property, Ariake clay refers to [20][21][22], Singapore marine clay to [23,24], Hong Kong marine clay to [25], Kuantan clay to [26,27], and Bangkok clay to [28][29][30][31][32], respectively.

Soil
In 2020, a collaborative study effort between Japanese and Thai academics, highlighted by the work of Shimada et al. [14], was conducted at the facilities of Chulalongkorn University to conduct both experiments and simulations on an operational GSHP system.The investigation included a Thermal Response Test (TRT) to validate the simulation model within a pilot plant setup.This GSHP configuration included two boreholes 6 m apart, each 50 m deep and 0.23 m in diameter.The first borehole was equipped with two sets of U-tubes at 10 and 15 m depths connected in series, while the second borehole had U-tubes installed at the 50 m mark.Notably, this GSHP system was integrated into the air-conditioning infrastructure of a compact building on the Bangkok campus known as "Parot Racha".The results of this study showed that the GSHP system was more efficient than traditional air source heat pump systems, achieving a reduction in electricity consumption of approximately 30% until the heat pump inlet temperature reached 5 • C above ambient.This result underscores the significant potential of GSHP technology to reduce long-term energy consumption in buildings despite Bangkok's elevated daytime temperatures.Further progress in 2020 has been documented by Jotisankasa et al. [13], who conducted a comprehensive energy pile study in a residential building.The study used a 17.5 m-long, 152 mm diameter steel pipe pile equipped with a U-tube heat exchanger in a parallel configuration.The GSHP setup consisted of two different loops, labeled number 1 and number 2, with the former consisting of six interconnected piles and the latter consisting of seven.Preliminary evaluations of the GSHP efficiency, as measured by the thermal conductivity of the system, yielded values of 2.01 W/m-K for loop number 1 and 1.49 W/m-K for loop number 2. These evaluations were conducted with the piles positioned in the stiff clay layers of Bangkok, highlighting the significant impact of groundwater infiltration on the performance of the GSHP system.This study underscores the critical role of assessing thermal dynamics within the soil for the effective implementation of GSHP technologies in smart city developments, particularly in energy pile designs.
Understanding the role of thermal movement in the soil is critical to building smart cities and integrating energy piles.Thermal conductivity, a determinant of how well heat moves through the soil, varies among different soil types, even those with similar engineering properties.The main factors affecting thermal conductivity are moisture content and density.In Southeast Asia, natural soft clays at shallow depths are often completely saturated with water.High moisture content and saturation levels reduce the presence of air, which acts as an insulating layer in the voids between soil particles.In addition, the movement of water through these voids enhances heat convection, resulting in increased thermal conductivity.Similarly, the density of the soil plays an important role; denser soil allows for better contact between particles, which increases heat conduction and, consequently, thermal conductivity.The application of GSHP and energy pile technologies in Southeast Asia is still in its infancy, indicating a significant area for study and exploration.The specific effects of variations in the effective thermal conductivity of GSHP systems, as well as the thermal conductivity and density of soft clays, on temperature modulation and heat dissipation in the soil, remain poorly understood.Thailand, which is at the forefront of the GSHP study and has extensive data on the characteristics of Bangkok soft clay, provides an excellent opportunity for detailed case studies.Such a study could provide valuable insights applicable to other Southeast Asian regions with similar geotechnical characteristics, thereby advancing GSHP feasibility studies, system designs, and performance evaluations.This effort is essential to realize the full potential of GSHP systems in supporting the sustainable infrastructure of smart cities.

Fundamental Theory of Heat Transfer in Soil Using Finite Difference Method
Typically, preliminary design heat transfer calculations from the energy foundation to the surrounding soil mass are evaluated based on the magnitude of temperature changes in the soil at various distances from the energy foundation.These calculations are critical for predicting the performance of ground source heat pump (GSHP) systems before the foundation design is finalized.
The basic theories needed to understand the finite difference method for studying soil heat transfer can be summarized as follows: 1.The definition and advantages of the finite difference method; 2. The model simplification and boundary conditions; 3. The model computational domain and computational nodes (meshing); 4. The transformation from differential equation to algebraic equation solution.
The analysis technique using the finite difference method (FDM) offers several advantages, including a relatively simple calculation procedure.It uses algebraic equations instead of the complex heat transfer differential equations typically used in other methods.In addition, this analysis can be performed using a simple spreadsheet, using a mathematical equation derived from a more complex differential equation that can be solved by entering specific parameters.Finite difference analysis allows for the calculation of the time-dependent temperature evolution of the soil around the energy pile, with the scale of the time step adjustable to seconds, minutes, hours, or days, depending on the study requirements.Given the time constraints of this study, the finite difference method is identified as the most appropriate and direct approach.
In this study, the finite difference model of the energy piles is conceptualized as elongated cylindrical shapes within an infinite soil mass, simplifying the three-dimensional calculations to a one-dimensional analysis focused on radial temperature changes from the center of the energy pile outward.The general domain of the model is a vertical plane bisecting the energy pile, containing a heat exchanger that maintains a constant temperature, simulating the heating mode of the GSHP system.Heat transfer is radial, moving outward from the heat exchanger to the infinite groundmass on three different soft marine clay densities.Before starting the calculations, it is necessary to define the reference domain of the model to outline all boundaries and define all required parameters and variables.Figure 4 illustrates the definition of a constant temperature boundary, also known as a "Dirichlet boundary", which assumes that the soil mass around the energy pile or GSHP system acts as a large, finite heat reservoir.The temperature along certain edges or surfaces of the model domain is fixed at a predetermined soil mass temperature, ignoring time-dependent temperature variations.This boundary condition ensures the stability and accuracy of the calculations by preventing instability.The bottom and side boundaries of the soil mass domain (as shown in Figure 4, No. 1) are set to the natural background temperature of the soil type.Conversely, the surface of the energy piles (as shown in Figure 4, No. 2) represents the operating scenario of the GSHP by injecting a constant temperature similar to the temperature of the heat exchanger.The next step is to define the calculation node (as shown in Figure 4, No. 3) as a reference point within the domain, overlaying a grid array.The distance ∆r represents the distance between computation nodes, with ∆r = 0 at the center of the energy pile.The time iterations, or time step ∆t, indicate the interval between calculation iterations, which is essential for analyzing the time-dependent changes in the soil temperature.The temperature of the neighboring calculation nodes T 0 , T 1 , T 2 , T 3 , and T 4 within the mesh facilitates the process description.
overlaying a grid array.The distance ∆ represents the distance between computat nodes, with ∆ = 0 at the center of the energy pile.The time iterations, or time step indicate the interval between calculation iterations, which is essential for analyzing time-dependent changes in the soil temperature.The temperature of the neighboring culation nodes , , , , and within the mesh facilitates the process descripti The temperature rise ( ) within any domain element can be calculated using the o dimensional transient heat balance differential equation, as shown in Equation ( 1): where the thermal diffusivity is the product of the thermal conductivity (λ) divi by the soil density ( ) and specific heat capacity (c), is the time, and is the radial tance.
Break down the finite differential equation (Equation ( 1)) into Equations ( 2)-( 4) s arately and solve for the explicit form of the functions individually.The named calculat nodes to referred to in Figure 4 are used to translate the differential term into ordinary algebraic equation.The temperature rise (T) within any domain element can be calculated using the one-dimensional transient heat balance differential equation, as shown in Equation ( 1): where the thermal diffusivity (α d ) is the product of the thermal conductivity (λ) divided by the soil density (ρ) and specific heat capacity (c), t is the time, and r is the radial distance.Break down the finite differential equation (Equation ( 1)) into Equations ( 2)-( 4) separately and solve for the explicit form of the functions individually.The named calculation nodes T 1 to T 4 referred to in Figure 4 are used to translate the differential term into an ordinary algebraic equation.
∂T ∂t ∂T ∂r By substituting Equation (2) to Equation (3) into the heat balance equation, the final explicit solution for the temperature of point T 1 can be derived as given in Equation (5): where the required parameters for the calculation of the changing temperature based on Equation ( 5) are the temperature of the surrounding points (T 0 , T 2 , and T 4 ), the radius of each node from the center line (r), the changing radius (∆r), the thermal diffusivity (α d ), and the thermal time step (∆t).

Methodology
When building smart cities, the integration of energy piles into the urban infrastructure and the evaluation of the resulting heat movement on the ground is of paramount importance.These considerations are critical for optimizing the efficiency and sustainability of GSHP systems, which play a significant role in the energy framework of smart cities.In this study, the predictive power of a one-dimensional FDM was used to estimate and project the soil temperatures around an energy pile over various distances and time periods.The input parameters, which include the detailed physical and thermal characteristics of the materials and soil, were carefully collected from an extensive literature review of GSHP systems in Southeast Asia, design manuals, and reputable textbooks [12][13][14][15][16]33].The model formulation is further enhanced by an in-depth study conducted at a physical GSHP test site and an energy pile installation at Kasetsart University, Thailand [13].
The focus of our analysis will primarily be placed on two critical parameters: the apparent thermal conductivity of the soil and its density.The model assumes that the thermal conductivity increases proportionally with the temperature increase in the heat exchanger connected to the GSHP system.Variations in soil density will reflect the different geographical locations in Southeast Asia where GSHP and energy pile systems could be deployed.This approach allows for a comprehensive examination of how soil density affects temperature rise, with further details on the different soil densities and their underlying rationale to be provided in a future section.
A distinctive feature of the finite difference model is its ability to distill complex real-world heat transfer processes into a manageable mathematical model.This simplification facilitates the analysis of the vertical heat exchanger integrated into the energy pile, underscoring the model's utility in advancing the development of smart city infrastructure.

Model Definition
As urban areas in Southeast Asian regions expand, there is an increase in the construction of small to medium-sized buildings to meet the growing demand.These structures often require pile foundations to be installed in firmer soil layers to minimize settlement, especially in the soft clayey soil layer (see Table 1).The development of a one-dimensional FDM model requires simplified inputs.
The concept of energy piles and their foundations as applied in FDM modeling is shown in Figure 5.The model assumes that the energy pile is in a soft clay layer with an average thickness of 10 to 15 m below the ground surface.For the integration of a GSHP system into the pile, it is assumed that a suitable pile type is a spun pile with a diameter of 0.3 m and an opening width of 0.15 m along the depth of the pile.After the installation of the GSHP system, the pile is filled with concrete.The length of the pile is determined to be 12 m, with a working depth of 10 m for the heat exchanger, considering the soil plug at the top of the pile.The U-tube heat exchanger, which is essential to the system, has a total length of 20 m, equally divided between the inlet and the outlet.The 10 m heat exchanger's inlet, illustrated in red on the right side of the pile's figure, provides the downward flow path of the high temperature-heat medium.In contrast, the 10 m heat exchanger's outlet, illustrated in blue on the left side of the pile's figure, provides the upward flow path.Thermal input parameters for all components of the GSHP system and energy pile, as referenced in the American Society of Heating, Refrigerating and Air-Conditioning Engineers (ASHRAE) [34] and the Analysis and Design of Energy Geostructures [6], are provided in Table 2 for use in the numerical model.
Conversely, a coarser mesh reduces computational time but may sacrifice accuracy 4 provides an explanation of the mesh used in this study; each computation point i at the horizontal spacing of 0.15 m, and the spacing distance can be referred t increment r.In addition, the vertical spacing between each horizontal roller is set to accommodate the varying temperature of the heat exchanger at depth.The mesh is the grid or mesh of nodes used to compute the numerical solution of the governing equations.It divides the entire domain into small rectangular sections in 2D.The number of nodes can vary across the domain, with finer meshes used in areas where the solution needs to be finer or change rapidly.At each node, the governing equations are approximated by finite difference formulas that are replaced by references to the values of the neighboring computational nodes.The location and setup of the computing nodes has a significant impact on the accuracy and efficiency of the FDM result.
A finer mesh generally results in higher accuracy and resolution, but higher computational cost.Conversely, a coarser mesh reduces computational time but may sacrifice accuracy.Figure 4 provides an explanation of the mesh used in this study; each computation point is placed at the horizontal spacing of 0.15 m, and the spacing distance can be referred to as the increment r.In addition, the vertical spacing between each horizontal roller is set to 1.0 m to accommodate the varying temperature of the heat exchanger at depth.

Constant-Temperature Boundary Condition
To simulate the initial operating condition of the energy pile system, the temperature of a 20 m-long heat exchanger is set as a constant temperature boundary within the analysis spreadsheet.This approach simplifies the model by ignoring the radial heat convection within the heat exchanger and the fluid medium it contains.It is important to consider the different heat contributions of the medium flowing downward through the 10 m-long inlet and upward through the 10 m-long outlet of the heat exchanger and treat them as independent heat sources.
Previous studies [12,13,16,33] have monitored the average temperatures at the inlet and outlet of the heat exchanger, showing that the inlet water temperature is 37 • C, while the outlet water temperature is measured at 33 • C. Direct interpolation is used to calculate the temperature at each meter along the heat exchanger, which indicates a decrease of 0.2 • C/m along its length, as shown in Table 3.The operation of the GSHP is often tied to the demand of the building's air-conditioning system, which requires 12 h of operation to maintain a constant temperature range of 33 to 37 • C, depending on the depth of the heat exchanger.This is followed by a 12 h shutdown of the GSHP system, during which the temperature at the constant temperature boundary is adjusted to match the ambient ground temperature.Studies conducted in the Bangkok area [12][13][14], using Thermal Response Tests (TRT) or operational GSHP air-conditioning systems, have found that the background temperature of the Bangkok clay varies between 29 • C and 31 • C during different seasons.On average, this is 5 • C cooler than the ambient air temperature [15].For the FDM temperature boundary calculations, the background temperature of Bangkok clay is set to 29 • C.
To help understand the concept of a constant temperature boundary, Figure 6 illustrates the asymmetric temperature bulb surrounding the energy pile in the soil.This is based on direct interpolation of the constant temperature boundary 1 year after the simulation, using a finite element model (FEM) through PLAXIS 2D CONNECT V21.01.00.479 software [35,36], within a 3 m radius of the pile centers.
ture boundary calculations, the background temperature of Bangkok clay is set to 29 °C.
To help understand the concept of a constant temperature boundary, Figure 6 illustrates the asymmetric temperature bulb surrounding the energy pile in the soil.This is based on direct interpolation of the constant temperature boundary 1 year after the simulation, using a finite element model (FEM) through PLAXIS 2D CONNECT V21.01.00.479 software [35,36], within a 3 m radius of the pile centers.

Variable Thermal Conductivities and Soil Densities
The thermal conductivity of the GSHP system, which is integral to the performance of energy piles, is affected by various physical factors of the system components.These include the mineralogy, dry density, water content, and gradation of the soil materials as identified in previous studies [37,38].Experimental evidence from lateritic clay samples [39] indicates that thermal conductivity values for samples from the Guilin and Liuzhou areas in the Guangxi region of China are temperature-dependent.In laboratory tests, at temperatures ranging from 5 to 90 °C, clay samples at saturation levels of 4%, 44.90%, and 83.90% exhibited a linear increase in thermal conductivity, with values ranging from 0.2 to 0.4 W/m•°C, 0.6 to 1.6 W/m•°C, and 0.9 to 1.3 W/m•°C, respectively, as shown in Figure 7.
Assuming a linear relationship between temperature and thermal conductivity for clay samples, this relationship is applied in the FDM model.The thermal conductivity of soft clay from different locations is expressed as a linear function of the average inletoutlet temperature of the GSHP system.Bangkok soft clay, which has been extensively studied, serves as a representative of other soft clays.Its thermal and engineering

Variable Thermal Conductivities and Soil Densities
The thermal conductivity of the GSHP system, which is integral to the performance of energy piles, is affected by various physical factors of the system components.These include the mineralogy, dry density, water content, and gradation of the soil materials as identified in previous studies [37,38].Experimental evidence from lateritic clay samples [39] indicates that thermal conductivity values for samples from the Guilin and Liuzhou areas in the Guangxi region of China are temperature-dependent.In laboratory tests, at temperatures ranging from 5 to 90 • C, clay samples at saturation levels of 4%, 44.90%, and 83.90% exhibited a linear increase in thermal conductivity, with values ranging from 0.2 to 0.4 W/m• • C, 0.6 to 1.6 W/m• • C, and 0.9 to 1.3 W/m• • C, respectively, as shown in Figure 7.
Assuming a linear relationship between temperature and thermal conductivity for clay samples, this relationship is applied in the FDM model.The thermal conductivity of soft clay from different locations is expressed as a linear function of the average inlet-outlet temperature of the GSHP system.Bangkok soft clay, which has been extensively studied, serves as a representative of other soft clays.Its thermal and engineering properties are similar to those of other East Asian soft marine clays, such as Ariake clay, Singapore marine clay, Hong Kong marine clay, Kuantan clay, and Ho Chi Minh City marine clay.They may provide insight into the thermal performance of similar soil types.Formed during the Holocene by sediments from the Chao Phraya River delta, Bangkok clay is characterized by its high moisture content, low strength, and high compressibility, which cause it to pose geotechnical challenges.Notable studies include those by Abuel-Naga et al. [28,29], who used temperature-controlled triaxial equipment at 20 • C to 90 • C under undrained conditions and found thermal conductivities of 0.825, 1.1, and 1.3 W/m• • C. Jotisankasa et al. [13], who implemented a full-scale energy pile system in a residential setting, found thermal conductivities for steel and spun piles of 1.49, 2.01, and 3. Formed during the Holocene by sediments from the Chao Phraya River delta, Bangkok clay is characterized by its high moisture content, low strength, and high compressibility, which cause it to pose geotechnical challenges.Notable studies include those by Abuel-Naga et al. [28,29], who used temperature-controlled triaxial equipment at 20 °C to 90 °C under undrained conditions and found thermal conductivities of 0.825, 1.1, and 1.3 W/m•°C.Jotisankasa et al. [13], who implemented a full-scale energy pile system in a residential setting, found thermal conductivities for steel and spun piles of 1.49, 2.01, and 3.3 W/m•°C.Another study in Hanoi, Vietnam [40] recorded heat exchanger temperatures ranging from 37 °C to 40 °C with resulting thermal conductivities of 1.82 and 1.97 W/m•°C.The study derives variable thermal conductivities for the GSHP system from these results, which are summarized in Table 4, which provides the data for constructing four variable thermal conductivity equations.These linear equations correlate the temperature of the heat exchanger to its thermal conductivity, with cases no. 2 and no. 4 based on Jotisankasa et al. [13] and Widiatmojo et al. [40], and cases no. 1 and no. 3 offering adjusted slopes to provide intermediate ranges of thermal conductivity for the model.For clarity, all four thermal conductivity equations are shown in Figure 8, starting at a temperature of 32 °C and an initial thermal conductivity of 1.49 W/m•°C.The upper limit of the average temperature of the heat exchanger will not exceed 40 °C, according to reviews of in situ GSHP experiments.
The variation in soil density is another crucial variable that correlates with the geolocation of the GSHP system.This study considers the density of soft clays varying from 1400 kg/m 3 (representing Singapore marine clay and Ariake clay), through 1600 kg/m 3 (Bangkok clay), to 1800 kg/m 3 (Hong Kong clay).Additional variables and their effect on system performance are detailed in Table 5.The study derives variable thermal conductivities for the GSHP system from these results, which are summarized in Table 4, which provides the data for constructing four variable thermal conductivity equations.These linear equations correlate the temperature of the heat exchanger to its thermal conductivity, with cases no. 2 and no. 4 based on Jotisankasa et al. [13] and Widiatmojo et al. [40], and cases no. 1 and no.The variation in soil density is another crucial variable that correlates with the geolocation of the GSHP system.This study considers the density of soft clays varying from 1400 kg/m 3 (representing Singapore marine clay and Ariake clay), through 1600 kg/m 3 (Bangkok clay), to 1800 kg/m 3 (Hong Kong clay).Additional variables and their effect on system performance are detailed in Table 5.

Soil Name
Density (kg/m 3 ) Singapore marine clay and Ariake clay 1400 Bangkok clay 1600 Hong Kong clay 1800

Evaluation of Temperature Distribution around Energy Piles
The operation of GSHP and energy piles is crucial for the efficient performance of smart cities, with typical operating times of up to 12 h per day.This study meticulously examines the temperature dynamics surrounding the soil adjacent to energy piles, using a detailed temporal analysis ranging from short-term (initial 7 days) to long-term (annual overview) perspectives.
During the first week of observation, temperature variations around the energy pile exhibit a cyclical pattern reminiscent of a sine wave, as shown in Figure 9.This cycle is attributed to the alternating operating phases of the GSHP system.Despite the cyclic nature, a gradual upward trend in temperature is observed, indicating an accumulation of heat within the soil mass over time.This phenomenon is primarily governed by the mechanisms of heat conduction and convection, which become more pronounced with increasing distance from the core of the energy pile.Specifically, at distances of 0.15 m, 0.30 m, 0.45 m, and 0.60 m from the pile center, the temperature variations progressively decrease, indicating a decrease in the thermal gradient due to the operating cycles of the energy pile.For example, at a 0.15 m distance, the temperature variations range from 31.25 °C to 30.25 °C, showing a variation of 1 °C.This variation narrows significantly at distances of 0.3 m (0.25 °C fluctuation) and 0.45 m (0.17 °C fluctuation), illustrating thermal diffusion and distribution patterns within the soil.
In contrast, the annual perspective shows a relatively stable temperature distribution with negligible long-term increases.At the end of the year, the maximum recorded soil temperature near the energy pile (at a radius of 0.15 m) stabilizes around 31.5 °C, as shown in Figure 10.This stability underscores the equilibrium achieved between the heat input from the GSHP system and the natural heat dissipation mechanisms within the soil over an extended period.
Figure 11 further elaborates on the temperature dynamics by comparing the effects of alternating 12 h operating cycles of the GSHP system over 1 month and 1 year.Within

Evaluation of Temperature Distribution around Energy Piles
The operation of GSHP and energy piles is crucial for the efficient performance of smart cities, with typical operating times of up to 12 h per day.This study meticulously examines the temperature dynamics surrounding the soil adjacent to energy piles, using a detailed temporal analysis ranging from short-term (initial 7 days) to long-term (annual overview) perspectives.
During the first week of observation, temperature variations around the energy pile exhibit a cyclical pattern reminiscent of a sine wave, as shown in Figure 9.This cycle is attributed to the alternating operating phases of the GSHP system.Despite the cyclic nature, a gradual upward trend in temperature is observed, indicating an accumulation of heat within the soil mass over time.This phenomenon is primarily governed by the mechanisms of heat conduction and convection, which become more pronounced with increasing distance from the core of the energy pile.Specifically, at distances of 0.15 m, 0.30 m, 0.45 m, and 0.60 m from the pile center, the temperature variations progressively decrease, indicating a decrease in the thermal gradient due to the operating cycles of the energy pile.For example, at a 0.15 m distance, the temperature variations range from   In contrast, the annual perspective shows a relatively stable temperature distribution with negligible long-term increases.At the end of the year, the maximum recorded soil temperature near the energy pile (at a radius of 0.15 m) stabilizes around 31.5 • C, as shown in Figure 10.This stability underscores the equilibrium achieved between the heat input from the GSHP system and the natural heat dissipation mechanisms within the soil over an extended period.et al. [28,29] on thermally induced volume change in soft Bangkok clay.The methodological framework of this study employs FDM, recognizing its inher limitations, particularly in the spatial resolution of the temperature nodes, which placed at 0.15 m intervals from the pile centerline.This limitation ensures model accur but limits the granularity of the temperature distribution analysis.The study contrasts results of the finite difference method with potentially more accurate modeling te The methodological framework of this study employs FDM, recognizing its inherent limitations, particularly in the spatial resolution of the temperature nodes, which are placed at 0.15 m intervals from the pile centerline.This limitation ensures model accuracy but limits the granularity of the temperature distribution analysis.The study contrasts the results of the finite difference method with potentially more accurate modeling techniques, such as the FEM, to highlight areas for future study improvements.Specifically, the study highlights the potential benefits of including heat energy as a boundary condition within FEM analyses to capture a more accurate and realistic representation of GSHP system operation.
In addition, this study identifies a significant data gap regarding the temperaturedependent thermal conductivity of clayey soils, particularly for Bangkok clay, which is critical for refining the accuracy of temperature distribution models around energy piles.Future study directions include conducting laboratory experiments to determine specific thermal conductivity values for Bangkok clay and other similar soils to improve the predictability and reliability of temperature distribution models in the context of smart city infrastructure development.

Evaluation of Thermal Conductivity Variations on Temperature Distribution
Exploring the role of thermal conductivity variations within different time frames, namely 1 month, 3 months, and 1 year, provides crucial insights into the thermal dynamics around energy piles, which are an integral part of smart city infrastructures.The meticulously detailed analysis, summarized in Figure 12, provides a clear visualization of the temperature profiles around energy piles, meticulously examining the potential shifts in temperature distribution attributable to changes in the thermal conductivity of the surrounding soil medium.This analysis is critical because thermal conductivity is a fundamental property affecting heat transfer efficiency within the soil.It directly affects the operational effectiveness of GSHP systems and energy piles, which are critical components in the sustainable energy framework of smart cities.Therefore, understanding the sensitivity of temperature distribution to variations in thermal conductivity is essential for predicting the long-term performance of these systems.
Smart Cities 2024, 7, FOR PEER REVIEW 18 condition within FEM analyses to capture a more accurate and realistic representation of GSHP system operation.
In addition, this study identifies a significant data gap regarding the temperaturedependent thermal conductivity of clayey soils, particularly for Bangkok clay, which is critical for refining the accuracy of temperature distribution models around energy piles.Future study directions include conducting laboratory experiments to determine specific thermal conductivity values for Bangkok clay and other similar soils to improve the predictability and reliability of temperature distribution models in the context of smart city infrastructure development.

Evaluation of Thermal Conductivity Variations on Temperature Distribution
Exploring the role of thermal conductivity variations within different time frames, namely 1 month, 3 months, and 1 year, provides crucial insights into the thermal dynamics around energy piles, which are an integral part of smart city infrastructures.The meticulously detailed analysis, summarized in Figure 12, provides a clear visualization of the temperature profiles around energy piles, meticulously examining the potential shifts in temperature distribution attributable to changes in the thermal conductivity of the surrounding soil medium.This analysis is critical because thermal conductivity is a fundamental property affecting heat transfer efficiency within the soil.It directly affects the operational effectiveness of GSHP systems and energy piles, which are critical components in the sustainable energy framework of smart cities.Therefore, understanding the sensitivity of temperature distribution to variations in thermal conductivity is essential for predicting the long-term performance of these systems.The results of the study are quite revealing, indicating that within the observed range of thermal conductivity variations, there is remarkable stability in the temperature distribution around the energy piles over the time periods studied.This observation suggests that the thermal behavior of the system exhibits a degree of resilience to changes in soil thermal conductivity.This resilience is particularly noteworthy in the context of smart city planning and development, where energy efficiency and sustainability are paramount.Such stability in temperature distribution, despite variations in thermal conductivity, underscores a critical aspect of the thermal dynamics involved in the operation of energy piles.It suggests that the design of these systems can afford a degree of flexibility with respect to the thermal properties of the soil.This finding is invaluable to engineers and The results of the study are quite revealing, indicating that within the observed range of thermal conductivity variations, there is remarkable stability in the temperature distribution around the energy piles over the time periods studied.This observation suggests that the thermal behavior of the system exhibits a degree of resilience to changes in soil thermal conductivity.This resilience is particularly noteworthy in the context of smart city planning and development, where energy efficiency and sustainability are paramount.Such stability in temperature distribution, despite variations in thermal conductivity, underscores a critical aspect of the thermal dynamics involved in the operation of energy piles.It suggests that the design of these systems can afford a degree of flexibility with respect to the thermal properties of the soil.This finding is invaluable to engineers and urban planners, as it implies that minor variations in soil thermal conductivity may not require adjustments to the design or operation of GSHP systems.
However, it is important to contextualize these findings within the limitations of the study.The analysis, while comprehensive, is based on the range of thermal conductivity variations observed within the scope of the study.This range is inherently related to the specific soil types, moisture content, and other environmental conditions pertinent to the study area.Consequently, extending these findings to different regions or soil conditions requires careful consideration of local soil properties and environmental factors.Furthermore, the implications of this stability go beyond the immediate operational considerations of energy piles.They touch on broader issues of sustainability and resilience in urban planning.In a world where climate variability and environmental change are of increasing concern, the ability to design energy systems that are robust to changes in underlying environmental characteristics is invaluable.The findings of this study contribute to a growing body of knowledge that supports the development of more sustainable and resilient urban environments.

Impact of Clay Density Variation on Temperature Predictions
This study explores the intricate dynamics of soil density variation and its subsequent impact on the thermal behavior around energy piles, a key component in the thermal management systems of smart cities.Using FDMs, this study categorizes soft clay into three distinct density types: Singapore marine clay and Ariake clay with a density of 1400 kg/m 3 , soft Bangkok clay with a density of 1600 kg/m 3 , and Hong Kong clay with a density of 2200 kg/m 3 , as contextually implied in subsequent discussions.These categorizations are critical to understanding the nuanced impact of soil density on the thermal properties and temperature regulation capabilities of the soil surrounding energy piles.Figure 13 provides a visual representation of how variations in soil density affect temperature predictions over a one-year observation period.Notably, the graph illustrates a convergence of temperature differences over time, suggesting a diminishing effect of density variations on temperature fluctuations with long-term system operation.This phenomenon underscores the complex interplay between soil physical properties and thermal dynamics and highlights the critical role of soil density in determining the thermal efficiency of subsurface energy systems.
Figure 13.One-year predicted temperatures with different densitie Further analysis, as shown in Figures 14-16, systematically evaluates the temperature response to increasing soil density at different radial distances from the energy pile (0.3 m, 0.6 m, and 1.0 m) and over different time periods (3 months, 6 months, and 1 year).This multidimensional analysis reveals a consistent trend: as soil density increases from 1400 to 1800 kg/m 3 , a 28.57% increase, the resulting effect on soil temperature shows a decreasing trajectory.Specifically, temperature decreases of 0.02 • C, 0.01 • C, and 0.0025 • C are observed with each density increase over the specified time periods.This trend is further quantified by noting that for every 10% increase in soil density, there is a corresponding average temperature decrease of 0.007 • C at a 0.3 m radius from the energy pile.The granularity of this analysis is further expanded in Figures 14-16, where the temperature effects of increased soil density are assessed at radial distances of 0.6 m and 1 m from the energy pile.These results provide a nuanced understanding of the thermal dynamics of the soil, revealing a proportional relationship between soil density and temperature modulation.For example, at a 0.6 m radius, each 10% increase in soil density correlates with an average temperature decrease of 0.003 • C, while at a 1 m radius, the average temperature decrease is 0.0009 • C for the same density increase.
Smart Cities 2024, 7, FOR PEER REVIEW 20 corresponding average temperature decrease of 0.007 °C at a 0.3 m radius from the energy pile.The granularity of this analysis is further expanded in Figures 14-16, where the temperature effects of increased soil density are assessed at radial distances of 0.6 m and 1 m from the energy pile.These results provide a nuanced understanding of the thermal dynamics of the soil, revealing a proportional relationship between soil density and temperature modulation.For example, at a 0.6 m radius, each 10% increase in soil density correlates with an average temperature decrease of 0.003 °C, while at a 1 m radius, the average temperature decrease is 0.0009 °C for the same density increase.corresponding average temperature decrease of 0.007 °C at a 0.3 m radius from the energy pile.The granularity of this analysis is further expanded in Figures 14-16, where the temperature effects of increased soil density are assessed at radial distances of 0.6 m and 1 m from the energy pile.These results provide a nuanced understanding of the thermal dynamics of the soil, revealing a proportional relationship between soil density and temperature modulation.For example, at a 0.6 m radius, each 10% increase in soil density correlates with an average temperature decrease of 0.003 °C, while at a 1 m radius, the average temperature decrease is 0.0009 °C for the same density increase.This comprehensive study of the effects of soil density on the thermal behavior around energy piles provides valuable insights for the design and operation of geothermal systems in smart cities.It highlights the importance of the physical properties of soil, particularly density, in influencing thermal interactions within the soil mass, which in turn affects the efficiency and sustainability of energy pile systems.The results suggest that careful consideration of soil density can improve the predictive accuracy of thermal models, thereby facilitating optimized design strategies for subsurface energy solutions.In addition, the results of the study underscore the importance of incorporating detailed soil property analyses into the planning and implementation phases of geothermal projects.Understanding the specific effects of soil density variations on temperature regulation can lead to more effective thermal management practices, contributing to the overall resilience and sustainability of urban infrastructure.This study not only advances the scientific understanding of soil-thermal interactions but also provides a pragmatic guide for engineers and urban planners seeking to integrate efficient geothermal technologies into the fabric of smart cities.

Conclusions
The study on the use of GSHP systems and energy piles in smart city development has yielded several significant findings: (1) GSHP systems, especially when integrated with energy piles, have proven to be efficient and sustainable solutions for heating and cooling in urban environments, significantly reducing electricity consumption and greenhouse gas emissions.(2) The study emphasized the critical role of understanding thermal movement within the soil, particularly in soft marine clays common in Southeast Asia, for optimal implementation of GSHP systems.Soil thermal conductivity and density were identified as key factors affecting system performance.(3) Using a finite difference method (FDM), the study successfully modeled temperature distribution around energy piles and evaluated the impact of variations in soil thermal conductivity and density on system efficiency.(4) GSHP systems offer a promising solution to mitigate the urban heat island effect, especially in densely populated areas where conventional HVAC systems are less efficient.
Based on the results of this study, several avenues for future study were identified: (1) While FDM provided valuable insights, future studies could explore more sophisti- This comprehensive study of the effects of soil density on the thermal behavior around energy piles provides valuable insights for the design and operation of geothermal systems in smart cities.It highlights the importance of the physical properties of soil, particularly density, in influencing thermal interactions within the soil mass, which in turn affects the efficiency and sustainability of energy pile systems.The results suggest that careful consideration of soil density can improve the predictive accuracy of thermal models, thereby facilitating optimized design strategies for subsurface energy solutions.In addition, the results of the study underscore the importance of incorporating detailed soil property analyses into the planning and implementation phases of geothermal projects.Understanding the specific effects of soil density variations on temperature regulation can lead to more effective thermal management practices, contributing to the overall resilience and sustainability of urban infrastructure.This study not only advances the scientific understanding of soil-thermal interactions but also provides a pragmatic guide for engineers and urban planners seeking to integrate efficient geothermal technologies into the fabric of smart cities.

Conclusions
The study on the use of GSHP systems and energy piles in smart city development has yielded several significant findings: (1) GSHP systems, especially when integrated with energy piles, have proven to be efficient and sustainable solutions for heating and cooling in urban environments, significantly reducing electricity consumption and greenhouse gas emissions.(2) The study emphasized the critical role of understanding thermal movement within the soil, particularly in soft marine clays common in Southeast Asia, for optimal implementation of GSHP systems.Soil thermal conductivity and density were identified as key factors affecting system performance.(3) Using a finite difference method (FDM), the study successfully modeled temperature distribution around energy piles and evaluated the impact of variations in soil thermal conductivity and density on system efficiency.(4) GSHP systems offer a promising solution to mitigate the urban heat island effect, especially in densely populated areas where conventional HVAC systems are less efficient.
Based on the results of this study, several avenues for future study were identified: (1) While FDM provided valuable insights, future studies could explore more sophisticated modeling techniques, such as the finite element method (FEM), to achieve greater accuracy in predicting ground temperature distribution and system performance.(2) Additional study is needed to determine the temperature-dependent thermal conductivity of different types of clayey soils, particularly Bangkok clay, to refine the accuracy of temperature distribution models around energy piles.(3) Investigating the effects of rapid urbanization and the subsequent increase in the urban heat island effect on the efficiency of GSHP systems could provide valuable insights for urban planning and smart city development.(4) Future work should focus on optimizing GSHP system configurations and operating parameters to maximize the use of geothermal energy in urban development.(5) Further studies could explore the integration of GSHP systems within the broader context of sustainable and resilient urban infrastructure, considering climate variability and environmental change.
This study contributes significantly to the growing body of knowledge on sustainable urban development and highlights the potential of GSHP systems and energy piles in smart cities.By addressing the outlined future study directions, the full potential of geothermal energy in urban environments can be realized, supporting the transition to more sustainable and energy-efficient urban habitats.
This study has several limitations that should be acknowledged.First, the short period of data collection may not adequately represent the long-term performance and seasonal variations of ground source heat pump (GSHP) systems.Second, the lack of temperature data from thermal response tests (TRTs) and current GSHP installations in Southeast Asia hinders comprehensive regional studies and comparisons.Third, the one-dimensional finite difference model, while computationally efficient, simplifies thermal behavior by ignoring vertical temperature changes and possibly omitting localized effects.Finally, the study emphasizes the need to use heat energy as a boundary condition in finite element method (FEM) calculations to improve accuracy.This suggests that the existing model could be improved using more complex methods.These limitations highlight potential areas for future research to improve GSHP modeling and assessment in the region.

Figure 2 .
Figure 2. Schematic mechanism of a heat exchanger in an energy pile.

Figure 2 .
Figure 2. Schematic mechanism of a heat exchanger in an energy pile.

Figure 4 .
Figure 4. Analysis domain for FDM model with calculation node and boundary conditions.

Figure 4 .
Figure 4. Analysis domain for FDM model with calculation node and boundary conditions.

Figure 5 .
Figure 5.General definition of energy pile in FDM model.

Figure 6 .
Figure 6.Asymmetrical soil temperature around energy pile after 1 year, simulated with FEM model.

Figure 6 .
Figure 6.Asymmetrical soil temperature around energy pile after 1 year, simulated with FEM model.
3 W/m• • C. Another study in Hanoi, Vietnam [40] recorded heat exchanger temperatures ranging from 37 • C to 40 • C with resulting thermal conductivities of 1.82 and 1.97 W/m• • C.
3 offering adjusted slopes to provide intermediate ranges of thermal conductivity for the model.For clarity, all four thermal conductivity equations are shown in Figure 8, starting at a temperature of 32 • C and an initial thermal conductivity of 1.49 W/m• • C. The upper limit of the average temperature of the heat exchanger will not exceed 40 • C, according to reviews of in situ GSHP experiments.
31.25 • C to 30.25 • C, showing a variation of 1 • C.This variation narrows significantly at distances of 0.3 m (0.25 • C fluctuation) and 0.45 m (0.17 • C fluctuation), illustrating thermal diffusion and distribution patterns within the soil.Smart Cities 2024, 7, FOR PEER REVIEW 16 a radius of 1 m from the center of the pile, temperatures are observed to range from 29.5 to 32 °C, with a broader dispersion observed within a radius of 1 to 2 m, where temperatures range from 29.0 to 30 °C.These observations suggest that the effects of thermal consolidation, both primary and secondary, have minimal influence on soil volume change within the temperature range studied (25 to 90 °C), referencing the work of Abuel-Naga et al. [28,29] on thermally induced volume change in soft Bangkok clay.

Figure 9 .
Figure 9. Seven-day temperature cycle at various radii.Figure 9. Seven-day temperature cycle at various radii.

Figure 9 .
Figure 9. Seven-day temperature cycle at various radii.Figure 9. Seven-day temperature cycle at various radii.

Figure 9 .
Seven-day temperature cycle at various radii.

Figure 10 .
Figure 10.One-year temperature cycle at various radii.

Figure 10 .
Figure 10.One-year temperature cycle at various radii.

Figure 11
Figure 11 further elaborates on the temperature dynamics by comparing the effects of alternating 12 h operating cycles of the GSHP system over 1 month and 1 year.Within a radius of 1 m from the center of the pile, temperatures are observed to range from 29.5 to 32 • C, with a broader dispersion observed within a radius of 1 to 2 m, where temperatures range from 29.0 to 30 • C.These observations suggest that the effects of thermal consolidation, both primary and secondary, have minimal influence on soil volume change within the temperature range studied (25 to 90 • C), referencing the work of Abuel-Naga et al. [28,29] on thermally induced volume change in soft Bangkok clay.

Figure 9 .
Figure 9. Seven-day temperature cycle at various radii.

Figure 10 .
Figure 10.One-year temperature cycle at various radii.

Figure 12 .
Figure 12.One-year predicted temperatures with different thermal conductivities.

Figure 12 .
Figure 12.One-year predicted temperatures with different thermal conductivities.

Figure 13 .
Figure 13.One-year predicted temperatures with different densities of Bangkok clay.

Figure 14 .
Figure 14.Temperature prediction in 1 month, 6 months, and 1 year and Bangkok clay density correlation over time at 0.3 m radius.

Figure 15 .
Figure 15.Temperature prediction in 3 months, 6 months, and 1 year and Bangkok clay density correlation over time at 0.6 m radius.

Figure 14 .
Figure 14.Temperature prediction in 1 month, 6 months, and 1 year and Bangkok clay density correlation over time at 0.3 m radius.

Figure 14 .
Figure 14.Temperature prediction in 1 month, 6 months, and 1 year and Bangkok clay density correlation over time at 0.3 m radius.

Figure 15 .
Figure 15.Temperature prediction in 3 months, 6 months, and 1 year and Bangkok clay density correlation over time at 0.6 m radius.

Figure 15 .
Figure 15.Temperature prediction in 3 months, 6 months, and 1 year and Bangkok clay density correlation over time at 0.6 m radius.

Figure 16 .
Figure 16.Temperature prediction in 3 months, 6 months, and 1 year and Bangkok clay density correlation over time at 0.3 m radius.

Figure 16 .
Figure 16.Temperature prediction in 3 months, 6 months, and 1 year and Bangkok clay density correlation over time at 0.3 m radius.

Table 1 .
Basic engineering properties of clay sedimentary formations in different countries.

Table 1 .
Basic engineering properties of clay sedimentary formations in different countries.

Table 2 .
Input parameters for FDM analysis.
Figure 5.General definition of energy pile in FDM model.

Table 2 .
Input parameters for FDM analysis.

Table 3 .
The temperature of the HDPE heat exchanger at different depths of energy pile.

Table 4 .
Summary of deriving variable thermal conductivities in GSHP systems for constructing four variable thermal conductivity equations.

Table 4 .
Summary of deriving variable thermal conductivities in GSHP systems for constructing four variable thermal conductivity equations.

Table 5 .
The difference in density for clay in different countries.

Table 5 .
The difference in density for clay in different countries.