Three-dimensional Numerical Model Research Articles

Hydrodynamic Model of the Area of the Żelazny Most Mining Waste Storage Facility to Reconstruct the Migration of Saline Groundwater

This paper presents the construction of a numerical three-dimensional model of the area of the Żelazny Most Mining Waste Storage Facility (MWSF). In the study area, the difficult geological conditions associated with glaciotectonics are accompanied by a complex hydrotechnical system of sediment deposition and sedimentary water drainage. In order to effectively reflect the water flow paths, a detailed schematization was carried out, using 700,000 boreholes and more than 300 hydrogeological cross-sections. In addition, numerous drainage sections, streams, and ditches were included to reliably assess the amount of saline water entering the underlying aquifers. This research was supported by magnetic resonance sounding (MRS) studies of the reservoir’s sediments. The MWSF is currently being expanded, so the work primarily focuses on illustrating changes in the hydrodynamic field resulting from the inclusion of the new southern section. Models of similar facilities have been implemented before, but in the current one, the combination of meticulous analysis of the hydro-structural system, the water balance, a significant amount of data, the size of the facility, and the use of an unstructured discretization grid in the calculations is undoubtedly innovative and will be an important contribution to the development of analogous solutions around the world.

Evaluation Method of the Impact of Twin Shield Tunneling Construction on Elevated Bridges: Case Study

In urban metro construction, shield tunneling often needs to pass through building and bridge pile foundations, potentially affecting the stability of existing structures. Therefore, accurately assessing the impact of shield tunneling on bridges and buildings is crucial. This study presents a comprehensive prediction method combining numerical simulation and empirical formulas, taking the underpass project of the Shijiazhuang–Wuhan High-Speed Railway Bridge by Zhengzhou Metro Line 5 as a case study. Three-dimensional numerical model calculations were performed using finite element software to analyze the displacement and stress changes of buildings and tunnel structures at different construction stages, revealing the deformation patterns of buildings adjacent to the tunnel during shield tunneling. In particular, the ground settlement caused by twin-tunnel excavation was compared with Peck’s empirical formula to verify the reliability of the numerical simulation. The results show that twin-tunnel excavation exacerbates the horizontal displacement, uplift, and settlement of the ground, with maximum deformation rates increasing by 7.10%, 20%, and 11.4%, respectively. Comparing the ground deformation results of Peck’s empirical formulas with numerical calculations revealed similar trends in the settlement curves, with a maximum deviation of 6.67%. It can be concluded that using Peck’s empirical formula to calculate ground deformation characteristics complements the limitations of numerical simulations, making the assessment results more reliable. The findings of this study demonstrate that integrating numerical simulation with empirical formulas significantly enhances the reliability of deformation predictions in complex tunneling scenarios. This research not only offers a comprehensive safety assessment method for shield tunneling construction but also provides valuable guidance for the design and construction of similar projects, serving as a theoretical reference for future engineering endeavors.

Mass transfer of bilirubin and bovine serum albumin in hollow fiber membrane module of artificial liver

Understanding the mass transfer behaviors in hollow fiber membrane module of artificial liver is important for improving toxin removal efficiency. A three-dimensional numerical model was established to study the mass transfer of small molecule bilirubin and macromolecule bovine serum albumin (BSA) in the hollow fiber membrane module. Effects of tube-side flow rate, shell-side flow rate, and hollow fiber length on the mass transfer of bilirubin and BSA were discussed. The simulation results showed that the clearance of bilirubin was significantly affected by both convective and diffusive solute transport, while the clearance of macromolecule BSA was dominated by convective solute transport. The clearance rates of bilirubin and BSA increasd with the increase of tube-side flow rate and hollow fiber length. With the increase of shell-side flow rate, the clearance rate of bilirubin first rose rapidly, then slowly rose to an asymptotic value, while the clearance rate of BSA gradually decreased. The results can provide help for designing structures of hollow fiber membrane module and operation parameters of clinical treatment.

3D numerical modelling and laboratory study of flow field induced by a group of submerged vegetations

