1. Introduction
In order to maintain sustainability, intelligent agriculture and precision agriculture have become an important research direction of modern agriculture [
1,
2]. Variable adjustment technology of flow rate in irrigation pipelines helps the irrigation and fertilization operation match the actual needs of crops and, therefore, achieves the goal of water saving, fertilizer saving, and production increment [
3,
4,
5].
Research on flow regulation in the field of agricultural irrigation focuses on three main areas. (1) PWM (Pulse Width Modulation) variable technique: Mo Jinqiu [
6] and Ramon Salcedo [
7] analyze pressure pulses and propose mitigation measures for large PWM variable sprinkler physical systems. Li Jianian [
8] and Chang Chen [
9] apply PWM variable technique combined with Venturi fertilizer applicator to fertilizer application control of the fertilizer applicator. (2) H Zhu [
10], Wang Lixai [
11] regulate the flow rate by adjusting the line pressure with a frequency converter and apply it to water-fertilizer-pharmacy integration systems and variable spraying systems. (3) Applying flow control by flow control valves, Sun Caizhen [
12] designs a valve-regulated proportional valve fertilizer applicator to achieve control of the fertilizer supply flow.
Using flow control valves to regulate flow may make the regulation more stable and continuous. Flow control and valve control techniques are categorized as follows. (1) An intelligent control strategy based on PID control combined with multiple algorithms. For example, Song Lepeng [
13,
14] uses the variable theory domain adaptive fuzzy PID control algorithm to regulate electromechanical flow control valves, and Young-Rae Ko [
15,
16] directly uses PID closed-loop controllers as the control strategy. (2) Combined with the mathematical model control strategy of valve regulation law, Li Jinyang [
17] uses BP neural networks to establish the relationship between valve opening and flow, Do-Hun Chin [
18] and Zaryankin, A. E. [
19] use simulation to obtain the operating characteristics of the regulating valve, and Rezazadeh [
20] utilizes the flow equation for the flow and differential pressure relationship of MCOP, and uses for flow control. Additionally, the flow rate can also be estimated by modeling the regulation law of the valve body. Zhang [
21] uses parameters including oil temperature, differential pressure, and spool displacement as input and flow rate as output with the help of neural networks to achieve the prediction of flow rate as a feedback signal to the flow controller.
IoT-based remote control enhances system automation and production efficiency [
22]. Remote flow control technology is becoming increasingly used in modern agriculture [
23]. Tu lv [
24] and H. Ruan [
25] have implemented remote control of the flow with the help of PWM technology and variable frequency pump technology. However, its implementation is relatively complex, which has led to most of the research on remove control methods for regulating valves staying in the laboratory and making it difficult to promote them on a large scale. Compared with other developed agricultural economies, the penetration rate of such technology in China is much lower [
26]. Therefore, achieving acceptable flow control accuracy and responsiveness with simpler structures and control methods has become a crucial issue.
This paper (1) designed an electric flow regulator, (2) established mathematical relationships between valve parameters through performance tests, and (3) developed a control program and flow measurement program to operate at a remote terminal. The aim is to achieve remote and accurate control and measurement of the flow rate.
2. Actuating Component Design
The integrated flow rate regulation device consists of a local actuating component and a remote control system. The actuating component consists of two parts: the valve body structure and the driver module.
2.1. Valve Body Structure
The valve body is angularly structured. Its flow path is simple and resistant to precipitation blockage. The valve body adapts well to the agricultural irrigation system to transport fertilizer, liquid, and other viscous media and contains suspended matter and other conditions. Its vertical installation also has a certain self-cleaning function. The spool is made of elastic material so that the valve has the ability to close when the opening degree is zero, eliminating leakage losses.
The structure of the body is shown in
Figure 1. The main components are plunger, rubber cushion, inlet shrinkage tube, and valve body shell; four parts. The plunger is used to drive the radial action of the adjustment and the outlet section of the inlet shrinkage tube to control the shape of the flow channel together. The rubber gasket is used to seal the gap between the plunger and the inlet shrinkage outlet when the valve body is closed and fixed by the nut on the water ward side of the plunger. The inlet shrinkage tube consists of two parts, the inlet shrinkage section and the outlet section, which are used to gather the fluid entering the valve body with an inner diameter of 32 mm at the outlet. The valve body shell is a three-way structure, and consists of three parts: water inlet, chamber, and water outlet. The valve body fixes the internal structure of the valve body. The inlet shrinkage section and the valve body shell are made of PVC (polyvinyl chloride).
2.2. Driver Module
The driver module of the device is powered by a stepper motor. The power output shaft of the stepper motor is connected to the screw through a coupling to transmit the torque, and the radial movement of the plunger is ensured by the limitation of the plunger’s degrees of freedom by the plunger guide to adjust the valve opening, as shown in
Figure 1.
To ensure that the motor drives the torque in the working pressure range of the valve body, it is necessary to calculate the motor output shaft torque and the screw torque, and then compare. Screw the driven plunger with a screw drive, which is mainly subject to resistance to the action of the plunger face level before and after the pressure difference. The required screw torque Tn can be obtained from the following formula [
27].
where
D is the diameter of the plunger facing the horizontal surface in mm, Δ
p is the maximum pressure difference of the valve standard working condition in MPa,
n is the number of threads,
S is the thread lead in mm,
d is the thread middle distance in mm,
f is the friction coefficient of the bolt (it is also a dimensionless number and can be checked in the spiral sub by material),
α is the thread tooth angle in °.
Because of the existence of mechanical losses in the device, the motor torque selected should be greater than the requested torque Tn. It follows that the final choice of 42HS6315B4 type stepper motor with DV542C two-phase digital stepper motor driver. Both devices can be driven by 12 to 24 V DC.
Commands are sent to the stepper motor driver by the controller, including the rotational direction, number of pulses, and pulse frequency. The stepper motor’s angular displacement is controlled by the number of pulses, and the frequency of the pulses determines how quickly the motor rotates. This enables precise positioning of the opening of the regulating valve and control of the regulating speed.
3. Experimental Study
The development of the device control system is based on the hydraulic performance relationship of the valve body, and the opening-pulse relationship of the drive module, so the performance relationship between these two modules should be investigated through experimental studies.
3.1. Testing Device
The structure of the test device is shown in
Figure 2a. The test device consists of a constant pressure water supply module, a flow regulation module, and a field control device. The constant pressure water supply module provides constant water pressure using variable frequency pumps, the flow regulation module includes the regulating valve body and its driving module, and the control device consists of electromagnetic flow meter, pressure transmitter, gate valve, field controller and field control terminal. A physical diagram of the flow regulation module and pressure transmitter installation is shown in
Figure 2b.
3.2. Experimental Design
Characteristic flow coefficient test: The differential pressure Δp before and after the valve body is maintained at an average of 15 levels in 0.02–0.30 Mpa, and the pipeline flow q, pre-valve pressure p1 with the valve body fully open is measured.
Performance parameter test: Keep the differential pressure Δp before and after the valve body at 5 levels in 0.05–0.25 MPa. The valve opening is selected uniformly at 25 levels in 0–50% and 12 levels uniformly in 50–100%. The flow rate, pre-valve pressure, and back-end pressure are q, p1 and p2, respectively.
Opening and pulse relationship test: The valve opening is selected uniformly at 25 levels in 0–50% and 12 levels uniformly in 50–100%. Then, measure the actual opening of the valve body la and record the cumulative number of pulses np.
During the test, the test level is adjusted continuously, and the pressure and flow rate under the current test level is stabilized for 15 s before switching to the next test level, the data is collected every 5 s, the measured data is obtained by averaging the 3 data before switching the level after the test, and each group is repeated for 3 times.
3.3. Test Index
Pressure head: Pressure transmitter (model: PCM300, accuracy 0.5%) is used before and after the flow adjustment module to monitor.
Flow rate: The pipeline flow rate is monitored by the electromagnetic flow meter (EMF5000, accuracy 0.5%).
Cumulative number of pulses: The sum of the number of pulses of the stepper motor action is monitored and recorded by the controller.
Valve body relative opening: Expressed by the ratio of the accumulated number of pulses of the stepper motor action and the total number of pulses of its full stroke, which is monitored and recorded by the controller.
The actual opening of the valve body: Measured by the vernier calipers, the spool is in the fully closed position and recorded as 0 mm.
3.4. Analysis of Test Results
3.4.1. Valve Body Hydraulic Characteristics
The flow path of the valve body will make the medium subject to local pressure loss. By adjusting the opening of the valve body, one can change this loss, and thus regulate the flow. When the regulating valve is in operation, the relationship between the lead pressure
p1, the relative opening
lre, and the flow rate
q through the valve body can be obtained from the following formula [
28].
where
q is the flow rate through the valve body, m
3/h;
KV is the dimensionless characteristic flow coefficient obtained through tests;
CPP is the post-valve pressure prediction coefficient, a dimensionless number;
kL is the relative flow coefficient, which is only related to the relative opening lre;
p1 is the pre-valve pressure, MPa;
ρ is the density of the medium in kg∙L
−1 (the test uses clear water, the density is about 1 kg∙L
−1)
The valve body characteristic flow coefficient
KV is determined by using the relationship between flow and pressure at the full start of the valve body, and is calculated by using the following formula [
28].
where Δ
p is the differential pressure Δ
p before and after the valve body, MPa;
n is the number of data groups, (the number of levels in this study is 15, with 3 replication groups, so
n is 45, and its calculation table (
Table 1) is as follows). Then,
KV = 84.61.
The post-valve pressure prediction coefficient CPP is the proportionality coefficient between the pre-valve pressure p1 and the post-valve pressure p2,100 under the same p1 when the valve body is fully open, which is a constant in the same pipeline. The CPP under this device is 0.375, and the measured data and the calculated data will use Pearson correlation analysis to get R2 of 0.999, indicating an ideal fitting performance.
The relative flow coefficient
kL is measured experimentally by performance parameters. The relative flow coefficients at different pressures with relative openings are shown in
Figure 3. It is observed that the relative flow coefficient
kL varies less at different pressures, and increases with the relative opening
lre at the same pressure. In the case of a small relative opening (0% to 5%), the relative flow coefficient at different pressures does not change immediately with the increase in opening, and the greater the pressure, the faster the relative flow coefficient changes. The main reason is the anti-leakage design of the valve body. The existence of the pre-closing volume when the spool is closed and greater water pressure can make the spool and the valve seat separate more quickly.
A polynomial was fitted to the test data to establish the relationship between opening and flow. The relationship established is as follows.
The correlation coefficient R2 between the fitted and measured values at all five pressure levels is above 0.996, which is a good fit.
3.4.2. The Relationship between the Valve Opening and the Number of Pulses of the Stepper Motor
In order to verify the feasibility of using the ratio of the accumulated number of pulses from the stepper motor action to the total number of pulses from its full stroke as the relative opening of the valve body, a relationship between the actual opening of the valve and the number of pulses from the stepper motor is established.
The opening degree
l of the control valve is proportional to the number of driving pulses
np of the stepper motor, and its proportionality factor
k can be calculated by the number of threads
n, thread guide
S, and the number of steps divided by the driver
nc. Its relationship is shown in Equation (5) [
27].
If the driver divides one revolution into 25,600 steps, substituting the parameters yields a proportionality factor k of approximately 1.142 × 10−4.
The theoretically calculated value is close to the measured value,
R2 is 0.999, and the proposed relationship is reliable. It is shown that the method of characterizing the relative opening by the ratio of the number of pulses is reliable (
Figure 4).
4. Control System Design
4.1. Hardware Design
The control system is composed of a field acquisition and execution subsystem, IOT cloud box, and remote monitoring subsystem. The overall framework is shown in
Figure 5. The field acquisition and execution subsystem takes a programmable logic controller as the core, and connects pressure transmitter, motor driver, and other devices respectively. The remote monitoring subsystem takes cloud platform as the core, which can realize the monitoring of the acquisition data through remote management software and process the acquired data and output the valve opening size synchronously by running the opening decision algorithm, as well as the control command issuance. As the core hub, the IOT connects the site monitoring subsystem and remote decision monitoring subsystem through RS-232 serial port and 4G network, respectively, and realize the communication between data information and control signal at the local end and the cloud platform.
4.2. Software Design
The device control software is developed based on the relationship among the pressure in front of the valve, the valve opening, the pipeline flow, the accumulated pulses, and the actual opening. Its operation flow is shown in
Figure 6. The steps are as follows.
(1) Configure initialization parameters such as target flow rate and allowable error. (2) Read the pressure before the valve and the number of opening pulses to calculate the current flow rate. (3) Judge whether the current flow rate is in the range of a reasonable flow rate interval; if it is within the range, the program ends; if it is not within the range, calculate the target opening pulse number. (4) Calculate a difference between the target opening pulse number and the current pulse number, and adjust the opening to the target opening based on its positive or negative value to drive the motor rotation direction. (5) Update the adjusted number of pulses to the current number of pulses. (6) The program jumps back to step 2 until the end of the program.
The measurement of the flow rate can be calculated directly into Equation (2). The calculation of its opening is the reverse of the former: The target flow rate as q with other parameters into Equation (2) to get the corresponding flow coefficient kL, then its inverse function lre = f(kL) into Equation (3) to obtain the relative opening lre, multiplied by the number of pulses corresponding to the full stroke of the stepper motor to obtain the number of pulses np1 corresponding to the target flow rate q0.
5. Performance Test
The validation test was conducted with a remote control system attached to the performance test platform (
Figure 2). The sensors used are consistent with the performance test platform. The control software was run on the remote control terminal to make decisions and give regulation commands. The test data were collected through the local control module, and the data is collected in the same way as above. Through several trial adjustments, the pulse frequency of the stepper motor driver is set to 400 HZ. This opening adjustment speed can make the flow adjustment take into account the adjustment speed and stability, and at the same time, the water hammer phenomenon will not occur. Our goal is to test the regulation performance of the flow rate, and the flow measurement performance, respectively.
The test is set up with multiple levels of test pressure p1, target flow q0, and target opening lre. The flow measurement software output data qm and the actual data qa measured by the electromagnetic flowmeter are recorded once per second. The regulation time ts is used to evaluate the regulation rapidity, the overshoot Δh to evaluate the regulation stability, and the steady-state error ess to evaluate the regulation accuracy. The relative error δ of qa and qm evaluate measurement accuracy. All tests are repeated 5 times.
The average values of
qs,
ts, Δ
h,
ess, the mean values of
qa and
qm and
δ of the measurement test for each working condition of the regulation experiment are calculated in
Table 2.
From the table, we can obtain the following information. In the regulation performance tests, the final steady-state flow rates achieved were all similar to the target flow rates, with steady-state errors within 6%, and in most cases within 2%. The higher the pressure and flow rate, the higher the regulation accuracy. This meets the requirements of flow regulation in production. The rapidity and stability of the flow rate regulation is also good. The regulation time decreases as the pressure increases, and increases as the target flow rate increases, with a maximum of 78 s. Overshoot is within 5%, with less overshoot at moderate pressures.
In the test performance tests, the measured flow rates were compared at multiple opening levels at multiple pressures, and it can be seen that the measured flow rates match the actual flow rates with an error of 3% or less. The accuracy of the measurement is higher at higher opening and flow levels, but the effect of pressure on the accuracy of the measurement is not significant.
6. Conclusions
- (1)
A remote control-based integrated device for flow regulation and measurement is designed. The device has a stepper motor as the power input, an angular regulating valve as the actuating component, a pressure transmitter as the collecting component, a PLC and a cloud box as the signal transmission relay, and a remote control terminal as the data processing and human-computer interaction platform, enabling remote pipeline flow regulation and measurement.
- (2)
The device’s regulation properties were studied. The flow coefficient KV of the valve body is 84.61, and the regulation curve is in accordance with the law of fast-opening characteristics. There is no leakage loss when closed within 0.3 MPa of the pressure difference between the front and rear of the valve body. This paper also created the connection between the pressure before valve p1, valve body opening lre, pipe flow q, and the relationship between the valve opening degree lre and the number of driving pulses np of the stepper motor. The correlation of the measured data under test circumstances is more than 0.99, indicating excellent fitting performance.
- (3)
A control program for adjustment and the flow measurement is designed based on the adjustment characteristics of the valve body. The test shows that the flow adjustment effect is excellent, with a maximum adjustment time ts of 78 s, an overshoot Δh of less than 5%, a steady-state error ess of less than 6%, and a better adjustment effect with increased pressure. The measured flow and the actual flow are also well matched, with an error of less than 3%.
- (4)
Compared to other relevant devices on the market, the device is simple in construction, less costly, and the regulation meets production needs. Subsequent research will be based on the application of this valve body in specific irrigation scenarios.
Author Contributions
Conceptualization, D.Z., M.L., H.Z. and X.N.; methodology, M.L., C.S. and H.T.; software, M.L., H.T. and C.S.; validation, B.J. and H.F.; resources, D.Z. and H.Z.; data curation, B.J. and H.F.; writing—original draft preparation, M.L.; writing—review and editing, D.Z., H.T. and H.Z.; visualization, X.N.; project administration, D.Z.; funding acquisition, D.Z. All authors have read and agreed to the published version of the manuscript.
Funding
This work was funded by Ningxia Hui Autonomous Region Key R&D Program, grant number 2022BBF02026; National Key R&D Program Project, grant number 2021YFE0103000 and Shaanxi Province Key R&D Program, grant number 2020ZDNY01-01.
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Informed consent was obtained from all subjects involved in the study.
Data Availability Statement
Raw data are not publicly available, though the data may be made available on request from the corresponding author.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Abioye, E.A.; Abidin, M.S.Z.; Mahmud, M.S.A.; Buyamin, S.; Ishak, M.H.I.; Rahman, M.K.I.A.; Otuoze, A.O.; Onotu, P.; Ramli, M.S.A. A review on monitoring and advanced control strategies for precision irrigation. Comput. Electron. Agric. 2020, 173, 105441. [Google Scholar] [CrossRef]
- Yang, M. Development Status of Agricultural Machinery Automation in my country. China Collect. Econ. 2021, 26, 21–22. [Google Scholar]
- Li, H.; Tang, P.; Chen, C.; Zhang, Z.; Xia, H. Research status and development trend of fertilization equipment used in fertigation in China. J. Drain. Irrig. Mach. Eng. 2021, 39, 200–209. [Google Scholar] [CrossRef]
- Sui, R.; O’Shaughnessy, S.; Evett, S.R.; Andrade-Rodriguez, A.; Baggard, J. Evaluation of a decision support system for variable-rate irrigation in a humid region. Trans. Asabe 2020, 63, 1207–1215. [Google Scholar] [CrossRef]
- Sui, R.; Yan, H. Field Study of Variable Rate Irrigation Management in Humid Climates. Irrig. Drain. 2017, 66, 327–339. [Google Scholar] [CrossRef]
- Mo, J.; Huang, X.; Li, W.; Li, Y. Research and optimization of hydraulic characteristics of large-scale variable sprinkler irrigation machine based on PWM technology. Trans. Chin. Soc. Agric. Eng. 2020, 36, 10. [Google Scholar] [CrossRef]
- Salcedo, R.; Zhu, H.; Zhang, Z.; Wei, Z.; Chen, L.; Ozkan, E.; Falchieri, D. Foliar deposition and coverage on young apple trees with PWM-controlled spray systems. Comput. Electron. Agric. 2020, 178, 105794. [Google Scholar] [CrossRef]
- Li, J.; Hong, T.; Feng, R.; Yue, X.; Luo, Y. Design and experiment of Venturi variable fertilizer apparatus based on pulse width modulation. Trans. Chin. Soc. Agric. Eng. 2012, 28, 6. [Google Scholar] [CrossRef]
- Chen, C.; He, P.; Zhang, J.; Li, X.; Ren, Z.; Zhao, J.; He, J.; Wang, Y.; Liu, H.; Kang, J. A fixed-amount and variable-rate fertilizer applicator based on pulse width modulation. Comput. Electron. Agric. 2018, 148, 330–336. [Google Scholar] [CrossRef]
- Zhu, H.; Sorensen, R.; Butts, C.; Lamb, M.; Blankenship, P. A pressure regulating system for variable irrigation flow controls. Appl. Eng. Agric. 2002, 18, 533–540. [Google Scholar]
- Wang, L.; Zhang, S.; Ma, C.; Xu, Y.; Qi, J.; Wang, W. Design of variable spraying system based on ARM. Trans. Chin. Soc. Agric. Eng. 2010, 26, 6. [Google Scholar] [CrossRef]
- Sun, C. Theoretical and Experimental Research on Valve Regulating Proportional Ferilizer Injector. Master’s Thesis, Jiangsu University, Zhenjiang, China, 2018. [Google Scholar]
- Song, L.; Dong, Z.; Xiang, L.; Xing, S. Variable universe adaptive fuzzy PID control of spray flow valve. Trans. Chin. Soc. Agric. Eng. 2010, 26, 114–118. [Google Scholar] [CrossRef]
- Wang, H.; Wang, X.; Huang, J.; Quan, L. Flow Control for a Two-Stage Proportional Valve with Hydraulic Position Feedback. Chin. J. Mech. Eng. 2020, 33, 1–13. [Google Scholar] [CrossRef]
- Zhuang, Y.; Zhang, N.; Ma, L.; Huang, Y. Design of flowrate control system for liquid flow standard device. China Meas. Test 2020, 46, 126–131. [Google Scholar] [CrossRef]
- Ko, Y.; Kim, T. Feedforward Plus Feedback Control of an Electro-Hydraulic Valve System Using a Proportional Control Valve. Actuators 2020, 9, 45. [Google Scholar] [CrossRef]
- Li, J.; Jia, W.; Wei, X. On-line Mixing Pesticide Device Based on Flow Control Valve and Neural Network. Trans. Chin. Soc. Agric. Mach. 2014, 45, 98–103. [Google Scholar] [CrossRef]
- Do-Hun, C. A Study on the Numerical Analysis of Internal Flow in a Cone Type Valve. J. Korean Soc. Ind. Converg. 2020, 23, 199–207. [Google Scholar]
- Zaryankin, A.E.; Zaryankin, V.A.; Akatov, A.S.; Zonov, A.S. Development and Investigation of a New Rotary Valve for Power Steam Turbines. Therm. Eng. 2020, 67, 249–255. [Google Scholar] [CrossRef]
- Rezazadeh, P.; Bijankhan, M.; Mahdavi Mazdeh, A. An experimental study on a flow control device applicable in pressurized networks. Flow Meas. Instrum. 2019, 68, 101533. [Google Scholar] [CrossRef]
- Zhang, J.; Wang, D.; Xu, B.; Su, Q.; Lu, Z.; Wang, W. Flow control of a proportional directional valve without the flow meter. Flow Meas. Instrum. 2019, 67, 131–141. [Google Scholar] [CrossRef]
- Li, J.; Chen, S.; Xiao, S.; Wang, C.; Li, M.; Xiong, W.; Li, K.; Liu, C. Application of Internet of Things Technology in Precision Agriculture. China Sci. Technol. Inf. 2022, 4, 84–86. [Google Scholar] [CrossRef]
- Masseroni, D.; Arbat, G.; de Lima, I.P. Editorial-Managing and Planning Water Resources for Irrigation: Smart-Irrigation Systems for Providing Sustainable Agriculture and Maintaining Ecosystem Services. Water-Sui 2020, 12, 263. [Google Scholar] [CrossRef] [Green Version]
- Lu, T. Research on Intelligent Water and Fertilizer Integration Irrigation System Based on Internet of Things. Master’s Thesis, North China University of Water Resources and Hydropower, Zhengzhou, China, 2019. [Google Scholar]
- Zhu, D.; Ruan, H.; Wu, P.; Li, J.; Lu, L. Remote fuzzy PID control strategy for fertilizer conductivity of water and fertilizer machine. J. Agric. Mach. 2022, 53, 186–191. [Google Scholar]
- Li, W.; Clark, B.; Taylor, J.A.; Kendall, H.; Jones, G.; Li, Z.; Jin, S.; Zhao, C.; Yang, G.; Shuai, C.; et al. A hybrid modelling approach to understanding adoption of precision agriculture technologies in Chinese cropping systems. Comput. Electron. Agric. 2020, 172, 105305. [Google Scholar] [CrossRef]
- Feng, X.; Li, B.; Han, S. Mechanical Principles and Mechanical Design, 2nd ed.; Higher Education Press: Beijing, China, 2014; pp. 147–175. [Google Scholar]
- Lu, P. Practical Technology of Control Valve, 2nd ed.; Mechanical Industry Press: Beijing, China, 2017; pp. 123–125. [Google Scholar]
| Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).