Chinese  |  English
Home Table of Contents

15 December 2024, Volume 64 Issue 12
    

  • Select all
    |
    SPECIAL SECTION: BIG DATA
  • WU Houyue, LI Xianwei, ZHANG Shunxiang, ZHU Honghao, WANG Ting
    Journal of Tsinghua University(Science and Technology). 2024, 64(12): 1997-2006. https://doi.org/10.16511/j.cnki.qhdxxb.2024.21.027
    Abstract ( ) Download PDF ( ) HTML ( )   Knowledge map   Save
    [Objective] The generation of adversarial samples in text represents a significant area of research in natural language processing. The process is employed to test the robustness of machine learning models and has gained widespread attention from scholars. Owing to the complex nature of Chinese semantics, generating Chinese adversarial samples remains a major challenge. Traditional methods for generating Chinese adversarial samples mainly involve word replacement, deletion/insertion, and word order adjustment. These methods often produce samples that are easily detectable and have low attack success rates, and thus, the methods struggle to balance attack effectiveness and semantic coherence. To address these limitations, this study introduces DiffuAdv, a novel method for generating Chinese adversarial samples. This approach enhances the generation process by simulating the data distribution during the adversarial attack phase. The gradient changes between adversarial and original samples are used as guiding conditions during the model's reverse diffusion phase in pre-training, resulting in the generation of more natural and effective adversarial samples. [Methods] DiffuAdv entails the introduction of diffusion models into the generation of adversarial samples to improve attack success rates while ensuring the naturalness of the generated text. This method utilizes a gradient-guided diffusion process, leveraging gradient information between original and adversarial samples as guiding conditions. It consists of two stages: forward diffusion and reverse diffusion. In the forward diffusion stage, noise is progressively added to the original data until a noise-dominated state is achieved. The reverse diffusion stage involves the reconstruction of samples, in which the gradient changes between adversarial and original samples are leveraged to maximize the adversarial objective. During the pre-training phase, data capture and feature learning occur under gradient guidance, with the aim of learning the data distribution of original samples and analyzing the deviations from adversarial samples. In the reverse diffusion generation phase, adversarial perturbations are constructed using gradients and integrated into the reverse diffusion process, ensuring that at each step of reverse diffusion, samples evolve toward greater adversarial effectiveness. To validate the effectiveness of the proposed method, extensive experiments are conducted across multiple datasets and various natural language processing tasks, and the performance of the method is compared with those of seven existing state-of-the-art methods. [Results] Compared with existing methods for generating Chinese adversarial samples, DiffuAdv demonstrates higher attack success rates across three tasks: text sentiment classification, causal relation extraction, and sentiment cause extraction. Ablation experiments confirm the effectiveness of using gradient changes between original and adversarial samples to guide the generation of adversarial samples and improve their quality. Perplexity (PPL) measurements indicate that the adversarial samples generated by DiffuAdv have an average PPL value of only 0.518, demonstrating that these samples are superior in rationality and readability compared with the samples generated by other methods. [Conclusions] DiffuAdv effectively generates high-quality adversarial samples that closely resemble real text in terms of fluency and naturalness. The adversarial samples produced by this method not only achieve high attack success rates but also exhibit strong robustness. The introduction of DiffuAdv enhances the research perspective on generating adversarial text samples and broadens the approaches for tasks such as text sentiment classification, causal relationship extraction, and emotion-cause pair extraction.
  • LI Jiayi, HUANG Ruizhang, CHEN Yanping, LIN Chuan, QIN Yongbin
    Journal of Tsinghua University(Science and Technology). 2024, 64(12): 2007-2018. https://doi.org/10.16511/j.cnki.qhdxxb.2024.21.028
    Abstract ( ) Download PDF ( ) HTML ( )   Knowledge map   Save
    [Objective] The increasing maturity of large language model technology has facilitated its widespread application in downstream tasks across various vertical fields. Large language models have exhibited beneficial performance in text summarization tasks in general fields, such as news and art. However, the highly specific language style in the judicial field and the unique complexity of judicial documents in terms of structure and logic make it difficult for large language models to generate judicial document summaries. This study aims to combine prompt learning with large language models to explore their performance in summarizing judicial documents. Prompt templates containing structural information and judicial documents are used as inputs for fine-tuning large language models. As a result, large language models can generate judicial document summaries that adhere to judicial language styles and the structural and logical complexities of judicial documents. [Methods] This study proposes a judicial document summary method that combines prompt learning and the Qwen large language model. Judicial document data are used as the input for fine-tuning a large language model using supervised fine-tuning technology to enhance its applicability in the judicial field. Simultaneously, prompt templates that incorporate structural information and role instructions are designed to optimize summary generation to more accurately reflect the structural characteristics and logical relationships of documents. According to the characteristics of the pretraining data format of the large language model, the fine-tuning data were constructed in the form of question-answer pairs. [Results] The experimental results show that the proposed method improves the F1 of the baseline model by 21.44%, 28.50%, and 28.97% in ROUGE-1, ROUGE-2, and ROUGE-L, respectively, and exceeds all of the comparison models. The ablation experiment demonstrated that the summary generation method using prompt learning was superior to the method without prompt learning for all indicators, and the performance of summarization generated by the large language model utilizing prompt learning was significantly enhanced. The case demonstration reveals that after prompt learning is used to enhance the perception of structural information in the judgment document by the large language model, the judgment document summary generated by this model can better capture and retain key information in the judgment document. Moreover, the language style of this model is closer to that of a real judgment document summary, which further illustrates the effectiveness of the proposed method. [Conclusions] This study integrates the structural information of a judgment document into the task of generating a judgment document summary using a large language model in the form of prompt templates. Prompt templates containing structural information are used to assist the large language model in summarization generation. Therefore, the model can focus on the key information in the judgment document and capture deeper semantic logical relationships. The results demonstrate that after fine-tuning the large language model with judicial document data and introducing structural information, the model demonstrated excellent performance and great application potential in the judicial document summary task. The proposed method can effectively enhance the capability of a large language model in the field of judicial document summaries.
  • SPECIAL SECTION: NUCLEAR WASTE WATER AND GAS
  • YANG Li, ZHANG Yujie, FANG Sheng, SONG Jiayue, LI Xinpeng, CHEN Yixue
    Journal of Tsinghua University(Science and Technology). 2024, 64(12): 2019-2030. https://doi.org/10.16511/j.cnki.qhdxxb.2024.27.029
    Abstract ( ) Download PDF ( ) HTML ( )   Knowledge map   Save
    [Objective] Local-scale atmospheric dispersion modeling of radionuclides is crucial for nuclear emergency response during the early phase. The Lagrangian puff dispersion model excels in accurately and rapidly reproducing radioactive fields at this scale by accounting for natural turbulence and integrating wind fields with spatial and temporal variations. Given that nuclear power plants (NPPs), especially Chinese NPPs, are often located in heterogeneous terrains, which lead to channeling and slope flows, puff splitting in the puff dispersion model is necessary to accurately represent the phenomenon of plume splitting and layer decoupling phenomena. Despite its importance, the threshold values for puff splitting have not been adequately studied. In addition, the complex terrain around NPP sites generates highly complicated flows, necessitating the use of a diagnostic wind field model coupled with the atmospheric dispersion model to improve the accuracy of dispersion simulations. [Methods] To further provide an effective atmospheric dispersion modeling and establish threshold values of puff splitting for the Lagrangian puff dispersion model, the local-scale Lagrangian splitting puff dispersion model (SPUFF) was developed and fully integrated with the California meteorological model (CALMET). Two local-scale dispersion simulations were conducted using the CALMET to drive the SPUFF: one against the Sanmen NPP wind tunnel experiments with east (E) and northeast (NE) wind directions and another to simulate the Fukushima Daiichi nuclear accident. These simulations aimed to validate SPUFF's performance and practicality. Furthermore, a comprehensive sensitivity analysis was performed to determine the credible range of horizontal threshold values for puff splitting. The dispersion results were evaluated using multiple statistical metrics: the fraction of simulations within a factor of two/five/ten of the observations (FAC2/5/10), fractional mean bias (FB), normalized mean-square error (NMSE), normalized absolute difference (NAD), and geometric mean bias (MG). [Results] Validation results indicated that plumes generated by SPUFF effectively covered the majority of measurement sites, with coverage rates reaching 99.60% and 97.54% in the E and NE directions, respectively. All four crucial statistical metrics for SPUFF met acceptable criteria (FAC2: 0.52, FB: -0.17; NMSE: 0.75, NAD: 0.31 in the E direction; FAC2: 0.48, FB: 0.37; NMSE: 1.28, NAD: 0.39 in the NE direction), indicating remarkable performance. Practical evaluations demonstrated that SPUFF can reproduce more measurements in the Futaba station compared to the Lagrangian particle model (LAPMOD). SPUFF also successfully captured the concentration peak effects resulting from the reactor events during the Fukushima nuclear accident. Sensitivity analysis suggested that applying no puff splitting module might be sufficient for complex terrains with constant meteorological conditions (constant wind fields). However, puff splitting becomes crucial in complex terrains with variable meteorological conditions. For local-scale dispersion scenarios involving NPPs, the recommended threshold values for puff splitting range between 700 m and 1 100 m. [Conclusions] This paper provides a comprehensive evaluation of the Lagrangian splitting puff dispersion model (SPUFF) and demonstrates its practical application. The results strongly indicate that SPUFF is a valuable tool for future nuclear emergency responses. Additionally, this paper proposes a credible range of threshold values for puff splitting, offering guidelines for applying the puff model in local-scale dispersion scenarios at NPP sites.
  • WANG Xuebin, WU Dong, YIN Yuguo, LIU Yaming, RUAN Hao
    Journal of Tsinghua University(Science and Technology). 2024, 64(12): 2031-2044. https://doi.org/10.16511/j.cnki.qhdxxb.2024.27.019
    Abstract ( ) Download PDF ( ) HTML ( )   Knowledge map   Save
    [Significance] In nuclear reactors and spent fuel reprocessing plants, the production of tritiated light water is unavoidable, amounting to thousands of tons annually. The direct discharge of this byproduct into the environment poses significant ecological risks. Consequently, strict tritium emission standards have been established worldwide, propelling the development of detritiation technologies. Among these, combined electrolysis and catalytic exchange (CECE) technology stands out because of its high detritiation factor and mild operating conditions, positioning it as a focal point in global research. This study explores the current state of CECE technology, highlights the three key technologies that underpin it, and addresses the challenges faced in its engineering application, thereby promoting its practical implementation. [Progress] CECE technology comprises liquid-phase catalytic exchange (LPCE), electrolysis, and hydrogen-oxygen recombination processes. LPCE technology is instrumental in the operation of CECE technology. The LPCE column, a critical component, operates on complicated principles, and its efficiency is influenced by various factors such as temperature, pressure, and packing material. Research conducted over the years has shed light on the effect of these elements on the performance of LPCE columns. Electrolysis technology serves as the bottom reflux mechanism within the CECE, with alkaline electrolyzers and proton exchange membrane (PEM) electrolyzers as the main devices. Alkaline electrolyzers, characterized by their limited liquid inventory, good tritium radiation resistance, and high operational stability, are widely regarded as mature technologies. Efforts are currently being directed toward increasing gas production capabilities. PEM electrolyzers represent a new area of development. Compared to alkaline electrolyzers, their notable advantage lies in the absence of alkaline electrolytes. However, their susceptibility to tritium poses significant challenges to their widespread application. Hydrogen-oxygen recombination technology is the top reflux technology of the CECE technology, with the recombiner device playing a pivotal role. Recent advancements have seen a transition from hydrophilic to hydrophobic catalysts within the recombiners, coupled with a reduction in the reaction temperature by over 100℃ while maintaining an efficiency rate exceeding 99.9%. Concurrently, theoretical simulations of CECE technology have evolved with the development of models such as two-film mass transfer and three-fluid models, alongside simulation programs such as FLOSHEET and EVIO. These tools have been instrumental in guiding the design of the CECE process by combining theoretical simulations and experimental analyses. With the development of theoretical simulations and key technologies underpinning CECE, several countries have designed processing schemes to remove tritium from tritiated light water using CECE technology. This study details the process proposed by Canada and Japan. [Conclusions and Prospects] Advances in the key technologies of CECE demonstrate significant advantages in removing tritium from tritiated light water. Moreover, there is substantial potential for further development in engineering applications. Furthermore, the efficiency and cost-effectiveness of CECE technology can be further improved in several ways: the development of more efficient and economical catalysts, the enhancement of PEM electrolyzers to offer better resistance to tritium irradiation, increased gas production, and advancements in hydrogen fuel cell technology.
  • FENG Qingliang, LONG Zeyu, CAO Shaofei, GONG Wenjing, SONG QinNan
    Journal of Tsinghua University(Science and Technology). 2024, 64(12): 2045-2052. https://doi.org/10.16511/j.cnki.qhdxxb.2024.27.030
    Abstract ( ) Download PDF ( ) HTML ( )   Knowledge map   Save
    [Objective] Although nuclear energy provides economic benefits, it releases substantial radioactive substances into nearby waters, thereby affecting the ecological environment. Aquatic organisms in radioactive waters accumulate nuclides, reaching a stable concentration within their bodies over time. At this point, the concentration factor, combined with the nuclide concentration in aquatic organisms, can be used to estimate the radioactivity level in target waters. Radiation doses at other trophic levels can be estimated through the food chain. During the transmission process of the food chain, different degrees of radionuclide accumulation were observed at various fish trophic levels. In addition, at the same trophic level, radionuclide accumulation usually differs in fish individuals with various feeding and metabolic rates. Therefore, the concentration factor of different fish species, especially those with large differences in body size, varies considerably. Direct application of the concentration factor recommended by the International Atomic Energy Agency may lead to large deviations in the calculated results. Zebrafish are commonly used aquatic model organisms in the laboratory. This typically small fish has a length of 3—5 cm and is sensitive to changes in the water environment, making it an ideal indicator organism for environmental monitoring. Therefore, this study considered zebrafish as the research subject to explore the accumulation and discharge kinetics of 137Cs, which is one of the main radionuclides in the liquid effluents of nuclear power plants, to provide a reference for assessing the impact of nuclear power development on aquatic organisms. [Methods] Zebrafish were bred in radioactive water via a specially designed breeding method and equipment, and 137Cs activity in live zebrafish was measured periodically through high-purity germanium gamma spectrometry. The reliability of the living efficiency scale was verified through the gray treatment of zebrafish postmortem. [Results] The results showed a 8.16% average absolute relative error between measurements obtained from living and deceased zebrafish, which confirmed the scale's reliability. Accumulation and discharge experiments revealed that under approximately 1 000 Bq/L radioactivity, 137Cs after being rapidly accumulated in the zebrafish during the initiation of the experiment and 137Cs activity in zebrafish increases linearly. Then, after approximately 50 days, the value remained unchanged and reached equilibrium. At the beginning of the discharge experiment, the 137Cs activity in zebrafish decreased rapidly. However, from around the 60th day, the decreased trend of 137Cs activity became insignificant, and the curve flattened in later stages. [Conclusions] Zebrafish exhibited relatively low 137Cs accumulation with a concentration factor of (5.81±0.23)L/kg at equilibrium. The 137Cs discharged by zebrafish exhibited an initial rapid change, followed by a slow rate of change. When all the zebrafish died on the 198th day, the specific activity of 137Cs remained at approximately 2 370 Bq/kg. This finding suggests the prolonged discharge of 137Cs in zebrafish, indicating a relatively low metabolic rate and difficult discharge. Zebrafish are sensitive to changes in the water environment, and the accumulation and discharge characteristics of 137Cs can provide a reference for investigating the influence of liquid effluents from nuclear power plants on other aquatic organisms and assessing the water environment.
  • CHEN Jiachen, CHEN Hailong, LIAN Bing, WANG Yan
    Journal of Tsinghua University(Science and Technology). 2024, 64(12): 2053-2058. https://doi.org/10.16511/j.cnki.qhdxxb.2024.27.033
    Abstract ( ) Download PDF ( ) HTML ( )   Knowledge map   Save
    [Objective] Radioactive materials may enter the soil through atmospheric deposition, surface water, and groundwater pathways, resulting in radioactive contamination. Physical disturbances caused by natural winds or anthropogenic activities may cause particulate matter in the soil to be resuspended from the ground into the air. Some nuclides (e.g., uranium) that are difficult to transfer to the human body through the food chain are likely to be inhaled in the form of resuspended particulate matter and cause internal irradiation. Therefore, the resuspension behavior of radioactive particulate matter must be investigated for the control of suspended particles. CFD(computational fluid dynamics) studies are mostly carried out on the resuspension of particulate matter in the pipe wall; only a few CFD numerical simulation studies have been conducted on the resuspension release of particulate matter from contaminated soil. [Methods] In this work, a CFD simulation of the resuspension release of radioactive aerosol particles from a contaminated site is carried out. First, the resuspension of radioactive aerosol particles caused by the air inlet is achieved, and the influence of the air inlet and exhaust outlet settings on the distribution characteristics of the airflow field is then analyzed. In the simulation, the air supply port is set as the velocity inlet, and the air velocity of the air inlet is set to 10 m/s. The air exhaust port is set as the outflow boundary, the wall is set as the solid-wall boundary condition, and the ground surface is set as the surface source for generating aerosol particles. The generated aerosol particles have a size of 2.5 μm and a density of 1.65 g/cm3. The discrete phase model is chosen to solve the gas-solid two-phase flow, and the transient state is used to calculate the motion-diffusion characteristics of the resuspended aerosol particles. To determine the state of particulate matter and the trajectory of the motion with time and space, the coupled SIMPLE algorithm of pressure and velocity is selected with the standard pressure and the second-order on-wind discrete format, and the turbulence model is adopted as the standard k-ε model. According to the resuspension motion of the aerosol particles, the ventilation of the air inlet violently disturbs the airflow in the box, forming an approximate left-right symmetric airflow field in the box. [Results] Simulation results show that: 1) when the inlet wind speed is 10 m/s, the airflow inside the resuspension box is uniformly distributed, and the vortex flow is sufficient to resuspend the radioactive aerosol particles on the ground. In addition, the amount of resuspension of particles is large, which is convenient for aerosol sampling and measurements. 2) The simulated wind speed of the resuspension box near the ground is within the range of the ground surface wind speed in the area and thus can be applied to characterize the effect of surface wind speed on the resuspension of particles in this area. [Conclusions] This simulation provides a basis for the experimental design of the deposition and resuspension of radioactive aerosol particles. The findings provide guidance for the radiation protection of people working or moving in radiation-contaminated areas.
  • MECHANICAL ENGINEERING
  • YAN Bo, HU Jinchun, ZHU Yu, WEN Tingrui, XU Dengfeng
    Journal of Tsinghua University(Science and Technology). 2024, 64(12): 2059-2067. https://doi.org/10.16511/j.cnki.qhdxxb.2024.21.025
    Abstract ( ) Download PDF ( ) HTML ( )   Knowledge map   Save
    [Objective] Accurate pose measurement is crucial for precise motion control in multi-degree-of-freedom motors. Traditional methods for pose measurement often rely on external sensors such as encoders, inertial measurement units, and optical sensors, which can be complex and less reliable. This paper introduces a three-degree-of-freedom pose measurement technique that uses redundant magnetic field information. By utilizing the motor's intrinsic magnetic field, this approach aims to simplify the system design and improve its robustness. [Methods] The proposed method employs an array of redundant magnetic sensors to detect the magnetic field generated by the motor's rotor permanent magnets. A measurement algorithm is developed to quickly calculate the motor's multi-degree-of-freedom pose from the detected magnetic flux density data. In addition, a calibration algorithm is introduced to ensure real-time, precise measurement of the motor's pose by accurately aligning the magnetic field model parameters with the actual magnetic field generated by the motor. This calibration process leverages simulation data and a linear magnetic field model to achieve accurate alignment. To further enhance accuracy, spatial transformation matrices are used to construct a mathematical model of the magnetic sensor signals. This allows the system to effectively map the detected magnetic field information to the motor's pose. The Gauss-Newton method is then employed to solve the overdetermined equations arising from the redundant information provided by multiple sensors, and multi-degree-of-freedom pose measurement and position parameter calibration are realized. [Results] A significant advantage of this method is its reliance on the motor's intrinsic magnetic field for pose measurement. This intrinsic measurement approach simplifies the overall system structure and enhances its robustness. By measuring the pose based on the motor's magnetic field, the system can simultaneously determine the installation position of the sensor array and calibrate the parameter offsets of the permanent magnets. This dual functionality is achieved without any special requirements for the sensor arrangement, providing versatility and adaptability to different motor configurations. In addition, this method is less susceptible to external environmental factors such as vibrational noise and dust contamination. Traditional external sensors can be significantly affected by such factors, resulting in measurement errors and reduced reliability. In contrast, the proposed method provides a more reliable and precise measurement system that can be easily integrated into various motor configurations and withstand demanding environmental conditions. [Conclusions] To validate the proposed method, experiments were conducted on a self-built physical system that included the Real-Time eXtension (RTX) operating system, a Windows software platform, a magnetic sensor array, and an FPGA hardware platform. These experiments covered permanent magnet parameter offset calibration, stability tests, and comparative performance evaluations. The results from both simulations and experiments show that the proposed method can rapidly and accurately measure the pose of a three-degree-of-freedom rotary motor. The root-mean-square errors in orientation measurement for the three rotational axes were 0.031, 0.025, and 0.056 rad, respectively, with a resolution of approximately 1.77?0-4 rad. These results confirm the method's capability to provide high-precision pose measurements, positioning it as a promising solution for advanced motion control applications in multi-degree-of-freedom motors. Both simulation and experimental results validate the high-precision pose measurement capability of the proposed method, highlighting its potential for such advanced motion control applications.
  • XIE Jiaqi, ZHANG Han, ZHU Zhiming
    Journal of Tsinghua University(Science and Technology). 2024, 64(12): 2068-2083. https://doi.org/10.16511/j.cnki.qhdxxb.2024.22.043
    Abstract ( ) Download PDF ( ) HTML ( )   Knowledge map   Save
    [Significance] Rotary friction welding (RFW) is an extensively studied and applied solid-state welding method that can achieve high-quality welding between similar or dissimilar materials. RFW involves complex thermal, mechanical, and metallurgical processes, with heat generation occurring due to intense thermomechanical coupling. The stress and temperature histories at the friction interface considerably affect the microstructure evolution and mechanical properties of welded joints. Consequently, a primary focus of RFW research is deeply exploring the mechanisms and processes of thermomechanical coupling, obtaining the stress and temperature histories during welding to guide process-parameter selection and welded-joint-microstructure regulation. However, because of high-speed rotation and substantial plastic deformation during RFW, the temperature evolution and plastic deformation at the welding interface cannot be directly measured experimentally. Thus, mathematical models must be developed to study RFW. Currently, numerical simulation has become the predominant method for RFW theoretical research. With the development of computational technology, various numerical simulation methods have emerged, further elucidating the evolution laws of various physical fields during RFW and supporting theoretical research on the thermomechanical coupling behavior of RFW. [Progress] This paper reviews the research progress on the thermomechanical coupling behavior and numerical simulation technologies of RFW. It encompasses the theoretical underpinnings of the friction behavior of RFW, the development of heat-generation models, and the discussion of prevalent analytical and numerical methods for calculating temperature and stress fields during RFW. The proposal and research on RFW have a long history, resulting in the establishment of three friction-behavior theories: slide, stick, and slide-tick friction theories. These theories have informed the development of various thermomechanical coupling heat-generation models as well as material models. Analytical methods directly employ the thermomechanical coupling model to compute analytical solutions for the temperature and stress fields. These methods offer high computational efficiency and provide intuitive insights into heat generation and transfer processes, as well as material flow and deformation characteristics during RFW. However, analytical methods have challenges, such as their reliance on one-dimensional assumptions and simplified boundary conditions. Different stages of friction require distinct mathematical physics equations, complicating the achievement of coherent calculations for the entire welding process. Meanwhile, numerical simulation methods are more various, mainly including thermal conduction numerical models and the finite element method (FEM). As a mainstream numerical simulation method, the FEM can simulate material flow models and friction models, extending from two-dimensional to three-dimensional analyses. This allows for the obtainment of detailed information on temperature, stress and strain fields, residual stress distribution after welding, interface contact, and joint formation during RFW. In addition, the FEM can be effectively integrated with other simulation and prediction methods, such as microstructure evolution simulations and neural networks, offering comprehensive guidance for welded-joint-microstructure regulation and process-parameter selection. [Conclusions and Prospects] Presently, the simulation accuracy of RFW highly depends on material parameters and boundary condition settings, and the prediction capability of simulation models remains limited. Therefore, further research on welding mechanisms and the introduction of various computational methods to enhance the efficiency and accuracy of numerical simulation technologies while reducing computational costs represent current challenges and developmental directions for RFW simulations.
  • LI Zheng, ZHENG Shigang, DANG Xiaoyong, JI Wen, LIU Qu, CAI Zhipeng
    Journal of Tsinghua University(Science and Technology). 2024, 64(12): 2084-2091. https://doi.org/10.16511/j.cnki.qhdxxb.2024.22.030
    Abstract ( ) Download PDF ( ) HTML ( )   Knowledge map   Save
    [Objective] The primary objective of this research is to meticulously examine how pulse magnetic field assisted deep cryogenic (MDC) treatment affects the transformation and stabilization of retained austenite in Cr4Mo4V bearing steel. This study aims to elucidate the underlying mechanisms by which the pulse magnetic field influences the microstructural changes in bearing steel, particularly focusing on the stabilization of retained austenite, which plays a crucial role in determining the mechanical properties and overall performance of the steel. [Methods] To achieve a comprehensive understanding of how retained austenite transformed under various treatment conditions, this study utilized several material characterization techniques, including X-ray diffraction (XRD), vibrating sample magnetometry (VSM), and electron backscatter diffraction (EBSD). The use of EBSD analysis allows for a detailed comparison of variations in the dislocation density among samples processed under different conditions. For comparative analysis, the experimental set-up was divided into two distinct treatment processes: the conventional deep cryogenic (DC) treatment and the MDC treatment. Following these treatments, the samples were subjected to high-temperature tempering to evaluate the thermal stability of the retained austenite. [Results] The XRD analysis revealed a reduction in the volume fraction of retained austenite from (23.8%?.6)% to (21.5%?.9)% following the DC process. A relatively smaller reduction to (22.5%?.5)% was observed with the MDC process. These results, supported by VSM and EBSD analyses, highlight the capacity of the pulse magnetic field to partially inhibit the transformation of retained austenite. Further examination of the high-temperature stability of austenite in samples treated with DC and MDC revealed that MDC samples demonstrated improved retention, maintaining 7.1% of retained austenite after high-temperature tempering, compared to 4.9% in DC-treated samples. This indicates that the retained austenite in Cr4Mo4V bearing steel exhibits improved high-temperature stability following treatment with the MDC process. Furthermore, the dislocation density analysis revealed that the DC process led to a 9.8% increase in the dislocation density, whereas the MDC process moderated this increase to only 6.5%. This difference suggests the magnetic field's role in inhibiting dislocation diffusion, which in turn reduces martensite nucleation sites, thereby stabilizing retained austenite. The dislocation density change of the samples treated with DC and MDC after a high-temperature tempering validates this point. The dislocation density in DC-treated samples was approximately 1.23?015 m-2, while it decreased to 1.13?015 m-2 in MDC-treated samples. The dislocation density change reflects the extent of phase transformation. [Conclusions] This study provides a thorough analysis that clearly demonstrates the significant impact of applying a pulse magnetic field during deep cryogenic treatment on the microstructural evolution of Cr4Mo4V bearing steel. The magnetic field not only moderates the increase in the dislocation density but also enhances the mobility of dislocations. This contributes to the stabilization of retained austenite, which is crucial for improving the mechanical properties and performance of bearing steel. The findings of this research lay a solid foundation for optimizing heat treatment processes using the magnetic field assisted deep cryogenic treatment.
  • DAI Jingzhou, TIAN Ling, HAN Tianlin
    Journal of Tsinghua University(Science and Technology). 2024, 64(12): 2092-2104. https://doi.org/10.16511/j.cnki.qhdxxb.2024.27.025
    Abstract ( ) Download PDF ( ) HTML ( )   Knowledge map   Save
    [Objective] Sliding bearings are critical components of modern machinery, and their proper function is critical. However, inadequate lubrication can cause significant wear and degradation of the bearing contact surfaces, posing significant safety risks. Monitoring the wear state of sliding bearings during operation and predicting their remaining useful life (RUL) is crucial for ensuring equipment safety and reducing maintenance costs. Despite this need, the current online diagnostic and prognostic methods for sliding bearings are lacking. To address this issue, this study proposes an intelligent diagnostic and prognostic method for sliding bearing wear based on multidomain features and relevance-vector-based iterative exponential degradation (RV-IED). [Methods] This paper designed and built a test rig to simulate real-world operating conditions of sliding bearings and collected data under different wear conditions. Wear degradation and vibration measurement tests were conducted to measure the maximum wear depth (MWD) and vibration signals during the tests. For feature engineering, multidomain features combining the time domain, frequency domain, and nonlinear characteristics were constructed. From the energy of the vibration signals, time-domain features were derived. A Fourier transform was then applied to these signals to obtain frequency-domain waveforms, which are decomposed into multiple normal distributions. Calculating the intensity values of the top ten peaks' $\widetilde{\mu} \pm \widetilde{\sigma}$ provided a ten-dimensional frequency-domain vector, which was then reduced to one-dimensional frequency-domain features using principal component analysis. To handle the high-dimensionality of feature spaces, dynamic time warping was used to compute distances between different spectra as nonlinear features. These multidomain features served as input vectors for diagnosis, with the corresponding MWDs used as labels for training the relevance vector machine (RVM). New samples were diagnosed by outputting the current MWD. After each diagnosis, another RVM extracted sparse relevance vectors and corresponding weights from historical diagnosis data, fitting an exponential model using nonlinear least squares. This model predicts the bearing RUL by extending the trend to a preset threshold. [Results] The proposed method was evaluated against traditional time domain and frequency domain features combined with SVM/RVM methods as a control group. The experimental results showed that: (1) In diagnosing wear depth, the proposed method achieved a diagnostic error within 10%, outperforming the control group; (2) Although the diagnostic error increased as the training set size decreased, the changes were minimal beyond a reduction of 30%, making the method suitable for small sample sizes. We recommend a dataset size that does not exceed 1?04; (3) For RUL prediction, the proposed method's cumulative relative accuracy is 0.59, compared to 0.36 for the control group. [Conclusions] By leveraging the constructed sliding bearing wear test rig and monitoring data, multidomain features were created to accurately reflect bearing wear degradation. A diagnostic and prognostic method based on RV-IED for sliding bearing wear provides accurate diagnostics and RUL predictions, even for small sample sizes. This method surpasses traditional approaches and effectively supports intelligent diagnostics and predictive maintenance of sliding bearing wear states. This innovative approach holds promise for further advancement in fault prediction and health management in mechanical systems.
  • HAO Yuxing, WU Hongyu, ZHANG Yuling, WU Qingjian, TAN Lijie, YAN Shaoze
    Journal of Tsinghua University(Science and Technology). 2024, 64(12): 2105-2114. https://doi.org/10.16511/j.cnki.qhdxxb.2024.27.026
    Abstract ( ) Download PDF ( ) HTML ( )   Knowledge map   Save
    [Objective] Underwater glider is a low-energy and buoyancy-driven exploration robot that is used in long-term ocean exploration tasks, such as anticyclonic eddy observations, marine noise detection, and biogeochemical analyses. To enhance its overall performance, multidisciplinary optimization during the design and application phases is crucial. An effective method for enhancing hydrodynamic performance involves optimizing the surface micromorphology of underwater vehicles. Drawing inspiration from this micromorphology design concept, this paper conducts a multiobjective design optimization of the riblet parameters on the wing surface of an underwater glider to enhance its overall performance. [Methods] This study begins with the “Petrel-II” underwater glider prototype, proposing an initial riblet design scheme based on the movement characteristics of the underwater glider. Subsequently, a computational fluid dynamics simulation model of the entire glider was developed. An optimized Latin hypercube experimental design was used to obtain sample riblet parameters and attack angles, followed by computational fluid dynamics simulations under various operating conditions. Hydrodynamic equations, including riblet parameters, were established through polynomial fitting and computational fluid dynamics results. These hydrodynamic equations were integrated into the dynamic model of the glider, allowing for a comprehensive dynamic analysis that includes riblet effects. Using the entire vehicle dynamic model, a theoretical derivation of the underwater glider’s speed, static stability, and energy consumption evaluation was achieved, enabling a quantitative and comprehensive performance analysis. A surrogate model for performance evaluation was established using optimized Latin hypercube experimental designs, dynamic simulations, and performance models, considerably improving computational efficiency. The optimization objective function for the riblet parameters of the wing surface was determined. Riblet spacing, depth, and direction were selected as optimization design variables. The surrogate model and second-generation nondominated sorting genetic algorithm were used for optimization calculations. The optimization aimed to simultaneously improve the horizontal speed, vertical speed, energy efficiency, and static stability of the underwater glider. As a result, four single-objective optimal solutions and one multiobjective compromise solution were obtained, leading to the final riblet design scheme. [Results] Numerical examples demonstrate that the smaller direction and the larger depth of the riblet can effectively and simultaneously improve the glider’s speed, static stability, and energy efficiency. The effect of riblet spacing on the glider’s performance indicators shows certain contradictory and nonlinear characteristics. Compared with the performance analysis result of the glider without riblet on the wing surface, the optimization result considerably enhances the glider’s overall performance, proving the effectiveness of the proposed design method. In addition, 3D printing technology was used to produce a prototype of the wing. The surface of the obtained wing prototype is smooth, with the clear and undeformed riblet profile, which confirms the machinability of the riblet-enhanced surface. [Conclusions] The riblet on the wing surface can effectively improve the underwater glider’s overall performance, and the proposed optimization method can make the above improvement more significant. This research provides theoretical guidance and a reference for the actual optimization design of underwater gliders. The next step is to consider improving the micromorphology of the hull’s surface to further enhance the glider’s overall performance and carrying out the corresponding parameter optimization.
  • LU Yangping, TAN Lei
    Journal of Tsinghua University(Science and Technology). 2024, 64(12): 2115-2121. https://doi.org/10.16511/j.cnki.qhdxxb.2024.27.032
    Abstract ( ) Download PDF ( ) HTML ( )   Knowledge map   Save
    [Objective] Centrifugal pumps are a common type of hydraulic machinery widely used in water transfer, energy storage, agricultural irrigation, oil production, and more. Owing to system regulation demands, these pumps often need to switch operating conditions and, therefore, operate for long periods in off-design states. This requires a wide efficient operating range and the development of advanced design methods. To address the issues of low precision in load control design methods and narrow efficient operating ranges, this paper proposes a centrifugal pump impeller design method based on the three-segment quadratic function controlling blade load distribution. [Methods] The proposed method divides the blade load from the impeller inlet to the outlet into three segments: the inlet segment, the intermediate segment, and the outlet segment. Each segment is given a quadratic function-type load distribution to construct the impeller design method based on three-segment quadratic function load control. By assigning values to 13 independent control parameters, the specific form of load distribution is determined, resulting in the final structure of the impeller. Using this design method, a centrifugal pump with a specific speed of 102 is designed and numerically simulated to analyze its energy and flow-field characteristics. This study also establishes a comprehensive test bench for closed centrifugal pump hydraulic model testing, which consists of the test pump, pipeline system, and sensors that measure inlet pressure, outlet pressure, flow rate, rotational speed, and torque. The designed pump is manufactured and tested on this bench. [Results] The simulation results demonstrate that the optimized impeller achieves a 4.11% increase in efficiency at the designed point Qd and a 5.35% increase in the average efficiency under multiple operating conditions ranging from 0.7 Qd to 1.3 Qd, with an almost unchanged flow-head curve. The energy characteristics indicate that the optimized impeller has a wide efficient operating range. Internal flow-field analysis reveals a decrease in the area and strength of vortex and flow separation at off-design points. At the 50% span of the blades, the pressure distribution from the impeller inlet to the outlet is more uniform, and the streamline distribution is more reasonable. Furthermore, changing the blade load in the inlet section from a fast rise to a slow increase improves pressure distribution on the blade surface, reduces local pressure gradients, and achieves a uniform variation. The tested results show that the maximum efficiency of the designed pump is 78.20%, which meets the design requirements. Meanwhile, the simulated results are compared with experimental results. The error between the simulated and experimental results is lower than 5%, validating the accuracy of the simulation method. [Conclusions] Both simulated and experimental results confirm that the centrifugal pump impeller design method based on the three-segment quadratic function load control significantly enhances pump efficiency, improves flow patterns, and provides a reference for the development of high-efficiency and wide-operating-condition impellers.
  • HYDRAULIC ENGINEERING
  • LI Dan, QIN Chao, XUE Yuan, XU Mengzhen, LIU Yingjie
    Journal of Tsinghua University(Science and Technology). 2024, 64(12): 2122-2131. https://doi.org/10.16511/j.cnki.qhdxxb.2024.21.021
    Abstract ( ) Download PDF ( ) HTML ( )   Knowledge map   Save
    [Objective] Rivers carry water and material transport within certain boundary forms. Due to the difficulties of in-site measurement, limited hydrological observation stations, and the precision constraints of digital elevation model (DEM), there is a significant scarcity of information on small river cross-sections distributed in river source areas, mountainous regions, and remote areas, which hinders research on river hydrology and hydraulic processes. Since the beginning of this century, the International Association of Hydrological Sciences (IAHS) has been advocating for solutions to the challenges of hydrological prediction in ungauged basins (PUB). Despite large-scale remote sensing technology is increasingly applied to the extraction of river hydraulic parameters, the spatial resolution of satellite altimetry data is too low for small rivers (river width less that 150 m), which account for a high proportion of river networks. The National Aeronautics and Space Administration (NASA) launched the ICESat-2 (Ice, Cloud and land Elevation Satellite-2) satellite in 2018. This satellite was equipped with the Advanced Topographic Laser Altimeter System (ATLAS), a photon-counting LiDAR system for the first time. The light spot (footprint) it projects onto the Earth's surface has a diameter of about 17 m, with a center-to-center distance between spots of only 0.7 m. This allows for the acquisition of photon point cloud data with smaller spots and higher density along the track. The high-density photon point cloud provided by ICESat-2 offers the possibility of extracting hydrological parameters of narrow rivers with high precision. This study focuses on the small rivers with a river width of less than 10 m in the Huangfu River basin, a first-order tributary of the middle reaches of the Yellow River, which is a data-scarce region. [Methods] A method for extracting the cross-sectional morphology of small rivers using ICESat-2 ATL03 data has been proposed. First, photons with medium to high confidence levels are selected to eliminate most of the noise. Then, a smoothing filter is applied for precise de-noising. Finally, the point cloud that has been precisely de-noised is manually edited to generate a DEM. The cross-sectional morphology of three different locations of small rivers were extracted based on DEM and compared with unmanned aerial vehicle (UAV) in-situ measurement results. [Results] The results show that: (1) the method proposed in this study, which combines the selection of medium to high confidence point clouds with filtering denoising, can effectively remove the noise from photon point clouds, with a denoising rate above 63%; (2) the completeness and richness of ground points extracted based on ATL03 data are superior to those obtained by reclassifying ATL03 using ATL08 product; (3) the results of river cross-sections extracted based on ATL03 data are basically consistent with the UAV in-situ measurement results (R2>0.96, RMSE=0.69 m). [Conclusions] The research results preliminarily demonstrate the feasibility of using ICESat-2 altimetry data for extracting cross-sections of small rivers in data-scarce areas, partially supplying the three-dimensional spatiotemporal of small rivers in data-deficient areas, providing technical support for construction of three-dimensional river network across entire basins and the simulation of hydrological and hydraulic processes. This also indicates that ICESat-2 altimetry data have research prospects in obtaining hydrological parameters of rivers in data-scarce areas.
  • YANG Fang, HU Yuying, SONG Lixiang, ZHAO Jianshi
    Journal of Tsinghua University(Science and Technology). 2024, 64(12): 2132-2143. https://doi.org/10.16511/j.cnki.qhdxxb.2024.22.032
    Abstract ( ) Download PDF ( ) HTML ( )   Knowledge map   Save
    [Objective] The two-dimensional (2D) hydrodynamic model can simulate the process of flood inundation and evolution. This model is widely used in flood forecasting. The number of grids and time steps considerably affect the computational efficiency of this model. Various methods have been proposed to improve the computational efficiency of the 2D hydrodynamic model. Some researchers use the dynamic grid strategy to calculate only effective grid cells to reduce the impact of the increased number of grids. Others use the local time step (LTS) technology to decrease the time consumption caused by small time steps. Whether the efficiency of model computation can be improved by combining the advantages of the two strategies requires further research. [Methods] Herein, a hybrid algorithm that combines the dynamic grid strategy and LTS technology is proposed to further improve the model performance based on the self-developed flood simulation model, HydroMPM. Technically, the grid cells that actually assist in the flux calculation are first selected as effective cells. Then, the LTS technology is applied to these cells to further optimize the flux calculation and update strategy. The calculation accuracy and efficiency of the hybrid algorithm are compared and analyzed using an ideal dam break case and a typical flood simulation scenario in the Nangang River basin, Guangdong Province, China. [Results] The dynamic grid strategy can accelerate model computation by computing only the effective grid cells. However, the effective cells actually contain all computation grids when the computation area is completely submerged. In this case, the dynamic grid strategy may lead to a high computation amount and low computation efficiency owing to the dynamic update mechanism on all grid cells. The LTS technology can improve the average time step by hierarchically updating the grid cells. However, the performance of this technology barely depends on the difference in grid scale and flow velocity distribution. The urban flood process often has a scattered distribution of the local inundation area, which is suitable for the application of the dynamic grid strategy. At the same time, local mesh refinement is also required in urban flood simulation. This refinement enables the model to better describe the topographic variation in local waterlogging-vulnerable areas. However, it also leads to a large difference in the maximum time step allowed by different grid scales, which requires the application of the LTS technology. By combining the two strategies, the hybrid dynamic grid and LTS technology can further enhance the model performance. In the ideal dam break case, the hybrid algorithm can reduce computation time by 13.1%-64.8%. In the practical application case of the Nangang River basin, the calculation time can be saved by approximately 60%. [Conclusions] The hybrid algorithm successfully combines the advantages of the original dynamic grid strategy and LTS technology to further accelerate the efficiency of the model computation. However, the performance of this hybrid algorithm varies depending on the application scenarios. If the 2D hydrodynamic model is used to simulate only river floods, estuaries, and other areas, the effect of the dynamic grid strategy may not be obvious. The LTS technology may have a general effect when the grid distribution is uniform and the flow state is stable. The hybrid algorithm can combine the advantages of the two abovementioned algorithms. However, this hybrid algorithm has the inherent shortcomings of the two algorithms. In practical applications, a suitable algorithm should be chosen depending on the specific application scenario to achieve good simulation effect and computing efficiency.
  • ZHANG Hui, SHEN Xiaopeng, SHEN Chen, MO Weibin, HUYAN Bin, WANG Rui
    Journal of Tsinghua University(Science and Technology). 2024, 64(12): 2144-2154. https://doi.org/10.16511/j.cnki.qhdxxb.2024.22.035
    Abstract ( ) Download PDF ( ) HTML ( )   Knowledge map   Save
    [Objective] Vibration-sinking caissons are increasingly popular in coastal engineering because of their low cost, ease of installation, and high stability. However, verifying excitation force and vibration amplitude before installation is essential. Current specifications require extensive geological surveys and geotechnical tests, which are costly and time-consuming. A rapid and reliable assessment of vibration-sinking feasibility with simple on-site tests is essential, especially for large-diameter caissons. [Methods] This study proposes a method for quickly evaluating the vibration-sinking feasibility of steel cylindrical caissons based on standard penetration test (SPT) blow counts. The proposed method correlates soil shear strength with SPT blow counts through empirical relationships to calculate installation resistance during the vibrating penetration process of steel cylindrical caissons. The accuracy of the method is verified by comparing results with current specifications and on-site records. In addition, a quantitative analysis assesses the influence of soil properties and interlayers, represented by SPT blow counts, on the maximum installable diameter of the caisson. [Results] The calculated installation resistance at different depths using the proposed method generally aligns with current specifications. For a steel caisson with a 24-m diameter, the maximum installation resistance calculated during the vibration sinking matches the maximum vibratory force recorded on-site with an 8% error, demonstrating the reliability of the method. When the SPT blow counts are small, installation resistance in clay is greater than that in sand for the same SPT blow counts, whereas higher blow counts yield greater resistance in sand than in clay. For shallow installation depths, vibration amplitude requirements dominate the feasibility analysis, whereas excitation force requirements become primary for deeper sinking depths. The consideration of interlayers is important for feasibility analysis, with a critical number of SPT blow counts identified. The interlayer has little influence on the vibration-installation capacity unless the SPT blow counts exceed the critical number. Beyond this point, the vibration-installation capacity rapidly decreases as SPT blow counts of the interlayer increase. The critical number of SPT blow counts is related to the interlayer thickness and bottom position; deeper and thicker interlayers have a greater impact and lower critical number of SPT blow counts. The impact of the bottom position becomes considerably greater than the thickness when the blow counts exceed the critical number. [Conclusions] The proposed method provides a quick and reliable estimation of the vibration-sinking feasibility of steel caissons based on SPT blow counts. The method effectively predicts the maximum installable diameter of offshore cylindrical caisson foundations based on SPT blow counts and offers valuable insights into the effects of soil conditions and interlayers, thereby enhancing the feasibility and optimization of large-diameter steel caissons.
  • CIVIL ENGINEERING
  • ZHAO Huanshuai, PAN Yongtai, YU Chao, QIAO Xin, CAO Xingjian, NIU Xuechao
    Journal of Tsinghua University(Science and Technology). 2024, 64(12): 2155-2165. https://doi.org/10.16511/j.cnki.qhdxxb.2024.22.029
    Abstract ( ) Download PDF ( ) HTML ( )   Knowledge map   Save
    [Objective] In rock-crushing processes, external loading methods are important factors affecting the mechanical properties and fracture behavior of rocks. Among these loading methods, vibration and impact methods are the most common ones. However, previous research has mainly focused on macroscopic failure features and energy dissipation properties under the singular loading of vibration or impact. Research on the composite loading of vibration and impact is relatively scarce, and few studies have investigated the influence of vibration loading on the microscopic fracture characteristics and energy evolution during rock impacts. In particular, quantitative characterization studies are lacking. The research on the influence of vibration loading on the propagation of impact cracks and the energy utilization efficiency in rocks has significant academic and engineering applications to fully adapt to the needs of modern mine construction and high efficiency, energy saving, and green production. [Methods] The quasi-brittle green sandstone material, commonly used in rock-crushing operations, was taken as the research object. The macro/micromechanical response relationship of green sandstone was established by integrating indoor experiments with microscopic parameter calibration. The parallel bonding model was adopted, and two loading methods — impact and composite loading of vibration and impact — were compared and analyzed to investigate the influence of vibration loading on the propagation of impact cracks and the energy utilization efficiency in the failure process of green sandstone. The analysis was conducted using the particle flow code (PFC). [Results] The research results indicate that under the same impact velocity, increasing the frequency or amplitude of vibration leads to an increasing trend in the number of cracks in green sandstone. Under the two loading methods, the maximum number of cracks in green sandstone shows a nearly linear increase as the impact velocity increases, with the majority being tensile cracks. The distribution characteristics of cracks exhibit the X-shaped conjugate slope. However, the growth rate of cracks is relatively high under composite loading of vibration and impact. The quantitative characterization of the increase in the number of cracks and impact velocity under vibration loading is established. Under equivalent impact velocity, as the frequency and amplitude increase, there is a corresponding increase in both the proportion of vibration input energy and the energy utilization efficiency in green sandstone. However, as the impact velocity increases, the proportion of vibration input energy within the total input energy in green sandstone decreases. Concurrently, the maximum energy utilization efficiency shows a trend of rapid increase followed by a decrease, with the maximum increase reaching 0.725%. [Conclusions] In practical rock-crushing applications, appropriately increasing vibration loading can exacerbate the damage and deterioration of rocks. This process significantly enhances the energy utilization efficiency with lower energy input. This study preliminarily explores the impact of vibration loading on the propagation of impact cracks and the energy utilization efficiency in green sandstone to provide a reference for the rational selection of parameters in rock-crushing processes.
  • ZHOU Jiaxing, WANG Jin-an, LI Fei
    Journal of Tsinghua University(Science and Technology). 2024, 64(12): 2166-2176. https://doi.org/10.16511/j.cnki.qhdxxb.2024.21.018
    Abstract ( ) Download PDF ( ) HTML ( )   Knowledge map   Save
    [Objective] As the depth of coal mining increases, the activation of discontinuous structures, such as faults, poses a significant risk to the safe and efficient mining of coal seams. Therefore, acquiring precise knowledge of the distribution of in-situ stress is paramount for the design, construction, and disaster prevention of mining engineering. [Methods] This study proposes an inversion method for in-situ fields applicable to discontinuous zones of deep coal seams. (1) Given the discontinuity characteristics of the deep in-situ stress field, stability discriminants for normal faults and stability discriminant equations for positive faults, reverse faults, and strike-slip fault zones are derived based on the lateral pressure coefficients of in-situ stress. (2) A long short-term memory neural network algorithm is adopted to optimize the learning of the in-situ stress field data formed in different periods sequentially to effectively solve the nonlinearity, discreteness, and multi-noise problems of the measured deep in-situ stress data and to ensure that the excellent in-situ stress data information is remembered for a long time and that the inferior in-situ stress data information is forgotten in time. [Results] This study considers the main and auxiliary well areas of Yingjun's second mining area in Shanghai Miao as the research background and establishes an algorithm model for long short-term memory neural networks. Given the distribution characteristics of the in-situ stress field in fault areas at different scales, an inversion calculation of the in-situ stress field in discontinuous areas of deep coal seams was conducted. [Conclusions] The correlation coefficient between inverted and measured stress fields was 0.945, with an average error of 12.897%. The standard deviation of the stress difference is 2.000. The amount and direction of the regional stress field of the fault will also change. Compared with the regional stress field, the in-situ stress field in the DF15 and SF15 large-scale fault zones is approximately 5 MPa lower, and the counterclockwise direction is deflected. The in situ stress in the surrounding rock area at the top and bottom of the coal seams adjacent to DF15 and SF15 large-scale fault zones is relatively small, and no stress concentration area is detected. The eighth overlying coal seam was tilted toward horizontal in-situ stress by extrusion during deposition, and a concentration of in-situ stress was detected on the top rock. Therefore, for deep coal seam mining in the well field, the protective-layer mining method can be adopted, i.e., the lower 15 coals can be mined first to provide protection and pressure relief for the mining of the upper 8 coals to ensure the safety and reliability of deep-mining in the well field. Folds mainly control the distribution of the horizontal in-situ stress field, and the in-situ stress of the rock mass in the axial part of the backslope and the inner arc increases, while the in-situ stress of the outer arc of the obliquity is relatively small. Therefore, the inversion method proposed in this study can offer a new perspective for reconstructing the stress fields in deep discontinuous areas.
  • PUBLIC SAFETY
  • ZHANG Ying, Wu Lingli, Zhang Bo
    Journal of Tsinghua University(Science and Technology). 2024, 64(12): 2177-2184. https://doi.org/10.16511/j.cnki.qhdxxb.2024.27.020
    Abstract ( ) Download PDF ( ) HTML ( )   Knowledge map   Save
    [Objective] To investigate the impact of thermal aging on the fire propagation characteristics of conductors, the fire propagation characteristics of conductors under different equivalent aging ages were experimentally studied. The effects of different equivalent aging years on the fire propagation behavior of polyethylene wire, including the spread and dripping behavior, were studied. [Methods] First, the service life of the expected test wire in an actual scene was derived from the empirical formula of the aging model investigated in previous studies, and the required test conditions were accordingly calculated. [Results] The experimental results show that with the extension of thermal aging time, the flame height, width, and spread speed decrease. During the equivalent aging time, the frequency and mass of the droplets in the melt decreased with increasing equivalent aging time, whereas the mass of individual droplets increased with increasing equivalent aging time. According to the experimental results and theoretical analysis, a thermal aging law for the wire fire propagation characteristics was established. An analysis of the effect mechanism of thermal aging on fire spread behavior reveals that an increase in the equivalent aging age mainly affects the change in the elongation at break and melting point of the thermal insulation layer. With the increase in the equivalent aging time, the elongation at break of the polyethylene insulation layer decreases, which affects the increase in the surface tension of the layer. Under the influence of high temperature, the melting speed of the solid insulation layer was accelerated, and the layer can more easily transition into a liquid state. The liquid melt on the upper surface of the insulation layer flows down the surface under the action of gravity. Once the self-gravity of the accumulated molten material overcomes the constraint of its surface tension, individual droplets are formed for dripping, and this dynamic constraint process becomes longer as the equivalent aging age increases. This reduces the frequency of dripping molten droplets. With the accumulation of the melt on the surface of the insulation layer, the volume of the melt accumulation increases, which affects the increase in the heat loss of the melt in the combustion zone. The flame located above the insulation layer reduces the width, height, and width of the flame because of the loss of fuel on the surface of the insulation layer, thus affecting the reduction of the fire spread speed. The self-constructed fire spread and propagation model reveals that the change in flame width directly affects the convective heat transfer of flame to the insulation layer of the pyrolysis zone. Moreover, less heat is transferred via convection, which results in a slower spread of the fire. [Conclusions] Overall, this study revealed the internal mechanism of fire spread and dripping behavior of heat-aged wire, which can inspire engineers to develop more flame-retardant materials, promote the innovation of refractory materials, enhance the fire safety of various electrical equipment, and mitigate the harm caused by aging wire fires to human life and property.