Numerical Simulation Study Research Articles

Understanding Low-Speed Streaks and Their Function and Control through Movable Shark Scales Acting as a Passive Separation Control Mechanism.

The passive bristling mechanism of the scales on the shortfin mako shark (Isurus oxyrinchus) is hypothesized to play a crucial role in controlling flow separation. In the hypothesized mechanism, the scales are triggered in response to patches of reversed flow at the onset of separation occurring in the low-speed streaks that form in a turbulent boundary layer. The two goals of this investigation were as follows: (1) to measure the reversing flow occurring within the low-speed streaks in a separating turbulent boundary layer; (2) to understand the passive flow control mechanism of movable shark skin scales that inhibit reversing flow within the low-speed streaks. Experiments were conducted using digital particle image velocimetry (DPIV). DPIV was used to analyze the flow in a turbulent boundary layer subjected to an adverse pressure gradient formation over both a smooth flat plate and a flat plate on which shark skin specimens were affixed. The experimental analysis of the flow over the smooth flat plate corroborated the findings of previous direct numerical simulation studies, which indicated that the average spanwise spacing of the low-speed streaks increases in the presence of adverse pressure gradients upstream of the point of separation. However, the characteristics of the flow over the shark skin specimen more closely resemble that of a zero-pressure gradient turbulent boundary layer. A comparative analysis of the width and velocity of the reversed streaks between flat plate and shark skin cases reveals that the mean spanwise spacing decreases, and thus, the number of streaks increases over the shark skin. Additionally, the reversed streaks observed over shark scales are thinner and the highest negative velocity within the streaks falls within the range required to bristle the scales.

Experiment and numerical simulation study of polycarboxylate superplasticizer modified cemented ultrafine tailings filling slurry: Rheology, fluidity, and flow properties in pipeline

Understanding the rheological properties of cemented paste backfill is crucial for optimizing its flowability and ensuring the efficient transportation in mine site. Due to the superfine size, the high-concentration filling slurry prepared by ultrafine tailings (the mean size < 0.019 mm) usually exhibits poor flow performance. In this study, the polycarboxylate superplasticizer (PCS) was adopted as an additive to modify cemented ultrafine tailings filling slurry. Experiment tests and numerical simulation were conducted to investigate effects of PCS dosages, solid and binder contents on rheological, fluidity, and flow performance. Results indicate that the PCS modified filling slurry conforms to Bingham model and shows shear-thinning behavior. As PCS dosage increases from 0 to 0.30 wt%, the yield stress gradually decreases, whereas the plastic viscosity initially increases then decreases. The electrostatic repulsion and steric hindrance generated by PCS can disintegrate the floc structures formed by ultrafine particles, leading to more release of free water, and reducing interparticle friction, correspondingly, improving flowability of filling slurry. The optimal PCS dosage is determined as 0.15–0.25 wt% considering transportation requirement. For a given PCS dosage, the yield stress and plastic viscosity exhibit an exponential growth with solid content, while a continuous decline with binder content. Furthermore, the fluidity increases linearly and exponentially with PCS dosage for filling slurry at solid content of 59–62 wt% and 63 wt%, respectively. An exponential model is established to quantify the relationship between yield stress and fluidity, which provides a method for predicting rheological parameters. With addition of PCS dosage, the flow performance of filling slurry is enhanced significantly in transportation pipeline. The maximal flow velocities at the pipeline center increase, and those near the pipeline wall decrease. The findings afford a practical guidance for designing and optimizing materials proportions of PCS modified filling slurry to achieve the desired flowability in mine site.

Comparative analysis of tube designs in heat exchangers: A numerical simulation study for enhanced thermal-flow efficiency, "optimizing wavy tubes bundle geometries for enhanced heat transfer in underwater applications"

With the growing challenges of heat transfer in marine vehicles and underwater power stations, the demand for an effective tube heat exchanger has surged. Numerical simulations were conducted on the constant, laminar flow of an incompressible viscous fluid inside a wavy tube while keeping the wall temperature consistent. The laminar flow was modeled in two dimensions using ANSYS FLUENT 19.2. For the water fluid model, an integrated algorithm was used to link the pressure and velocity fields. Primarily examined were circular tube banks organized in a staggered arrangement. The primary equations were approached via a finite volume method based on a combined technique. This study aimed to evaluate and contrast the heat exchange and flow traits of four distinct tube designs: asteroid (ASTB), hypotrochoid (HTB), Wavy edge circle (WCTB), and mixed bundle (MBT). The influence of tube design on a heat exchanger’s thermal-flow efficiency was measured across Reynolds numbers ranging from 100 to 400. The results indicated an improved Colburn factor and a decrease in both the friction factor and tube performance for all scenarios examined. Three key metrics determined the optimal tube design: the performance evaluation criterion (PEC), the global performance criterion (GPC), and the average Nusselt number (Nu). Regardless of the Reynolds number, the conventional HC had the least efficient PEC, while the STB registered the peak PEC value. The highest GPC was attained with the STB. Nu values remained relatively consistent for ASTB, HTB, WCTB, and MBT. In evaluating the average Nusselt number, the ASTB outperformed the STB (92–106% higher), HTB (36–48% higher), WCTB (21–32% higher), and MBT (57–64% higher).

Experimental and Analytical Study on Shear Lag Effect of T-Shaped Reinforced Concrete Shear Walls

Because of the flanges on T-shaped shear walls (TSSWs), the shear force acting on such walls results in a shear lag effect, making it impossible to forecast with accuracy the normal stresses of the flanges using the Bernoulli–Euler assumption. Shear lag (SL) in flanged walls has, however, received less attention from researchers, particularly in experimental studies. Understanding the SL in T-shaped reinforced concrete shear walls under shear and axial force is the main goal of this work. First, a SL model is suggested for TSSWs. In this model, the SL deflection is considered to be the generalized displacement and the SL warping deformation, and it is assumed to be a quadratic nonlinear function. Then, experimental and numerical simulation studies are, respectively, conducted to investigate SL effect of TSSWs, and also to evaluate the accuracy of the SL method. Finally, the parameter analysis is conducted to investigate the influence of axial load, shear force, and flange length on the SL effect of TSSWs. The results show that the SL of the TSSW is significant, the normal stress distribution (NSD) of the flange is uneven, and the normal stresses near the web are higher, according to the results of the analytical, simulated, and experimental results. The SL model can accurately predict the normal stresses of the flange of TSSWs, and the quadratic parabola assumption of the SL warp displacement of TSSWs is reasonable. Parameter analysis shows that axial force has little effect on the SL effect of TSSWs. The TSSWs under larger shear force have the more obvious SL effect. A more obvious SL effect occurs in the TSSWs with longer flanges.

Numerical Simulation Study on the Impact of Deep Foundation Pit Excavation on Adjacent Rail Transit Structures—A Case Study

Excavation in foundation pits can result in serious issues for nearby tunnel structures like deformation, differential settlement, and seepage damage, which profoundly impact project timelines and potentially endanger life and property safety. Therefore, it is imperative to investigate these impacts before and after construction and to facilitate timely adjustments of construction measures and reinforcement where possible. In this study, a foundation pit construction project near a rail transit line is employed as a case to comprehensive study the impact of on-site deep foundation pit excavation on adjacent rail transit structures by numerical simulation. A three-dimensional finite-element model of the foundation pit based on site geological characteristics and construction procedures is established to study the excavation and maintenance processes. Through analysis of key parameters including soil deformation, displacement, shear force, and bending moment of the tunnel structures, the designed protective structure is found to have effectively mitigated soil deformation, ensuring the stability of the foundation pit. As excavation progresses, lateral soil deformation and vertical uplift gradually increase but remain within specified control values. During various excavation stages, the maximum displacement of the tunnel structure gradually increases, with the increase rates of maximum settlement being 29.09%, 20.51%, and 6.45%, respectively. This indicates a gradual enhancement of the stability of the tunnel structure. Additionally, excavation of the foundation pit has a significant impact on the bending moment distribution of the tunnel structure but does not affect the axial force and shear force of the tunnel structure. The findings of this study offer crucial scientific insights for evaluating the safety and stability of construction near tunnel structures.

Numerical Simulation Study on the Deformation Patterns of Surrounding Rock in Deeply Buried Roadways under Seepage Action

To reveal the deformation patterns of the surrounding rock in deeply buried straight-wall arch-shaped roadways under seepage action, this study, based on an FLAC3D numerical simulation and classic elastoplastic theory, investigates the influences of surrounding rock classification, roadway burial depth, pore water pressure, and roadway cross-sectional dimensions on the deformation of surrounding rock. A multivariate regression prediction model for rock deformation was established based on the numerical simulation conclusions, and the correctness of the conclusions was verified through comparative analysis. Correlation analysis of various factors with rock deformation was conducted, ranking their impact as follows: pore water pressure &gt; roadway burial depth &gt; surrounding rock classification &gt; roadway height &gt; roadway width. The research results can provide guidance for the construction and support of deeply buried roadways under seepage action.

Numerical Simulation Study on the Influence of Fracture on Borehole Wave Modes: Insights from Fracture Width, Filling Condition, and Acoustic Frequency.

Unconventional reservoirs, such as shale and tight formations, have become increasingly vital contributors to oil and gas production. In these reservoirs, fractures serve as crucial spaces for fluid migration and storage, making their precise assessment essential. Array acoustic logging stands out as a pivotal method for evaluating fractures. To investigate the impact of fracture width, fracture-filling conditions, and acoustic frequency on compressional and shear waves, a three-dimensional variable mesh finite difference program was employed for acoustic logging numerical simulation. Firstly, numerical models representing fractured formations with varying fracture widths and distinct fluid-filling conditions were established, and array acoustic logging numerical simulations were conducted at different frequencies. Subsequently, the waveform data were processed to extract acoustic characteristic parameters, such as velocities and amplitude attenuations of compressional and shear waves. Finally, a quantitative analysis was conducted to examine the variation patterns of characteristic parameters of refracted compressional and shear waves in relation to fracture properties. The research results indicate that amplitude attenuation information derived from borehole wave modes is particularly sensitive to the changes in fracture properties. As fracture width increased, we observed a significant amplitude attenuation in both compressional and shear waves, proportional to the logarithm of the attenuation coefficients. Furthermore, when the fracture width was constant, gas-filled fractures exhibited more prominent amplitude attenuation than water-filled fractures, with shear wave attenuation being more sensitive to the filling material. Moreover, from a quantitative perspective, the analysis revealed that the attenuation coefficients of refracted compressional and shear waves exhibited an exponential variation with gas saturation. Notably, once fracture width and filling conditions were established, the amplitudes of compressional and shear waves at the dominant frequency of 40 kHz were significantly reduced compared to those at 8 kHz, accompanied by increased attenuation. Subsequent quantitative analysis revealed that, when the product of fracture width and dominant frequency remains constant, the corresponding attenuation coefficient ratios approach 1. This indicates that the attenuation process of acoustic propagation in fractured media follows the principle of acoustic similarity. The findings of this study provide reference for further research on fracture property evaluation methods based on array acoustic logging data.

Exploring thermal buoyant flow in urban street canyons: Influence of approaching turbulent boundary layer

Turbulent boundary layer inflow is a critical factor in urban climate research, shaping canyon flow dynamics, air ventilation patterns, and heat flux distribution. In numerical simulation studies, it serves as a fundamental inflow boundary condition, profoundly influencing overall results. In this study, simultaneous Particle Image Velocimetry and Laser-Induced Fluorescence (PIV-LIF) measurements are utilized within a large closed-circuit water tunnel. This approach allows comprehensive flow data to be gathered under varied flow and thermal conditions, encompassing a spectrum of Richardson numbers ranging from 0.01 to 1.34. The investigation aims to elucidate the effects of turbulent boundary layer flows on heat transfer mechanisms and flow behaviours within a two-dimensional street canyon model with a unit aspect ratio. The analysis reveals distinct heat and fluid flow characteristics, highlighting the interplay between thermal conditions and flow dynamics. The three chosen turbulent boundary layer flows demonstrate unique influences on flow characteristics and heat removal capacity. Significant variations in ventilation rates are observed, with a maximum difference of 80% among the tested boundary layer flows. Additionally, the most pronounced variation in heat removal capacity is approximately 45%. Thicker boundary layers with lower velocities near the canyon exhibit reduced ventilation and heat removal capabilities. Furthermore, the investigation reveals that varied turbulence inlet profiles result in diverse fluctuating features at the canyon roof level, with a comparatively lesser impact on the deeper regions of the canyon.

An efficient and robust real-time longitudinal ventilation control method for unpredictable tunnel fire scenarios

A novel real-time control method, named SF-B control method, was proposed and a comprehensive numerical simulation study was conducted for feasibility analysis. The simulation utilized a combination of Python scripts and FDS. The SF-B control method incorporates automatic detection and localization of the fire source, and the existing ventilation velocity is corrected based on the real-time state of motion of smoke fronts and the back-layering length. The control logic of the SF-B control method was described in detail, and the criterion for determining stabilization was provided. After a comparison with PID control method, scenarios involving random constant HHRs, random fire locations and variable HRR were set up for separate simulations, and the results were compared with previous experimental findings. The outcomes demonstrate that the fluctuation of the SF-B control method is only about 1/10 of that of the PID control method, and the control performance is not significantly influenced by the HRR or fire location. It shows a more powerful control ability on smoke back-layering in scenarios involving variable HRR. This method exhibits excellent compatibility with different fire conditions, showcasing high ventilation control efficiency, robust control capabilities, and delivering reliable control results.

Study on the optimization of heating exchanger in electric heating solid energy storage heating system

This paper addresses prevalent issues of suboptimal compatibility between the heating exchanger and the thermal storage unit, poor safety performance, and overall insufficient heat exchange efficiency within the application of heating exchangers in electric heating solid energy storage heating systems. Experimental testing and numerical simulation studies were conducted. The research investigates the effect of the temperature of inlet air as well as velocity on the heat exchange performance of the heating exchanger as well as temperature variations of single-row heat pipes. Drawing upon these change patterns, an optimized heating exchanger structure is proposed and subsequently investigated through optimization simulation studies. The study results indicate that the best overall optimization effect is achieved with a heating exchanger arranged with finned tube combinations of 4 mm in two rows, 6 mm in two rows, 8 mm in two rows, 10 mm in two rows, 12 mm in two rows, and 14 mm in ten rows arranged successively from front to back. When the heating exchanger’s inlet air speed is relatively high, this combination’s heat exchange capacity surpasses the original structure. Additionally, the uniformity of air-side temperature drop improved by 44.89%, while the finned area was reduced by 25.62%.

Numerical Simulation and Experimental Study of Gas–Solid Two-Phase Spraying of Dry Powder Fire-Extinguishing System Based on Fire-Extinguishing Inspection Robot

In order to solve the problem where the traditional intelligent inspection robot only has a single inspection function, we studied the use of a dry powder (including an ultra-fine dry powder) as a fire-extinguishing medium for the first time. In fire-extinguishing robots, the spray pressure is difficult to control, and there are several other issues. For integrated inspection, an intelligent, nitrogen-driven fire-extinguishing robot using a dry powder in a pressure-controlled spray was developed. On this basis, in order to investigate nitrogen-driven dry powder particle spraying as a gas–solid two-phase mechanism, as well as the flow characteristics and the influence of relevant parameters on the spraying effect, a nitrogen-driven dry powder particle spraying system was established as part of a gas–solid two-phase computational fluid dynamics model. The flow field of the spraying system and the particle motion characteristics were analyzed to explore the micro-mechanisms of the influence of different driving pressures, pipe diameters, and nozzle configurations on the spraying of the dry powder. In order to investigate the macroscopic effect of dry powder spraying where the gas–solid two-phase micro-mechanisms could not be revealed, an experimental platform was set up, and the experiments verified the accuracy of the numerical simulation results. We also investigated the dry powder spraying effect under different driving pressures, pipe diameters, nozzle configurations, and loading ratios. Finally, an orthogonal test was designed based on the results of the single-factor experiments to find the best combination of parameters required to achieve the optimal spraying effect. The research results can provide a theoretical and technical reference for the design and development of nitrogen-driven dry powder spraying systems.

Numerical simulation and experimental study on effect of cooling rate on microstructure and strength of nanostructured materials

In this study, by numerical simulation (finite element method, FEM) and experimental, the cooling rate was investigated by changing the product thickness (20, 3, 2, 1, 0.5 and 0.3) mm of Al based two-phase nanostructured materials casted through a copper mold. The effect of cooling rate on the microstructure and strength of the alloy was studied. The Experimental results showed that the precipitated intermetallic phases have a decreasing size corresponding to the increasing cooling rates by simulation from ~102 K/s to ~104 K/s. The results show that an appropriate cooling rate can improve the microstructure and properties of the alloy. The Abaqus/Standard capability for uncoupled heat transfer analysis was intended to model solid body heat conduction with general, temperature-dependent conductivity; internal energy (including latent heat effects); and quite general convection and radiation boundary conditions. This study describes the basic energy balance, constitutive models, boundary conditions, finite element discretization, and time integration procedures used. The time step used an automatic algorithm through the smallest tolerance. The maximum temperature change was allowed over a period and the increment was adjusted for this parameter, as was the rate of convergence in the non-linear cases. First-order heat transfer elements used the rule of numerical integration with integrated stations located at the corners of the element for thermal capacitance terms. (Jacobian terminology). This approach is particularly effective when there is a strong latent thermal effect. Thus, first - order elements were used in the case of latent heat. The HEATCAP element is available for single - point pooled thermal capacitance modeling. Centralized film loading options between the mold and the casting were specified by the user.

Numerical Simulation and Crack Suppression Study on High-Speed Laser Cladding of ZL101 Aluminum Alloy

Numerical Simulation and Crack Suppression Study on High-Speed Laser Cladding of ZL101 Aluminum Alloy

Effect of oxy-fuel combustion with different O2/CO2 fractions on the pulverized coal combustion mechanism and heat transfer characteristics in rotary kilns: A numerical simulation study

Oxy-fuel combustion combined with carbon capture and storage is the most effective and direct technology to achieve carbon neutrality in cement industry. In this work, numerical simulation was employed to investigate the effects of air combustion and oxy-fuel combustion with different O2/CO2 fractions on the airflow field, temperature field, pulverized coal combustion mechanism, heat transfer characteristics and NOx distribution in rotary kiln. The results show that the airflow field under different conditions has no significant change. When transitioning from the air condition to the oxy-fuel condition, the ignition point of pulverized coal extends backward by approximately 0.2 m. The maximum temperature decreases from 2306 K to 2052 K, and the temperature becomes more uniform. As the oxygen concentration in the primary air increases, both the maximum and average temperatures gradually rise, and the flame length shortens from 25.4 m to 24.16 m. In the rotary kiln, the locations of char oxidation and gasification reactions differ, but these reaction locations remain consistent under different conditions. Under oxy-fuel conditions, the oxidation rate and oxidation proportion of char decrease, while the gasification rate and gasification proportion of char increase, resulting in a shortening of the burnout distance of coal from 34.5 m to 31.3 m. With increasing oxygen concentration in the primary air, the rates of char oxidation and gasification gradually rise, further shortening the burnout distance. The decrease in temperature at the front end of the rotary kiln under oxy-fuel conditions reduces the total heat transfer in the burning zone by 15 %. When the oxygen content of the primary air is 70 %, it ensures stable calcination of clinker. Additionally, under oxy-fuel conditions, the NOx generation is significantly reduced. With the increase of oxygen concentration in the primary air, the NOx concentration at the outlet decreases slightly from 451 ppm to 406 ppm, due to the increase in kiln temperature that favors the reduction of fuel NOx.

Desulfurization characteristics of slaked lime and regulation optimization of circulating fluidized bed flue gas desulfurization process—A combined experimental and numerical simulation study

Circulating fluidized bed flue gas desulfurization (CFB-FGD) process has been widely applied in recent years. However, high cost caused by the use of high-quality slaked lime and difficult operation due to the complex flow field are two issues which have received great attention. Accordingly, a laboratory-scale fluidized bed reactor was constructed to investigate the effects of physical properties and external conditions on desulfurization performance of slaked lime, and the conclusions were tried out in an industrial-scale CFB-FGD tower. After that, a numerical model of the tower was established based on computational particle fluid dynamics (CPFD) and two-film theory. After comparison and validation with actual operation data, the effects of operating parameters on gas–solid distribution and desulfurization characteristics were investigated. The results of experiments and industrial trials showed that the use of slaked lime with a calcium hydroxide content of approximately 80% and particle size greater than 40 μm could significantly reduce the cost of desulfurizer. Simulation results showed that the flow field in the desulfurization tower was skewed under the influence of circulating ash. We obtained optimal operating conditions of 7.5 kg·s−1 for the atomized water flow, 70 kg·s−1 for circulating ash flow, and 0.56 kg·s−1 for slaked lime flow, with desulfurization efficiency reaching 98.19% and the exit flue gas meeting the ultraclean emission and safety requirements. All parameters selected in the simulation were based on engineering examples and had certain application reference significance.

Numerical simulation study on the effect of underground drainage pipe network in typical urban flood

Rapid urbanization and climate change have heightened urban flood risks. An in-depth examination of the underground drainage pipe network’s effect on urban flood inundation can enhance urban stormwater management and mitigate disaster risks. This study presents a full hydrodynamic urban flood model (FHUM) coupling the pipe transport module of the Storm Water Management Model (SWMM) with a shallow water model. Based on unstructured grid cells, runoff yield and confluence calculations are carried out to effectively separate the underground drainage pipe network from other factors. The model is applied to Minzhi Area and Haidian Island to quantitatively assess the influence of the underground drainage system on urban flooding. Results indicate that FHUM performs well, accurately reflecting the process from rainfall to inundation. The underground drainage pipe network reduces surface inundation in the study area and improves the city’s stormwater drainage capacity. With increasing rainstorm intensity, the impact of the pipe network does not exhibit a consistent trend. During heavy rainstorms, the pipe network consistently plays a significant and stable role, particularly in reducing high inundation depths. However, attention should be paid to the potential flood risk transfer associated with using pipe networks in urban flood management. This work provides a scientific basis for urban stormwater management, holding significant importance in improving flood control and disaster reduction capabilities.

Numerical simulation and experimental study on axial compression electromagnetic bulging of aluminum alloy tube

Numerical simulation and experimental study on axial compression electromagnetic bulging of aluminum alloy tube

Numerical simulation and theoretical study on the impact of wind-sand flow of high-speed trains in long tunnel space

