1. Introduction
Microsatellites have been actively developed [
1,
2]. However, except in the case that a large number of microsatellites are launched into the same orbit for constellation flight, these satellites are launched piggyback riding along with large satellites or released from the International Space Station; thus, they have no choice of schedule or orbit. To launch these microsatellites quickly and at low cost, the authors have been developing of a micro hybrid rocket since 2014 and have conducted launch experiments every year until 2018. Hybrid rocket is consisted by liquid or gaseous oxidizer and a solid fuel like plastics or paraffins. Because the physical partition of fuel and oxidizer, hybrid rockets have high resistance to explosion hazards and can achieve safer and more reliable space transport systems compared to conventional liquid and solid rockets. Moreover, the safer character can realize low costs on ground manipulation on fuels and rockets, such as storage, production, transportation. Thus, many studies on hybrid rocket were conducted.
Because of a plastic grain like Acrylonitrile-Butadiene-Styrene (ABS) can be used for the fuel of hybrid rocket, many studies for 3-D Printed plastic grain. Lyne et al. developed powder fuel contained 3-D printed fuels to improve regression rate of the grain [
30]. McFarland and Antunes conducted static combustion test 3D-printed grain made from ABS, acrylonitrile styrene acrylate, polylactic acid (PLA), polypropylene, polyethylene terephthalate glycol, Nylon, and PLA with aluminum particles to investigate the regression rate of the grain [
31]. Wang et al. investigated combustion performance of fuels comprised a paraffin-based fuel embedded on a 3-D printed nested helical ABS grain [
32].
Paraffins can also be used for the grain. Nakagawa et al. investigate a regression rate of five-type cylindrical port paraffin fuels with oxidizer swirl flow and concluded that the regression rate can be increased with oxidizer swirl flow [
33]. Tang et al. investigated mechanical modified paraffin-based grain on combustion test to obtain an excellent paraffin-based grain with both good mechanical properties and combustion performance [
34]. In 2017, the authors developed the 3D-printed solid fuel grain with a star-fractal port to improve the regression rate and used in the launch experiment. In 2018, the authors added swirl shape to a star-fractal port. The vehicle used star-fractal swirl port hybrid rocket engine launched in that year reached an altitude of 6.2 km [
8]. In 2019, a 50 kNs-class hybrid rocket engine was developed with the goal of reaching a maximum altitude of 15 km; however, an engine rupture [
9] occurred in combustion test hence the launch was abandoned. In 2020, COVID-19 prevented sufficient experiments from being conducted to determine the cause of the rupture failure, and the launch was abandoned again. In March 2021, the cause of the rupture was identified. An engine equipped with countermeasures successfully completed five combustion tests from May to July 2021, and a rocket named “Amanuma Maru” was launched at 6:00 am on September 19, 2021 from Noshiro City, Akita Prefecture. Atmospheric pressure telemetry data indicated that it reached a maximum altitude of 10.1 km. This paper reports a summary of its development and launch test.
2. Outline of the Rocket and the Launcher
Figure 1 is an overview of the rocket. Hybrid rocket propulsion has the propellant in different phase, i.e., liquid oxidizer and solid fuel, and has an essentially non-explosive nature due to its boundary layer combustion mechanism [
10]. The reason for employing a hybrid rocket propulsion as our rocket engine is its high safety. The rocket was designed and developed in 2019 [
11] and the launch was scheduled for 2019. However, several ruptures of the engine prevented the launch.
An improved engine was subsequently developed and fitted.
Table 1 shows its mass properties in comparison with those of CAMUI-500P [
12], which held the altitude record for hybrid rockets in Japan, and our FY2018 model, “Fukushima-Maru”. The object of the development was to launch microsatellites quickly and at low cost; thus, the method of manufacture used typical processing methods such as lathe processing and the vehicle was designed to be lightweight.
2.1. Engine
There is a technical issue in hybrid rocket propulsion – low fuel regression rate. This can cause low levels of thrust. To overcome this issue, many research and developments were conducted. For example, Fuller et al. developed and studied a grain with a star swirl port [
13,
14]. Whitmore et al. developed and studied a grain with a circular helical port [
15,
16]. Yuasa et al. studied the improvement of the regression rate of hybrid rockets by a swirling oxidizer flow [
17,
18]. The authors also developed a 3D-printed solid fuel grain with a grain with a star-fractal [
19], cherry-blossom-fractal [
19] and a star-fractal swirl port was developed [
20].
Figure 2 shows the cross-section of the port has star-fractal geometry, and it is swept while rotating axially. The purposes of the star-fractal swirl port are to have a large surface area and to generate a swirling oxidizer flow. The effectiveness of this port was demonstrated in ground combustion experiments [
21], the improvement of the regression-rate of the fuel was evaluated [
22], and hence this port was adopted for the launch experiment.
The fuel grain was made by a fused-deposition-modeling-type, commercial-grade 3D printer. The filament material of the 3D printer was an acrylonitrile butadiene styrene (ABS) resin. Because 3D printing was a time-consuming process, only portions near the port were 3D-printed. Commercially available ABS pipes were used as other portions of the fuel grain. This also led to cost reductions.
The fuel grain was covered by a glass-fiber-reinforced-plastics (GFRP) ablator and inserted into a motor case. There were two types of motor cases, such as the external thread type and the welded flange type. These motor cases were different only for the fastening method between the motor case and the injector, and between the motor case and the nozzle housing. For these motor case, however, their function as a combustion chamber was the same by design. The oxidizer, nitrous oxide (N2O), was injected into the combustion chamber through the injector. The bell nozzle was made by graphite and had a 31 mm throat.
The development of the engine was started in 2019. The preliminary model was designed, developed, and tested. However, several ruptures occurred in the combustion test. Their cause was determined to be self-decomposition of the oxidizer, nitrous oxide [
23] and insufficient strength of the resin housing of the nozzle. Improvements made to both the oxidizer and housing prevented the rupture.
Figure 3 shows combustion test facility, which was also improved to have enough strength to prevent failure on engine rapture.
Figure 4 shows the thrust and time curve of the five successful combustion tests, and
Table 2 confirms the reproducibility of the total impulse (integration of thrust and time curve), which was within -5.0% to +5.1% of the average.
2.2. Oxdizer Tank
The oxidizer tank is based on the one developed in 2018 [
24]. The total length was extended to 2 m to increase the carrying capacity. Due to extended length and the difference in the coefficient of thermal expansion between the aluminum alloy body and the carbon fiber reinforced plastic (CFRP) that reinforces it, the CFRP reinforcement part cracked as a result of the temperature drop during the filling of the tank with the oxidizer. The part was subsequently improved by changing the lamination sequence. As a result, the weight index
pV/(
mg) (maximum internal pressure × volume / tank weight) was improved 1.7 times, i.e., from 1.4 to 2.4.
2.3. Separation Mechanism
As shown in
Figure 1, the rocket had separation mechanisms on the first and second stages. The purpose the separation mechanism of the first stage is to release a streamer to slightly decelerate the rocket at the top of its flight, while that of the second stage is to release the parachute and deployable floats to further decelerate the rocket and land on the sea. The first-stage separation mechanism used a Nichrome wire to cut Reny (high performance polyamide based on mainly polyamide MXD6) bolts, which is a type of polyamide resin. This mechanism had been used in launches in 2016, 2017, and 2018. The second-stage separation mechanism used a non-explosive low-shock separation mechanism [
25] with a separation nut.
2.4. Telemetory & Datalogger
Telemetry measured air pressure and GPS position coordinates during flight, and it downlinked to the ground by using LoRa (Long Range) communication protocol, a type of LPWA (Low Power Wide Area) using 920 MHz-band. The data logger recorded information, such as acceleration, geomagnetism, and gyro, and was picked up when the rocket was recovered. Both units were mounted on the measurement module shown in
Figure 1. Because the altitude of the rocket can be estimated from the air pressure during flight, it can be measured in real time from the downlink signal. The GPS signal was used to identify the position of the rocket when it landed on sea. The telemetry and data logger had initially been developed at Kanagawa University and flew in 2017 [
26] and 2018 [
27]; they were made more compact, lightweight and reliable by Fullheart Japan Company [
28].
2.5. Body
The strength design was taken into account for a compression load by thrust load of rocket engine, bending load by wind gust, and shock load by parachute opening. To achieve ultra-lightweight, the body was constructed by thin-walled, composite material cylinders. The critical failure load of the thin-walled cylinders is buckling due to compression load. It is well-known that the experiment buckling load has large scatter and is smaller than the theoretical buckling load. The ratio, the experimental buckling load divide by the theoretical buckling load, is called “Knockdown factor”. Many efforts to find the cause of the difference between the experimental buckling load and the theoretical buckling load have been conducted. The cause, however, is not revealed. One of the authors, Takano, derived the theoretical buckling load with less approximation for composite cylinder [
29,
30] and, obtained the statistical knockdown factor with confidence bound [
31]. The method by the theoretical buckling load with the statistical knockdown factor was applied for the design of the body hence it achieved ultra-light weight body. Divergence (static divergence due to aeroelasticity) of the fuselage was also taken into consideration in the design of the bending stiffness fuselage and the joints (couplers) between the stages. CFRP was used for the body material that achieves high-bending stiffness, and glass fiber reinforced plastic (GFRP) was used for the body of the equipment that generated radio waves, such as the measurement module. The coupler was fixed to the fuselage with shear bolts. However, because it is difficult to predict the rigidity reduction of the joint, the bending stiffness of the fuselage was measured to determine the stiffness reduction of the joint and it was used to evaluate the divergence. An internal resonance phenomenon, which was studied by Wang and Chan [
32], was not considered. Thus, the consideration is the future work.
2.6. Fin & Finstay
Flutter, which is an aeroelastic phenomenon, was also considered in the design of the fins and fin stays (L-clip that fixes the fins to the body). The fins had a diamond airfoil cross section and an aluminum honeycomb sandwich panel with a CFRP skin to improve stiffness. To avoid a decrease in stiffness in the fin stay, which is the joint between the fin and the body, an L-clip made of CFRP was used and both adhesion and peeling bolt joints were used.
2.7. Launcher [33]
The launcher had a glide length of 12 m with elevation angle from 60° to 86° and was designed for the maximum rocket weight of 100 kg and rocket length of 5 m. The maximum elevation angle error due to elastic deformation by the 150kg weight of rocket was designed less than 0.55°. It was designed to be divided into pieces of 3.6 m or less in length to enable it to be moved to the launch site, with most parts being 20 kg or less in weight; thus it could be assembled without the need for heavy machinery. Design load for the launcher was not only the rocket weight but also wind loads. Design load by wind for a crane design standard [
34] was applied. For the operational state, which was elevation angle with 80° and set the rocket, 16m/s was applied and the minimum margin of safety (MS) was 5.5. For non-operational state, which was elevation angle 0° and without the rocket, 56m/s was applied and the minimum MS was 0.45.
3. Launch Test and Results
A launch test was conducted at 6:00 am on September 19, 2021 at the former Ochiaihama in Noshiro City, Akita Prefecture, for the purpose of a flight demonstration and attempting to break the Japanese altitude record for hybrid rockets.
Figure 5 shows the state of launch.
The oxidizer filled by GSE (Ground Support Equipment) shown in
Figure 3. without freight container, engine stand and linear motion guide.
From the visual observation by telescope camera and the altitude data described later, it was judged that the combustion of the engine was normally complete. Altitude data from the atmospheric pressure telemetry was also successfully received. The launch was planned to be carried out on any day with good weather from September 18th to 20th, 2021, including the backup date. However, due to the effects of Typhoon No. 14, which took a peculiar course of going from western Japan to eastern Japan after staying in the East China Sea, the launch was postponed on the 18th as it was judged that the sea swells were high and the navigation and recovery work of the recovery ship would be difficult. On the 19th, it was forecast that swells would remain on the sea the day before, but there was a risk that the launch on the 20th would fail due to problems with the rocket and GSE. Therefore, a final decision was made based on the weather at 4:00 AM on the 19th. Although the land was covered with fog, the surface of the sea could be seen from the vicinity of the launch point, and so it was decided to proceed with the launch. After inserting the rocket into the launcher, it became possible to turn on the power of the on-board camera; however, there was a problem that recording could not be started with the remote control.
Condensation on the day of the flight was severe and the video cameras on the ground also had condensation inside the camera lens. Thus, the condensation was judged to be the cause of the problem. If the recording of the camera mounted on the rocket could not start, it would not be possible to shoot videos during flight. However, it was judged that the purpose of the experiment could be achieved even if the camera recording was not possible, so it was decided to launch the rocket. As shown in
Table 3, just before launch, one of the two atmospheric pressure telemetry units and one GPS telemetry unit were out of communication. However, previous experience had indicated that the position of the antenna and the position of the receiving station could be improved depending on the attitude of the rocket, so the launch proceeded. An investigation after the flight indicated a possible cause was that the GPS and atmospheric pressure telemetry had been unable to be received before the launch because of the carrier sense function, which stops radio waves when radio waves are emitted in the same frequency band near the transmitter to avoid interference.
GPS reception became impossible after the launch, but this seems to be due to restricted GPS reception for high-speed moving objects. However, there was also a problem that communication with the forced separation circuit became impossible, and it was not possible to confirm the voltage supply signal to the Nichrome wire for separation or receive the separation signal.
Investigations conducted after flight determined the cause to be forced separation circuit failure due to condensation or low temperatures under 0 ℃ at high altitude.
However, only atmospheric pressure telemetry was received.
Figure 6 shows the result of estimating the altitude from the atmospheric pressure data.
Altitude was estimated using the height measurement formula (Equation (1)) derived from the ICAO (International Civil Aviation Organization) standard atmosphere. Note that the subscript 0 means the ground value. As a result of calculation using the pre-launch pressure
p0 =1014.42 hPa, temperature
T0 =13 ℃, and the minimum measured atmospheric pressure
p0 =256.85 hPa, the maximum attainable altitude was calculated to be 10124 m. This means the flight broke the previous altitude record of 8.3 km (estimated value) held by CAMUI-500P.
Figure 5 also shows the altitude simulation results. The data after the apex show ballistic fall (fall without streamer and parachute deployed) and deceleration fall (fall with streamer and parachute deployed). Considering this result and the fact that the separation signal could not be received, it is thought that the first stage did not separate normally and fell, landing on the water before the second stage separation. As a result of subsequent investigations, the cause was determined that a forced separation circuit failure had occurred due to dew condensation and a drop in temperature at high altitudes.
Velocities also evaluated by the altitude simulation results and numerical differentiation of measured height. The numerical differentiation is second-ordered, using three measured point shown in Equation (2).
Figure 7 shows the measured and simulated velocity. The maximum velocities of height direction are 309 m/s (measured) and 290 m/s (simulated). The possibility cause of the difference is prediction accuracy of drag coefficient (
Cd) of the rocket.
Cd was predicted by the following step. Calculating
Cd on the subsonic (Mach=0.8) and then predicting
Cd change by Mach number (Mach-
Cd) was using data of similar rocket. Thus, the
Cd around Mach=1 has less accuracy. Improvement of the Mach-
Cd is the future work.
4. Conclusion
A hybrid rocket was launched and achieved an altitude of 10.1 km, setting a new domestic altitude record for hybrid rockets. However, the confirmation of the separation signal, the reception of the GPS signal, and the recovery of the rocket failed, and the data of the on-board data logger could not be recovered. Furthermore, it is presumed that the first stage did not separate normally either. Although the domestic altitude record for hybrid rockets was broken, many problems remained. An investigation conducted after flight revealed the cause of the failure of the telemetry and forced separation circuits. In the future, we would like to take countermeasures against these problems and improve the attainable altitude.
Author Contributions
Conceptualization, At.T.; methodology, At.T. and F.Y.; software, R.K., Y.K. and M.M.; investigation, At.T., K.Y. Y.F. and R.K.; data curation, X.X.; writing—original draft preparation, At.T.; writing—review and editing, At.T., K.T., Ak.T., T.M. and S.U.; project administration, At.T. All authors have read and agreed to the published version of the manuscript.
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Not applicable.
Data Availability Statement
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the fact they are not open data.
Conflicts of Interest
Author Yoshihiko Kunihiro and Makoto Miyake were employed by the company Fullheart Japan Corporation. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
- Microsatellites from UNISEC, UNISEC homepage, http://unisec.jp/unisec/satellites, (accessed on 13 Feberuary 2024).
- Shou, H.N. Micro-Satellite Configuration of Discoid and Asymmetrical, Gyroless with Thrusters Three-Axis Robust Control and Stability Analysis, International Journal of Engineering and Technology Innovation, 2014, 1, 1–17.
- Lyne, J. E.; Brigham, A.; Savery, R.; Karcher, K.; Pyron, J.; Adams, L.; Reagan, G.; Furches, H.; Sola, D.; Melendez, L.; Keck, C. The Use of a 3-D Printed, Polymer Matrix Containing Pulverized Fuel in a Hybrid Rocket, AIAA Propulsion and Energy Forum. 2018. [Google Scholar] [CrossRef]
- McFarland, M.; Antunes, E. Small-Scale Static Fire Tests of 3D Printing Hybrid Rocket Fuel Grains Produced from Different Materials, Aerospace, 2019, 6, 81.
- Wang, Z.; Lin, X.; Li, F.; Yu, X. Combustion performance of a novel hybrid rocket fuel grain with a nested helical structure, Aerospace Science and Technology, 2020, 97, 105613.
- Nakagawa, I.; Hikone, S.; Ishiguro, T.; Maruyama, S. A Study on Paraffin-based Fuel Hybrid Rockets, Trans. JSASS Aerospace Tech. Japan, 2012, 10, ists28, To_1_7–To_1_10. [Google Scholar]
- Tang, Y.; Chen, S.; Zhang, W.; Shen, R.; DeLuca, L. T.; Ye, Y. Mechanical Modifications of Paraffin-based Fuels and the Effects on Combustion Performance. Propellants Explosives Pyrotechnics, 2017, 42, 1268–1277. [Google Scholar] [CrossRef]
- Takano, A.; Funami, Y.; Nishino, S. Conceptual design and launch of hybrid rocket with star fractal swirl fuel grain. UNISEC Space Takumi Journal, 2020, 9, 1–16. (in Japanese). [Google Scholar]
- Igarashi, H.; Takano, A. ; Kitamura, R, Funami, Y. Burst case of hybrid rocket engine with nitrous oxide and countermeasures. 10th UNISEC Space Takumi Conference, 2020, UNISEC 2020-004 (in Japanese).
- Takahashi, A.; Shimada, T. Essentially non-explosive propulsion paving a way for fail-safe space transportation. Transactions of the JSASS, Aerospace Technology Japan, 2018, 16, 1–8. [Google Scholar] [CrossRef]
- Bao, J.; Kitamura, R. ; Takano, A, Funami Y, “Conceptual design of a hybrid rocket toward an altitude of 15km. 2nd Hybrid Rocket Symposium, 2019, HR-2019-006 (in Japanese).
- Nagata, H.; Wakita, M.; Totani, T.; Uematsu, T. Development and Flight Demonstration of 5kN Thrust Class CAMUI Type Hybrid Rocket. Transactions of the JSASS, Aerospace Technology Japan, 2014; 12, Ta 1–Ta 4. [Google Scholar]
- 13. Fuller, J-K.; Ehrlich, D.A.; Lu, P. C<bold>.,</bold> Jansen, R.P.; Hoffman, J. D. Advantages of Rapid Prototyping for Hybrid Rocket Motor Fuel Grain Fabrication. AIAA Paper, 2011-5821. 2011.
- Armold, D.; Boyer, J. E.; Kuo, K. K.; DeSain, J. D.; Curtiss, T. J.; Fuller, J. K. Test of Hybrid Rocket Fuel Grains with Swirl Patterns Fabricated Using Rapid Prototyping Technology. AIAA Paper, 2013-4141. 2013. [Google Scholar]
- Whitmore, S.A.; Walker, S.D.; Merkley, D.P.; Sobbi, M. High Regression Rate Hybrid Rocket Fuel Grains with Helical Port Structures. J. Propul. Power, 2015, 32, 1727–1738. [Google Scholar] [CrossRef]
- Whitmore, S.A.; Walker, S. D. Engineering Model for Hybrid Fuel Regression Rate Amplification Using Helical Ports. J. Propul. Power, 2017, 33, 398–407. [Google Scholar] [CrossRef]
- Yuasa, S.; Shimada, O.; Imamura, T.; Tamura, T.; Yamamoto, K. A Technique for Improving the Performance of Hybrid Rocket Engines. AIAA Paper, 1999-2322. 1999. [Google Scholar]
- Yuasa, S.; Shiraishi, N.; Hirata, K. Controlling Parameters for Fuel Regression Rate of Swirling-Oxidizer-Flow-Type Hybrid Rocket Engine. AIAA Paper, 2012-4106. 2012. [Google Scholar]
- Tateyama, T.; Takano, A. Development of CFRP Reinforced Lightweight Hybrid Rocket Engine. Aerospace Technology Japan, 2018, 17, 237–243. (in Japanese). [Google Scholar] [CrossRef]
- Shizume, M.; Takano, A.; Funami, Y.; Morohoshi, H.; Tahara, K.; Terada, T. Development of Hybrid Rocket Engine with Star Fractal Swirl Grain Made by 3D Printer. Space Transportation Symposium FY2017, Sagamihara, Japan, STPC-2017-012, 2018 (in Japanese).
- Funami, Y.; Takano, A. Ground combustion experiments of hybrid rocket fuel grain with a star fractal swirl port. 63rd Space Sciences and Technology Conference, 2E07 (JSASS-2019-4316), 2019 (in Japanese).
- Funami, Y.; Takano, A. Regression-Rate Evaluation of Hybrid-Rocket Fuel Grain with a Star-Fractal Swirl Port. Trans. Japan Soc. Aero. Space Sci. 2023, 66, 61–69. [Google Scholar] [CrossRef]
- Hamazaki, A.; Kaneyori, H.; Masu, K.; Kumada, K.; Funami, Y.; Kitamura, R.; Takano, A. Investigation of decomposition on nitrous oxide for hybrid rocket engine. Proceedings of Space Transportation Symposium FY2020, STCP-2019-013, 2021 (in Japanese).
- Tateyama, T.; Takano, A. Development of Low Cost/Light Weight Oxidizer Tank in Hybrid Rocket. Aerospace Technology Japan, 2020, 19, 19–26. (in Japanese). [Google Scholar] [CrossRef]
- Takano, A.; Nishino, S. Development of Separation Mechanism Using Non-Explosive Separation Nuts. Aerospace Technology Japan, 2020, 19, 179–185. (in Japanese). [Google Scholar] [CrossRef]
- Shimazaki, T.; Takano, A. Development of Telemetry device and Launch result at Izu Oshima. 8th UNISEC Space Takumi Conference, UNISEC 2018-003, 2018(in Japanese).
- Yoshida, M.; Shimazaki, T.; Takano, A. Development and flight test result of telemeters for hybrid rocket. Space Transportation Symposium FY2018, STCP-2018-024, 2019 (in Japanese).
- Hoshi, T.; Takano, A.; Miyake, M.; Kunihiro, Y. Development of independent mounting type Telemetry device for micro launch vehicle. 10th UNISEC Space Takumi Conference, UNISEC 2020-003, 2020 (in Japanese).
- Takano, A. Improvement of Flügge’s Equations for Buckling of Moderately Thick Anisotropic Cylindrical Shells, AIAA Journal, 2008, 46, 903–911.
- Takano, A. Simple Closed-Form Solution for the Buckling of Moderately Thick Anisotropic Cylinders, Trans. JSASS, Aerospace Technology Japan, 2012, 10, 17–26. [Google Scholar] [CrossRef]
- Takano, A. Statistical Knockdown Factors of Buckling Anisotropic Cylinders under Axial Compression, ASME Journal of Applied Mechanics. 2012; 79. [Google Scholar] [CrossRef]
- Wang, Y.R.; Chang, Y.S. Study of Primary and Internal Resonance on 3D Free-Free Double-Section Beam. Adv. technol. innov., 2020, 5, 270–291. [Google Scholar] [CrossRef]
- Tachibana, Y.; Yamaguchi, Y.; Takano, A. Development of 12 m launcher for 5 m hybrid rocket launch. Proceedings of Space Transportation Symposium FY2018, TCP-2018-025, 2019 (in Japanese).
- Japanese Industrial Standards Committee, Cranes-Wind load assessment, JIS B 8830, 2001. +5.
|
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. |
© 2024 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/).