Mathematical Modeling of Heat and Mass Transfer in Energy Science and Engineering 2014
View this Special IssueResearch Article  Open Access
Kewei Song, Wenkai Li, Yang Zhou, Yuanru Lu, "Numerical Study of Buoyancy Convection of Air under Permanent Magnetic Field and Comparison with That under Gravity Field", Mathematical Problems in Engineering, vol. 2014, Article ID 494585, 13 pages, 2014. https://doi.org/10.1155/2014/494585
Numerical Study of Buoyancy Convection of Air under Permanent Magnetic Field and Comparison with That under Gravity Field
Abstract
Magnetothermal free convection of air in a square enclosure under a nonuniform magnetic field provided by a permanent neodymiumironboron magnet is numerically studied. The natural convection under the gravity field alone is also studied for comparison. The physical fields of magnetizing force, velocity, and temperature as well as the local distribution characteristic of Nusselt number are all presented in this paper. The results show that the buoyancy convection of air in the square enclosure under magnetic field is quite different from that under the gravity field. The local value of Nusselt number under the magnetic field supplied by a permanent magnet with a residual magnetic flux density of about 4.5 Tesla can reach a high value of about three times larger than the maximum local value of Nusselt number under the gravity field. Relatively uniform distributions of temperature gradient and Nusselt number can be obtained along the cold wall of the enclosure under the magnetic field. A permanent magnet with high magnetic energy product with reaching to 3.5 Tesla can play a comparative role on the averaged Nusselt number compared with that under the gravity environment.
1. Introduction
The heat transfer and flow characteristics of natural convection in enclosures have attracted much research over the years due to their many practical engineering applications, such as in building insulation, growing crystals, solar energy collection, cooling of electrical industries, and flows in rooms due to thermal energy sources [1, 2]. Natural convection flows in a vertical cavity with two vertical walls at different temperatures and with adiabatic horizontal surfaces are the most considered configuration in the studies of natural convection because of their relative simplicity and practical importance [3–5].
Fast development of superconducting materials at high temperature has enabled us to utilize commercial superconducting magnets that produce a magnetic flux density up to 10 Tesla or more. With the advent of superconducting magnets, magnetic convection can be studied on Earth, and various new findings have been reported in the last decades [6–10]. Wakayama [8–10] has been active in finding new and notable effects of a strong magnetic field in fluid convection. In the fields of fluid mechanics and heat transfer, the application of such high magnetic field on the control of natural convection has received considerable attention [11–21]. The heat transfer control of natural convection is based on the principle that the magnetizing force acts at a molecular level depending on the values of magnetic susceptibility of fluid. Therefore, the utilization of the magnetizing force is equivalent to the local gravity control within the region where the high gradient magnetic field exists. Braithwaite et al. [11] reported enhancement or cancellation of gravitational convection due to a magnetic field for a solution of gadolinium nitrate in a shadow layer heated from below and cooled from above. Tagawa et al. [12–20] carried out studies about the effects of magnetic fields on the natural convection. These studies have shown experimentally and numerically how to enhance or suppress heat transfer by the application of a magnetic field supplied by superconducting magnet.
Most of the previous studies about the magnetothermal natural convection in a vertical cavity focus on the effect of magnetic field supplied by superconducting magnet together with gravity field on natural convection. Bednarz et al. [21] studied the transient magnetic convection boundary layer in a microgravity environment. But the magnetic field supplied by the superconducting magnet is often assumed to be uniform; that is, the magnetic field gradient is constant. In this case, the problem under the uniform magnetic field is similar to that under gravity field. As the magnetic field supplied by permanent magnet is nonuniform and the gravity field is uniform, the natural convection characteristics under such two fields must be quite different. Recently, along with the development of the technology and improvement of manufacturing process, the residual magnetic flux density of a permanent magnet increases quickly. The maximum value of the residual magnetic flux density of a neodymiumironboron magnet can reache to 1.5 Tesla or more. It is believed that the residual magnetic flux density of a permanent magnet will reach a larger value relying on the fast development of technology. The permanent magnet with large residual magnetic flux density can supply a nonuniform magnetic field and drive the gas free convection easily under a gradient temperature field. Thus, the permanent magnet which has high magnetic energy product can be widely applied in future for control of natural convection without any power supply from outside like a superconducting magnet.
The present study focuses on the heat transfer and flowing characteristic of the magnetothermal convection of air in a square enclosure driven by a permanent magnet field alone. Comparison is also carried out for the convection driven by magnetic field and gravity field, respectively.
2. Physical Model
The physical model, twodimensional square enclosure, considered in this paper is shown in Figure 1. A cubic permanent magnet with the residual magnetic flux density in the direction is put on the bottom of the square enclosure. The length of both the magnet and square enclosure is 0.05 m. There is a gap of 1 mm between the magnet and the square enclosure. Figures 1(a) and 1(b) show the physical model in which the left wall is heated and right wall is cooled isothermally, and the other two walls are thermally insulated. Figure 1(a) shows the physical model under gravity field; Figure 1(b) shows the model under magnetic field. The direction of magnetizing force and the gravity acceleration are also shown in Figure 1.
(a)
(b)
3. Mathematical Equations
In a single phase with isothermal state, both the magnetizing force and gravitational force are the conservative force and hence convection does not occur by itself no matter how strong it is. In order to arouse the convection, both the magnetic field gradient and the temperature gradient are necessary. In a normal expression, the magnetizing force is written by the equation as follows [22]:
The momentum equation including the magnetizing force alone is as follows:
And the momentum equation including the gravity buoyancy force alone is
We assume that the air is an incompressible Newtonian fluid and the Boussinesq approximation is employed. When we consider the isothermal state, the fluid density and susceptibility would be constant, and convection does not arise in both the gravity and magnetic fields. Parameters under this state are , , and at reference temperature . Hence we get
When there is a temperature difference, the magnetic susceptibility and density also change with temperature. Subtracting (4) from (2) and (3) gives
The density can be expressed by a Taylor expansion (keeping two terms only) around a reference state as follows:
On the other hand, the thermal volume coefficient of expansion is defined for fluid as follows [23]:
Then, (6) can be written as follows:
The difference of can be written as
Magnetic susceptibility of air is a function of temperature; this follows the Curie law:
Here, is a constant; then
The volumetric coefficient of expansion for an ideal gas can be given from the ideal gas law:
Then, from (7) so that
Given , the momentum equation (5) can be written as follows:
Thus, the government equations for the natural convection caused by the magnetic buoyancy force or gravity buoyancy force are summarized as follows.
Continuity equation is
Momentum equation, for magnetic field alone is
And for gravity field alone
Energy equation is
In all the subsequent computation, the reference temperature is selected equal to the temperature of the cold wall .
The initial condition is
Boundary conditions are
The Rayleigh number is
The local heat transfer coefficient is defined as
The local Nusselt number Nu_{local} represents the ratio of heat transfer rate by convection to that by conduction in the fluid. Nu_{local} is given by
The overall averaged Nusselt number is obtained by averaging Nu_{local} on the wall involved in heat transfer:
In this study, the magnetic field is assumed to be provided by a cubic permanent neodymiumironboron (NdFeB) magnet which has a high magnetic energy product and could drive the gas free convection easily under a gradient temperature field. The magnetic field is solved using the vector potential method [24].
According to Maxwell’s equations, the divergences of is always zero:
The fact that the divergence of is zero means that there cannot exist sources of . Because the divergence of the curl of a vector is identically equal to zero. Thus the magnetic induction can be expressed as the curl of a vector .
is the vector potential. The vector potential generated by the permanent magnet at the calculated point is as follows: Here, is the volume of the magnet, is the distance from the element volume of the magnet to the point where the field is calculated, and is the unit vector along . When the calculated field point is fixed and the source point, that is the center of element volume , is moving, the gradient of taken at the source point is Equation (28) can be rewritten as According to the curl calculus identities, the curl of the product of a scalar function and a vector function is And one of the properties of vector product is that the vector product is anticommutative: Then, Equation (30) can be rewritten as The curl calculus identity has where is the surface bounding the volume ; is element area of the surface , is the direction normal to the surface . The magnetic vector potential can be rewritten as is the equivalent surface current density and is the equivalent volume current density; the unit vector points, from the source, which is the loop, to the point where the magnetic field is calculated. Here using the BiotSavart law we can write down the value of : Thus, if we know the permanent magnet’s magnetization , we can find the equivalent current densities. Then we can use the equivalent currents to calculate by the upper equation.
In this study, the permanent magnet is applied near the bottom wall, and the poles of the magnet are lined parallel to the direction. The square enclosure is placed on top of the cubic magnet in the middle of the side and parallel to the  plane. The computed magnetic field lines supplied by the cubic permanent magnet in the enclosure are shown in Figure 2.
The SIMPLE algorithm developed by Patankar [25] is used to solve the coupled heat transfer and fluid flow problem. The powerlaw scheme is used in the finite difference formulation of convection terms and a fine grid system is selected to raise the simulation precision. The gridindependence test for the solutions is carried out with three grid systems under the gravity field only. The averaged Nusselt numbers at different grid systems are presented in Table 1. Since the difference of Nusselt number between three grid systems is less than 0.5%, the grid system 82 × 82 is selected for all computations. The physical properties of air and other constants used in this paper are summarized in Table 2.


4. Results and Discussion
All the computations are carried out from to , and the residual flux densities of permanent magnet changes from 0.5 T to 4.5 T. The flow field, temperature field, and the distribution of local Nusselt number along the walls are obtained under the magnetic field with ; the difference of the effect of magnetic field and gravity field on natural convection is also studied.
4.1. Temperature Fields
Figure 3 shows the temperature fields of the studied physical models at and 10^{5}. For the physical model under the gravity field, as shown in Figure 3(a), the temperature field looks rotationally symmetric about the center of the enclosure. The temperature gradient normal to the hot wall decreases and the temperature gradient near the cold wall increases gently from the bottom wall to the top wall. The largest temperature gradient near the hot wall locates at the position close to the bottom wall, but near the cold wall it is close to the top wall. The temperature field structure under the magnetic field, as shown in Figure 3(b), is quite different from the temperature field under the gravity field. The temperature gradient near the hot and cold walls decreases along direction. The maximum gradient of the temperature near the hot wall locates at the bottom of the hot wall, and the maximum temperature gradient on the cold wall locates at the position close to the bottom wall. The temperature gradient normal to the cold wall changes a little from top wall to the position close to the bottom wall, and then it decreases quickly till to the bottom wall in a small range.
(a)
(b)
4.2. Magnetizing Force Fields
Figure 4 shows the gravity buoyancy force and magnetizing buoyancy force fields in the enclosure. From these figures we can find that the magnetizing buoyancy force is different from the gravity buoyancy force. The gravity buoyancy force parallels along direction, but the magnetizing buoyancy force’s direction changes in the enclosure. The magnetizing buoyancy force increased with increase of magnet strength. There is large magnetizing force near the bottom wall of the enclosure and the force near the hot wall is larger than that near the cold wall. The magnetizing force decreases quickly along direction. But for the gravity field, the distribution of gravity buoyancy force is quite different. The larger gravity buoyancy force locates in the region near the top and hot wall, and the gravity buoyancy force decreases from hot wall to cold wall and increases from bottom wall to top wall. The maximum value of magnetizing buoyancy force is larger than the maximum gravity buoyancy force.
(a)
(b)
4.3. Velocity Fields
Figure 5 shows the velocity fields under gravity field and magnetic field. When there is gravity field only, as shown in Figure 5(a), the large value of velocity locates near the walls and the value of velocity decreases from the walls to the center of the enclosure. There is one vortex in the center of the enclosure at . When increases to 10^{5}, the core of the vortex is divided into two cores, one is located near the top and hot walls, and the other one locates close to the cold and bottom walls. When there is magnetic field only, the velocity field is different from that under the gravity field. A vortex is also formed in the enclosure, but the core of the vortex locates close to the bottom wall at . When increases to 10^{5}, there is still only one core of the vortex and the core moves much closer to the hot wall. Figures 6(a) and 6(b) show the streamlines in the enclosure under gravity field and magnetic field, respectively. The streamlines related to gravity field look rotationally symmetric about the center of the enclosure; large velocity gradient normal to the wall locates near the bottom segment of the hot wall and the top segment of the cold wall. But under the magnetic field, the streamlines are not symmetric and there are dense streamlines near the corner of the hot and bottom walls where there is large value of magnetizing force as shown in Figure 4(b). In the region near the top wall, the velocity gradient normal to the wall is significantly less than that in other regions.
(a)
(b)
(a)
(b)
4.4. Distribution of Temperature along Lines
The distribution curves of temperature at along the lines parallel to axis are plotted in Figure 7; these three lines , , and are at the positions , 0.5, and 0.8, respectively. In the region near the hot wall, the temperature under the magnetic field along the line is lower than that under the gravity field, but along the lines and , the temperature under the magnetic field is higher than the temperature under the gravity field. In the region near the cold wall, the temperature related to magnetic field is higher than the temperature related to gravity field along the lines and , and the temperature related to magnetic field is lower than the temperature related to the gravity field along line . The high temperature near the cold wall means a high temperature gradient and high temperature near the hot wall means a low temperature gradient. For the magnetic field, the largest temperature gradient near the cold wall is along the line and the largest temperature gradient near the hot wall is along the line . We can obtain the same conclusion as in Figure 7 that the distribution of the temperature under the magnetic field is quite different from that under the gravity field.
4.5. Heat Transfer on the Walls
The distribution curves of the local Nusselt number along the hot and cold walls are plotted in Figure 8. The distribution of the local Nusselt number under gravity field is shown in Figure 8(a). The peak value along the hot wall locates near the bottom wall and decreases nearly linearly to the top wall. The distribution of local Nusselt number along the cold wall is symmetric about the central line of the enclosure with the distribution of local Nusselt number along the hot wall. The value of the Nusselt number increases nearly linearly from bottom wall and gets the peak value near the top wall; then it decreases again till to the top wall. Figure 8(b) shows the distribution curves of local Nusselt number under the magnetic field with . We can see that the distribution of the local Nusselt number for the model under magnetic field is different from that shown in Figure 8(a). The local Nusselt number gets a high peak value at the bottom of the hot wall and it decreases quickly for a short distance and then continues decreasing at a relatively slow rate till it reaches the top wall. The change of the local Nusselt number along the cold wall is much gentler compared with that along the hot wall. The local Nusselt number increases slowly from the top wall to the position near the bottom wall and then it decreases quickly till it reaches the bottom wall.
(a)
(b)
The comparison of the distribution curves of the local Nusselt number under gravity field and magnetic field are presented in Figure 9. The value of local Nusselt number along the hot wall is larger than the value under the gravity filed in a short line segment from the bottom wall; the peak value of local Nusselt number related to magnetic field is more than three times larger than the peak value related to the gravity field. On the cold wall, the local Nusselt number related to magnetic field changes much more gently and has larger value on the lower part of the wall compared with that related to gravity field.
(a)
(b)
Figure 10 shows the distribution of the average Nusselt number with ranging from 10^{4} to 10^{6} and the residual magnetic flux density of a permanent magnet ranging from 1.5 T to 4.5 T. The average under the magnetic field has a similar distribution curve compared with that under the gravity field. The averaged increases with increase of . The averaged Nusselt numbers on the hot wall with are summarized in Table 3. When is about 3.5 T, the averaged is nearly the same as that under gravity field; this means that the magnetizing buoyancy force has a comparative effect on in the enclosure compared with the gravity buoyancy force.

5. Conclusion
Natural convection of air in a twodimensional square enclosure under a nonuniform magnetic field provided by a permanent magnet is carried out. The results are compared with that under the gravity field. As the magnetizing buoyancy force of air in the magnetic field is quite different from the gravitational buoyancy force, the natural convection in the enclosure under a magnetic field provided by a permanent magnet has quite different characteristics. For the model studied in this paper, the velocity and temperature fields are rotationally symmetric about the center of the enclosure. The heat transfer Nusselt numbers on the hot wall and the cold wall are symmetrical about the center line of the enclosure. But the velocity and temperature fields under the magnetic field do not have such characteristics. Heat transfer Nusselt number is different on the hot and cold walls. The value of local Nusselt number has a peak value on the bottom of the hot wall and decreases quickly along the hot wall, but the local Nusselt number changes with a relatively gentle rate on the cold wall. The peak value on the hot wall is three times larger than that under the gravity field when and 2 times larger when . When is about 3.5 T, the value of averaged is nearly the same with that under gravity field. This means that in the zero gravity environment, a permanent magnet with can play a comparative role on the averaged Nusselt number compared with that under the gravity environment.
Nomenclature
:  Vector potential 
:  Magnetic flux density (T) 
:  Residual magnetic flux density of a permanent magnet (T) 
:  Magnetizing force (N/m^{3}) 
:  Acceleration coefficient of gravity (m/s^{2}) 
:  Magnetic field strength (A/m) 
:  Magnetization 
:  Nusselt number 
:  Pressure (N/m^{2}) 
:  Reference pressure without convection of gas (N/m^{2}) 
:  Perturbed pressure due to convection 
:  Prandtl number 
:  Unit vector 
:  Rayleigh number 
:  Time (s) 
:  Temperature (K) 
:  Velocity vector (u, w)(m/s) 
:  Coordinate (m) 
:  Coordinate (m). 
:  Thermal diffusivity of gas (m^{2}/s) 
:  Volume coefficient of expansion of gas due to temperature difference (1/K) 
:  Viscosity of gas (Pa·s) 
:  Magnetic permeability of free space (H/m) 
:  Kinematic viscosity of gas (m^{2}/s) 
:  Density of gas (kg/m^{3}) 
:  Density of gas at a reference state of no convection (kg/m^{3}) 
:  Mass magnetic susceptibility of air (m^{3}/kg) 
:  Conduction coefficient (W/(m·K)). 
0:  Reference state 
:  That of a cold wall 
:  That of a hot wall. 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work is supported by Gansu Provincial Funds for Distinguished Young Scientists (no. 145RJDA324), the National Natural Science Foundation of China (no. 51376086), and the College Students’ Innovative Experiment Project of Lanzhou Jiaotong University (20140048, 20140050).
References
 B. Gebhart, Y. Yaluria, R. L. Mahajan, and B. Sammakia, Buoyancy Induced Flows and Transport, Hemisphere, Washington, DC, USA, 1988.
 K. T. Yang, “Natural convection in enclosures,” in Handbook of Single Phase Convection Heat Transfer, Wiley, New York, NY, USA, 1987. View at: Google Scholar
 G. de Vahl Davis, “Natural convection of air in a square cavity: a benchmark numerical solution,” International Journal for Numerical Methods in Fluids, vol. 3, no. 3, pp. 249–264, 1983. View at: Publisher Site  Google Scholar
 C. J. Freitas, R. L. Street, A. N. Findikakis, and J. R. Koseff, “Numerical simulation of three dimensional flow in a cavity,” International Journal for Numerical Methods in Fluids, vol. 5, no. 6, pp. 561–575, 1985. View at: Publisher Site  Google Scholar
 T. Fusegi, J. M. Hyun, K. Kuwahara, and B. Farouk, “A numerical study of threedimensional natural convection in a differentially heated cubical enclosure,” International Journal of Heat and Mass Transfer, vol. 34, no. 6, pp. 1543–1557, 1991. View at: Publisher Site  Google Scholar
 D. Braithwaite, E. Beaugnon, and R. Tournier, “Magnetically controlled convection in a paramagnetic fluid,” Nature, vol. 354, no. 6349, pp. 134–136, 1991. View at: Publisher Site  Google Scholar
 N. I. Wakayama, “Behavior of gas flow under gradient magnetic fields,” Journal of Applied Physics, vol. 69, no. 4, pp. 2734–2736, 1991. View at: Publisher Site  Google Scholar
 N. I. Wakayama, “Magnetic promotion of combustion in diffusion flames,” Combustion and Flame, vol. 93, no. 3, pp. 207–214, 1993. View at: Publisher Site  Google Scholar
 N. I. Wakayama, “Magnetic acceleration and deceleration of O_{2} gas streams injected into air,” IEEE Transactions on Magnetics, vol. 31, no. 1, pp. 897–901, 1995. View at: Publisher Site  Google Scholar
 N. I. Wakayama and M. Wakayama, “Magnetic acceleration of inhaled and exhaled flows in breathing,” Japanese Journal of Applied Physics, vol. 39, no. 3, pp. L262–L264, 2000. View at: Publisher Site  Google Scholar
 O. Braithwaite, E. Beaugnon, and R. Tournier, “Magnetically controlled convection in a paramagnetic fluid,” Nature, vol. 354, no. 6349, pp. 134–136, 1991. View at: Publisher Site  Google Scholar
 T. Tagawa, H. Ozoe, K. Inoue, M. Ito, K. Sassa, and S. Asai, “Transient characteristics of convection and diffusion of oxygen gas in an open vertical cylinder under magnetizing and gravitational forces,” Chemical Engineering Science, vol. 56, no. 14, pp. 4217–4223, 2001. View at: Publisher Site  Google Scholar
 T. Tagawa, R. Shigemitsu, and H. Ozoe, “Magnetizing force modeled and numerically solved for natural convection of air in a cubic enclosure: effect of the direction of the magnetic field,” International Journal of Heat and Mass Transfer, vol. 45, no. 2, pp. 267–277, 2002. View at: Publisher Site  Google Scholar
 M. Akamatsu, M. Higano, Y. Takahashi, and H. Ozoe, “Numerical computation on the control of aerial flow by the magnetizing force in gravitational and nongravitational fields,” Numerical Heat Transfer A: Applications, vol. 43, no. 1, pp. 9–29, 2003. View at: Publisher Site  Google Scholar
 T. Bednarz, T. Tagawa, M. Kaneda, H. Ozoe, and J. S. Szmyd, “Convection of air in a cubic enclosure with an electric coil inclined in general orientations,” Fluid Dynamics Research, vol. 36, no. 2, pp. 91–106, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 P. Filar, E. Fornalik, M. Kaneda, T. Tagawa, H. Ozoe, and J. S. Szmyd, “Threedimensional numerical computation for magnetic convection of air inside a cylinder heated and cooled isothermally from a side wall,” International Journal of Heat and Mass Transfer, vol. 48, no. 9, pp. 1858–1867, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 T. Bednarz, E. Fornalik, T. Tagawa, H. Ozoe, and J. S. Szmyd, “Experimental and numerical analyses of magnetic convection of paramagnetic fluid in a cube heated and cooled from opposing verticals walls,” International Journal of Thermal Sciences, vol. 44, no. 10, pp. 933–943, 2005. View at: Publisher Site  Google Scholar
 T. P. Bednarz, C. Lei, J. C. Patterson, and H. Ozoe, “Effects of a transverse, horizontal magnetic field on natural convection of a paramagnetic fluid in a cube,” International Journal of Thermal Sciences, vol. 48, no. 1, pp. 26–33, 2009. View at: Publisher Site  Google Scholar
 T. Bednarz, E. Fornalik, H. Ozoe, J. S. Szmyd, J. C. Patterson, and C. Lei, “Influence of a horizontal magnetic field on the natural convection of paramagnetic fluid in a cube heated and cooled from two vertical side walls,” International Journal of Thermal Sciences, vol. 47, no. 6, pp. 668–679, 2008. View at: Publisher Site  Google Scholar
 T. Bednarz, T. Tagawa, M. Kaneda, H. Ozoe, and J. S. Szmyd, “Magnetic and gravitational convection of air with a coil inclined around the x axis,” Numerical Heat Transfer A: Applications, vol. 46, no. 1, pp. 99–113, 2004. View at: Publisher Site  Google Scholar
 T. P. Bednarz, W. Lin, J. C. Patterson, C. Lei, and S. W. Armfield, “Scaling for unsteady thermomagnetic convection boundary layer of paramagnetic fluids of $Pr>\text{1}$ in microgravity conditions,” International Journal of Heat and Fluid Flow, vol. 30, no. 6, pp. 1157–1170, 2009. View at: Publisher Site  Google Scholar
 B. Bai, A. Yabe, J. Qi, and N. I. Wakayama, “Quantitative analysis of air convection caused by magneticfluid coupling,” AIAA Journal, vol. 37, no. 12, pp. 1538–1543, 1999. View at: Publisher Site  Google Scholar
 R. B. Bird, W. E. Stewart, and E. N. Lightfoot, Transport Phenomena, John Wiley & Sons, New York, NY, USA, 2007.
 P. Lorrain and D. R. Corson, Electromagnetic Fields and Waves, W.H. Freeman and Company, San Francisco, Calif , USA, 1970.
 S. V. Patankar, Numerical Heat Transfer and Fluid Flow, Hemisphere, New York, NY, USA, 1980.
Copyright
Copyright © 2014 Kewei Song et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.