Review

Mathematical methods for modeling the microcirculation

  • Received: 02 March 2017 Accepted: 09 May 2017 Published: 07 June 2017
  • The microcirculation plays a major role in maintaining homeostasis in the body. Alterations or dysfunctions of the microcirculation can lead to several types of serious diseases. It is not surprising, then, that the microcirculation has been an object of intense theoretical and experimental study over the past few decades. Mathematical approaches offer a valuable method for quantifying the relationships between various mechanical, hemodynamic, and regulatory factors of the microcirculation and the pathophysiology of numerous diseases. This work provides an overview of several mathematical models that describe and investigate the many different aspects of the microcirculation, including geometry of the vascular bed, blood flow in the vascular networks, solute transport and delivery to the surrounding tissue, and vessel wall mechanics under passive and active stimuli. Representing relevant phenomena across multiple spatial scales remains a major challenge in modeling the microcirculation. Nevertheless, the depth and breadth of mathematical modeling with applications in the microcirculation is demonstrated in this work. A special emphasis is placed on models of the retinal circulation, including models that predict the influence of ocular hemodynamic alterations with the progression of ocular diseases such as glaucoma.

    Citation: Julia C. Arciero, Paola Causin, Francesca Malgaroli. Mathematical methods for modeling the microcirculation[J]. AIMS Biophysics, 2017, 4(3): 362-399. doi: 10.3934/biophy.2017.3.362

    Related Papers:

    [1] Waraporn Apiwatanapiwat, Prapassorn Rugthaworn, Pilanee Vaithanomsat, Warunee Thanapase, Akihiko Kosugi, Takamitsu Arai, Yutaka Mori, Yoshinori Murata . Ethanol production at high temperature from cassava pulp by a newly isolated Kluyveromyces marxianus strain, TISTR 5925. AIMS Energy, 2013, 1(1): 3-16. doi: 10.3934/energy.2013.1.3
    [2] S.O Adegoke, O.A Falode, A.A Adeleke, P.P Ikubanni, O.O Agboola . Production of an alternative fuel from a blend of polypropylene wastes and Jatropha distillates for Oil Well drilling. AIMS Energy, 2020, 8(6): 1127-1142. doi: 10.3934/energy.2020.6.1127
    [3] Hans Mattila, Dina Kačar, Tuulia Mali, Taina Lundell . Lignocellulose bioconversion to ethanol by a fungal single-step consolidated method tested with waste substrates and co-culture experiments. AIMS Energy, 2018, 6(5): 866-879. doi: 10.3934/energy.2018.5.866
    [4] Rumana Islam, Charushi Panditharatne, John Schellenberg, Richard Sparling, Nazim Cicek, David B. Levin . Potential of thin stillage as a low-cost nutrient source for direct cellulose fermentation by Clostridium thermocellum. AIMS Energy, 2015, 3(4): 711-727. doi: 10.3934/energy.2015.4.711
    [5] Fitriani Tupa R. Silalahi, Togar M. Simatupang, Manahan P. Siallagan . A system dynamics approach to biodiesel fund management in Indonesia. AIMS Energy, 2020, 8(6): 1173-1198. doi: 10.3934/energy.2020.6.1173
    [6] Gustavo Aguilar, Pranjali D. Muley, Charles Henkel, Dorin Boldor . Effects of biomass particle size on yield and composition of pyrolysis bio-oil derived from Chinese tallow tree (Triadica Sebifera L.) and energy cane (Saccharum complex) in an inductively heated reactor. AIMS Energy, 2015, 3(4): 838-850. doi: 10.3934/energy.2015.4.838
    [7] Michaela Weissgram, Thomas Ters, Hedda K. Weber, Christoph Herwig . Investigating the potential of thermophilic species for ethanol production from industrial spent sulfite liquor. AIMS Energy, 2015, 3(4): 592-611. doi: 10.3934/energy.2015.4.592
    [8] Lifita N. Tande, Valerie Dupont . Autothermal reforming of palm empty fruit bunch bio-oil: thermodynamic modelling. AIMS Energy, 2016, 4(1): 68-92. doi: 10.3934/energy.2016.1.68
    [9] Divakar Dahiya, Poonam Singh Nigam . Bioethanol synthesis for fuel or beverages from the processing of agri-food by-products and natural biomass using economical and purposely modified biocatalytic systems. AIMS Energy, 2018, 6(6): 979-992. doi: 10.3934/energy.2018.6.979
    [10] Yadessa Gonfa Keneni, Jorge Mario Marchetti . Oil extraction from plant seeds for biodiesel production. AIMS Energy, 2017, 5(2): 316-340. doi: 10.3934/energy.2017.2.316
  • The microcirculation plays a major role in maintaining homeostasis in the body. Alterations or dysfunctions of the microcirculation can lead to several types of serious diseases. It is not surprising, then, that the microcirculation has been an object of intense theoretical and experimental study over the past few decades. Mathematical approaches offer a valuable method for quantifying the relationships between various mechanical, hemodynamic, and regulatory factors of the microcirculation and the pathophysiology of numerous diseases. This work provides an overview of several mathematical models that describe and investigate the many different aspects of the microcirculation, including geometry of the vascular bed, blood flow in the vascular networks, solute transport and delivery to the surrounding tissue, and vessel wall mechanics under passive and active stimuli. Representing relevant phenomena across multiple spatial scales remains a major challenge in modeling the microcirculation. Nevertheless, the depth and breadth of mathematical modeling with applications in the microcirculation is demonstrated in this work. A special emphasis is placed on models of the retinal circulation, including models that predict the influence of ocular hemodynamic alterations with the progression of ocular diseases such as glaucoma.


    1. Introduction

    The oil palm (Elaeis guineensis) is one of the most rapidly expanding equatorial crops globally. Plantations are replanted every 25 years due to decreased oil productivity [1]. At least 120,000 hectares of oil palm plantation was estimated to be replanted annually in Malaysia from 2006 to 2010 [2], and 450,000 hectares of oil palm plantation are expected to be replanted annually in Indonesia in the next 25 years. Palm trunks are not used as lumber due to their high moisture content (70-80% based on total mass), resulting in significant warping after drying [3]. The only current practical use for felled oil palm trunks is as plywood [4]. Recently, palm sap squeezed from old oil palm trunks (OPTs) was shown to contain large quantities of sugar, and ethanol was produced from this sugar by fermentation by Saccharomyces cerevisiae [5,6]. Sap containing a high content of available sugars (glucose 85.2 g/L) was found in the inner part of the OPT and found to also contain a significant amount of amino acids (serine and glutamic acid), minerals such as calcium and manganese, and several organic acids (such as malic acid and maleic acid). This sap was therefore used directly as a fermentable medium for ethanol or lactic acid production, without supplementation with nitrogen or minerals [6], and could be used as feedstock for bioethanol production without the need for complicated pretreatment such as the saccharification of lignocellulose. However, some palm saps contained a low content of sugar, depending on the species of oil palm or the physiology of the plant. Here we condensed sugar in low-sugar-content sap using flat membrane filtration in order to efficiently utilize the sap.

    S. cerevisiaeis used as an fermentation yeast for ethanol production in the brewing industry. However, several unique yeasts have been isolated from natural habitats that exhibit unique characteristics for fuel ethanol production using various bioresources as substrates under various conditions [7,8,9]. Kluyveromyces marxianusis is an thermotolerant yeast for bioethanol production at elevated temperatures [10]. The advantages of thermotolerant yeast for industrial ethanol production are: 1. The high fermentation activity of thermotolerant yeasts at elevated temperature decreases the chances of bioreactor contamination. 2. The optimum temperature for S. cerevisiae is 28-32°C whereas K. marxianus has a wide optimum temperature range, from 28 to 45 °C. 3. Cooling and distillation costs are decreased when fermentation is conducted at 40 °C. 4. Increased enzyme activity at higher temperature should decrease the need for added enzyme and thus decrease the overall costs for exogenous enzyme [11].

    We isolated the thermotolerant yeast, K. marxianus TISTR5925, from rotting fruit in Thailand. This yeast strain can produce ethanol by hydrolyzing cassava pulp, the waste product from starch factories [12]. In this study, we compared the growth of three thermotolerant yeast strains, K. marxianus NCYC578, NCYC2791, and TISTR5925 with three non-thermotolerant yeast strains, S. cerevisiae kyokai7, NCYC3233, and S288c at elevated temperature. In addition, we studied the fermentation of condensed sap from oil palm trunk by K. marxinanus TISTR5925 and S. cerevisiae Kyokai7 as representative strains of the two species at temperatures from 37 °C to 50 °C. The energy required to concentrate the sugar in the sap by flat membrane filtration was measured. To evaluate ethanol production from palm sap, the total energy required to prepare and concentrate the palm sap was taken as the input energy, and the potential energy obtained from the produced ethanol and the effective energy obtained from the solid residue after squeezing sap from oil palm trunk were taken as the output energy. The Net Energy Ratio (NER) is the proportion of the output energy to the input energy; the NER for different contents of sugar in the sap were compared to identify the sugar content required for economically-viable ethanol production.

    2. Materials and Method

    2.1. Preparation of sap from old oil palm trunk

    OPTs 23-25 years old were harvested from several plantations in Johor province (1o 28′0″ North, 103°45′0″ East), Malaysia. In initial studies, 26 oil palm (E. guinesis) trees were cut in Malaysia and their sap was squeezed [13]. At this time, logged trunks were cut into 1.2 m lengths, and then debarked to the trunk core (1.2 m long × 20 cm diameter). Cores from three trees were prepared. The trunk cores (approx. 30 kg) were shredded and squeezed and palm sap was obtained using the method of Murata et al. [13]; a total of 17 kg sap and 10.1 kg solid residue after extraction of sap was obtained from one trunk core [13]. The squeezed sap and the residues of oil palm trunk were stored at −20°C to avoid contamination of microbes. After thawing and analysis sugar content, the sap was used for ethanol fermentation and the effective energy obtained from combustion of the residue was estimated.

    The content of various sugars in the sap was measured by HPLC (TOSOH, Osaka, Japan) using a CARBOSep CHO-682 LEAD column (Transgenomic Inc., Omaha, NE, USA), a temperature of 80 °C, a flow rate of 0.4 mL/min, and a refractive index (RI) detector (TOSOH RI-8020). The sugar components were identified by comparing their retention times with those of authentic standards.

    2.2. Membrane filtration of the palm sap with thin sugar

    When the total sugar in the sap was less than 5%, the sap was condensed to at least 10% total sugar using a flat membrane filtration system built by Sanko Ltd. (Kanagawa, Japan) and consisting of a membrane filtration device (LabStak M20; Alfa Laval, Lund, Sweden) and a high pressure pump (Wanner Hydracell G10-X; Wanner Engineering, Minneapolis, MN, USA). The palm sap was pumped into the membrane filtration device using a high pressure pump and then condensed by the reverse osmosis membrane (FILMTEC membrane, BW30-365 high rejection brackish water RO element, Dow Chemical Co., Minneapolis, MN, USA) in the LabStak M20.

    2.3. Growth of the yeasts on YPD plates and ethanol production from oil palm sap

    K. marxianus TISTR5925 was isolated from rotting fruit from a local market in Thailand [12]. S. cerevisiae Kyokai7 was purchased from the National Research Institute of Brewing (Hiroshima, Japan), and S. cerevisiae NCYC3233, and K. marxiaus NCYC587 and NCYC2791, were purchased from the National Collection of Yeast Cultures (Norwich, United Kingdom). Yeast extract, peptone, and agar were purchased from Difco Laboratories (Detroit, MI, USA). Yeast extract peptone dextrose (YPD) medium comprised 1% yeast extract, 2% peptone, and 2% glucose. Liquid YPD medium was used to prepare plates by adding 2% agar. 10% YPD (1% yeast extract, 2% peptone, 10% glucose) was used as a control fermentation medium for comparison with palm sap and for pre-cultivation. All chemicals, except for the medium, were purchased from Wako Pure Chemical Industries Ltd. (Osaka, Japan).

    K. marxianus (NCYC587, NCYC2791, TISTR5925) and S. cerevisiae (Kyokai7, NCYC3233, S288c) were precultured in liquid 10% YPD medium at 30 °C for 24 h, then streaked on YPD plates and incubated at 30 °C and 45 °C. Of the six strains, TISTR5925 and Kyokai7 were used as representative strains of K. marxianus and S. cerevisiae, respectively, to ferment palm sap. Precultures of these strains (3 mL) were centrifuged for 5 min at 4,000 x g and the pellets were washed twice with sterile water. The cells were resuspended in 3 mL palm sap (without added nitrogen sources), and ethanol fermentation was conducted at temperatures between 37 °C to 50 °C. Each fermentation was repeated at least three times. The ethanol produced and glucose consumed were monitored simultaneously by HPLC using an Aminex HPX87-H column (300 mm × 7.8 mm; Bio-Rad, Hercules, CA, USA). The ethanol and glucose peaks were identified by comparing the retention times with those of standards [14]. The column was operated at 65 °C and a flow rate of 0.6 ml/min using sulfuric acid (5.0 mM) as the mobile phase. The data were processed using Chromato-PRO (Runtime Instruments, Kanagawa, Japan).

    2.4. Measurement of energy consumed for flat membrane filtration

    The power (kW) required to concentrate the sap was measured using an electric power meter (Clamp on Power Hi Tester 3169-01, HIOKI E.E. Corp. Ueda, Nagano, Japan). The input energy was then converted into thermal terms using equation 1:

    The energy required for sap content (mega joule: MJ)

    = (Average power (kW) x running time (h) x 3.6 MJ)/0.4 (Equation 1)

    Units and factors:

    *a Average power (kW): Average of the power required to concentrate the palm sap

    *b kWh: = 3.6 MJ = 3.6 × 106 joule

    *c 0.4: the conversion coefficient used in Japan for electric power produced from petroleum oil

    The output energy was estimated from the energy in the ethanol produced from the palm sap. The effective heating value for ethanol was obtained by subtracting the heat of vaporization of water from the heating value of pure ethanol, and used to estimate the net energy ratio. The effective heating value of ethanol is 21.2 MJ L-1. The power required to drive the other devises (rotary lathe, shredder, and mill) to shred and squeeze the oil palm trunk were taken from Murata et al. [13].

    3. Results and Discussion

    3.1. Condensation of palm sap with thin sugar by flat membrane filtration

    We developed a system consisting of a shredder and mills for squeezing palm sap from debarked old oil palm trunk cores, and the saps were extracted from several oil palm trunks [13]. The sugar components in the sap were identified by their HPLC peaks as sucrose, glucose and fructose (Figure 1A Untreated sap). The major sugar was glucose (28.8 g/L), followed by fructose (around 2.4 g/L); little sucrose was detected (Table 1 Untreated sap). In contrast, the major sugar in unfermented sap from oil palm tree grown in Nigeria is sucrose, with both glucose and fructose present as minor sugars [15,16]. In our earlier research, glucose, sucrose and fructose were contained in the sap from oil palm trunks [10,20]. A significant amount of sucrose was detected in oil palm sap from one trunk (Figure 1 and Table 1, Condensed sap 2) but little (Figure 1 and Table 1, Condensed sap 1) or none (Table 1, Condensed sap 3) in oil palm sap from other trunks, indicating that the content of specific sugars depends on the species of oil palm, its physiological state, and cultivation conditions.

    Figure 1. Sugar peaks in sap squeezed from oil palm trunks. HPLC analysis of sugar components in sap squeezed from oil palm trunks. untreated sap, condensed sap 1, condensed sap 2, condensed sap3. Peak 1: sucrose, Peak 2: glucose, Peak 3: fructose. No peak corresponding to sucrose was observed in condensed sap 1 and condensed sap3 but a peak corresponding to sucrose was observed in condensed sap 2. The X axis shows the refractive index, and the Y axis shows the retention time of the sugar.
    Table 1. The sugar components in palm sap.
    g/L
    Sucrose Glucose Fructose Total
    Untreated sap -28.82.4 31.2
    Condensed sap 12.574.716.8 94.0
    Condensed sap 250.7 65.5 17.8134.0
    Condensed sap 3-42.616.6 59.2
     | Show Table
    DownLoad: CSV

    Sucrose and starch are the major products of photosynthetic CO2 fixation. Triose phosphate (e.g., glycerol aldehyde-3-phosphate; 3PGA) is produced during CO2 fixation and is transported from the chloroplast to the cytosol, where it enters the sucrose biosynthetic pathway. Sucrose is biosynthesized by the phosphatase reaction via formation of sucrose-6-phosphate from fructose-6-phosphate and uridine diphosphate glucose (UDP-glucose) [17]. In contrast, starch is synthesized via an elongation reaction of adenine diphosphate glucose (ADP-glucose) from triose phosphate remaining in the chloroplast. Starch and sucrose synthesis is regulated by the inorganic phosphate (Pi) content in the cytosol [18]. Triose phosphate transporter, located in the chloroplast envelope membrane, exchanges triose phosphate for Pi between the stroma and cytoplasm. If the level of Pi in the cytosol is decreased, triose phosphate is not exported from the chloroplast stroma to the cytosol. Consequently, triose phosphate remains in the stroma and enters the starch biosynthesis pathway, and starch accumulates in the chloroplasts [17,18]. Starch synthesis is also favored by high light conditions, high temperature, and by the presence of triose phosphate [18]. Therefore, starch and sucrose synthesis are regulated in a complex fashion by the physiological condition of the oil palm tree. Facilitation of starch degradation results in high-glucose sap, whereas facilitation of sucrose synthesis produces sap rich in sucrose. Since sucrose is degraded into glucose and fructose by intrinsic invertase, high sucrose palm sap would also contain glucose and fructose. We found that the sap extracted from old palm trunks grown in Johor province contained almost glucose, but the few trunks had the sap contained sucrose except for glucose (Table 1 Condensed sap2).

    Several squeezed saps had low contents of sugar (31.2 g/L), and the ethanol produced from these saps would be low content (less than 15.6 g/L). Therefore, the sap with thin sugar showed the insufficient yield on ethanol production, because of waste of much energy than the stored in ethanol for distillation. We used a flat membrane filtration approach to condense the sugars and allow efficient ethanol production. Three palm saps containing less than 30-40 g/L sugar were condensed 3 to 4-fold (from 18 L to 5.4 L), then the sucrose, glucose, and fructose were quantified (Table 1). Condensed sap 2 exhibited a large sucrose peak (50.7 g/L, Table 1 Condensed sap 2), whereas condensed sap1 contained 2.5 g/L sucrose, and condensed sap 3 contained too little sucrose to quantify (Table 1 Condensed sap 3). Sucrose is a disaccharide comprising glucose and fructose joined by a glycosidic linkage, and is the major transportable carbohydrate in plants. Since saps squeezed from different trunks contained very different amounts of sucrose, the amount of sucrose in palm sap depends on plant growth conditions.

    3.2. Selection and Evaluation of thermotolerant yeasts for ethanol production from the sap

    We incubated K. marxianus (TISTR5925, NCYC587, NCYC2791) and S. cerevisiae (Kyokai7, S288c, NCYC3233) on YPD plates at 30 °C and 45°Cand compared their growth. As shown in Figure 2, the K. marxianus strains showed good growth at 45 °C, whereas the S. cerevisiaestrains essentially did not grow (Figure 2).

    Figure 2. Growth of S. cerevisiae and K. marxianus on YPD plates at 30 °C and 45°C. S. cerevisiae strains (Kyokai7, S288c, and NCYC3233) and K. marxianus (TISTR5925, NCYC587, and NCYC2791) were streaked on YPD plates, and incubated at 30 °C (A) and 45 °C (B).

    There are several reports regarding bioethanol production by yeasts other than the representative species, S. cerevisiae [7,8,9]. We isolated the thermotolerant yeast, K. marxianus TISTR5925, in Thailand, and observed good ethanol production due to its high fermentation activity at various temperatures and its ability to ferment various sugars [12].

    Condensed palm sap no.1 (Table 1) was fermented by the representative yeast strains, TISTR5925 and Kyokai7, at various temperatures. TISTR5925 consumed all the sugars and produced ethanol at temperatures between 37 °C and 45°C, but not at 50°C(Figure 3a and b), whereas Kyokai7 did not produce ethanol efficiently at 42 °C and 45°C, and sugars remained in the condensed sap (Figure 3c and d). K. marxianus TISTR5925 showed high fermentation activity from 37 °C to 45°Cbecause of its thermotolerance, whereas S. cerevisiae Kyokai7 showed little fermentation at 42 °C and 45 °C due to heat shock stress. The ethanol yields for K. marxianus TISTR5925 and S. cerevisiae Kyokai7 at each temperature were calculated and compared. The theoretical value for the complete conversion of the sugar in the palm sap (96.3 g/L) to ethanol is 49.2 g/L (Table 2). The ethanol yield at each temperature was calculated as the ratio with this theoretical value (Table 2a and b). K. marxianus TISTR5925 provided high ethanol yields, from 74.8% to 92.2%, at temperatures up to 45 °C, whereas S. cerevisiae Kyokai7 provided low ethanol yields (from 49.6% to 68.0%) above 42 °C. Yeast cells are killed by temperature fluctuations result in decreased ethanol production [11]. The thermotolerant yeast, K. marxianus TISTR5925, holds promise for ethanol production from oil palm sap at high temperature because of its fermentation activity over a range of temperatures.

    Figure 3. Fermentation of palm sap by TISTR5925 and Kyokai7. Palm sap was condensed to 10% sugars by flat membrane filtration. Both strains were precultured on 10% YPD medium, resuspended in condensed palm sap, and incubated at various temperatures. The experiments were repeated three times and the average values were plotted.
    Table 2. Fermentation of palm sap.
    a. TISTR5925
    Temperature, °C Time, h Max ethanol, g/L Ethanol yield*1, % Ethanol productivity*2, g/L/h
    0 7 18 24 48 72
    *1 Glucose (1.00 g/L) is converted to a theoretical value of 0.511g/L ethanol. If the sugar content of palm sap was 96.26 g/L, it was converted to 49.23 g/L ethanol (theoretical value). The ethanol yield at each temperature was calculated as the ratio with the theoretical value (49.23 g/L). *2 Ethanol productivity from 0 h to 7 h was calculated as initial rate.
    37 0.0 9.8 41.7 42.7 33.7 32.5 42.7 86.7 1.40
    40 0.0 10.7 45.4 42.6 34.6 31.5 45.4 92.2 1.52
    42 0.0 12.9 41.6 37.9 30.6 28.0 41.6 84.5 1.84
    45 0.0 5.9 31.6 36.8 35.6 35.0 36.8 74.8 0.84
    50 0.0 2.5 1.9 2.1 1.2 1.5 2.5 5.1 0.36
    a. TISTR5925
    Temperature, °C Time, h Max ethanol, g/L Ethanol yield*1, % Ethanol productivity*2, g/L/h
    0 7 18 24 48 72
    37 0.0 8.4 44.5 40.8 33.4 32.0 44.5 90.4 1.20
    40 0.0 7.4 43.0 42.7 38.2 35.8 43.0 87.3 1.06
    42 0.0 3.5 23.7 33.2 33.5 36.3 33.5 68.0 0.49
    45 0.0 3.4 9.4 10.3 23.7 24.4 24.4 49.6 0.49
    50 0.0 2.7 2.2 2.1 1.6 1.9 2.7 5.5 0.39
     | Show Table
    DownLoad: CSV

    Besides, the ethanol productivity from 0 h to 7 h of K. marxianus TISTR5925 was faster maximum 3.7 folds at 42 °C than ones by S. cerevisiae Kyokai7 (Table 2). This quick fermentation showed the advantage on the prevention of contamination and downsizing of fermentation tank, therefore, the fermentation by thermotolerant yeast would lead to cost saving in plant construction.

    3.3. Net e nergy ratio for the ethanol production using the flat membrane filtration and the fermentation with the thermotolerant yeast

    The energy required to condense the sap was estimated. The content of 18 L (sugar content: 31.2 g/L) of sap to 5.4 L (sugar content: 96.3 g/L) took 1.8 h and required 1.22 kWh of electricity (calculated from equation 1, Materials and Methods), or 10.9 MJ (Table 3a). The consumed power in order to condense sugar in sap depended on the initial content of sugar. If the initial sugar content in sap was too low, the much more energy could be wasted for condensation of palm sap.

    Table 3. Energy required to produce EtOH from condensed sap (Input energy).
    Sap squeezing *1, MJ Initial sugar content *2, g/L Final sugar content, g/L Difference between initial and final sugar content *3 Sap condensation *4, MJ EtOH fermentation *5, MJ Total input energy *6, MJ
    *1 The energy required to squeeze sap is cited from Murata et al., 2013.
    *2 Sugar content before condensation by flat membrane filtration
    *3 Difference between initial and final sugar content in sap: Final sugar content, 96.3g/L − Initial sugar, g/L
    *4 The energy required to condense sap, 10.9MJ was used for different content of sugar 65.1 g/L (= 96.3g/L − 31.2 g/L). Energy required to condense the palm saps were calculated as follows. Energy for sap condensation, MJ = 10.9 MJ × (Difference from b to g)/65.1g/L
    *5 The energy required for EtOH production is cited from Murata et al., 2013.
    *6 Total input energy = (Sap squeezing, MJ) + (Sap condensation, MJ) + (EtOH production, MJ)
    a 5.8 31.2 96.3 65.1 10.9 0.85 17.6
    b 5.8 41.2 96.3 55.1 9.3 0.85 15.9
    c 5.8 51.2 96.3 45.1 7.6 0.85 14.2
    d 5.8 61.2 96.3 35.1 5.9 0.85 12.5
    e 5.8 71.2 96.3 25.1 4.2 0.85 10.9
    f 5.8 81.2 96.3 15.1 2.5 0.85 9.2
    g 5.8 91.2 96.3 5.1 0.9 0.85 7.5
     | Show Table
    DownLoad: CSV

    We estimated the initial content of sugar in the saps to range from 41.2 g/L to 91.2 g/L to evaluate profitability for ethanol production from palm sap, and we calculated the amount of electricity required to concentrate the sap sugars to 96.3 g/L (Table 3 Initial sugar). The required energy was found to range from 0.9 MJ to 10.9 MJ, depending on the initial content of sugar in the sap (Table 3 Sap condensation). The energy required for squeezing the sap and for ethanol fermentation was estimated as 5.8 MJ and 0.85 MJ, respectively (Table 3 Sap squeezing and EtOH fermentation). The total input energy was calculated as the sum of the joules required for squeezing and condensing the sap and for ethanol production (Table 3 Total input energy) and was found to depend on the initial content of sugar in the sap.

    The sugars in sap squeezed from oil palm trunks were analyzed by HPLC as described in Materials and Methods. Saps with low sugar content (31.2 g/L) were condensed, from 17 L to 5.4 L, using a flat membrane filter. Fermentation of condensed palm sap (96.3 g/L; Table 2) by K. marxianus (Figure 3) provided ethanol (45.4 g/L); the volume of ethanol was calculated as 0.31 L according to the equation in Table 4 (Table 4 EtOH). The stored energy in the ethanol was calculated as 6.7 MJ from the energy content of ethanol (21.2 MJ/L) and the volume (0.31 L; Table 4 Energy from EtOH). If the residue obtained after squeezing sap from the oil palm trunk was used as a fuel for the distillation boiler, the amount of residue after extracting sap (17 kg) from an oil palm log (30 kg) was 10.1 kg [13] with a potential energy of 75.8 MJ (Table 4 Potential energy). We estimated that 20% of this energy could be utilized, giving an effective energy of 15.2 MJ (Table 4 Effective energy) [13].

    Table 4. The energy from ethanol and the residue (Output), and Net Energy Ratio (NER) of output to input energy.
    Folds of sugar content in sap *1 Sap*2, L EtOH *3, L Energy from EtOH *4, MJ Effective energy from the residues*5, MJ Total output energy *6, MJ Total input energy *7, MJ Net Energy Ratio *8
    *1 Fold of sugar content in sap = 96.3g/L /(Initial sugar, g/L; Citation From Table 3)
    *2 Sap, L = 17 L/(folds of sugar content)
    *3 Ethanol was quantified by HPLC as 45.4 g/L from the fermentation data in Table 2 Ethanol (L) was calculated as follows; EtOH, L = (0.0454kg/L × Sap volume, L)/0.789 kg/L
    *4 Energy in ethanol; 21.2 MJ L−1 Energy from EtOH, MJ = (21.2 MJ/L) × (each EtOH, L)
    *5 Average weight of residue after squeezing sap was 10.1 kg. Low heat value (LHV) of residue with 50% moisture was estimated as 7.5 MJ kg−1. Potential heat, MJ = (Residue after squeezing; 10.10 kg) x 7.5 = 75.8 MJ Effective energy was estimated as 20% of potential heat of residue. Effective energy of the residues = 75.8 MJ × 0.2 = 15.2 MJ
    *6 Total energy = (Energy from EtOH, MJ) + (Effective heat of residues, MJ)
    *7 Values cited from Table 3
    *8 Net Energy Ratio = (Total output energy) / (Total input energy)
    a 3.1 5.5 0.31 6.7 15.2 21.9 17.6 1.2
    b 2.3 7.3 0.41 8.8 15.2 24.0 15.9 1.5
    c 1.9 9.0 0.52 10.9 15.2 26.1 14.2 1.8
    d 1.6 10.8 0.62 13.1 15.2 28.3 12.5 2.3
    e 1.4 12.6 0.72 15.2 15.2 30.4 10.9 2.8
    f 1.2 14.3 0.82 17.3 15.2 32.5 9.2 3.5
    g 1.1 16.1 0.92 19.5 15.2 34.7 7.5 4.6
     | Show Table
    DownLoad: CSV

    Consequently, total output energy was calculated as 21.8 MJ (the sum of the energy in the produced ethanol and the effective energy from the residue; Table 4 Total output energy). The output energy depended on the volume of the sap and EtOH (Figure 4 Sap and EtOH), which in turn were associated with the initial content of sugar (Table 3 Initial sugar). Therefore, the output energy depended on the initial sugar content.

    The ratio of output energy to input energy was calculated (Table 4 Net Energy Ratio). A high initial sugar content provided a high ethanol output energy. An initial sugar content of 31.2 g/L required a total input energy of 17.6 MJ (Table 3a), whereas a sugar content of 96.3 g/L (5.5 L) provided 0.31 L of ethanol (6.7 MJ; Table 4a). This value is less than the 17.6 MJ required as total input energy for sap content and ethanol production. A sustainable process requires that the input energy required for ethanol production should not exceed the output energy that can be extracted from the biomass. An initial sugar content in sap of 61.2 g/L (Table 3c, Initial sugar) required a total input energy of 12.5 MJ to squeeze and concentrate the sap, and to produce ethanol. The output energy from the 0.62 L of ethanol produced was 13.1 MJ (Table 4d Energy from EtOH), providing an NER of output energy to input energy of 1.04 (data not shown). Therefore, positive energy efficiency requires a sugar content in the raw sap above 61.2 g/L (Table 3d and Table 4d).

    The biomass residue should be used for energy production in order to improve the energy efficiency of the process [19]. Squeezing an oil palm log (20 cm diameter, 1.2 m long, approx. 30 kg) provided 17 L of sap and 10.1 kg of residue. The potential energy of this residue was calculated to be 75.8 MJ from the low heat value (LHV) of residue containing 50% water (Table 4 Potential energy), and the effective energy was estimated to be 20% of the potential energy of the residue (Table 4 Effective energy).

    An effective energy of 15.2 MJ was estimated, and this energy was added as output energy (Table 4 Total output energy). Adding the effective energy (15.2 MJ) using the residue as fuel during distillation provides a total output energy of 13.1 MJ + 15.2 MJ = 28.3 MJ (Table 4c total output energy). Consequently, the energy efficiency of ethanol production from oil palm sap was improved by using the residue as fuel (Table 4 Net Energy Ratio). In this study the themotolerant yeast, K. marxianus TISTR5925, was used to ferment palm sap and was shown to have several advantageous features compared to S. cerevisiae due to its thermotolerance. Ethanol fermentation conducted between 35 °C to 40 °C reduces costs by an estimated 3, 000, 000 yen per year, for 30, 000 kL ethanol produced per year [2]. Therefore, less energy is required to produce ethanol using thermotolerant yeast due to decreased the size of fermentation tank and distillation tower cooling costs, thus improving the energy efficiency of the process.

    4. Conclusion

    Sap with a low content of sugar was obtained by squeezing old oil palm trunks. The sugar in the sap was condensed and fermented into ethanol using K. marxianus TISTR5925. The input energy required to condense the sap and for ethanol production, and the output energy contained in the heat energy of the ethanol, were calculated and compared. At a sugar content above 61.2 g/L, the output energy exceeded the input energy.

    We estimated the energy from combustion of the residue after squeezing sap from oil palm trunk. Using the residue as fuel during distillation and improving fermentation by using thermotolerant yeast can increase the efficiency of ethanol production from oil palm sap.

    Acknowledgments

    This study was supported by the Japan International Research Center for Agricultural Sciences (JIRCAS).

    Conflict of Interests

    The authors declare no conflict of interests.

    [1] Tuma R, Duran WN, Ley K (2008) Microcirculation. 2Eds., Elsevier.
    [2] Gutterman DD, Chabowski DS, Kadlec AO, et al. (2016) The human microcirculation: regulation of flow and beyond. Circ Res 118: 157–172. doi: 10.1161/CIRCRESAHA.115.305364
    [3] Mortillaro N (1983) The physiology and pharmacology of the microcirculation. Elsevier.
    [4] Dalkara T (2015) Cerebral Microcirculation: An Introduction, In: Lanzer P, Editor, Springer Berlin Heidelberg, 655–680.
    [5] Navar LG, Arendshorst WJ, Pallone TL, et al. (2011) The renal microcirculation. Compr Physiol 9: 550–683.
    [6] Kvietys PR (2010) The gastrointestinal circulation. Morgan & Claypool Life Sciences.
    [7] Ivanov KP (2013) Circulation in the lungs and microcirculation in the alveoli. Respir Physiol Neurobiol 187: 26–30. doi: 10.1016/j.resp.2013.02.022
    [8] Murray A, Dinsdale G (2016) Imaging the microcirculation. Microcirculation 23: 335–336. doi: 10.1111/micc.12282
    [9] Leahy M (2012) Microcirculation Imaging. Wiley-Blackwell.
    [10] Eriksson S, Nilsson J, Sturesson C (2014) Non-invasive imaging of microcirculation: a technology review. Med Devices 7: 445–452.
    [11] Secomb TW, Beard DA, Frisbee JC, et al. (2008) The role of theoretical modeling in microcirculation research. Microcirculation 15: 693–698. doi: 10.1080/10739680802349734
    [12] Lee J, Smith NP (2008) Theoretical modeling in hemodynamics of microcirculation. Microcirculation 15: 699–714. doi: 10.1080/10739680802229589
    [13] Popel AS, Johnson PC (2005) Microcirculation and hemorheology. Annu Rev Fluid Mech 37: 43–69. doi: 10.1146/annurev.fluid.37.042604.133933
    [14] Gompper G, Fedosov DA (2016) Modeling microcirculatory blood flow: current state and future perspectives. Wiley Interdiscip Rev Syst Biol Med 8: 157–168. doi: 10.1002/wsbm.1326
    [15] Secomb TW (2017) Blood flow in the microcirculation. Annu Rev Fluid Mech 49: 443–461. doi: 10.1146/annurev-fluid-010816-060302
    [16] Liu D, Wood NB, Witt N, et al. (2009) Computational analysis of oxygen transport in the retinal arterial network. Curr Eye Res 34: 945–956. doi: 10.3109/02713680903230079
    [17] Gould IG, Linninger AA (2015) Hematocrit distribution and tissue oxygenation in large microcirculatory networks. Microcirculation 22: 1–18. doi: 10.1111/micc.12156
    [18] Reichold J, Stampanoni M, Keller AL, et al. (2009) Vascular graph model to simulate the cerebral blood flow in realistic vascular networks. J Cereb Blood Flow Metab 29: 1429–1443. doi: 10.1038/jcbfm.2009.58
    [19] Guibert R, Fonta C, Plouraboué F (2010) A new approach to model confined suspensions flows in complex networks: application to blood flow. Transport Porous Med 83: 171–194. doi: 10.1007/s11242-009-9492-0
    [20] Muller LO, Toro EF (2014) Enhanced global mathematical model for studying cerebral venous blood flow. J Biomech 47: 3361–3372. doi: 10.1016/j.jbiomech.2014.08.005
    [21] Ursino M, Lodi CA (1997) A simple mathematical model of the interaction between intracranial pressure and cerebral hemodynamics. J Appl Physiol 82: 1256–1269.
    [22] Arciero JC, Carlson BE, Secomb TW (2008) Theoretical model of metabolic blood flow regulation: roles of ATP release by red blood cells and conducted responses. Am J Physiol Heart Circ Physiol 295: H1562–H1571. doi: 10.1152/ajpheart.00261.2008
    [23] Takahashi T, Nagaoka T, Yanagida H, et al. (2009) A mathematical model for the distribution of hemodynamic parameters in the human retinal microvascular network. J Biorheol 23: 77–86. doi: 10.1007/s12573-009-0012-1
    [24] Gabryś E, Rybaczuk M, Kędzia A (2005) Fractal models of circulatory system. Symmetrical and asymmetrical approach comparison. Chaos, Solitons & Fractals 24: 707–715.
    [25] Olufsen MS, Peskin CS, Kim WY, et al. (2000) Numerical simulation and experimental validation of blood flow in arteries with structured-tree outflow conditions. Ann Biomed Eng 28: 1281–1299. doi: 10.1114/1.1326031
    [26] David T, Alzaidi S, Farr H (2009) Coupled autoregulation models in the cerebro-vasculature. J Eng Math 64: 403–415. doi: 10.1007/s10665-009-9274-2
    [27] Ganesan P, He S, Xu H (2011) Development of an image-based model for capillary vasculature of retina. Comput Meth Prog Bio 102: 35–46. doi: 10.1016/j.cmpb.2010.12.009
    [28] Schröder S, Schmid-Schönbein GW, Schmid-Schönbein H, et al. (1990) A method for recording the network topology of human retinal vessels. Klin Monbl Augenheilkd 197: 33–39. doi: 10.1055/s-2008-1046240
    [29] Causin P, Guidoboni G, Malgaroli F, et al. (2016) Blood flow mechanics and oxygen transport and delivery in the retinal microcirculation: multiscale mathematical modeling and numerical simulation. Biomech Model Mechanobiol 15: 525–542. doi: 10.1007/s10237-015-0708-7
    [30] Pan Q, Wang R, Reglin B, et al. (2014) A one-dimensional mathematical model for studying the pulsatile flow in microvascular networks. J Biomech Eng 136: 011009.
    [31] Linninger AA, Gould IG, Marinnan T, et al. (2013) Cerebral microcirculation and oxygen tension in the human secondary cortex. Ann Biomed Eng 41: 2264–2284. doi: 10.1007/s10439-013-0828-0
    [32] Gould I, Tsai P, Kleinfeld D, et al. (2017) The capillary bed offers the largest hemodynamic resistance to the cortical blood supply. J Cereb Blood Flow Metab 37: 52–68. doi: 10.1177/0271678X16671146
    [33] Fry BC, Roy TK, Secomb TW (2013) Capillary recruitment in a theoretical model for blood flow regulation in heterogeneous microvessel networks. Physiol Rep 1: e00050.
    [34] Su SW, Catherall M, Payne S (2012) The influence of network structure on the transport of blood in the human cerebral microvasculature. Microcirculation 19: 175–187. doi: 10.1111/j.1549-8719.2011.00148.x
    [35] Safaeian N, David T (2013) A computational model of oxygen transport in the cerebrocapillary levels for normal and pathologic brain function. J Cereb Blood Flow Metab 33: 1633–1641. doi: 10.1038/jcbfm.2013.119
    [36] Soltani M, Chen P (2013) Numerical modeling of interstitial fluid flow coupled with blood flow through a remodeled solid tumor microvascular network. PLoS One 8: e67025. doi: 10.1371/journal.pone.0067025
    [37] Sefidgar M, Soltania M, Raahemifarc K, et al. (2015) Numerical modeling of drug delivery in a dynamic solid tumor microvasculature. Microvasc Res 99: 43–56. doi: 10.1016/j.mvr.2015.02.007
    [38] Levine H, Pamuk S, Sleeman B, et al. (2002) Mathematical modeling of capillary formation and development in tumor angiogenesis: penetration into the stroma. Bull Math Biol 63: 801–863.
    [39] Rieger H, Welter M (2015) Integrative models of vascular remodeling during tumor growth. Wiley Interdiscip Rev Syst Biol Med 7: 113–129. doi: 10.1002/wsbm.1295
    [40] Erbertseder K, Reichold J, Flemisch B, et al. (2012) A coupled discrete/continuum model for describing cancer-therapeutic transport in the lung. PLoS One 7: e31966. doi: 10.1371/journal.pone.0031966
    [41] Dash RK, Bassingthwaighte JB (2006) Simultaneous blood-tissue exchange of oxygen, carbon dioxide, bicarbonate, and hydrogen ion. Ann Biomed Eng 34: 1129–1148. doi: 10.1007/s10439-005-9066-4
    [42] Moschandreou TE, Ellis CG, Goldman D (2011) Influence of tissue metabolism and capillary oxygen supply on arteriolar oxygen transport: a computational model. Math Biosci 232: 1–10. doi: 10.1016/j.mbs.2011.03.010
    [43] Welter M, Rieger H (2016) Computer simulations of the tumor vasculature: applications to interstitial fluid flow, drug delivery, and Oxygen supply. Adv Exp Med Biol 936: 31–72. doi: 10.1007/978-3-319-42023-3_3
    [44] Preziosi L, Tosin A (2009) Multiphase and multiscale trends in cancer modelling. Math Model Nat Phenom 4: 1–11.
    [45] Cornelis F, Sautd O, Cumsillea P, et al. (2013) In vivo mathematical modeling of tumor growth from imaging data: Soon to come in the future? Diagn Interv Imaging 94: 593–600. doi: 10.1016/j.diii.2013.03.001
    [46] Wu J, Xu S, Long Q, et al. (2008) Coupled modeling of blood perfusion in intravascular, interstitial spaces in tumor microvasculature. J Biomech 41: 996–1004. doi: 10.1016/j.jbiomech.2007.12.008
    [47] Ye GF, Moore TW, Jaron D (1993) Contributions of oxygen dissociation and convection to the behavior of a compartmental oxygen transport model. Microvasc Res 46: 1–18. doi: 10.1006/mvre.1993.1031
    [48] Piechnik SK, Chiarelli PA, Jezzard P (2008) Modelling vascular reactivity to investigate the basis of the relationship between cerebral blood volume and flow under CO2 manipulation. Neuroimage 39: 107–118. doi: 10.1016/j.neuroimage.2007.08.022
    [49] Fantini S (2014) Dynamic model for the tissue concentration and oxygen saturation of hemoglobin in relation to blood volume, flow velocity, and oxygen consumption: Implications for functional neuroimaging and coherent hemodynamics spectroscopy (CHS). Neuroimage 85: 202–221. doi: 10.1016/j.neuroimage.2013.03.065
    [50] Ursino M, Lodi CA (1998) Interaction among autoregulation, CO2 reactivity, and intracranial pressure: a mathematical model. Am J Physiol 274: H1715–H1728.
    [51] Spronck B, Martens EG, Gommer ED, et al. (2012) A lumped parameter model of cerebral blood flow control combining cerebral autoregulation and neurovascular coupling. Am J Physiol Heart Circ Physiol 303: H1143–H1153. doi: 10.1152/ajpheart.00303.2012
    [52] Payne S (2006) A model of the interaction between autoregulation and neural activation in the brain. Math Biosci 204: 260–281. doi: 10.1016/j.mbs.2006.08.006
    [53] Diamond SG, Perdue KL, Boas DA (2009) A cerebrovascular response model for functional neuroimaging including dynamic cerebral autoregulation. Math Biosci 220: 102–117. doi: 10.1016/j.mbs.2009.05.002
    [54] Muller LO, Toro EF (2014) A global multiscale mathematical model for the human circulation with emphasis on the venous system. Int J Numer Method Biomed Eng 30: 681–725. doi: 10.1002/cnm.2622
    [55] Gutierrez G (2004) A mathematical model of tissue-blood carbon dioxide exchange during hypoxia. Am J Respir Crit Care Med 169: 525–533. doi: 10.1164/rccm.200305-702OC
    [56] Vazquez AL, Masamoto K, Kim SG (2008) Dynamics of oxygen delivery and consumption during evoked neural stimulation using a compartment model and CBF and tissue P(O2) measurements. Neuroimage 42: 49–59. doi: 10.1016/j.neuroimage.2008.04.024
    [57] Barrett MJ, Suresh V (2013) Extra permeability is required to model dynamic oxygen measurements: evidence for functional recruitment? J Cereb Blood Flow Metab 33: 1402–1411. doi: 10.1038/jcbfm.2013.74
    [58] Fang Q, Sakadžić S, Ruvinskaya L, et al. (2008) Oxygen advection and diffusion in a three dimensional vascular anatomical network. Opt Express 16: 17530–17541. doi: 10.1364/OE.16.017530
    [59] Lorthois S, Cassot F, Lauwers F (2011) Simulation study of brain blood flow regulation by intra-cortical arterioles in an anatomically accurate large human vascular network: Part I: methodology and baseline flow. Neuroimage 54: 1031–1042. doi: 10.1016/j.neuroimage.2010.09.032
    [60] Boas DA, Jones SR, Devor A, et al. (2008) A vascular anatomical network model of the spatio-temporal response to brain activation. Neuroimage 40: 1116–1129. doi: 10.1016/j.neuroimage.2007.12.061
    [61] Takahashi T (2014) Microcirculation in fractal branching networks. Springer Japan.
    [62] Ganesan P, He S, Xu H (2010) Analysis of retinal circulation using an image-based network model of retinal vasculature. Microvasc Res 80: 99–109. doi: 10.1016/j.mvr.2010.02.005
    [63] El-Bouri WK, Payne SJ (2015) Multi-scale homogenization of blood flow in 3-dimensional human cerebral microvascular networks. J Theor Biol 380: 40–47. doi: 10.1016/j.jtbi.2015.05.011
    [64] Gagnon L, Sakadžić S, Lesage F, et al. (2015) Multimodal reconstruction of microvascular-flow distributions using combined two-photon microscopy and Doppler optical coherence tomography. Neurophotonics 2: 015008–015008. doi: 10.1117/1.NPh.2.1.015008
    [65] Tsoukias NM, Goldman D, Vadapalli A, et al. (2007) A computational model of oxygen delivery by hemoglobin-based oxygen carriers in three-dimensional microvascular networks. J Theor Biol 248: 657–674. doi: 10.1016/j.jtbi.2007.06.012
    [66] Park CS, Payne SJ (2016) Modelling the effects of cerebral microvasculature morphology on oxygen transport. Med Eng Phys 38: 41–47. doi: 10.1016/j.medengphy.2015.09.004
    [67] Gorodnova NO, Kolobov AV, Mynbaev OA, et al. (2016) Mathematical modeling of blood flow alteration in microcirculatory network due to angiogenesis. Lobachevskii J Math 37: 541–549. doi: 10.1134/S199508021605005X
    [68] Cristini V, Kassab GS (2005) Computer modeling of red blood cell rheology in the microcirculation: a brief overview. Ann Biomed Eng 33: 1724–1727. doi: 10.1007/s10439-005-8776-y
    [69] Sriram K, Intaglietta M, Tartakovsky DM (2014) Non-Newtonian flow of blood in arterioles: consequences for wall shear stress measurements. Microcirculation 21: 628–639. doi: 10.1111/micc.12141
    [70] Kim S, Ong PK, Yalcin O, et al. (2009) The cell-free layer in microvascular blood flow. Biorheology 46: 181–189.
    [71] Pries AR, Secomb TW (2003) Rheology of the microcirculation. Clin Hemorheol Microcirc 29: 143–148.
    [72] Ganesan P, He S, Xu H (2011) Modelling of pulsatile blood flow in arterial trees of retinal vasculature. Med Eng Phys 33: 810–823. doi: 10.1016/j.medengphy.2010.10.004
    [73] Secomb TW, Pries AR (2013) Blood viscosity in microvessels: experiment and theory. C R Phys 14: 470–478. doi: 10.1016/j.crhy.2013.04.002
    [74] Pries AR, Secomb TW (2005) Microvascular blood viscosity in vivo and the endothelial surface layer. Am J Physiol Heart Circ Physiol 289: H2657–H2664. doi: 10.1152/ajpheart.00297.2005
    [75] Lipowsky HH, Usami S, Chien S (1980) In vivo measurements of "apparent viscosity" and microvessel hematocrit in the mesentery of the cat. Microvasc Res 19: 297–319. doi: 10.1016/0026-2862(80)90050-3
    [76] Haynes RH (1960) Physical basis of the dependence of blood viscosity on tube radius. Am J Physiol 198: 1193–1200.
    [77] Neofytou P (2004) Comparison of blood rheological models for physiological flow simulation. Biorheology 41: 693–714.
    [78] Hellums JD, Nair PK, Huang NS, et al. (1996) Simulation of intraluminal gas transport processes in the microcirculation. Ann Biomed Eng 24: 1–24. doi: 10.1007/BF02649700
    [79] Hellums JD (1977) The resistance to oxygen transport in the capillaries relative to that in the surrounding tissue. Microvasc Res 13: 131–136. doi: 10.1016/0026-2862(77)90122-4
    [80] Formaggia L, Quarteroni A (2009) Cardiovascular mathematics: modeling and simulation of the circulatory system. Milan: Springer-Verlag.
    [81] Pries AR, Ley K, Claassen M, et al. (1989) Red cell distribution at microvascular bifurcations. Microvasc Res 38: 81–101. doi: 10.1016/0026-2862(89)90018-6
    [82] Canic S, Kim EH (2003) Mathematical analysis of the quasilinear effects in a hyperbolic model blood flow through compliant axi-symmetric vessels. Math Method Appl Sci 26: 1161–1186. doi: 10.1002/mma.407
    [83] Gagnon L, Smith AF, Boas DA, et al. (2016) Modeling of cerebral oxygen transport based on in vivo microscopic imaging of microvascular network structure, blood flow, and oxygenation. Front Comput Neurosci 1.
    [84] Nair PK, Huang NS, Hellums JD, et al. (1990) A simple model for prediction of oxygen transport rates by flowing blood in large capillaries. Microvasc Res 39: 203–211. doi: 10.1016/0026-2862(90)90070-8
    [85] Pandey HM, Negi DS, Bisht MS (2011) The study of mathematical modelling of human blood circulatory system. Int J Math Scis Appl 1: 425–434.
    [86] Sharan M, Popel AS (2001) A two-phase model for flow of blood in narrow tubes with increased effective viscosity near the wall. Biorheology 38: 415–428.
    [87] Namgung B, Ju M, Cabrales P, et al. (2013) Two-phase model for prediction of cell-free layer width in blood flow. Microvasc Res 85: 68–76. doi: 10.1016/j.mvr.2012.10.006
    [88] Chebbi R (2015) Dynamics of blood flow: modeling of the Fåhræus-Lindqvist effect. J Biol Phys 41: 313–326. doi: 10.1007/s10867-015-9376-1
    [89] Verma SR, Srivastava A (2014) Analytical study of a two-phase model for steady flow of blood in a circular tube. IJERA 12: 1–10.
    [90] Das B, Johnson PC, Popel AS (1998) Effect of nonaxisymmetric hematocrit distribution on non-Newtonian blood flow in small tubes. Biorheology 35: 69–87. doi: 10.1016/S0006-355X(98)00018-3
    [91] Gupta BB, Nigam KM, Jaffrin MY (1982) A three-layer semi-empirical model for flow of blood and other particular suspensions through narrow tubes. J Biomech Eng 104: 129. doi: 10.1115/1.3138326
    [92] Walburn FJ, Schneck DJ (1976) A constitutive equation for whole human blood. Biorheology 13: 201–210.
    [93] Zaman A, Ali N, Sajid M, et al. (2016) Numerical and analytical study of two-layered unsteady blood flow through catheterized artery. PLoS One 11: e0161377. doi: 10.1371/journal.pone.0161377
    [94] Secomb TW, Pries AR (2011) The microcirculation: physiology at the mesoscale. J Physiol 589: 1047–1052. doi: 10.1113/jphysiol.2010.201541
    [95] AlMomani TD, Vigmostad SC, Chivukula VK, et al. (2012) Red blood cell flow in the cardiovascular system: a fluid dynamics perspective. Crit Rev Biomed Eng 40: 427–440. doi: 10.1615/CritRevBiomedEng.v40.i5.30
    [96] Chung B, Johnson PC, Popel AS (2007) Application of Chimera grid to modelling cell motion and aggregation in a narrow tube. Int J Numer Meth Fluids 53: 105–128. doi: 10.1002/fld.1251
    [97] Ye T, Phan-Thien N, Lim CT (2016) Particle-based simulations of red blood cells-A review. J Biomech 49: 2255–2266. doi: 10.1016/j.jbiomech.2015.11.050
    [98] Tsubota K, Wada S (2010) Effect of the natural state of an elastic cellular membrane on tank-treading and tumbling motions of a single red blood cell. Phys Rev E Stat Nonlin Soft Matter Phys 81: 011910. doi: 10.1103/PhysRevE.81.011910
    [99] Nakamura M, Bessho S, Wada S (2013) Spring-network-based model of a red blood cell for simulating mesoscopic blood flow. Int J Numer Method Biomed Eng 29: 114–128. doi: 10.1002/cnm.2501
    [100] Hariprasad DS, Secomb TW (2014) Two-dimensional simulation of red blood cell motion near a wall under a lateral force. Phys Rev E Stat Nonlin Soft Matter Phys 90: 053014. doi: 10.1103/PhysRevE.90.053014
    [101] Pozrikidis C (2003) Numerical simulation of the flow-induced deformation of red blood cells. Ann Biomed Eng 31: 1194–1205. doi: 10.1114/1.1617985
    [102] Gounley J, Peng Y (2015) Computational modeling of membrane viscosity of red blood cells. Commun Comput Phys 17: 1073–1087. doi: 10.4208/cicp.2014.m355
    [103] Eggleton CD, Popel AS (1998) Large deformation of red blood cell ghosts in a simple shear flow. Phys Fluids 10: 1834–1845. doi: 10.1063/1.869703
    [104] Bagchi P (2007) Mesoscale simulation of blood flow in small vessels. Biophys J 92: 1858–1877. doi: 10.1529/biophysj.106.095042
    [105] Sun C, Munn LL (2005) Particulate nature of blood determines macroscopic rheology: A 2-D lattice Boltzmann analysis. Biophys J 88: 1635–1645. doi: 10.1529/biophysj.104.051151
    [106] Li X, Sarkar K (2008) Front tracking simulation of deformation and buckling instability of a liquid capsule enclosed by an elastic membrane. 227: 4998–5018.
    [107] Pan TW, Wang T (2009) Dynamical simulation of red blood cell rheology in microvessels. Int J Numer Anal Mod 6: 455–473.
    [108] Bessonov N, Crauste F, Fischer S, et al. (2011) Application of hybrid models to blood cell production in the bone marrow. Math Model of Nat Phenom 6: 2–12.
    [109] Volpert V, Bessonov N, Sequeira A, et al. (2015) Methods of blood flow modelling. Math Model Nat Phenom 11: 1–25.
    [110] Munn LL, Dupin MM (2008) Blood cell interactions and segregation in flow. Ann Biomed Eng 36: 534–544. doi: 10.1007/s10439-007-9429-0
    [111] Volpert V, Pujo-Menjouet L (2016) Blood cell dynamics: half of a century of modelling. Math Model Nat Phenom 11: 92–115. doi: 10.1051/mmnp/201611106
    [112] Imai Y, Omori T, Shimogonya Y, et al. (2016) Numerical methods for simulating blood flow at macro, micro, and multi scales. J Biomech 49: 2221–2228. doi: 10.1016/j.jbiomech.2015.11.047
    [113] Liu Y, Liu WK (2006) Rheology of red blood cell aggregation by computer simulation. J Comput Phys 220: 139–154. doi: 10.1016/j.jcp.2006.05.010
    [114] Liu Y, Zhang L, Wang X, et al. (2004) Coupling of Navier-Stokes equations with protein molecular dynamics and its application to hemodynamics. Int J Numer Meth Fluids 46: 1237–1252. doi: 10.1002/fld.798
    [115] Bagchi P, Johnson PC, Popel AS (2005) Computational fluid dynamic simulation of aggregation of deformable cells in a shear flow. J Biomech Eng 127: 1070–1080. doi: 10.1115/1.2112907
    [116] Le DV, White J, Peraire J, et al. (2009) An implicit immersed boundary method for three-dimensional fluid-membrane interactions. J Comput Phys 228: 8427–8445. doi: 10.1016/j.jcp.2009.08.018
    [117] Wu T, Feng JJ (2013) Simulation of malaria-infected red blood cells in microfluidic channels: Passage and blockage. Biomicrofluidics 7: 44115. doi: 10.1063/1.4817959
    [118] Ju M, Ye SS, Namgung B, et al. (2015) A review of numerical methods for red blood cell flow simulation. Comput Methods Biomech Biomed Engin 18: 130–140. doi: 10.1080/10255842.2013.783574
    [119] Shaw S, Ganguly S, Sibanda P, et al. (2014) Dispersion characteristics of blood during nanoparticle assisted drug delivery process through a permeable microvessel. Microvasc Res 92: 25–33. doi: 10.1016/j.mvr.2013.12.007
    [120] Shipley RJ, Chapman SJ (2010) Multiscale modelling of fluid and drug transport in vascular tumours. Bull Math Biol 72: 1464–1491. doi: 10.1007/s11538-010-9504-9
    [121] Pittman RN (2011) Regulation of Tissue Oxygenation. Morgan & Claypool Life Sciences.
    [122] Popel AS (1989) Theory of oxygen transport to tissue. Crit Rev Biomed Eng 17: 257–321.
    [123] Goldman D (2008) Theoretical models of microvascular oxygen transport to tissue. Microcirculation 15: 795–811. doi: 10.1080/10739680801938289
    [124] Krogh A (1919) The supply of oxygen to the tissues and the regulation of the capillary circulation. J Physiol 52: 457–474. doi: 10.1113/jphysiol.1919.sp001844
    [125] Buxton RB, Frank LR (1997) A model for the coupling between cerebral blood flow and oxygen metabolism during neural stimulation. J Cereb Blood Flow Metab 17: 64–72.
    [126] Buxton RB, Uludag K, Dubowitz DJ, et al. (2004) Modeling the hemodynamic response to brain activation. Neuroimage 23 (Suppl 1): S220–S233.
    [127] Secomb TW (2015) Krogh-cylinder and infinite-domain models for washout of an inert diffusible solute from tissue. Microcirculation 22: 91–98. doi: 10.1111/micc.12180
    [128] Schumacker PT, Samsel RW (1989) Analysis of oxygen delivery and uptake relationships in the Krogh tissue model. J Appl Physiol 67: 1234–1244.
    [129] Lagerlund TD, Low PA (1991) Axial diffusion and Michaelis-Menten kinetics in oxygen delivery in rat peripheral nerve. Am J Physiol 260: R430–R440.
    [130] Piiper J, Scheid P (1991) Diffusion limitation of O2 supply to tissue in homogeneous and heterogeneous models. Respir Physiol 85: 127–136. doi: 10.1016/0034-5687(91)90011-7
    [131] Whiteley JP, Gavaghan DJ, Hahn CE (2003) Mathematical modelling of pulmonary gas transport. J Math Biol 47: 79–99. doi: 10.1007/s00285-003-0196-8
    [132] McGuire BJ, Secomb TW (2001) A theoretical model for oxygen transport in skeletal muscle under conditions of high oxygen demand. J Appl Physiol 91: 2255–2265.
    [133] McGuire BJ, Secomb TW (2003) Estimation of capillary density in human skeletal muscle based on maximal oxygen consumption rates. Am J Physiol Heart Circ Physiol 285: H2382–H2391. doi: 10.1152/ajpheart.00559.2003
    [134] Federspiel WJ, Popel AS (1986) A theoretical analysis of the effect of the particulate nature of blood on oxygen release in capillaries. Microvasc Res 32: 164–189. doi: 10.1016/0026-2862(86)90052-X
    [135] Page TC, Light WR, Hellums JD (1998) Prediction of microcirculatory oxygen transport by erythrocyte/hemoglobin solution mixtures. Microvasc Res 56: 113–126. doi: 10.1006/mvre.1998.2088
    [136] Stathopoulos NA, Nair PK, Hellums JD (1987) Oxygen transport studies of normal and sickle red cell suspensions in artificial capillaries. Microvasc Res 34: 200–210. doi: 10.1016/0026-2862(87)90053-7
    [137] Severns ML, Adams JM (1982) The relation between Krogh and compartmental transport models. J Theor Biol 97: 239–249. doi: 10.1016/0022-5193(82)90101-1
    [138] Chen X, Buerk DG, Barbee KA, et al. (2007) A model of NO/O2 transport in capillary-perfused tissue containing an arteriole and venule pair. Ann Biomed Eng 35: 517–529. doi: 10.1007/s10439-006-9236-z
    [139] Ye GF, Moore TW, Buerk DG, et al. (1994) A compartmental model for oxygen-carbon dioxide coupled transport in the microcirculation. Ann Biomed Eng 22: 464–479. doi: 10.1007/BF02367083
    [140] Hayashi T, Watabe H, Kudomi N, et al. (2003) A theoretical model of oxygen delivery and metabolism for physiologic interpretation of quantitative cerebral blood flow and metabolic rate of oxygen. J Cereb Blood Flow Metab 23: 1314–1323.
    [141] Vadapalli A, Pittman RN, Popel AS (2000) Estimating oxygen transport resistance of the microvascular wall. Am J Physiol Heart Circ Physiol 279: H657–H671.
    [142] d'Angelo C (2007) Multiscale modelling of metabolism and transport phenomena in living tissues. Lausanne: PhD Thesis, EPFL.
    [143] Eggleton CD, Vadapalli A, Roy TK, et al. (2000) Calculations of intracapillary oxygen tension distributions in muscle. Math Biosci 167: 123–143. doi: 10.1016/S0025-5564(00)00038-9
    [144] Lucker A, Weber B, Jenny P (2015) A dynamic model of oxygen transport from capillaries to tissue with moving red blood cells. Am J Physiol Heart Circ Physiol 308: H206–H216. doi: 10.1152/ajpheart.00447.2014
    [145] Rees SE, Klaestrup E, Handy J, et al. (2010) Mathematical modelling of the acid-base chemistry and oxygenation of blood: a mass balance, mass action approach including plasma and red blood cells. Eur J Appl Physiol 108: 483–494.
    [146] Dash RK, Bassingthwaighte JB (2004) Blood HbO2 and HbCO2 dissociation curves at varied O2, CO2, pH, 2,3-DPG and temperature levels. Ann Biomed Eng 32: 1676–1693. doi: 10.1007/s10439-004-7821-6
    [147] Dash RK, Bassingthwaighte JB (2010) Erratum to: blood HbO2 and HbCO2 dissociation curves at varied O2, CO2, pH, 2,3-DPG and temperature levels. Ann Biomed Eng 38: 1683–1701. doi: 10.1007/s10439-010-9948-y
    [148] Schacterle RS, Adams JM, Ribando RJ (1991) A theoretical model of gas transport between arterioles and tissue. Microvasc Res 41: 210–228. doi: 10.1016/0026-2862(91)90023-5
    [149] Li Z, Yipintsoi T, Bassingthwaighte JB (1997) Nonlinear model for capillary-tissue oxygen transport and metabolism. Ann Biomed Eng 25: 604–619. doi: 10.1007/BF02684839
    [150] Vadapalli A, Goldman D, Popel AS (2002) Calculations of oxygen transport by red blood cells and hemoglobin solutions in capillaries. Artif Cells Blood Substit Immobil Biotechnol 30: 157–188. doi: 10.1081/BIO-120004338
    [151] Cabrera ME, Saidel GM, Kalhan SC (1998) Role of O2 in regulation of lactate dynamics during hypoxia: mathematical model and analysis. Ann Biomed Eng 26: 1–27. doi: 10.1114/1.28
    [152] Chen HS, Gross JF (1979) Estimation of tissue-to-plasma partition coefficients used in physiological pharmacokinetic models. J Pharmacokinet Biopharm 7: 117–125. doi: 10.1007/BF01059446
    [153] Sharan M, Popel AS (2002) A compartmental model for oxygen transport in brain microcirculation in the presence of blood substitutes. J Theor Biol 216: 479–500. doi: 10.1006/jtbi.2002.3001
    [154] Hsu R, Secomb TW (1989) A Green's function method for analysis of oxygen delivery to tissue by microvascular networks. Math Biosci 96: 61–78. doi: 10.1016/0025-5564(89)90083-7
    [155] Secomb TW, Hsu R, Park EY, et al. (2004) Green's function methods for analysis of oxygen delivery to tissue by microvascular networks. Ann Biomed Eng 32: 1519–1529. doi: 10.1114/B:ABME.0000049036.08817.44
    [156] Secomb TW (2016) A Green's function method for simulation of time-dependent solute transport and reaction in realistic microvascular geometries. Math Med Biol 33: 475–494. doi: 10.1093/imammb/dqv031
    [157] Lamkin-Kennard KA, Buerk DG, Jaron D (2004) Interactions between NO and O2 in the microcirculation: a mathematical analysis. Microvasc Res 68: 38–50. doi: 10.1016/j.mvr.2004.03.001
    [158] Sriram K, Vazquez BY, Yalcin O, et al. (2011) The effect of small changes in hematocrit on nitric oxide transport in arterioles. Antioxid Redox Signal 14: 175–185. doi: 10.1089/ars.2010.3266
    [159] Gundersen SI, Chen G, Palmer AF (2009) Mathematical model of NO and O2 transport in an arteriole facilitated by hemoglobin based O2 carriers. Biophys Chem 143: 1–17. doi: 10.1016/j.bpc.2009.02.005
    [160] Vullo V (2014) Circular Cylinders and Pressure Vessels. Springer.
    [161] Gonzalez-Fernandez JM, Ermentrout B (1994) On the origin and dynamics of the vasomotion of small arteries. Math Biosci 119: 127–167. doi: 10.1016/0025-5564(94)90074-4
    [162] Achakri H, Rachev A, Stergiopulos N, et al. (1994) A theoretical investigation of low frequency diameter oscillations of muscular arteries. Ann Biomed Eng 22: 253–263. doi: 10.1007/BF02368232
    [163] Carlson BE, Arciero JC, Secomb TW (2008) Theoretical model of blood flow autoregulation: roles of myogenic, shear-dependent, and metabolic responses. Am J Physiol Heart Circ Physiol 295: H1572–H1579. doi: 10.1152/ajpheart.00262.2008
    [164] Ford Versypt AN, Makrides E, Arciero JC, et al. (2015) Bifurcation study of blood flow control in the kidney. Math Biosci 263: 169–179. doi: 10.1016/j.mbs.2015.02.015
    [165] Carlson BE, Secomb TW (2005) A theoretical model for the myogenic response based on the length-tension characteristics of vascular smooth muscle. Microcirculation 12: 327–338. doi: 10.1080/10739680590934745
    [166] Arciero J, Harris A, Siesky B, et al. (2013) Theoretical analysis of vascular regulatory mechanisms contributing to retinal blood flow autoregulation. Invest Ophthalmol Vis Sci 54: 5584–5593. doi: 10.1167/iovs.12-11543
    [167] Ursino M, Di Giammarco P (1991) A mathematical model of the relationship between cerebral blood volume and intracranial pressure changes: the generation of plateau waves. Ann Biomed Eng 19: 15–42. doi: 10.1007/BF02368459
    [168] Ursino M (1988) A mathematical study of human intracranial hydrodynamics part 1-The cerebrospinal fluid pulse pressure. Ann Biomed Eng 16: 379–401. doi: 10.1007/BF02364625
    [169] Banaji M, Tachtsidis I, Delpy D, et al. (2005) A physiological model of cerebral blood flow control. Math Biosci 194: 125–173. doi: 10.1016/j.mbs.2004.10.005
    [170] Beard DA, Bassingthwaighte JB (2001) Modeling advection and diffusion of oxygen in complex vascular networks. Ann Biomed Eng 29: 298–310. doi: 10.1114/1.1359450
    [171] Jespersen SN, Ostergaard L (2012) The roles of cerebral blood flow, capillary transit time heterogeneity, and oxygen tension in brain oxygenation and metabolism. J Cereb Blood Flow Metab 32: 264–277. doi: 10.1038/jcbfm.2011.153
    [172] Lemon DD, Nair PK, Boland EJ, et al. (1987) Physiological factors affecting O2 transport by hemoglobin in an in vitro capillary system. J Appl Physiol 62: 798–806.
    [173] Goldman D, Popel AS (2000) A computational study of the effect of capillary network anastomoses and tortuosity on oxygen transport. J Theor Biol 206: 181–194. doi: 10.1006/jtbi.2000.2113
    [174] Tsoukias NM, Popel AS (2003) A model of nitric oxide capillary exchange. Microcirculation 10: 479–495. doi: 10.1038/sj.mn.7800210
    [175] Delp MD, Armstrong RB (1988) Blood flow in normal and denervated muscle during exercise in conscious rats. Am J Physiol 255: H1509–H1515.
    [176] Pohl U, De Wit C, Gloe T (2000) Large arterioles in the control of blood flow: role of endothelium-dependent dilation. Acta Physiol Scand 168: 505–510. doi: 10.1046/j.1365-201x.2000.00702.x
    [177] Segal SS (2005) Regulation of blood flow in the microcirculation. Microcirculation 12: 33–45. doi: 10.1080/10739680590895028
    [178] Harris A, Jonescu-Cuypers CP, Kagemann L, et al. (2003) Atlas of ocular blood flow: vascular anatomy, pathophysiology, and metabolism. Butterworth-Heinemann/Elsevier.
    [179] Aletti M, Gerbeau JF, Lombardi D (2016) A simplified fluid-structure model for arterial flow. Application to retinal hemodynamics. Comput Methods Appl Mech Engrg 306: 77–94. doi: 10.1016/j.cma.2016.03.044
    [180] Causin P, Malgaroli F (2016) Mathematical modeling of local perfusion in large distensible microvascular networks. arXiv:161002292 [q-bioTO].
    [181] Carichino L, Guidoboni G, Siesky B, et al. (2013) Effect of intraocular pressure and cerebrospinal fluid pressure on the blood flow in the central retinal vessels, In: Causin P, Guidoboni G, Sacco R et al., Editors, Integrated Multidisciplinary Approaches in the Study and Care of the Human Eye, Kugler Publications, 59–66.
    [182] Cassani S, Arciero J, Guidoboni G, et al. (2016) Theoretical predictions of metabolic flow regulation in the retina. J Model Ophthalmol 1.
    [183] Yu DY, Cringle SJ (2001) Oxygen distribution and consumption within the retina in vascularised and avascular retinas and in animal models of retinal disease. Prog Retin Eye Res 20: 175–208. doi: 10.1016/S1350-9462(00)00027-6
    [184] Yu DY, Cringle SJ, Su EN (2005) Intraretinal oxygen distribution in the monkey retina and the response to systemic hyperoxia. Invest Ophthalmol Vis Sci 46: 4728–4733. doi: 10.1167/iovs.05-0694
    [185] Roos MW (2004) Theoretical estimation of retinal oxygenation during retinal artery occlusion. Physiol Meas 25: 1523–1532. doi: 10.1088/0967-3334/25/6/016
    [186] Haugh LM, Linsenmeier RA, Goldstick TK (1990) Mathematical models of the spatial distribution of retinal oxygen tension and consumption, including changes upon illumination. Ann Biomed Eng 18: 19–36. doi: 10.1007/BF02368415
    [187] Cornelissen AJ, Dankelman J, VanBavel E, et al. (2002) Balance between myogenic, flow-dependent, and metabolic flow control in coronary arterial tree: a model study. Am J Physiol Heart Circ Physiol 282: H2224–H2237. doi: 10.1152/ajpheart.00491.2001
    [188] He Z, Lim JK, Nguyen CT, et al. (2013) Coupling blood flow and neural function in the retina: a model for homeostatic responses to ocular perfusion pressure challenge. Physiol Rep 1: e00055.
    [189] Ursino M (1991) Mechanisms of cerebral blood flow regulation. Crit Rev Biomed Eng 18: 255–288.
    [190] Ursino M (1991) A mathematical model of overall cerebral blood flow regulation in the rat. IEEE Trans Biomed Eng 38: 795–807. doi: 10.1109/10.83592
    [191] Nobrega AC, O'Leary D, Silva BM, et al. (2014) Neural regulation of cardiovascular response to exercise: role of central command and peripheral afferents. Biomed Res Int 2014: 478965.
    [192] Sigal IA, Ethier CR (2009) Biomechanics of the optic nerve head. Exp Eye Res 88: 799–807. doi: 10.1016/j.exer.2009.02.003
    [193] Downs JC, Roberts MD, Burgoyne CF, et al. (2009) Multiscale finite element modeling of the lamina cribrosa microarchitecture in the eye. Conf Proc IEEE Eng Med Biol Soc 2009: 4277–4280.
    [194] Grytz R, Girkin CA, Libertiaux V, et al. (2012) Perspectives on biomechanical growth and remodeling mechanisms in glaucoma. Mech Res Commun 42: 92–106. doi: 10.1016/j.mechrescom.2012.01.007
    [195] Newson T, El-Sheikh A (2006) Mathematical modeling of the biomechanics of the lamina cribrosa under elevated intraocular pressures. J Biomech Eng 128: 496–504. doi: 10.1115/1.2205372
    [196] Harris A, Guidoboni G, Arciero JC, et al. (2013) Ocular hemodynamics and glaucoma: the role of mathematical modeling. Eur J Ophthalmol 23: 139–146.
    [197] Guidoboni G, Harris A, Carichino L, et al. (2014) Effect of intraocular pressure on the hemodynamics of the central retinal artery: a mathematical model. Math Biosci Eng 11: 523–546. doi: 10.3934/mbe.2014.11.523
    [198] Cassani S, Harris A, Siesky B, et al. (2015) Theoretical analysis of the relationship between changes in retinal blood flow and ocular perfusion pressure. J Coupled Syst Multiscale Dyn 3: 38–46. doi: 10.1166/jcsmd.2015.1063
    [199] Causin P, Guidoboni G, Harris A, et al. (2014) A poroelastic model for the perfusion of the lamina cribrosa in the optic nerve head. Math Biosci 257: 33–41. doi: 10.1016/j.mbs.2014.08.002
    [200] Biesecker KR, Srienc AI, Shimoda AM, et al. (2016) Glial cell Calcium signaling mediates capillary regulation of blood flow in the retina. J Neurosci 36: 9435–9445. doi: 10.1523/JNEUROSCI.1782-16.2016
  • This article has been cited by:

    1. Jairam Choudhary, Surender Singh, Lata Nain, Thermotolerant fermenting yeasts for simultaneous saccharification fermentation of lignocellulosic biomass, 2016, 21, 07173458, 82, 10.1016/j.ejbt.2016.02.007
    2. Maria Priscila Lacerda, Eun Joong Oh, Carrie Eckert, The Model System Saccharomyces cerevisiae Versus Emerging Non-Model Yeasts for the Production of Biofuels, 2020, 10, 2075-1729, 299, 10.3390/life10110299
    3. Shweta Priya, Rameshwar Tiwari, Sarika Rana, Mohanram Saritha, Surender Singh, Anju Arora, Lata Nain, Saccharification of biopretreated paddy straw with indigenous holocellulase and fermentation with Saccharomyces cerevisiae LN1 under optimized conditions, 2016, 1, 2363-7692, 419, 10.1007/s40974-016-0021-z
    4. Rameshwar Avchar, Marizeth Groenewald, Abhishek Baghela, Wickerhamiella shivajii sp. nov., a thermotolerant yeast isolated from distillery effluent, 2019, 69, 1466-5026, 3262, 10.1099/ijsem.0.003616
    5. Vinay Sharma, A Review on Current Technological Advancement of Lignocellulosic Bioethanol Production, 2016, 1, 25728466, 10.15406/jabb.2016.01.00011
    6. E. Raja Sathendra, R. Praveenkumar, Baskar Gurunathan, S. Chozhavendhan, Mani Jayakumar, 2022, 9780323852692, 87, 10.1016/B978-0-323-85269-2.00016-2
    7. Yoshinori Murata, Charles O. Nwuche, Julius E. Nweze, Ifeanyi A. Ndubuisi, James C. Ogbonna, Potentials of multi-stress tolerant yeasts, Saccharomyces cerevisiae and Pichia kudriavzevii for fuel ethanol production from industrial cassava wastes, 2021, 111, 13595113, 305, 10.1016/j.procbio.2021.11.014
    8. K H A Rahman, K S K Ismail, N Najimudin, Growth of thermotolerant Pichia kudriavzevii UniMAP 3-1 strain for ethanol production using xylose and glucose at different fermentation temperatures, 2021, 765, 1755-1307, 012107, 10.1088/1755-1315/765/1/012107
    9. Estifanos Hawaz, Mesfin Tafesse, Anteneh Tesfaye, Dereje Beyene, Solomon Kiros, Gessese Kebede, Teun Boekhout, Bart Theelen, Marizeth Groenewald, Ayantu Degefe, Sissay Degu, Alene Admas, Diriba Muleta, Isolation and characterization of bioethanol producing wild yeasts from bio-wastes and co-products of sugar factories, 2022, 72, 1590-4261, 10.1186/s13213-022-01695-3
    10. Teun Boekhout, Isolation and Characterization of Bioethanol Producing Wild Yeasts from Bio-Wastes and Co-Products of Sugar Factories, 2022, 1556-5068, 10.2139/ssrn.4048132
    11. N. Kusumawati, S. H. Sumarlan, E. Zubaidah, A. K. Wardani, Isolation of xylose-utilizing yeasts from oil palm waste for xylitol and ethanol production, 2023, 10, 2197-4365, 10.1186/s40643-023-00691-y
    12. Pratima Bajpai, 2024, 9780443214929, 133, 10.1016/B978-0-443-21492-9.00006-9
    13. Islem Dammak, Nourelhouda Abdelkefi, Imen Ben Atitallah, Michael Brysch-Herzberg, Abdulrahman H. Alessa, Salma Lasram, Hela Zouari-Mechichi, Tahar Mechichi, Characterization and biocontrol potential of Wickerhamomyces subpelliculosus yeasts isolated from dates: volatile compounds-mediated antifungal activity against mycotoxigenic Penicillium strains, 2024, 24058440, e39504, 10.1016/j.heliyon.2024.e39504
  • 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(11300) PDF downloads(3247) Cited by(22)

Article outline

Figures and Tables

Figures(3)  /  Tables(5)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog