Friction losses of Newtonian and non-Newtonian fluids flowing in laminar regime in a helical coil

This study aimed to carry out experimental work to determine, for Newtonian and non-Newtonian fluids, the friction factor (fc) with simultaneous heat transfer, at constant wall temperature as boundary condition, in fully developed laminar flow inside a vertical helical coil. The Newtonian fluids studied were aqueous solutions of glycerol, 25%, 36%, 43%, 59% and 78% (w/w). The non-Newtonian fluids were aqueous solutions of carboxymethylcellulose (CMC), a polymer, with concentrations of 0.2%, 0.3%, 0.4% and 0.6% (w/w) and aqueous solutions of xanthan gum (XG), another polymer, with concentrations of 0.1% and 0.2% (w/w). According to the rheological study done, the polymer solutions had shear-thinning behavior and different values of viscoelasticity. The helical coil used has an internal diameter, curvature ratio, length and pitch, respectively: 0.00483 m, 0.0263, 5.0 m and 11.34 mm. It was concluded that the friction factors, with simultaneous heat transfer, for Newtonian fluids can be calculated using expressions from literature for isothermal flows. The friction factors for CMC and XG solutions are similar to those for Newtonian fluids when the Dean number, based in a generalized Reynolds number, is less than 80. For Dean numbers higher than 80, the friction factors of the CMC solutions are lower those of the XG solutions and of the Newtonian fluids. In this range the friction factors decrease with the increase of the viscometric component of the solution and increase for increasing elastic component. The change of behavior at Dean number 80, for Newtonian and non-Newtonian fluids, is in accordance with the study of Ali [4]. There is a change of behavior at Dean number 80, for Newtonian and non-Newtonian fluids, which is in according to previous studies. The data also showed that the use of the bulk temperature or of the film temperature to calculate the physical properties of the fluid has a residual effect in the friction factor values.

Heat transfer coefficients from Newtonian and non-Newtonian fluids flowing in laminar regime in a helical coil

This study aimed to carry out experimental work to obtain, for Newtonian and non-Newtonian fluids, heat transfer coefficients, at constant wall temperature as boundary condition, in fully developed laminar flow inside a helical coil. The Newtonian fluids studied were aqueous solutions of glycerol, 25%, 36%, 43%, 59% and 78% (w/w) and the non-Newtonian fluids aqueous solutions of carboxymethylcellulose (CMC), a polymer, with concentrations 0.1%, 0.2%, 0.3%, 0.4% and 0.6% (w/w) and aqueous solutions of xanthan gum (XG), another polymer, with concentrations 0.1% and 0.2% (w/w). According to the rheological study performed, the polymer solutions had shear thinning behavior and different values of elasticity. The helical coil used has internal diameter, curvature ratio, length and pitch, respectively: 0.004575 m, 0.0263, 5.0 m and 11.34 mm. The Nusselt numbers for the CMC solutions are, on average, slightly higher than those for Newtonian fluids, for identical Prandtl and generalized Dean numbers. As outcome, the viscous component of the shear thinning polymer tends to potentiate the mixing effect of the Dean cells. The Nusselt numbers of the XG solutions are significant lower than those of the Newtonian solutions, for identical Prandtl and generalized Dean numbers. Therefore, the elastic component of the polymer tends to diminish the mixing effect of the Dean cells. A global correlation, for Nusselt number as a function of Péclet, generalized Dean and Weissenberg numbers for all Newtonian and non-Newtonian solutions studied, is presented.

Convex optimization framework for intermediate deadline assignment in soft and hard real-time distributed systems

It is generally challenging to determine end-to-end delays of applications for maximizing the aggregate system utility subject to timing constraints. Many practical approaches suggest the use of intermediate deadline of tasks in order to control and upper-bound their end-to-end delays. This paper proposes a unified framework for different time-sensitive, global optimization problems, and solves them in a distributed manner using Lagrangian duality. The framework uses global viewpoints to assign intermediate deadlines, taking resource contention among tasks into consideration. For soft real-time tasks, the proposed framework effectively addresses the deadline assignment problem while maximizing the aggregate quality of service. For hard real-time tasks, we show that existing heuristic solutions to the deadline assignment problem can be incorporated into the proposed framework, enriching their mathematical interpretation.

A survey of technologies and applications for climbing robots locomotion and adhesion

The interest in the development of climbing robots has grown rapidly in the last years. Climbing robots are useful devices that can be adopted in a variety of applications, such as maintenance and inspection in the process and construction industries. These systems are mainly adopted in places where direct access by a human operator is very expensive, because of the need for scaffolding, or very dangerous, due to the presence of an hostile environment. The main motivations are to increase the operation efficiency, by eliminating the costly assembly of scaffolding, or to protect human health and safety in hazardous tasks. Several climbing robots have already been developed, and other are under development, for applications ranging from cleaning to inspection of difficult to reach constructions. A wall climbing robot should not only be light, but also have large payload, so that it may reduce excessive adhesion forces and carry instrumentations during navigation. These machines should be capable of travelling over different types of surfaces, with different inclinations, such as floors, walls, or ceilings, and to walk between such surfaces (Elliot et al. (2006); Sattar et al. (2002)). Furthermore, they should be able of adapting and reconfiguring for various environment conditions and to be self-contained. Up to now, considerable research was devoted to these machines and various types of experimental models were already proposed (according to Chen et al. (2006), over 200 prototypes aimed at such applications had been developed in the world by the year 2006). However, we have to notice that the application of climbing robots is still limited. Apart from a couple successful industrialized products, most are only prototypes and few of them can be found in common use due to unsatisfactory performance in on-site tests (regarding aspects such as their speed, cost and reliability). Chen et al. (2006) present the main design problems affecting the system performance of climbing robots and also suggest solutions to these problems. The major two issues in the design of wall climbing robots are their locomotion and adhesion methods. With respect to the locomotion type, four types are often considered: the crawler, the wheeled, the legged and the propulsion robots. Although the crawler type is able to move relatively faster, it is not adequate to be applied in rough environments. On the other hand, the legged type easily copes with obstacles found in the environment, whereas generally its speed is lower and requires complex control systems. Regarding the adhesion to the surface, the robots should be able to produce a secure gripping force using a light-weight mechanism. The adhesion method is generally classified into four groups: suction force, magnetic, gripping to the surface and thrust force type. Nevertheless, recently new methods for assuring the adhesion, based in biological findings, were proposed. The vacuum type principle is light and easy to control though it presents the problem of supplying compressed air. An alternative, with costs in terms of weight, is the adoption of a vacuum pump. The magnetic type principle implies heavy actuators and is used only for ferromagnetic surfaces. The thrust force type robots make use of the forces developed by thrusters to adhere to the surfaces, but are used in very restricted and specific applications. Bearing these facts in mind, this chapter presents a survey of different applications and technologies adopted for the implementation of climbing robots locomotion and adhesion to surfaces, focusing on the new technologies that are recently being developed to fulfill these objectives. The chapter is organized as follows. Section two presents several applications of climbing robots. Sections three and four present the main locomotion principles, and the main "conventional" technologies for adhering to surfaces, respectively. Section five describes recent biological inspired technologies for robot adhesion to surfaces. Section six introduces several new architectures for climbing robots. Finally, section seven outlines the main conclusions.

The Neogene of the Lower Tagus Basin (Portugal)

Revista Española de Paleontologia 19 (2), 229-242

Electricity markets modeling considering complex contracts and aggregators

All over the world, the liberalization of electricity markets, which follows different paradigms, has created new challenges for those involved in this sector. In order to respond to these challenges, electric power systems suffered a significant restructuring in its mode of operation and planning. This restructuring resulted in a considerable increase of the electric sector competitiveness. Particularly, the Ancillary Services (AS) market has been target of constant renovations in its operation mode as it is a targeted market for the trading of services, which have as main objective to ensure the operation of electric power systems with appropriate levels of stability, safety, quality, equity and competitiveness. In this way, with the increasing penetration of distributed energy resources including distributed generation, demand response, storage units and electric vehicles, it is essential to develop new smarter and hierarchical methods of operation of electric power systems. As these resources are mostly connected to the distribution network, it is important to consider the introduction of this kind of resources in AS delivery in order to achieve greater reliability and cost efficiency of electrical power systems operation. The main contribution of this work is the design and development of mechanisms and methodologies of AS market and for energy and AS joint market, considering different management entities of transmission and distribution networks. Several models developed in this work consider the most common AS in the liberalized market environment: Regulation Down; Regulation Up; Spinning Reserve and Non-Spinning Reserve. The presented models consider different rules and ways of operation, such as the division of market by network areas, which allows the congestion management of interconnections between areas; or the ancillary service cascading process, which allows the replacement of AS of superior quality by lower quality of AS, ensuring a better economic performance of the market. A major contribution of this work is the development an innovative methodology of market clearing process to be used in the energy and AS joint market, able to ensure viable and feasible solutions in markets, where there are technical constraints in the transmission network involving its division into areas or regions. The proposed method is based on the determination of Bialek topological factors and considers the contribution of the dispatch for all services of increase of generation (energy, Regulation Up, Spinning and Non-Spinning reserves) in network congestion. The use of Bialek factors in each iteration of the proposed methodology allows limiting the bids in the market while ensuring that the solution is feasible in any context of system operation. Another important contribution of this work is the model of the contribution of distributed energy resources in the ancillary services. In this way, a Virtual Power Player (VPP) is considered in order to aggregate, manage and interact with distributed energy resources. The VPP manages all the agents aggregated, being able to supply AS to the system operator, with the main purpose of participation in electricity market. In order to ensure their participation in the AS, the VPP should have a set of contracts with the agents that include a set of diversified and adapted rules to each kind of distributed resource. All methodologies developed and implemented in this work have been integrated into the MASCEM simulator, which is a simulator based on a multi-agent system that allows to study complex operation of electricity markets. In this way, the developed methodologies allow the simulator to cover more operation contexts of the present and future of the electricity market. In this way, this dissertation offers a huge contribution to the AS market simulation, based on models and mechanisms currently used in several real markets, as well as the introduction of innovative methodologies of market clearing process on the energy and AS joint market. This dissertation presents five case studies; each one consists of multiple scenarios. The first case study illustrates the application of AS market simulation considering several bids of market players. The energy and ancillary services joint market simulation is exposed in the second case study. In the third case study it is developed a comparison between the simulation of the joint market methodology, in which the player bids to the ancillary services is considered by network areas and a reference methodology. The fourth case study presents the simulation of joint market methodology based on Bialek topological distribution factors applied to transmission network with 7 buses managed by a TSO. The last case study presents a joint market model simulation which considers the aggregation of small players to a VPP, as well as complex contracts related to these entities. The case study comprises a distribution network with 33 buses managed by VPP, which comprises several kinds of distributed resources, such as photovoltaic, CHP, fuel cells, wind turbines, biomass, small hydro, municipal solid waste, demand response, and storage units.

Perinatal Outcome and Change in Body Mass Index in Mothers of Dichorionic Twins: a Longitudinal Cohort Study

We used a prospective cohort to analyze the effect of change in BMI rather than change in weight, in mothers carrying dichorionic twins from a population that did not receive any dietary intervention. A total of 269 mothers (150 nulliparas and 119 multiparas) were evaluated. The average change (%) from the pre-gravid BMI was 7.2+/-6.1, 17.4+/-8.2, and 28.7+/-10.8, at 12-14, 22-25, and 30-34 weeks, respectively, without difference between nulliparas and multiparas. The comparison between maternities below or above the average change from the pregravid BMI failed to demonstrate an advantage (in terms of total twin birthweight and gestational age) of an above average change from the pregravid BMI, even when the lower versus upper quartiles were compared. Our observations reached different conclusions regarding the recommended universal dietary intervention in twin gestations. A cautious approach is advocated towards seemingly harmless excess weight gain, as normal weight women may turn overweight, or even obese, by the end of pregnancy, and be exposed to the untoward effects of obesity on future health and body image.

Analysis and treatment of a nineteenth century oil painting

