Experimental Study on the Thermal Conductivity of Improved Graphite Composite Insulation Boards

: The thermal conductivity of thermal insulation materials directly affects the building energy consumption. The types and constituents of thermal insulation materials in thermal insulation boards are the key to determining the insulation performance. By optimizing the material constituents and ratios, this paper proposes an improved graphite composite insulation board (GCIB), which has lower thermal conductivity and good ﬁre resistance. Through theoretical derivation, it is found that the limit range of the thermal conductivity of the new GCIB is 0.042–0.064 W/(m · K). Combined with the results of theoretical value analysis, and according to the ratios of material components, the random distribution function of each material component is constructed, and the numerical model of GCIB is established. Through numerical analysis, the range of thermal conductivity of the new composite insulation board is 0.046–0.050 W/(m · K). Finally, we establish an experimental model of the new GCIB. Through the model test of six GCIBs, the thermal conductivity of the new GCIB is obtained as 0.046 W/(m · K), which is in good agreement with the results of theoretical analysis and numerical simulation. Through theoretical analysis, numerical simulation and a sample test, this paper veriﬁes the better thermal insulation performance of the improved GCIB, providing theoretical and numerical simulation methods for the new GCIB, as well as a theoretical reference for the promotion and application of the GCIB.


Introduction
In recent years, the global energy expenditure in industrial and residential construction has become one of the most important concerns [1]; due to industrial development and population growth, the amount of housing, industry, commerce and urban construction will increase significantly.It is estimated that the global energy consumption will increase by 64% by 2040 [2].Therefore, in order to achieve the sustainable development of the country and society, it is necessary to use thermal insulation materials to better save energy and strengthen the sustainable energy strategy of the construction sector.The construction industry improves its energy efficiency by continuously developing new thermal insulation materials.The thermal insulation efficiency of thermal insulation materials mainly depends on the thermal conductivity and their ability to maintain their thermal characteristics for a period of time.Thermal conductivity is one of the main characteristics of thermal insulation materials in the building industry.The most accurate method to obtain the thermal conductivity of composite insulation board is to measure it according to the standard test method [3].The thermal insulation of the building envelope is very important for energy conservation and a comfortable indoor environment.For the envelope of a building structure, the lower the thermal conductivity is, the better the thermal insulation performance of the materials is, and the higher the energy conservation rate of the building is.However, the thermal insulation materials are mostly flammable materials, and the improvement of the thermal insulation performance will reduce the fire resistance of the building structure.Therefore, on the premise of ensuring the fire protection requirements of building structure maintenance component materials, it is crucial to improve the thermal insulation performance of building structures by adjusting and improving the thermal conductivity of thermal insulation components in maintenance components.
As energy becomes increasingly scarce, thermal insulation materials are being forcibly used in buildings.Thermal insulation materials are a type of material or a combination of materials, usually composed of solid matrix materials and gas materials, and gas materials are randomly or regularly scattered in cells, pores or gaps [4,5].Thermal insulation can be applied to certain structural components, including walls, roofs, ceilings, windows and floors [6].Abu Jdayil et al. [7] reviewed the different types, manufacturing methods and characteristics of the traditional and most advanced thermal insulation materials in recent decades.Most of the available thermal insulation materials can be divided into four categories, including inorganic materials, organic materials, composite materials and advanced materials [8].The consumption of inorganic materials accounts for 60% of the market, while that of organic insulating materials accounts for 27% [1].According to the selection of thermal insulation materials for the exterior walls of public buildings in Jiangxi Province in 2019, the use of inorganic thermal insulation mortar accounted for 52.8%, rock wool board accounted for 28.0%, foamed cement accounted for 9.4%, inorganic nano silicon thermal insulation board accounted for 4.9%, vitrified micro bead board accounted for 4.6%, reflective thermal insulation coating for 0.1% and foam glass thermal insulation board for 0.2%, and most of the thermal insulation materials were inorganic materials.Lakatos et al. [9] used different methods to measure the thermal insulation performance of vacuum insulation board with an EPS protective layer, and they performed a comparative study.Berardi [10] studied the influence of environmental factors on the equivalent thermal conductivity of several foam materials.Kumar et al. [11] studied the performance of a variety of building insulation materials and gave a comparison of their performance in different climatic regions.Yang et al. [12] conducted research on the thermal conductivity of aerogel thermal insulation board for building energy conservation through numerical simulation and tests.König et al. [13] quantified the contribution of the effective thermal conductivity of open cell foam glass.Lamy-Mendes et al. [14] introduced the remarkable thermal conductivity of silica aerogel in the application of building thermal insulation.Li et al. [15] carried out a series of experimental studies on the thermal performance of vacuum insulation board and extruded polystyrene foam building walls.Wang et al. [16] studied the change in thermal conductivity of eight common building insulation materials (i.e., glass wool, rock wool, silicone blanket, foamed polystyrene, extruded polystyrene, phenolic foam, foam ceramics and foam glass) with temperature and relative humidity through experiments.Under different external conditions, many scholars have carried out a series of studies on the thermal conductivity of thermal insulation materials [6,7,[17][18][19].Petrosyan [20] studied the influence of the presence of thermal insulation in the structure on the cold load required for cooling, and revealed a pattern of cost changes in the case of insulating materials.Khoukhi [21] elucidated the combined impact of heat and humidity transfer on the thermal conductivity of polystyrene in building insulation materials.Berardi et al. [22] obtained the effects of accelerated aging processes in laboratory conditions over the thermal performance of aerogel-enhanced insulating materials.Berardi et al. [23] presented an experimental analysis of the thermal conductivity of four materials, namely rock wool, fiberglass, extruded polystyrene and polyisocyanurate, and explained how assumptions about thermal conductivity in simulations affect performance estimates.Al-Homoud [24] presented an overview of the performance characteristics and the main features of common building thermal insulating materials and their applications.Berardi et al. [25] took common insulating materials as examples (e.g., fiberglass, rock wool, polyisocyanurate and extruded polystyrene), and quantified the impact of the temperature dependency of the thermal conductivity in exterior walls and flat roofs.Wu et al. [26] investigated the thermal conductivity of polyurethane (PU) foams theoretically and experimentally.Hoseini et al. [27] presented a theoretical and experimental study on the effective thermal conductivity of aerogel composites.Majumder et al. [28] provided a detailed analysis for the thermal characterization of recycled materials for building insulation.While thermal insulation boards are mostly composed of a variety of materials, the ratios of thermal insulation materials will directly affect the thermal conductivity [1] and fire resistance [29] of thermal insulation boards.The ratios of each material component in building thermal insulation boards should comprehensively consider the thermal conductivity, fire resistance and other factors, which is relatively lacking in research.The GCIB uses highly flame-retardant graphite composite polystyrene particles as a filler, high-strength inorganic polymer materials as a binder and a proper amount of functional additives.It is processed by mixing, equipment pressing, curing, cutting and other manufacturing processes.It has good fire resistance (combustion performance can reach Grade A), low thermal conductivity, a good thermal insulation effect, easy construction and a long service life, good security and other advantages, and the GCIB will not reduce its combustion performance and mechanical properties when reducing its density and thermal conductivity [30].
By the end of 2019, China's total urban and rural buildings had reached 62 billion m 2 , and building energy consumption accounted for around 22% of the total social energy consumption.Improving building energy efficiency has become one of the most important means to achieve China's carbon peak and neutrality goals in the construction field [31][32][33].This paper takes the improved GCIB as the research object.In Section 3, according to the ratios of each material component in the insulation board, the series and parallel theoretical analysis models of the heat conduction theory of the graphite composite insulation board are proposed, and the thermal conductivity limit of the graphite composite insulation board is obtained by using the Fourier heat conduction theory.In Section 4, according to the ratios of each material component, we establish a numerical model of the GCIB with random material distribution.Through numerical simulation, we obtain the average heat flow density distribution relationship of the GCIB, and calculate the thermal conductivity of the GCIB.In Section 5, we configure samples of the GCIB corresponding to the material ratios, and further verify the correctness of the theoretical analysis and numerical simulation through heat conduction tests.This paper studies the thermal conductivity of graphite composite insulation board from three aspects of theoretical analysis, numerical simulation and model tests,; puts forward corresponding theoretical and numerical simulation calculations and analysis methods; verifies the excellent thermal insulation performance of the graphite composite insulation board, and provides a theoretical reference for the promotion and application of graphite composite insulation board.

Components of the GCIB
The interior of the GCIB is mainly cemented by aggregates and cementitious materials, mainly consisting of cement, vitrified microspheres, graphite polystyrene particles and silica fume.In order to meet the requirements of thermal conductivity and fire protection in the building code, a DRCD3030 intelligent thermal conductivity tester, electronic balance and other major equipment are used to measure the thermal conductivity according to the lightweight aggregate standard [34] and the guarded hot plate method [35]; the component ratios and thermal conductivity of each material of the improved GCIB are shown in Table 1.Among them, vitrified microspheres and graphite polystyrene particles are the main thermal insulation materials.Cement, as the main cementing material, forms a slurry after adding water to firmly cement the vitrified microspheres and graphite polystyrene particles together.Silica fume, as an additive, can improve the tensile strength, compressive strength, wear resistance, corrosion resistance and flame retardancy of GCIBs, and also reduce the cost of products.(1) Cement The test mainly adopts P • O 42.5 ordinary Portland cement [36], with a bulk density of 1200 kg/m 3 , a specific surface area of 320 m 2 /kg and thermal conductivity of around 0.5 W/(m • K).Its chemical components are shown in Table 2. (2) Vitrified microspheres A vitrified microsphere is a type of lightweight filler aggregate and thermal insulation material, which has very stable physical and chemical properties, strong anti-aging and weather resistance properties and excellent thermal insulation, fire prevention and sound absorption properties.As a type of green and environmentally friendly building thermal insulation material, the vitrified microsphere has a porous internal structure, a glassy closed external surface and a fully spherical body, with good thermal insulation performance.Compared with the traditional expanded perlite, the vitrified microsphere has weak water absorption capacity and good aging resistance [38].
As one of the main aggregates of the GCIB, the main chemical components are SiO 2 , Al 2 O 3 , CaO.There are different types of vitrified microspheres, which are granular and different in appearance.The size range of vitrified microspheres is 0.1-2 mm, which can be divided into 20-30 mesh, 30-50 mesh, 50-70 mesh, 70-90 mesh and 90-120 mesh according to the size; the appearance of 20-30 mesh vitrified microspheres is as shown in Figure 1.For different sizes of vitrified microspheres, the bulk density is 50-100 kg/m 3 , the thermal conductivity is 0.028-0.048W/m • K, the floating rate is greater than 95%, the vitrification rate of the ball is greater than 95%, the water absorption rate is less than 50% and the melting temperature is 1200 °C.The density of vitrified microspheres is far lower than that of mortar, and the thermal conductivity is significantly lower than that of mortar.
(3) Graphite polystyrene particles As shown in Figure 2, graphene polyphenylene particles are organic materials with good sphericity, low water absorption, low thermal conductivity (generally less than 0.041 W/(m • K)), excellent thermal insulation performance and a packing density of approximately 1/10 that of the vitrified microsphere [39].However, they have poor particle grading, a large void ratio, low cementitious strength with cement and poor hightemperature resistance and aging resistance.
Due to the poor grading of graphite polystyrene particles, the double mixed aggregate technology can be used to mix vitrified microspheres and polystyrene particles in an appropriate proportion, which can not only overcome defects such as the high water absorption of the vitrified microsphere and the lower thermal insulation performance of the product in the later period, but also improve the aggregate grading, reduce the cement consumption, meet the dry density requirements and achieve a good thermal insulation effect.(4) Silica fume In the process of smelting industrial silicon and ferrosilicon in the industrial electric furnace at a high temperature, silica fume (Figure 3) can be obtained through the collection and treatment of waste gas and soot by the capture device [40].SiO 2 content accounts for around 90% of the total smoke and dust in the escaping smoke and dust, so the main component of silica fume is SiO 2 .
Silica fume is not easy to react with other substances, and does not react chemically with most acids and alkalis.Its particles are evenly covered on the surface of the object, with strong corrosion resistance.The particle grading of silica fume is reasonable, which can reduce and eliminate sedimentation and delamination, improve the grading of vitri-fied microspheres and graphite polystyrene particles and improve the strength, friction resistance, corrosion resistance and other properties of GCIBs.0 10 20 mm

Methods
In this paper, theoretical analysis, numerical simulation and sample tests are used to study the thermal conductivity of GCIBs.Through comparative analysis, the correctness of the test results is verified, and the theoretical and numerical simulation methods that can be used for the thermal conductivity of GCIBs are given.The research technology route is shown in Figure 4. First, according to the material components of the GCIB, a theoretical analysis model is established.Through series and parallel models, the theoretical calculation expression of the upper and lower limits of the thermal conductivity of the GCIB is obtained.Secondly, the numerical calculation model of the GCIB is established, and the simulation results of the thermal conductivity of the improved GCIB are obtained by using ANSYS for numerical analysis, and the influence law of material components on thermal conductivity is further studied.Next, the sample of the improved GCIB is obtained through the steps of weighing -> mixing -> molding -> curing -> molding.The thermal conductivity of the improved GCIB is measured by the steady-state method.Finally, by comparing the results of the theoretical analysis, numerical simulation and experiment, the final conclusion is drawn.

Theoretical Analysis of Heat Conduction in the GCIB
For thermal insulation materials with different material components and structural compositions, the calculation methods of thermal conductivity are different.Song et al. [41] gave a summary of the fundamentals, constituents, constructions and performance of the vacuum insulation panel (VIP).The Simpson [42] model considers the influence of the time effect and obtains the theoretical calculation expression of expanded polystyrene (EPS) foam boards at a given time.
where λ is the thermal conductivity at a given time, the subscript i denotes the initial thermal conductivity and t is the time in days.∆λ represents the incremental increase in thermal conductivity from the initial value to the final plateau value.C is a time constant that determines the aging rate.
The Wei [43] model studies the thermal conductivity of silica aerogel and proposes a unit cell model for the conductive heat transfer of the xonotlite-aerogel composite insulation material, by simplification; xonotlite-type calcium silicate is considered as a periodical array of hollow cubic structures with connecting bars, and the final heat transfer expression can be written as where γ a = a/l, γ b = 2h/a and γ c = c/a.a, l, c and h are the structural parameters in the unit cell; β 1 = 1 − k g /k cq , k g and k cq are the effective thermal conductivities of gas in xonotlite-type calcium silicate and the shell of xonotlite-type calcium silicate, respectively.ψ is the filling coefficient of the aerogel, and the expression of k ae is where a 1 = d/D, a 2 = a/d.k g is the gaseous thermal conductivity in the aerogel, which is apparently different from that of the gas in free space.
where f 1 and f 2 are the stem density of fibers and the material density; λ s,g is the the described equation for serial systems and D is the shape-related parameter.
A large number of studies have been conducted on various types of thermal insulation materials.Based on these studies, this section mainly uses series and parallel models to study the thermal conductivity of GCIBs.
In combination with the comprehensive requirements of thermal insulation materials for thermal conductivity and fire resistance, the improved GCIB is composed of four materials.Assuming that the density, mass percentage and thermal conductivity of the various materials constituting the insulation board are known, the volume ratio of each corresponding material of the insulation board is where Ω i is the volume ratio of the ith material in the GCIB, ω i and ω k are the mass percentages of the ith constituent material, ρ i is the density of the ith constituent material (unit: kg/m 3 ), and i is the various materials that make up the GCIB.The local Fourier heat conduction expression of the four materials constituting the GCIB is [48-50] where Q i is the heat transferred from the inside and outside of the ith material in the GCIB (unit: W), λ i is the thermal conductivity of the ith material (unit: W/(m • K)), A i is the heat transfer area of the ith material (unit: m 2 ), ∆T i is the temperature difference between the inside and outside surfaces of the ith material (unit: K), and ∆d i is the thickness of the ith material (unit: m).The macro Fourier heat conduction expression of the GCIB is [48-50] where Q is the heat transferred from the inside and outside of the GCIB (unit: W), λ is the equivalent thermal conductivity of the insulation board (unit: W/(m • K)), A is the heat transfer area of the insulation board (unit: m 2 ), ∆T is the temperature difference between the internal and external surfaces of the insulation board (unit: K), and ∆d is the thickness of the insulation board (unit: m).Since the GCIB is composed of four component materials, in order to obtain the upper and lower limit values of the equivalent thermal conductivity, we assume that the combination of the four materials has two theoretical models, parallel and series, as shown in Figure 5, in which T h and T c represent the temperature loads on both sides of the insulation board, respectively.
(1) Upper limit of thermal conductivity in parallel model of insulation board In the insulation board, the upper bound λ ub of thermal conductivity in the parallel model is the volume average of the constituent materials' thermal conductivity.According to the ratios of each component of the GCIB, the parallel analysis model of the insulation board is established, which meets the basic assumptions.
By ( 9) and ( 10), we can obtain the volume ratio of various materials in the insulation board: Using the relationship in (10), and combined with ( 7)- (11), there are or Substituting ( 11) into (13), we obtain the upper bound equation for calculating the thermal conductivity of the insulation board: where λ ub is the upper limit value of the thermal conductivity of the GCIB.
(2) Lower limit of thermal conductivity in series model of insulation board In the insulation board, in order to obtain the lower bound λ lb of the thermal conductivity coefficient in the series model, we first obtain the reciprocal of the constituent materials' thermal conductivity, take the volume average and then calculate the reciprocal.Similarly, the series analysis model of the GCIB is established to meet the basic assumptions By ( 15) and ( 16), we can obtain the volume ratio of various components of the insulation board: Combining ( 7), ( 8) and ( 15)-( 17), we obtain By (18), we obtain i.e., ∆d i By ( 17) and (20), we obtain Substituting ( 16) into (22), we obtain the lower bound equation for calculating the thermal conductivity of the GCIB, i.e., where λ lb is the lower limit value of the thermal conductivity in the GCIB (unit: W/(m • K)).
(3) Equivalent thermal conductivity of the GCIB Assuming that the particles in the GCIB are evenly distributed in all directions, the true value of the equivalent thermal conductivity in the insulation board should be close to the average of its upper and lower limits, hence

Numerical Simulation Analysis of Heat Conduction of the GCIB 4.1. Finite Element Model
The GCIB is composed of four materials randomly, and its interior is relatively complex.In order to facilitate numerical simulation, we propose the following assumptions:

•
The GCIB is continuous inside, and the material is evenly distributed and dense, without cracks and gaps.

•
The four materials inside the insulation board are closely bonded, and the materials are in a binding state during the analysis.

•
The thermal conductivity of each constituent material is constant and independent of size and temperature [51].
In order to conduct the numerical simulation on the GCIB, a representative thermal conductivity unit is constructed as shown in Figure 6, and its geometric structure is set as 1 unit [51,52]-that is, the side length of the planar thermal conductivity unit L = 1.
The thermal conductivity of the GCIB is analyzed numerically by using the higher steady state method [51][52][53].Planar thermal unit PLANE55 is selected to simulate various materials in the GCIB.The density, mass ratio and thermal conductivity of each component material are shown in Table 1.According to the assumption, since the four component materials are uniformly distributed and dense in the GCIB, the randomly distributed unit numbering program is compiled using Python to establish the calculation and analysis model of the four randomly distributed and dense materials, as shown in Figure 6.According to Figure 6, the finite element model of the GCIB is established by using ANSYS, as shown in Figure 7, which includes the four component materials shown in Table 1, where L = 1 is the size of the heat conduction unit.Since the steady state method is used for numerical simulation, the upper and lower boundary conditions of the finite element model are set to the ideal thermal insulation state, and the left and right ends are, respectively, constant temperature boundaries (T h = 50 K and T c = −50 K) along the X axis.We construct different temperature loads on the left and right sides of the finite element model, respectively, and set the upper and lower sides of the Y axis as the insulation boundary state.
In the process of steady state heat conduction analysis, the heat exchange caused by convection and radiation (i.e., ideal heat conduction state) is not considered, and the heat conduction in the plane meets where T is the temperature (unit: K), q is the heat flux vector (unit: W/m 2 ) and λ is the thermal conductivity (unit: W/(m • K)).
By ( 26)-( 29), we can obtain the heat flow density distribution and temperature field of the thermal conductivity calculation unit of the GCIB, further calculate the average heat flow density of the thermal conductivity unit and substitute it into Fourier equation ( 30); thus, we can obtain the thermal conductivity coefficient of the thermal conductivity unit in the insulation board.
where λ is the equivalent thermal conductivity of the thermal conductivity calculation unit in the GCIB (unit: W/(m • K)), ∆T is the temperature difference in the direction of heat conduction (i.e., T h − T c , unit: K) and q avg is the average heat flow density of the thermal conductivity calculation unit of the GCIB (unit: W/m 2 ).

Numerical Simulation Results and Analysis
(1) Thermal conductivity of the improved GCIB Through the steady state analysis of heat flow, Figure 8 shows the analysis results of the heat flow density of the thermal conductivity calculation unit in the improved GCIB.Through parametric programming, we extract and analyze the average heat flow density of the thermal conductivity calculation unit as q avg400 = 4.603/(W/m 2 ), q avg10,000 = 4.593/(W/m 2 ).
According to the calculation results of the average heat flow density obtained from the numerical simulation, using (30), we obtain the thermal conductivity of the two models as λ 400 ≈ λ 10,000 = 0.046 W/(m • K), which is within the range of the theoretical calculation results.

Influence of Material Component Ratios
In order to explore the influence of the volume ratio of graphite polystyrene particles on the thermal conductivity in the GCIB, and keep Ω 1 and Ω 2 unchanged (i.e., Ω 1 = 5%, Ω 2 = 25%), through numerical simulation, we obtain the influence of the change in the Ω 3 and Ω 4 ratio on the thermal conductivity of the insulation board.Similarly, we keep Ω 2 and Ω 4 unchanged (i.e., Ω 2 = 25%, Ω 4 = 2%), and we obtain the effect of the change in the ratios of Ω 3 and Ω 1 on the thermal conductivity of the insulation board.We keep Ω 1 and Ω 4 unchanged (i.e., Ω 1 = 5%, Ω 4 = 2%), and we obtain the effect of the change in the ratios of Ω 3 and Ω 2 on the thermal conductivity of the insulation board.Through calculation and analysis, the influence law curve is obtained, as shown in Figure 9.
According to the influence curve shown in Figure 9, we can find that the thermal conductivity of the GCIB decreases rapidly at the initial stage with the increase in the volume ratio of the graphite polystyrene particles.When the volume ratio of the material components approaches 10 (i.e., Ω 3 /Ω 4 → 10, Ω 3 /Ω 1 → 10, Ω 3 /Ω 2 → 10), the thermal conductivity slowly tends to stabilize, and the ratios of material components in the improved GCIB are in a stable range.

Test Verification of Thermal Conductivity in the GCIB
There are many testing methods for measuring the thermal conductivity of thermal insulation materials [54,55].The thermal conductivity of thermal insulation materials is affected by various factors (e.g., material types, material components, ambient temperature and humidity) [56], which brings more difficulties to the measurement of the materials' thermal conductivity under actual working conditions.Under laboratory conditions, we use the thermal conductivity meter to build a relatively ideal thermal conductivity environment; the steady state method was used to measure the thermal conductivity of the improved GCIB.

Measuring Equipment and Principle
The main equipment for measuring the thermal conductivity of the GCIB by the steady state method is the intelligent thermal conductivity tester (device model: DRCD-3030), which is mainly composed of an electric heater, heating disk and radiator, as shown in Figure 10.The measurement range of the thermal conductivity is 0.010-3.000W/m • K, and the test standard adopted is "Thermosetting composite polystyrene foam insulation board" [57].
The temperature of the heating plate is controlled and measured by the temperature control and measurement sensor.The accuracy of the two temperature measurement sensors is 0.1 °C.The upper and lower end surfaces of the GCIB sample are fully contacted with the heating plate and the cooling plate, respectively, similar to a "sandwich" ' structure.During the test, we ignore the horizontal heat transfer, and the heat is only transferred vertically along the upper and lower directions.After the heat transfer, heat dissipation and heat conduction, the system finally reaches a steady state of heat transfer.When the system reaches the steady state, the heating rate, heat transfer rate and heat dissipation rate of the system are all equal.According to the Fourier heat conduction equation, the thermal conductivity can be expressed as [58] where T 1 and T 2 are the temperatures of the upper and lower surfaces of the test piece in the steady state; is the mass of the heat sink, is the specific heat capacity of the heat sink (3.80 × 10 2 J • kg −1 • K −1 ), a P is the side length of the heat sink, h P is the thickness of the heat sink, a B is the side length of the sample, h B is the thickness of the sample, and dT dt T=T 2 is the heat dissipation coefficient when the temperature of the heat sink is T 2 .

Test Sample Preparation and Test Process
According to the requirements of the material components shown in Table 1, we weigh the materials in the corresponding ratios, add water to mix them to the flow plastic state and prepare 6 samples with the size of through the setting template, and we ensure that the surface roughness is less than ±2% of the thickness.The preparation process is shown in Figure 11.
We placed the 6 test samples in an electric blast drying oven, slowly raised the temperature to (65 ± 5) °C, dried them to a constant mass (the change rate of the sample mass is 0.2% twice at a constant temperature of 3 h) and then moved them to a dryer to cool them to a normal atmospheric temperature to obtain the test samples, as shown in Figure 12.
We used the thermal conductivity tester to measure the thermal conductivity of the six samples, placed the samples between the heating disk and the cooling disk and used the equipment to start the measurement, as shown in Figure 13.When the cold plate temperature, hot plate temperature and protection temperature are equal to the set temperature, the test reaches a stable state, and the time is counted in seconds after the temperature is stable.If the temperature changes and becomes stable again, the time will be counted again.When the power is stable, the thermal conductivity is also stable (three decimal places are stable, and the thermal conductivity of the insulation board can be determined when the stability time is 10,000 s or more).

Test Result
Through the tests on 6 groups of samples, we obtained the process change curve of the thermal conductivity in the GCIB, as shown in Figure 14.When the thermal conductivity test reached a steady state, the temperature of the cold zone was kept at 15 °C and the temperature of the hot zone was kept at 35 °C.The test results show that the temperature of the insulation board basically enters the equilibrium and stable stage in the cold and hot areas in around 4000 s, and the system reaches the stable state in around 12,000 s.We take the average of the thermal conductivity test results of 3000 s in the last stage as the final thermal conductivity of the insulation board, and the test results are shown in Table 3.According to the test results in Table 3, the thermal conductivity of the GCIB is within the range of 0.044-0.050W/(m • K).The test results of 1# − 6# insulation boards were taken as the analysis object.The average thermal conductivity of the thermal insulation board is λ s = 0.047 W/(m • K), and the standard deviation is σ s = 0.002 W/(m • K).  1 into (6), we can obtain the volume ratio of each component material: According to the thermal conductivity of each material in Table 1, substituting the volume ratio in (32) into ( 14), we can obtain the upper limit value of the thermal conductivity in the parallel model of the improved GCIB: By (22), we can obtain the lower limit value of the thermal conductivity in the series model of the improved GCIB: Assuming that the constituent particles of the improved GCIB are uniformly distributed in all directions, and the true value of the thermal conductivity in the insulation board is between λ lb and λ ub , we can approximately take the equivalent thermal conductivity as the average of the upper and lower limit values: (2) Comparison of results Through the analysis of series and parallel theoretical models, the maximum thermal conductivity of the improved GCIB is λ ub = 0.064 W/(m • K), the minimum value is λ lb = 0.042 W/(m • K), and the equivalent value is λ = 0.053 W/(m • K).According to Section 4.2, the calculation results of the thermal conductivity of the two numerical models are λ 400 ≈ λ 10,000 = 0.046 W/(m • K).Through the sample test, we obtain the average thermal conductivity of the improved GCIB as λ s = 0.047 W/(m • K).
Through the comparison of the theoretical analysis, numerical simulation and test results, we found that the sample test results are in good agreement with the numerical simulation analysis results, and the error of the two results is within 0.5%, which is closer to the calculation results of the series theoretical model, between λ ub and λ eq .

Comparison of the Influence of Material Component Ratios
According to (14), (24) and (25), the theoretical thermal conductivity of GCIBs under different material ratios can be obtained.Then, according to Sections 4.3 and 5.3, the theoretical results, numerical simulation results and experimental results are compared, as shown in Figure 15.
According to Figure 15, it is found that the thermal conductivity of the GCIB decreases in two stages as the proportion of graphite polystyrene particles increases.In the first stage, Ω 3 /Ω 4 ≤ 10 (Ω 3 /Ω 1 ≤ 10 or Ω 3 /Ω 2 ≤ 10); as the Ω 3 /Ω 4 (Ω 3 /Ω 1 or Ω 3 /Ω 2 ) ratio increases, the thermal conductivity decreases obviously.In the second stage, Ω 3 /Ω 4 > 10 (Ω 3 /Ω 1 > 10 or Ω 3 /Ω 2 > 10); as the Ω 3 /Ω 4 (Ω 3 /Ω 1 or Ω 3 /Ω 2 ) ratio increases, the thermal conductivity decreases slowly.The results of the numerical simulation and sample test have strong agreement, and the analysis results of the numerical simulation are closer to the lower limit value of the theoretical results.The change curve of the numerical simulation of thermal conductivity is between the theoretical lower limit value curve and the equivalent value curve.

Conclusions
Through theoretical analysis, numerical simulation and experimental tests on the improved GCIB, several conclusions are derived as follows.

1.
The series and parallel models of the GCIB were proposed, which were used to calculate the thermal conductivity of the improved GCIB.The theoretical range of the thermal conductivity of the improved GCIB was obtained.

2.
According to the ratios of each material component in the improved composite insulation board, we established a numerical analysis model of the insulation board with a random distribution of each material component.Through analysis, we obtained the average heat flow density of the thermal conductivity calculation unit.According to the Fourier heat conduction calculation formula, we further obtained the thermal conductivity of the insulation board, which was within the reasonable range of theoretical calculation.

3.
Through numerical simulation, we studied the influence of the volume fraction of graphite polystyrene particles on the thermal conductivity of the GCIB.With the increase in the volume ratio of graphite polystyrene particles, the thermal conductivity of the GCIB decreased rapidly at the initial stage.When the volume ratio approached 10 (i.e., Ω 3 /Ω 4 → 10, Ω 3 /Ω 1 → 10, Ω 3 /Ω 2 → 10), the thermal conductivity value tended to be stable.

4.
The thermal conductivity of the improved GCIB was obtained through the testing of the samples.The test results were within the range of the theoretical calculation results, and the numerical simulation results were in good agreement with the error within 0.5%.
In this paper, the theoretical calculation model of the GCIB was first given.Through numerical simulation and test research, we verified the thermal conductivity of the improved GCIB.The research results have high theoretical and application value.The ratios of each component material in the improved GCIB are within a reasonable range (i.e., Ω 3 /Ω 4 → 10, Ω 3 /Ω 1 → 10, Ω 3 /Ω 2 → 10); however, there is still considerable space for optimization and improvement.We hope to achieve better improvements in future research.

Figure 5 .
Figure 5. Theoretical model of equivalent thermal conductivity of the GCIB: (a) the parallel model; (b) the series model.

Figure 6 .
Figure 6.Geometric model of heat conduction unit of the GCIB.

Figure 7 .
Figure 7. Finite element model of heat conduction unit in the GCIB: (a) 400 element model; (b) 10,000 element model.

Figure 8 .
Figure 8. Calculation results of thermal conductivity of the improved GCIB: (a) the heat flow density cloud chart (400 units); (b) the heat flow density vector chart (400 units); (c) the heat flow density cloud chart (10,000 units); (d) the heat flow density vector chart (10,000 units).

Figure 9 .
Figure 9. Influence curve of thermal conductivity of the GCIB under different component material ratios.

Table 1 .
The constituents of the GCIB.

SiO 2 CaO Al 2 O 3 Fe 2 O 3 MgO SO 3 K 2 O Na 2 O TiO 2 Other Surplus
k g is the gaseous thermal conductivity in the aerogel, and k s is the solid thermal conductivity in the nanospheres.Based on Xie et al.'s [44] model and Karamanos et al.'s [45,46] model, Csanády et al. [47] calculated the thermal conductivity factor of the stem in straw-based thermal insulating materials.

Table 3 .
Thermal conductivity of the test sample.