When high-speed trains (HST) run in enclosed spaces such as long tunnels, the thermal accumulation of their suspension devices is continuous and cannot be effectively dissipated. In addition, previous experiments or simulations for the heat dissipation of HST in tunnel spaces did not consider the impact of sand. To clarify the impact of HWS-LT on the heat accumulation of HST equipment cabin, this study used the CFD method to numerically simulate the impact of different wind-sand flow concentrations or no-sand wind on the cooling of equipment in the long tunnel space. Firstly, the sand particles in the wind-sand flow gather at the tunnel entrance and enter the equipment cabin with the train as it enters the tunnel. This boundary condition is more in line with actual engineering situations. Secondly, both flows show asymmetric intrusion into the cabin due to the asymmetrical tunnel arrangement, but the sand particles in the wind-sand flow are affected by the vortices and tunnel walls, resulting in more asymmetric flow and some particles being trapped in the grids or filters, leading to outflow ρQ < inflow ρQ. Under the wind-sand flow condition, the temperature of some equipment surfaces shows more significant increases than under the no-sand wind. Finally, contrary to popular perception, the wind-sand flow carrying sand particles can dissipate heat more effectively than no-sand wind, and the higher the volume fraction φ within a certain concentration range, the better the heat dissipation effect. This is because the wind-sand flow has a higher specific heat capacity, which can remove some heat from the contact point between the sand particles and the equipment wall upon contact. The higher sand particle concentration increases the contact frequency and contact area between the sand particles and the equipment wall, and the heat transfer pathway and heat dissipation efficiency are improved.

Numerical simulation study of combined shot peening and laser shock peening on surface integrity of Ti-6Al-4V titanium alloy

Combined shot peening and laser shock peening (CSP) can introduce high surface compressive residual stress (CRS) and deep CRS layer, and reduce the high surface roughness caused by shot peening. The relationship between CSP sequences and surface integrity, including CRS and surface roughness, is far from satisfying the actual demand. In this paper, we construct a finite element model considering initial surface profile features. The experimental results agree with the simulated results. Based on the simulation results, axial and longitudinal stress wave transfer effects on surface CRS and depth CRS are investigated. Compared with shot peening (SP) and laser shock peening (LSP) alone, CSP treatment can increase CRS depth and CRS magnitude by 1060 μm and 50 ∼ 80Mpa, respectively. In addition, the roughness following CSP treatment is lower than that observed after SP treatment alone, and the roughness varies depending on the sequence of CSP treatment. Specifically, the roughness values for SP, SP + LSP, and LSP + SP are 0.938 ± 0.031 μm, 0.982 ± 0.11 μm, and 1.168 ± 0.093 μm, respectively. The surface morphology, hardness, and simulation results are combined to reveal the reasons for the reduction of surface roughness after CSP treatment. The changes of peaks and valleys after CSP treatment compared with single SP treatment were studied, and the variation trend of surface roughness was further explained.

Numerical simulation and experimental study of three-phase distribution characteristics of leaked light non-aqueous phase liquid from buried pipelines in soils containing groundwater and gas.

Leakage accidents of buried pipelines have become increasingly common due to the prolonged service of some pipelines which have been in use for more than 150years. Therefore, there is an urgent need for accurate prediction of pollution scope to aid in the development of emergency remediation strategies. This study investigated the distribution of a light non-aqueous phase liquid in soils containing gas and water through numerical simulations and laboratory experiments. Firstly, a three-dimensional porous medium model was established using ANSYS FLUENT, and for the first time, the distribution of gas and groundwater in soil environments was simulated in the model. Subsequently, the distribution of the three phases of diesel, gas, and water in soil was studied with different leakage velocities and it was found that the leakage velocity played a significant role in the distribution. The areas of diesel in soils at 60min were 0.112m2, 0.194m2, 0.217m2, and 0.252m2, with corresponding volumes of 0.028m3, 0.070m3, 0.086m3, and 0.106m3, respectively, for leakage velocities of 1.3m/s, 3.4m/s, 4.6m/s, and 4.9m/s. Calculation formulas for distribution areas and volumes were also developed to aid in future prevention and control strategies under different leakage velocities. The study also compared the distribution areas and volumes of diesel in soils with and without groundwater, and it was found that distribution scopes were larger in soils containing groundwater due to capillary force. In order to validate the accuracy of the numerical simulation, laboratory experiments were conducted to study the diffusion of oil, gas, and water under different leakage velocities. The results showed good agreement between the experiments and the simulations. The research findings are of great significance for preventing soil pollution and provide a theoretical basis for developing scientifically sound soil remediation strategies.