Two Phase Heat Transfer on Porous Media Base Termite Nest Structure with Finite Element Method

ABSTRACT


INTRODUCTION
Macrotermes sp.termites are one of the termite species found in Indonesia.With their tendency to build nests in the soil, including the species Macrotermes gilvus, they are known for being quite active in constructing their nests, often comprising complex structures with sufficiently wide air passages within [1] The composition of Macrotermes sp.termite nests, consisting of a mixture of soil, saliva, and other organic materials, is known to maintain a stable temperature even as the temperature outside fluctuates [2].This makes the termite nest characteristic applicable to buildings in Indonesia.One common application is the use of environmentally friendly building techniques inspired by the structure of termite nests.For example, environmentally friendly construction concepts that utilize a mixture of natural and recycled materials, as well as building designs that consider natural air circulation, can reflect the principles found in the structure of Macrotermes sp.termite nests [3].By applying these principles, it is hoped that buildings in Indonesia can become more energy-efficient and environmentally friendly, while benefiting from the uniqueness and sustainability of the local environment.
Thermal insulation is an essential aspect of constructing buildings that can affect their energy efficiency and comfort levels.Termite nests serve as a great inspiration for thermal insulation structures due to their shaft structures that maintain a stable temperature inside the building, both during the day and night [2].However, further research is required to optimize the use of the termite nest model for thermal insulation in buildings [4].
Porous media is a commonly occurring structure in nature that consists of solid objects with many empty spaces.The Diffusion Equation is not only limited to fluid flow but also heat transfer.It is used for modeling fluid flow in porous media [5].The heat transfer mechanism is of particular interest in porous media, and termite nests serve as the primary inspiration for this heat flow model.According to research from Harvard, the primary heat flow process in termite nests is convection, which maintains temperature stability.The heat transfer model is a partial differential equation [4].Partial differential equations can be solved in two ways: analytical solutions and numerical solutions.Numerical solutions are widely discussed through various methods, such as the Finite Element Method, Finite Difference Method, and Finite Volume Method [6][7][8].The finite element method is widely used in solving solutions of hot equations [9][10][11].
The thermal behavior of termite nests cannot be accurately predicted with traditional approaches such as thermal conduction models.Therefore, further research is needed to develop mathematical models that can predict the thermal behavior of termite nests.A study by Lu et al. [12] in 2024 has developed a porous media-based mathematical model to simulate fluid flow and heat transfer.
The integration of energy conservation concepts and fluid flow equations through porous media in porous media mathematical models will provide precise and accurate solutions for modeling the thermal behavior of termite nests as thermal insulation.A study from Feng et al. [13] has demonstrated that porous media mathematical models can provide more accurate results in predicting the thermal behavior of porous insulation materials compared to thermal conduction models.
Research on the development of thermal insulation models in termite nests using a porous media mathematical approach needs to be conducted.This will help us understand how heat is transferred in such structures.The model is expected to provide significant benefits for the development of thermal insulation technology in buildings.This will result in buildings that are more efficient and environmentally friendly.Recent research published in the journal Science Advances revealed that termite nests consist of many interconnected microscopic pores [14].These pores enable termite nests to regulate thermal, humidity, and internal carbon dioxide gas concentrations.This allows termites to move around without any interference.Essentially, termites work like lungs, inhaling and exhaling air automatically.The microscopic structure of building walls can optimize thermoregulation and ventilation.This latest finding could also lead to the development of building designs that have thermal regulation systems and automatic air circulation.
Darcy's Law is the basic law of fluid motion in porous media Darcy's Law defines the flow rate of water through a porous medium, assuming laminar flow.By combining the concepts of permeability and pressure in determining the flow rate of fluid, the Darcy equation is obtained as follows [15,16]: where, Q=flow discharge (m 3 /s); A=the area of the crosssectional area passed by the fluid (m 2 ); k=hydraulic conductivity (m/s); L=the length of cross-sectional distance (m); ∆ℎ = ℎ  − ℎ  changes in water level or hydraulic head (m).
A negative sign indicates that the flow is a combination of increasing gradients, in other words the flow will occur from the high end to the low end.
Darcy's Law defines permeability as the measure of a porous material's ability to allow fluid to flow through it.In other words, the rate of flow of a homogeneous fluid through a porous material is directly proportional to the pressure gradient and the normal direction to the cross-sectional area of the flow direction and inversely proportional to viscosity.The value of permeability depends on the size of the pores and their distribution within the porous material [3,17,18].
The part of the pore space that is filled with fluid is referred to as fluid saturation or fluid saturation.Saturation is defined as the ratio between the pore volume of a porous medium occupied by a particular fluid and the total pore volume of a porous medium [14].

𝑆 = 𝑉 𝑓 𝑉 𝑡
where,  = Saturation;   = Pore volume containing fluid;   =Total interconnected pore volume.By S is saturation, Vf the volume of the pore containing the fluid, and Vt the total volume of the corresponding pore.The mechanism of fluid movement in a porous medium is determined by fluid transport.The diffuse fluid in porous media is caused by dispersion.In practice, dispersions are used in the fluid flowing state.A measure of the ability of a porous medium to spread fluid in different directions is called the dispersion coefficient [19].

𝐽 = −𝐷 ( 𝑑𝐶 𝑑𝑥 )
J is the fluid transfer rate, D is the disperse coefficient, is the concentration gradient, and-means the direction of diffusion follows the direction of the concentration gradient.The amount of fluid mass in a porous medium will forever remain unchanged.This is following the law of conservation of mass which states that mass can change its shape, but cannot be made or destroyed.This makes fluid movement in porous media easy to predict for a certain time [20].
The numerical method is a solution technique that is formulated mathematically utilizing calculation/arithmetic operations and is done repeatedly with the help of a computer or manually.The Finite Element Method is a numerical method used to solve problems of a continuous system, where in solving the problem, the system was first idealized as an arrangement of discrete elements that make up the distribution of a voltage and displacement.So that the complete solution is produced from the combination of stress and displacement distributions that meet the equation of equilibrium of force and suitability of displacement at the connection of each element.Finite element methods are used in finding numerical solutions of partial differential equations, as in porous media.Based on the laws used, the exact research parameters to obtain mathematical modeling use parameters such as permeability, porosity, compressibility, viscosity, and thermal conductivity of termite nest objects.

Data
The research utilizes heat transfer information obtained through observations.The observations were carried out on a termite mound belonging to Captotermes curviganathus.Data collection took place in Boyolali regency, Indonesia, between September 11th and 26th, 2023.A total of 11 data points for location and 24 data points for time were gathered for this study.

Methods
Data is collected with thermocouple in from outside termite nest into inside termite nest.The obtained data is then used for determining the boundary conditions of the equation.Determination of the boundary conditions employs polynomial interpolation.The obtained data is used for validating parameters.

MATHEMATICAL MODELING
In termite nests, thermals propagate into the nest by going through the soil first.Thermals propagate on the ground which is the solid part, so that thermals propagate conductively.In this section thermal insulation occurs.After going through the conduction phase, thermals propagate on termite nests which are part of the void.Thermal propagation occurs by convection in termite nests as shown in Figure 1.

Figure 1. Termite nest structure
In this phase, the thermals are no longer isolated so that the temperature in the termite nest tends to be stable.This shows if the thermal transfer in termite nests is multiphase, namely from the conduction phase to the convection phase.In this research, we limit mathematical modeling into only 1dimensional conditions, where heat flows through the soil surface in the porous media of termite nests.
The isolation model is constructed by several laws of physics, the conduction model is built.The law of conservation of energy and the law of Fourier.In the conduction phase, thermal uses Fourier's law because there are conduction variables.Fourier's law can be used to explain how thermal is transferred among the parts of a termite nest.This allows the termite nest to maintain a suitable temperature for their activities inside the nest.First of all, we use the law of energy conversation to illustrate heat transfer.
The Law of Energy Conservation is as follows: Heat changes with time, so that the thermal and heat parts change into derivatives of time as follows: Fourier's law is formulated as follows: The total heat of a system is a combination of exhaustive heat and incoming heat.
Press ( 3) is substituted in the section of Eq. ( 4) so that the following equation is obtained: The heat conduction equation in porous media is as follows: The equation is then written in derived form as follows: E is energy,  is mass, c is specific thermal capacity, u is thermal, ρ is density, μ is viscosity, A is cross-sectional area, and u is discharge.
Convection constructed with Darcy's Law and Continuity Equation, Darcy's Law in press (1), can be rewritten as u.
In the above equation, is the density of the fluid,  is porosity, and   is the source of mass.
In the above equation,   and   are the compressibility of the fluid and the compressibility of the rock, respectively.  refers to total compressibility.Substitute Eq. ( 11) becomes Eq. ( 9): In the equation above,  0 is the initial permeability and  0 is the initial porosity.
According to Gay-Lussac's law: where,  is thermal,  is a constant, or it can also be written  ∝ .

𝑃 = 𝑟𝑢
So the heat equation in a porous medium is: where, (, ) is the heat of the solution while r is the ratio of pressure and thermal (constant).
So that the equation of 1-dimensional heat convection in porous media is: Or in general it can be written as: In the conduction phase using boundary conditions with boundary conditions, then the convection part of the boundary conditions is:

NUMERICAL SIMULATION
The system employs 11 element partitions in the x position section, employing a linear partitioning approach.Time discretization utilizes the Runge-Kutta method with a 1-hour time unit for single day.The simulation results are as follows.
Point 0 is the time that shows t=6 (06:00).The simulation shows the highest thermal of termite nests at 32℃, the lowest thermal at 20℃, and thermal stable at 26℃.The solution produces a multiphase heat equation, that is, phase conduction and convection phase.The thermal is isolated in the conduction phase to a position x<35cm.Therefore, the thermal is stable in the convection phase when x>35cm.In phase convection, its thermal conductivity is minimal, hindering the propagation of heat.As a result, the thermal is always stable at 26℃.
Next, a thermal comparison between the simulation outcomes and real data is conducted (Figure 2).We depict 24 time points per hour derived from measurements over the course of a single day.
Based on Figure 3, it is evident that the real data closely aligns with the simulation results, indicating the predictive capability of the simulation data.During the conduction phase, thermal levels exhibit an ascent at t=1 (01:00) and t=19 (19:00), while experiencing a decline at t=7 (07:00) and t=13 (13:00).In contrast, during the convection phase, thermals maintain a consistent movement throughout each time unit.This observation suggests that during the conduction phase, thermal insulation is maximized, resulting in the cessation of heat propagation during the convection phase.
After obtaining a thermal insulation model, a model evaluation was carried out using MAPE and R 2 Score.From the Table 1 obtained the value of MAPE in the conduction phase of 3% with its R 2 Score is 0.74.In the convection phase, the MAPE value is 0.2% and its R 2 Score is 0.85.Based on MAPE criteria and research methods, the results of model evaluation show good preparation.This is evidenced by the MAPE <20% and 0.7≤R 2 Score<0.9.
Errors obtained in model evaluation are required to validate the model.Good isolation model is indicated by a small error value.Furthermore, such models can be used for prediction and analysis of thermal distributions.In the model there are parameters that affect it.Therefore, parameter analysis is needed to see the distribution of the model.
The parameter k in this conduction section is k1=kρc.Four parameter schemes were used in this study are, k1=1, k2=0.494,k3=0.0055, and k4=0.00275.Optimal k is a parameter obtained from field data.In the conduction phase, a given k shows the resulting thermal difference at each point (position).This means that parameters affect the conduction phase.While in the convection phase, the thermal remains stable at every point (position).This shows that the parameter has no effect on the convection phase.
The following thermal patterns in the parameter analysis during the conduction phase.
From Figure 4, the given k values result in different positions.The parameters k4 are isolated at point (position) 10cm, k3 is isolated at point 1cm, k is optimal at point 20cm, k2 and k1 are squeezed at point 35cm.So, the position length at k4<k3<k is optimal <k1.This shows that the smaller the k value, the closer it is to 0 and the smaller the length (Figure 5).
In the convection phase there are combined parameters of . When these parameters are raised to the maximum or lowered to an infinite minimum, the thermal remains stable.
In the Figure 4, 12 parameters k are given in the convection phase.The value 0.17 ≤  ≤ 1.00, obtained thermally that is always stable.Experiments when given k close to 0, i.e. k1 and k2, show unstable thermals.This can be affected by its thermal conductivity.Thus, broadly speaking, the given parameters k does not have a significant influence on the convection phase.As a result, the thermals are always stable.
Based on Figure 4, the thermal will be more quickly isolated in the conduction section if the parameter is smaller, in the graph of the parameter value below 0.011, it is quite difficult to find a material with that much thermal conductivity, so it can be modified by reducing the density or compressibility of the material.Materials with a conductivity of 0.011 can maintain thermal by carrying out heat storage.The thermal will be isolated more slowly if the total parameter is close to 1 or reaches the limit of the insulator material, because the thermal has been isolated in the conduction section, the size of the volume in the convection section is also affected in thermal stability.

CONCLUSION
In termite nests, there is a change in heat propagation from conduction to convection, so that heat undergoes 2 phases in heat propagation.In the conduction phase, the thermal magnitude is still influenced by parameters in the termite nest.While in the convection phase, the thermal is always stable and not too influenced by any parameters.Based on the results of this research, it can be understood how changes in parameters affect thermal distribution.Permeability, porosity, compressibility, viscosity, and thermal conductivity can serve as a basis for the construction of building compositions with natural thermal insulation.The conducted simulation demonstrates that the constructed model adeptly captures the heat distribution within termite nests during both the Conduction and Convection Phases, achieving an MAPE (Mean Absolute Percentage Error) of 0.0349 with an R 2 Score of 0.742 for the Conduction Phase, and an MAPE of 0.0023 with an R 2 Score of 0.85 for the Convection Phase.

Figure 3 .
Figure 3.Comparison of simulation results and real data

Figure 5 .
Figure 5. Temperature isolation position Substitution of the equation of Darcy's Law by velocity in the continuity equation yields:

Table 1 .
Evaluation score