<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.1 20151215//EN" "http://jats.nlm.nih.gov/publishing/1.1/JATS-journalpublishing1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="1.1">
<front>
<journal-meta>
<journal-id journal-id-type="pmc">FDMP</journal-id>
<journal-id journal-id-type="nlm-ta">FDMP</journal-id>
<journal-id journal-id-type="publisher-id">FDMP</journal-id>
<journal-title-group>
<journal-title>Fluid Dynamics &#x0026; Materials Processing</journal-title>
</journal-title-group>
<issn pub-type="epub">1555-2578</issn>
<issn pub-type="ppub">1555-256X</issn>
<publisher>
<publisher-name>Tech Science Press</publisher-name>
<publisher-loc>USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">16692</article-id>
<article-id pub-id-type="doi">10.32604/fdmp.2021.016692</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Article</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Experimental and Numerical Analysis of the Relationship between Pressure and Drip Rate in a Vertical Polypropylene Infusion Bag</article-title><alt-title alt-title-type="left-running-head">Experimental and Numerical Analysis of the Relationship between Pressure and Drip Rate in a Vertical Polypropylene Infusion Bag</alt-title><alt-title alt-title-type="right-running-head">Experimental and Numerical Analysis of the Relationship between Pressure and Drip Rate in a Vertical Polypropylene Infusion Bag</alt-title>
</title-group>
<contrib-group content-type="authors">
<contrib id="author-1" contrib-type="author" corresp="yes">
<name name-style="western">
<surname>Duan</surname>
<given-names>Weiwei</given-names>
</name>
<xref ref-type="aff" rid="aff-1"/>
<email>2190130134@stu.ahpu.edu.cn</email>
</contrib>
<contrib id="author-2" contrib-type="author">
<name name-style="western">
<surname>Tang</surname>
<given-names>Lingfeng</given-names>
</name>
<xref ref-type="aff" rid="aff-1"/>
</contrib>
<aff id="aff-1">
<institution>School of Mechanical Engineering, Anhui Polytechnic University</institution>, <addr-line>Wuhu, 241000</addr-line>, <country>China</country></aff>
</contrib-group><author-notes><corresp id="cor1">&#x002A;Corresponding Author: Weiwei Duan. Email: <email>2190130134@stu.ahpu.edu.cn</email></corresp></author-notes>
<pub-date pub-type="epub" date-type="pub" iso-8601-date="2021-09-04">
<day>4</day>
<month>9</month>
<year>2021</year>
</pub-date>
<volume>17</volume>
<issue>6</issue>
<fpage>1143</fpage>
<lpage>1164</lpage>
<history>
<date date-type="received">
<day>20</day>
<month>4</month>
<year>2021</year>
</date>
<date date-type="accepted">
<day>13</day>
<month>7</month>
<year>2021</year>
</date>
</history>
<permissions>
<copyright-statement>&#x00A9; 2021 Duan and Tang</copyright-statement>
<copyright-year>2021</copyright-year>
<copyright-holder>Duan and Tang</copyright-holder>
<license xlink:href="https://creativecommons.org/licenses/by/4.0/">
<license-p>This work is licensed under a <ext-link ext-link-type="uri" xlink:type="simple" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution 4.0 International License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p>
</license>
</permissions>
<self-uri content-type="pdf" xlink:href="TSP_FDMP_16692.pdf"></self-uri>
<abstract>
<p>Vertical infusion (self-emptying) bags used for Intravenous infusion are typically obtained by moulding a soft envelope of polypropylene. In normal conditions a continuous flow of liquid can be obtained with no need to use a pump. In the present study, the relationship between air pressure effects and the drug drip rate have been investigated experimentally and numerically. After determining relevant experimental data about the descending height of liquid level, the dropping speed and pressure, the ordinary least square method and MATLAB have been used to reconstruct the related variation and interrelation laws. Numerical simulations have been performed to determine the best gas-liquid volume ratio and improve the overall performances of these bags. According to these results, that the biggest effect on the drip rate is produced by the diameters of the used needles.</p>
</abstract>
<kwd-group kwd-group-type="author">
<kwd>Polypropylene infusion bag</kwd>
<kwd>numerical simulations</kwd>
<kwd>gas-liquid volume ratio</kwd>
<kwd>overall performances</kwd>
<kwd>drip rate</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1">
<label>1</label>
<title>Introduction</title>
<p>Intravenous infusion is commonly used in clinical medicine and infusion bags, whose standard capacity is 250, are widely used for various annual demands [<xref ref-type="bibr" rid="ref-1">1</xref>]. At present, most hospitals lack of that equipment for air purification and hygiene cleaning. Therefore, it is particularly paramount to make infusion bags self-emptying in barometric pressure in the closed infusion process, rather than air circuit. In addition, the right way to regulate the rate of infusion is effective treatment and to control the rate of infusion is an important link to ensure the quality of treatment and care [<xref ref-type="bibr" rid="ref-2">2</xref>]. The study on the pressure and rules of drip rate in the infusion process of vertical polypropylene infusion bags could be practical to guide new structures of infusion bags and to design more reasonable infusion bags [<xref ref-type="bibr" rid="ref-3">3</xref>&#x2013;<xref ref-type="bibr" rid="ref-6">6</xref>].</p>
<p>With the progress of medical science, more stringent requirements for airtightness, costs in economy, suitability of drug containers, vertical infusion bags have been on the stage, which attracts scholars to conduct numerous researches on vertical infusion bags. Liang et al. [<xref ref-type="bibr" rid="ref-7">7</xref>] have discussed the application of pressurized infusion bag in cardiac arrest. Using the cardioplegia perfusion to protect myocardium is safe, easy to operate the pressure of perfusion and the controllability of capacity are in good controllability. However, they did not simulate the pattern of pressure and simulation within infusion bags was not simulated. Song [<xref ref-type="bibr" rid="ref-8">8</xref>] has analyzed the rapid method to determinate the rate of intravenous drip, which improved the working efficiency. However, according to long-term clinical practices, it could not be used to accurately describe the drip rate of infusion process. Xiao et al. [<xref ref-type="bibr" rid="ref-9">9</xref>] and other researchers have proposed that real-time monitoring of the rate of infusion drip needs to monitor liquid droplets and the liquid level in most cases. In addition, there are few mistakes of measurements to monitor infusion drips, alarms are timely and reliable, due to various factors in the environment, when to call the nurse to replace the liquid still needs further consideration, and timely replacement of liquid still needs further consideration. Most of these aforementioned discussions are based on experience rather than scientific experiments, which demonstrate that this research is not scientific enough.</p>
<p>Scientific experiments are generally based on the theories while the results are often affected by the factors which should be found in comprehensive experiments. In addition, the differences of individuals in experiments often lead to the variety and instability of the results. Therefore, it is difficult to apply data in statistics for medical treatment. In this paper, a kind of simulator of human blood vessel was used to build devices of infusion and make experimental environment closer to real hospitals, then the pressure and drip rate in self-emptying were obtained by using the principle of fluid mechanics. The simulation of flow field was carried out by CFD to obtain the pressure and velocity field. The experimental data were fitted with MATLAB to require condition that the pressure and the drawdown rate of liquid level are in curve relationships and constitute a quartic function; lastly, the biggest factor that affects the infusion rate through experiments would be fingered out [<xref ref-type="bibr" rid="ref-10">10</xref>]. This article provides a theoretical basis for obtaining the design method of the infusion bag and obtaining the drip rate law of intravenous infusion.</p>
</sec>
<sec id="s2">
<label>2</label>
<title>Theoretical Calculations</title>
<p>As the actual liquid transport speed is comparatively slow (less than 90 drops), which conforms to the laminar flow law, it is possible to calculate the total pipe pressure loss by water head losing and the local pressure loss [<xref ref-type="bibr" rid="ref-11">11</xref>].</p>
<p>The pressure loss of the fluid in the pipe flow is defined as the <xref ref-type="disp-formula" rid="eqn-1">Eq. (1)</xref>:</p>
<p><disp-formula id="eqn-1">
<label>(1)</label>
<!--<alternatives>
<graphic mimetype="image" mime-subtype="png" xlink:href="eqn-1.png"/><tex-math id="tex-eqn-1"><![CDATA[{\rm \Delta p = }{\rm \lambda }\displaystyle{l \over d}{\rm }\displaystyle{{{\rm \rho }{{\rm \nu }^2}} \over {{\rm 2 }}}]]></tex-math>--><mml:math id="mml-eqn-1" display="block"><mml:mrow><mml:mi mathvariant="normal">&#x0394;</mml:mi><mml:mi mathvariant="normal">p</mml:mi><mml:mo>&#x003D;</mml:mo></mml:mrow><mml:mrow><mml:mi>&#x03BB;</mml:mi></mml:mrow><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mi>l</mml:mi><mml:mi>d</mml:mi></mml:mfrac></mml:mrow><mml:mrow></mml:mrow><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mrow><mml:mrow><mml:mi>&#x03C1;</mml:mi></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mi>&#x03BD;</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:mrow></mml:mfrac></mml:mrow></mml:mstyle></mml:mstyle></mml:math>
<!--</alternatives>--></disp-formula></p>
<p>where, l is pipe length, d is pipe diameter, &#x03C1; is fluid density, v is fluid speed, and &#x03BB; is the loss of resistance along the project.</p>
<p>When fluid flows through various local obstruction devices, due to different influential factors, the flow is extremely complex and not easy to be calculated theoretically; hence the local pressure loss is generally calculated by the following <xref ref-type="disp-formula" rid="eqn-2">Eq. (2)</xref>:</p>
<p><disp-formula id="eqn-2">
<label>(2)</label>
<!--<alternatives>
<graphic mimetype="image" mime-subtype="png" xlink:href="eqn-2.png"/><tex-math id="tex-eqn-2"><![CDATA[{\rm \Delta }{{\rm p}_{\rm \zeta }}{\rm = }{\rm \zeta }\displaystyle{{{\rm \nu }{{\rm }^{\rm 2}}} \over {{\rm 2 g}}}]]></tex-math>--><mml:math id="mml-eqn-2" display="block"><mml:mrow><mml:mi mathvariant="normal">&#x0394;</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x03B6;</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mo>&#x003D;</mml:mo></mml:mrow><mml:mrow><mml:mi>&#x03B6;</mml:mi></mml:mrow><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mrow><mml:mrow><mml:mi>&#x03BD;</mml:mi></mml:mrow><mml:mrow><mml:msup><mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mi mathvariant="normal">g</mml:mi></mml:mrow></mml:mrow></mml:mfrac></mml:mrow></mml:mstyle></mml:math>
<!--</alternatives>--></disp-formula></p>
<p>where, &#x03B6; is the coefficient of local resistance.</p>
<p>Therefore, the total pressure loss in the pipeline is the <xref ref-type="disp-formula" rid="eqn-3">Eq. (3)</xref>:</p>
<p><disp-formula id="eqn-3">
<label>(3)</label>
<!--<alternatives>
<graphic mimetype="image" mime-subtype="png" xlink:href="eqn-3.png"/><tex-math id="tex-eqn-3"><![CDATA[{\rm \Delta p = }\sum {\rm \lambda } {\rm }\displaystyle{l \over d}{\rm }\displaystyle{{{\rm \rho }{{\rm \nu }^2}} \over {{\rm 2 }}}{\rm + }\sum {\rm \zeta } {\rm }\displaystyle{{{{\rm \nu }^2}} \over {2g}}]]></tex-math>--><mml:math id="mml-eqn-3" display="block"><mml:mrow><mml:mi mathvariant="normal">&#x0394;</mml:mi><mml:mi mathvariant="normal">p</mml:mi><mml:mo>&#x003D;</mml:mo></mml:mrow><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>&#x03BB;</mml:mi></mml:mrow><mml:mrow></mml:mrow><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mi>l</mml:mi><mml:mi>d</mml:mi></mml:mfrac></mml:mrow><mml:mrow></mml:mrow><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mrow><mml:mrow><mml:mi>&#x03C1;</mml:mi></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mi>&#x03BD;</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:mrow></mml:mfrac></mml:mrow><mml:mrow><mml:mo>&#x002B;</mml:mo></mml:mrow><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>&#x03B6;</mml:mi></mml:mrow><mml:mrow></mml:mrow><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mi>&#x03BD;</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mi>g</mml:mi></mml:mrow></mml:mfrac></mml:mrow></mml:mstyle></mml:mstyle></mml:mstyle></mml:math>
<!--</alternatives>--></disp-formula></p>
<p>Infusion process diagram is shown in <xref ref-type="fig" rid="fig-1">Fig. 1</xref>, and the <xref ref-type="disp-formula" rid="eqn-4">Eq. (4)</xref> can be obtained by Bernoulli equation in the initial state.</p>
<p><disp-formula id="eqn-4">
<label>(4)</label>
<!--<alternatives>
<graphic mimetype="image" mime-subtype="png" xlink:href="eqn-4.png"/><tex-math id="tex-eqn-4"><![CDATA[\displaystyle{{{{\rm \nu }_0}^2} \over {{\rm 2g }}}{\rm + }{z_0}{\rm + }\displaystyle{{{p_0}} \over {{\rm \rho }g}}{\rm = (}\displaystyle{{{{\rm \nu }_1}^2} \over {2g}}{\rm + }{z_1}{\rm + }\displaystyle{{{p_1}} \over {{\rm \rho }g}}{\rm ) + \Delta }{{\rm p}_1}]]></tex-math>--><mml:math id="mml-eqn-4" display="block"><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x03BD;</mml:mi></mml:mrow><mml:mn>0</mml:mn></mml:msub></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mi mathvariant="normal">g</mml:mi></mml:mrow></mml:mrow></mml:mfrac></mml:mrow><mml:mrow><mml:mo>&#x002B;</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:mo>&#x002B;</mml:mo></mml:mrow><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mrow><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mi>&#x03C1;</mml:mi></mml:mrow><mml:mi>g</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mrow><mml:mo>&#x003D;</mml:mo><mml:mo stretchy="false">(</mml:mo></mml:mrow><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x03BD;</mml:mi></mml:mrow><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mi>g</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mrow><mml:mo>&#x002B;</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:mo>&#x002B;</mml:mo></mml:mrow><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mrow><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mi>&#x03C1;</mml:mi></mml:mrow><mml:mi>g</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo>&#x002B;</mml:mo><mml:mi mathvariant="normal">&#x0394;</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:mstyle></mml:mstyle></mml:mstyle></mml:mstyle></mml:math>
<!--</alternatives>--></disp-formula></p>
<fig id="fig-1">
<label>Figure 1</label>
<caption>
<title>Schematic diagram of infusion process, 1&#x2014;gas zone, 2&#x2014;murphy tube, 3&#x2014;water column of 24 mm height</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-1.png"/>
</fig>
<p>At the end of the state, <xref ref-type="disp-formula" rid="eqn-5">Eq. (5)</xref> can be obtained as follows:</p>
<p><disp-formula id="eqn-5">
<label>(5)</label>
<!--<alternatives>
<graphic mimetype="image" mime-subtype="png" xlink:href="eqn-5.png"/><tex-math id="tex-eqn-5"><![CDATA[\displaystyle{{{{\rm \nu }_2}^2} \over {{\rm 2g }}}{\rm + }{z_2}{\rm + }\displaystyle{{{p_2}} \over {{\rm \rho }g}}{\rm = (}\displaystyle{{{{\rm \nu }_1}^2} \over {2g}}{\rm + }{z_1}{\rm + }\displaystyle{{{p_1}} \over {{\rm \rho }g}}{\rm ) + \Delta }{{\rm p}_2}]]></tex-math>--><mml:math id="mml-eqn-5" display="block"><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x03BD;</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mi mathvariant="normal">g</mml:mi></mml:mrow></mml:mrow></mml:mfrac></mml:mrow><mml:mrow><mml:mo>&#x002B;</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:mo>&#x002B;</mml:mo></mml:mrow><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mrow><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mi>&#x03C1;</mml:mi></mml:mrow><mml:mi>g</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mrow><mml:mo>&#x003D;</mml:mo><mml:mo stretchy="false">(</mml:mo></mml:mrow><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x03BD;</mml:mi></mml:mrow><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mi>g</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mrow><mml:mo>&#x002B;</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mi>z</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:mo>&#x002B;</mml:mo></mml:mrow><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mrow><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mi>&#x03C1;</mml:mi></mml:mrow><mml:mi>g</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo>&#x002B;</mml:mo><mml:mi mathvariant="normal">&#x0394;</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi mathvariant="normal">p</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:mstyle></mml:mstyle></mml:mstyle></mml:mstyle></mml:math>
<!--</alternatives>--></disp-formula></p>
<p>At initial condition, v<sub>0</sub> is the speed of the liquid in the bag. p<sub>0</sub> is the air pressure, v<sub>1</sub> is the speed of the needle in the initial state. v<sub>2</sub> and p<sub>2</sub> is the speed and pressure of the liquid at a certain state. p<sub>1</sub> and z<sub>2</sub> are the pressure and height of the venous blood vessel simulator. z<sub>0</sub> and z<sub>2</sub> are the height from the liquid in bag to the arm of the venous blood vessel simulation device height in the initial state and the state of the infusion bag in any state, &#x0394;p<sub>1</sub> and &#x0394;p<sub>2</sub> are the total pipeline loss in the infusion process.</p>
<p>In Murphy tube, one hundred drops of liquid is about 6.5 ml, which can be converted to the velocity of 15 drops/ml. The output of the local loss of the needle accounts for a large proportion of water head losing, and the proportion of the input needle takes up a small proportion in all the energy loss, hence the latter is negligible. As the local pressure loss becoming larger with the higher speed, and reaching the maximum speed at the initial state, the corresponding local pressure loss in the initial state is selected.</p>
<p>When the needle external diameter is 0.5 mm, the internal diameter is 0.3 mm, and the initial speed is 57 drops/min, the flow rate could be calculated as 0.83 m/s. And when the final drip rate is 40 drops/min, the flow rate could be of 0.58 m/s. Confirming the table to get the local energy loss coefficient &#x03B6; of the needle with external diameter of 0.5 mm, it can be measured that the value of z<sub>0</sub>, z<sub>1</sub>, and z<sub>2</sub> of the 250 ml infusion bag are 1.14 m, 24 cm and 1.025 m. Substituting z<sub>0</sub>, z<sub>1</sub>, v<sub>0</sub>, v<sub>1</sub>, p<sub>1</sub>, p<sub>2</sub> into <xref ref-type="disp-formula" rid="eqn-1">Eqs. (1)</xref>&#x2013;<xref ref-type="disp-formula" rid="eqn-5">(5)</xref> to acquire the minimum value of P<sub>2</sub> is &#x2013;2640 Pa when dripping at 40 drops/min.</p>
</sec>
<sec id="s3">
<label>3</label>
<title>Flow Field Simulation of Infusion Process</title>
<p>In order to find the appropriate modeling method and to determine the reasonable boundary conditions and flow patterns, the infusion bag with the capacity of 250 ml was taken as an example, the FLUENT [<xref ref-type="bibr" rid="ref-12">12</xref>,<xref ref-type="bibr" rid="ref-13">13</xref>] software was used to simulate the unsteady flow of the infusion process in the initial state and near the evacuating state [<xref ref-type="bibr" rid="ref-14">14</xref>&#x2013;<xref ref-type="bibr" rid="ref-16">16</xref>].</p>
<sec id="s3_1">
<label>3.1</label>
<title>The Implementation of 3D Models</title>
<p>In this chapter, a 250 ml infusion bag was taken as an example and the simulation process of which is similar to other polypropylene infusion bags with different volumes. The theoretical data of three vertical polypropylene infusion bags and infusion tubes were obtained from measurements of Vernier calipers and rulers. Combined with the height of the infusion tube from the venous vessels described in Chapter 2, UG software was adopted to establish 3D models of infusion bottles and tubes. Finally, the model was implemented as follows: the liquid volume of the infusion bag was 250 ml; the initial liquid height in the bottle (the distance between the liquid surface and the bottle mouth) was 40.5 mm; the inner diameter of the output needle was 0.3 mm. The model is shown in <xref ref-type="fig" rid="fig-2">Fig. 2</xref>.</p>
<fig id="fig-2">
<label>Figure 2</label>
<caption>
<title>Infusion bag model</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-2.png"/>
</fig>
</sec>
<sec id="s3_2">
<label>3.2</label>
<title>Meshing</title>
<p>In general, the calculation of denser grids is extremely accurate and time-consuming. The tubes were meshed by a high-density grid while the polypropylene infusion bag was meshed by a low-density grid. This method not only guaranteed the calculation accuracy but also took into account the calculation economy. The grid was divided by hexahedral structured mesh, and tetrahedral unstructured mesh according to the structural characteristics of each part. Before the grid was specifically divided, the model needed to be processed accordingly.</p>
<list list-type="order">
<list-item><p>Firstly, all the independent parts such as the bottle, the bottom of the bottle and the tube were combined as a whole by Boolean operation, and then the overall model was divided according to the meshing needs. The operation of the merged model is shown in <xref ref-type="fig" rid="fig-3">Fig. 3</xref>.</p></list-item>
<list-item><p>The infusion bottle model was divided into several parts according to the structure and size. A surface was created to divide the body. <xref ref-type="fig" rid="fig-3">Fig. 3</xref> shows the segmentation operation of the model. There was a big difference between the size of the tube and the bottom of the bottle. Thus, the mesh needed to be divided separately. The rest of the body was also divided in the same way.</p></list-item>
<list-item><p>For parts such as the bottle and the tube, the relevant surface should be meshed before dividing the body grids. <xref ref-type="fig" rid="fig-4">Fig. 4</xref> demonstrates the meshing operation of the tube section by surface meshing. <xref ref-type="fig" rid="fig-5">Fig. 5</xref> illustrates the grid of the middle section of the infusion bottle, which adopts quadrilateral structure. The surface meshing size of each part is determined according to the structure size.</p></list-item>
<list-item><p>The body meshing was performed on the basis of surface meshing. Considering the structural characteristics of infusion bottles and tubes, the cooper method was adopted. <xref ref-type="fig" rid="fig-6">Fig. 6</xref> is the operation of structural body meshing. The unstructured mesh was used for the lower part of the bottle and the common surfaces of the adjacent body share a set of grids, which was convenient for later analysis and calculation. The body meshing of infusion bottle is shown in <xref ref-type="fig" rid="fig-7">Fig. 7</xref> and the total number of grids is about 1.7 million.</p></list-item>
</list>
<fig id="fig-3">
<label>Figure 3</label>
<caption>
<title>Segmentation operation of model</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-3.png"/>
</fig>
<fig id="fig-4">
<label>Figure 4</label>
<caption>
<title>Surface meshing operation</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-4.png"/>
</fig>
<fig id="fig-5">
<label>Figure 5</label>
<caption>
<title>Middle section mesh of infusion bottle</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-5.png"/>
</fig>
<fig id="fig-6">
<label>Figure 6</label>
<caption>
<title>Body meshing operation</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-6.png"/>
</fig>
<fig id="fig-7">
<label>Figure 7</label>
<caption>
<title>Body meshing of infusion bottle</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-7.png"/>
</fig>
<p>Two speed monitoring points were set at the exit of the infusion bag and the middle of the tube respectively, as shown in <xref ref-type="fig" rid="fig-8">Fig. 8</xref>. Grid-independent verification was shown in <xref ref-type="table" rid="table-1">Tab. 1</xref>.</p>
<fig id="fig-8">
<label>Figure 8</label>
<caption>
<title>Speed monitoring points</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-8.png"/>
</fig>

<table-wrap id="table-1">
<label>Table 1</label>
<caption>
<title>Grid-independent verification</title>
</caption>
<table>
<colgroup>
<col/>
<col/>
<col/>
<col/>
</colgroup>
<thead>
<tr>
<th>Number</th>
<th>Grids</th>
<th>Velocity of point 1 (m/s)</th>
<th>Velocity of point 2 (m/s)</th>
</tr>
</thead>
<tbody>
<tr>
<td>1</td>
<td>1034299</td>
<td>0.0321</td>
<td>0.0524</td>
</tr>
<tr>
<td>2</td>
<td>1193090</td>
<td>0.0334</td>
<td>0.0537</td>
</tr>
<tr>
<td>3</td>
<td>1323427</td>
<td>0.0353</td>
<td>0.0546</td>
</tr>
<tr>
<td>4</td>
<td>1783267</td>
<td>0.0372</td>
<td>0.0550</td>
</tr>
<tr>
<td>5</td>
<td>2263253</td>
<td>0.0374</td>
<td>0.0551</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3_3">
<label>3.3</label>
<title>Flow Field Calculation</title>
<p>For the flow in the infusion bottle, the liquid in the bottle was continuously reduced and the liquid level was continuously decreasing as time goes by, which was an unsteady flow. Besides, the infusion bottle was constantly deformed and the gas phase pressure in the bottle was constantly changing with the liquid level. The flow field analysis cannot analyze the deformation of the bottle and the pressure change in the bottle due to the deformation of the bottle body could not be obtained. Therefore, the deformation of the bottle was not taken into consideration in this calculation. A fixed position was taken for steady-state analysis in order to analyze the flow rate at the outlet under the conditions of a given bottle pressure and the height of the liquid level in the bottle.</p>
</sec>
<sec id="s3_4">
<label>3.4</label>
<title>Setting Boundary Conditions</title>
<p>(1) Calculating Model and Boundary Conditions</p>
<p>Calculation time was Steady. Viscous Model was Laminar. The VOF [<xref ref-type="bibr" rid="ref-17">17</xref>,<xref ref-type="bibr" rid="ref-18">18</xref>] gas-liquid two-phase flow model was adopted with the gas phase as the main term and the liquid phase as the second term. The calculation type was Steady state; the viscous flow model used a laminar flow equation. The outlet was a pressure outlet with a pressure of 2,500 Pa, shown in <xref ref-type="fig" rid="fig-9">Fig. 9</xref>. All the remaining wall conditions were set as WALL. The VOF model was set as shown in <xref ref-type="fig" rid="fig-10">Fig. 10</xref>. Solution Methods were PISO, and Under-Relaxation factors were all set to 0.2. Other parameters were default. Because it is a steady state calculation, the under-relaxation factors were adjusted to 0.2 in order to improve the stability of the calculation.</p>
<fig id="fig-9">
<label>Figure 9</label>
<caption>
<title>Outlet boundary condition setting</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-9.png"/>
</fig>
<fig id="fig-10">
<label>Figure 10</label>
<caption>
<title>VOF model setting</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-10.png"/>
</fig>
<p>The VOF method determine the free face by studying the volume ratio function of fluid and mesh in the grid cell, and track the change of fluid rather than the motion of particles on the free face. The definition function for the existence of the phase flow point on the grid is as follows (<xref ref-type="disp-formula" rid="eqn-6">Eq. (6)</xref>):</p>
<p><disp-formula id="eqn-6">
<label>(6)</label>
<!--<alternatives>
<graphic mimetype="image" mime-subtype="png" xlink:href="eqn-6.png"/><tex-math id="tex-eqn-6"><![CDATA[f\left( {x,y,t} \right) = \left\{ {\matrix{ 1 \cr 0 \cr } } \right.]]></tex-math>--><mml:math id="mml-eqn-6" display="block"><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>x</mml:mi><mml:mo>,</mml:mo><mml:mi>y</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x003D;</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mtable columnspacing="1em" rowspacing="4pt"><mml:mtr><mml:mtd><mml:mn>1</mml:mn></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn>0</mml:mn></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo stretchy="true" symmetric="true" fence="true"></mml:mo></mml:mrow></mml:math>
<!--</alternatives>--></disp-formula></p>
<p><disp-formula id="eqn-7">
<label>(7)</label>
<!--<alternatives>
<graphic mimetype="image" mime-subtype="png" xlink:href="eqn-7.png"/><tex-math id="tex-eqn-7"><![CDATA[\displaystyle{{\partial f} \over {\partial t}} + \displaystyle{{\partial uf} \over {\partial x}} + \displaystyle{{\partial vf} \over {\partial y}} = 0]]></tex-math>--><mml:math id="mml-eqn-7" display="block"><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mo>&#x002B;</mml:mo><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>u</mml:mi><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>x</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mo>&#x002B;</mml:mo><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>v</mml:mi><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mi>y</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mstyle></mml:mstyle></mml:mstyle></mml:math>
<!--</alternatives>--></disp-formula></p>
<p><disp-formula id="eqn-8">
<label>(8)</label>
<!--<alternatives>
<graphic mimetype="image" mime-subtype="png" xlink:href="eqn-8.png"/><tex-math id="tex-eqn-8"><![CDATA[\displaystyle{{F_{i,\;j}^{n + 1} - F_{i,\;j}^n} \over {\Delta t}} + \displaystyle{{\delta {F_{i + \displaystyle{1 \over 2},j}} - \delta {F_{i - \displaystyle{1 \over 2},j}}} \over {\Delta {x_i}}} + \displaystyle{{\delta {F_{i,\;j + \displaystyle{1 \over 2}}} - \delta {F_{i,\;j - \displaystyle{1 \over 2}}}} \over {\Delta {y_i}}} = 0]]></tex-math>--><mml:math id="mml-eqn-8" display="block"><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mrow><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>&#x002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x2212;</mml:mo><mml:msubsup><mml:mi>F</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi></mml:mrow><mml:mi>n</mml:mi></mml:msubsup></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x0394;</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mo>&#x002B;</mml:mo><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mrow><mml:mi>&#x03B4;</mml:mi><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x002B;</mml:mo><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac></mml:mrow><mml:mo>,</mml:mo><mml:mi>j</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mi>&#x03B4;</mml:mi><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac></mml:mrow><mml:mo>,</mml:mo><mml:mi>j</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x0394;</mml:mi><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mfrac></mml:mrow><mml:mo>&#x002B;</mml:mo><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mrow><mml:mi>&#x03B4;</mml:mi><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi><mml:mo>&#x002B;</mml:mo><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac></mml:mrow></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mi>&#x03B4;</mml:mi><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mstyle scriptlevel="0" displaystyle="true"><mml:mrow><mml:mfrac><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac></mml:mrow></mml:mstyle></mml:mrow></mml:msub></mml:mrow></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x0394;</mml:mi><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mfrac></mml:mrow><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mstyle></mml:mstyle></mml:mstyle></mml:math>
<!--</alternatives>--></disp-formula></p>
<p>The mass conservation formulation is shown in <xref ref-type="disp-formula" rid="eqn-7">Eq. (7)</xref>, the difference equation is shown in <xref ref-type="disp-formula" rid="eqn-8">Eq. (8)</xref>, and the difference grid is shown in <xref ref-type="fig" rid="fig-11">Fig. 11</xref>.</p>
<fig id="fig-11">
<label>Figure 11</label>
<caption>
<title>Variables on grid cell</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-11.png"/>
</fig>
<p>(2) Initialization</p>
<p>After the relevant settings were completed, and the flow field initialized, the gas phase area was defined first, and then calculation was performed. Then the gas phase area was under &#x201C;patch&#x201D; process. The gas percentage in gas phase area (The gas phase was located in the upper part of the infusion bottle.) was set as 1. <xref ref-type="fig" rid="fig-12">Fig. 12</xref> is the &#x201C;patch&#x201D; operation. When calculating, the residual graph was opened and the exit velocity value was monitored to judge whether the calculation was over.</p>
<fig id="fig-12">
<label>Figure 12</label>
<caption>
<title>Initial flow field patch settings</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-12.png"/>
</fig>
<p>In the initial stage of infusion, the model was implemented as follows: the volume of the infusion bottle was 250 ml; the height of the upper gas was 45 mm; the diameter of needle was 0.3 mm. When the air was approached, the initial height of the liquid in the bottle was 40.5 mm, and the needle&#x2019;s diameter was 0.3 mm. Thus, a two-part fluid model and meshing were implemented. The grid division with the total number of 567,778 is shown in <xref ref-type="fig" rid="fig-13">Fig. 13</xref>.</p>
<fig id="fig-13">
<label>Figure 13</label>
<caption>
<title>Meshing</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-13.png"/>
</fig>
<p>At the end of the infusion, the model was used as the calculation area. The pressure on the liquid surface was set at &#x2013;2,980 Pa, the pressure at the outlet of the needle was set at 2,400 Pa, and the flow medium was simplified as water. The infusion process simulated the flow condition within 3 s from this state. In order to compare with the experimental data and to verify the reliability of the simulation, the speed point P was set at the center of needle outlet section in the calculation process.</p>
</sec>
<sec id="s3_5">
<label>3.5</label>
<title>Simulation Analysis</title>
<p>The contour of velocity of the infusion bag outlet at the beginning of infusion is showed in <xref ref-type="fig" rid="fig-14">Fig. 14</xref>. It can be seen that the speed of the fluid is maximum at the center and the speed of fluid near the tube was 0, which is consistent with the theoretical analysis. Through the &#x201C;report&#x201D; function of FLUENT the average value of speed taken out of the outlet is converted into the dripping rate. The dripping rate at the beginning of infusion is about 57 drops. When the pressure in the bottle is &#x2013;2,000 Pa, the dripping rate is about 44 drops at the end.</p>
<fig id="fig-14">
<label>Figure 14</label>
<caption>
<title>Contour of velocity of the infusion bag outlet</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-14.png"/>
</fig>
<p><xref ref-type="fig" rid="fig-15">Fig. 15</xref> shows contour of velocity of the infusion bag inlet. It can be seen that the speed of the fluid is maximum at needle axis, which is approximately 0.5 m/s, while the speed declines gradually to zero near the tube wall.</p>
<fig id="fig-15">
<label>Figure 15</label>
<caption>
<title>Contour of velocity of the infusion bag inlet</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-15.png"/>
</fig>
<p><xref ref-type="fig" rid="fig-16">Fig. 16</xref> shows the pressure distribution in the bag. It can be seen that the pressure increases with the increase of the height. This is because the height of tube is 1 m, and the pressure difference is the manifestation of static pressure. The outlet boundary condition is the pressure outlet and the pressure value is 2,500 Pa and the pressure difference between it and the negative pressure is the static pressure of 1 m, which is reasonable.</p>
<fig id="fig-16">
<label>Figure 16</label>
<caption>
<title>Contour of static pressure in the bottle</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-16.png"/>
</fig>
<p>In the calculation process, a pressure measurement point P1 and a velocity measurement point P2 are respectively set, where P1 is located above the liquid inside the bottle (this point is always air), and P2 is the center point of the needle outlet section.</p>
<p><xref ref-type="fig" rid="fig-17">Fig. 17</xref> shows pressure change at a distance of 20 mm from the top of the bag. As time goes by, the pressure in the bottle gradually decreases. As can be seen from <xref ref-type="fig" rid="fig-18">Fig. 18</xref>, the pressure at P1 above the bottle was gradually decreasing linearly with time in the first 3 s. In about 1 s, the negative pressure reaches &#x2013;2,993 Pa. This is caused by the increase in the air volume in the bottle as the decrease of the liquid level, which is consistent with the actual physical change.</p>
<fig id="fig-17">
<label>Figure 17</label>
<caption>
<title>Pressure change with time</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-17.png"/>
</fig>
<fig id="fig-18">
<label>Figure 18</label>
<caption>
<title>The pressure at P1 point with time</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-18.png"/>
</fig>
<p>The outlet velocity curve at a distance of 20 mm from the top of the bag with time can be seen from <xref ref-type="fig" rid="fig-19">Fig. 19</xref>, and <xref ref-type="fig" rid="fig-19">Fig. 19(b)</xref> is the partial enlarged view of <xref ref-type="fig" rid="fig-19">Fig. 19(a)</xref>. It indicates that the speed fluctuates slightly at first and then gradually stabilizes.</p>
<fig id="fig-19">
<label>Figure 19</label>
<caption>
<title>Outlet peed curve with time curve</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-19.png"/>
</fig>
<p>It can be seen from <xref ref-type="fig" rid="fig-20">Fig. 20(a)</xref> that except for the unstable in the first 0.5 s, the speed gradually tends to be stable and decreases with time. The speed value is 0.499 m/s at 1 s, 0.498 m/s at 2 s and 0.496 m/s at 3 s. <xref ref-type="fig" rid="fig-20">Fig. 20(b)</xref> is a partial enlarged view of <xref ref-type="fig" rid="fig-20">Fig. 20(a)</xref>, from which the variation of the speed is more clearly.</p>
<fig id="fig-20">
<label>Figure 20</label>
<caption>
<title>(a) Velocity of point P2 changes with time, (b) Speed change curve within 2&#x2013;3 s</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-20.png"/>
</fig>
</sec>
</sec>
<sec id="s4">
<label>4</label>
<title>Experimental Verification</title>
<p>The previous part simulated the internal flow field of infusion bag during human intravenous infusion. However, some boundary conditions may be set inaccurately due to the deformation of infusion bag, leaving the simulation results be different from the actual data. In this paper, a simulated experimental device for human intravenous injection was implemented. Through experiments, a lot of data such as the pressure and droplet rate of infusion bag liquid during human intravenous injection were obtained. According to numerical analysis of the data, an accurate variation rule of the pressure and droplet rate in human intravenous infusion bag was implemented.</p>
<sec id="s4_1">
<label>4.1</label>
<title>Human Intravenous Simulator</title>
<p>Human beings&#x2019; central venous pressure is normally 5 to 12 cm of water column, which is less than the peripheral venous pressure. In order to be safe and reliable, the static pressure of 24 cm was designed for the intravenous simulation device pressure which is twice the maximum normal value. If the total mechanical energy of the needle at the needle exit was always larger than the static pressure, the liquid would pass through the device smoothly, so that the infusion can be completed normally. Therefore, it is reasonable to judge whether there is a blood return by confirming if the liquid could flow out of the simulation device.</p>
</sec>
<sec id="s4_2">
<label>4.2</label>
<title>Infusion Device</title>
<p>The entire infusion process test equipment included: a thickened clawed bag stopper with adjustable frame, multiple of disposable infusion sets, multiple medical infusion tubes (needle size is 0.5 mm), an infusion tube with small needle (for U-tube to measure pressure), stopwatch, paper, pen, tape.</p>
<p>The infusion bag test bench device was mainly composed of the infusion of soft bags, a Murphy tube, connecting hose and designed fixed frame. Infusion soft bag was fixed on the soft bag holder; one end of the Philippine tube was fixed on the frame with one end connected to the bottom of the tube and the other hanging on the liquid recovery box. One side of the hose was connected to the bottom of the infusion soft bag and the other side to the fixed frame. Both the main mark and the secondary bar were equipped with the scale marks; the main benchmark bar was installed on the main scale while the sub-bar was installed on the sub-standard in order to read the infusion of the soft bag level and the height of the hose level easily. The main scale and the sub-scale that could slide up and down to facilitate the liquid level descent height calibration were equipped with a guide sleeve fitted with fastening nut and a pointer, which will help to read the data.</p>
<p>In airtight infusion, with the liquid flowing out of the bag, the volume inside the bag was getting smaller and the bag started deforming gradually and the drip rate became slowly, there is a certain relationship between the output of liquid volume, pressure of infusion bag and the drop velocity. Because it was difficult to measure the volume of the output liquid and there was a part of the liquid remained in the tube, the mathematical model with the volume as an independent variable and drip as the dependent variable was inappropriate. Due to the height of the liquid drop can be observed without interrupting the experimental observation, the relationship between the height drop and the titer can be better observed by viewing the drop level of the liquid as an independent variable [<xref ref-type="bibr" rid="ref-19">19</xref>].</p>
<p>In the infusion process, the pressure inside the bag is real-time changing with the continuous flow of liquid. The conventional infusion device usually directly hangs the infusion bag at a certain height and the infusion was performed by inserting an infusion tube into the rubber stopper. Because there was no way to penetrate the internal measurement of pressure, the pressure changes cannot be measured by the infusion of soft bottles in the infusion process. In order to obtain real-time infusion process infusion bottle pressure and the pressure of the changing pattern, put forward a pressure measurement device, as shown in <xref ref-type="fig" rid="fig-21">Fig. 21</xref>.</p>
<fig id="fig-21">
<label>Figure 21</label>
<caption>
<title>Infusion test bench device diagram. 1&#x2014;Polypropylene infusion bag, 2&#x2014;Murphy&#x2019;s tube, 3&#x2014;Connected hose, 4&#x2014;Block, 5&#x2014;Vinyl holder, 6&#x2014;Fluoride retainer, 7&#x2014;Murphy&#x2019;s holder, 8&#x2014;Hose holder, 9&#x2014;Main benchmark, 10&#x2014;Sub-bar, 11&#x2014;Main scale, 12&#x2014;Sub-scale, 13&#x2014;Sleeve, 14&#x2014;Pointer, 15&#x2014;Fastening nut, 16&#x2014;Liquid recovery box</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-21.png"/>
</fig>
<p>Dripping rate refers to the droplets in per unit time of the Murphy tube. In the initial state, the liquid level in the Murphy tube should not be higher than half of the height of the Murphy, which facilitated the statistical drip rate. Using cylinder to measure 100 drops, and calculating each drop&#x2019;s volume, then converted into meters per second. Through the scale of the ruler that connected to infusion bag, the liquid declining level can be measured. The drip rate measurement was shown in <xref ref-type="fig" rid="fig-22">Fig. 22</xref>. The principle of measuring the liquid height was shown in <xref ref-type="fig" rid="fig-23">Fig. 23</xref> and the measurement of the 0.55 mm needle was shown in <xref ref-type="table" rid="table-2">Tab. 2</xref>.</p>
<fig id="fig-22">
<label>Figure 22</label>
<caption>
<title>Schematic diagram of the pressure test device. 1&#x2014;Polypropylene infusion bag; 2&#x2014;infusion tube; 3&#x2014;Murphy tube; 4&#x2014;U-tube, A&#x2014;liquid level</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-22.png"/>
</fig>
<fig id="fig-23">
<label>Figure 23</label>
<caption>
<title>Schematic diagram of drip rate measurement</title>
</caption>
<graphic mimetype="image" mime-subtype="png" xlink:href="FDMP_16692-fig-23.png"/>
</fig>

<table-wrap id="table-2">
<label>Table 2</label>
<caption>
<title>Pressure and drip rate data in the infusion process of 250 ml infusion bag (0.5 mm needle)</title>
</caption>
<table>
<colgroup>
<col/>
<col/>
<col/>
<col/>
<col/>
<col/>
<col/>
<col/>
<col/>
<col/>
<col/>
</colgroup>
<tbody>
<tr>
<td></td>
<td colspan="2">Group<break/>one</td>
<td colspan="2">Group<break/>two</td>
<td colspan="2">Group<break/>three</td>
<td colspan="2">Group<break/>four</td>
<td colspan="2">Average<break/>value</td>
</tr><tr>
<td>Height<break/>(cm)</td>
<td>Pressure<break/>(Pa)</td>
<td>Drip<break/>(drop)</td>
<td>Pressure<break/>(Pa)</td>
<td>Drip<break/>(drop)</td>
<td>Pressure<break/>(Pa)</td>
<td>Drip<break/>(drop)</td>
<td>Pressure<break/>(Pa)</td>
<td>Drip<break/>(drop)</td>
<td>Pressure<break/>(Pa)</td>
<td>Drip<break/>(drop)</td>
</tr>
<tr>
<td>0</td>
<td>&#x2013;450</td>
<td>57</td>
<td>&#x2013;290</td>
<td>57</td>
<td>&#x2013;540</td>
<td>63</td>
<td>&#x2013;490</td>
<td>52</td>
<td>&#x2013;442</td>
<td>57</td>
</tr>
<tr>
<td>0.5</td>
<td>&#x2013;900</td>
<td>53</td>
<td>&#x2013;800</td>
<td>57</td>
<td>&#x2013;1030</td>
<td>60</td>
<td>&#x2013;950</td>
<td>49</td>
<td>&#x2013;920</td>
<td>55</td>
</tr>
<tr>
<td>1</td>
<td>&#x2013;890</td>
<td>53</td>
<td>&#x2013;930</td>
<td>55</td>
<td>&#x2013;1060</td>
<td>59</td>
<td>&#x2013;1010</td>
<td>51</td>
<td>&#x2013;920</td>
<td>55</td>
</tr>
<tr>
<td>1.5</td>
<td>&#x2013;890</td>
<td>53</td>
<td>&#x2013;790</td>
<td>57</td>
<td>&#x2013;1100</td>
<td>57</td>
<td>&#x2013;980</td>
<td>51</td>
<td>&#x2013;940</td>
<td>55</td>
</tr>
<tr>
<td>2</td>
<td>&#x2013;840</td>
<td>53</td>
<td>&#x2013;820</td>
<td>56</td>
<td>&#x2013;890</td>
<td>59</td>
<td>&#x2013;1010</td>
<td>51</td>
<td>&#x2013;890</td>
<td>55</td>
</tr>
<tr>
<td>2.5</td>
<td>&#x2013;820</td>
<td>53</td>
<td>&#x2013;890</td>
<td>55</td>
<td>&#x2013;1110</td>
<td>56</td>
<td>&#x2013;1040</td>
<td>50</td>
<td>&#x2013;965</td>
<td>54</td>
</tr>
<tr>
<td>3</td>
<td>&#x2013;800</td>
<td>53</td>
<td>&#x2013;820</td>
<td>55</td>
<td>&#x2013;1080</td>
<td>56</td>
<td>&#x2013;1030</td>
<td>50</td>
<td>&#x2013;932</td>
<td>54</td>
</tr>
<tr>
<td>3.5</td>
<td>&#x2013;850</td>
<td>52</td>
<td>&#x2013;860</td>
<td>55</td>
<td>&#x2013;1060</td>
<td>56</td>
<td>&#x2013;960</td>
<td>50</td>
<td>&#x2013;932</td>
<td>54</td>
</tr>
<tr>
<td>4</td>
<td>&#x2013;900</td>
<td>52</td>
<td>&#x2013;880</td>
<td>54</td>
<td>&#x2013;1110</td>
<td>55</td>
<td>&#x2013;973</td>
<td>49</td>
<td>&#x2013;900</td>
<td>53</td>
</tr>
<tr>
<td>4.5</td>
<td>&#x2013;1010</td>
<td>51</td>
<td>&#x2013;700</td>
<td>54</td>
<td>&#x2013;1130</td>
<td>55</td>
<td>&#x2013;1040</td>
<td>49</td>
<td>&#x2013;970</td>
<td>50</td>
</tr>
<tr>
<td>5</td>
<td>&#x2013;1100</td>
<td>50</td>
<td>&#x2013;780</td>
<td>53</td>
<td>&#x2013;1210</td>
<td>54</td>
<td>&#x2013;1030</td>
<td>48</td>
<td>&#x2013;1030</td>
<td>51</td>
</tr>
<tr>
<td>5.5</td>
<td>&#x2013;1200</td>
<td>50</td>
<td>&#x2013;830</td>
<td>53</td>
<td>&#x2013;1290</td>
<td>53</td>
<td>&#x2013;1020</td>
<td>48</td>
<td>&#x2013;1085</td>
<td>51</td>
</tr>
<tr>
<td>6</td>
<td>&#x2013;1400</td>
<td>48</td>
<td>&#x2013;890</td>
<td>53</td>
<td>&#x2013;1320</td>
<td>53</td>
<td>&#x2013;1100</td>
<td>47</td>
<td>&#x2013;1178</td>
<td>50</td>
</tr>
<tr>
<td>6.5</td>
<td>&#x2013;1510</td>
<td>46</td>
<td>&#x2013;1040</td>
<td>53</td>
<td>&#x2013;1330</td>
<td>52</td>
<td>&#x2013;1170</td>
<td>47</td>
<td>&#x2013;1263</td>
<td>49</td>
</tr>
<tr>
<td>7</td>
<td>&#x2013;1400</td>
<td>46</td>
<td>&#x2013;1250</td>
<td>51</td>
<td>&#x2013;1420</td>
<td>52</td>
<td>&#x2013;1380</td>
<td>45</td>
<td>&#x2013;1363</td>
<td>49</td>
</tr>
<tr>
<td>7.5</td>
<td>&#x2013;1570</td>
<td>45</td>
<td>&#x2013;1490</td>
<td>50</td>
<td>&#x2013;1470</td>
<td>50</td>
<td>&#x2013;1500</td>
<td>43</td>
<td>&#x2013;1507</td>
<td>47</td>
</tr>
<tr>
<td>8</td>
<td>&#x2013;1930</td>
<td>42</td>
<td>&#x2013;1710</td>
<td>47</td>
<td>&#x2013;1950</td>
<td>49</td>
<td>&#x2013;1800</td>
<td>42</td>
<td>&#x2013;1847</td>
<td>45</td>
</tr>
<tr>
<td>8.5</td>
<td>&#x2013;2130</td>
<td>42</td>
<td>&#x2013;1680</td>
<td>47</td>
<td>&#x2013;2220</td>
<td>45</td>
<td>&#x2013;2030</td>
<td>40</td>
<td>&#x2013;2015</td>
<td>44</td>
</tr>
<tr>
<td>9</td>
<td>&#x2013;2640</td>
<td>40</td>
<td>&#x2013;1920</td>
<td>45</td>
<td>&#x2013;2700</td>
<td>43</td>
<td>&#x2013;1970</td>
<td>38</td>
<td>&#x2013;2095</td>
<td>42</td>
</tr>
<tr>
<td>9.5</td>
<td>&#x2013;2980</td>
<td>35</td>
<td>&#x2013;2150</td>
<td>43</td>
<td>&#x2013;3020</td>
<td>40</td>
<td>&#x2013;1920</td>
<td>37</td>
<td>&#x2013;2220</td>
<td>39</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s4_3">
<label>4.3</label>
<title>The Implementation of Mathematical Model</title>
<p>After the implementation of an experimental platform, according to the above description, the device was used to measure the droplet velocity. An infusion experience was made with the specification of 250 ml to detect the variation of droplet velocity and the liquid drop height in unit time. The relationship between the droplet velocity and liquid drop height in the tube was observed. With 8&#x2013;10 groups of experiment, the most interfering factors were removed, and the remaining four groups were averaged. An experimental model of pressure in different volumes of bags between the droplet velocity and the drop height was conducted by means of numeric analysis to conclude the law of drip velocity. The data of droplet velocity can be seen in <xref ref-type="table" rid="table-2">Tab. 2</xref>. The surface drop height value h, pressure P and the drip value V measured in above experiment were input into the MATLAB [<xref ref-type="bibr" rid="ref-20">20</xref>] software for analyzing the obtaining of the change curve and the function formula of drop rate as well as the height when the liquid dropped to h; the coefficients of the fitted curve was gained with least square fitting of a polynomial. Finally, the quartic mathematical model was implemented [<xref ref-type="bibr" rid="ref-21">21</xref>].</p>
<p>The function formulas between droplet velocity or pressure in the infusion process and the liquid surface are as following <xref ref-type="disp-formula" rid="eqn-6">Eqs. (6)</xref> and <xref ref-type="disp-formula" rid="eqn-7">(7)</xref> show:</p>
<p><disp-formula id="eqn-9">
<label>(9)</label>
<!--<alternatives>
<graphic mimetype="image" mime-subtype="png" xlink:href="eqn-9.png"/><tex-math id="tex-eqn-9"><![CDATA[V = {p_1} \times {h^4} + {p_2} \times {h^3} + {p_3} \times {h^2} + {p_4} \times h + {p_5}]]></tex-math>--><mml:math id="mml-eqn-9" display="block"><mml:mi>V</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mo>&#x00D7;</mml:mo><mml:mrow><mml:msup><mml:mi>h</mml:mi><mml:mn>4</mml:mn></mml:msup></mml:mrow><mml:mo>&#x002B;</mml:mo><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mo>&#x00D7;</mml:mo><mml:mrow><mml:msup><mml:mi>h</mml:mi><mml:mn>3</mml:mn></mml:msup></mml:mrow><mml:mo>&#x002B;</mml:mo><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow><mml:mo>&#x00D7;</mml:mo><mml:mrow><mml:msup><mml:mi>h</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mo>&#x002B;</mml:mo><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mn>4</mml:mn></mml:msub></mml:mrow><mml:mo>&#x00D7;</mml:mo><mml:mi>h</mml:mi><mml:mo>&#x002B;</mml:mo><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mn>5</mml:mn></mml:msub></mml:mrow></mml:math>
<!--</alternatives>--></disp-formula></p>
<p><disp-formula id="eqn-10">
<label>(10)</label>
<!--<alternatives>
<graphic mimetype="image" mime-subtype="png" xlink:href="eqn-10.png"/><tex-math id="tex-eqn-10"><![CDATA[P = {P_1} * {h^3} + {P_2} * {h^2} + {P_3}h + {P_4}]]></tex-math>--><mml:math id="mml-eqn-10" display="block"><mml:mi>P</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mo>&#x2217;</mml:mo><mml:mrow><mml:msup><mml:mi>h</mml:mi><mml:mn>3</mml:mn></mml:msup></mml:mrow><mml:mo>&#x002B;</mml:mo><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mo>&#x2217;</mml:mo><mml:mrow><mml:msup><mml:mi>h</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mo>&#x002B;</mml:mo><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow><mml:mi>h</mml:mi><mml:mo>&#x002B;</mml:mo><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mn>4</mml:mn></mml:msub></mml:mrow></mml:math>
<!--</alternatives>--></disp-formula></p>
<p>h is dropping height, V is droplet velocity, P is pressure in soft bags.</p>
<p>The formal <xref ref-type="disp-formula" rid="eqn-8">Eqs. (8)</xref> and <xref ref-type="disp-formula" rid="eqn-9">(9)</xref> were obtained by least squares polynomial fitting:</p>
<p><disp-formula id="eqn-11">
<label>(11)</label>
<!--<alternatives>
<graphic mimetype="image" mime-subtype="png" xlink:href="eqn-11.png"/><tex-math id="tex-eqn-11"><![CDATA[\left[ {\matrix{ {20} &#9; {\sum\nolimits_{i = 0}^{19} {{h_i}} } &#9; {\sum\nolimits_{i = 0}^{19} {{h_i}^2} } &#9; {\sum\nolimits_{i = 0}^{19} {{h_i}^3} } &#9; {{{\sum\nolimits_{i = 0}^{19} {{h_i}} }^4}} \cr {\sum\nolimits_{i = 0}^{19} {{h_i}} } &#9; {\sum\nolimits_{i = 0}^{19} {{h_i}^2} } &#9; {\sum\nolimits_{i = 0}^{19} {{h_i}^3} } &#9; {{{\sum\nolimits_{i = 0}^{19} {{h_i}} }^4}} &#9; {{{\sum\nolimits_{i = 0}^{19} {{h_i}} }^5}} \cr {\sum\nolimits_{i = 0}^{19} {{h_i}^2} } &#9; {\sum\nolimits_{i = 0}^{19} {{h_i}^3} } &#9; {{{\sum\nolimits_{i = 0}^{19} {{h_i}} }^4}} &#9; {{{\sum\nolimits_{i = 0}^{19} {{h_i}} }^5}} &#9; {{{\sum\nolimits_{i = 0}^{19} {{h_i}} }^6}} \cr {\sum\nolimits_{i = 0}^{19} {{h_i}^3} } &#9; {{{\sum\nolimits_{i = 0}^{19} {{h_i}} }^4}} &#9; {{{\sum\nolimits_{i = 0}^{19} {{h_i}} }^5}} &#9; {{{\sum\nolimits_{i = 0}^{19} {{h_i}} }^6}} &#9; {{{\sum\nolimits_{i = 0}^{19} {{h_i}} }^7}} \cr {{{\sum\nolimits_{i = 0}^{19} {{h_i}} }^4}} &#9; {{{\sum\nolimits_{i = 0}^{19} {{h_i}} }^5}} &#9; {{{\sum\nolimits_{i = 0}^{19} {{h_i}} }^6}} &#9; {{{\sum\nolimits_{i = 0}^{19} {{h_i}} }^7}} &#9; {{{\sum\nolimits_{i = 0}^m {{h_i}} }^8}} \cr } } \right] \times \left( {\matrix{ {{P_5}} \cr {{P_4}} \cr {{P_3}} \cr {{P_2}} \cr {{P_1}} \cr } } \right) = \left( {\matrix{ {\sum\nolimits_{i = 0}^{19} {{v_i}} } \cr {\sum\nolimits_{i = 0}^{19} {{h_i}{v_i}} } \cr {\sum\nolimits_{i = 0}^{19} {{h_i}^2{v_i}} } \cr {\sum\nolimits_{i = 0}^{19} {{h_i}^3{v_i}} } \cr {\sum\nolimits_{i = 0}^{19} {{h_i}^4{v_i}} } \cr } } \right)]]></tex-math>--><mml:math id="mml-eqn-11" display="block"><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mtable columnspacing="1em" rowspacing="4pt"><mml:mtr><mml:mtd><mml:mrow><mml:mn>20</mml:mn></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>3</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow><mml:mn>4</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>3</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow><mml:mn>4</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow><mml:mn>5</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>3</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow><mml:mn>4</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow><mml:mn>5</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow><mml:mn>6</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>3</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow><mml:mn>4</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow><mml:mn>5</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow><mml:mn>6</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow><mml:mn>7</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow><mml:mn>4</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow><mml:mn>5</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow><mml:mn>6</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow><mml:mn>7</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>m</mml:mi></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow><mml:mn>8</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>&#x00D7;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtable columnspacing="1em" rowspacing="4pt"><mml:mtr><mml:mtd><mml:mrow><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mn>5</mml:mn></mml:msub></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mn>4</mml:mn></mml:msub></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x003D;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtable columnspacing="1em" rowspacing="4pt"><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>v</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>v</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>2</mml:mn></mml:msup><mml:mrow><mml:msub><mml:mi>v</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>3</mml:mn></mml:msup><mml:mrow><mml:msub><mml:mi>v</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>4</mml:mn></mml:msup><mml:mrow><mml:msub><mml:mi>v</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:math>
<!--</alternatives>--></disp-formula></p>
<p><disp-formula id="eqn-12">
<label>(12)</label>
<!--<alternatives>
<graphic mimetype="image" mime-subtype="png" xlink:href="eqn-12.png"/><tex-math id="tex-eqn-12"><![CDATA[\left[ {\matrix{ {20} &#9; {\sum\nolimits_{i = 0}^{19} {{h_i}} } &#9; {\sum\nolimits_{i = 0}^{19} {{h_i}^2} } &#9; {\sum\nolimits_{i = 0}^{19} {{h_i}^3} } \cr {\sum\nolimits_{i = 0}^{19} {{h_i}} } &#9; {\sum\nolimits_{i = 0}^{19} {{h_i}^2} } &#9; {\sum\nolimits_{i = 0}^{19} {{h_i}^3} } &#9; {\sum\nolimits_{i = 0}^{19} {{h_i}^4} } \cr {\sum\nolimits_{i = 0}^{19} {{h_i}^2} } &#9; {\sum\nolimits_{i = 0}^{19} {{h_i}^3} } &#9; {\sum\nolimits_{i = 0}^{19} {{h_i}^4} } &#9; {\sum\nolimits_{i = 0}^{19} {{h_i}^5} } \cr {\sum\nolimits_{i = 0}^{19} {{h_i}^3} } &#9; {\sum\nolimits_{i = 0}^{19} {{h_i}^4} } &#9; {\sum\nolimits_{i = 0}^{19} {{h_i}^5} } &#9; {\sum\nolimits_{i = 0}^{19} {{h_i}^6} } \cr } } \right] \times \left( {\matrix{ {{b_4}} \cr {{b_3}} \cr {{b_2}} \cr {{b_1}} \cr } } \right) = \left( {\matrix{ {\sum\nolimits_{i = 0}^{19} {{P_i}} } \cr {\sum\nolimits_{i = 0}^{19} {{h_i}{P_i}} } \cr {\sum\nolimits_{i = 0}^{19} {{h_i}{P_i}^2} } \cr {\sum\nolimits_{i = 0}^{19} {{h_i}{P_i}^3} } \cr } } \right)]]></tex-math>--><mml:math id="mml-eqn-12" display="block"><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mtable columnspacing="1em" rowspacing="4pt"><mml:mtr><mml:mtd><mml:mrow><mml:mn>20</mml:mn></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>3</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>3</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>4</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>3</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>4</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>5</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>3</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>4</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>5</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>6</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>&#x00D7;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtable columnspacing="1em" rowspacing="4pt"><mml:mtr><mml:mtd><mml:mrow><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mn>4</mml:mn></mml:msub></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x003D;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtable columnspacing="1em" rowspacing="4pt"><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>19</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow><mml:mn>3</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:math>
<!--</alternatives>--></disp-formula></p>
<p>The coefficients p<sub>1</sub>, p<sub>2</sub>, p<sub>3</sub>, p<sub>4</sub>, p<sub>5</sub> and b<sub>0</sub>, b<sub>1</sub>, b<sub>2</sub>, b<sub>3</sub> can be obtained by solving the above-mentioned normal equations.</p>
<p>MATLAB was applied to build the mathematical model between the 250 ml polypropylene infusion bag droplet velocity, pressure and liquid level:</p>
<p>(1) The level of drop heights, pressure values and drop rate values were put into MATLAB software:</p>
<p>Height: <inline-formula id="ieqn-1">
<!--<alternatives><inline-graphic xlink:href="ieqn-1.tif"/><tex-math id="tex-ieqn-1"><![CDATA[h = \left[ {0,0.5,1,1.5,2,2.5,3,3.5,4,4.5,5,5.5,6,6.5,7,7.5,8,8.5,9,9.5} \right]]]></tex-math>--><mml:math id="mml-ieqn-1"><mml:mi>h</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mn>0.5</mml:mn><mml:mo>,</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>1.5</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mn>2.5</mml:mn><mml:mo>,</mml:mo><mml:mn>3</mml:mn><mml:mo>,</mml:mo><mml:mn>3.5</mml:mn><mml:mo>,</mml:mo><mml:mn>4</mml:mn><mml:mo>,</mml:mo><mml:mn>4.5</mml:mn><mml:mo>,</mml:mo><mml:mn>5</mml:mn><mml:mo>,</mml:mo><mml:mn>5.5</mml:mn><mml:mo>,</mml:mo><mml:mn>6</mml:mn><mml:mo>,</mml:mo><mml:mn>6.5</mml:mn><mml:mo>,</mml:mo><mml:mn>7</mml:mn><mml:mo>,</mml:mo><mml:mn>7.5</mml:mn><mml:mo>,</mml:mo><mml:mn>8</mml:mn><mml:mo>,</mml:mo><mml:mn>8.5</mml:mn><mml:mo>,</mml:mo><mml:mn>9</mml:mn><mml:mo>,</mml:mo><mml:mn>9.5</mml:mn></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:math>
<!--</alternatives>--></inline-formula></p>
<p>Drip: <inline-formula id="ieqn-2">
<!--<alternatives><inline-graphic xlink:href="ieqn-2.tif"/><tex-math id="tex-ieqn-2"><![CDATA[V = \left[ {57,55,55,55,55,54,54,54,53,50,51,51,50,49,49,47,45,44,42,39} \right]]]></tex-math>--><mml:math id="mml-ieqn-2"><mml:mi>V</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mn>57</mml:mn><mml:mo>,</mml:mo><mml:mn>55</mml:mn><mml:mo>,</mml:mo><mml:mn>55</mml:mn><mml:mo>,</mml:mo><mml:mn>55</mml:mn><mml:mo>,</mml:mo><mml:mn>55</mml:mn><mml:mo>,</mml:mo><mml:mn>54</mml:mn><mml:mo>,</mml:mo><mml:mn>54</mml:mn><mml:mo>,</mml:mo><mml:mn>54</mml:mn><mml:mo>,</mml:mo><mml:mn>53</mml:mn><mml:mo>,</mml:mo><mml:mn>50</mml:mn><mml:mo>,</mml:mo><mml:mn>51</mml:mn><mml:mo>,</mml:mo><mml:mn>51</mml:mn><mml:mo>,</mml:mo><mml:mn>50</mml:mn><mml:mo>,</mml:mo><mml:mn>49</mml:mn><mml:mo>,</mml:mo><mml:mn>49</mml:mn><mml:mo>,</mml:mo><mml:mn>47</mml:mn><mml:mo>,</mml:mo><mml:mn>45</mml:mn><mml:mo>,</mml:mo><mml:mn>44</mml:mn><mml:mo>,</mml:mo><mml:mn>42</mml:mn><mml:mo>,</mml:mo><mml:mn>39</mml:mn></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:math>
<!--</alternatives>--></inline-formula></p>
<p>Pressure:</p>
<p><disp-formula id="eqn-112">
<!--<alternatives>
<graphic mimetype="image" mime-subtype="png" xlink:href="eqn-112.png"/><tex-math id="tex-eqn-112"><![CDATA[P = \left[ {\matrix{ { - 442,} &#9; { - 920,} &#9; { - 920,} &#9; { - 940,} &#9; { - 890,} &#9; { - 965,} &#9; { - 932,} &#9; { - 932,} &#9; { - 973,} &#9; { - 970,} &#9; { - 1030,} \cr { - 1085,} &#9; { - 1178,} &#9; { - 1263,} &#9; { - 1363,} &#9; { - 1507,} &#9; { - 1847,} &#9; { - 2015,} &#9; { - 2095,} &#9; { - 2220} &#9; {} &#9; {} \cr } } \right]]]></tex-math>--><mml:math id="mml-eqn-112"><mml:mi>P</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mtable columnspacing="1em" rowspacing="4pt"><mml:mtr><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>442</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>920</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>920</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>940</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>890</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>965</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>932</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>932</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>973</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>970</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1030</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1085</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1178</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1263</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1363</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1507</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1847</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>2015</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>2095</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>2220</mml:mn></mml:mrow></mml:mtd><mml:mtd><mml:mrow></mml:mrow></mml:mtd><mml:mtd><mml:mrow></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:math>
<!--</alternatives>--></disp-formula></p>
<p>(2) The correct function was selected to guarantee the corresponding point on the function image as much as possible. So the formula for the 4 power function of velocity and liquid drop height is as the following:</p>
<p><disp-formula id="eqn-13">
<label>(13)</label>
<!--<alternatives>
<graphic mimetype="image" mime-subtype="png" xlink:href="eqn-13.png"/><tex-math id="tex-eqn-13"><![CDATA[V = {p_1} \times {h^4} + {p_2} \times {h^3} + {p_3} \times {h^2} + {p_4} \times h + {p_5}]]></tex-math>--><mml:math id="mml-eqn-13" display="block"><mml:mi>V</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mo>&#x00D7;</mml:mo><mml:mrow><mml:msup><mml:mi>h</mml:mi><mml:mn>4</mml:mn></mml:msup></mml:mrow><mml:mo>&#x002B;</mml:mo><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mo>&#x00D7;</mml:mo><mml:mrow><mml:msup><mml:mi>h</mml:mi><mml:mn>3</mml:mn></mml:msup></mml:mrow><mml:mo>&#x002B;</mml:mo><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow><mml:mo>&#x00D7;</mml:mo><mml:mrow><mml:msup><mml:mi>h</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mo>&#x002B;</mml:mo><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mn>4</mml:mn></mml:msub></mml:mrow><mml:mo>&#x00D7;</mml:mo><mml:mi>h</mml:mi><mml:mo>&#x002B;</mml:mo><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mn>5</mml:mn></mml:msub></mml:mrow></mml:math>
<!--</alternatives>--></disp-formula></p>
<p>And the cubic function of pressure and liquid drop height is as the following:</p>
<p><disp-formula id="eqn-14">
<label>(14)</label>
<!--<alternatives>
<graphic mimetype="image" mime-subtype="png" xlink:href="eqn-14.png"/><tex-math id="tex-eqn-14"><![CDATA[P = {P_1} * {h^3} + {P_2} * {h^2} + {P_3}h + {P_4}]]></tex-math>--><mml:math id="mml-eqn-14" display="block"><mml:mi>P</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mo>&#x2217;</mml:mo><mml:mrow><mml:msup><mml:mi>h</mml:mi><mml:mn>3</mml:mn></mml:msup></mml:mrow><mml:mo>&#x002B;</mml:mo><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mo>&#x2217;</mml:mo><mml:mrow><mml:msup><mml:mi>h</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mo>&#x002B;</mml:mo><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow><mml:mi>h</mml:mi><mml:mo>&#x002B;</mml:mo><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mn>4</mml:mn></mml:msub></mml:mrow></mml:math>
<!--</alternatives>--></disp-formula></p>
<p>(3) The functions are output to obtain the mathematical relationship of the height drop change, which are as the following:</p>
<p><disp-formula id="eqn-15">
<label>(15)</label>
<!--<alternatives>
<graphic mimetype="image" mime-subtype="png" xlink:href="eqn-15.png"/><tex-math id="tex-eqn-15"><![CDATA[P = - 4.4{h^3} + 41.83{h^2} - 178.7h - 674.1]]></tex-math>--><mml:math id="mml-eqn-15" display="block"><mml:mi>P</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mo>&#x2212;</mml:mo><mml:mn>4.4</mml:mn><mml:mrow><mml:msup><mml:mi>h</mml:mi><mml:mn>3</mml:mn></mml:msup></mml:mrow><mml:mo>&#x002B;</mml:mo><mml:mn>41.83</mml:mn><mml:mrow><mml:msup><mml:mi>h</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>178.7</mml:mn><mml:mi>h</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>674.1</mml:mn></mml:math>
<!--</alternatives>--></disp-formula></p>
<p><disp-formula id="eqn-16">
<label>(16)</label>
<!--<alternatives>
<graphic mimetype="image" mime-subtype="png" xlink:href="eqn-16.png"/><tex-math id="tex-eqn-16"><![CDATA[V = - 0.8668{h^4} - 0.7227{h^3} + 2.739{h^2} + 5.983h + 49.45]]></tex-math>--><mml:math id="mml-eqn-16" display="block"><mml:mi>V</mml:mi><mml:mo>&#x003D;</mml:mo><mml:mo>&#x2212;</mml:mo><mml:mn>0.8668</mml:mn><mml:mrow><mml:msup><mml:mi>h</mml:mi><mml:mn>4</mml:mn></mml:msup></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>0.7227</mml:mn><mml:mrow><mml:msup><mml:mi>h</mml:mi><mml:mn>3</mml:mn></mml:msup></mml:mrow><mml:mo>&#x002B;</mml:mo><mml:mn>2.739</mml:mn><mml:mrow><mml:msup><mml:mi>h</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mo>&#x002B;</mml:mo><mml:mn>5.983</mml:mn><mml:mi>h</mml:mi><mml:mo>&#x002B;</mml:mo><mml:mn>49.45</mml:mn></mml:math>
<!--</alternatives>--></disp-formula></p>
</sec>
<sec id="s4_4">
<label>4.4</label>
<title>Effect of Needle Diameter on Droplet Velocity</title>
<p>At the beginning of the experiment, the same infusion tube was used in the first several groups until that the initial droplet velocity was found diverse. Because the needle diameter was so small that long glucose infusion will leave the glucose crystal remain in the needle diameter, leading to various experimental data. Therefore, new tubes should be replaced every time. There are various factors that can impact the droplet velocity such as the hanging height, level height, material, structure, deformation of infusion bag and the pipeline loss along infusion process, but according to the experimental analysis, the most influential factor is the needle diameter, which is the needle local piping loss.</p>
<p>It can be seen from <xref ref-type="table" rid="table-3">Tab. 3</xref>, in the range of the errors permitted, when the needle diameter was same, the droplet velocity increased with the increase of the infusion bag capacity; for the same specifications of the infusion bag, while the needle diameter increases, drip significantly changes with the increase in the diameter of the needle. Therefore, the relationship between needle diameter capacity and droplet velocity will be concluded from <xref ref-type="table" rid="table-3">Tab. 3</xref>.</p>

<table-wrap id="table-3">
<label>Table 3</label>
<caption>
<title>Effect of needle diameter on droplet velocity</title>
</caption>
<table>
<colgroup>
<col/>
<col/>
<col/>
<col/>
<col/>
</colgroup>
<thead>
<tr>
<th>Needle diameter (mm)</th>
<th>100 ml</th>
<th>250 ml</th>
<th>500 ml</th>
<th>Average</th>
</tr>
</thead>
<tbody>
<tr>
<td>0.5 mm</td>
<td>47</td>
<td>55</td>
<td>57</td>
<td>52</td>
</tr>
<tr>
<td>0.6 mm</td>
<td>69</td>
<td>71</td>
<td>73</td>
<td>71</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn id="table-3fn1" fn-type="other">
<p>Note: The statistics of the droplet velocity are for the initial droplet velocity, unit: drop/min.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
</sec>
<sec id="s5">
<label>5</label>
<title>Discussion</title>
<p>With the development of medicine, sealing and economy, the applicability of medicines&#x2019; storage containers have become higher and higher, and the vertical infusion bags have been put on the medical stage. Therefore, the analysis of the vertical infusion bag has become the focus of many scholars and drives them to conduct more researches on the thickness profile of the vertical infusion bag. Under the condition of meeting the required drip rate, there is a lack of research on optimal control of the gas-liquid ratio in the bottle. Certain pharmaceutical companies also companies are manual controlling the gas-ratio, so that some medicine bottles are not used up. Therefore, through the model analysis and the flow field analysis of the vertical infusion bag, the gas-liquid ratio in the vertical infusion bag can be simulated, which can effectively achieve the energy-saving effects and effectively reduce its unusable volume. The possibility of pre-simulation makes it easy to provide pharmaceutical companies with certain reliable data and put them into use. The fully enclosed infusion has been proven effectively in avoiding contamination of infusions by bacterial, viruses, and dust in the air, improving the safety and effectiveness of infusion, and gradually becoming the mainstream of the infusion industry [<xref ref-type="bibr" rid="ref-22">22</xref>,<xref ref-type="bibr" rid="ref-23">23</xref>].</p>
<p>There are a large number of researches and practical application cases on the vertical polypropylene infusion bag. For example, Liu Gexin, a bioengineering and clinical researcher, proposed that the wall thickness of the infusion bag could affect the infusion performance. The contour of Von Mises stress distribution showed that under the same capacity condition and VON Mises stress change, the Von Mises stress changes insignificantly when the wall thickness of 500 mL bag is 0.2 mm and 0.3 mm. When the wall thickness is 0.3 mm and 0.5 mm, the maximum Von Mises stress reduces by 28%. The maximum Von Mises stress area of each bag is on the side wall of the bag. The actual dosing space has a relationship with the thickness of the vertical infusion bag wall. The wall thickness of the 100 ml bag type is from 0.2 mm to 0.5 mm, and the actual dosing space is reduced by 27%, respectively. The actual dosing space of the 250 mL and 500 mL bags is reduced by 54% and 49%, respectively. The greater the wall thickness is, the larger the actual dosing space of the vertical infusion bag. The wall thickness also has a certain influence on the emptying rate. The vertical infusion bag indicates that when the wall thickness of the bag becomes from 0.2 mm to 0.5 mm, the emptying rate is 11%. Under the same wall thickness, the emptying rate of high-standard volume vertical infusion bag is significantly higher than that of low-standard one [<xref ref-type="bibr" rid="ref-24">24</xref>,<xref ref-type="bibr" rid="ref-25">25</xref>]. The final determination of the wall thickness also needs to consider factors such as the air permeability, moisture permeability and puncture resistance of the vertical infusion bag. Liu Gexin and other bioengineering and clinical researchers did not consider it.</p>
<p>Wang Huasheng explained in the Chinese Journal of Infectious Diseases that the fully closed infusion method used the pressure of atmospheric pressure to directly act on the wall of the bag. The bag body shrank automatically to allow the liquid enter the human body by itself, avoiding the contact between the infusion and the air, and preventing the contamination of the infusion by bacteria, viruses, and dust in the air. The vertical polypropylene infusion bag had almost no adsorption on insulin nitroglycerin and paclitaxel, which ensured that the infusion was safe and effective.</p>
<p>Chen Wei, from the Department of Pharmacy, Tongji Hospital affiliated to Huazhong University of Science and Technology, examined the clinical performance of the vertical polypropylene infusion bag, and explained the vertical polypropylene infusion bag was completely transparent and had a high light transmittance which was convenient for clinical inspection to ensure the safety of medication. The maximum dosing space was greatly increased, which was convenient for clinical use [<xref ref-type="bibr" rid="ref-26">26</xref>]. The maximum dosing space of different infusion devices was different, the maximum dosing space of the vertical polypropylene infusion bag is the largest, the maximum dosing space of the soft bag is second and the maximum dosing space of the glass bottle is the smallest. The different packages of each infusion specification are compared in <xref ref-type="table" rid="table-4">Tab. 4</xref>, and the difference is statistically significant.</p>

<table-wrap id="table-4">
<label>Table 4</label>
<caption>
<title>Dosing space of three different packaging infusion devices</title>
</caption>
<table>
<colgroup>
<col/>
<col/>
<col/>
<col/>
</colgroup>
<thead>
<tr><th rowspan="2">Package style</th><th colspan="3">Volume</th>
</tr>
<tr>
<th>100 ml</th>
<th>250 ml</th>
<th>500 ml</th>
</tr>
</thead>
<tbody>
<tr>
<td>Soft bags</td>
<td>53.41 &#x00B1; 1.11</td>
<td>86.95 &#x00B1; 0.72</td>
<td>134.89 &#x00B1; 2.11</td>
</tr>
<tr>
<td>Glass bottles</td>
<td>19.08 &#x00B1; 0.55</td>
<td>44.60 &#x00B1; 0.66</td>
<td>77.27 &#x00B1; 0.85</td>
</tr>
<tr>
<td>Vertical infusion bags</td>
<td>106.95 &#x00B1; 1.63</td>
<td>124.20 &#x00B1; 1.23</td>
<td>181.09 &#x00B1; 2.23</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The vertical infusion bag is especially convenient to add large-volume traditional Chinese medicine injections; its packing material is extremely flexible to withstand temperature changes, and can withstand a high maximum pressure, which is convenient for transportation and storage. The material of the infusion bag can be upright like a glass bottle, which is in line with the clinical operation habits of nursing staff and improves the work efficiency. The vertical polypropylene infusion bag integrates all the advantages of glass bottles, plastic bottle, and non-soft bags and it has extremely high clinical value. Compared with soft bags, the maximum pressure that the vertical polypropylene infusion bag can withstand is significantly higher. Pressure bearing capacity of three different packaging infusion devices are compared in <xref ref-type="table" rid="table-5">Tab. 5</xref>, and the difference is statistically significant.</p>

<table-wrap id="table-5">
<label>Table 5</label>
<caption>
<title>Pressure bearing capacity of three different packaging infusion devices</title>
</caption>
<table>
<colgroup>
<col/>
<col/>
<col/>
<col/>
</colgroup>
<thead>
<tr><th rowspan="2">Package style</th><th colspan="3">Volume</th>
</tr>
<tr>
<th>100 ml</th>
<th>250 ml</th>
<th>500 ml</th>
</tr>
</thead>
<tbody>
<tr>
<td>Soft bags</td>
<td>0.251 &#x00B1; 0.01</td>
<td>0.16 &#x00B1; 0.007</td>
<td>0.13 &#x00B1; 0.005</td>
</tr>
<tr>
<td>Vertical infusion bags</td>
<td>0.78 &#x00B1; 0.09</td>
<td>1.01 &#x00B1; 0.027</td>
<td>0.68 &#x00B1; 0.052</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The quality of the vertical polypropylene infusion bag was investigated in the Journal of Hunan University of Traditional Chinese Medicine. Taking the vertical polypropylene infusion bag as the research object, comparing the multi-layer co-extruded infusion bag and the polypropylene infusion bottle, the light transmittance, adaptability, and adsorption of them were studied. Among them, the vertical polypropylene infusion bag and the polypropylene infusion bottle are mainly made of polypropylene; the material of the multi-layer co-extruded infusion bag is composed of polyethylene, polyamine, and polypropylene. In clinical applications, polypropylene infusion bottles, like glass infusion bottles, are semi-open infusion containers, and vertical polypropylene infusion bags and multi-layer co-extruded infusion bags are fully enclosed infusion containers. In the light transmittance test, the vertical polypropylene infusion bag, the multi-layer co-extruded infusion bag and the polypropylene infusion bottle were detected by an ultraviolet spectrophotometer. It is generally believed that the transmittance of the vertical polypropylene infusion bags and the multi-layer co-extruded infusion bags is greater than 75.0%, and the test samples are of good quality, all of which are higher than the standard. The average transmittance of the three was (95.20 &#x00B1; 0.55)%, (91.96 &#x00B1; 0.84)%, and (93.78 &#x00B1; 0.37)% (<italic>P</italic> &#x003C; 0.05). The transmittance of the vertical polypropylene infusion bags is better than that of the multi-layer co-extruded infusion bags. Comparing the multi-layer co-extruded infusion bags and the polypropylene infusion bottles, the light transmittance, adaptability and adsorption of them were investigated. The vertical polypropylene infusion bag is a new type of patented infusion packaging container. Compared with the multi-layer co-extruded infusion bag, it is not easy to be broken and easy to transport; it has good transparency, small deformation, and good adaptability; its material is chemically stable and is suitable for drugs; it can be placed upright, which is in line with operating habits. In addition, compared with the polypropylene infusion bottle, it does not need to introduce air during use, effectively avoiding secondary pollution [<xref ref-type="bibr" rid="ref-27">27</xref>&#x2013;<xref ref-type="bibr" rid="ref-30">30</xref>].</p>
</sec>
<sec id="s6">
<label>6</label>
<title>Conclusions and Prospect</title>
<list list-type="order">
<list-item><p>According to the principle of fluid mechanics, the theoretical exhaust pressure at the end of the drip rate of 40 drops/min was &#x2013;2,640 Pa.</p></list-item>
<list-item><p>The results simulated by CFD showed the maximum speed inside the needle was at the position of axis, approximately 0.5 m/s, and the speed near the wall gradually decreased to zero. However, from the actual situation, the material of the infusion bag and the deformation of the force also have an effect on the drip rate. The necessary coupling analysis of these parameters can make the graphics and data more comprehensive and accurate.</p></list-item>
<list-item><p>A human blood vessel vein simulator and measuring devices were used to obtain relevant data about pressure and drop rate. The least squares method and MATLAB software were used to simulate the relation of inside pressure, fluid speed and liquid height. The inside pressure and liquid height were in the cubic function relation while the drop rate and the liquid height were in the quadric function relation.</p></list-item>
<list-item><p>Through the experiments, it could be concluded that the needle diameter was one of the most effective factors among all influencing factors. However, the force analysis of the infusion bag has not been carried out by ANSYS software. If the analysis by ANYSY software is combined with experiments and simulations, the transient changes of the infusion process can be better reflected, the data can be more accurate, and the simulation results can be more convincing.</p></list-item>
</list>
</sec>
</body>
<back><fn-group>
<fn fn-type="other">
<p><bold>Funding Statement:</bold> This article belongs to the project of &#x201C;The University Synergy Innovation Program of Anhui Province (GXXT-2019-004)&#x201D;; &#x201C;Teaching Research Project of Anhui Education Department (2019jyxm0229)&#x201D;.</p>
</fn>
<fn fn-type="conflict">
<p><bold>Conflicts of Interest:</bold> The authors declare that they have no conflicts of interest to report regarding the present study.</p>
</fn>
</fn-group>
<ref-list content-type="authoryear">
<title>References</title>
<ref id="ref-1">
<label>1.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Liu</surname>, 
<given-names>G. H.</given-names></string-name>
</person-group> (
<year>2009</year>). 
<article-title>Factors affecting the infusion rate and its adjustment method</article-title>. 
<source>Nursing Practice and Research</source><italic>,</italic> 
<volume>6</volume><italic>(</italic>
<issue>1</issue><italic>),</italic> 
<fpage>88</fpage>&#x2013;
<lpage>90</lpage>.</mixed-citation>
</ref>
<ref id="ref-2">
<label>2.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Kinna</surname>, 
<given-names>S.</given-names></string-name>, <string-name>
<surname>Jayshree</surname>, 
<given-names>T.</given-names></string-name>, <string-name>
<surname>Dushyant</surname>, 
<given-names>V.</given-names></string-name>
</person-group> (
<year>2020</year>). 
<article-title>Comparison of glottis views with fixed-height pillow <italic>vs.</italic> adjustable-pillow height by pressure infusion bag for successful intubation</article-title>. 
<source>Bali Journal of Anesthesiology</source><italic>,</italic> 
<volume>4</volume><italic>(</italic>
<issue>2</issue><italic>),</italic> <fpage>35</fpage>.</mixed-citation>
</ref>
<ref id="ref-3">
<label>3.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Ranjana</surname>, 
<given-names>B.</given-names></string-name>, <string-name>
<surname>Jyoti</surname>, 
<given-names>P.</given-names></string-name>, <string-name>
<surname>Anurag</surname>, 
<given-names>K.</given-names></string-name>
</person-group> (
<year>2019</year>). 
<article-title>Pharmacokinetic-pharmacodynamic (PK-PD) modelling of effect of naringenin &#x0026; its surface modified nanocarriers on associated &#x0026; core behaviours of autism spectrum disorders (ASD)</article-title>. 
<source>Planta Medica International Open</source><italic>,</italic> 
<volume>6</volume><italic>(</italic>
<issue>2</issue><italic>),</italic> <fpage>e41</fpage>&#x2013;<lpage>e49</lpage>.</mixed-citation>
</ref>
<ref id="ref-4">
<label>4.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Matey</surname>, 
<given-names>J. M.</given-names></string-name>, <string-name>
<surname>Moreno de Simon</surname>, 
<given-names>M. D.</given-names></string-name>, <string-name>
<surname>Garc&#x00ED;a-Ruiz</surname>, 
<given-names>C.</given-names></string-name>, <string-name>
<surname>Montalvo</surname>, 
<given-names>G.</given-names></string-name>
</person-group> (
<year>2019</year>). 
<article-title>A validated GC-MS method for ketamine and norketamine in hair and its use in authentic cases</article-title>. 
<source>Forensic Science International</source><italic>,</italic> 
<volume>301</volume><italic>,</italic> <fpage>447</fpage>&#x2013;<lpage>454</lpage>.</mixed-citation>
</ref>
<ref id="ref-5">
<label>5.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Chen</surname>, 
<given-names>X.</given-names></string-name>, <string-name>
<surname>Min</surname>, 
<given-names>L. H.</given-names></string-name>
</person-group> (
<year>2020</year>). 
<article-title>Design and application of a new type of double-headed infusion set</article-title>. 
<source>Chinese Clinical Nursing</source><italic>,</italic> 
<volume>12</volume><italic>(</italic>
<issue>1</issue><italic>),</italic> 
<fpage>92</fpage>.</mixed-citation>
</ref>
<ref id="ref-6">
<label>6.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Dar&#x00E1;s Orenga</surname>, 
<given-names>M. &#x00C1;.</given-names></string-name>, <string-name>
<surname>Gellida Vilarroig</surname>, 
<given-names>M.</given-names></string-name>, <string-name>
<surname>Vives Piqueres</surname>, 
<given-names>L.</given-names></string-name>, <string-name>
<surname>Sanz Garc&#x00ED;a</surname>, 
<given-names>M.</given-names></string-name>, <string-name>
<surname>Inoges</surname>, 
<given-names>R. T.</given-names></string-name> <etal>et al.</etal>
</person-group> (
<year>2019</year>). 
<article-title>Comparative study of different epidural infusion sets at maximum speeds for labor analgesia</article-title>. 
<source>Revista Espa&#x00F1;ola de Anestesiolog&#x00ED;a y Reanimaci&#x00F3;n</source><italic>,</italic> 
<volume>67</volume><italic>(</italic>
<issue>4</issue><italic>),</italic> <fpage>179</fpage>&#x2013;<lpage>186</lpage>.</mixed-citation>
</ref>
<ref id="ref-7">
<label>7.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Liang</surname>, 
<given-names>Y. Y.</given-names></string-name>, <string-name>
<surname>Wu</surname>, 
<given-names>J. T.</given-names></string-name>, <string-name>
<surname>Lan</surname>, 
<given-names>Y. N.</given-names></string-name>, <string-name>
<surname>Lan</surname>, <given-names>Q.</given-names></string-name>
</person-group> (
<year>2012</year>). 
<article-title>Application of pressurized infusion bag in pediatric cardiac arrest fluid infusion</article-title>. 
<source>Nursing Practice and Research</source><italic>,</italic> 
<volume>9</volume><italic>(</italic>
<issue>19</issue><italic>),</italic> 
<fpage>46</fpage>&#x2013;
<lpage>47</lpage>.</mixed-citation>
</ref>
<ref id="ref-8">
<label>8.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Song</surname>, 
<given-names>Z. F.</given-names></string-name>
</person-group> (
<year>2010</year>). 
<article-title>Intravenous infusion rate rapid determination method</article-title>. 
<source>Chinese Journal of Traumatology</source><italic>,</italic> 
<volume>18</volume><italic>(</italic>
<issue>2</issue><italic>),</italic> 
<fpage>126</fpage>.</mixed-citation>
</ref>
<ref id="ref-9">
<label>9.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Xiao</surname>, 
<given-names>Z.</given-names></string-name>, <string-name>
<surname>Liu</surname>, 
<given-names>J. S.</given-names></string-name>, <string-name>
<surname>Wang</surname>, <given-names>L.</given-names></string-name>, <string-name>
<surname>Zhang</surname>, <given-names>T. N.</given-names></string-name>
</person-group> (
<year>2011</year>). 
<article-title>Design and implementation of real-time monitoring system for infusion drip rate</article-title>. 
<source>Studies</source><italic>,</italic> 
<volume>26</volume><italic>(</italic>
<issue>11</issue><italic>),</italic> 
<fpage>19</fpage>&#x2013;
<lpage>21</lpage>.</mixed-citation>
</ref>
<ref id="ref-10">
<label>10.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Wang</surname>, 
<given-names>G. S.</given-names></string-name>
</person-group> (
<year>2019</year>). 
<article-title>Fluid-structure coupling simulation and experimental verification of compression process of multi-chamber infusion soft bag</article-title>. 
<source>Journal of Medical Biomechanics</source><italic>,</italic> 
<volume>34</volume><italic>(</italic>
<issue>S1</issue><italic>),</italic> 
<fpage>183</fpage>.</mixed-citation>
</ref>
<ref id="ref-11">
<label>11.</label><mixed-citation publication-type="book">
<person-group person-group-type="author"><string-name>
<surname>Long</surname>, 
<given-names>T. Y.</given-names></string-name>, <string-name>
<surname>Cai</surname>, 
<given-names>Z. J.</given-names></string-name>
</person-group> (
<year>2019</year>). 
<source>Fluid mechanics</source>. 
<publisher-loc>Beijing</publisher-loc>: 
<publisher-name>China Architecture &#x0026; Building Press</publisher-name>.</mixed-citation>
</ref>
<ref id="ref-12">
<label>12.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Zhao</surname>, 
<given-names>H. W.</given-names></string-name>, <string-name>
<surname>Chen</surname>, 
<given-names>Y. H.</given-names></string-name>, <string-name>
<surname>Liu</surname>, 
<given-names>X. G.</given-names></string-name>
</person-group> (
<year>2019</year>). 
<article-title>Optimization of the heat dissipation structure and temperature distribution in an electric vehicle power battery</article-title>. 
<source>Fluid Dynamics &#x0026; Materials Processing</source><italic>,</italic> 
<volume>15</volume><italic>(</italic>
<issue>4</issue><italic>),</italic> 
<fpage>293</fpage>&#x2013;
<lpage>305</lpage>. DOI 
<pub-id pub-id-type="doi">10.32604/fdmp.2019.08361</pub-id>.</mixed-citation>
</ref>
<ref id="ref-13">
<label>13.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Du</surname>, 
<given-names>J.</given-names></string-name>, <string-name>
<surname>Wu</surname>, 
<given-names>X.</given-names></string-name>, <string-name>
<surname>Li</surname>, 
<given-names>R. N.</given-names></string-name>, <string-name>
<surname>Cheng</surname>, 
<given-names>R. R.</given-names></string-name>
</person-group> (
<year>2019</year>). 
<article-title>Numerical simulation and optimization of a mid-temperature heat pipe exchanger</article-title>. 
<source>Fluid Dynamics &#x0026; Materials Processing</source><italic>,</italic> 
<volume>15</volume><italic>(</italic>
<issue>4</issue><italic>),</italic> 
<fpage>293</fpage>&#x2013;
<lpage>305</lpage>. DOI 
<pub-id pub-id-type="doi">10.32604/fdmp.2019.08361</pub-id>.</mixed-citation>
</ref>
<ref id="ref-14">
<label>14.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Fang</surname>, 
<given-names>J.</given-names></string-name>, <string-name>
<surname>Wang</surname>, 
<given-names>L.</given-names></string-name>, <string-name>
<surname>Hong</surname>, 
<given-names>Y.</given-names></string-name>, <string-name>
<surname>Zhao</surname>, 
<given-names>J. D.</given-names></string-name>
</person-group> (
<year>2021</year>). 
<article-title>Influence of solid-fluid interaction on impact dynamics against rigid barrier: CFD-DEM modelling</article-title>. 
<source>G&#x00E9;otechnique</source><italic>,</italic> <fpage>1</fpage>&#x2013;<lpage>16</lpage>.</mixed-citation>
</ref>
<ref id="ref-15">
<label>15.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Alekseeenko</surname>, 
<given-names>V. V.</given-names></string-name>, <string-name>
<surname>Yashchenko</surname>, 
<given-names>V. P.</given-names></string-name>, <string-name>
<surname>Sukhanova</surname>, 
<given-names>K. M.</given-names></string-name>, <string-name>
<surname>Shestakov</surname>, 
<given-names>A. N.</given-names></string-name>
</person-group> (
<year>2021</year>). 
<article-title>Determination of elasticity modulus of flexible reinforced PVC hoses material at temperature change based on experimental studies and finite element modeling</article-title>. 
<source>IOP Conference Series: Earth and Environmental Science</source><italic>,</italic> 
<volume>751</volume><italic>(</italic>
<issue>1</issue><italic>),</italic> <fpage>012071</fpage>.</mixed-citation>
</ref>
<ref id="ref-16">
<label>16.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Guo</surname>, 
<given-names>H.</given-names></string-name>, <string-name>
<surname>Zhou</surname>, 
<given-names>Y.</given-names></string-name>, <string-name>
<surname>Li</surname>, 
<given-names>Y.</given-names></string-name>, <string-name>
<surname>Zhang</surname>, 
<given-names>Z.</given-names></string-name>, <string-name>
<surname>Li</surname>, 
<given-names>H.</given-names></string-name>
</person-group> (
<year>2021</year>). 
<article-title>An experiment-based simplified method for the model of building groups in CFD simulation</article-title>. 
<source>Advances in Civil Engineering</source><italic>,</italic> <fpage>1</fpage>&#x2013;<lpage>13</lpage>.</mixed-citation>
</ref>
<ref id="ref-17">
<label>17.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Yu</surname>, 
<given-names>D. L.</given-names></string-name>, <string-name>
<surname>Zhang</surname>, 
<given-names>H. L.</given-names></string-name>, <string-name>
<surname>Zeng</surname>, 
<given-names>X.</given-names></string-name>, <string-name>
<surname>Liao</surname>, 
<given-names>D. H.</given-names></string-name>, <string-name>
<surname>Wu</surname>, 
<given-names>N. X.</given-names></string-name>
</person-group> (
<year>2021</year>). 
<article-title>Influence of nozzle orifice shape on the atomization process of Si<sub>3</sub>N<sub>4</sub> in a dry granulation process</article-title>. 
<source>Fluid Dynamics &#x0026; Materials Processing</source><italic>,</italic> 
<volume>17</volume><italic>(</italic>
<issue>3</issue><italic>),</italic> 
<fpage>569</fpage>&#x2013;
<lpage>586</lpage>. DOI 
<pub-id pub-id-type="doi">10.32604/fdmp.2021.014711</pub-id>.</mixed-citation>
</ref>
<ref id="ref-18">
<label>18.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Tawfiq</surname>, 
<given-names>C.</given-names></string-name>
</person-group> (
<year>2019</year>). 
<article-title>Droplet breakup regime in a cross-junction device with lateral obstacles</article-title>. 
<source>Fluid Dynamics &#x0026; Materials Processing</source><italic>,</italic> 
<volume>15</volume><italic>(</italic>
<issue>5</issue><italic>),</italic> 
<fpage>545</fpage>&#x2013;
<lpage>555</lpage>. DOI 
<pub-id pub-id-type="doi">10.32604/fdmp.2019.01793</pub-id>.</mixed-citation>
</ref>
<ref id="ref-19">
<label>19.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Philipp</surname>, 
<given-names>E.</given-names></string-name>, <string-name>
<surname>Braitsch</surname>, 
<given-names>M.</given-names></string-name>, <string-name>
<surname>Bichsel</surname>, 
<given-names>T.</given-names></string-name>, <string-name>
<surname>M&#x00FC;hlebach</surname>, <given-names>S.</given-names></string-name>
</person-group> (
<year>2016</year>). 
<article-title>Diluting ferric carboxymaltose in sodium chloride infusion solution (0.9% w/v) in polypropylene bottles and bags: Effects on chemical stability</article-title>. 
<source>European Journal of Hospital Pharmacy&#x2013;Science and Practice</source><italic>,</italic> 
<volume>23</volume><italic>(</italic>
<issue>1</issue><italic>),</italic> 
<fpage>22</fpage>&#x2013;
<lpage>27</lpage>.</mixed-citation>
</ref>
<ref id="ref-20">
<label>20.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Chen</surname>, 
<given-names>P. X.</given-names></string-name>, <string-name>
<surname>Huang</surname>, 
<given-names>K.</given-names></string-name>, <string-name>
<surname>Wang</surname>, 
<given-names>F. H.</given-names></string-name>, <string-name>
<surname>Xie</surname>, 
<given-names>W. J.</given-names></string-name>, <string-name>
<surname>Wei</surname>, 
<given-names>S.</given-names></string-name> <etal>et al.</etal>
</person-group> (
<year>2019</year>). 
<article-title>Simulation of heat and mass transfer in a grain pile on the basis of a 2D irregular pore network</article-title>. 
<source>Fluid Dynamics &#x0026; Materials Processing</source><italic>,</italic> 
<volume>15</volume><italic>(</italic>
<issue>4</issue><italic>),</italic> 
<fpage>367</fpage>&#x2013;
<lpage>389</lpage>. DOI 
<pub-id pub-id-type="doi">10.32604/fdmp.2019.07762</pub-id>.</mixed-citation>
</ref>
<ref id="ref-21">
<label>21.</label><mixed-citation publication-type="book">
<person-group person-group-type="author"><string-name>
<surname>Cai</surname>, 
<given-names>G. C.</given-names></string-name>
</person-group> (
<year>2019</year>). 
<source>Modern numerical analysis</source>. 
<publisher-loc>Beijing</publisher-loc>: 
<publisher-name>Science Press</publisher-name>.</mixed-citation>
</ref>
<ref id="ref-22">
<label>22.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Carly</surname>, 
<given-names>W.</given-names></string-name>, <string-name>
<surname>Dominic</surname>, 
<given-names>F.</given-names></string-name>, <string-name>
<surname>Galal</surname>, 
<given-names>H. G.</given-names></string-name>, <string-name>
<surname>Ann</surname>, 
<given-names>B.</given-names></string-name>, <string-name>
<surname>Bryony</surname>, 
<given-names>D. F.</given-names></string-name>
</person-group> (
<year>2020</year>). 
<article-title>Patients&#x2019; perspectives on the quality and safety of intravenous infusions: A qualitative study</article-title>. 
<source>Journal of Patient Experience</source><italic>,</italic> 
<volume>7</volume><italic>(</italic>
<issue>3</issue><italic>),</italic> <fpage>380</fpage>&#x2013;<lpage>385</lpage>.</mixed-citation>
</ref>
<ref id="ref-23">
<label>23.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Zeng</surname>, 
<given-names>F. Q. C.</given-names></string-name>
</person-group> (
<year>2019</year>). 
<article-title>Discussion on the clinical effect of nursing safety management on reducing the risk of intravenous infusion</article-title>. 
<source>The Medical Forum</source><italic>,</italic> 
<volume>23</volume><italic>(</italic>
<issue>9</issue><italic>),</italic> 
<fpage>1301</fpage>&#x2013;
<lpage>1302</lpage>.</mixed-citation>
</ref>
<ref id="ref-24">
<label>24.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Le</surname>, 
<given-names>G. G.</given-names></string-name>, <string-name>
<surname>Vieillard</surname>, 
<given-names>V.</given-names></string-name>, <string-name>
<surname>Mouraud</surname>, 
<given-names>S.</given-names></string-name>, <string-name>
<surname>Do</surname>, 
<given-names>B.</given-names></string-name>, <string-name>
<surname>Marabelle</surname>, 
<given-names>A.</given-names></string-name> <etal>et al.</etal>
</person-group> (
<year>2020</year>). 
<article-title>Stability of nivolumab in its original vials after opening and handing in normal saline bag for intravenous infusion</article-title>. 
<source>European Journal of Cancer</source><italic>,</italic> 
<volume>135</volume><italic>,</italic> <fpage>192</fpage>&#x2013;<lpage>202</lpage>.</mixed-citation>
</ref>
<ref id="ref-25">
<label>25.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Chambers</surname>, 
<given-names>D. J.</given-names></string-name>
</person-group> (
<year>2018</year>). 
<article-title>Principles of intravenous drug infusion</article-title>. 
<source>Anaesthesia &#x0026; Intensive Care Medicine</source><italic>,</italic> <volume>20</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>61</fpage>&#x2013;<lpage>64</lpage>.</mixed-citation>
</ref>
<ref id="ref-26">
<label>26.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Liu</surname>, 
<given-names>X. H.</given-names></string-name>, <string-name>
<surname>Cheng</surname>, 
<given-names>Z. P.</given-names></string-name>, <string-name>
<surname>Chen</surname>, 
<given-names>D. G.</given-names></string-name>
</person-group> (
<year>2013</year>). 
<article-title>The influence of wall thickness of vertical polypropylene infusion bag on infusion performance</article-title>. 
<source>Bioengineering and Clinical</source><italic>,</italic> 
<volume>9</volume><italic>(</italic>
<issue>2</issue><italic>),</italic> 
<fpage>3</fpage>&#x2013;
<lpage>6</lpage>.</mixed-citation>
</ref>
<ref id="ref-27">
<label>27.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Maxime</surname>, 
<given-names>P.</given-names></string-name>, <string-name>
<surname>Bertrand</surname>, 
<given-names>D.</given-names></string-name>, <string-name>
<surname>Wadih</surname>, 
<given-names>A. C.</given-names></string-name>, <string-name>
<surname>Brigitte</surname>, 
<given-names>N.</given-names></string-name>, <string-name>
<surname>Laurent</surname>, 
<given-names>S.</given-names></string-name> <etal>et al.</etal>
</person-group> (
<year>2018</year>). 
<article-title>Effectiveness of in-line filters to completely remove particulate contamination during a pediatric multidrug infusion protocol</article-title>. 
<source>Scientific Reports</source><italic>,</italic> 
<volume>8</volume><italic>(</italic>
<issue>1</issue><italic>),</italic> <fpage>1</fpage>&#x2013;<lpage>8</lpage>.</mixed-citation>
</ref>
<ref id="ref-28">
<label>28.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Marzie</surname>, 
<given-names>F.</given-names></string-name>, <string-name>
<surname>Abdolreza</surname>, 
<given-names>M.</given-names></string-name>, <string-name>
<surname>Vahid</surname>, 
<given-names>V.</given-names></string-name>
</person-group> (
<year>2021</year>). 
<article-title>Influence of solvent type on the ability and properties of the boehmite nanoparticles embedded emulsion polyvinyl chloride nanocomposite ultrafiltration membranes</article-title>. 
<source>Iranian Polymer Journal</source><italic>,</italic> <fpage>1</fpage>&#x2013;<lpage>15</lpage>.</mixed-citation>
</ref>
<ref id="ref-29">
<label>29.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Jain</surname>, 
<given-names>D.</given-names></string-name>, <string-name>
<surname>Bhagat</surname>, 
<given-names>H.</given-names></string-name>, <string-name>
<surname>Jain</surname>, 
<given-names>D.</given-names></string-name>
</person-group> (
<year>2020</year>). 
<article-title>Effect of intravenous lignocaine infusion on the quality of emergence in patients undergoing transsphenoidal resection of pituitary tumors&#x2014;A prospective, randomized controlled trial</article-title>. 
<source>Surgical Neurology International</source><italic>,</italic> 
<volume>11</volume>.</mixed-citation>
</ref>
<ref id="ref-30">
<label>30.</label><mixed-citation publication-type="journal">
<person-group person-group-type="author"><string-name>
<surname>Lv</surname>, 
<given-names>J.</given-names></string-name>, <string-name>
<surname>Zhang</surname>, 
<given-names>M. L.</given-names></string-name>, <string-name>
<surname>Liu</surname>, 
<given-names>H. Y.</given-names></string-name>, <string-name>
<surname>Wang</surname>, 
<given-names>Y.</given-names></string-name>, <string-name>
<surname>Ma</surname>, 
<given-names>Y. F.</given-names></string-name>
</person-group> (
<year>2018</year>). 
<article-title>Study on the influence of reserved air volume in Non-PVC double valve soft bag on the quality of intravenous infusion</article-title>. 
<source>Journal of Nurses Training</source><italic>,</italic> 
<volume>33</volume><italic>(</italic>
<issue>17</issue><italic>),</italic> 
<fpage>1614</fpage>&#x2013;
<lpage>1616</lpage>.</mixed-citation>
</ref>
</ref-list>
</back>
</article>