A surface model of nonlinear, non-steady-state phloem transport

  • Received: 01 April 2016 Accepted: 01 February 2017 Published: 01 August 2017
  • MSC : Primary: 35K61, 92C80; Secondary: 65M60

  • Phloem transport is the process by which carbohydrates produced by photosynthesis in the leaves get distributed in a plant. According to Münch, the osmotically generated hydrostatic phloem pressure is the force driving the long-distance transport of photoassimilates. Following Thompson and Holbrook[35]'s approach, we develop a mathematical model of coupled water-carbohydrate transport. It is first proven that the model presented here preserves the positivity. The model is then applied to simulate the flow of phloem sap for an organic tree shape, on a 3D surface and in a channel with orthotropic hydraulic properties. Those features represent an significant advance in modelling the pathway for carbohydrate transport in trees.

    Citation: Youcef Mammeri, Damien Sellier. A surface model of nonlinear, non-steady-state phloem transport[J]. Mathematical Biosciences and Engineering, 2017, 14(4): 1055-1069. doi: 10.3934/mbe.2017055

    Related Papers:

    [1] Liwei Chen, Jaslyn Lee, Wei Ning Chen . The use of metabolic engineering to produce fatty acid-derived biofuel and chemicals in Saccharomyces cerevisiae: a review. AIMS Bioengineering, 2016, 3(4): 468-492. doi: 10.3934/bioeng.2016.4.468
    [2] Clara Navarrete, José L. Martínez . Non-conventional yeasts as superior production platforms for sustainable fermentation based bio-manufacturing processes. AIMS Bioengineering, 2020, 7(4): 289-305. doi: 10.3934/bioeng.2020024
    [3] George H McArthur IV, Pooja P Nanjannavar, Emily H Miller, Stephen S Fong . Integrative metabolic engineering. AIMS Bioengineering, 2015, 2(3): 93-103. doi: 10.3934/bioeng.2015.3.93
    [4] Pawel Jajesniak, Tuck Seng Wong . From genetic circuits to industrial-scale biomanufacturing: bacterial promoters as a cornerstone of biotechnology. AIMS Bioengineering, 2015, 2(3): 277-296. doi: 10.3934/bioeng.2015.3.277
    [5] Murtaza Shabbir Hussain, Gabriel M Rodriguez, Difeng Gao, Michael Spagnuolo, Lauren Gambill, Mark Blenner . Recent advances in bioengineering of the oleaginous yeast Yarrowia lipolytica. AIMS Bioengineering, 2016, 3(4): 493-514. doi: 10.3934/bioeng.2016.4.493
    [6] Edwige Arnold, Iness Hammami, Jingkui Chen, Sachin Gupte, Yves Durocher, Mario Jolicoeur . Overexpression of G6PDH does not affect the behavior of HEK-293 clones stably expressing interferon-α2b. AIMS Bioengineering, 2016, 3(3): 319-336. doi: 10.3934/bioeng.2016.3.319
    [7] Ozgun Firat Duzenli, Sezer Okay . Promoter engineering for the recombinant protein production in prokaryotic systems. AIMS Bioengineering, 2020, 7(2): 62-81. doi: 10.3934/bioeng.2020007
    [8] Flávia da Silva Fernandes, Luan Reis Honorato da Silva, Érica Simplício de Souza, Lívia Melo Carneiro, João Paulo Alves Silva, Steven Zelski, João Vicente Braga de Souza, Jacqueline da Silva Batista . Isolation, genetic identification of Amazonian yeasts and analysis of thermotolerance and alcohol tolerance of Saccharomyces cerevisiae from Theobroma grandiflorum and Eugenia stipitata. AIMS Bioengineering, 2024, 11(1): 24-43. doi: 10.3934/bioeng.2024003
    [9] Luan Reis Honorato da Silva, Flávia da Silva Fernandes, Jocélia Pinheiro Santos, Érica Simplício de Souza, Lívia Melo Carneiro, João Paulo Alves Silva, Jacqueline da Silva Batista, João Vicente Braga de Souza . Bioprospecting Saccharomyces cerevisiae in fruits from Amazonian region for beer brewing. AIMS Bioengineering, 2024, 11(2): 130-146. doi: 10.3934/bioeng.2024008
    [10] Bashir Sajo Mienda, Mohd Shahir Shamsir . An overview of pathway prediction tools for synthetic design of microbial chemical factories. AIMS Bioengineering, 2015, 2(1): 1-14. doi: 10.3934/bioeng.2015.1.1
  • Phloem transport is the process by which carbohydrates produced by photosynthesis in the leaves get distributed in a plant. According to Münch, the osmotically generated hydrostatic phloem pressure is the force driving the long-distance transport of photoassimilates. Following Thompson and Holbrook[35]'s approach, we develop a mathematical model of coupled water-carbohydrate transport. It is first proven that the model presented here preserves the positivity. The model is then applied to simulate the flow of phloem sap for an organic tree shape, on a 3D surface and in a channel with orthotropic hydraulic properties. Those features represent an significant advance in modelling the pathway for carbohydrate transport in trees.


    SiC ceramics are often used in electronic devices as a semiconductor component in the manufacture of MOSFET (Metal Oxide Semiconductor Field Effect Transistor), IGBT (Insulated Gate Bipolar Transistor), DMOS (Double diffused Metal Oxid Semiconductor), SiC–DMOS (Silicon Carbid Double diffused Metal Oxid Semiconductor), SiC–SBD (Silicon Carbid transistors0 and other components [1,2]. The semiconductor is usually connected to THE conductive part of a packaged unit [3].

    Soldering is the most used technology employed for bond formations between SiC semiconductors and conductive parts. Great attention is at present devoted to soldering with ultrasound assistance without the use of flux [4,5,6,7,8]. Many authors at present deal with the issues of soldering SiC ceramics with Cu [9,10,11].

    A group of scientists [12] has dealt with the research of a suitable solder for applications at higher service temperatures in the automotive industry. They solved the application of a power SiC device in electronic hybrid vehicles with a service temperature around 200 ℃. The capability to operate SiC-based power electronic modules at higher temperatures significantly simplifies the control system, which in turn reduces total costs. Most of the lead-free, Sn or Bi-based solders exert a melting point within the temperature range of 170–229 ℃, although these are not suitable for permanent service at higher application temperatures. Therefore the BiAgPd and BiAgSn (AgBiX®, Indium Corporation) solders were developed as an alternative to a fluxless solder with a melting point of 265 ℃ [13,14,15]. This solder was developed with the aim to improve wettability, which enables a more intense formation of intermetallic phases in the joint boundary.

    An alternative is solder with a matrix containing Ag particles [16,17,18]. This solder, which is destined for higher application temperatures, is often used for reactive soldering of SiC with aluminium oxide [19,20]. However, this effort is focused on the development of new universal solders that can join metal and ceramic materials in a direct way (without coatings or interlayers) by a flux-less method. The design of new solder alloys must also be based on the use of lead-free filler materials.

    In this paper, research is directed to the use of a new solder alloy for higher application temperatures based on Bi–Ag–Ti. The solder was used to solder the SiC semiconductor substrate with Cu. This combination is directly focused upon power semiconductor applications. In this paper, the suitability of this alloy for soldering of this combination was evaluated on the basis of standard research methods. The joint microstructure and analyses of interactions on the boundary were studied and the shear strength of the joint was measured.

    Ultrasonic soldering was employed for soldering SiC ceramics with Cu. The CERAN 500 hot plate with thermostatic control was used as the heat source for melting solder. Surface oxides were removed by application of Hanuz UZP 2 ultrasonic equipment at a frequency of 40 kHz. Ultrasonic vibrations were applied for 3 to 5 s. The parameters of ultrasonic soldering are given in Table 1.

    Table 1.  Soldering parameters.
    Parameters
    Ultrasound power 400 W
    Working frequency 40 kHz
    Amplitude 2 μm
    Soldering temperature 380 ℃
    Time of ultrasound activation 5 s

     | Show Table
    DownLoad: CSV

    Substrates of the following dimensions were used for soldering:

    ●  SiC ceramics in the form of a disc in diameter of Φ 15 mm × 3 mm

    ●  Copper in the form of a square in dimension of 10 × 10 × 2 mm

    An active Bi–Ag–Ti-based solder was applied for joint fabrication. The chemical composition of the soldering alloy is given in Table 2.

    Table 2.  Chemical composition of soldering alloy.
    Element Quantity (hm%)
    Bi 86
    Ag 11
    Ti 3

     | Show Table
    DownLoad: CSV

    Figure 1 shows a schematic procedure of ultrasonic soldering. The melting temperature of the soldering alloy was measured by DTA analysis. The microstructure of the soldered joint was observed by use of light microscopy on NEOPHOT 32 equipment. The microstructure of joints was analysed by electron scanning microscopy (SEM). EDX analysis was applied for analysis of SiC/Bi11Ag3Ti and Cu/Bi11Ag3Ti joint boundaries. A LabTest 5.250SP1-VM device was used for shear strength measurement. A jig with a defined specimen shape was used to change the direction of application of the tensile axial loading force on the specimen. This special jig provides uniform shear loading to the specimen in the plane of the interface of the solder and substrate.

    Figure 1.  Shows a schematic procedure of ultrasonic soldering.

    Figure 2 shows the result of DTA analysis of soldering alloy type Bi11Ag3Ti. The analysis follows the course of solder heating from room temperature to melting point. The first reaction took place at a temperature of 258 ℃. This temperature characterizes the eutectic transformation between Bi and Ag. It is the temperature at which the binary mixture of Bi–Ag begins to melt. The second reaction has occurred at a temperature of 366 ℃, which is when total fusion of Bi–Ag mixture occurs. The third reaction took place at 399 ℃. This reaction characterizes the fusion of the last phase of Ag–Ti. The melting point of Bi11Ag3Ti solder is 411 ℃. Based on DTA analysis, the soldering temperature was determined to be 380 ℃.

    Figure 2.  DTA analysis of soldering alloy type Bi11Ag3Ti.

    Figure 3 shows the microstructure of the SiC/Bi11Ag3Ti/Cu soldered joint. The lower part of the joint is formed of Cu and the upper part consists of SiC. It is obvious from the figure that the bond is sound and free from cracks or pores. An undulated boundary may be seen between the solder and copper, which is caused by the effect of ultrasonic vibrations. The results of light microscopy suggest the presence Ti phases, which are concentrated not only in solder volume but also on the boundary with the copper substrate.

    Figure 3.  Microstructure of SiC/Bi11Ag3Ti/Cu joint.

    EDX microanalysis was applied for the identification of constituents and phases. Figure 4 shows the results of line EDX microanalysis of the SiC/Bi11Ag3Ti joint.

    Figure 4.  The line EDX microanalysis of boundary of SiC/Bi11Ag3Ti joint.

    The close boundary of the SiC/Bi11Ag3Ti joint was mainly formed of Bi and Ti. The minor presence of Ag was also proven. The bond was formed owing to an increased concentration of Ti, which is obvious from Figure 4—the green curve. Ti is an active component of solder. The high affinity of Ti to elements of the ceramics component ensured the reaction between Ti and ceramics at the formation of a reaction band. Based on line analysis, the diffusion of Si to the volume of soldering alloy was identified.

    The diffused silicon is concentrated in the bond zone. Thus the reaction band is formed of the Ti–Si phase. The thickness of the reaction band at the point of measurement is approximately 1 µm. Figure 5 shows the map of elements in the boundary of SiC/Bi11Ag3Ti joint. The yellow color identifies the presence of Ti in the boundary. As is obvious from Figure 5, Ti was distributed to the boundary with SiC during the soldering process.

    Figure 5.  (a) Planar distribution of elements in the boundary of SiC/Bi11Ag3Ti joint (b) Si, (c) Bi, (d) Ag, and (e) Ti.

    The planar distribution of elements also proved the presence of Si in the boundary. The presence of Si in the boundary is depicted in Figure 5 by green color. In several sections of the joint the presence of silver was also identified in the boundary—Figure 5d. Silver occurs in the solder in the form of an intermetallic phase.

    Another analyzed boundary comprised the Cu/Bi11Ag3Ti joint. The EDX line analysis was performed in this boundary and its results are shown in Figure 6. No metallurgical bond was formed in the boundary between Cu and Bi11Ag3Ti solder. This finding was also proved by the authors in works [21,22]. Based on their observation, only a mechanical bond exists in the boundary. The presence of Ag and/or Ti was not identified in this boundary. Based on the line analysis it may be stated that the bond in a close boundary is formed by the mechanical mixture of Cu and Bi. The thickness of the reaction band in Bi11Ag3Ti/Cu joint was approximately 0.8 μm.

    Figure 6.  Line EDX microanalysis in the boundary of Cu/Bi11Ag3Ti joint.

    Finally, the shear strength of the joints soldered with Bi–Ag–Ti solder was measured. The shear strength was graphically evaluated and is shown in Figure 7. An important value for this work was the result of the average shear strength of the SiC/Cu joint, which is marked in red in the graph. The average value reached 40 MPa.

    Figure 7.  Shear strength of the joints soldered with Bi–Ag–Ti solder.

    For comparison, shear strength was also measured at the Cu/Cu joint. The average strength of this joint was 33.5 MPa. The SiC/Cu combination is a main combination for a large number of applications in modern electronics. In this study, this bond achieved a higher strength than the bond strength of metallic (Cu/Cu) substrates. The bond strength for applications in power semiconductor devices must be at least 40 MPa. A new solder alloy fulfills this condition.

    The aim of this study was to assess the suitability of a newly developed active solder type Bi–Ag–Ti for soldering the SiC/Cu combination. This Bi–Ag–Ti based solder belongs to solders for higher application temperatures. Its selection is oriented for application in power semiconductor devices. The achieved results may be characterised as follows:

    -The melting temperature of Bi11Ag3Ti solder and reactions during heating was determined on the basis of DTA analysis. Three reactions were observed. First reaction occurs at the temperature of 258 ℃, when the binary mixture of Bi–Ag begins to melt. A second reaction at 366 ℃ represents the end of melting of Bi–Ag mixture. The third reaction—melting of Ag–Ti phase, occurs at 399 ℃. The solder was completely melted at 411 ℃.

    -The microstructure of soldered joints was free from cracks, pores or other imperfections,

    -The boundary in SiC/Bi11Ag3Ti joint was predominantly formed of Si–Ti based phase. The thickness of the reaction layer in the boundary was approximately 1 μm. The presence of an active element Ti ensured the wetting of SiC substrate during soldering process,

    -The bond in the boundary of Bi11Ag3Ti/Cu joint was formed on the basis of formation of a mechanical mixture of Bi + Cu. The thickness of reaction layer in the boundary of SiC/Bi11Ag3Ti joint was approximately 0.8 μm. Titanium at the interface was not observed,

    -The average shear strength of SiC/Cu joint was 40 MPa. This value was higher than the shear strength of the Cu/Cu joint (33.5 MPa).

    The research of soldering the combination of SiC/Cu by use of Bi11Ag3Ti solder proved the suitability of application of the selected soldering alloy. Regarding the results of analyses of joint boundaries, the formation of a sound joint was proved. The average shear strength of SiC/Cu joint reached 40 MPa. This joint was even stronger than the Cu/Cu joint using the same solder. This value (40 MPa) of shear strength is suitable for applications in modern power semiconductor devices.

    This work was supported by the Slovak Research and Development Agency under the contract no. APVV-17-0025. The paper was also prepared with the support of the VEGA 1/0089/17 project: Research of new alloys for direct soldering of metallic and ceramic materials and project for supporting excellent teams of young researchers in conditions of Slovak Technical University-Project Young Researcher-excellent teams-no. 1392.

    All authors declare no conflicts of interest in this paper.

    [1] [ P. Cabrita,M. Thorpe,G. Huber, Hydrodynamics of steady state phloem transport with radial leakage of solute, Frontiers Plant Sci., 4 (2013): 531-543.
    [2] [ A. L. Christy,J. M. Ferrier, A mathematical treatment of Münch's pressure-flow hypothesis of phloem translocation, Plant Physio., 52 (1973): 531-538.
    [3] [ T. K. Dey,J. A. Levine, Delaunay meshing of isosurfaces, Visual Comput., 24 (2008): 411-422.
    [4] [ J. M. Ferrier, Further theoretical analysis of concentration-pressure-flux waves in phloem transport systems, Can. J. Bot., 56 (1978): 1086-1090.
    [5] [ F. G. Feugier,A. Satake, Dynamical feedback between circadian clock and sucrose availability explains adaptive response of starch metabolism to various photoperiods, Frontiers Plant Sci., 305 (2013): 1-11.
    [6] [ D. B. Fisher,C. Cash-Clark, Sieve tube unloading and post-phloemtransport of fluorescent tracers and proteins injected into sieve tubes via severed aphid stylets, Plant Physio., 123 (2000): 125-137.
    [7] [ J. D. Goeschl,C. E. Magnuson, Physiological implications of the Münch--Horwitz theory of phloem transport: effect of loading rates, Plant Cell Env., 9 (1986): 95-102.
    [8] [ J. Gričar,L. Krže,K. Čufar, Number of cells in xylem, phloem and dormant cambium in silver fir (Abies alba), in trees of different vitality, IAWA Journal, 30 (2009): 121-133.
    [9] [ J. Hansen,E. Beck, The fate and path of assimilation products in the stem of 8-year-old {Scots} pine (Pinus sylvestris {L}.) trees, Trees, 4 (1990): 16-21.
    [10] [ F. Hecht, New Developments in Freefem++, J. Num. Math., 20 (2012): 251-265.
    [11] [ L. Horwitz, Some simplified mathematical treatments of translocation in plants, Plant Physio., 33 (1958): 81-93.
    [12] [ T. Hölttä,M. Mencuccini,E. Nikinmaa, Linking phloem function to structure: Analysis with a coupled xylem-phloem transport model, J. Theo. Bio., 259 (2009): 325-337.
    [13] [ T. Hölttä,T. Vesala,S. Sevanto,M. Perämäki,E. Nikinmaa, Modeling xylem and phloem water flows in trees according to cohesion theory and Münch hypothesis, Trees, 20 (2006): 67-78.
    [14] [ W. Hundsdorfer and J. G. Verwer, Numerical Solution of Time-Dependent Advection-Diffusion-Reaction Equations Springer Series in Comput. Math., 33, Springer, 2003.
    [15] [ K. H. Jensen,J. Lee,T. Bohr,H. Bruus,N. M. Holbrook,M. A. Zwieniecki, Optimality of the Münch mechanism for translocation of sugars in plants, J. R. Soc. Interface, 8 (2011): 1155-1165.
    [16] [ A. Kagawa,A. Sugimoto,T. C. Maximov, CO 2 pulse-labelling of photoassimilates reveals carbon allocation within and between tree rings, Plant Cell Env., 29 (2006): 1571-1584.
    [17] [ E. M. Kramer, Wood grain pattern formation: A brief review, J. Plant Growth Reg., 25 (2006): 290-301.
    [18] [ H.-O. Kreiss and J. Lorenz, Initial-Boundary Value Problems and the Navier-Stokes Equations Classics in Applied Mathematics, SIAM, 2004.
    [19] [ A. Lacointe,P. E. H. Minchin, Modelling phloem and xylem transport within a complex architecture, Funct. Plant Bio., 35 (2008): 772-780.
    [20] [ A. Lang, A model of mass flow in the phloem, Funct. Plant Bio., 5 (1978): 535-546.
    [21] [ P. E. H. Minchin,M. R. Thorpe,J. F. Farrar, A simple mechanistic model of phloem transport which explains sink priority, Journal of Experimental Botany, 44 (1993): 947-955.
    [22] [ E. Münch, Die Stoffbewegungen in der Pflanze Jena, Gustav Fischer, 1930.
    [23] [ K. A. Nagel,B. Kastenholz,S. Jahnke,D. van Dusschoten,T. Aach,M. Mühlich,D. Truhn,H. Scharr,S. Terjung,A. Walter,U. Schurr, Temperature responses of roots: Impact on growth, root system architecture and implications for phenotyping, Funct. Plant Bio., 36 (2009): 947-959.
    [24] [ E. M. Ouhabaz, Analysis of Heat Equations on Domains London Math. Soc. Monographs Series, Princeton University Press, 2005.
    [25] [ S. Payvandi,K. R. Daly,K. C. Zygalakis,T. Roose, Mathematical modelling of the phloem: The importance of diffusion on sugar transport at osmotic equilibrium, Bull. Math Biol., 76 (2014): 2834-2865.
    [26] [ S. Pfautsch,J. Renard,M. G. Tjoelker,A. Salih, Phloem as capacitor: Radial transfer of water into xylem of tree stems occurs via symplastic transport in ray parenchyma, Plant Physio., 167 (2015): 963-971.
    [27] [ O. Pironneau,M. Tabata, Stability and convergence of a Galerkin-characteristics finite element scheme of lumped mass type, Int. J. Num. Meth. Fluids, 64 (2000): 1240-1253.
    [28] [ G. E. Phillips,J. Bodig,J. Goodman, Flow grain analogy, Wood Sci., 14 (1981): 55-64.
    [29] [ R. J. Phillips,S. R. Dungan, Asymptotic analysis of flow in sieve tubes with semi-permeable walls, J. Theor. Biol., 162 (1993): 465-485.
    [30] [ D. Rotsch, T. Brossard, S. Bihmidine, W. Ying, V. Gaddam, M. Harmata, J. D. Robertson, M. Swyers, S. S. Jurisson and D. M. Braun, Radiosynthesis of 6'-Deoxy-6'[18F]Fluorosucrose via automated synthesis and its utility to study in vivo sucrose transport in maize (Zea mays) leaves PLoS ONE 10 (2015), e0128989.
    [31] [ D. Sellier,J. J. Harrington, Phloem transport in trees: A generic surface model, Eco. Mod., 290 (2014): 102-109.
    [32] [ D. Sellier,M. J. Plank,J. J. Harrington, A mathematical framework for modelling cambial surface evolution using a level set method, Annals Bot., 108 (2011): 1001-1011.
    [33] [ R. Spicer, Symplasmic networks in secondary vascular tissues: Parenchyma distribution and activity supporting long-distance transport, J. Exp. Bot., 65 (2014): 1829-1848.
    [34] [ J. F. Swindells, C. F. Snyder, R. C. Hardy and P. E. Golden, Viscosities of sucrose solutions at various temperatures: Tables of recalculated values, NBS Circular 440 (1958).
    [35] [ M. V. Thompson,N. M. Holbrook, Application of a single-solute non-steady-state phloem model to the study of long-distance assimilate transport, J. Theo. Bio., 220 (2003): 419-455.
  • This article has been cited by:

    1. Felix Weber, Markus Rettenmayr, Joining of SiO2 glass and 316L stainless steel using Bi–Ag-based active solders, 2021, 56, 0022-2461, 3444, 10.1007/s10853-020-05426-4
    2. Rudimylla Septimio, Clarissa Cruz, Bismarck Silva, Amauri Garcia, José Eduardo Spinelli, Noé Cheung, Microstructural and segregation effects affecting the corrosion behavior of a high‐temperature Bi‐Ag solder alloy in dilute chloride solution, 2021, 0021-891X, 10.1007/s10800-021-01533-5
    3. Tomáš Meluš, Roman Koleňák, Jaromír Drápala, Paulína Babincová, Matej Pašák, Ultrasonic soldering of Al2O3 ceramics and Ni-SiC composite by use of Bi-based active solder, 2023, 10, 2372-0484, 213, 10.3934/matersci.2023012
    4. Tatsuya Kobayashi, Tetsuya Ando, Recent Development of Joining and Conductive Materials for Electronic Components, 2021, 62, 1345-9678, 1270, 10.2320/matertrans.MT-M2021060
  • Reader Comments
  • © 2017 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4282) PDF downloads(536) Cited by(6)

Figures and Tables

Figures(8)  /  Tables(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog