[go: up one dir, main page]

 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (1,727)

Search Parameters:
Keywords = discrete element method

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
22 pages, 1885 KiB  
Article
A Study on the Bearing Performance of an RC Axial Compression Shear Wall Strengthened by a Replacement Method Using Local Reinforcement with an Unsupported Roof
by Yuanwen Liu, Jie Deng, Yigang Jia, Guangyu Wu, Naiwen Ke and Xianglan Wei
Buildings 2024, 14(9), 2926; https://doi.org/10.3390/buildings14092926 (registering DOI) - 15 Sep 2024
Abstract
When compared with conventional replacement reinforcement methods, the method of replacement using a local reinforcement with an unsupported roof has the advantages of shortening the reinforcement cycle and reducing material loss, and many scholars have carried out useful explorations thereof. At present, the [...] Read more.
When compared with conventional replacement reinforcement methods, the method of replacement using a local reinforcement with an unsupported roof has the advantages of shortening the reinforcement cycle and reducing material loss, and many scholars have carried out useful explorations thereof. At present, the formula for the bearing capacity of reinforcement by replacing concrete in the Code does not consider the effect of stress hysteresis on the parts of the reinforcement; so when the initial stress level is greater than 0.4, the Code’s strength utilization coefficient of 0.8 for the new concrete in the replaced area is on the side of insecurity. In this study, we are trying to improve and supplement the formula in the Code through the following work. Firstly, 18 groups of shear wall models were constructed using a VFEAP finite element analysis program to analyze the bearing performances of the shear walls after the replacement. The results showed that the replacement concrete strength, the initial stress level and the size of the replaced area had a significant influence on the bearing capacity of the shear wall after the replacement. Secondly, utilizing the replacement concrete strength, the initial stress level and the size of the replacement area as key parameters, then introducing the strength improvement coefficient considering the constraints of the stirrups, the modified strength utilization coefficient of new concrete in the replaced area was formulated. Finally, based on the modified strength utilization coefficient, the replacement bearing capacity formulas for the one-batch, two-batch, and three-batch replacements were derived, and an N-batch replacement bearing capacity formula was regressed and fitted on the basis of these equations, which are less discrete and more secure than the Code’s formula. Full article
(This article belongs to the Special Issue Advanced Studies on Steel Structures)
18 pages, 3521 KiB  
Article
Design Optimization and Experimentation of Triple-Head Gradually Reducing Spiral Precision Fertilizer Apparatus
by Guoqiang Dun, Quanbao Sheng, Haitian Sun, Xinxin Ji, Zhenzhen Yu, Hongxuan Wang, Xingpeng Wu, Yuhan Wei, Chaoxia Zhang, Shang Gao and Hailiang Li
Processes 2024, 12(9), 1988; https://doi.org/10.3390/pr12091988 (registering DOI) - 14 Sep 2024
Viewed by 208
Abstract
In order to solve the problems of existing spiral fertilizer apparatuses, such as the variation in cavity filling rate with rotational speed, fluctuation of fertilizer discharge flow, and inability to discharge fertilizer precisely, a triple-head gradually reducing spiral fertilizer apparatus is designed, which [...] Read more.
In order to solve the problems of existing spiral fertilizer apparatuses, such as the variation in cavity filling rate with rotational speed, fluctuation of fertilizer discharge flow, and inability to discharge fertilizer precisely, a triple-head gradually reducing spiral fertilizer apparatus is designed, which gradually compresses fertilizer particles through the triple-head reducing fertilizer spiral structure to achieve complete cavity filling and uniform fertilizer discharge. The main factors that affect the particle motion state and the structural size of the spiral fertilizer through theoretical analysis are determined, and its theoretical fertilizer discharge amount and rotational speed are calculated. Using EDEM (Discrete Element Method Software 2022) to establish a simulation model of a single-head gradually reducing fertilizer apparatus, the spiral lead reduction percentage x1, spiral diameter reduction percentage x2, and rotational speed x3 are determined as experimental factors, and the filling rate μ and spiral torque Yaverage are used as experimental indicators to conduct a simulation study on the secondary universal rotation combination design experiment. The results show that when the rotational speed is 95 r/min, the spiral lead reduction percentage is 60.00~73.21%, the spiral diameter reduction percentage is 86.55~97.05%, the filling rate μ is greater than 95%, and the spiral torque Yaverage is less than 16 N·m. In order to further improve the uniformity of fertilizer discharge and ensure the controllable accuracy of fertilizer discharge, comparative verification experiments are conducted on single-, double-, and triple-head gradually reducing spiral fertilizer discharge devices and ordinary spiral fertilizer discharge devices. The results show that the precision of the gradually reducing spiral fertilizer apparatus is better than that of the ordinary spiral fertilizer apparatus. Moreover, it is determined that the three-head style performed best. The triple-head gradually reducing spiral fertilizer apparatus is also validated by randomly adjusting six rotational speeds. The experiment results show that the average deviation of the fertilizer discharge flow rate of the fertilizer apparatus from the preset value is 3.16%. The two have a minor deviation, and the fertilizer precision is high. Precise control of the amount of fertilizer discharged can be achieved by adjusting the rotational speed, and the research can provide a specific reference for the improved design and precise control of the spiral fertilizer apparatus. Full article
(This article belongs to the Section Advanced Digital and Other Processes)
16 pages, 6104 KiB  
Article
Engineering Discrete Element Method-Based Design and Optimization of the Key Components of a Spoon-Wheel Spinach Seed-Metering Device
by Gang Zheng, Bing Qi, Wenyi Zhang, Weixing Shao, Lei Zhang, Yunxia Wang and Youqiang Ding
Agronomy 2024, 14(9), 2096; https://doi.org/10.3390/agronomy14092096 (registering DOI) - 14 Sep 2024
Viewed by 163
Abstract
In view of issues affecting manual on-demand sowing of spinach, such as high working intensity, low seeding efficiency, and high cost, a mechanized precision seeding device for spinach seeds is proposed, based on a spoon-wheel including a combined cell with two intersecting surfaces. [...] Read more.
In view of issues affecting manual on-demand sowing of spinach, such as high working intensity, low seeding efficiency, and high cost, a mechanized precision seeding device for spinach seeds is proposed, based on a spoon-wheel including a combined cell with two intersecting surfaces. The key structure was designed, and motion analysis was conducted employing SolidWorks software (2022). Additionally, the seeding process was simulated using EDEM simulation software (2022). The number of spoons, the radius of the spoons, and the speed of the seed wheel were chosen as the influencing factors, while the qualified index, reseeding index, and missing index were utilized as the evaluation indexes. A response surface method test based on the Box–Behnken design was carried out, and the test results were analyzed to obtain the optimal parameter combination for the seed-metering device. The field test results showed that when the rotation speed was 18 rpm, the radius of each seed spoon was 2.5 mm, and the number of seed spoons was 50, the average qualified index was 90.89%, the reseeding index was 8.22%, and the missed seeding index was 0.89%. The test results satisfy the technical production requirements for spinach seed sowing. Full article
Show Figures

Figure 1

Figure 1
<p>Structure diagram of seed-metering device: (<b>a</b>) assembled model; (<b>b</b>) exploded views: 1. shaft; 2. rear shell; 3. brush; 4. seed-carrying wheel; 5. adjustment baffle; 6. seed-picking wheel; 7. front shell.</p>
Full article ">Figure 2
<p>Work area division diagram of seed-metering device: I—seed picking; II—seed clearing; III—seed carrying; IV—seed dropping.</p>
Full article ">Figure 3
<p>Parameter of the hole and force analysis diagram of seeds.</p>
Full article ">Figure 4
<p>Simulation models (<b>a</b>). The model of seeds (<b>b</b>) The model of seeding device.</p>
Full article ">Figure 5
<p>Effect of test factors on the qualified index: (<b>a</b>) response surface diagram of interaction between the rotational speed and the radius of the spoons; (<b>b</b>) response surface diagram of interaction between the rotational speed and the number of the spoons; (<b>c</b>) response surface diagram of interaction between the radius and the number of spoons.</p>
Full article ">Figure 6
<p>Effect of test factors on the missing index: (<b>a</b>) response surface diagram of interaction between the rotational speed and the radius of the spoons; (<b>b</b>) response surface diagram of interaction between the rotational speed and the number of the spoons; (<b>c</b>) response surface diagram of interaction between the radius and the number of spoons.</p>
Full article ">Figure 7
<p>Bench test of the seed-metering device.</p>
Full article ">Figure 8
<p>Field testing of planter: (<b>A</b>) photographs of field trials; (<b>B</b>) field trial sowing effect.</p>
Full article ">
24 pages, 5903 KiB  
Article
Design and Testing of Soybean Double-Row Seed-Metering Device with Double-Beveled Seed Guide Groove
by Huajiang Zhu, Sihao Zhang, Wenjun Wang, Hongqian Lv, Yulong Chen, Long Zhou, Mingwei Li and Jinhui Zhao
Agriculture 2024, 14(9), 1595; https://doi.org/10.3390/agriculture14091595 - 13 Sep 2024
Viewed by 174
Abstract
During the operation of a shaped hole seed-metering device, poor seed-filling quality and inconsistent seed-casting points lead to poor seed spacing uniformity, especially in a one-chamber double-row seed-metering device. To solve this problem, a soybean double-row seed-metering device with double-beveled seed guide groove [...] Read more.
During the operation of a shaped hole seed-metering device, poor seed-filling quality and inconsistent seed-casting points lead to poor seed spacing uniformity, especially in a one-chamber double-row seed-metering device. To solve this problem, a soybean double-row seed-metering device with double-beveled seed guide groove was designed to ensure a high single-seed rate and seed-casting point consistency. Through the theoretical analysis of the working process of the seed-metering device, dynamic and kinematic models of the seeds were established, and the main structural parameters of the seed discharge ring, triage convex ridge, shaped hole, and seed guide groove were determined. The main factors affecting the seeding performance were obtained as the following: the inclination angle of the triage convex ridge, the radius of the shaped hole, and the depth of the seed guide groove. A single-factor test was carried out by discrete element simulation to obtain the inclination angle of the triage convex ridge α3 = 29°, the radius of the shaped hole r1 = 4.16–4.5 mm, and the depth of the seed guide groove l1 = 0.49–1.89 mm. A two-factor, five-level, second-order, orthogonal rotation combination test was conducted to further optimize the structural parameters of the seed-metering device. The two test factors were the radius of the shaped hole and the depth of the seed guide groove, and the evaluation indices were the qualified rate, replay rate, and missed seeding rate. The results showed that the optimal combinations of the structural parameters were the radius of the shaped hole r1 = 4.33 mm and the depth of the seed guide groove l1 = 1.20 mm. Subsequent bench testing demonstrated that the seed discharge’s qualified rate was above 94% at operating speeds of 6–10 km/h, and the seeding performance was stable. The final results of the soil trench test showed that the seed-metering device exhibited a qualified rate of 93.31%, replay rate of 2.04%, and missed seeding rate of 4.65% at an operating speed of 8 km/h. This research outcome may serve as a valuable reference and source of inspiration for the innovative design of precision seed-metering devices. Full article
Show Figures

Figure 1

Figure 1
<p>Diagram of the structural principles of a soybean double-row seed-metering device with double-beveled seed guide groove, where 1 is the front cover; 2 is the front seed chamber; 3 is the seed discharge ring; 4 is the rear seed chamber; 5 is the brush wheel; 6 is the drive disc; 7 is the seed guard; 8 is the rear shell; 9 is the drive shaft; and 10 is the seeding wheel.</p>
Full article ">Figure 2
<p>Schematic diagram of the division of the working process area of the seed-metering device, where I is seed-filling area; II is seed-clearing area; III is seed transport area; IV is seed-casting area.</p>
Full article ">Figure 3
<p>Schematic diagram of the structure of the seed discharge ring.</p>
Full article ">Figure 4
<p>Schematic diagram of the forces on the seeds during seed-filling.</p>
Full article ">Figure 5
<p>Schematic diagram of seed force and speed in the process of seed-casting.</p>
Full article ">Figure 6
<p>Schematic diagram of the structure of the triage convex ridge.</p>
Full article ">Figure 7
<p>Schematic diagram of the structure of the shaped hole.</p>
Full article ">Figure 8
<p>Schematic diagram of distribution distance of the shaped hole.</p>
Full article ">Figure 9
<p>Schematic diagram of the curve equation of the guide groove on the seed-metering ring, where 1 is the shaped hole; 2 is the seed guide groove; 3 is seeds in the process of moving; 4 is the position the soybean seed needs to reach at time <span class="html-italic">t</span><sub>4</sub>; the dotted line <span class="html-italic">A</span><sub>1</sub><span class="html-italic">A</span><sub>3</sub> segment is the absolute motion trajectory of the soybean seed; and the Arc <span class="html-italic">A</span><sub>2</sub><span class="html-italic">A</span><sub>3</sub> segment is the relative motion trajectory of the soybean seed.</p>
Full article ">Figure 10
<p>Seed-metering device model and discrete element simulation process.</p>
Full article ">Figure 11
<p>Bench validation test. (<b>a</b>) Bench test device; (<b>b</b>) Test process; (<b>c</b>) Test effect.</p>
Full article ">Figure 12
<p>Soil trench validation test. (<b>a</b>) Soil trench test device; (<b>b</b>) Test process; (<b>c</b>) Test effect. where the red circles in (<b>b</b>,<b>c</b>) indicate the position of the seeds after landing.</p>
Full article ">Figure 13
<p>Effect curve of the inclination angle of the triage convex ridge on the test index.</p>
Full article ">Figure 14
<p>Effect curve of the radius of the shaped hole on the test index. (<b>a</b>) Pre-test results; (<b>b</b>) Test results.</p>
Full article ">Figure 15
<p>Effect curve of the depth of the seed guide groove on the test index.</p>
Full article ">Figure 16
<p>Response surface diagram: (<b>a</b>) the radius of the shaped hole and the depth of the seed guide groove on the qualified rate; (<b>b</b>) the radius of the shaped hole and the depth of the seed guide groove on the replay rate; and (<b>c</b>) the radius of the shaped hole and the depth of the seed guide groove on the missed seeding rate.</p>
Full article ">Figure 17
<p>Bench test results.</p>
Full article ">
22 pages, 10121 KiB  
Article
Simulation Analysis and Test of a Cleaning Device for a Fresh-Peanut-Picking Combine Harvester Based on Computational Fluid Dynamics–Discrete Element Method Coupling
by Jie Ling, Man Gu, Weiwen Luo, Haiyang Shen, Zhichao Hu, Fengwei Gu, Feng Wu, Peng Zhang and Hongbo Xu
Agriculture 2024, 14(9), 1594; https://doi.org/10.3390/agriculture14091594 - 13 Sep 2024
Viewed by 175
Abstract
In order to solve the problems of high impurity rate and large loss rate in the whole fresh peanut harvesting and production process in hilly areas of China, the method of computational fluid dynamics (CFD) and discrete element (DEM) coupling is used to [...] Read more.
In order to solve the problems of high impurity rate and large loss rate in the whole fresh peanut harvesting and production process in hilly areas of China, the method of computational fluid dynamics (CFD) and discrete element (DEM) coupling is used to examine the gas–solid two-phase simulation of the cleaning device of the crawler fresh-peanut-picking combine harvester. In addition, the Box–Behnken test method is used to analyze the influence of different parameters on the impurity content and loss rate of materials in the cleaning process by taking the fan speed, feed amount, and air inlet angle as the test factors, and the optimal combination of working parameters is sought. The simulation results show that when the fan speed is 2905.07 r/min, the feed rate is 0.80 kg/s, the air inlet angle is 43.14°, the impurity content is 7.32%, and the loss rate is 4.78%. Compared with the simulation test results, the impurity content is increased by 0.68%, and the loss rate is increased by 1.24%, which verifies the reliability of the simulation model, and the research results provide some technical support for the improvement of the cleaning device in the later stages. Full article
(This article belongs to the Section Agricultural Technology)
Show Figures

Figure 1

Figure 1
<p>Crawler fresh peanut-picking combine: 1. chassis walking system; 2. dynamical system; 3. peanut collection device; 4. cleaning device; 5. inclined scraper conveyor; 6. transmission system; 7. control console; 8. peanut picking device; 9. picking platform.</p>
Full article ">Figure 2
<p>Cleaning device. 1. Miscellaneous outlets l; 2. Cleaning chamber 2; 3. Centrifugal fan; 4. Cleaning chamber 1.</p>
Full article ">Figure 3
<p>The working principle of the cleaning device.</p>
Full article ">Figure 4
<p>Centrifugal fan.</p>
Full article ">Figure 5
<p>Main and side views of fan blades.</p>
Full article ">Figure 6
<p>Material discrete element model.</p>
Full article ">Figure 7
<p>Cleaning device model.</p>
Full article ">Figure 8
<p>Meshing of cleaning devices.</p>
Full article ">Figure 9
<p>Simulation of material movement at different fan speeds: (<b>a</b>) 2500 r/min; (<b>b</b>) 2750 r/min; (<b>c</b>) 3000 r/min.</p>
Full article ">Figure 10
<p>Validation test.</p>
Full article ">Figure 11
<p>The influence of interaction factors on impurity rate. (<b>a</b>) The effect of fan speed and feeding amount on impurity rate. (<b>b</b>) The effect of fan speed and air inlet angle on impurity rate. (<b>c</b>) The effect of feeding amount and air inlet angle on impurity rate.</p>
Full article ">Figure 12
<p>The influence of interaction factors on loss rate. (<b>a</b>) The effect of fan speed and feeding amount on loss rate. (<b>b</b>) The effect of fan speed and air inlet angle on loss rate. (<b>c</b>) The effect of feeding amount and air inlet angle on loss rate.</p>
Full article ">Figure 13
<p>Analysis contour of the flow field inside the cleaning chamber. (<b>a</b>) Cleaning chamber trace distribution. (<b>b</b>) Cleaning chamber vector distribution. (<b>c</b>) Air velocity cloud picture in cleaning chamber 1. (<b>d</b>) Air velocity cloud picture in cleaning chamber 2.</p>
Full article ">Figure 13 Cont.
<p>Analysis contour of the flow field inside the cleaning chamber. (<b>a</b>) Cleaning chamber trace distribution. (<b>b</b>) Cleaning chamber vector distribution. (<b>c</b>) Air velocity cloud picture in cleaning chamber 1. (<b>d</b>) Air velocity cloud picture in cleaning chamber 2.</p>
Full article ">Figure 14
<p>Airflow velocity distribution at the observation surface of the cleaning chamber. (<b>a</b>) Transverse airflow velocity distribution. (<b>b</b>) Longitudinal airflow velocity distribution.</p>
Full article ">Figure 15
<p>Material cleaning process in the cleaning chamber. (<b>a</b>) t = 0.08 s. (<b>b</b>) t = 0.56 s. (<b>c</b>) t = 1.64 s. (<b>d</b>) t = 2.38 s.</p>
Full article ">Figure 16
<p>Field test.</p>
Full article ">
21 pages, 13910 KiB  
Article
Sensitivity Analysis on Influential Factors of Strain Rockburst in Deep Tunnel
by Jiheng Gu, Jiaqi Guo, Zihui Zhu, Feiyue Sun, Benguo He and Hengyuan Zhang
Buildings 2024, 14(9), 2886; https://doi.org/10.3390/buildings14092886 - 12 Sep 2024
Viewed by 270
Abstract
Strain rockburst is a severe failure phenomenon caused by the release of elastic strain energy in intact rocks under high-stress conditions. They frequently occur in deep tunnels, causing significant economic losses, casualties, and construction delays. Understanding the factors influencing this disaster is of [...] Read more.
Strain rockburst is a severe failure phenomenon caused by the release of elastic strain energy in intact rocks under high-stress conditions. They frequently occur in deep tunnels, causing significant economic losses, casualties, and construction delays. Understanding the factors influencing this disaster is of significance for tunnel construction. This paper first proposes a novel three-dimensional (3D) discrete element numerical analysis method for rockburst numerical analysis considering the full stress state energy based on the bonded block model and the mechanics, brittleness, integrity, and energy storage of the surrounding rock. This numerical method is first validated via laboratory tests and engineering-scale applications and then is applied to study the effects of compressive and tensile strengths of rock mass, tunnel depth, and lateral pressure coefficient on strain rockburst. Meanwhile, sensitivity analyses of these influencing factors are conducted using numerical results and systematic analysis methods, and the influence degree of each factor on the rockburst tendency is explored and ranked. The results reveal that laboratory tests and actual engineering conditions are consistent with numerical simulation results, which validates the rationality and applicability of the novel rockburst analysis method proposed in this paper. With the increase in compressive strength, the stress concentration degree, energy accumulation level, maximum stress difference, and maximum elastic strain energy within the rock mass all increase, leading to a stronger rockburst tendency. Tunnel depth and the lateral stress coefficient are positively correlated with rockburst tendency. As the lateral pressure coefficient and tunnel depth increase, rockburst tendency exponentially increases, while the maximum stress difference and maximum elastic strain energy within the rock mass also increase. The influence degree of each factor is ranked from highest to lowest as follows: tensile strength, lateral pressure coefficient, compressive strength, and tunnel depth. The research results provide theoretical support and technical guidance for the effective prediction, prevention, and control of rock burst disasters in deep tunnels. Full article
Show Figures

Figure 1

Figure 1
<p>Process of BBM-3DEC coupled numerical simulation method for strain rockburst.</p>
Full article ">Figure 2
<p>Comparison of laboratory tests and numerical simulation results on rockburst failure of rock samples. (<b>a</b>) laboratory rock sample structure and size; (<b>b</b>) high voltage servo true three-axis testing machine; (<b>c</b>) rock specimen; (<b>d</b>) high-speed camera system; (<b>e</b>) numerical experimental model (<b>f</b>) failure situation of laboratory test rock samples; (<b>g</b>) lateral failure situation of laboratory test rock samples (“张拉破坏” means tensile failure; “剪切破坏” means shear failure); (<b>h</b>) distribution of crack displacement (Unit: m); (<b>i</b>) lateral distribution of crack displacement (Unit: m); (<b>j</b>) distribution of displacement vector; (<b>k</b>) lateral distribution of displacement vector; (<b>l</b>) principal stress difference (Pa); (<b>m</b>) lateral principal stress difference (Pa); (<b>n</b>) elastic strain energy density (J/m<sup>3</sup>); (<b>o</b>) lateral elastic strain energy density (J/m<sup>3</sup>); (<b>p</b>) rockburst tendency <span class="html-italic">REC</span>; (<b>q</b>) lateral rockburst tendency <span class="html-italic">REC</span>.</p>
Full article ">Figure 3
<p>Comparison between engineering rockburst situation and numerical simulation results.</p>
Full article ">Figure 4
<p>Three-dimensional discrete element numerical model for strain-type rockburst: (<b>a</b>) model size and boundary conditions; (<b>b</b>) tunnel size and monitoring points.</p>
Full article ">Figure 5
<p>Multivariate information response and rockburst tendency of surrounding rock under different tunnel depths: (<b>a</b>) <span class="html-italic">σ</span><sub>c</sub> = 60 MPa, principal stress difference (Pa); (<b>b</b>) <span class="html-italic">σ</span><sub>c</sub> = 60 MPa, elastic strain energy (J/m<sup>3</sup>); (<b>c</b>) <span class="html-italic">σ</span><sub>c</sub> = 60 MPa, <span class="html-italic">REC</span>; (<b>d</b>) <span class="html-italic">σ</span><sub>c</sub> = 80 MPa, principal stress difference (Pa); (<b>e</b>) <span class="html-italic">σ</span><sub>c</sub> = 80 MPa, elastic strain energy (J/m<sup>3</sup>); (<b>f</b>) <span class="html-italic">σ</span><sub>c</sub> = 80 MPa, <span class="html-italic">REC</span>; (<b>g</b>) <span class="html-italic">σ</span><sub>c</sub> = 100 MPa, principal stress difference (Pa); (<b>h</b>) <span class="html-italic">σ</span><sub>c</sub> = 100 MPa, elastic strain energy (J/m<sup>3</sup>); (<b>i</b>) <span class="html-italic">σ</span><sub>c</sub> = 100 MPa, <span class="html-italic">REC</span>.</p>
Full article ">Figure 5 Cont.
<p>Multivariate information response and rockburst tendency of surrounding rock under different tunnel depths: (<b>a</b>) <span class="html-italic">σ</span><sub>c</sub> = 60 MPa, principal stress difference (Pa); (<b>b</b>) <span class="html-italic">σ</span><sub>c</sub> = 60 MPa, elastic strain energy (J/m<sup>3</sup>); (<b>c</b>) <span class="html-italic">σ</span><sub>c</sub> = 60 MPa, <span class="html-italic">REC</span>; (<b>d</b>) <span class="html-italic">σ</span><sub>c</sub> = 80 MPa, principal stress difference (Pa); (<b>e</b>) <span class="html-italic">σ</span><sub>c</sub> = 80 MPa, elastic strain energy (J/m<sup>3</sup>); (<b>f</b>) <span class="html-italic">σ</span><sub>c</sub> = 80 MPa, <span class="html-italic">REC</span>; (<b>g</b>) <span class="html-italic">σ</span><sub>c</sub> = 100 MPa, principal stress difference (Pa); (<b>h</b>) <span class="html-italic">σ</span><sub>c</sub> = 100 MPa, elastic strain energy (J/m<sup>3</sup>); (<b>i</b>) <span class="html-italic">σ</span><sub>c</sub> = 100 MPa, <span class="html-italic">REC</span>.</p>
Full article ">Figure 6
<p>Multivariate information response and rockburst tendency of surrounding rock under different tunnel depths: (<b>a</b>) <span class="html-italic">σ</span><sub>t</sub> = 4 MPa, principal stress difference (Pa); (<b>b</b>) <span class="html-italic">σ</span><sub>t</sub> = 4 MPa, elastic strain energy (J/m<sup>3</sup>); (<b>c</b>) <span class="html-italic">σ</span><sub>t</sub> = 4 MPa, <span class="html-italic">REC</span>; (<b>d</b>) <span class="html-italic">σ</span><sub>t</sub> = 8 MPa, principal stress difference (Pa); (<b>e</b>) <span class="html-italic">σ</span><sub>t</sub> = 8 MPa, elastic strain energy (J/m<sup>3</sup>); (<b>f</b>) <span class="html-italic">σ</span><sub>t</sub> = 8 MPa, <span class="html-italic">REC</span>; (<b>g</b>) <span class="html-italic">σ</span><sub>t</sub> = 12 MPa, principal stress difference (Pa); (<b>h</b>) <span class="html-italic">σ</span><sub>t</sub> = 12 MPa, elastic strain energy (J/m<sup>3</sup>); (<b>i</b>) <span class="html-italic">σ</span><sub>t</sub> = 12 MPa, <span class="html-italic">REC</span>.</p>
Full article ">Figure 6 Cont.
<p>Multivariate information response and rockburst tendency of surrounding rock under different tunnel depths: (<b>a</b>) <span class="html-italic">σ</span><sub>t</sub> = 4 MPa, principal stress difference (Pa); (<b>b</b>) <span class="html-italic">σ</span><sub>t</sub> = 4 MPa, elastic strain energy (J/m<sup>3</sup>); (<b>c</b>) <span class="html-italic">σ</span><sub>t</sub> = 4 MPa, <span class="html-italic">REC</span>; (<b>d</b>) <span class="html-italic">σ</span><sub>t</sub> = 8 MPa, principal stress difference (Pa); (<b>e</b>) <span class="html-italic">σ</span><sub>t</sub> = 8 MPa, elastic strain energy (J/m<sup>3</sup>); (<b>f</b>) <span class="html-italic">σ</span><sub>t</sub> = 8 MPa, <span class="html-italic">REC</span>; (<b>g</b>) <span class="html-italic">σ</span><sub>t</sub> = 12 MPa, principal stress difference (Pa); (<b>h</b>) <span class="html-italic">σ</span><sub>t</sub> = 12 MPa, elastic strain energy (J/m<sup>3</sup>); (<b>i</b>) <span class="html-italic">σ</span><sub>t</sub> = 12 MPa, <span class="html-italic">REC</span>.</p>
Full article ">Figure 7
<p>Multivariate information response and rockburst tendency of surrounding rock under different tunnel depth: (<b>a</b>) <span class="html-italic">D</span> = 500 m, principal stress difference (Pa); (<b>b</b>) <span class="html-italic">D</span> = 500 m, elastic strain energy (J/m<sup>3</sup>); (<b>c</b>) <span class="html-italic">D</span> = 500 m, <span class="html-italic">REC</span>; (<b>d</b>) <span class="html-italic">D</span> = 1500 m, principal stress difference (Pa); (<b>e</b>) <span class="html-italic">D</span> = 1500 m, elastic strain energy (J/m<sup>3</sup>); (<b>f</b>) <span class="html-italic">D</span> = 1500 m, <span class="html-italic">REC</span>; (<b>g</b>) <span class="html-italic">D</span> = 2500 m, principal stress difference (Pa); (<b>h</b>) <span class="html-italic">D</span> = 2500 m, elastic strain energy (J/m<sup>3</sup>); (<b>i</b>) <span class="html-italic">D</span> = 2500 m, <span class="html-italic">REC</span>.</p>
Full article ">Figure 8
<p>Multivariate information response and rockburst tendency of surrounding rock under different lateral pressure coefficient: (<b>a</b>) <span class="html-italic">λ</span> = 0.9, principal stress difference (Pa); (<b>b</b>) <span class="html-italic">λ</span> = 0.9, elastic strain energy (J/m<sup>3</sup>); (<b>c</b>) <span class="html-italic">λ</span> = 0.9, <span class="html-italic">REC</span>; (<b>d</b>) <span class="html-italic">λ</span> = 1.5, principal stress difference (Pa); (<b>e</b>) <span class="html-italic">λ</span> = 1.5, elastic strain energy (J/m<sup>3</sup>); (<b>f</b>) <span class="html-italic">λ</span> = 1.5, <span class="html-italic">REC</span>; (<b>g</b>) <span class="html-italic">λ</span> = 2.1, principal stress difference (Pa); (<b>h</b>) <span class="html-italic">λ</span> = 2.1, elastic strain energy (J/m<sup>3</sup>); (<b>i</b>) <span class="html-italic">λ</span> = 2.1, <span class="html-italic">REC</span>.</p>
Full article ">Figure 9
<p>Fitting curves of <span class="html-italic">REC</span> changes under different influencing factors: (<b>a</b>) compressive strength; (<b>b</b>) tensile strength; (<b>c</b>) tunnel depth; (<b>d</b>) coefficient of lateral pressure.</p>
Full article ">Figure 10
<p>Steps of grey relational analysis method.</p>
Full article ">Figure 11
<p>Correlation coefficient of rockburst influencing factors.</p>
Full article ">Figure 12
<p>Average correlation degree of influencing factors of strain rockburst.</p>
Full article ">
16 pages, 9808 KiB  
Article
Influence of Fine Content and Mean Diameter Ratio on the Minimum and Maximum Void Ratios of Sand–Fine Mixtures: A Discrete Element Method Study
by Huaqiao Zhong, Zhehao Zhu, Jiajin Zhao, Lanyi Wei, Yanyan Zhang, Jiayu Li, Jiajun Wang and Wenguo Yao
Buildings 2024, 14(9), 2877; https://doi.org/10.3390/buildings14092877 - 11 Sep 2024
Viewed by 260
Abstract
As urbanization accelerates and surface space becomes increasingly scarce, the development and utilization of urban underground space have become more critical. The sand–fine mixture soils commonly found in river-adjacent and coastal areas pose significant challenges to the design and construction of underground structures [...] Read more.
As urbanization accelerates and surface space becomes increasingly scarce, the development and utilization of urban underground space have become more critical. The sand–fine mixture soils commonly found in river-adjacent and coastal areas pose significant challenges to the design and construction of underground structures due to their unique mechanical properties. In soil mechanics, the minimum and maximum void ratios are crucial indicators for assessing soil compressibility, permeability, and shear strength. This study employed the discrete element method (DEM) to simulate the minimum and maximum void ratios of sand–fine mixtures under various conditions by setting six fine contents and three mean diameter ratios. The results indicate that as the fine content increases, these void ratios exhibit a trend of initially decreasing and then increasing, which can be effectively modelled using a single-parameter quadratic function. Additionally, the initial shear modulus was closely related to the uniformity of contact distribution at the microscopic level within the specimens. This study also introduced a dimensionless parameter that simultaneously described changes in contact distribution and initial shear modulus. Full article
Show Figures

Figure 1

Figure 1
<p>Grain-size distribution curve of the mixture with varying mean diameter ratios and fine contents.</p>
Full article ">Figure 2
<p>Determination of the minimum void ratio in DEM.</p>
Full article ">Figure 3
<p>Determination of the maximum void ratio in DEM.</p>
Full article ">Figure 4
<p>The triaxial shear test in DEM.</p>
Full article ">Figure 5
<p>Comparison of effective stress paths obtained from the laboratory test and the DEM simulation [<a href="#B38-buildings-14-02877" class="html-bibr">38</a>].</p>
Full article ">Figure 6
<p>The impact of mean diameter ratio and fine content on minimum void ratio.</p>
Full article ">Figure 7
<p>The correlation between the transitional minimum void ratio and mean diameter ratio.</p>
Full article ">Figure 8
<p>The functional relationship between the quadratic team ratio and mean diameter ratio.</p>
Full article ">Figure 9
<p>The impact of mean diameter ratio and fine content on maximum void ratio.</p>
Full article ">Figure 10
<p>Correlation between the maximum void ratio and minimum void ratio.</p>
Full article ">Figure 11
<p>The impact of mean diameter ratio and fine content on initial shear modulus.</p>
Full article ">Figure 12
<p>The coefficient of variation in contact force with varying fine contents.</p>
Full article ">Figure 13
<p>The variation in force chain structure with fine content.</p>
Full article ">Figure 14
<p>The variation in force chain structure with mean diameter ratio.</p>
Full article ">Figure 15
<p>Correlation between the initial shear modulus and the deviation of minimum void ratio.</p>
Full article ">
23 pages, 9821 KiB  
Article
Packing Characteristics and Heat Transfer Performance of Non-Spherical Particles for Concentrated Solar Power Applications
by Aidana Boribayeva, Xeniya Gvozdeva and Boris Golman
Energies 2024, 17(18), 4552; https://doi.org/10.3390/en17184552 - 11 Sep 2024
Viewed by 268
Abstract
Concentrated solar power (CSP) technology relies on thermal energy storage to extend operating hours, making the selection of heat storage media crucial for system efficiency. Bauxite powder, known for its availability and high-temperature stability, emerges as a potential alternative to conventional materials in [...] Read more.
Concentrated solar power (CSP) technology relies on thermal energy storage to extend operating hours, making the selection of heat storage media crucial for system efficiency. Bauxite powder, known for its availability and high-temperature stability, emerges as a potential alternative to conventional materials in CSP systems. This study employed the discrete element method to investigate the influence of particle shape on the packing and heat transfer characteristics of non-spherical particles. The research focused on assessing the impact of particle sphericity by comparing spherical particles with non-spherical shapes, including ellipsoids and cylinders, and exploring the effect of varying the aspect ratio (AR) of the cylindrical particles. Particle sphericity significantly influenced packing morphology, with the cylindrical particles exhibiting distinct structural patterns that were absent in the ellipsoidal particles, and strongly affected heat transfer, as observed in the average temperature variations within the packed bed over time. The cylinders with higher aspect ratios demonstrated enhanced heat transfer rates, driven by the increased contact area and coordination numbers, despite their predominant misalignment with the heat flux direction. These insights are valuable for optimizing thermal energy storage media in CSP systems. Full article
(This article belongs to the Special Issue Highly Efficient Thermal Energy Storage (TES) Technologies)
Show Figures

Figure 1

Figure 1
<p>Schematic of cylindrical particle packing used for modeling heat transfer in packed beds.</p>
Full article ">Figure 2
<p>Schematic representation of a five-layer packed bed.</p>
Full article ">Figure 3
<p>Illustration of (<b>a</b>) particle orientation angle and (<b>b</b>) cylindrical particle sides.</p>
Full article ">Figure 4
<p>Temperature distribution in packed beds of (<b>a</b>) spheres, (<b>b</b>) ellipsoids, (<b>c</b>) and cylinders at 1 s and 5 s.</p>
Full article ">Figure 5
<p>Temporal evolution of average layer temperature: (<b>a</b>) layer 1, (<b>b</b>) layer 2, and (<b>c</b>) layer 3 for spherical, ellipsoidal, and cylindrical particles.</p>
Full article ">Figure 6
<p>Total contact area in each layer for compacts of spherical, ellipsoidal, and cylindrical particles.</p>
Full article ">Figure 7
<p>The force chain of (<b>a</b>) spherical, (<b>b</b>) ellipsoidal, and (<b>c</b>) cylindrical particles.</p>
Full article ">Figure 8
<p>Distribution of particle orientation angle for cylindrical and ellipsoidal particles.</p>
Full article ">Figure 9
<p>(<b>a</b>) Coordination number distribution and (<b>b</b>) arithmetic average of coordination numbers for cylindrical particles of various aspect ratios.</p>
Full article ">Figure 10
<p>RDF of cylindrical particle compacts normalized by particle diameter.</p>
Full article ">Figure 11
<p>Illustration of contact types between cylinders: (<b>a</b>) face–band, (<b>b</b>) parallel band–band, (<b>c</b>) band–edge, (<b>d</b>) face–edge.</p>
Full article ">Figure 12
<p>Variation of local packing fraction in (<b>a</b>) <span class="html-italic">x</span>, (<b>b</b>) <span class="html-italic">y</span>, and (<b>c</b>) <span class="html-italic">z</span> planes.</p>
Full article ">Figure 13
<p>Temperature distribution in packed beds of (<b>a</b>) cylinders with all ARs in temperature range from 300 K to 450 K and (<b>b</b>) cylinders with AR = 2.5 and AR = 3.0 in temperature range from 300 K to 400 K.</p>
Full article ">Figure 14
<p>Temporal evolution of average layer temperature for cylindrical particles with various ARs: (<b>a</b>) layer 1, (<b>b</b>) layer 2, and (<b>c</b>) layer 3.</p>
Full article ">Figure 15
<p>Distribution of orientation angles for cylindrical particles with various ARs.</p>
Full article ">Figure 16
<p>Total contact area in each layer for compacts of cylindrical particles with various ARs.</p>
Full article ">Figure 17
<p>Distribution of contact types for cylindrical particles with various ARs.</p>
Full article ">Figure 18
<p>Comparative visualization of cylinder pairs with (<b>a</b>) band–band skewed contact, AR = 1.0; (<b>b</b>) band–band skewed contact, AR = 3.0; (<b>c</b>) band–edge contact, AR = 1.0; (<b>d</b>) band–edge contact, AR = 3.0.</p>
Full article ">
21 pages, 15945 KiB  
Article
Mechanisms of Proppant Transport in Rough Fractures of Offshore Unconventional Reservoirs: Shale and Tight Sandstone
by Biao Yin, Yishan Lou, Shanyong Liu and Peng Xu
J. Mar. Sci. Eng. 2024, 12(9), 1582; https://doi.org/10.3390/jmse12091582 - 7 Sep 2024
Viewed by 459
Abstract
After hydraulic fracturing, unconventional reservoirs frequently encounter challenges related to limited effective proppant support distance and suboptimal proppant placement. Due to the strong heterogeneity of offshore reservoirs, which causes varying fracture roughnesses depending on different lithologies, a systematic study of the relationship between [...] Read more.
After hydraulic fracturing, unconventional reservoirs frequently encounter challenges related to limited effective proppant support distance and suboptimal proppant placement. Due to the strong heterogeneity of offshore reservoirs, which causes varying fracture roughnesses depending on different lithologies, a systematic study of the relationship between roughness and proppant transport could optimize operational parameters. This study incorporates the box dimension method for fractal dimension analysis to quantify roughness in auto-correlated Gaussian distributed surfaces created by true triaxial tests. Combined with the numerical analysis of (computational fluid dynamics) CFD-DEM (discrete element method) for bidirectional coupling, the laws of proppant deposition and transport processes within fractures with different roughnesses are obtained through comparative verification simulations. The results show that for rougher fractures of shale, the proppants are transported farther, but at JRC_52, (joint roughness coefficient), where there may be plugging in curved areas, there is a risk of near-well blockages. Compared to the smooth model, fluctuations in JRC_28 (tight sandstone) drastically increase turbulent kinetic energy within the fracture, altering particle transport dynamics. Moreover, smaller proppants (d/w ≤ 0.3) exhibit better transport capacity due to gravity, but the conductivity of the proppant is limited when the particles are too small. A d/w of 0.4 is recommended to guarantee transport capacity and proppant efficiency near the well. Additionally, proppants injected sequentially from small to large in shale fractures offer optimal propping effects, and can take advantage of the better transport capacity of smaller proppants in rough fractures. The large proppant (d/w = 0.8) is primarily deposited by gravity and forms a sloping sand bed, which subsequently ensures the aperture of the fractures. This research provides a fresh perspective on the influence of fracture roughness on proppant transport in offshore unconventional reservoirs and offers valuable considerations for the order of proppant injection. Full article
(This article belongs to the Section Ocean Engineering)
Show Figures

Figure 1

Figure 1
<p>CFD-DEM coupling process (<b>a</b>); schematic diagram of CFD-DEM coupling [<a href="#B26-jmse-12-01582" class="html-bibr">26</a>] (<b>b</b>).</p>
Full article ">Figure 2
<p>True triaxial test system and post-fracturing results of rock sample #1.</p>
Full article ">Figure 3
<p>Pump–pressure curves of the rock sample #1 (<b>a</b>); pump–pressure curves and result of the rock sample #3 after fracturing (<b>b</b>).</p>
Full article ">Figure 4
<p>Scanned fracture surfaces of two different rocks (<b>a</b>); calculation of the fractal dimensions (<b>b</b>).</p>
Full article ">Figure 5
<p>Histograms of fracture aperture and corresponding probability density fitted with a Gaussian curve.</p>
Full article ">Figure 6
<p>Comparison of typical values and equation-predicted JRC values [<a href="#B31-jmse-12-01582" class="html-bibr">31</a>] (<b>a</b>); the box dimension method calculates the fractal dimension D (<b>b</b>).</p>
Full article ">Figure 7
<p>Post-fracturing roughness identification results for two different lithologies.</p>
Full article ">Figure 8
<p>Establishment of 3D surfaces with different JRC (slight aperture range variation and pronounced undulations on rougher surfaces).</p>
Full article ">Figure 9
<p>The process of establishing the rough channel model of JRC_21.</p>
Full article ">Figure 10
<p>Flow channel model with five different roughnesses.</p>
Full article ">Figure 11
<p>Simulation of boundary conditions during proppant transportation.</p>
Full article ">Figure 12
<p>Verification of grid independence.</p>
Full article ">Figure 13
<p>Average settling velocity of the proppant (<b>a</b>); comparison of the experiment and simulation (<b>b</b>).</p>
Full article ">Figure 14
<p>Distribution curves of wall shear and turbulent kinetic energy for two different roughnesses.</p>
Full article ">Figure 15
<p>Contours of the proppant volume fraction under different roughness conditions (<b>a</b>); shear–stress distribution on P90 for different roughnesses (<b>b</b>).</p>
Full article ">Figure 16
<p>Deposition of particles in the DEM at different moments for tight sandstone (JRC_28).</p>
Full article ">Figure 17
<p>Deposition of particles in the DEM at different moments for shale (JRC_52).</p>
Full article ">Figure 18
<p>Transport and deposition curves of particles with different roughnesses.</p>
Full article ">Figure 19
<p>Distribution of proppant on two surfaces at different times (<b>a</b>); relationship between the number of collisions and roughness (<b>b</b>).</p>
Full article ">Figure 20
<p>Volume fraction distribution of the proppant for different particle sizes.</p>
Full article ">Figure 21
<p>Schematic diagram of the arrangement of different particle sizes in rough fractures.</p>
Full article ">Figure 22
<p>Velocity contours of mixed particle size transport in rough fractures (<b>a</b>); location and distribution of mixed particle size in rough fractures (<b>b</b>).</p>
Full article ">
26 pages, 12522 KiB  
Article
A Vision–Language Model-Based Traffic Sign Detection Method for High-Resolution Drone Images: A Case Study in Guyuan, China
by Jianqun Yao, Jinming Li, Yuxuan Li, Mingzhu Zhang, Chen Zuo, Shi Dong and Zhe Dai
Sensors 2024, 24(17), 5800; https://doi.org/10.3390/s24175800 - 6 Sep 2024
Viewed by 327
Abstract
As a fundamental element of the transportation system, traffic signs are widely used to guide traffic behaviors. In recent years, drones have emerged as an important tool for monitoring the conditions of traffic signs. However, the existing image processing technique is heavily reliant [...] Read more.
As a fundamental element of the transportation system, traffic signs are widely used to guide traffic behaviors. In recent years, drones have emerged as an important tool for monitoring the conditions of traffic signs. However, the existing image processing technique is heavily reliant on image annotations. It is time consuming to build a high-quality dataset with diverse training images and human annotations. In this paper, we introduce the utilization of Vision–language Models (VLMs) in the traffic sign detection task. Without the need for discrete image labels, the rapid deployment is fulfilled by the multi-modal learning and large-scale pretrained networks. First, we compile a keyword dictionary to explain traffic signs. The Chinese national standard is used to suggest the shape and color information. Our program conducts Bootstrapping Language-image Pretraining v2 (BLIPv2) to translate representative images into text descriptions. Second, a Contrastive Language-image Pretraining (CLIP) framework is applied to characterize not only drone images but also text descriptions. Our method utilizes the pretrained encoder network to create visual features and word embeddings. Third, the category of each traffic sign is predicted according to the similarity between drone images and keywords. Cosine distance and softmax function are performed to calculate the class probability distribution. To evaluate the performance, we apply the proposed method in a practical application. The drone images captured from Guyuan, China, are employed to record the conditions of traffic signs. Further experiments include two widely used public datasets. The calculation results indicate that our vision–language model-based method has an acceptable prediction accuracy and low training cost. Full article
Show Figures

Figure 1

Figure 1
<p>Basic workflow of YOLO-based traffic sign recognition.</p>
Full article ">Figure 2
<p>Basic workflow of the vision–language model-based traffic sign detection approach. The Chinese national standard of traffic sign can be found in Reference [<a href="#B37-sensors-24-05800" class="html-bibr">37</a>].</p>
Full article ">Figure 3
<p>Six categories of Chinese traffic signs.</p>
Full article ">Figure 4
<p>Keyword dictionary constructed by Chinese national standards.</p>
Full article ">Figure 5
<p>Keyword dictionary constructed by representative images and BLIPv2 model. The mark * indicates that the frozen encoder networks are used by BLIPv2.</p>
Full article ">Figure 6
<p>High-dimensional feature extraction with CLIP network. The mark * indicates that we adopt the frozen encoder networks to realize feature extraction.</p>
Full article ">Figure 7
<p>Locations of Guyuan in China.</p>
Full article ">Figure 8
<p>Cosine distance between traffic sign images and keywords.</p>
Full article ">Figure 9
<p>Predicted probability distribution of traffic sign images.</p>
Full article ">Figure 10
<p>Classification accuracy of traffic sign images based on the proposed method.</p>
Full article ">Figure 11
<p>Time consumption of each component in the vision–language model-based method.</p>
Full article ">Figure 12
<p>Classification accuracy of traffic sign images based on the YOLO programs.</p>
Full article ">Figure 13
<p>Distributions of four evaluation metrics in the bootstrapping test.</p>
Full article ">Figure 14
<p>Example images in CCTSDB2021 and TT100K.</p>
Full article ">Figure 15
<p>Quantitative evaluation of the CCTSDB2021 dataset.</p>
Full article ">Figure 16
<p>Examples of the false detections in two public datasets.</p>
Full article ">
22 pages, 6901 KiB  
Article
Imaging Pressure Distribution in Geological Reservoirs from Surface Deformation Data
by Reza Abdollahi, Sirous Hosseinzadeh, Abbas Movassagh, Dane Kasperczyk and Manouchehr Haghighi
Sustainability 2024, 16(17), 7710; https://doi.org/10.3390/su16177710 - 5 Sep 2024
Viewed by 379
Abstract
Geological reservoirs are widely used for storing or disposing of various fluids and gases, including groundwater, wastewater, carbon dioxide, air, gas, and hydrogen. Monitoring these sites is essential due to the stored assets’ economic value and the disposed materials’ hazardous nature. Reservoir pressure [...] Read more.
Geological reservoirs are widely used for storing or disposing of various fluids and gases, including groundwater, wastewater, carbon dioxide, air, gas, and hydrogen. Monitoring these sites is essential due to the stored assets’ economic value and the disposed materials’ hazardous nature. Reservoir pressure monitoring is vital for ensuring operational success and detecting integrity issues, but it presents challenges due to the difficulty of obtaining comprehensive pressure distribution data. While direct pressure measurement methods are costly and localized, indirect techniques offer a viable alternative, such as inferring reservoir pressure from surface deformation data. This inversion approach integrates a forward model that links pressure distribution to deformation with an optimization algorithm to account for the ill-posed nature of the inversion. The application of forward models for predicting subsidence, uplift, and seismicity is well-established, but using deformation data for monitoring underground activity through inversion has yet to be explored. Previous studies have used various analytical, semi-analytical, and numerical models integrated with optimization tools to perform efficient inversions. However, analytical or semi-analytical solutions are impractical for complex reservoirs, and advanced numerical models are computationally expensive. These studies often rely on prior information, which may only sometimes be available, highlighting the need for innovative approaches. This study addresses these challenges by leveraging advanced numerical models and genetic algorithms to estimate pressure distribution from surface deformation data without needing prior information. The forward model is based on a discrete Green matrix constructed by integrating the finite element method with Python scripting. This matrix encapsulates the influence of reservoir properties and geometry on the displacement field, allowing for the rapid evaluation of displacement due to arbitrary pressure distributions. Precomputing Green’s matrix reduces computational load, making it feasible to apply advanced optimization methods like GA, which are effective for solving ill-posed problems with fewer observation points than unknown parameters. Testing on complex reservoir cases with synthetic data showed less than 5% error in predicted pressure distribution, demonstrating the approach’s reliability. Full article
Show Figures

Figure 1

Figure 1
<p>Inversion workflow for imaging reservoir pressure distribution from deformation data.</p>
Full article ">Figure 2
<p>Schematic of Green matrix construction from observation points and reservoir elements.</p>
Full article ">Figure 3
<p>Geomechanical properties and geometry of Geertsma analytical model.</p>
Full article ">Figure 4
<p>Geometry and boundary conditions of circular shape reservoir in an infinite media.</p>
Full article ">Figure 5
<p>Numerical result stability test vs. different mesh sizes and selecting the optimum size.</p>
Full article ">Figure 6
<p>Vertical (<b>a</b>) and radial (<b>b</b>) displacements resulting from a cylindrical constant pressure reservoir (the conventional numerical method).</p>
Full article ">Figure 7
<p>Vertical (<b>a</b>) and radial (<b>b</b>) displacements resulting from a cylindrical constant pressure reservoir (analytical method).</p>
Full article ">Figure 8
<p>Comparison of radial displacement at the surface resulting from a cylindrical constant pressure reservoir using different forward models.</p>
Full article ">Figure 9
<p>Comparison of vertical displacement at the surface resulting from a cylindrical constant pressure reservoir using different forward models.</p>
Full article ">Figure 10
<p>Arrangement of reservoir elements and pressure distribution throughout the cylindrical-shaped reservoir.</p>
Full article ">Figure 11
<p>Vertical (<b>a</b>) and radial (<b>b</b>) displacements resulting from cylindrical variable pressure reservoir (conventional numerical method).</p>
Full article ">Figure 12
<p>Comparison of radial displacement at the surface resulting from a cylindrical reservoir (distributed pore pressure) using conventional numerical model and discrete Green matrix.</p>
Full article ">Figure 13
<p>Comparison of vertical displacement at the surface resulting from a cylindrical reservoir (distributed pore pressure) using conventional numerical model and discrete Green matrix.</p>
Full article ">Figure 14
<p>Comparison of estimated pressure distribution (inversion) with real pressure distribution for cylindrical-shaped reservoir.</p>
Full article ">Figure 15
<p>Schematic of reservoir and overlaying layers in scenario 3.</p>
Full article ">Figure 16
<p>Vertical (<b>a</b>) and horizontal (<b>b</b>) displacements resulting from constant pressure irregular shape reservoir (conventional numerical method).</p>
Full article ">Figure 17
<p>Comparison of horizontal displacement at surface estimated from conventional forward model and discrete Green matrix.</p>
Full article ">Figure 18
<p>Comparison of vertical displacement at surface estimated from conventional forward model and discrete Green matrix.</p>
Full article ">Figure 19
<p>Arrangement of reservoir elements and pressure distribution throughout the irregular-shaped reservoir.</p>
Full article ">Figure 20
<p>Vertical (<b>a</b>) and horizontal (<b>b</b>) displacements resulting from distributed pressure in the irregular-shaped reservoir (conventional numerical method).</p>
Full article ">Figure 21
<p>Comparison of estimated pressure distribution (inversion) with real pressure distribution in irregular-shaped reservoir.</p>
Full article ">Figure 22
<p>Error distribution across reservoir elements for pore pressure prediction.</p>
Full article ">Figure 23
<p>Sensitivity analysis of inversion approach concerning geomechanical properties.</p>
Full article ">
17 pages, 16231 KiB  
Article
Probing Internal Damage in Grey Cast Iron Compression Based on Acoustic Emission and Particle Flow
by Zhen Li, Zhao Lei, Sheng Xu, Hengyang Sun, Bin Li and Zhizhong Qiao
Processes 2024, 12(9), 1893; https://doi.org/10.3390/pr12091893 - 4 Sep 2024
Viewed by 365
Abstract
Grey cast iron releases energy in the form of stress waves when damaged. To analyse the evolution of the physical and mechanical properties and acoustic emission characteristics of grey cast iron under uniaxial compression, acoustic emission signals were collected at different rates (0.5, [...] Read more.
Grey cast iron releases energy in the form of stress waves when damaged. To analyse the evolution of the physical and mechanical properties and acoustic emission characteristics of grey cast iron under uniaxial compression, acoustic emission signals were collected at different rates (0.5, 1, and 2 mm/s). Combined with load-time curves, damage modes were identified and classified using the parametric RA-AF correlation analysis method. The results indicate the loading rate effects on the strength, deformation, acoustic emission (AE), and energy evolution of grey cast iron specimens. The acoustic emission counts align with the engineering stress–strain response. To better illustrate the entire failure process of grey cast iron, from its internal microstructure to its macroscopic appearance, X-ray diffraction (XRD) and optical microscopy (OM) were employed for qualitative and quantitative analyses of the material’s internal microstructural characteristics. The equivalent crystal model of grey cast iron was constructed using a Particle Flow Software PFC2D 6.00.30 grain-based model (GBM) to simulate uniaxial compression acoustic emission tests. The calibration of fine parameters with indoor test results ensured good agreement with numerical simulation results. Acoustic emission dynamically monitors the compression process, while discrete element particle flow software further analyses the entire damage process from the inside to the outside. It provides a new research method and idea for the study of crack extension in some metal materials such as grey cast iron. Full article
(This article belongs to the Section Particle Processes)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Experimental process.</p>
Full article ">Figure 2
<p>(<b>a</b>) Microstructure XRD of grey cast iron; (<b>b</b>) the proportion of each organisation in the material; (<b>c</b>) grey cast iron metallographic micrograph.</p>
Full article ">Figure 3
<p>Schematic diagram of the behaviour and components of (<b>a</b>) PBM and (<b>b</b>) SJM.</p>
Full article ">Figure 4
<p>(<b>a</b>) Initial model; (<b>b</b>) the particles are divided into polygonal regions using the Voronoi algorithm; (<b>c</b>) particle-swarm-based grouping of material micro parameters, with OM physical map of the material on the far left and PFC comparison simulation in the middle.</p>
Full article ">Figure 5
<p>The engineering stress–strain curve of cast iron HT150 under a uniaxial compression test.</p>
Full article ">Figure 6
<p>Engineering stress–strain curves and AE ring-down count relationships at various com-pression rates: (<b>a</b>) 0.5 mm/s; (<b>b</b>) 1 mm/s; (<b>c</b>) 2 mm/s.</p>
Full article ">Figure 7
<p>Distribution of AF and RA correlations for cast iron in different stages (I, II, III and IV) at different loading rates: (<b>a</b>) 0.5 mm/s; (<b>b</b>) 1 mm/s; (<b>c</b>) 2 mm/s.</p>
Full article ">Figure 8
<p>RA-AF correlation density map with different loading rates: (<b>a</b>) 0.5 mm/s; (<b>b</b>) 1 mm/s; (<b>c</b>) 2 mm/s. (<b>d</b>) The proportion of cracks in specimens under different loading rates is different.</p>
Full article ">Figure 9
<p>Comparison of the results of laboratory test and numerical simulation: (<b>a</b>) stress–strain curve; (<b>b</b>) failure mode and crack distribution.</p>
Full article ">Figure 10
<p>Numerical simulations of acoustic emission counts for different phases (I, II, III and IV) at different rates: (<b>a</b>) 0.5; (<b>b</b>) 1; and (<b>c</b>) 2 mm/s.To study the spatio-temporal evolution of acoustic emissions and cracks, we divided the simulated stress–strain curve into four stages, corresponding to the four stages observed in the experiment, and recorded the evolution process at each stage. The acoustic emission counts and crack formation in the model were mainly caused by bond breakdown between the particles inside the crystal (see the red line segment in <a href="#processes-12-01893-f011" class="html-fig">Figure 11</a>). As seen in <a href="#processes-12-01893-f011" class="html-fig">Figure 11</a>, when the axial strain is 0.04%, cracks can be observed between particles within various grains. However, these cracks are widely dispersed and sparse, resulting in a low count of acoustic emissions. This simulates the situation when compressing the grey cast iron material in the experiment, where a small number of acoustic emission counts and tiny cracks are present at the beginning of loading due to the presence of microscopic cracks and pores within the material. Loading up to 0.08%, it can be observed from <a href="#processes-12-01893-f011" class="html-fig">Figure 11</a> that the cracks appear mainly within the graphite flake crystals, at which point the acoustic emission counts start to slowly rise and become more dense. This corresponds to the experimental observation that after a period of loading, the graphite phase within the cast iron material begins to be compressed, generating a large number of stress waves that are detected by the acoustic emission equipment. As loading continues in the model, reaching a strain of 0.12%, cracks and acoustic emissions still predominantly occur within the graphite flake crystals. However, cracks begin to appear within the boundaries between pearlite crystals and graphite flakes, indicating the emergence of crack bands. At this stage, the acoustic emission count reaches its peak. Upon reaching the maximum load during experimentation, grey cast iron experiences a proliferation of quasi-cleaved surfaces and cracks propagating along grain boundaries, leading to plastic deformation and the rupture of the surrounding matrix. Final phase, extensive cracking and acoustic emission were observed in all crystals within the model, resulting in the formation of a 45-degree tilt crack band. Experimentally, the material eventually developed macroscopic cracks with a similar inclination following the penetration of minute cracks.</p>
Full article ">Figure 11
<p>The spatio-temporal evolution of force chains, cracks, and acoustic emissions during the loading process is as follows: At 0.04% strain, cracks (indicated by red lines in the figure) are sparsely dispersed throughout the various crystals. At 0.08% strain, a large number of cracks appears within the graphite flake crystals, with minimal variation in the force chains observed in the model up to this point. By 0.12% strain, cracks emerge at the grain boundaries of pearlite and graphite flake crystals, indicating a trend towards crack band formation across the entire model, accompanied by a significant reduction in both strong and weak force chains. At 0.13% strain, clear crack bands appear throughout the model, with strong force chains concentrating in their vicinity.</p>
Full article ">Figure 11 Cont.
<p>The spatio-temporal evolution of force chains, cracks, and acoustic emissions during the loading process is as follows: At 0.04% strain, cracks (indicated by red lines in the figure) are sparsely dispersed throughout the various crystals. At 0.08% strain, a large number of cracks appears within the graphite flake crystals, with minimal variation in the force chains observed in the model up to this point. By 0.12% strain, cracks emerge at the grain boundaries of pearlite and graphite flake crystals, indicating a trend towards crack band formation across the entire model, accompanied by a significant reduction in both strong and weak force chains. At 0.13% strain, clear crack bands appear throughout the model, with strong force chains concentrating in their vicinity.</p>
Full article ">Figure 12
<p>Displacement field distribution after damage of a single cracked specimen under uniaxial compression: (<b>a</b>) tension; (<b>b</b>) shearing; (<b>c</b>) tension–shearing.</p>
Full article ">Figure 13
<p>Numerical simulation of shear and tensile crack trends.</p>
Full article ">Figure 14
<p>Comparison of energy trend change.</p>
Full article ">
19 pages, 20977 KiB  
Article
Research on the Collection Characteristics of a Hydraulic Collector for Seafloor Massive Sulfides
by Huan Dai and Yan Li
J. Mar. Sci. Eng. 2024, 12(9), 1534; https://doi.org/10.3390/jmse12091534 - 3 Sep 2024
Viewed by 352
Abstract
Ore collection is very important in deep-sea mining for seafloor massive sulfide (SMS). In view of the characteristics of SMS ores produced by mechanical crushing, which contain coarse particles and a wide particle size distribution, in-depth research on the collection process with a [...] Read more.
Ore collection is very important in deep-sea mining for seafloor massive sulfide (SMS). In view of the characteristics of SMS ores produced by mechanical crushing, which contain coarse particles and a wide particle size distribution, in-depth research on the collection process with a device combining a rotary crushing head and a flat suction mouth was conducted. In this paper, solid–liquid two-phase flow in the hydraulic collection process with a drum rotation is carried out using the computational fluid dynamics-discrete element method (CFD-DEM), and the flow field characteristics and particle motion characteristics are analyzed. The results indicate that particles with a maximum diameter of 20 mm can be effectively collected when the suction velocity is 3 m/s. The collection process of SMS mainly goes through three stages: particle disturbance start-up, partial particle influx, and stable collection. In addition, the appropriate drum speed facilitates the collection of SMS ore. Finally, the correctness of the numerical method was assessed using similarity experiments. This work can be used to guide the design of underwater mining equipment. Full article
(This article belongs to the Section Ocean Engineering)
Show Figures

Figure 1

Figure 1
<p>Configure diagram of a typical SMS mining system.</p>
Full article ">Figure 2
<p>Typical mining machines: (<b>a</b>) dredging collection; (<b>b</b>) 2-drum counter rotation collection [<a href="#B18-jmse-12-01534" class="html-bibr">18</a>]; (<b>c</b>) milling + suction collection [<a href="#B19-jmse-12-01534" class="html-bibr">19</a>]; (<b>d</b>) jet punching collection [<a href="#B17-jmse-12-01534" class="html-bibr">17</a>].</p>
Full article ">Figure 3
<p>Particle–particle collision model.</p>
Full article ">Figure 4
<p>CFD-DEM coupling schematic.</p>
Full article ">Figure 5
<p>Scheme model of the integrated collecting device.</p>
Full article ">Figure 6
<p>Main structural parameters of the suction head.</p>
Full article ">Figure 7
<p>(<b>a</b>) Basin model and (<b>b</b>) grid model.</p>
Full article ">Figure 8
<p>The formation process of the ore deposit in EDEM: (<b>a</b>) particle generation; (<b>b</b>) particle sedimentation; (<b>c</b>) cleaning.</p>
Full article ">Figure 9
<p>Detection settings in EDEM.</p>
Full article ">Figure 10
<p>Static pressure contour of the hydraulic collection.</p>
Full article ">Figure 11
<p>Velocity contour of the hydraulic collection.</p>
Full article ">Figure 12
<p>Track diagram of the selected particles: (<b>a</b>) start-up; (<b>b</b>) partial inflow of particles; (<b>c</b>) stable collection.</p>
Full article ">Figure 13
<p>Velocity vector diagram at different speeds: (<b>a</b>) 15 r/min; (<b>b</b>) 30 r/min; (<b>c</b>) 45 r/min; (<b>d</b>) 60 r/min.</p>
Full article ">Figure 14
<p>Instantaneous state of particle lifting at 10 s: (<b>a</b>)15 r/min; (<b>b</b>) 30 r/min; (<b>c</b>) 45 r/min; (<b>d</b>) 60 r/min.</p>
Full article ">Figure 14 Cont.
<p>Instantaneous state of particle lifting at 10 s: (<b>a</b>)15 r/min; (<b>b</b>) 30 r/min; (<b>c</b>) 45 r/min; (<b>d</b>) 60 r/min.</p>
Full article ">Figure 15
<p>Top view of the instantaneous state of particles at different drum speeds: (<b>a</b>) 15 r/min; (<b>b</b>) 30 r/min; (<b>c</b>) 45 r/min; (<b>d</b>) 60 r/min.</p>
Full article ">Figure 16
<p>Mass flow rate changes over time.</p>
Full article ">Figure 17
<p>Changes in mass flow rate with speed.</p>
Full article ">Figure 18
<p>Variation in outlet velocity with RPM.</p>
Full article ">Figure 19
<p>Experimental configuration.</p>
Full article ">Figure 20
<p>Experimental system: (<b>a</b>) experimental platform; (<b>b</b>) simulated ore.</p>
Full article ">Figure 21
<p>Experimental result: (<b>a</b>) the grooves formed on the particle bed; (<b>b</b>) collected particles (most of the water was removed).</p>
Full article ">
18 pages, 4887 KiB  
Article
Finite Element Modeling of the Dynamic Response of Plywood
by Arkadiusz Charuk, Katarzyna Gawdzińska and Paweł Dunaj
Materials 2024, 17(17), 4358; https://doi.org/10.3390/ma17174358 - 3 Sep 2024
Viewed by 279
Abstract
Modeling the dynamic properties of wood and wood-based composites is a challenging task due to naturally growing structure and moisture-dependent material properties. This paper presents the finite element modeling of plywood panels’ dynamic properties. Two panels differing in thickness were analyzed: (i) 18 [...] Read more.
Modeling the dynamic properties of wood and wood-based composites is a challenging task due to naturally growing structure and moisture-dependent material properties. This paper presents the finite element modeling of plywood panels’ dynamic properties. Two panels differing in thickness were analyzed: (i) 18 mm and (ii) 27 mm. The developed models consisted of individual layers of wood, which were discretized using three-dimensional finite elements formulated using an orthotropic material model. The models were subjected to an updating procedure based on experimentally determined frequency response functions. As a result of a model updating relative errors for natural frequencies obtained numerically and experimentally were not exceeding 2.0%, on average 0.7% for 18 mm thick panel and not exceeding 2.6%, on average 1.5% for 27 mm thick panel. To prove the utility of the method and at the same time to validate it, a model of a cabinet was built, which was then subjected to experimental verification. In this case, average relative differences for natural frequencies of 6.6% were obtained. Full article
Show Figures

Figure 1

Figure 1
<p>Analyzed plywood structure with sheets of veneers (<b>a</b>) and the wood structure (<b>b</b>).</p>
Full article ">Figure 2
<p>Finite element model of an 18 mm thick plywood panel.</p>
Full article ">Figure 3
<p>Modal analysis test stand with measurement points arrangement: schematic representation (<b>a</b>) and actual stand (<b>b</b>).</p>
Full article ">Figure 4
<p>Selected mode shapes comparison—18 mm poplar plywood.</p>
Full article ">Figure 5
<p>Comparison of calculated and experimentally determined accelerance functions for plywood 18 mm (<b>a</b>) and 27 mm (<b>b</b>).</p>
Full article ">Figure 6
<p>Schematic representation of P<sub>6</sub> parameter definition.</p>
Full article ">Figure 7
<p>Sensitivity analysis results for 18 mm (<b>a</b>) and for 27 mm (<b>b</b>) plywood panels.</p>
Full article ">Figure 8
<p>Mode shapes distortion for 18 mm thick (<b>a</b>) and 27 mm thick (<b>b</b>) plywood panels.</p>
Full article ">Figure 9
<p>Comparison of frequency response functions before and after the model updating for 18 mm (<b>a</b>) and 27 mm (<b>b</b>) thick plywood panel.</p>
Full article ">Figure 10
<p>Geometrical structure (<b>a</b>) and finite element model (<b>b</b>) of cabinet.</p>
Full article ">Figure 11
<p>Comparison of frequency response functions for cabinet model MES and experimentally determined.</p>
Full article ">
20 pages, 7686 KiB  
Review
A Review of Discrete Element Method Applications in Soil–Plant Interactions: Challenges and Opportunities
by Yuyuan Tian, Zhiwei Zeng and Yuan Xing
Agriculture 2024, 14(9), 1486; https://doi.org/10.3390/agriculture14091486 - 1 Sep 2024
Viewed by 604
Abstract
The discrete-element method (DEM) has become a pivotal tool for investigating soil–plant interactions in agricultural and environmental engineering. This review examines recent advancements in DEM applications, focusing on both the challenges and opportunities that shape future research in this field. This paper first [...] Read more.
The discrete-element method (DEM) has become a pivotal tool for investigating soil–plant interactions in agricultural and environmental engineering. This review examines recent advancements in DEM applications, focusing on both the challenges and opportunities that shape future research in this field. This paper first explores the effectiveness of DEM in simulating soil and plant materials, including seeds, roots, and residues, highlighting its role in understanding interactions that affect agricultural practices. Challenges such as long computation times and the complexity of determining accurate contact parameters are discussed, alongside emerging methods like machine learning that offer potential solutions. Notable advancements include the integration of machine learning algorithms for contact parameter estimation, the use of expanded particle models for dynamic processes, and the development of new techniques for detailed post-processing of DEM simulations. The review also identifies key future research directions, including the incorporation of environmental factors such as air and water, and the exploration of residue management for carbon storage and erosion prevention. By addressing these challenges and seizing these opportunities, future research can enhance the accuracy and applicability of DEM models, advancing our understanding of soil–plant interactions and contributing to more sustainable agricultural and environmental practices. Full article
(This article belongs to the Section Agricultural Technology)
Show Figures

Figure 1

Figure 1
<p>Common soil–plant interaction systems simulated by discrete-element method (DEM).</p>
Full article ">
Back to TopTop