The three-dimensional (3D) numerical modelling in an open channel flow field of a group of submerged vegetations using computational fluid dynamics (CFD) platform of FLOW-3D HYDRO was performed in this study. A set of acoustic Doppler velocimetry (ADV) measurements have been conducted as benchmark to validate the numerical model. A quantitative comparison was performed on several hydrodynamic variables that impacted the vegetated open channel flow, such as flow depth, streamwise water velocity, turbulent intensity, and Reynolds shear stress. In the numerical analysis, the flow turbulence was treated using the RANS approach (within RNG k-ε); while the Volume Of Fluid (VOF) method was used to track the air-water interface. Structured meshes with hexahedral elements were used to discretize the channel geometry. In the findings, the numerical model reasonably reproduced the flow field and presented corresponding agreement with the experimental turbulent structures. This study showed that the differences in results between various analyses were all less than 10% and concludes that the presented numerical approach can be utilised as an efficient tool for simulations of the flow field within a vegetation patch (i.e. by using the simplified RANS approach).

Study on the damage cause of tunnel junction under different component earthquakes based on on-site investigation and numerical simulation

The Shengli tunnel sustained significant damage during the 2022 Luding earthquake (Mw 6.6), with numerous circumferential and inclined cracks observed in the tunnel lining. This paper develops a three-dimensional numerical model to examine the tunnel junction's dynamic responses to three seismic components: east-west (EW), north-south (NS), and up-down (UD). It also analyzes the causes of tunnel damage. A comparison between the numerical simulation results and the on-site investigation findings shows a high degree of agreement, with damages predominantly concentrated at the junction between the right tunnel and the branch tunnel. The tunnel cross sections' plastic damage areas, damage parameters, K values, and convergent deformation are significantly larger during horizontal vibrations (EW and NS components) than during vertical vibrations (UD component). The stiffness differential between the right and branch tunnels results in varying dynamic responses, leading to geometric incompatibility at the junction as the primary cause of damage. Under horizontal earthquakes, tunnels aligned perpendicular to the vibration direction exhibit significant lateral deformations. Consequently, the extent of damage development in the right tunnel is greater than in the branch tunnel under an EW component earthquake, and the extent of damage development in the branch tunnel is greater than that in the right tunnel under an NS component earthquake. In contrast, damage under the UD component is confined near the junction, without extending into the right or branch tunnels. From an energy perspective, the tunnel system accumulates the least inelastic strain energy during UD component earthquakes, resulting in minimal damage despite the highest Peak Ground Acceleration in the UD direction.

Dynamic heat transfer characteristics of printed circuit precooler in supercritical CO2 Brayton cycle

The dynamic characteristics of heat exchangers are critical in supercritical CO2 Brayton cycle, especially in the precooler where significant thermophysical property variations occur. However, few studies have comprehensively analyzed the dynamic heat transfer characteristics in precooler, particularly with respect to the detailed flow and thermal field evolution. A three-dimensional numerical model was developed to study the dynamic performance of a printed circuit precooler, with an emphasis on the impact of thermophysical properties. The study examines responses to abrupt changes in CO2 and water inlet temperature, as well as mass flux variations. Results show that the precooler exhibits rapid response characteristics with mild fluctuations in heat transfer parameters and flow structure. Notably, the large specific heat near the critical point results in minimal variation in CO2 outlet temperature despite abrupt increases in CO2 inlet temperature, thereby enhancing system stability. Conversely, a decrease in CO2 mass flux causes a transition from forced convection to mixed flow, with strengthened buoyancy effects and decreased entropy production. These variations lead to significant changes in CO2 outlet temperature, posing potential challenges to system stability under CO2 mass flux step change conditions. The heat transfer shows a relatively lower response rate to parameter changes of water compared to CO2. Although variations in water parameters have a negligible effect on CO2 heat transfer, they influence the overall heat transfer coefficient. The dynamic heat transfer characteristics of the precooler presented in this paper address research gaps left by 1D numerical and experimental studies, which have provided limited information on flow and thermal fields during dynamic processes. These insights are essential for understanding the dynamic behavior of the entire system.

Utilizing high-resolution 3D Voronoi meshing to analyze field data from the Brine Availability Test in Salt (BATS)

Salt is an attractive disposal medium for radioactive waste because intact salt is essentially impermeable and non-porous. However, upon drift or borehole excavation a damaged region develops surrounding the excavation which causes increased permeability and porosity creating potential flow paths for brine. Brine leads to corrosion of waste forms and waste packages and is a possible transport vector for radionuclides, so it is important to better understand the early-time behavior and evolution of brine flow in a salt. As a result, this study is part of Task E of DECOVALEX-2023 which focuses on understanding the evolution of thermal, two-phase hydrological, and mechanical processes in the excavation damaged zone in salt. Field measurements from The Brine Availability Test in Salt (BATS) 1a heater experiment are analyzed by implementing a high-resolution three-dimensional numerical model. This salt heater experiment consists of 28 days of heating and 13 days of cooling in a central borehole within bedded salt at the Waste Isolation Pilot Plant (WIPP). Here, the flow simulator PFLOTRAN is utilized; simulations are run on a Voronoi mesh, with temperature-dependent thermal conductivity, permeability and porosity decay away from excavations. The temperature-dependency of permeability is done to match field measurements. Results from the simulation match temperature measured in the field within + /- 0.1 °C and the total brine inflow over the 41-day experiment. This study illustrates that the accuracy of the temperature evolution within salt is critically important when analyzing and modeling experimental data by simulating three heating scenarios of the BATS 1a experiment showing that temperature has a direct effect on total brine inflow.

An assessment of controlled source EM for monitoring subsurface CO2 injection at the wyoming carbonSAFE geologic carbon storage site

We evaluate if electromagnetic (EM) geophysical methods for monitoring geologic carbon storage (GCS) efforts at the Wyoming CarbonSAFE project adjacent to the Dry Fork Station power plant near Gillette, Wyoming. This first involved acquiring both electric and magnetic fields at eleven different locations ranging in distance from immediately adjacent to 4 km from the plant. Passive EM measurements were made to provide spectral EM noise measurements generated by electricity production at the plant and to determine if useful magnetotelluric (MT) data can be successfully collected in the region. The processed data indicate that useful MT data can be collected as long as the site is located more than 2km away from the power plant as well as active roads and rail lines. Controlled source EM data were collected using three different source configurations, two of which connected to steel casings used to complete the injection wells. Comparing the EM noise measurements to the CSEM data show measurable electric and magnetic field signals at all sites. Next a series of three-dimensional (3D) numerical models were built that simulate resistivity changes caused by the proposed CO2 injection at depths ranging from 2.4 to 3.0km. These models were used to simulate various EM measurement configurations. The modeling shows that casing-source CSEM monitoring can provide sensitivity to the injected CO2 if source electrodes are connected to the bottom of one or both of the injection wells.

Analysis of the behavior of piled foundations in unconventional soil

Building foundations using pile groups generally assume that the piles absorb the entire load of the structure. However, it is known that pile caps rest not only on the piles but also on the surrounding soil. These foundations involve great complexity in the phenomenon of load sharing, especially when they are situated on unconventional soils, such as tropical climate soils. To study this behavior, static load tests were conducted on real foundations. These tests were also used to calibrate a three-dimensional finite element numerical model (FEA-3D) using Rocscience's RS3 software. The results showed that the settlement corresponding to the allowable load of the pile groups was significantly lower, ranging from 22 % to 73 %, compared to pile caps that exert contact pressure on the soil. For pile groups, the efficiency ranged between 92 % and 109 %, while for pile caps with contact, it ranged from 132 % to 178 %. In this sense, it becomes evident that estimating the settlement of the pile group based on the conventional design, which does not consider the contact between the cap and the soil, can significantly underestimate the settlement, and provide an overly optimistic analysis of the soil-structure interaction.

Study on the Diffusion Law of Grouting Slurry at the Pile Tip of Bored Piles in Gravel Pebble Layers

The post-construction improvement technology of bored pile tips has been widely used in the construction of gravel pebble layers, but the diffusion law of the grouting slurry remains to be revealed. To study the diffusion law of the grouting slurry at the tip of bored piles in gravel pebble layers, an equivalent relationship between the seepage continuity equation and the mass diffusion continuity equation was established based on an in-depth comparative analysis; further relying on the Diluted Species Transport Module in Porous Media of Comsol Multiphysics, a three-dimensional numerical model of slurry diffusion was established to systematically study the impact of key factors such as different porosities, grouting times, grouting pressures, slurry diffusion coefficients, and the ratio of vertical to horizontal slurry diffusion coefficients on the diffusion radius of the slurry. The calculation results show that with the continuation of grouting, the increase in porosity, grouting pressure, and slurry diffusion coefficient, or the decrease in the ratio of vertical to horizontal slurry diffusion coefficients, the diffusion range of the slurry increases accordingly. Furthermore, comparatively speaking, the impact of porosity is the smallest, while the impact of the slurry diffusion coefficient is the greatest. The research findings reflect the diffusion trend of the slurry in the gravel pebble layer, which has guiding significance for the quality control of the actual design and construction of grouting.

Effect of the pulsed current of on-chip thermoelectric cooler on temperature response during thermal shock in localized cooling

On-chip thermoelectric cooling (TEC) emerges as a promising solution to address the issue of localized chip overheating, effectively mitigating thermal shock induced by local heat flow conduction and enhancing both reliability and performance of chips. In this study, a three-dimensional numerical model is developed to investigate the effectiveness of on-chip thermoelectric cooler’s pulse current in mitigating temperature fluctuations caused by thermal shock at hot spots. The influence of the magnitude and lag time of the pulse current on the temperature fluctuation of Si chip during thermal shock is analyzed. The calculation results show that the implementation of on-chip TEC can effectively mitigate the issue of localized temperature overheating caused by thermal shock. By adjusting the maximum and minimum current of the on-chip thermoelectric cooler, it is possible to mitigate the fluctuation amplitude of chip temperature induced by thermal shock at the hot spot. The findings also indicate that the lag time between input current and thermal shock can lead to an increase in the maximum temperature of the chip. In addition, experimental results using commercially thermoelectric cooler also show that when the input current of TEC is pulse current, the temperature fluctuation caused by thermal shock can be reduced better.

Energy consumption and dilution of toxic gases in underground infrastructures: A case study in a railway tunnel under forced ventilation

Toxic fumes are released after blasting in underground excavations. The analysis of safe re-entry times into the work areas is essential for addressing safety and productivity concerns. In this paper, a railway tunnel 2.1 km long under forced ventilation system is selected as a case study to assess the energy consumption and re-entry times after blasting under different ventilation conditions. First, three-dimensional CFD numerical models were conducted in order to predict the migration time to reduce the levels of noxious gases below the limits established in safety regulations. Field measurements were conducted to verify the accuracy of the results obtained in the simulations. Then, additional analytical models were conducted to evaluate the energy consumed by the ventilation system. The results obtained show that the accumulated re-entry time reaches 414 h for the excavation of the tunnel employing an airflow rate of 40 m3 s−1. This time increases by 33 % when the airflow rate decreases to 30 m3 s−1. However, even though the time decreases when the ventilation volume increases, energy consumption increases significantly, reaching 1,306 MWh for a ventilation volume of 40 m3 s−1 compared to 551 MWh using an airflow of 30 m3 s−1.

Investigation on the heat extraction performance considering a rough fracture network in heterogeneous reservoirs

In this paper, a new three-dimensional numerical model of geothermal extraction is established. The model utilizes spatial frequency to describe the heterogeneity of multiple physical parameters of the reservoir, and rough discrete fractures are constructed based on the Weierstrass-Mandelbrot fractal function. Finally, the study examines the impact of geological-engineering parameters on the heat extraction performance. The results show that after the fracture number increases to 7–9, the influence of fracture number on the low-temperature area gradually diminishes. With the increase in fracture width, the thermal breakthrough phenomenon becomes more pronounced, leading to a reduction in the heat extraction performance. The increase in fracture roughness coefficient can delay thermal breakthrough and improve heat extraction performance. When the fracture roughness coefficient increases from 1.6 to 11.6, the effective electric power increases by 19.44 %. The increase in injection rate intensifies thermal breakthrough, but it also enhances the heat extraction performance. When the injection rate increases from 60 l/s to 160 l/s, the effective electric power increases by 55.7 %. Increasing the injection temperature can delay thermal breakthrough and reduce the heat extraction performance. Abnormal high temperature reservoir (>0.03 K/m) has better heat extraction performance.

Numerical analysis of running resistance of a completion string with packers and sliding sleeves in horizontal wells

The installation of a completion string in small-diameter horizontal wells often presents challenges such as instability, bending, and high frictional resistance. To accurately estimate the running resistance during deployment and ensure safe descent to the bottom of the well, a three-dimensional numerical model has been developed incorporates real-time contact interaction between the completion string and the wellbore, while maintaining precision and field applicability. The analysis evaluates several factors, including buckling behavior, average and dynamics transient value of contact force, axial forces, and frictional forces. Results reveal that the upper pipe string exerting thrust on the lower pipe string causes bending and spiral buckling along with resulting contact with wellbore. Maximum hook load and maximum negative axial force occur at a inclination angle of 65°. Contact strength is more pronounced in inclined and horizontal section compared to vertical section. Leveraging extensive data, the study has yielded a quick prediction formula for quickly calculating the resistance experienced by a completion string of field area. The numerical results offer valuable insights into comprehending the resistance encountered during the descent of a completion string in horizontal wells, and provide robust theoretical support for predicting descent resistance.