This thesis focused on the study and treatment of a 19th century female portrait in oil from ECOMUSEU Municipal do Seixal, Portugal. The portrait, which depicts Isabel Maria Lourenço Affonso was in poor condition and a large strip of paint and canvas was missing (approximately 9cm by 66cm, almost 11% of the total surface area). The portrait is a companion piece to a male portrait (the relationship was established as part of this study), therefore a technical study of both paintings was considered essential to support the choices made during the treatment. The project involved three main areas: - The study of the history, condition, materials and techniques of both paintings. This allowed their comparison and understanding of their relationship; - The treatment of Isabel Maria Lourenço Affonso. The choices made and problems encountered are described. - The production of a replacement for the missing strip of paint and canvas. The practical solution developed to overcome such an unusual challenge is described along with the creative and flexible thinking required. Because not all traditional infill materials cope well on a mechanical level with thin layers over a very large surface (many are too brittle), strict criteria had to be employed to choose the appropriate material. The primary goal was to find a fill which would remain flexible and be capable of accepting surface texture, such that there would be a good visual match with the painting. Analysis and testing was carried out to evaluate the physical properties of the fill material chosen, BEVA® Gesso-P. The successful creation of the replacement strip has resulted in two publications and one presentation: Publication pending in The Picture Restorer, Leslie Carlyle, Raquel Marques, Isabel Pombo Cardoso and Sara Babo, “Creating a Textured Replacement Strip for the Missing Lower Portion of an Oil Portrait: Problem Solving and Practical Solutions”. Abstract accepted for presentation and publication, International Meeting on Retouching of Cultural Heritage (2RECH), Raquel Marques, Leslie Carlyle and Isabel Pombo Cardoso, “Textured Replacement Strip for a Missing Portion of a Portrait: Problem Solving and Practical Solutions”.

Population dynamics of 0-group flounder (Platichthys flesus L.) in Galway Bay, West of Ireland and the value of flounder from shore angling tourism in Ireland

The economic value of flounder from shore angling around Ireland was assessed. Flounder catches from shore angling tournaments around Ireland were related to domestic and overseas shore angling expenditure in order to determine an economic value for the species. Temporal trends in flounder angling catches, and specimen (trophy) flounder reports were also investigated. Flounder was found to be the most caught shore angling species in competitions around Ireland constituting roughly one third of the shore angling competition catch although this did vary by area. The total value of flounder from shore angling tourism was estimated to be of the order of €8.4 million. No significant temporal trends in flounder angling catches and specimen reports were found. Thus there is no evidence from the current study for any decline in flounder stocks. The population dynamics of 0-group flounder during the early benthic stage was investigated at estuarine sites in Galway Bay, west of Ireland. Information was analysed from the March to June sampling period over five years (2002 to 2006). Spatial and temporal variations in settlement and population length structure were analysed between beach and river habitats and sites. Settlement of flounder began from late March to early May of each year, most commonly in April. Peak settlement was usually in April or early May. Settlement was recorded earlier than elsewhere, although most commonly was similar to the southern part of the UK and northern France. Settlement was generally later in tidal rivers than on sandy beaches. Abundance of 0-group flounder in Galway Bay did not exhibit significant inter -annual variability. 0-group flounder were observed in dense aggregations of up to 105 m'2, which were patchy in distribution. Highest densities of 0-group flounder were recorded in limnetic and oligohaline areas as compared with the lower densities in polyhaline and to a lesser extent mesohaline areas. Measurements to of salinity allowed the classification of beaches, and tidal river sections near the mouth, into a salinity based scheme for length comparisons. Beaches were classified as polyhaline,the lower section of rivers as mesohaline, and the middle and upper sections as oligohaline. Over the March to June sampling period 0-group flounder utilised different sections at different length ranges and were significantly larger in more upstream sections. During initial settlement in April, 0-group flounder of 8-10 mm (standard length, SL) were present in abundance on polyhaline sandy beaches. By about 10mm (SL), flounder were present in all polyhaline, mesohaline and (oligohaline) sections. 0-group flounder became absent or in insignificant numbers in polyhaline and mesohaline sections in a matter of weeks after first appearance. From April to June, 0-group flounder of 12-30mm (SL) were found in more upstream locations in the oligohaline sections. About one month (May or June) after initial settlement, 0-group flounder became absent from the oligohaline sections. Concurrently, flounder start to reappear in mesohaline and polyhaline areas at approximately 30mm (SL) in June. The results indicate 0-group flounder in the early benthic stage are associated with low salinity areas, but as they grow, this association diminishes. Results strongly suggest that migration of 0-group flounder between habitats takes place during the early benthic phase.

Relay-assisted transmission and radio resource management for wireless networks

The following paper presents an overview of the Ph.D Thesis1 presented in [1], which compiles all the research done during the period of time between 2004-2007. In that dissertation the relay-assisted transmission with half-duplex relays is analyzed from different points of view. This study is motivated by the necessity of finding innovative solutions to cope with the requirements of next generation wireless services, and with current radio technology. The use of relayed communications represents a change of paradigm of conventional communications, and requires the definition and evaluation of protocols to be applied to single or multiple-user relay communication. With the two fold goal of enhancing spectral efficiency and homogenize service in cellular communications, system design is investigated at physical (type of transmissions of the relay, decoding mode, ..) and upper layers (resource allocation, dynamic link control).

