Next Article in Journal
Mango Fruit Fly Trap Detection Using Different Wireless Communications
Next Article in Special Issue
Weed Identification in Maize Fields Based on Improved Swin-Unet
Previous Article in Journal
The Effects of Different Sowing Density and Nitrogen Topdressing on Wheat Were Investigated under the Cultivation Mode of Hole Sowing
Previous Article in Special Issue
Performance of an Automatic Variable-Rate Fertilization System Subject to Different Initial Field Water Conditions and Fertilizer Doses in Paddy Fields
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Measurement and CFD-DEM Simulation of Suspension Velocity of Peanut and Clay-Heavy Soil at Harvest Time

1
Nanjing Institute of Agricultural Mechanization, Ministry of Agriculture and Rural Affairs, Nanjing 210014, China
2
Graduate School of Chinese Academy of Agricultural Sciences, Beijing 100083, China
*
Authors to whom correspondence should be addressed.
Submission received: 6 May 2023 / Revised: 23 June 2023 / Accepted: 26 June 2023 / Published: 28 June 2023
(This article belongs to the Special Issue Precision Operation Technology and Intelligent Equipment in Farmland)

Abstract

:
The suspension velocity is the core of the cleaning and sorting mechanisms that utilize a combination of a fan and vibrating sieve. To investigate this, various experimental subjects, such as peanuts with different kernels and clay-heavy clods in different states, were used. The experiment involved simulating the suspension velocity of materials through numerical calculations using fluid dynamics and particle discrete element coupling. The Eularian model was employed to study the coupled gas-solid two-phase flow. The experiment measured the suspension velocities of single and double kernel peanuts, which were found to be 8.34~9.40 m/s and 8.13~9.51 m/s, respectively. Under 20.4% water content and lumpy conditions, the suspension velocities of smaller clods, side by side clods, and larger clods were 12.61~14.30 m/s, 14.16~15.76 m/s and 16.44~18.72 m/s, respectively; under 20.4% water content and smaller clods, the suspension velocities of lumpy and strip of clods were 12.61~14.30 m/s, 11.90~14.13 m/s, respectively; under lumpy and smaller clods, the suspension velocity at 17.6%, 20.4%, and 23.9% water content ranged from 12.38 to 14.20 m/s, 12.61 to 14.30 m/s, and 12.62 to 14.49 m/s, respectively. The simulations showed that the suspension velocity for different types of peanuts, clod sizes, shapes, and water contents was less different from the actual experiments. Specifically, the relative errors in suspension velocity for single-kernel peanuts, double-kernel peanuts, smaller clods, side-by-side clods, larger clods, lumpy clods, strips of clods, and clods with 17.3%, 20.4%, and 23.9% water content were 1.2%, 4.1%, 0.4%, 2.0%, 4.4%, 0.4%, 5.1%, 5.4%, 0.4%, and 1.9%, respectively, compared to actual experiment measurements. The results indicate a significant difference in the suspension velocity between peanuts and clay-heavy clods, which can be distinguished from each other based on this difference. Furthermore, the simulation results have been found to be consistent with the experimental results, thus verifying the feasibility of measuring the material suspension velocity using CFD-DEM gas-solid coupling.

1. Introduction

Currently, both domestic and foreign peanut harvesting machinery utilize cleaning and sorting mechanisms that combine a fan and vibrating sieve. This principle is based on the difference in suspension velocity caused by varying material densities, allowing for the selection of appropriate air speeds to achieve material separation [1,2]. In addition, pneumatic conveying is also widely used in mechanical devices such as sowing, fertilizing, picking up, and harvesting [3,4,5,6]. As can be seen, the suspension velocity of the material is an important parameter in the design of agricultural machinery, such as scavenging and conveying devices. China’s peanut planting area ranks second in the world. Nearly one-third is planted in hilly mountainous areas with mostly clay-heavy soil [7]. Due to its small fields, sloping arable land, complex and diverse planting agronomy, and other multiple factors, peanut production mechanization and industrial development is severely restricted. Current peanut combine harvesting machines are designed for sandy loam soil planting, making them poorly adaptable for direct application to clay-heavy soil planting. This results in low harvesting efficiency, significant impurity rates, and high loss rates. These defects highlight the need for improved peanut harvesting machinery specifically designed for clay soil planting. Therefore, there is an urgent need for the research and development of combined harvesting machinery specifically designed for peanut cultivation on clay-heavy soils. This development is of great significance as it can solve the problems of low harvesting efficiency and increase farmers’ income. Furthermore, cleaning is a crucial step in peanut combined harvesting machinery that directly affects the rate of impurities and losses. Therefore, the separation and clearing of peanuts grown on clay-heavy soils for combined harvesting is of great importance to fill the gaps in the industry, while suspension velocity plays a central role in this process.
So far, many scholars have conducted suspension velocity measurement experiments and numerical simulation analyses related to different crop particles and debris. B.Y. Gorial et al. [8] numerically solved a mathematical model of particle motion to simulate the trajectory of wheat and straw in an air stream, suggesting the feasibility of achieving their separation based on particle suspension velocity. El-Sayed A.S. et al. [9] studied the physical and aerodynamic characteristics of different peanut varieties (American, Chinese, and Egyptian) and derived optimum airflow velocities for peanut shell separation, broken seed kernels, and intact seed kernels by means of suspension velocity experiments. Dai Fei et al. [10] measured the suspension velocities of flax grains, capsules, short stems, and capsule husks at different stages to derive the overall interval distribution and separation characteristics of each component of the flax threshing material, and used the gas-solid coupling method to simulate the suspension velocities of each component of the flax threshing material to verify the feasibility of measuring suspension velocities based on CFD-DEM simulation experiments. Yaquan Liang et al. [11] used a DFPF-25 material suspension velocity measurement device to measure the suspension velocity of each component of rice material, and verified through fluid simulation of the cleaning and sorting chamber that a dual fan cleaning structure consisting of a cross-flow fan and a centrifugal fan can significantly improve the cleaning efficiency. Lianxing Gao et al. [12] conducted aerodynamic experiments on the main components of a peanut hulling machine and impurities, and obtained the floating coefficients of broken peanut pods, unpurified peanut pods requiring secondary shelling, peanut rice, and peanut shells and stones, which provided a reference basis for the design of peanut hulling machine sorting devices and pneumatic conveying devices. Wang Yifei [13] measured the suspension velocity of peanut threshing material, including pods, long stalks, short stalks, stalks, and miscellaneous leaves, and used coupled CFD-DEM to simulate the material transport process and flow field simulation for the cleaning and sorting device of the peanut combine harvester.
Existing studies have shown the value of the coupled CFD-DEM-based theory as a guide and practical application for the analysis of the aerodynamic characteristics, study separation, and clearing of peanuts grown in clay-heavy soils and side-by-side clods, which can save costs and improve research efficiency.

2. Materials and Methods

2.1. Structure and Working Principle of the Levitation Wind Tunnel

The system composition of the device for measuring the suspension velocity of peanuts and clay-heavy clods is shown in Figure 1 and mainly consists of a holding net, an inlet, a stable section, a test section, a damping net, a buzzer, a motor, and a frequency converter. The working wind speed of the wind tunnel is adjusted by controlling the speed of the motor through the frequency converter.
The main technical parameters of the peanut clay suspension velocity measuring device are shown in Table 1.
The velocity measurement system consists of a differential pressure transmitter, an industrial control computer, and velocity measurement using the drop method, that is, by measuring the static pressure difference between the test section of the wind tunnel and the upstream stable section according to the fluid dynamics “Bernoulli” theorem to convert it to the test section’s working wind velocity. The differential pressure is first converted by the differential pressure transmitter into a voltage value, then the computer collects the voltage value as a digital quantity, which is calculated into the wind tunnel working wind speed value (unit: m/s) by the program cured in the computer and displayed to the digital tube in real time.
To ensure the accuracy of the measurement results, corrections should be made after the measurements are completed compared to the theoretically calculated values. First, set the material in the conical observation tube at a relatively stable suspension height, read the height range value and inverter frequency, and close the power supply to remove the material. Then, re-open the power supply and allow the fan speed to reach the same value through the high-precision rotating blade type anemometer, and read the airflow speed at the small end of the conical observation tube, according to Formula (1), to calculate the suspension velocity of the material [14].
V h = V 0 r r + 2 H tan θ 2
where, V h is the airflow velocity at any section of the conical observation tube, m/s; V 0 is the air velocity at the small end of the conical observation tube, m/s; r is the diameter of the small end of the conical observation tube, m; H is the suspension height of the material in the conical tube, m; θ is the angle of inclination of the conical observation tube °.

2.2. Experimental Materials and Suspension Velocity Measurement Methods

The experiment material consisted of Wanhua 17 peanuts and clay-heavy clods, all taken from the semi-fed peanut combine harvester in the field after cleaning and sorting into the fruit collection box section, as shown in Figure 2. The water content of the peanuts was measured at 24.7% and the clods at 20.4%. After sorting, the average geometry of the single-kernel peanuts was 25.04 × 16.95 × 16.11 mm, with an average mass of 1.92 g. The average geometry of the double-kernel peanuts was 39.24 × 16.17 × 17.02 mm, with an average mass of 3.06 g. Ten random samples of each were taken.
The clay-heavy clods were divided into three size ranges based on the actual harvest: smaller clods, side-by-side clods, and larger clods. The size of the side-by-side clods is essentially the same as the peanut pods, with a mass of 7.71 ± 0.50 g. The smaller clods are 1/2 the size of the side-by-side clods and are the smallest clods that can remain on the sieve surface under the action of the vibrating sieve, with a mass of 3.80 ± 0.50 g. The larger clods are three times the size of the smaller clods, with a mass of 11.74 ± 0.50 g. The shapes of the clods are divided into two types: lumps and strips. The soil water content is divided into three classes: 17.6%, 20.4%, and 23.9%. The water content of the soils in this sample was 20.4%, which was scaled appropriately to consider factors such as clay-heavy soils and harvesting weather in different areas. Ten random samples were taken from each.
The velocity of suspension of each classified material was tested 10 times, with each test repeated 3 times in a single-factor test. In the test, to measure the suspension velocity of different masses of clods, the shape of the clods was taken to be lumpy based on the number and mass share of the two shapes in the fruiting soil mixture, and a water content of 20.4% was taken. When conducting experiments to measure the suspension velocity of different shapes of clods, the smaller clods were taken based on the number and mass share of clods in the three size zones and the ease of performing clod removal. Then, a water content of 20.4% was taken. When carrying out experiments to measure the suspension velocity of clods with different water contents, smaller and lumpy clods were selected.
The experimental tools include height measuring scales, the HHF142 high-precision rotary vane anemometer (measuring range 0.2 to 40 MPS), vernier calipers, and electronic measuring scales, etc.
First, turn on the wind tunnel power supply, warm up for half an hour, and then place the material to be tested on the holding grid from the inlet, adjust the wind tunnel wind speed by inverter, and observe the suspension of the material in the conical tube (Figure 3). As the materials to be tested are irregularly shaped and unevenly distributed in mass, a fixed and precise suspension height value will not be obtained when the suspension velocity is measured. When the material suspension height is stabilized in a certain height range in the conical observation tube, stop pressing the wind speed adjustment button on the inverter panel and record the suspension height range [H1, H2] and the frequency of the inverter. Calculate the average suspension height range for the three replicate tests [h1, h2], and calculate the material suspension velocity range according to Equation (1). After recording, turn off the power to the wind tunnel, remove the material, turn the power back on, and adjust the fan speed to the recorded value. The air velocity at the small end of the conical tube is measured with the “OMEGA” HHF142 high-precision rotary vane anemometer and the arithmetic average of six measurement points is taken. The six measuring points are determined in accordance with the equal area method for circular duct sections, i.e., the small end of the conical observation tube section is divided into three concentric circles of equal area, and then each circle is divided into two parts of equal area, and the measuring points are chosen on the dividing line of these two parts [15].

3. Numerical Simulation Modelling

3.1. Controlling Equations for the Fluid and Particle Phases

In the CFD-DEM coupling process, the fluid phase is controlled by the Navier–Stokes equations for incompressible fluids:
α ρ t + Δ · α ρ u = 0
where, ρ is the is the density of the gas, kg/m3; α   is the volume fraction of the gas; u is the rate of gas flow, m/s; t is the time, s; Δ is the Hamiltonian differential operator.
The equation of conservation of momentum is
α ρ u t + Δ · α ρ u 2 = α Δ τ + α ρ g α Δ p i = 1 N F i V c
where, p is the is pressure on the micro-element, Pa; g   is the acceleration of gravity, m/s2; τ is the gas kinetic viscosity, Pa · s; F i is the resistance of particle i to the gas, N; V c is the volume of the cell grid, m3.
The particle phase control equation is as follows:
m d v p d t = m g + F p
d ω p d t = M p I p
where, m is the mass of the particle, kg; v p   is the velocity of the particle, m/s; F p is the combined force on the particle other than gravity, such as buoyancy, drag and additional mass force, N; ω P is the angular velocity of the particle, rad/s; M p is the combined moment of rotation of the particles, I p is the rotational inertia of the particle.

