Next Article in Journal
Natural and Anthropogenic Controls of Groundwater Quality in Sri Lanka: Implications for Chronic Kidney Disease of Unknown Etiology (CKDu)
Previous Article in Journal
Impact of Water Rights Transaction in the Beijing-Tianjin-Hebei Region in China Based on an Improved Computable General Equilibrium Model
Previous Article in Special Issue
Numerical Analysis on the Heat Transfer Characteristics of Supercritical Water in Vertically Upward Internally Ribbed Tubes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Modeling of Microfluid Dynamics in Xylem Vessels of Khaya grandifoliola

by
Daniel Somma de Araujo
1,
Diogo Henrique Morato de Moraes
2,
Marcio Mesquita
2,*,
Rilner Alves Flores
2,
Rafael Battisti
2,
Glenio Guimarães Santos
2,
Fabio Ponciano de Deus
3 and
Rhuanito Soranz Ferrarezi
4
1
School of Agronomy, Federal University of Goiás, Goiânia 74690-900, GO, Brazil
2
Postgraduate Program in Agronomy, Federal University of Goiás, Goiânia 74690-900, GO, Brazil
3
Water Resources and Sanitation Department, Federal University of Lavras, Lavras 37200-900, MG, Brazil
4
Department of Horticulture, University of Georgia, Athens, GA 30602, USA
*
Author to whom correspondence should be addressed.
Water 2021, 13(19), 2723; https://doi.org/10.3390/w13192723
Submission received: 18 August 2021 / Revised: 24 September 2021 / Accepted: 27 September 2021 / Published: 1 October 2021

Abstract

:
Computational fluid dynamic (CFD) can be used to quantify the internal flow variables of xylem conducting vessels. This study aims to analyze through numerical simulations the xylem water ascent of African mahogany (Khaya grandifoliola) cultivated under different irrigation regimes. We determined a geometric model, defined through the variability of the anatomical structures of the species, observing characteristics of the xylem vessels such as diameter, length, number of pits, and average surface area of the pits. Then we applied numerical simulation through an Eulerian mathematical model with the discretization of volumes via CFD. Compared to other models, we observed that numerical simulation using CFD represented the xylem microstructures in a greater level of detail, contributing to the understanding of the flow of xylem vessels and the interference of its various structures. Analyzing the micrographs, we observed the non-irrigated vessels had a higher number of pits in the secondary wall thickening when compared to the irrigated treatments. This trend influenced the variability of the radial flow of the xylem vessels, causing greater fluid movement in this region and decreasing the influence of the smooth part of the wall, resulting in a lower total resistance of these vessels.

1. Introduction

The water transport and flow characteristics of plant xylem are explored widely through botanical experimentation [1], mainly using two-dimensional mathematical models of the xylem structure [2,3]. Although these analyses help us to understand the partitioning of xylem resistance, they ignore the complex nature of these vessel structures (secondary wall thickening, membranes, and perforations) and the effects of the three-dimensionality of arrangements in conduits and water transport [4]. To better assimilate the interference of these anatomical mechanisms in the microfluidic effects of plant vascular flow, computational fluid dynamics (CFD) can be used to provide a greater level of detail of the phenomena involved, in addition to enabling the simulation of natural conditions and, close to ideal, by using a high level of detail for the physical variables through the governing equations and their flow models [5].
The anatomical structures that make up the xylem vessels play fundamental roles in the resistance of water transport in plants [6,7,8,9,10]. Chen et al. [9], using the numerical model of turbulence for shear stress transport and k–ε to simulate the annular and helical thickening of vessels, verified that the width and inclination angle of the secondary wall thickening had a significant influence on the flow resistance. Schulte et al. [6] analyzed the dynamic effect of flow in the torus region on the conifer pit membrane. They determined that the membrane represented 6% and 18% of the flow resistance for the models exposed to the sun and shade, respectively. Xu et al. [10] found that the resistance of smooth vessels was higher in relation to other structural elements of the xylem, accounting for 66.2% of the total resistance, where the approach used was based on anatomical experiments and numerical simulations.
Characterizing these micro interferences in the vessel water ascent is fundamental, as abiotic factors can influence the vessels during the growth cycle. For example, the water regime influences the formation and development of xylem tissues, phloem, and cambium cell division [11,12,13]. In addition, water deficit can reduce the diameter of the xylem vessels, maintaining their density, hydraulic conductance, and conductive area, and modify the pattern of the secondary wall thickening (number, spatial distribution, and positioning of anatomical structures), influencing the resistances involved in water transport [12,14,15]. The understanding of the xylem’s anatomical variations and its potential interference in water ascent allows researchers to increase water use efficiency, as plants with well-developed vessels that are not collapsed by embolism/cavitation positively influence the transport of water and nutrients. This knowledge enables the study of plant tolerance of climatic stresses and the application of mathematical modeling in water ascent studies [9,14,15].
This study aims to analyze through numerical simulations the water ascent in the xylem vessels of African mahogany (Khaya grandifoliola) cultivated under different irrigation regimes, to define the variability of the anatomical structures of the species.

2. Materials and Methods

2.1. Description of Plant Material and Analysis of Anatomical Samples

Khaya grandifoliola plants were collected from a commercial field in the municipality of Engenheiro Navarro, Minas Gerais, Brazil (17°25′03.5′′ S, 44°04′22.4′′ W, 700 m). The climate of the region is classified by Köppen as Aw [13]. The average annual rainfall is between 950 and 1050 mm, with four to six months of dry season. The mean temperatures in summer range from 24.1 to 25.0 °C and in winter from 21.1 to 22.0 °C [13].
The trees were planted in a 6 m × 5 m plot with spacing in July 2008. Six Khaya grandifoliola plants were randomly selected, with three grown under irrigation and three rainfed. In the irrigated regime, 10 L of water daily per plant was added until the fifth year of planting; in subsequent years, the volume was increased to 50 L per day. In the rainfed plants, the irrigation was used during transplanting to increase the survival rate of seedlings in the field. After this period, there was no water supplementation, with watering dependent on natural rainfall.
Five disc-shaped cross-sections were removed from each plant in different longitudinal locations along the stem at 0%, 25%, 50%, 75%, and 100% of the total height (Figure 1a). From these sections, specimens with dimensions of 1.5 cm (tangential section) × 2.0 cm (radial section) × 3.0 cm (cross-section) were removed [11]. Subsequently, we polished the specimens with sandpapers with a 60-, 120-, 240-, 360-, and 600-grain variation. Samples were submitted to macroscopic analysis under a 25-fold increase to characterize the average diameter of the xylem vessels using the software Image-Pro Plus (Media Cybernetics, Silver Spring, MD, USA).
To measure the length of the xylem vessels, type, and area of interference of the secondary wall, small fragments of the heartwood were removed from the specimens in the radial position of 50%, for maceration, based on the methodology proposed by Trevisan et al. [17] (Figure 1c,d). The samples were then washed in distilled water, stained with 50% (p/v) hydroalcoholic Safranina dye, and kept in a flask until provisional slides were prepared. For a more detailed analysis, we used a scanning electron microscope (SEM) with energy-dispersive X-ray analysis (Jeol JSM-IT300LV; JEOL USA Inc, Peabody, MA, USA) [18]. With these analyses, we determined the vessels’ length, diameter, number of interferences and surface area.

2.2. Determination of the Anatomical Model and Numerical Simulation Parameterization

Three-dimensional calculations and the numerical mesh were elaborated using the software ICEM-CFD (ANSYS Inc., Canonsburg, PA, USA). Unstructured tetrahedrons and hexahedrons were used to discretize the vessel geometry. The control volumes near the walls have been refined enough to keep them within the logarithmic boundary layer. We performed mesh independence tests to ensure that the computed results were independent of refinement. The flux depends on the steady-state conservation equations for mass and momentum in a fluid, which are provided by Equations (1) and (2).
Continuity equation:
u x + v y + w z = 0 ,
Momentum equation:
{ ρ ( u u x + v u y + w u z ) = P x + μ ( 2 u x 2 + 2 u y 2 + 2 u z 2 ) ρ ( u v x + v v y + w v z ) = P y + μ ( 2 v x 2 + 2 v y 2 + 2 v z 2 ) ρ ( u w x + v w y + w w z ) = P z + μ ( 2 w x 2 + 2 w y 2 + 2 w z 2 )
where u, v, and w are the velocity vector components along the x, y, and z directions, respectively, ρ is the fluid density, P is the fluid pressure, and µ is the dynamic viscosity.
In the model, the fluid domain was considered as water, with an average turbulent intensity of 5%, a temperature of 25 °C, fluid density of 997 kg m−3, and fluid viscosity of 0.9 10−3 Pa/s. The shear stress transport (SST) model was adopted to obtain a solution for the governing equations. In this model, we used k-ω for the region close to the wall, where the flow is of low velocity and flow is less, and k-ε for the rest of the fluid, where there is a greater influence from the turbulent regime [19]. The number of interactions in each simulation ranged from 50 to 1000, which generated residual errors smaller than 10−8. The computational mesh used in the study had a total of 1.08 × 107 elements. The simulations were performed on one desktop with the help of a twelve-core AMD Ryzen 9 3900X processor (Advanced Micro Devices, Santa Clara, CA, USA) with 32 GB of RAM. The simulations took around 15 h to complete.
As boundary conditions, we defined zero pressure at the exit of the model and a flow velocity of 0.5 mm s−1 to analyze the resistive effects of anatomical structures on the flow. Xu et al. [10] determined that an extended smooth segment of 25 µm long should be added at the inlet and outlet to avoid backflow effects [20].
To estimate the strengths of the different structures of the components of the vessels, the methodology proposed by Xu et al. [10] was used to obtain data by numerical simulation, resulting in the pressure profile (ΔP) and the average flow (q). The total resistance of the vessels was calculated from the relationship between ΔP and q.

2.3. Statistical Analysis

Data variance analyses were performed using the F distribution test (p < 0.05); the Chi-square test will be applied to assess the distribution of the anatomical characteristics of the vessels, and the Shapiro-Wilk normality test for the frequency of values. The anatomical characteristics of the factorial analysis of variance and differences will be determined by the Tukey test (p < 0.05).

3. Results and Discussion

Figure 2 and Figure 3 show the micrographs of the plants tested in this study. The secondary wall thickening was pitted, with openings throughout the longitudinal section of the xylem vessels (Figure 3).
Figure 2 and Figure 3 show the micrographs of the plants tested in this study. The secondary wall thickening was pitted, with openings throughout the longitudinal section of the xylem vessels (Figure 3).
Table 1 shows the mean geometric characteristics of the xylem vessels of irrigated and non-irrigated plants. We observed a significant difference in the number of pits and the relationship between the total pit area and the vessel’s surface area. However, the vessel length, diameter, and surface area did not significantly differ.
Donkor [21] reported diameters ranging from 100 and 275 µm and lengths between 121 and 683 µm in their review. Soranso et al. [22] worked with Khaya ivorensis, a variety of the same species, and has diameters between 41 and 141 µm. Unfortunately, few studies verify the anatomical characteristics of the vessels of Khaya grandifoliola for comparison, but the results mentioned above corroborate the values found in the current study.
With the definition of the geometric characteristics, it was possible to define the geometric model (Figure 4). The secondary wall structures were created using polygons perpendicular to the wall by an extrusion process along the axis of the vessels.
The pressure profile in the xylem vessels of irrigated and non-irrigated plants is observed in Figure 5 and Figure 6. In the inner part of the vessel, it was possible to observe a loss of energy of 0.49 Pa in irrigated and 0.47 Pa in the non-irrigated vessels. There is a gradual decrease in pressure between the “inlet” and the “outlet”, along with the entire longitudinal extension. It is possible to observe a non-gradual pressure passage zone with lower pressure than in the surrounding vessel area in the pit. It is still possible to verify that the pressure on the scores that form near the “inlet” is higher than those of the “outlet”.
The numerical predictions of pressure in the transition region between the vessel wall and the pit are shown in Figure 7 and Figure 8. There was a pressure variation in this region and in the entire longitudinal profile in both water regimes. This variation occurs as a function of the determined boundary condition, in which we tried to represent a semipermeable membrane. The radial flow is driven by an osmotic difference in solute concentration between internal and external liquids, representing the forces of cellular absorption. Unlike what was observed by Jensen et al. [20], the axial flow decreases due to the observed radial efflux. Note that the non-irrigated vessel presented higher variability in the pressure profile, due to the higher number of outlets (scores) that were observed anatomically (Table 1).
This pressure drop effect is hydraulically described by Christiansen [23]; the total friction loss equals the sum of the losses between the radial outlets, concentrating on the first outlets that make up 75% of the horizontal axis. Thus, the pressure will be minimal at the distal end and gradually increase towards the initial source. Consequently, the discharge at the end will be minimal, and, as we move towards the beginning, each exit from the pit will, in turn, discharge a little more flow due to the increase in pressure (Figure 7b,c and Figure 8b,c).
The maximum velocity distribution profile presented in Figure 9 highlights a parabolic distribution profile, typical in Newtonian incompressible fluids. In the wall and structure of the secondary wall, the formation of the boundary layer is evident, characterizing an outer flow region in which the velocity varies slowly and a possible tangential zone of velocity increase in the outputs of the scores. In Figure 10, a velocity profile was traced in the longitudinal section of the thickening of the second wall to observe this region and its interaction with the rest of the vessel. In the initial section (25–35 µm) of the non-irrigated vessel, it was possible to observe a variability of the profile of 56.5%, with the irrigated vessel having 19.2%. There was less variability in the vessel’s final section (last 40 µm), around 12.7% and 9.9%, respectively, for non-irrigated and irrigated vessels (Figure 11).
Comparing the velocity and pressure profiles, it was possible to observe that in the areas with a decrease in pressure due to the existence of pits, there was an increase in the flow velocity, mainly due to the decrease in the fluid movement area, indicating the interference of the secondary wall thickness in the flow. Xu et al. and Schulte et al. [7,24] modeled these flow dynamics in the pit and verified the same behavior, increasing the velocity at the center of the structure and decreasing static pressure at the outlet. These structures presented a significant resistance to flow (Figure 12). Xu et al. [24] further reiterate that in this region, Re is relatively small and the opening of the pit, being a laminar flow without compromising other structures, comprises a direct relationship between the opening diameter of the pit and its resistivity.
Partitioning the resistances, the smooth wall represented 65.30% and 62.57% of the total resistance, respectively, for irrigated and non-irrigated vessels. The non-irrigated vessel presented a resistance that was 3.04% lower than the irrigated vessel. Water deficit can suppress exchange=cell division and inhibit cell growth caused by turgor [25]. Non-irrigated vessels had a higher number of pits, which increased the resistance of the secondary wall of non-irrigated vessels. Zwieniecki and Secchi [26] relate this change in resistance to adaptation to ensure the hydraulic safety of vessels, as drought can affect anatomical characteristics, such as the structure of the secondary wall, the size of the conduit, and the density of the wood. The literature shows that the anatomical characteristics of the xylem adjust to drought conditions in trees, as drier conditions result in thicker cell walls and a more significant number of lateral structural elements, increasing the efficiency of water transport [26,27].
The secondary wall thickening represented 34.70% and 37.43% of the total resistance, respectively, for irrigated and non-irrigated vessels. Non-irrigated vessels showed higher resistance to the secondary wall thickening due to the higher number of pits, causing a more significant radial movement of fluids in these regions (Figure 8), translating a higher fluid volume to these structures. Thus, the smooth part of the wall had a lower resistance, consequently affecting the total resistance (Figure 10). Zimmermann [28] states that these structures cause a variation in conductivities along the length of the plant, mainly as a function of the vessel diameter and branch junctions, called “hydraulic bottlenecks”. These cellular structures present in the xylem walls influence the radial movement of water inside the plants, resulting from the decrease in diameter along the length of the stem. Therefore, the vessel presents a significant radial flow and exerts forces on the flow resistance [6,28,29,30]. The main challenge of the cohesion-tension theory is to explain how plants, particularly trees, that carry water at negative pressures can avoid cavitation or embolisms that interrupt the flow and result in dry-out [31]. The pit pores add substantial credibility to this theory, as they play essential roles in vessel resistivity [20].
The flow model via CFD in Khaya grandifoliola xylem vessels allowed us to better understand the flow phenomena of different resistances, secondary wall thickening, and smooth vessels. However, when the model is simulated with the universally diffused reading vessels analogous to Ohm’s law, we can see that the resistance has been underestimated at 68.7% and 61.4% for irrigated and non-irrigated vessels, respectively (Table 2). This is because this model uses two-dimensional mathematical models of the xylem structure and simple botanical analyses that neglect the other structural components present in the vessel [2,3]. Three-dimensionalizing and including the use of CFD in these microfluidic analyses is necessary due to the small size, varied structure, and complex internal flow of xylem vessels in plant species, which makes botanical observation tests very complicated, and it is often difficult to obtain accurate results because of the test conditions. Furthermore, other water transmission models are performed by macro-simulation models, which is different from the actual plant water ascent. Therefore, it is difficult to analyze the complex flow pattern and flow mechanism phenomenon within the catheter through experiments alone [9]. Our study and the various other studies cited above [7,8,9,10] contribute to further research into xylem vessel flow and the role of multiple structures in vessels.

4. Conclusions

The xylem vessels in Khaya grandifoliola provide a path of low resistance, where the smooth vessels present higher resistance, followed by thickening of the secondary wall. The simulations demonstrated a relationship between the total vessel strength and the total number of scores. Between the irrigated and non-irrigated plants, there was a difference of 2.8% of the total resistance. The non-irrigated vessel had a lower total resistivity but higher resistance to secondary wall thickening. Thus, the total number of crosshairs had a significant influence on radial transmission efficiency. More pits allow the plant to have more xylem vessels connected. Such connectivity will, in the case of an embolism event, allow water to be easily rerouted without drastically reducing the hydraulic conductivity.

Author Contributions

Conceptualization, methodology and writing—original draft preparation, D.S.d.A., D.H.M.d.M. and M.M.; writing—review and checking, R.A.F., R.B., G.G.S., F.P.d.D. and R.S.F.; project administration and funding acquisition, M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Council for Scientific and Technological Development (CNPq), grant number 141706/2020-1 and 146774/2020-5.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The authors thank the Postgraduate Program in Agronomy—Federal University of Goiás, (UFG); to the National Council for Scientific and Technological Development (CNPq), to the Coordination for the Improvement of Higher Education Personnel (CAPES) for the financial and structural support for carrying out this study and to the Regional Center for Technological Development and Innovation (CRTI) for the analyses in the scanning electron microscope.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Koch, G.W.; Sillett, S.C.; Jennings, G.M.; Davis, S.D. The Limits to Tree Height. Nature 2004, 428, 851–854. [Google Scholar] [CrossRef]
  2. Hall, A.J.; Minchin, P.E.H. A Closed-Form Solution for Steady-State Coupled Phloem/Xylem Flow Using the Lambert-W Function. Plant Cell Environ. 2013, 36, 2150–2162. [Google Scholar] [CrossRef] [Green Version]
  3. Lacointe, A.; Minchin, P.E.H. Modelling Phloem and Xylem Transport within a Complex Architecture. Funct. Plant Biol. 2008, 35, 772–780. [Google Scholar] [CrossRef]
  4. Loepfe, L.; Martinez-Vilalta, J.; Piñol, J.; Mencuccini, M. The Relevance of Xylem Network Structure for Plant Hydraulic Efficiency and Safety. J. Theor. Biol. 2007, 247, 788–803. [Google Scholar] [CrossRef]
  5. Patankar, S.V. Numerical Heat Transfer and Fluid Flow; CRC Press: Boca Raton, FL, USA, 2018; ISBN 978-1-315-27513-0. [Google Scholar]
  6. Schulte, P.J. Computational Fluid Dynamics Models of Conifer Bordered Pits Show How Pit Structure Affects Flow. New Phytol. 2012, 193, 721–729. [Google Scholar] [CrossRef]
  7. Schulte, P.J.; Hacke, U.G.; Schoonmaker, A.L. Pit Membrane Structure Is Highly Variable and Accounts for a Major Resistance to Water Flow through Tracheid Pits in Stems and Roots of Two Boreal Conifer Species. New Phytol. 2015, 208, 102–113. [Google Scholar] [CrossRef] [PubMed]
  8. QingLin, A.; Fang, X.; Qi, C.; JiaoLiao, C.; Peng, W. Flow resistance characteristics of scalariform perforation plates in plant xylem vessels. Nongye Jixie Xuebao Trans. Chin. Soc. Agric. Mach. 2011, 42, 110–143. [Google Scholar]
  9. Chen, Q.; Xu, F.; Ai, Q.; Zhang, L. Hydrodynamic Model and Flow Resistance Characteristics of Plant Vessel Wall Thickenings. Trans. Chin. Soc. Agric. Eng. 2015, 31, 1–8. [Google Scholar]
  10. Xu, T.; Zhang, L.; Li, Z. Computational Fluid Dynamics Model and Flow Resistance Characteristics of Jatropha Curcas L Xylem Vessel. Sci. Rep. 2020, 10, 14728. [Google Scholar] [CrossRef]
  11. Motta, J.P.; da Silva Oliveira, J.T.; Braz, R.L.; Duarte, A.P.C.; Alves, R.C. Caracterização da madeira de quatro espécies florestais. Cienc. Rural 2014, 44, 2186–2192. [Google Scholar] [CrossRef] [Green Version]
  12. Pinto, C.A.; David, J.S.; Cochard, H.; Caldeira, M.C.; Henriques, M.O.; Quilhó, T.; Paço, T.A.; Pereira, J.S.; David, T.S. Drought-Induced Embolism in Current-Year Shoots of Two Mediterranean Evergreen Oaks. For. Ecol. Manag. 2012, 285, 1–10. [Google Scholar] [CrossRef]
  13. Reboita, M.S.; Rodrigues, M.; Silva, L.F.; Alves, M.A. Aspectos Climáticos Do Estado de Minas Gerais (Climate Aspects in Minas Gerais State). Rev. Bras. Climatol. 2015, 17. [Google Scholar] [CrossRef] [Green Version]
  14. Gea-Izquierdo, G.; Fonti, P.; Cherubini, P.; Martín-Benito, D.; Chaar, H.; Cañellas, I. Xylem Hydraulic Adjustment and Growth Response of Quercus Canariensis Willd. to Climatic Variability. Tree Physiol. 2012, 32, 401–413. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Bryukhanova, M.; Fonti, P. Xylem Plasticity Allows Rapid Hydraulic Adjustment to Annual Climatic Variability. Trees 2013, 27, 485–496. [Google Scholar] [CrossRef] [Green Version]
  16. Moraes, M.D.A. De Qualidade E Potencial Energético Da Madeira de Desbaste de Mogno Africano. Master’s Thesis, Federal University of Goiás, Goiânia, Brazil, 2017. [Google Scholar]
  17. Trevisan, R.; Rosa, M.; Haselein, C.R.; Santini, E.J.; Gatto, D.A. Dimensões Das Fibras E Sua Relação Com A Idade de Transição Entre Lenho Juvenil E Adulto de Eucalyptus grandis W. Hill ex Maiden. Ciênc. Florest. 2017, 27, 1385–1393. [Google Scholar] [CrossRef] [Green Version]
  18. Shi, C.-L.; Butenko, M.A. Visualizing Morphological Changes of Abscission Zone Cells in Arabidopsis by Scanning Electron Microscope. In Plant Senescence: Methods and Protocols, Methods in Molecular Biology; Guo, Y., Ed.; Springer: New York, NY, USA, 2018; pp. 321–328. ISBN 978-1-4939-7672-0. [Google Scholar]
  19. Menter, F.R. Review of the Shear-Stress Transport Turbulence Model Experience from an Industrial Perspective. Int. J. Comput. Fluid Dyn. 2009, 23, 305–316. [Google Scholar] [CrossRef]
  20. Jensen, K.H.; Berg-Sørensen, K.; Bruus, H.; Holbrook, N.M.; Liesche, J.; Schulz, A.; Zwieniecki, M.A.; Bohr, T. Sap Flow and Sugar Transport in Plants. Rev. Mod. Phys. 2016, 88, 035007. [Google Scholar] [CrossRef] [Green Version]
  21. Donkor, B.N. Stem Wood Structure of Four Ghanaian Khaya Species. Ph.D. Thesis, Lakehead University, Faculty of Forestry, Thunder Bay, ON, Canada, 1997. [Google Scholar]
  22. Soranso, D.R.; Vidaurre, G.B.; Chagas, M.P.; Oliveira, J.T.d.S.; Silva, J.G.M.d.; Latorraca, J.V.d.F. Radial Growth Dynamics of Khaya Ivorensis Trees from Experimental Plantation. Rev. Árvore 2018, 42. [Google Scholar] [CrossRef]
  23. Christiansen, J.E. Hydraulics of Sprinkling Systems for Irrigation. Trans. Am. Soc. Civ. Eng. 1942, 107, 221–239. [Google Scholar] [CrossRef]
  24. Xu, F.; Chen, J.L.; Ai, Q.L.; Chen, Q. Modeling Fluid Flow in Angiosperms Xylem Pits. Appl. Mech. Mater. 2012, 195–196, 577–582. [Google Scholar] [CrossRef]
  25. Swidrak, I.; Gruber, A.; Oberhuber, W. Xylem and Phloem Phenology in Co-Occurring Conifers Exposed to Drought. Trees 2014, 28, 1161–1171. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Zwieniecki, M.A.; Secchi, F. Threats to Xylem Hydraulic Function of Trees under ‘New Climate Normal’ Conditions. Plant Cell Environ. 2015, 38, 1713–1724. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Martín, J.A.; Esteban, L.G.; de Palacios, P.; Fernández, F.G. Variation in Wood Anatomical Traits of Pinus Sylvestris L. between Spanish Regions of Provenance. Trees 2010, 24, 1017–1028. [Google Scholar] [CrossRef]
  28. Zimmermann, M.H. Xylem Structure and the Ascent of Sap; Springer Science & Business Media: Berlin, Germany, 2013; ISBN 978-3-662-22627-8. [Google Scholar]
  29. Schulte, P.J.; Castle, A.L. Water Flow Through Vessel Perforation Plates—The Effects of Plate Angle and Thickness for Liriodendron Tulipifera. J. Exp. Bot. 1993, 44, 1143–1148. [Google Scholar] [CrossRef]
  30. Chen, Q.; Xu, F.; Ai, Q.L.; Chen, J.L. Numerical Simulation of Water Transport through Vessel Perforation Plates. Available online: https://doi.org/10.4028/www.scientific.net/AMM.195-196.645 (accessed on 9 March 2020).
  31. Lucas, W.J.; Groover, A.; Lichtenberger, R.; Furuta, K.; Yadav, S.-R.; Helariutta, Y.; He, X.-Q.; Fukuda, H.; Kang, J.; Brady, S.M.; et al. The Plant Vascular System: Evolution, Development and FunctionsF. J. Integr. Plant Biol. 2013, 55, 294–388. [Google Scholar] [CrossRef]
Figure 1. Representation of the sampling points of Khaya grandifoliola. (a) longitudinal positions of the stem, at 0, 25, 50, 75, and 100% of the total height; (b) radial positions of the tree transversals, in the form of disc-shaped; (c) specimen removal point; and (d) fragment used for maceration removed from the specimen. Adapted from Moraes [16].
Figure 1. Representation of the sampling points of Khaya grandifoliola. (a) longitudinal positions of the stem, at 0, 25, 50, 75, and 100% of the total height; (b) radial positions of the tree transversals, in the form of disc-shaped; (c) specimen removal point; and (d) fragment used for maceration removed from the specimen. Adapted from Moraes [16].
Water 13 02723 g001
Figure 2. Detail of secondary wall thickening on light electron microscopy, along the superficial area of the vessels of Khaya grandifoliola.
Figure 2. Detail of secondary wall thickening on light electron microscopy, along the superficial area of the vessels of Khaya grandifoliola.
Water 13 02723 g002
Figure 3. Field emission scanning electron microscopy of the vessel of Khaya grandifoliola. (a) Full-length xylem vessel; (b,c) details of secondary wall thickening.
Figure 3. Field emission scanning electron microscopy of the vessel of Khaya grandifoliola. (a) Full-length xylem vessel; (b,c) details of secondary wall thickening.
Water 13 02723 g003
Figure 4. Three-dimensional anatomical model of xylem vessels of Khaya grandifoliola.
Figure 4. Three-dimensional anatomical model of xylem vessels of Khaya grandifoliola.
Water 13 02723 g004
Figure 5. Flow field pressure distribution in section along the axis of irrigated vessels.
Figure 5. Flow field pressure distribution in section along the axis of irrigated vessels.
Water 13 02723 g005
Figure 6. Flow field pressure distribution in section along the axis of non-irrigated vessels.
Figure 6. Flow field pressure distribution in section along the axis of non-irrigated vessels.
Water 13 02723 g006
Figure 7. Numerical predictions of the pressure profile (Pa) generated along with the anatomical model in the longitudinal section of the irrigated Khaya grandifoliola xylem vessels. (a) The entire longitudinal axis; (b) input detail; and (c) output detail.
Figure 7. Numerical predictions of the pressure profile (Pa) generated along with the anatomical model in the longitudinal section of the irrigated Khaya grandifoliola xylem vessels. (a) The entire longitudinal axis; (b) input detail; and (c) output detail.
Water 13 02723 g007
Figure 8. Numerical predictions of the pressure profile (Pa) generated along with the anatomical model in the longitudinal section of the non-irrigated Khaya grandifoliola xylem vessels. (a) The entire longitudinal axis; (b) input detail and (c) output detail.
Figure 8. Numerical predictions of the pressure profile (Pa) generated along with the anatomical model in the longitudinal section of the non-irrigated Khaya grandifoliola xylem vessels. (a) The entire longitudinal axis; (b) input detail and (c) output detail.
Water 13 02723 g008
Figure 9. Flow field velocity distribution in section along the axis of irrigated and non-irrigated vessels.
Figure 9. Flow field velocity distribution in section along the axis of irrigated and non-irrigated vessels.
Water 13 02723 g009
Figure 10. Numerical predictions of flow velocity generated along with the anatomical model in the longitudinal section of irrigated and non-irrigated Khaya grandifoliola xylem vessels.
Figure 10. Numerical predictions of flow velocity generated along with the anatomical model in the longitudinal section of irrigated and non-irrigated Khaya grandifoliola xylem vessels.
Water 13 02723 g010
Figure 11. Flow field velocity distribution in the secondary vessel wall structure of the irrigated Khaya grandifoliola xylem.
Figure 11. Flow field velocity distribution in the secondary vessel wall structure of the irrigated Khaya grandifoliola xylem.
Water 13 02723 g011
Figure 12. Total and structural component strength of the xylem vessel model in Khaya grandifoliola.
Figure 12. Total and structural component strength of the xylem vessel model in Khaya grandifoliola.
Water 13 02723 g012
Table 1. The mean geometric pattern of xylem vessels in the stems of irrigated and non-irrigated Khaya grandifoliola plants.
Table 1. The mean geometric pattern of xylem vessels in the stems of irrigated and non-irrigated Khaya grandifoliola plants.
TreatmentLength (µm)Diameter (µm)Reticle Area (µm2)N. of PitsRA/SAV 1
Irrigated359.97 ns137.01 ns1.80 ns4362 b5.13% b
Non-irrigated392.65 ns143.26 ns1.78 ns5732 a5.63% a
1 RA/SAV is the ratio between the total area of the pit and the vessel’s surface area; ns: non-significant; Means signed by different letters differ statistically by the F test.
Table 2. Comparisons of the mean values of total vessel resistance between those simulated and calculated by Ohm’s law, and analysis of their percentage difference (%Dif).
Table 2. Comparisons of the mean values of total vessel resistance between those simulated and calculated by Ohm’s law, and analysis of their percentage difference (%Dif).
Total Resistance
TreatmentsSimulatedOhm’s Law%Diff
Irrigated9.86 × 1083.08 × 10868.72%
Non-irrigated9.56 × 1083.68 × 10861.44%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

de Araujo, D.S.; de Moraes, D.H.M.; Mesquita, M.; Flores, R.A.; Battisti, R.; Santos, G.G.; de Deus, F.P.; Ferrarezi, R.S. Numerical Modeling of Microfluid Dynamics in Xylem Vessels of Khaya grandifoliola. Water 2021, 13, 2723. https://doi.org/10.3390/w13192723

AMA Style

de Araujo DS, de Moraes DHM, Mesquita M, Flores RA, Battisti R, Santos GG, de Deus FP, Ferrarezi RS. Numerical Modeling of Microfluid Dynamics in Xylem Vessels of Khaya grandifoliola. Water. 2021; 13(19):2723. https://doi.org/10.3390/w13192723

Chicago/Turabian Style

de Araujo, Daniel Somma, Diogo Henrique Morato de Moraes, Marcio Mesquita, Rilner Alves Flores, Rafael Battisti, Glenio Guimarães Santos, Fabio Ponciano de Deus, and Rhuanito Soranz Ferrarezi. 2021. "Numerical Modeling of Microfluid Dynamics in Xylem Vessels of Khaya grandifoliola" Water 13, no. 19: 2723. https://doi.org/10.3390/w13192723

APA Style

de Araujo, D. S., de Moraes, D. H. M., Mesquita, M., Flores, R. A., Battisti, R., Santos, G. G., de Deus, F. P., & Ferrarezi, R. S. (2021). Numerical Modeling of Microfluid Dynamics in Xylem Vessels of Khaya grandifoliola. Water, 13(19), 2723. https://doi.org/10.3390/w13192723

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop