Design and Optimization Strategies of a High-Performance Vented Box

Podsumowanie

Here, we present the range analysis method to optimize the sample points generated by an orthogonal experimental design to ensure that fresh food can be stored in a vented box for a long time by regulating the airflow pattern.

Streszczenie

This study aims to solve the problems of air flow chaos and poor performance in a vented box caused by the heterogeneous distribution of airflow through the design of the internal structure of the vented box with constant energy consumption. The final goal is to evenly distribute the airflow inside the vented box. Sensitivity analysis was performed for three structural parameters: the number of pipes, the number of holes in the middle pipe, and the number of each increment from the inside to the outside pipe. A total of 16 random array sets of three structural parameters with four levels were determined using the orthogonal experimental design. Commercial software was used for the construction of a 3D model for the selected experimental points, and this data was used to obtain the airflow velocities, which were then used to obtain the standard deviation of each experimental point. According to the range analysis, the combination of the three structural parameters was optimized. In other words, an efficient and economical optimization method considering the performance of the vented box was established, and it could be widely used to extend the storage time of fresh food.

Wprowadzenie

Fresh vegetables and fruits occupy a high proportion of human food consumption, not only because they have good taste and an attractive form, but also because they are of great benefit for people to obtain nutrition and maintain health¹. Many studies have shown that fresh fruits and vegetables play a unique role in preventing many diseases^2,3. In the storage process of fresh fruits and vegetables, fungi, light, temperature, and relative humidity are the important reasons for their deterioration^4,5,6,7,8. These external conditions affect the quality of stored fresh fruits and vegetables by affecting the internal metabolism or chemical reactions⁹.

Common treatment technologies for fruits and vegetables include nonthermal and thermal preservation. Among them, thermal pretreatment has a positive effect on the drying process, but it can also have adverse effects on product quality, such as loss of nutrients, change of flavor and odor, and change of color^10,11. Therefore, in recent years, the nonthermal preservation of products has received attention from the research perspective to meet the demand of consumers for fresh products. At present, there are mainly radiation processing, pulsed electric field, ozone processing, edible coatings, dense phase carbon dioxide, and other nonthermal preservation technologies to store fruits and vegetables, but these technologies often have shortcomings, such as the requirement of large equipment, high price, and the cost of use¹². Therefore, the design of a simple structure, low cost, and convenient control of the preservation equipment is very meaningful to the food industry.

In the storage environment for fruits and vegetables, a proper air circulation system helps to eliminate the heat generated by the product itself, reduce the temperature gradient, and maintain the temperature and humidity in the space where it is located. Proper air circulation also prevents weight loss due to respiration and fungal infections^13,14,15. Numerous studies have been conducted on airflow within different structures. Praeger et al.^16,17 measured the wind speed at different positions under different fan operating powers in a warehouse through sensors and found that there could be as big as a sevenfold difference in air velocity due to different vertical heights, and the air velocity at each position was positively correlated with the fan operating power. Moreover, a study examined the effect of cargo arrangement and the number of fans on airflow, and it was concluded that increasing the distance of some fan positions and rationally choosing the number of fans was helpful in improving the effect. Berry et al.¹⁸ studied the effect of airflow in different fruit storage environments on stomata distribution in packing boxes. Using simulation software, Dehghannya et al.^19,20 studied the air flow state of forced pre-cold air in the package with different vent areas, quantities, and distribution positions on the packaging wall, and obtained the nonlinear influence of each parameter on the air flow state. Delele et al.²¹ applied a computational fluid dynamics model to study the influence of products randomly distributed in different forms of ventilation boxes on airflow. They found that the product size, porosity, and box hole ratio had a greater impact on airflow, whereas random filling had a smaller impact. Ilangovan et al.²² studied airflow patterns and thermal behavior between the three packaging structures and compared the results with reference structural models. The results showed that the heat distribution in the box was not uniform because of the different locations and designs of the vent. Gong et al.²³ optimized the width of the gap between the edge of the tray and the wall of the container.

The techniques used in this paper include simulation and optimization methods. The principle of the former is that the governing equations were discretized and numerically solved using the finite volume method²¹. The optimization method used in this paper is referred to as orthogonal optimization²⁴. The orthogonal test is a typical multifactor and multilevel analysis method. The orthogonal table built using this method contains representative points uniformly distributed in the design space, which can visually describe the entire design space and be examined. That is, fewer points represent the full factor test, greatly saving time, manpower, material, and financial resources. The orthogonal test has been widely used in the design of experiments in the fields of power systems, chemistry, civil engineering, etc²⁵.

The objective of this study is to design and optimize a high-performance vented box. A vented box can be defined as an original box including a gas control device that disperses the gas uniformly in the box. Velocity uniformity refers to how evenly air flows through the vented box. Yun-De et al.²⁶ have previously shown that the property of multiporous material has an important effect on the velocity uniformity of a fresh vegetable box. In some experiments, a plenum or modulated chamber was left at both the top and the bottom of the test chamber to guarantee a homogenous distribution of either forced or induced air²⁷. The vented box designed in this paper contains arrays of pipes with zigzag holes. Controlling the airflow distribution in the vented box is the main preservation strategy. There are two air inlets of equal size set parallelly at the left and right sides of the vented box, and an outlet is set at the upper side of the box. Designing the internal structure of a vented box is the key to this study. In other words, the number of pipes and holes is an important parameter for changing the internal structure of the vented box. The reference model has 10 pipes. The two middle pipes have 10 holes each, which are staggered across the pipes. The number of holes from the middle to the outer pipe increase by two at a time.

In other words, when we keep fresh vegetables, fruits, and other products, continuous and stable airflow can reduce the respiration of products, reduce ethylene and other harmful substances for product preservation, and reduce the temperature produced by the products themselves. Due to the different parameters of the vented box, it is not easy to obtain the required airflow state, which will affect the preservation property of the vented box. Therefore, the project takes the internal airflow velocity uniformity of the vented box as the control objective. A sensitivity analysis was conducted for the structural parameters of the vented box. The samples were selected by orthogonal experimental design. We used range analysis to optimize the combination of the three structural parameters. Meanwhile, we verify the desirability of the optimization results.

Protokół

1. Pre-simulation processing

NOTE: Considering the arrays of pipes, the three-dimensional bottom half and the top half of the vented box models are established by using three-dimensional software and saving them as X_T files, the overall dimensions are shown in Figure 1. Configurations are shown in the table of materials.

1. Run the simulation software and drag the Mesh component from the "Component Systems" to the "Project Schematic" window. Name it "the bottom". Right-click Geometry and click Browse to import the bottom X_T file.
2. Right-click Geometry, click New DesignModeler Geometry to enter into the "Mesh-Designmodeler" window, and click Generate to display the bottom model.
3. Right-click the upper surface and click Named Selection to rename it "vented box_upper". Select Selection Filter: Bodies. Right-click the bottom model to select the Named Selection and rename it as "the bottom".
4. Select Selection Filter: Faces and switch the "Select Mode" to Box Select. Select all the inner surfaces, right-click to select the Named Selection, and rename it as "inner surfaces_external", defined as mech interfaces later. Return to the initial window.
5. Double-click the bottom Mesh. Enter the "Meshing" window. Change the "Physical Preferences" from Mechanical to CFD. Click Update to generate the mesh model. Return to the initial window.
6. Drag the Meshcomponent from the "Component Systems" to the "Project Schematic" window. Name it "top". Right-click Geometry and click Browse to import the top X_T file.
7. Right-click Geometry and click New DesignModeler Geometry to enter into the "Mesh-Designmodeler" window. Click Generate to display the top model.
8. Right-click the lower surface and click Named Selection to rename it as "vented box_lower". Select Selection Filter: Bodies. Right-click the top model to select the Named Selection and rename it as "top".
9. Select Selection Filter: Faces. Right-click the upper surface and click Named Selection to rename it as the outlet. Return to the initial window.
10. Double-click the top's Mesh. Enter the "Meshing" window. Change the "Physical Preferences" from Mechanical to CFD. Right-click Mesh to select the Sizing in the "Insert". Select Selection Filter: Bodies. Select the top model and type 18 in "Element Size". Click Update. Return to the initial window.
11. Drag the Mesh component from the "Component Systems" to the "Project Schematic" window. Name it as the pipe. Import the pipe X_T file by clicking Geometry.
12. Enter the "Mesh-Designmodeler" window. Click Generate to display the pipe model displays again.
13. Select the two end faces of the pipe and label them as "inlet1" and "inlet2", and the pipe by body selecting and labeling as a pipe.
14. All inner surfaces by box selecting are labeled as "inner surfaces_internal", defined as mesh interfaces later. Return to the initial window.
15. DoubleClick the pipe's Mesh. Enter the "Meshing" window. Change the "Physical Preferences" from Mechanical to CFD. The mesh model can be generated by clicking "Update". Return to the initial window.
    NOTE: Figure 2A shows the grid of the bottom half of the vented box, Figure 2B shows the grid of the top half of the vented box, and Figure 2C shows the grid of the pipe. As shown in Figure 3, with the number of grids increasing from 4,137,724 to 5,490,081, the standard deviation changes are less than 0.0008. Considering the computation ability and accuracy, the following analysis is based on a grid model with 4,448,536 grids.

2. Simulation analysis

NOTE: The following operations are described based on the general sequence of simulation analysis from setup to solution to result.

1. Drag the Simulation component to the "Project Schematic" window. Link three Mesh components to the simulation component and update to enter.
2. Setup
   NOTE: The "Setup" consists of five steps: General, Models, Materials, Cell Zone Conditions, and Boundary Conditions.
   1. General: Verify the validity of the mesh model. Check whether the mesh has negative volume. Select Steady. For the relaxation, residual, and time scale factors, select the default values. If there is any problem with the divided grid or model settings, an error message will pop up.
   2. Models: Enter the setting interface of "Viscous Model" to select the K-epsilon Model.
   3. Materials: Set the "air" material.
   4. Cell Zone Conditions: Change the type of cell zone to Fluid.
3. Boundary Conditions
   1. Convert the type of vented box_upper, vented box_lower, inner surfaces_external, and inner surfaces_internal from the default "Wall" to "Interface".
      NOTE: Simulation software immediately generates the "Mesh Interfaces" after ending the above step.
   2. Open Mesh Interfaces and enter the "Create/Edit Mesh Interfaces" window. Match inner surfaces_external to inner surfaces_internal. Match vented box_upper to vented box_lower. Finally, the two mesh interfaces are created in the vented box and named interface1 and interface2, respectively.
   3. Set the airflow velocities of all inlets as 8.9525 m/s in the "Velocity Inlet" window and the gauge pressure of the outlet as zero in the "Pressure Outlet" window.
4. Solution
   1. Set the style of solution initialization as Standard Initialization before initializing.
   2. Set the Number of Iterations as 2,000.
   3. Click Calculate to start the simulation and return to the initial window until the simulation ends.
5. Results
   1. Click Results. Enter the "CFD Post" window and click the streamline icon of in the toolbox.
   2. Select outlet in "Start From" and Backward in "Direction". Click Apply to generate the internal flow diagram of the vented box.
   3. Click Plane in "Location", select ZX Plane in "Method", and select the input value as 0.6. Click Apply to generate the plane 0.6 m from the bottom surface.
   4. Click the Contour icon in the toolbox, select Plane 1 in "Locations", select Velocity in "Variable", and select Local in "Range". Click Apply to generate the velocity contour.
   5. Export the flow rate data for the plane generated above. Acquire the standard deviation of the flow rate in spreadsheet software (e.g., Excel).
      ​NOTE: Sensitivity analysis of three structural variables of the vented box was carried out. The number of pipes has four levels: eight, 10, 12, and 14. The number of holes in the middle pipes has four levels: eight, 10, 12, and 14. The number of each increment from the inside to the outside pipe has four levels: zero, two, four, and six. Modify the base model according to changes in these structural parameters. Repeat steps 1.1 through 2.5 10 times to obtain the data in Table 1. It can be seen from the table that the three structural parameters have certain effects on the standard deviation of the flow rate.

3. Orthogonal experiment design and range analysis

1. Run the statistical analysis software. Click Data and Generate in "Orthogonal Design".
2. Enter pipe_number into "Factor Name" and A in "Factor Label". Click Add and Define Values to set four levels for the number of pipes. Click Continue and return to the "Generate Orthogonal Design" window.
3. Enter whole_number into "Factor Name" and B in "Factor Label". Click Add and Define Values to set four levels for the number of holes. Click Continue and return to the "Generate Orthogonal Design" window.
4. Enter cumulative_number into "Factor Name" and C in "Factor Label". Click Add and Define Values to set four levels for the number of increments. Click Continue and return to the "Generate Orthogonal Design" window.
5. Click Create new data file to generate 16 array samples. Click Variable View to select Nominal in "Measure" and Input in "Role". Rename it as "standard_deviation×100000".
6. Repeat steps 1.1 to 2.5 with the sample points above; the resulting 16 standard deviations multiplied by 100,000 are filled into the sample list for later optimization.
7. Click Analyze and Univariate in "General Linear Model". Fill standard_deviation×100000 into "Dependent Variable" and pipe_number, hole_number, cumulative_number into "Fixed Factor(s)". Click Model and Build terms. Change Interaction to Main effects. Fill A, B, C into "Model". Click Continue and return to the "Univariate" window.
8. Click EM Means and fill A, B, C into "Display Means For". Click Continue and go back to the "Univariate" window.
9. Click OK and get the optimization result; the minimum value of the "mean" column in the table corresponds to the optimal variable. Double-click on the table, enter into the "Pivot Table" window, click Edit, and click Bar in "Create graph" to generate the histogram.

Wyniki

Following the protocol, the first three parts were the most important, which include modeling, meshing, and simulation, all in order to obtain the standard deviation of the flow rate. Then, we completed the structure optimization of the vented box through orthogonal experiments and range analysis. The model used in the protocol is the reference-vented box model, which is the initial model obtained from the reference. Figure 4 shows the result of the streamlined flow of the reference vented b...

Dyskusje

Due to its high performance and complex structure, in this study, we built a ventilated box based on modeling software. We analyzed the internal flow by simulation software. Simulation software is known for its advanced physics modeling capabilities, which include turbulence modeling, single- and multiphase flows, combustion, battery modeling, fluid-structure interaction, and much more. The sample selection method used in this paper is the orthogonal experimental design method, which is suitable for mechanical production...

Materiały

Name	Company	Catalog Number	Comments
Hardware
NVIDIA GPU	NVIDIA	N/A	An NVIDIA GPU is needed as some of the software frameworks below will not work otherwise. https://www.nvidia.com
Software
Ansys-Workbench	ANSYS	N/A	Multi-purpose finite element method computer design program software.https://www.ansys.com
SOLIDWORKS	Dassault Systemes	N/A	SolidWorks provides different design solutions, reduces errors in the design process, and improves product quality www.solidworks.com
SPSS	IBM	N/A	Software products for statistical analytical operations, data mining, predictive analysis, and decision support tasks software.https://www.ibm.com