3.2. Material Model and Parameter Setting

Based on the actual distribution of fruit and debris in the fruit collection box after completion of harvesting by the peanut combine in hilly mountainous areas, the simulation experiments mainly analyzed the suspension state of different kernels of peanuts and clay-heavy clods with different masses and shapes as well as different water contents. The experiments were based on the discrete element simulation software EDEM 2021.2 with the multi-sphere filling method to create parametric models of differently shaped particle entities (shown in Figure 4) [16], analyze particle movement of the peanut pods and clay-heavy clods, and analyze the fluid dynamic properties based on Ansys 2022 R1 [17].
The physical and relevant action parameters of the particle model, derived from measurements, and a review of the relevant literature, are shown in Table 2 [18,19].

3.3. CFD-DEM Numerical Simulation Experiments Based on DDPM Model

The structural model of the suspension velocity experimental device is created using the 3D solid modeling software SpaceClaim 2022R1. To improve the quality of the mesh and reduce the solution time (the computer’s solution time increases exponentially with the number of mesh cells), the model should be simplified as much as possible. The wind tunnel model is extracted from the fluid domain, grouped, and named boundaries are created and saved in .scdoc format before being imported into Fluent Meshing for meshing. The minimum sized cell is 1.4 mm, and the maximum sized cell is 10 mm. The hexahedron dominant method is used to fill the body mesh to improve the calculation accuracy. A total of 89,999 cells are generated after the division is completed, with a maximum skewness of 0.48 and a mesh quality of 0.2, achieving the desired effect. The 3D model and the divided body mesh are shown in Figure 5. Save as a .msh file and import into EDEM and Fluent Solution, respectively, to prepare the coupling simulation.
After setting the parameters and importing the model in EDEM, a virtual plane of the dynamic particle plant was created with a diameter of 80 mm, slightly smaller than the diameter of the small end of the conical tube. The plane is located 30 mm from the bottom of the small end of the conical tube to ensure that the particles are completely in the fluid region when they are generated. Setting the particle generation after 0.01 s allows the wind field to be in a stable phase when the particles appear, limiting the number of generated particles to 1. Set the Rayleigh time step to 3 × 10−5 s [20]. Turn on the coupling switch after setting.
Import the mesh file in Fluent Solution with the solver as pressure basis and time as transient. After setting up the gravity field, import the coupling interface file to connect Fluent and EDEM. Activate the dense discrete phase model (DDPM) in the Eulerian model, which is solved using a multiphase flow framework with the fluid equations incorporating a volume fraction term [21]. Set up the discrete phase, define the jet source and phase, and define the inlet air flow rate in the boundary conditions and the outlet as a pressure outlet. After initialization, enter the solution interface and set the time step to 1.5 × 10−3 s, which should be no less than the Rayleigh time step in EDEM and should be an integer multiple. Then, set the number of time steps to 10,000 and the simulation time to 15 s.
In the calculation process, the setting of the air flow rate has a significant effect on the suspension of the material. A low air flow rate will prevent the material from overcoming gravity and will thus rest on the holding grid, while a high flow rate will cause the material to fly out of the outlet. It is therefore necessary to carry out pre-calculations with reference to the theoretically calculated values of the suspension velocity of each material to find the boundary conditions that allow the material to achieve the desired suspension effect.
The simulations show that when the inlet wind velocity is set at 17.86 m/s, 24.38 m/s, 26.76 m/s, and 30.78 m/s for the peanut pods, smaller clods, side by side clods, and larger clods, each material can reach a stable suspension and is at one third of the conical observation tube. In the experiment, it was found that each material was relatively close to the wall of the conical observation tube when it reached a stable suspension state under the action of the vertical wind field, and the smaller the mass of the material, the closer it was to the wall of the tube, as shown in Figure 6. The reason for this phenomenon is that the velocity of the airflow in any cross-section of the conical observation tube decreases from the center of the tube to the wall, and the pressure difference on the surface of the material particles causes the particles to be pushed upwards by the airflow while also moving in the direction of the tube wall. When the velocity of the airflow near the pipe wall is less than the velocity required to overcome the gravitational work of the particles, the particles begin to fall, creating a horizontal acceleration from the pipe wall towards the center of the pipe, moving towards the center of the pipe. When the particle falls to a position where the velocity of the airflow is greater than the velocity required to overcome the work of gravity, it rises again, and so on, finally reaching a relatively stable state of suspension near the wall of the tube. The smaller the mass of the material found in the experiments, the closer it is to the wall of the pipe because the suspension velocity is significantly influenced by the mass factor; the smaller the mass of the material suspension velocity, the more likely it will reach a stable suspension state in the cross-section of the pipe where the airflow is smaller, that is, near the wall of the pipe. If the mass of the experimental material is further reduced to a certain level, it will cling to the wall of the tube to achieve a stable suspension. The study by Xiangyu Wen et al. [22] also verifies the conjecture of this study.

4. Results and Analysis

4.1. Experimental Results and Analysis

4.1.1. Suspension Velocity of Peanuts of Different Kernels

The suspension velocity measurements of peanuts are shown in Table 3, where the peanuts have a water content of 25% and the frequency converter is 28 Hz, and x ¯ is the arithmetic mean.
As can be seen from the results of the experiment in Table 2 and the comparative graphs of suspension velocities of peanuts of different kernels in Figure 7a, the suspension velocities of single-kernel peanuts ranged from 8.34 to 9.40 m/s and double-kernel peanuts from 8.13 to 9.51 m/s, with little difference between the two. The overall data show that the suspension velocities of the double-kernel peanut are even slightly lower than that of the single-kernel peanut, which is related to the difference in the windward area of the wind tunnel due to the difference in shape and size between the two. The single-kernel peanut pods have a small difference in length, width, and height dimensions, with a shape close to an ellipsoidal shape, relatively uniform forces, and a relatively small floating range in the wind field when they reach a steady state of suspension. On the other hand, the double-kernel peanut pods have a large area of force in the vertical airflow despite their large mass and shape dimensions, which results in the suspension velocities of the double-kernel peanut pods being slightly less than that of the single-kernel peanut pods. In the practical experiments, it was evident that, due to the heterogeneous shape and size of the double-kernel peanut pods, although the stability of the suspended state was better than that of the single-kernel peanut pods, their up and down range in the wind tunnel was significantly greater. The independent sample t-test yielded a significance level of 0.768, which is greater than 0.05, showing that the effect of single and double kernel peanuts on their suspension velocity did not reach a significant state. Thus, overall, it can be concluded that the number of kernels in the peanut pods has a small effect on their suspension velocity.

4.1.2. Suspension Velocity of Different Masses of Clods

The clods were taken in lump form and had a water content of 20.4%. The values of suspension velocity for different masses of clay-heavy clods are shown in Table 4. During the experiment, the wind velocity was adjusted to different levels to achieve a relatively stable condition for all three size ranges of clay-heavy clods and to be located at 1/3 of the conical observation tube. The frequency converter is 38 Hz for smaller clods, 42 Hz for side-by-side clods, and 48 Hz for larger clods.
As can be seen from Table 3, the suspension velocity of smaller clods ranged from 12.61 to 14.30 m/s, side-by-side clods from 14.16 to 15.76 m/s, and larger clods from 16.44 to 18.52 m/s. It is clear that the magnitude of suspension velocity is directly related to the mass of the clods. The one-sample K-S test yielded an asymptotic significance level p = 0.93 > 0.05, meaning that the sample data obeyed a normal distribution. The chi-square test yielded a significance level based on the mean p = 0.29 > 0.05, i.e., the variance is chi-square. A one-way ANOVA test showed that the significance level of the effect of the clods on the suspension velocity was less than 0.001 for all three size intervals, reaching a highly significant level. The level of significant difference in suspension velocity between the smaller clods and the side-by-side clods was 0.001, and the level of significant difference between the smaller clods and the larger clods, and the side-by-side clods and the larger clods was less than 0.001, with significant differences in suspension velocity between the three sizes. From the box line comparison plot in Figure 7b, it can be seen that the up and down range of the clods in the wind field for the three size ranges tends to be essentially the same. Therefore, it is considered that the mass of the clods is positively correlated with the suspension velocity values.

4.1.3. Suspension Velocity of Different Shapes of Clods

Smaller clods were taken for the experiment, with a water content of 20.4%. The values of suspension velocity for different shapes of clods are shown in Table 5. The frequency converter is 38 Hz.
The data show that the suspension velocities of 12.95~14.97 m/s for lumpy clods and 12.71~14.75 m/s for strip of clods are not very different. As can be seen from the box line comparison plot in Figure 7c, the median line of suspension velocity for the two shapes of clods is essentially at the same level. In practical trials, it was found that the suspension height range of the clods was significantly lower than that of the strips, similar to the suspension of single and double kernel peanuts. The reason for this is that the windward area of the strips in the vertical airflow is large and stable (this phenomenon is also illustrated by the fact that the value of the suspension velocity of the strips in the box line diagram is smaller than that of the lumpy), but the windward area decreases abruptly when the airflow force is encountered at both ends during the turning process, resulting in a greater height difference between the up and down floating of the strips in the conical observation tube than that of the lumpy. The independent samples T-test yielded a significance level of 0.165, which is greater than 0.05. It is therefore concluded that the shape of the clods has less influence on their suspension velocity values, given the same water content and mass size.

4.1.4. Suspension Velocity of Clods at Different Water Contents

Smaller, lumpy clods were taken for the experiment. The suspension velocity values of clods with different water contents are shown in Table 6, where the frequency converter frequency is 38 Hz.
As can be seen from Table 6, the suspension velocity ranged from 12.38 to 14.20 m/s at 17.6% water content, 12.61 to 14.30 m/s at 20.4%, and 12.61 to 14.49 m/s at 23.9%. As can be seen from the box line comparison plot in Figure 7d, the suspension velocity values and median line for the three water content clods show a slight increasing trend. The one-sample K-S test yielded an asymptotic significance p = 0.165 > 0.05, meaning that the sample data obeyed a normal distribution. The chi-square test yielded a significance based on the mean p = 0.519 > 0.05, i.e., the variance is chi-square. A one-way ANOVA test yielded a significance level of 0.005, which is less than 0.05, suggesting that water content has a significant effect on suspension velocity. It is therefore believed that the soil water content affects its suspension velocity value, i.e., the suspension velocity value increases slowly with increasing water content.
In the experiments, it was found that when the water content reached 23.9%, the soil sample was very easy to adhere to the inner wall of the conical observation tube, which greatly increased the difficulty of the measurement. This is because too much water content in the soil particles will form a water film, which has a hydrogel effect and will make the adhesion between the soil and the tube wall increase [23]. Therefore, harvesting should be carried out as much as possible when the soil water content is below 20%. This not only reduces the difficulty of the operation but also ensures the effectiveness of soil removal and reduces the rate of contamination.

4.2. Numerical Simulation Results and Analysis

4.2.1. Analysis of Simulation Results

Once the calculations were completed, the flow field data calculated by Fluent 2022R1 software and the particle data calculated by EDEM 2021.2 software were imported into the post-processing software EnSight 2022R1 for coupled gas-solid flow analysis when each material was in steady suspension in the conical observation tube [24]. The suspension position and velocity of each material in stable suspension in the conical observation tube and the corresponding velocity distribution of the flow field are shown in Figure 8.
As can be seen from Figure 8, the particles within the flow field are blue in color, which indicates that they are all in a stable suspension. The corresponding airflow velocities were obtained from the coordinates of the particles’ positions within the flow field. Suspension velocities were 8.52~9.43 m/s for single-kernel peanuts and 7.74~9.20 m/s for double-kernel peanuts. Under 20.4% water content and lumpy conditions, suspension velocity was 12.18~14.62 m/s for smaller clods, 14.30~16.23 m/s for side-by-side clods, and 16.99~19.80 m/s for larger clods; under 20.4% water content and smaller clods conditions, suspension velocity was 12.18~14.62 m/s for lumpy clods and 12.67~14.76 m/s for strips; under the conditions of lumpy and smaller clods, the suspension velocity was 11.36~13.87 m/s at 17.6% water content, 12.18~14.62 m/s at 20.4%, and 12.52~14.09 m/s at 23.9%.
A comparison of the simulation experiment with the actual experimental results is shown in Figure 9. As the suspension velocity of each material is a range value, the relative error is calculated by taking its central value. The calculation formula is:
ϵ = v s v e v e × 100 %
where, v s   is the center value of the simulated experiment levitation velocity, m/s; v e   is the center of the actual experiments’ levitation velocity, and m/s; ϵ is the relative error.
The relative errors between the simulation and the actual experimental results were calculated to be 1.2%, 4.1%, 0.4%, 2.0%, 4.4%, 0.4%, 5.1%, 5.4%, 0.4%, and 1.9% for single-kernel peanuts, double-kernel peanuts, smaller clods, side-by-side clods, larger clods, lumpy clods, strip of clods, and a water content of 17.6%, 20.4%, and 23.9% respectively. The simulation results are less inaccurate than the actual experimental results, and it can be concluded that the simulation based on the coupled gas-solid two-phase flow method is highly accurate in calculating the suspension velocity of peanuts and clay-heavy clods.

4.2.2. Validation of the Model

In order to further determine the accuracy of the coupled model simulation to calculate the suspension velocity of peanut pods and clay-heavy clods, the material suspension velocity value was used as the evaluation index by increasing the air flow rate at the entrance of the flow field so that the peanut pods and clay-heavy clods reached a stable suspension state at different height positions in the conical observation tube and comparing whether the suspension velocity at this time is consistent with the actual experiment and the previously simulated suspension velocity value. The validation results are shown in Figure 10.
As can be seen from Figure 10, the suspension height of single and double kernel peanuts and clods of different mass size, shape, and water content in the conical observation tube increased with increasing inlet air flow rate, but their suspension velocity values remained stable within the range of suspension velocities measured in the actual experiment. It can therefore be concluded that the peanuts and clay-heavy colds suspension velocity measurement process calculated based on the coupled CFD-DEM model simulations is reliable and the gas-solid two-phase flow model has been validated.

5. Discussion

The suspension speed is the central basis of the combination of a fan and vibrating sieve and is essential for the effective operation of the cleaning and sorting device on peanut combine harvesters. Zhichao Hu et al. [25] presented the working principle and structural design of a 4LH2-type, half-feed, and self-propelled peanut combine harvester scavenging system. Lianxing Gao et al. [26] proposed the principle of combined double suction air outlet and vibrating sieve type of cleaning and studied the effect of vibrating frequency, vibrating sieve, height of suction air outlet, and speed of the fan on the rate of cleaning loss and trash content of dialogue peanut through cleaning performance tests. Yifei Wang [13] measured the suspension velocity of peanut threshing material, including pods, long stalks, short stalks, stalks, and miscellaneous leaves, and used coupled CFD-DEM to simulate the material transport process and flow field simulation for the cleaning and sorting device of the peanut combine harvester. Previous research has primarily focused on mechanism design, structure optimization, and material testing for sandy loam soil peanut cleaning machinery, with little attention given to peanut cleaning machinery designed for hilly, mountainous clay-heavy land. Further research is needed in this area to develop effective and efficient cleaning machinery for this type of soil. Therefore, the research in this paper has certain significance in filling the gap in the industry.
In measurements and simulation experiments with different kernel peanut pods and different-shaped clods of clay-heavy soil, we found that there were significant differences in the suspension state of the differently shaped materials in the suspension wind tunnel. The floating range of double- kernel peanut pods and strips of clods within the wind field is significantly greater than that of single-kernel peanut pods and lump clods, which is directly related to the windward area of the different shapes of material within the wind field. However, there was no significant difference in the suspension velocity values for the different shapes of material when a steady-state suspension was reached. As the proportion of double-kernel peanut pods is much greater than that of single-kernel peanut pods at harvest, this suspension characteristic may be more favorable for cleaning and sorting peanuts grown on clay-heavy soils at harvest.
In the measurement experiments with clods of different water content, we found that the higher the water content, the more significant the adhesion effect of the clods in the conical observation tube. During the measurement experiments with a water content of 23.9%, the difficulty of the experiments increased considerably and had to be interrupted several times because the clods adhered to the wall of the tube. During harvesting, it is common for the soil water content to be close to or even exceed the experimental value, resulting in clods sticking to the machine. To avoid operational difficulties and improve sorting efficiency, it is recommended to harvest when the soil water content is below 20%.
This study is based on existing peanut half-feed combine harvesting machinery to propose cleaning and sorting ideas applicable to peanuts grown on clay-heavy soils and is limited to full-feed combine harvesting machinery. In addition, the final shape of the cleaning and sorting device will include the design and optimization of the vibrating sieve and the grouping of the vibrating sieve with the negative pressure soil separation mechanism, which will be the main issue addressed in our subsequent research.

6. Conclusions

  • The suspension velocity of the large number of clay-heavy clods present in the fruit collection box during the combined harvest and Wanhua 17 peanuts was measured by means of a suspension velocity measurement device. The suspension velocities of single and double kernel peanuts were 8.34~9.40 m/s and 8.13~9.51 m/s, respectively. Under 20.4% water content and lumpy conditions, the suspension velocities of smaller clods, side by side clods, and larger clods were 12.61~14.30 m/s, 14.16~15.76 m/s, and 16.44~18.72 m/s, respectively; under 20.4% water content and smaller clods, the suspension velocities of lumpy and strips of clods were 12.61~14.30 m/s, 11.90~14.13 m/s, respectively; under lumpy and smaller clods, the suspension velocity at 17.6%, 20.4%, and 23.9% water content ranged from 12.38 to 14.20 m/s, 12.61 to 14.30 m/s, and 12.62 to 14.49 m/s, respectively. Suspension velocity values and suspension characteristics were obtained for different kernelled peanut pods and clay-heavy clods under different conditions.
  • Numerical simulations of the experimental process of peanuts and clay-heavy clods in a suspension velocity measurement device based on CFD-DEM theory were carried out, and the model was validated. The coupling calculations gave suspension velocities of 8.52~9.43 m/s for single-kernel peanuts and 7.74~9.20 m/s for double-kernel peanuts; under 20.4% water content and lumpy conditions, suspension velocity was 12.18~14.62 m/s for smaller clods, 14.30~16.23 m/s for side-by-side clods, and 16.99 ~19.80 m/s for larger clods; under 20.4% water content and smaller clod conditions, suspension velocity was 12.18~14.62 m/s for lumpy clods and 12.67~14.76 m/s for strips; under the conditions of lumpy and smaller clods, the suspension velocity was 11.36~13.87 m/s at 17.6% water content, 12.18~14.62 m/s at 20.4%, and 12.52~14.09 m/s at 23.9%; The simulation results for the suspension velocity values of each material are within the ideal range of error from the experimental results.
  • The feasibility of the idea of fruit-soil separation through suspension velocity was verified by comparing and analyzing the difference in suspension velocity between peanuts and clay-heavy clods of various sizes, providing a basis for subsequent research.

Author Contributions

Conceptualization, M.Q. and Y.J.; methodology, M.Q., Y.J., M.C. and W.L.; software, M.Q. and Y.J.; validation, M.Q. and W.L.; formal analysis, M.Q., M.C. and L.S.; investigation, M.Q. and Y.J.; resources, F.W. and F.G.; data curation, M.Q. and Y.J.; writing—original draft preparation, M.Q.; writing—review and editing, M.Q. and W.L.; visualization, M.Q. and L.S.; supervision, F.W., F.G. and Z.H.; project administration, F.W. and F.G.; funding acquisition, F.W. and F.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the National Peanut Industry Technology System (grant No. CARS-13).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, Z.; Gao, L.; Chen, C.; Butts, C. Analysis on technology status and development of peanut harvest mechanization of China and the United States. Trans. CSAM 2017, 48, 1–21. [Google Scholar]
  2. Wang, Q.; Yu, Z.; Zhang, Y.; Li, J.; Peng, B.; Wang, B.; Hu, Z. Experimental Study of a 4HLB-4 Half-Feed Four-Row Peanut Combine Harvester. Agronomy 2022, 12, 3094. [Google Scholar] [CrossRef]
  3. Klinzing, G.E. A review of pneumatic conveying status, advances and projections. Powder Technol. 2018, 333, 78–90. [Google Scholar] [CrossRef]
  4. Li, J.; Yuan, W.; Zhai, H.; Zhu, S. Application of pneumatic conveying in agricultural machinery. J. Chin. Agric. Mech. 2023, 44, 85–92. [Google Scholar]
  5. Tang, H.; Xu, F.; Xu, C.; Zhao, J.; Wang, Y.-J. The influence of a seed drop tube of the inside-filling air-blowing precision seed-metering device on seeding quality. Comput. Electron. Agric. 2023, 204, 107555. [Google Scholar] [CrossRef]
  6. Tang, H.; Xu, C.; Wang, Z.; Wang, Q.; Wang, J. Optimized Design, Monitoring System Development and Experiment for a Long-Belt Finger-Clip Precision Corn Seed Metering Device. Front. Plant Sci. 2022, 13, 67. [Google Scholar] [CrossRef]
  7. Ji, W.; Lu, Y.; Yang, M.; Wang, J.; Zhang, X.; Zhao, C.; Xia, B.; Wu, Y.; Ying, R. Geochemical Characteristics of Typical Karst Soil Profiles in Anhui Province, Southeastern China. Agronomy 2023, 13, 1067. [Google Scholar] [CrossRef]
  8. Gorial, B.; O’callaghan, J. Separation of grain from straw in a vertical air stream. J. Agric. Eng. Res. 1991, 48, 111–122. [Google Scholar] [CrossRef]
  9. El-Sayed, A.S.; Yahaya, R.; Wacker, P.; Kutzbach, H.D. Characteristic attributes of the peanut (Arachis hypogaea L.) for its separation. Int. Agrophysics 2001, 15, 225–230. [Google Scholar]
  10. Dai, F.; Guo, W.; Song, X.; Shi, R.; Qu, J.; Zhao, W. Measurement and simulation of the suspension velocity of flax threshing material using CFD-DEM. Int. J. Agric. Biol. Eng. 2021, 14, 230–237. [Google Scholar] [CrossRef]
  11. Liang, Y.; Tang, Z.; Zhang, H.; Li, Y.; Ding, Z.; Su, Z. Cross-flow fan on multi-dimensional airflow field of air screen cleaning system for rice grain. Int. J. Agric. Biol. Eng. 2022, 15, 223–235. [Google Scholar] [CrossRef]
  12. Gao, L.; Zhang, W.; Du, X.; Liu, X.; Yang, J.; Liu, M. Experiment on aerodynamic characteristics of threshed mixtures of peanut shelling machine. Trans. Chin. Soc. Agric. Eng. 2012, 28, 289–292. [Google Scholar]
  13. Wang, Y. Design of the Cleaning Device for the Peanut Combine Harvester and Performance Studies. Master’s Thesis, Hebei Normal University of Science & Technology, Qinhuangdao, China, 2022. [Google Scholar]
  14. Hou, H.; Cui, Q.; Guo, Y.; Zhang, Y.; Sun, D.; Lai, S.; Liu, J. Design and test of air-blowing type grain and oil crop debris cleaning suspension speed measurement device. Trans. Chin. Soc. Agric. Eng. 2018, 34, 43–49. [Google Scholar]
  15. Jia, Z.; Weidong, Z.; Bin, C. The development of the testbed for stuff of agriculture in suspension speed. J. Heilongjiang August First Land Reclam. Univ. 1998, 10, 56–59. [Google Scholar]
  16. Liu, L.; Hao, S.; Zhang, M.; Liu, D.; Jia, F.; Quan, L. Numerical simulation and experiment on paddy ventilation resistance based on CFD-DEM. Trans. CSAM 2015, 46, 27–32. [Google Scholar]
  17. Dong, H.; Zhang, B.; Jiang, T.; Zhang, Y.; Qu, J.; Chen, C.; Xiao, Y.; Ding, Y.; Xi, X. Design and Optimization of Rice Grain Screening System Based on DEM–CFD Coupled Rice Seed Testing Platform. Agronomy 2022, 12, 3069. [Google Scholar]
  18. Li, J.; Tong, J.; Hu, B.; Wang, H.; Mao, C.; Ma, Y. Calibration of parameters of interaction between clayey black soil with different moisture content and soil-engaging component in northeast China. Trans. CSAE 2019, 35, 130–140. [Google Scholar]
  19. Xu, X.; Wei, H.; Yan, J.; Bao, G.; Du, Y.; Xie, H. Parameter calibration of peanut pods discrete element simulation. J. Chin. Agric. Mech. 2022, 43, 81–89. [Google Scholar]
  20. Xu, J.; Hou, J.; Wu, W.; Han, C.; Wang, X.; Tang, T.; Sun, S. Key Structure Design and Experiment of Air-Suction Vegetable Seed-Metering Device. Agronomy 2022, 12, 675. [Google Scholar] [CrossRef]
  21. Liu, X.; Shi, L.; Xu, X. Status of research on Euler-Lagrange method for dense gas-solid two-phase flow. Chin. J. Comput. Mech. 2007, 24, 166–172. [Google Scholar]
  22. Xiangyu, W.; Honglei, J.; Shengwei, Z.; Hongfang, Y.; Gang, W.; Tianyou, C. Test of Suspension Velocity of Granular Fertilizer Based on EDEM-Fluent Coupling. Nongye Jixie Xuebao/Trans. Chin. Soc. Agric. Mach. 2020, 51, 69–77. [Google Scholar]
  23. Fan, W. Design and Experimental Research on Key Components of Small Self-Walking White Radish Harvester. Master’s Thesis, Guizhou University, Guiyang, China, 2022. [Google Scholar]
  24. Tang, Z.; Gong, H.; Wu, S.; Zeng, Z.; Wang, Z.; Zhou, Y.; Fu, D.; Liu, C.; Cai, Y.; Qi, L. Modelling of paddy soil using the CFD-DEM coupling method. Soil Tillage Res. 2023, 226, 105591. [Google Scholar] [CrossRef]
  25. Hu, Z.; Peng, B.; Yin, W.; Wang, H.; Tian, L.; Hu, L. Design of 4LH2 type half-feed and self-propelled peanut combine. Trans. CSAE 2008, 24, 148–153. [Google Scholar]
  26. Tang, Z.; Gong, H.; Wu, S.; Zeng, Z.; Wang, Z.; Zhou, Y.; Fu, D.; Liu, C.; Cai, Y.; Qi, L. Design and test on cleaning device of peanut pods with double air-suction inlets with vibration screen. Nongye Jixie Xuebao/Trans. Chin. Soc. Agric. Mach. 2015, 46, 11. [Google Scholar]
Figure 1. Suspension velocity measurement device for peanuts and clay-heavy clods: (a) Structure of the device: 1. centrifugal fan; 2. buzzer; 3. stable section; 4. holding net; 5. test section; 6. flow stabilizer; 7. inlet; 8. damping net; 9. square and round section; 10. frame; (b) Experiment benches.
Figure 1. Suspension velocity measurement device for peanuts and clay-heavy clods: (a) Structure of the device: 1. centrifugal fan; 2. buzzer; 3. stable section; 4. holding net; 5. test section; 6. flow stabilizer; 7. inlet; 8. damping net; 9. square and round section; 10. frame; (b) Experiment benches.
Agronomy 13 01735 g001
Figure 2. Classification of peanuts and clay-heavy clods: (a) Field harvesting in action; (b) Classification of peanuts; (c) Smaller clods; (d) Side-by-side clods; (e) Larger clods.
Figure 2. Classification of peanuts and clay-heavy clods: (a) Field harvesting in action; (b) Classification of peanuts; (c) Smaller clods; (d) Side-by-side clods; (e) Larger clods.
Agronomy 13 01735 g002
Figure 3. Measurement process of suspension velocity of peanuts and clay-heavy clods: (a) Peanuts in suspension; (b) Clay-heavy clods in suspension.
Figure 3. Measurement process of suspension velocity of peanuts and clay-heavy clods: (a) Peanuts in suspension; (b) Clay-heavy clods in suspension.
Agronomy 13 01735 g003
Figure 4. Individual material particles and their discrete element models: (a) Peanut pod; (b) Smaller lumpy clod; (c) Smaller strips; (d) Side-by-side clod; (e) Larger clod.
Figure 4. Individual material particles and their discrete element models: (a) Peanut pod; (b) Smaller lumpy clod; (c) Smaller strips; (d) Side-by-side clod; (e) Larger clod.
Agronomy 13 01735 g004
Figure 5. 3D model of the suspended wind tunnel and the generated body mesh.
Figure 5. 3D model of the suspended wind tunnel and the generated body mesh.
Agronomy 13 01735 g005
Figure 6. Suspension positions for different material suspension velocity simulation measurements: (a) Peanut pod; (b) Smaller lumpy clod; (c) Smaller strips; (d) Side-by-side clod; (e) Larger clod.
Figure 6. Suspension positions for different material suspension velocity simulation measurements: (a) Peanut pod; (b) Smaller lumpy clod; (c) Smaller strips; (d) Side-by-side clod; (e) Larger clod.
Agronomy 13 01735 g006
Figure 7. Comparison of the suspension velocity of each material: (a) Single and double kernel peanuts; (b) Different masses of clods; (c) Different shapes of clods; (d) Clods at different water contents. (a, b, c are significance levels for each factor, in increasing order; *, **, *** are inter-factor significance levels, in increasing order).
Figure 7. Comparison of the suspension velocity of each material: (a) Single and double kernel peanuts; (b) Different masses of clods; (c) Different shapes of clods; (d) Clods at different water contents. (a, b, c are significance levels for each factor, in increasing order; *, **, *** are inter-factor significance levels, in increasing order).
Agronomy 13 01735 g007
Figure 8. Suspension position, velocity and flow field distribution of each material: (a) Peanut pod; (b) Smaller lumpy clod; (c) Smaller strips; (d) Side-by-side clod; (e) Larger clod.
Figure 8. Suspension position, velocity and flow field distribution of each material: (a) Peanut pod; (b) Smaller lumpy clod; (c) Smaller strips; (d) Side-by-side clod; (e) Larger clod.
Agronomy 13 01735 g008
Figure 9. Comparison of simulation and actual experimental results.
Figure 9. Comparison of simulation and actual experimental results.
Agronomy 13 01735 g009
Figure 10. Comparison chart of model validation results: (a) Single and double kernel peanuts; (b) Different masses of clods; (c) Different shapes of clods; (d) Clods at different water contents.
Figure 10. Comparison chart of model validation results: (a) Single and double kernel peanuts; (b) Different masses of clods; (c) Different shapes of clods; (d) Clods at different water contents.
Agronomy 13 01735 g010
Table 1. Technical parameters of the suspension velocity measuring device for peanuts and clay-heavy clods.
Table 1. Technical parameters of the suspension velocity measuring device for peanuts and clay-heavy clods.
Technical ParametersValues
External dimensions (L × W × H)/(mm × mm × mm)810 × 755 × 3359
Power/kW2.2
Flow rate/(m3·h−1)2450
Inner diameter of the flow stabilizer/mm90
Inner diameter of conical tube tip/mm190
Inner diameter of stable section/mm500
Conical tube inclination/(°)2.6
Holding and damping nets20-mesh nylon mesh
Table 2. Basic material parameters.
Table 2. Basic material parameters.
MaterialsParametersValues
Peanut podsPoisson’s ratio0.4
Density/kg·m−3465.08
Shear modulus/Pa6.5 × 106
Clay heavy soilPoisson’s ratio0.46
Density/kg·m−32054
Shear modulus/Pa1 × 106
Rigid PVCPoisson’s ratio0.32~0.47
Density/kg·m−31418
Shear modulus/Pa1.50 × 109~3.92 × 109
Peanut pods–rigid PVCCoefficient of restitution0.275
Coefficient of static friction0.583
Coefficient of rolling friction0.104
Clay heavy soil–rigid PVCCoefficient of restitution0.31
Coefficient of static friction0.522
Coefficient of rolling friction0.43
Table 3. Suspension velocity of peanuts of different kernels.
Table 3. Suspension velocity of peanuts of different kernels.
ClassificationTest NumberExternal DimensionsMass/gAverage Suspension Height Range/mmSuspension Velocity/m·s−1
Length/mmWidth/mmHeight/mm
Single-Kernel Peanuts123.5716.1815.341.8386~4568.37~9.24
224.5316.0615.541.74359~4578.35~9.61
324.2916.3715.542.09344~4138.89~9.83
424.9116.315.51.48416~5867.04~8.85
526.918.6218.612.17402~4698.22~9.03
625.3716.8816.912376~4248.75~9.38
726.7417.5215.362.12359~4518.42~9.61
824.1916.215.881.84356~4348.63~9.66
923.6717.4715.531.68372~4438.52~9.43
1026.2318.2216.912.28376~4728.18~9.38
x ¯
25.0416.9516.111.92374.6~460.58.34~9.40
Double-Kernel Peanuts139.9415.8516.033.24318~4578.35~10.23
237.1816.5117.373.06354~5047.84~9.69
344.0117.0818.053.44392~4718.20~9.16
436.515.4216.052.63364~4678.24~9.54
544.717.4118.173.26418~5367.51~8.82
636.3916.39172.86357~4568.37~9.64
741.0316.4517.843.74349~4488.46~9.76
833.815.2815.852.35352~4618.31~9.71
936.8115.1716.272.48389~5137.74~9.20
1041.9916.1717.153.52375~4648.27~9.40
x ¯
39.2416.1717.023.06366.8~477.78.13~9.51
Table 4. Suspension velocity of peanuts of different masses of clods.
Table 4. Suspension velocity of peanuts of different masses of clods.
Test NumberSmaller ClodsSide-by-Side ClodsLarger Clods
Average Suspension Height Range/mmSuspension Velocity/m·s−1Average Suspension Height Range/mmSuspension Velocity/m·s−1Average Suspension Height Range/mmSuspension Velocity/m·s−1
1278~39012.54~14.86288~37114.15~16.05319~39215.79~17.60
2277~36612.99~14.88296~36214.34~15.85261~32917.33~19.27
3313~41212.15~14.07277~33614.91~16.33192~28817.46~21.59
4254~31813.96~15.43305~37514.07~15.63244~34216.99~19.80
5282~34613.38~14.76296~37414.09~15.85271~34117.02~18.96
6276~33513.60~14.90324~38613.84~15.18351~42215.12~16.77
7356~45011.52~13.18323~38213.93~15.21284~36916.33~18.58
8384~43611.74~12.65322~39613.65~15.23326~38315.99~17.41
9314~39612.43~14.05281~36414.30~16.23319~39415.74~17.60
10288~41012.18~14.62287~35814.43~16.08294~38216.02~18.29
x ¯ 302~38612.61~14.30300~37014.16~15.76286~36416.44~18.52
Table 5. Suspension velocity of peanuts of different shapes of clods.
Table 5. Suspension velocity of peanuts of different shapes of clods.
Test NumberLumpy ClodsStrip of Clods
Average Suspension Height Range/mmSuspension Velocity/m·s−1Average Suspension Height Range/mmSuspension Velocity/m·s−1
1221~32913.73~16.29265~35013.30~15.16
2223~33913.52~16.23282~38512.63~14.76
3274~34113.48~14.95275~39012.54~14.93
4268~35213.26~15.09285~35513.20~14.69
5282~37412.84~14.76280~39512.45~14.81
6286~38112.71~14.67294~41012.18~14.49
7308~38912.56~14.18274~38612.61~14.95
8281~37812.76~14.78292~38512.63~12.53
9289~39612.43~14.60276~38912.67~14.76
10299~40312.31~14.37296~37412.84~14.44
x ¯ 273~36812.95~14.97282~38312.71~14.75
Table 6. Suspension velocity of clods at different water contents.
Table 6. Suspension velocity of clods at different water contents.
Test Number17.6%20.4%23.9%
Average Suspension Height Range/mmSuspension Velocity/m·s−1Average Suspension Height Range/mmSuspension Velocity/m·s−1Average Suspension Height Range/mmSuspension Velocity/m·s−1
1278~33613.58~14.86277~36612.99~14.88302~39612.43~14.31
2326~41112.17~13.79313~41212.15~14.07288~36912.93~14.62
3307~38912.56~14.20254~31813.96~15.43297~39412.47~14.42
4317~37512.82~13.98282~34613.38~14.76287~37612.80~14.65
5307~40612.25~14.20276~33513.60~14.90290~36413.03~14.58
6298~37212.87~14.40356~45011.52~13.18312~39112.52~14.09
7310~39412.47~14.13384~43611.74~12.65294~40612.25~14.49
8322~46011.36~13.87314~39612.43~14.05291~37412.84~14.55
9297~41512.10~14.92288~41012.18~14.62286~36712.97~14.67
10307~39912.38~14.20302~38612.61~14.30294~38612.61~14.49
x ¯ 306~42911.86~14.22278~39012.54~14.86293~42211.98~14.51
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Qin, M.; Jin, Y.; Luo, W.; Wu, F.; Shi, L.; Gu, F.; Cao, M.; Hu, Z. Measurement and CFD-DEM Simulation of Suspension Velocity of Peanut and Clay-Heavy Soil at Harvest Time. Agronomy 2023, 13, 1735. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy13071735

AMA Style

Qin M, Jin Y, Luo W, Wu F, Shi L, Gu F, Cao M, Hu Z. Measurement and CFD-DEM Simulation of Suspension Velocity of Peanut and Clay-Heavy Soil at Harvest Time. Agronomy. 2023; 13(7):1735. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy13071735

Chicago/Turabian Style

Qin, Mingyang, Yu Jin, Weiwen Luo, Feng Wu, Lili Shi, Fengwei Gu, Mingzhu Cao, and Zhichao Hu. 2023. "Measurement and CFD-DEM Simulation of Suspension Velocity of Peanut and Clay-Heavy Soil at Harvest Time" Agronomy 13, no. 7: 1735. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy13071735

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop