CFD continuous simulation method and device for pedestrian height microclimate and medium

Abstract

The invention provides a CFD continuous simulation method and device for pedestrian height microclimate and a medium, and the method comprises the steps: determining a target region and a research time period, and dividing the target region and the research time period into a plurality of continuous moments; obtaining meteorological measurement experiment data including boundary conditions and verification actual measurement data through a meteorological measurement experiment; establishing a regional three-dimensional geometric model and an air calculation domain, and performing grid division; on the basis of the boundary condition data, CFD simulation parameters of multiple continuous moments are set, high-spatial-resolution CFD continuous moment simulation is carried out, and a wind thermal field simulation result is obtained; and finally, comparing and verifying a simulation result with actually measured data. According to the invention, high-resolution and high-efficiency simulation of the urban microclimate dynamic process is realized, and the accuracy and feasibility of pedestrian height microclimate simulation are ensured.

Description

Technical Field

[0001] This document relates to the field of microclimate technology, and in particular to a CFD continuous simulation method, equipment and medium for pedestrian height microclimate. Background Technology

[0002] Urban microclimate monitoring methods include sensor network monitoring, satellite remote sensing monitoring, UAV remote sensing monitoring, and CFD numerical simulation monitoring. Sensor network monitoring involves deploying various sensors in urban areas to directly measure physical quantities related to urban microclimate, such as temperature, humidity, wind speed, and wind direction. Sensor network monitoring can achieve high-precision real-time climate monitoring and provide real-time response and early warning, but its coverage is limited, maintenance costs are high, and it is susceptible to interference. Satellite remote sensing monitoring uses sensors carried by satellites to acquire urban microclimate information from a macroscopic perspective. Satellite remote sensing can cover large urban areas and is not limited by ground conditions, but its spatial resolution is low, it requires complex data inversion, and vertical information acquisition is difficult. UAV monitoring uses sensors carried by UAVs to measure physical quantities related to urban microclimate at low altitudes. UAV monitoring is highly flexible, has high imaging resolution, and can achieve precise low-altitude detection, but its monitoring range is relatively small, its endurance is limited, and it is greatly affected by weather. High spatial resolution CFD numerical simulation monitoring is based on physical models and mathematical algorithms to simulate and predict urban microclimate. For numerical simulation, if the model simplification assumptions are made, the simulation parameters are determined, and the verification is reasonable, then it can not only delve into the physical mechanisms and have high resolution, but also be free from time and space limitations, and have low cost while also being able to make predictions and analyses, which is far superior to other urban microclimate monitoring methods.

[0003] In the field of CFD numerical simulation monitoring, there are two main simulation methods: steady-state simulation and transient simulation. These two methods face an inherent contradiction in urban microclimate simulation due to the time scale gap. While transient simulation can capture dynamic changes in the flow field, its time step is limited by physical stability, typically requiring microsecond to second-level steps. Simulating temperature changes over a period of time places extremely high demands on computing power, even exceeding the practical limits of supercomputing. Steady-state simulation provides single-moment results by converging to an equilibrium state, but it cannot reflect the dynamic processes of urban microclimate over a period of time, such as the alternation of day and night and the periodic changes in solar radiation, leading to the neglect of key phenomena such as daily temperature fluctuations and thermal hysteresis effects.

[0004] Pedestrian height is the primary living space for urban residents, therefore, the microclimate characteristics at this height directly affect their quality of life, comfort, and even health. Consequently, simulating and predicting urban microclimate, especially during periods involving dynamic processes such as day-night cycles and cyclical changes in solar radiation, is crucial for researchers' understanding and analysis of urban microclimate characteristics.

[0005] In summary, there is an urgent need for a high spatial resolution CFD simulation method that can simulate urban microclimate over a period of time with a simulation resolution of at least one meter, and can ensure the accuracy and feasibility of simulating urban microclimate at pedestrian height. Summary of the Invention

[0006] According to embodiments of the present invention, a CFD continuous simulation method, device, and medium for pedestrian height microclimate are provided, aiming to solve the above-mentioned problems.

[0007] According to an embodiment of the present invention, a CFD continuous simulation method for pedestrian height microclimate is provided, comprising: S1. Determine the target region and study time period for CFD continuous time simulation, and divide the study time period into a continuous time periods; S2. During the research period, meteorological measurement experiments are conducted to obtain meteorological measurement experimental data, which includes boundary condition data for CFD simulation and measured data for result verification. S3. Establish a three-dimensional geometric model of the city and an air computation domain for the target area, and perform spatial discretization and meshing on the air computation domain; S4. Based on the boundary condition data, set the CFD simulation boundary conditions and parameters for a consecutive time periods, and perform high spatial resolution CFD simulation at consecutive time periods to obtain the wind and heat field simulation results of the target area. S5. Compare and verify the simulated wind-thermal field results with the measured data; The boundary conditions used for CFD simulation include: inlet boundary conditions, outlet boundary conditions, and wall boundary conditions.

[0008] According to an embodiment of the present invention, an electronic device is provided, comprising: Processor; and, A memory is configured to store computer-executable instructions, which, when executed, cause the processor to perform the steps of the aforementioned CFD continuous simulation method for pedestrian height microclimate.

[0009] According to an embodiment of the present invention, a storage medium is provided for storing computer-executable instructions, which, when executed, implement the steps of the above-described CFD continuous simulation method for pedestrian height microclimate.

[0010] By employing the embodiments of this invention, the contradiction between the excessively high computational requirements of transient simulation and the inability of steady-state simulation to reflect temporal dynamics in traditional CFD simulations is effectively resolved. Through continuous steady-state simulation, the dynamic processes of urban microclimates are efficiently captured. The organic combination of field measurements and CFD simulation ensures both the reliability of the boundary conditions in the simulation process and the accuracy of the simulation results through system verification. This achieves meter-level high spatial resolution urban microclimate simulation, accurately reflecting the microclimate characteristics of the key activity space of pedestrian height. It provides a reliable technical means for urban heat island effect analysis, urban planning and design, and outdoor environmental assessment, and has significant engineering application value. Attached Figure Description

[0011] To more clearly illustrate the technical solutions in one or more embodiments of this specification or in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments recorded in this specification. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0012] Figure 1 This is a flowchart of a CFD continuous simulation method for pedestrian height microclimate according to an embodiment of the present invention; Figure 2 This is a schematic diagram of a CFD continuous simulation method for pedestrian height microclimate, which is a specific implementation of an embodiment of the present invention. Figure 3 This is a schematic diagram of the three-dimensional geometric model of the city in the study area according to an embodiment of the present invention; Figure 4 This is a distribution map of the experimental measuring points for the meteorological measurement experiment according to an embodiment of the present invention; Figure 5 This is a schematic diagram of the air computational domain of the research area in an embodiment of the present invention; Figure 6 This is a schematic diagram of mesh division according to an embodiment of the present invention; Figure 7 The temperature cloud map at a pedestrian height of 1.5m at each moment is the result of the high spatial resolution CFD continuous time simulation in this embodiment of the invention. Figure 8 The wind speed cloud map at a pedestrian height of 1.5m at each moment is the result of continuous time simulation by high spatial resolution CFD in an embodiment of the present invention. Figure 9 This is a comparison chart of the temperature error between experimental data and simulation results at a pedestrian height of 2m, as shown in this embodiment of the invention. Detailed Implementation

[0013] To enable those skilled in the art to better understand the technical solutions in one or more embodiments of this specification, the technical solutions in one or more embodiments of this specification will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of this specification, and not all of the embodiments. Based on one or more embodiments of this specification, all other embodiments obtained by those skilled in the art without creative effort should fall within the protection scope of this document.

[0014] Method Implementation Examples According to embodiments of the present invention, a CFD continuous simulation method for pedestrian height microclimate is provided. Figure 1 This is a flowchart of a CFD continuous simulation method for pedestrian height microclimate according to an embodiment of the present invention. Figure 2 This is a schematic diagram of a CFD continuous simulation method for pedestrian height microclimate, a specific implementation of an embodiment of the present invention. Figure 1 and Figure 2 As shown, the CFD continuous simulation method for pedestrian height microclimate according to an embodiment of the present invention specifically includes: S1. Determine the target region and study time period for CFD continuous time simulation, and divide the study time period into a continuous time periods; Specifically, select a micro-scale urban densely built-up area, where the horizontal dimension is less than 2 kilometers.

[0015] The study period for simulating pedestrian height in urban microclimate in the research area was determined to be ≥1 hour. The study period was divided into (a-1) smaller time periods, i.e. a (a≥3) consecutive moments.

[0016] S2. During the research period, meteorological measurement experiments are conducted to obtain meteorological measurement experimental data, which includes boundary condition data for CFD simulation and measured data for result verification; the boundary conditions for CFD simulation include: inlet boundary conditions, outlet boundary conditions and wall boundary conditions.

[0017] The meteorological measurement experiment was carried out through a deployed measurement system, which included m representative measuring points in the region and n incoming wind measuring points; The m representative measuring points in the region are used to collect air temperature and wind speed data at heights of 0 meters, 1.5 meters and 2 meters at a consecutive time, and simultaneously collect wall temperature data of building exterior surfaces, concrete floors and asphalt pavements at a consecutive time. The n incoming air measurement points are used to collect incoming air temperature and wind speed data at heights of 0 meters, 1.5 meters, and 2 meters at a continuous time.

[0018] Specifically, during the study period, meteorological measurement experiments were conducted. m (m≥4) representative experimental measurement points and n (n≥1) incoming wind measurement points were selected. m personnel were distributed across the m measurement points in the study area, using handheld anemometers and thermometers to measure air temperature and wind speed at heights of 0m, 1.5m, and 2m at a consecutive time intervals. n personnel were distributed across the n incoming wind measurement points, using handheld anemometers and thermometers to measure air temperature and wind speed at heights of 0m, 1.5m, and 2m at a consecutive time intervals. Simultaneously, m personnel were distributed across the m measurement points in the study area, using high-definition infrared thermal imagers to measure the temperature of common urban surfaces (building exteriors, concrete floors, and asphalt pavements) at a consecutive time intervals. During each measurement, multiple measurements were recorded continuously after the readings stabilized, and the average value was taken to avoid reading errors.

[0019] S3. Establish a three-dimensional geometric model of the city and an air computation domain for the target area, and perform spatial discretization and meshing on the air computation domain; Specifically, remote sensing technology was used to establish a three-dimensional urban geometric model of the study area based on OSM road network data and building geometry data. This model was then imported into the geometry software SpaceClaim, and an air computational domain was established for the study area based on the direction of airflow. The principle for establishing the air computational domain was that the building obstruction rate was less than 3%. The air computational domain was discretized using the mesh generation software Workbench meshing. The surface mesh element size (≤1m) was set according to the minimum geometric dimensions that the three-dimensional urban geometric model of the study area needed to capture, and the mesh was generated. The mesh quality was continuously optimized based on the principle that the maximum mesh skew coefficient should not exceed 0.85.

[0020] S4. Based on the boundary condition data, set the CFD simulation boundary conditions and parameters for a consecutive time periods, and perform high spatial resolution CFD simulation at consecutive time periods to obtain the wind and heat field simulation results of the target area. S4 specifically includes: establishing boundary condition models for high spatial resolution CFD continuous-time simulations. The boundary condition models include inlet boundary, outlet boundary, and wall boundary.

[0021] For the inlet boundary conditions of the air fluid domain, the velocity inlet boundary conditions are adopted. It is necessary to set the inlet air temperature and wind speed at a consecutive time, as shown in equations (1) and (2).

[0022] (1); (2); in, The inlet boundary temperature matrix; The wind speed matrix at the inlet boundary; The input air temperature at the inlet boundary at time a is in K, and its value is the average air temperature at heights of 0m, 1.5m, and 2m at n incoming wind measurement points at time a, obtained from meteorological measurement experiments. The input wind speed at the inlet boundary at time a is in m / s. The value is the average wind speed at heights of 0m, 1.5m, and 2m at n incoming wind measurement points at time a, obtained from meteorological measurement experiments.

[0023] For the outlet boundary conditions of the air fluid domain, the pressure outlet boundary conditions are adopted. It is necessary to set the flow outlet air temperature at a consecutive time, as shown in equation (3).

[0024] (3); in, This is the temperature matrix at the export boundary.

[0025] For the wall boundary of the air fluid domain, including three common wall surfaces, namely the building exterior surface, concrete ground and asphalt pavement, a mixed wall boundary condition is adopted. It is necessary to set the incoming air temperature, wall temperature, convective heat transfer coefficient and external emissivity of each wall surface at a continuous time.

[0026] The first type of wall boundary is the outer surface of the building, as shown in equations (4)-(7). The external emissivity is taken as the external emissivity value of the building's outer surface material.

[0027] (4); (5); (6); (7); in, The air temperature matrix for the incoming airflow to the building's exterior surface; This is the temperature matrix of the building's exterior wall surface. The thermal conductivity matrix of air on the building's exterior surface; The matrix of convective heat transfer coefficients on the exterior surface of the building; The distance from the first-layer grid node of the airflow domain on the building's outer surface to the building's outer surface; The air temperature at the building's exterior surface at time 'a' is expressed in K. The value is the average air temperature at 0m and 1.5m heights of m representative experimental measurement points at time 'a' obtained from meteorological measurement experiments. The value represents the building's exterior wall temperature at time 'a', expressed in K. It is the average of the building's exterior wall temperatures at m representative experimental points at time 'a' obtained from meteorological measurements.

[0028] The second type of wall boundary is a concrete floor, as shown in equations (8)-(11), where the external emissivity is taken as the external emissivity value of the concrete floor material.

[0029] (8); (9); (10); (11); in, The air temperature matrix for the inflow to the concrete floor; This is the temperature matrix of the concrete floor surface; The thermal conductivity matrix of air for concrete floors; The matrix represents the convective heat transfer coefficient of the concrete floor. The distance from the first layer of grid nodes in the airflow domain of the concrete floor to the wall surface; The air temperature at the concrete ground at time 'a' is expressed in K. The value is the average air temperature at 0m and 1.5m heights of m representative experimental measurement points at time 'a' obtained from meteorological measurement experiments. The value represents the concrete floor surface temperature at time 'a', expressed in K. It is the average concrete floor surface temperature at m representative experimental points at time 'a' obtained from meteorological measurements.

[0030] The third type of wall boundary is asphalt pavement, as shown in equations (12)-(15), where the external emissivity is taken as the external emissivity value of the asphalt pavement material.

[0031] (12); (13); (14); (15); in, An air temperature matrix for asphalt pavement; This is the temperature matrix of the asphalt pavement wall. The air thermal conductivity matrix of asphalt pavement; The matrix represents the convective heat transfer coefficient of asphalt pavement. The distance from the wall to the first layer of grid nodes in the airflow domain of the asphalt pavement; The air temperature at time a represents the air temperature at the asphalt pavement, in K, and is the average air temperature at 0m and 1.5m heights of m representative experimental measuring points at time a obtained from meteorological measurement experiments. The value represents the asphalt pavement wall temperature at time 'a', in K, and is the average wall temperature of the asphalt pavement at m representative experimental points at time 'a' obtained from meteorological measurement experiments.

[0032] Parameter settings for high spatial resolution CFD continuous-time simulations. This includes setting the solver and solution method, configuring the turbulence model, enabling the energy equation, setting material properties, setting symmetric boundary conditions for all air-fluid domain boundaries except for the inlet, outlet, and wall boundaries, setting the equation discretization scheme, setting the convergence residual criterion, performing mixed initialization, and setting the number of iterations.

[0033] A combined Python and Fluent secondary development approach was used to perform continuous time-series CFD simulations of urban microclimates over a given period. The pyansys and pyfluent packages were installed in Python, and a for loop was used to drive functions within these packages, thereby enabling Fluent to perform high spatial resolution continuous-time steady-state CFD simulations of the urban microclimate in the study area over a consecutive time period.

[0034] S5. Compare and verify the simulated wind-thermal field results with the measured data; Specifically, using CFD-post, simulated temperature data at a height of 2m from m representative experimental measurement points in the high spatial resolution CFD continuous time simulation results of urban microclimate are extracted and compared with temperature data at a height of 2m from m representative experimental measurement points in the meteorological measurement experiment. The absolute and relative errors between the simulated and experimental temperature values ​​are calculated, and the relative error of the temperature at each measurement point is required to be less than 10%, in order to verify the feasibility and accuracy of the high spatial resolution CFD continuous time simulation results of urban pedestrian height microclimate.

[0035] The following specific example illustrates the CFD continuous simulation method for pedestrian height microclimate according to this invention. The study period is from 11 PM on September 1, 2024 to 10 PM on September 2, 2024. This period is divided into 23 smaller time intervals, resulting in 24 consecutive moments. The meteorological measurement experiment included four representative regional measurement points and one incoming wind measurement point. The minimum surface grid element size was 0.5m. , , The values ​​were 1m, 1m, and 0.5m, the number of iterations was 80, and the for loop was used 24 times. When verifying the results, the maximum relative error of the temperature at each measuring point was less than 7%. Figure 3 To study the three-dimensional geometric model of the city in the region, Figure 4 This is a map showing the distribution of experimental measuring points for meteorological measurements. Figure 5 To establish the air computational domain for the research area, Figure 6 This is a diagram showing the results of the grid division. Figure 7 This is a temperature cloud map at a pedestrian height of 1.5m at each time step, based on the results of a high spatial resolution CFD continuous-time simulation. Figure 8This is a wind speed cloud map at a pedestrian height of 1.5m at each time step, based on the continuous time-series simulation results of high spatial resolution CFD. Figure 9 This is a comparison chart of the temperature error between experimental data and simulation results at a pedestrian height of 2m.

[0036] In summary, this invention can simulate the time-by-time urban microclimate of densely built-up urban areas at a microscale (horizontal dimension less than 2 kilometers) using a high spatial resolution CFD continuous time simulation method, and the relative error of the simulation of the urban microclimate at pedestrian height is less than 10%, demonstrating extremely high accuracy and feasibility.

[0037] By employing the embodiments of the present invention, the following beneficial effects are achieved: This method effectively resolves the contradiction between the excessively high computational requirements of transient simulations and the inability of steady-state simulations to reflect temporal dynamics in traditional CFD simulations. By performing steady-state simulations at continuous time intervals, it achieves efficient capture of urban microclimate dynamic processes. Through the organic combination of field measurements and CFD simulations, it ensures the reliability of boundary conditions in the simulation process and guarantees the accuracy of simulation results through system verification. It achieves urban microclimate simulation with meter-level high spatial resolution, accurately reflecting the microclimate characteristics of the key activity space of pedestrian height. It provides a reliable technical means for urban heat island effect analysis, urban planning and design, and outdoor environmental assessment, and has significant engineering application value.

[0038] Device Example 1 According to an embodiment of the present invention, an electronic device is provided, comprising: Processor; and, A memory is configured to store computer-executable instructions that, when executed, cause the processor to perform the steps of the method embodiments described above.

[0039] Device Example 2 According to an embodiment of the present invention, a storage medium is provided for storing computer-executable instructions, which, when executed, implement the steps of the method embodiments described above.

[0040] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some or all of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention.

Claims

1. A CFD continuous simulation method for pedestrian height microclimate, characterized in that, include: S1. Determine the target region and study time period for CFD continuous time simulation, and divide the study time period into a continuous time periods; S2. During the research period, meteorological measurement experiments are conducted to obtain meteorological measurement experimental data, which includes boundary condition data for CFD simulation and measured data for result verification. S3. Establish a three-dimensional geometric model of the city and an air computation domain for the target area, and perform spatial discretization and meshing on the air computation domain; S4. Based on the boundary condition data, set the CFD simulation boundary conditions and parameters for a consecutive time periods, and perform high spatial resolution CFD simulation at consecutive time periods to obtain the wind and heat field simulation results of the target area. S5. Compare and verify the simulated wind-thermal field results with the measured data; The boundary conditions used for CFD simulation include: inlet boundary conditions, outlet boundary conditions, and wall boundary conditions.

2. The method according to claim 1, characterized in that, The meteorological measurement experiment was carried out through a deployed measurement system, which included m representative measuring points in the region and n incoming wind measuring points. The m representative measuring points in the region are used to collect air temperature and wind speed data at heights of 0 meters, 1.5 meters and 2 meters at a consecutive time, and simultaneously collect wall temperature data of building exterior surfaces, concrete floors and asphalt pavements at a consecutive time. The n incoming air measurement points are used to collect incoming air temperature and wind speed data at heights of 0 meters, 1.5 meters, and 2 meters at a continuous time.

3. The method according to claim 1, characterized in that, The inlet boundary is configured with the inlet air temperature and wind speed at a continuous time interval. （1）； （2）； in, The inlet boundary temperature matrix, The wind speed matrix at the inlet boundary. Let be the input air temperature at the inlet boundary at time 'a', in K, and be the average air temperature at heights of 0m, 1.5m, and 2m at n incoming wind measurement points at time 'a' obtained from a meteorological measurement experiment. The input wind speed at the inlet boundary at time a is in m / s. The value is the average wind speed at heights of 0m, 1.5m, and 2m at n incoming wind measurement points at time a, obtained from meteorological measurement experiments.

4. The method according to claim 1, characterized in that, The outlet boundary condition is defined by setting the air temperature at the outlet for *a* consecutive time points. （3）； in, This is the temperature matrix at the export boundary.

5. The method according to claim 1, characterized in that, The wall boundary conditions include the boundary conditions of the building's outer surface, concrete floor, and asphalt pavement; For the building exterior surface, see Formula (4)-(7), the external emissivity is taken as the external emissivity value of the building exterior surface material; （4）； （5）； （6）； （7）； in, The air temperature matrix for the incoming airflow to the building's exterior surface; This is the temperature matrix of the building's exterior wall surface. The thermal conductivity matrix of air on the building's exterior surface; The matrix of convective heat transfer coefficients on the exterior surface of the building; The distance from the first-layer grid node of the airflow domain on the building's outer surface to the building's outer surface; The air temperature at the building's exterior surface at time 'a' is expressed in K. The value is the average air temperature at 0m and 1.5m heights of m representative experimental measurement points at time 'a' obtained from meteorological measurement experiments. The value represents the building exterior wall temperature at time 'a', in K, and is the average of the building exterior wall temperatures at m representative experimental points at time 'a' obtained from meteorological measurement experiments. For concrete floors, see formulas (8)-(11), where the external emissivity is taken as the external emissivity value of the concrete floor material. （8）； （9）； （10）； （11）； in, The air temperature matrix for the inflow to the concrete floor; This is the temperature matrix of the concrete floor surface; The thermal conductivity matrix of air for concrete floors; The matrix represents the convective heat transfer coefficient of the concrete floor. The distance from the first layer of grid nodes in the airflow domain of the concrete floor to the wall surface; The air temperature at the concrete ground at time 'a' is expressed in K. The value is the average air temperature at 0m and 1.5m heights of m representative experimental measurement points at time 'a' obtained from meteorological measurement experiments. The value represents the concrete floor surface temperature at time 'a', in K, and is the average of the concrete floor surface temperatures at m representative experimental points at time 'a' obtained from meteorological measurement experiments. For asphalt pavement, see formulas (12)-(15), the external emissivity is taken as the external emissivity value of the asphalt pavement material; （12）； （13）； （14）； （15）； in, An air temperature matrix for asphalt pavement; This is the temperature matrix of the asphalt pavement wall. The air thermal conductivity matrix of asphalt pavement; The matrix represents the convective heat transfer coefficient of asphalt pavement. The distance from the wall to the first layer of grid nodes in the airflow domain of the asphalt pavement; The air temperature at time a represents the air temperature at the asphalt pavement, in K, and is the average air temperature at 0m and 1.5m heights of m representative experimental measuring points at time a obtained from meteorological measurement experiments. The value represents the asphalt pavement wall temperature at time 'a', in K, and is the average wall temperature of the asphalt pavement at m representative experimental points at time 'a' obtained from meteorological measurement experiments.

6. The method according to claim 1, characterized in that, When establishing the air computation domain in S3, the principle is to ensure that the building's obstruction effect on the incoming airflow is at a low level. When performing mesh generation, the surface mesh resolution is set based on the minimum geometric feature size that needs to be captured in the city's three-dimensional geometric model, and mesh optimization is performed with the principle of controlling the mesh quality within the acceptable range of computation.

7. The method according to claim 1, characterized in that, When performing high spatial resolution CFD continuous-time simulation in S4, the CFD solver is driven by an automated process through the coupling interface between the external script program and the CFD solver. The CFD solver reads the preset boundary conditions for a consecutive time periods in sequence and performs the corresponding steady-state simulation calculations, thereby realizing continuous-time simulation of urban microclimate changes within the study period.

8. The method according to claim 1, characterized in that, In step S5, when comparing and verifying the simulation results with the meteorological measurement experimental data, wind and heat environment parameters within the pedestrian height range corresponding to the spatial location of the representative experimental measurement points in the region are extracted from the simulation results and quantitatively compared and analyzed with the measured data of the corresponding experimental measurement points. The accuracy of the high spatial resolution CFD continuous time simulation results is comprehensively evaluated based on the preset error evaluation criteria.

9. An electronic device, comprising: processor; as well as, A memory configured to store computer-executable instructions, which, when executed, cause the processor to perform the steps of the CFD continuous simulation method for pedestrian height microclimate as described in any one of claims 1-8.

10. A storage medium for storing computer-executable instructions, which, when executed, implement the steps of the CFD continuous simulation method for pedestrian height microclimate as described in any one of claims 1-8.

Images