Analysis of Longitudinal Cracking and Mold Flux Optimization in High-Speed Continuous Casting of Hyper-Peritectic Steel Thin Slabs

Longitudinal crack defects are a frequent occurrence on the surface of thin slabs during the high-speed continuous casting process. Therefore, this study undertakes a detailed analysis of the solidification characteristics of hyper-peritectic steel thin slabs. By establishing a three-dimensional heat transfer numerical model of the thin slab, the formation mechanism of longitudinal cracks caused by uneven growth of the initial shell is determined. Based on the mechanism of longitudinal crack formation, by adjusting the performance parameters of the mold flux, the contradiction between the heat transfer control and lubrication improvement of the mold flux is fully coordinated, further reducing the incidence of longitudinal cracks on the surface of the casting thin slab. The results show that, using the optimized mold flux, the basicity increases from 1.60 to 1.68, the F- mass fraction increases from 10.67% to 11.22%, the Na2O mass fraction increases from 4.35% to 5.28%, the Li2O mass fraction increases from 0.68% to 0.75%, and the carbon mass fraction reduces from 10.86% to 10.47%. The crystallization performance and rheological properties of the mold flux significantly improve, reducing the heat transfer performance while ensuring the lubrication ability of the molten slag. After optimizing the mold flux, a surface detection system was used to statistically analyze the longitudinal cracks on the surface of the casting thin slab. The proportion of longitudinal cracks (crack length/steel coil length, where each coil produced is about 32 m long) on the surface of the thin slab decreases from 0.056% to 0.031%, and the surface quality of the thin slab significantly improves.

Performance Prediction of GFRP-Reinforced Concrete Deep Beams Containing a Web Opening in the Shear Span

This study aimed to investigate the nonlinear structural behavior of concrete deep beams internally reinforced with glass fiber-reinforced polymer (GFRP) reinforcing bars and containing a web opening of various sizes and locations within the shear span. Three-dimensional (3D) numerical simulation models were developed for large-scale GFRP-reinforced concrete deep beams (300 mm × 1200 mm × 5000 mm) with a shear span-to-depth ratio (a/h) of 1.04. Predictions of the numerical models were validated against published experimental data. A parametric study was conducted to examine the effect of varying the opening size and location on the shear response. Results of the numerical analysis indicated that the strength of the deep beam models with an opening in the middle of the shear span decreased with an increase in either the opening width or height. The rate of the strength reduction caused by increasing the opening height was, however, more significant than that produced by increasing the opening width. Placing a web opening in the compression zone close to the load plate was very detrimental to the beam strength. Conversely, a negligible strength reduction was recorded when the web opening was placed in the tension side above the flexural reinforcement and away from the natural load path. Data of the parametric study were utilized to introduce simplified analytical formulas capable of predicting the shear capacity of GFRP-reinforced concrete deep beams with a web opening in the shear span.

Dynamic response of OWT tripod pile foundation in clay considering scour

Tripod-supported offshore wind turbine (OWT) foundation must withstand wind, wave, and seismic loads. In addition, it may be threatened by scouring during its service life, which could affect its bearing capacity and the overall OWT dynamic response. This study develops a three-dimensional numerical model of OWT tripod pile foundation to investigate scour effects. The numerical model incorporates a simplified single bounding surface model to simulate the dynamic stress–strain response of clay soil. The results revealed that scouring can reduce the first-order natural frequency of tripod pile-OWT system; however, the reduction usually does not lead to resonance within the 1P frequency range. Nonetheless, as the scour depth increases under wind and wave loads, the horizontal displacement of the wind turbine structure increases monotonously; however, its peak acceleration increases initially then decreases slightly as the scour depth continues to increase. Under seismic loading only, the OWT response exhibits complex evolution pattern as the scour depth increases. On the other hand, the foundation horizontal displacement and rotation angle increase significantly under combined wind-wave-seismic loads when the scouring depth reaches three times the pile diameter. Finally, the scouring depth has a relatively small impact on the foundation vertical displacement.

Mechanical mechanism analysis of rockburst in deep-buried tunnel with high in-situ stress

The Qinling water conveyance tunnel has a large buried depth and high in-situ stress level, and rockburst disasters frequently occurred during excavation. In order to find out the mechanical mechanism of rockburst, the research work in this paper is as follows: (1) In-situ three-dimensional hydraulic fracturing method was used to measure the in-situ stress of the deep buried tunnel crossing the ridge. (2) Based on the measured in-situ stress results, the stress distribution characteristics of the tunnel crossing the ridge were obtained by the multiple linear regression method, and the rockburst tendency during construction was predicted. (3) A three-dimensional numerical model of tunnel excavation was established to analyze the dynamic adjustment characteristics of the surrounding rock stress and elastic strain energy during TBM excavation, and to clarify the mechanical mechanism of rockburst. The research results show that the maximum principal stress of the deep-buried tunnel crossing the ridge of Qinling is 40–66 MPa, which belongs to extremely high in-situ stress level, and medium-strong rockburst may occur during excavation. In the process of TBM excavation, the stress of the surrounding rock in the range of 2.6 times the diameter of the tunnel before and after the working face is adjusted violently, and the concentrated zones after the stress redistribution are mainly distributed in the arch roof and arch bottom, and the stress concentration coefficient can reach 2.06. The arch roof, arch waist, and arch bottom are susceptible to immediate rockburst due to stress transient unloading at the moment of excavation. After the elastic strain energy of the surrounding rock at the arch roof and the arch bottom is released and accumulated, it is easy to cause time delayed rockburst, and the depth of the rockburst pit can reach 3.5 m, which is consistent with the rockburst phenomenon in the field.

Bayesian neural networks for solar power forecasts in advanced thermoelectric systems

This study presents a Bayesian neural network for the comprehensive performance forecast of an advanced thermoelectric system featuring segmented thermocouples with variable area shape for optimal solar energy conversion. The training data is obtained from a three-dimensional numerical model developed in ANSYS software to account for the energy, exergy, and entropy flows in the thermoelectric module to maximize the power output, energy, and exergy efficiencies while minimizing the entropy generation and system irreversibilities. The thermodynamic model optimizes key parameters, including geometry (height, cross-section, and percentage skutterudite content) across realistic diverse heat fluxes and cooling coefficients. To facilitate a faster and more efficient optimization pipeline, various neural networks, leveraging algorithms like Levenberg-Marquardt, scaled conjugate gradient, and Bayesian regularization, are utilized for efficient data management and accurate performance prediction. The Bayesian neural network emerges as the most accurate and efficient model providing the lowest training and testing loss of ∼10−9. Furthermore, key findings highlight the importance of skutterudite content and solar concentration ratio. For instance, an optimal skutterudite content of 6.2 % at 25 Suns produces maximum power of 29.7 W, while 37 % skutterudite at 75 Suns minimizes exergy destruction. This study's methodology and findings contribute to designing more efficient clean energy systems and advancing sustainable energy goals.

Parametric study and optimization of H-type finned tube heat exchangers with honeycomb arrangement using Taguchi method

Enhancing the efficiency of flue gas waste heat utilization is an important way to achieve carbon peaking and carbon neutrality goals. In this paper, a three-dimensional numerical model of an H-type finned tube heat exchanger with honeycomb arrangement is developed. Firstly, the effects of five parameters are analyzed using the control variable method, including the oblique tube pitch, the ratio of fin height to oblique tube pitch, fin pitch, fin thickness and slit width. Then, the Taguchi method was used to investigate the interaction effects of the five geometric parameters. Numerical simulations of 27 cases with different parameter combinations were carried out and the heat transfer and flow friction characteristics and their overall performance were discussed in detail. The results show that oblique tube pitch has the most significant effect on heat transfer capacity, fin pitch has the most significant effect on overall performance, and slit width has the least effect on both. Finally, the optimum parameter combinations were obtained. When inlet velocity ranges from 5 m/s to 10 m/s, heat transfer performance is increased by 46.3 %-52.1 %, pressure drop is increased by 5.7 % − 11.3 %, and the overall performance PEC is improved by 63.5 %-105.9 %. The results of the study contribute to the development and optimization of new flue gas heat exchangers.