Capsule permeability via polymer and protein ingress/egress.

Static incubation tests, where microcapsules and beads are contacted with polymer and protein solutions, have been developed for the characterization of permselective materials applied for bioartificial organs and drug delivery. A combination of polymer ingress, detected by size-exclusion chromatography, and protein ingress/ egress, assessed by gel electrophoresis, provides information regarding the diffusion kinetics, molar mass cutoff(MMCO) and permeability. This represents an improvement over existing permeability measurements that are based on the diffusion of a single type of solute. Specifically, the permeability of capsules based on alginate, cellulose sulfate, polymethylene-co-guanidine were characterized as a function of membrane thickness. Solid alginate beads were also evaluated. The MMCO of these capsules was estimated to be between 80 and 90 kDa using polymers, and between 116-150 kDa with proteins. Apparently, the globular shape of the proteins (radius of gyration (Rg) of 4.2-4.6 nm) facilitates their passage through the membrane, comparatively to the polysaccharide coil conformation (Rg of 6.5-8.3 nm). An increase of the capsule membrane thickness reduced these values. The MMCO of the beads, which do not have a membrane limiting their permselective properties, was higher, between 110 and 200 kDa with dextrans, and between 150 and 220 kDa with proteins. Therefore, although the permeability estimated with biologically relevant molecules is generally higher due to their lower radius of gyration, both the MMCO of synthetic and natural watersoluble polymers correlate well, and can be used as in vitro metrics for the immune protection ability of microcapsules and microbeads. This article shows, to the authors' knowledge, the first reported concordance between permeability measures based on model natural and biological macromolecules.

Nava enhances ventilatory variability and diaphragmatic activity/tidal volume coupling

INTRODUCTION. Neurally Adjusted Ventilatory Assist (NAVA) is a new ventilatory mode in which ventilator settings are adjusted based on the electrical activity detected in the diaphragm (Eadi). This mode offers significant advantages in mechanical ventilation over standard pressure support (PS) modes, since ventilator input is determined directly from patient ventilatory demand. Therefore, it is expected that tidal volume (Vt) under NAVA would show better correlation with Eadi compared with PS, and exhibit greater variability due to the variability in the Eadi input to the ventilator. OBJECTIVES. To compare tidal volume variability in PS and NAVA ventilation modes, and its correlation with patient ventilatory demand (as characterized by maximum Eadi). METHODS. Acomparative study of patient-ventilator interaction was performed for 22 patients during standard PS with clinician determined ventilator settings; and NAVA, with NAVA gain set to ensure the same peak airway pressure as the total pressure obtained in PS. A 20 min continuous recording was performed in each ventilator mode. Respiratory rate, Vt, and Eadi were recorded. Tidal volume variance and Pearson correlation coefficient between Vt and Eadi were calculated for each patient. A periodogram was plotted for each ventilator mode and each patient, showing spectral power as a function of frequency to assess variability. RESULTS. Median, lower quartile and upper quartile values for Vt variance and Vt/Eadi correlation are shown in Table 1. The NAVA cohort exhibits substantially greater correlation and variance than the PS cohort. Power spectrums for Vt and Eadi are shown in Fig. 1 (PS and NAVA) for a typical patient. The enlarged section highlights how changes in Eadi are highly synchronized with NAVA ventilation, but less so for PS. CONCLUSIONS. There is greater variability in tidal volume and correlation between tidal volume and diaphragmatic electrical activity with NAVA compared to PS. These results are consistent with the improved patient-ventilator synchrony reported in the literature.

A meta-analysis of the association of fracture risk and body mass index in women.

Several recent studies suggest that obesity may be a risk factor for fracture. The aim of this study was to investigate the association between body mass index (BMI) and future fracture risk at different skeletal sites. In prospective cohorts from more than 25 countries, baseline data on BMI were available in 398,610 women with an average age of 63 (range, 20-105) years and follow up of 2.2 million person-years during which 30,280 osteoporotic fractures (6457 hip fractures) occurred. Femoral neck BMD was measured in 108,267 of these women. Obesity (BMI ≥ 30 kg/m(2) ) was present in 22%. A majority of osteoporotic fractures (81%) and hip fractures (87%) arose in non-obese women. Compared to a BMI of 25 kg/m(2) , the hazard ratio (HR) for osteoporotic fracture at a BMI of 35 kg/m(2) was 0.87 (95% confidence interval [CI], 0.85-0.90). When adjusted for bone mineral density (BMD), however, the same comparison showed that the HR for osteoporotic fracture was increased (HR, 1.16; 95% CI, 1.09-1.23). Low BMI is a risk factor for hip and all osteoporotic fracture, but is a protective factor for lower leg fracture, whereas high BMI is a risk factor for upper arm (humerus and elbow) fracture. When adjusted for BMD, low BMI remained a risk factor for hip fracture but was protective for osteoporotic fracture, tibia and fibula fracture, distal forearm fracture, and upper arm fracture. When adjusted for BMD, high BMI remained a risk factor for upper arm fracture but was also a risk factor for all osteoporotic fractures. The association between BMI and fracture risk is complex, differs across skeletal sites, and is modified by the interaction between BMI and BMD. At a population level, high BMI remains a protective factor for most sites of fragility fracture. The contribution of increasing population rates of obesity to apparent decreases in fracture rates should be explored. © 2014 American Society for Bone and Mineral Research.

Multi-stage garnet in the internal Brianconnais basement (Ambin massif, Savoy): New petrological constraints on the blueschist-facies metamorphism in the Western Alps and tectonic implications

Three types of garnet have been distinguished in pelitic schists from an epidote-blueschist-facies unit of the Ambin and South Vanoise Brianconnais massifs on the basis of texture, chemical zoning and mineral inclusion characterization. Type-1 garnet cores with high Mn/Ca ratios are interpreted as pre-Alpine relicts, whereas Type-1 garnet rims, Type-2 inclusion-rich porphyroblasts and smaller Type-3 garnets are Alpine. The latter are all characterized by low Mn/Ca ratios and a coexisting mineral assemblage of blue amphibole, high-Si phengite, epidote and quartz. Prograde growth conditions during Alpine D-1 high-pressure (HP) metamorphism are recorded by a decrease in Mn and increase in Fe (+/-Ca) in the Type-2 garnets, culminating in peak P-T conditions of 14-16 kbar and 500degreesC in the deepest parts of the Ambin dome. The multistage growth history of Type-1 garnets indicates a polymetamorphic history for the Ambin and South Vanoise massifs; unfortunately, no age constraints are available. The new metamorphic constraints on the Alpine event in the massifs define a metamorphic T `gap' between them and their surrounding cover (Brianconnais and upper Schistes Lustres units), which experienced metamorphism only in the stability field of carpholite-lawsonite (T < 400degreesC). These data and supporting structural studies confirm that the Ambin and South Vanoise massifs are slices of `eclogitized' continental crust tectonically extruded within the Schistes Lustres units and Brianconnais covers. The corresponding tectonic contacts with top-to-east movement are responsible for the juxtaposition of lower-grade metamorphic units on the Ambin and South Vanoise massifs.

Phosphate forms in plant and their internal buffering in five soybean cultivars

Differences among plants in their ability to support nutritional stress periods may be caused by a differential vacuole capacity of ion storage and release and may also depend on the intensity of nutrient re-translocation under such conditions. In five soybean cultivars, submitted to eight days of P deprivation, the dry matter production and the contents of three phosphorus (P) forms - inorganic (Pi), organic (Po), and acid-soluble total (Pts) of different plant organs were determined. Pi release velocity (RSPi) was estimated as the tangent to the equations obtained for Pi f(t) at the point t = 2 days (the mean point in the period of greatest Pi decrease), considering that -deltaPi/deltat expresses the rate of Pi release. The internal Pi buffering capacity (IBCPi) was calculated as the inverse of the RSPi. Cultivars' differences in size of the non-metabolic Pi pool, RSPi, and the ability to transport Pi from less to more actively metabolizing regions were evaluated. The preferential Pi source and sink compartments under limited P absorption conditions were also evaluated. The cultivar Santa Rosa showed the highest Pi storage ability when the external supply was high, and a more intensive release under low P supply conditions than IAC8 and UFV1. The cultivar Uberaba was superior to Doko in its ability to store and use Pi. In all cultivars, upper leaves and roots were the main sink of Pi stored in the middle and lower leaves. Roots and upper leaves showed larger RSPi and lower IBCPi values than middle and lower leaves.