Quantification of the effects of architectural traits on dry mass production and light interception of tomato canopy under different temperature regimes using a dynamic functional-structural plant model

There is increasing interest in evaluating the environmental effects on crop architectural traits and yield improvement. However, crop models describing the dynamic changes in canopy structure with environmental conditions and the complex interactions between canopy structure, light interception, and dry mass production are only gradually emerging. Using tomato (Solanum lycopersicum L.) as a model crop, a dynamic functional-structural plant model (FSPM) was constructed, parameterized, and evaluated to analyse the effects of temperature on architectural traits, which strongly influence canopy light interception and shoot dry mass. The FSPM predicted the organ growth, organ size, and shoot dry mass over time with high accuracy (>85%). Analyses of this FSPM showed that, in comparison with the reference canopy, shoot dry mass may be affected by leaf angle by as much as 20%, leaf curvature by up to 7%, the leaf length: width ratio by up to 5%, internode length by up to 9%, and curvature ratios and leaf arrangement by up to 6%. Tomato canopies at low temperature had higher canopy density and were more clumped due to higher leaf area and shorter internodes. Interestingly, dry mass production and light interception of the clumped canopy were more sensitive to changes in architectural traits. The complex interactions between architectural traits, canopy light interception, dry mass production, and environmental conditions can be studied by the dynamic FSPM, which may serve as a tool for designing a canopy structure which is 'ideal' in a given environment.

Drag reduction phenomenon in viscous oil-water dispersed pipe flow: Experimental investigation and phenomenological modeling

An experimental study on drag-reduction phenomenon in dispersed oil-water flow has been performed in a 26-mm-i.d. Twelve meter long horizontal glass pipe. The flow was characterized using a novel wire-mesh sensor based on capacitance measurements and high-speed video recording. New two-phase pressure gradient, volume fraction, and phase distribution data have been used in the analysis. Drag reduction and slip ratio were detected at oil volume fractions between 10 and 45% and high mixture Reynolds numbers, and with water as the dominant phase. Phase-fraction distribution diagrams and cross-sectional imaging of the flow suggested the presence of a higher amount of water near to the pipe wall. Based on that, a phenomenology for explaining drag reduction in dispersed flow in a flow situation where slip ratio is significant is proposed. A simple phenomenological model is developed and the agreement between model predictions and data, including data from the literature, is encouraging. (c) 2011 American Institute of Chemical Engineers AIChE J, 2012

Modelagem numérica de escoamento em padrão anular com transferência de calor e massa

Annular flow is the prevailing pattern in transport and energy conversion systems and therefore, one of the most important patterns in multiphase flow in ducts. The correct prediction of the pressure gradient and heat transfer coefficient is essential for optimizing the system s capacity. The objective of this work is to develop and implement a numerical algorithm capable of predicting hydrodynamic and thermal characteristics for upflow, vertical, annular flow. The numerical algorithm is then complemented with the physical modeling of phenomena that occurs in this flow pattern. These are, turbulence, entrainment and deposition and phase change. For the development of the numerical model, axial diffusion of heat and momentum is neglected. In this way the time-averaged equations are solved in their parabolic form obtaining the velocity and temperature profiles for each axial step at a time, together with the global parameters, namely, pressure gradient, mean film thickness and heat transfer coefficient, as well as their variation in the axial direction. The model is validated for the following conditions: fully-developed laminar flow with no entrainment; fully developed laminar flow with heat transfer, fully-developed turbulent flow with entrained drops, developing turbulent annular flow with entrained drops, and turbulent flow with heat transfer and phase change

Analytical solution for transient one-dimensional couette flow considering constant and time-dependent pressure gradients

This paperaims to determine the velocity profile, in transient state, for a parallel incompressible flow known as Couette flow. The Navier-Stokes equations were applied upon this flow. Analytical solutions, based in Fourier series and integral transforms, were obtained for the one-dimensional transient Couette flow, taking into account constant and time-dependent pressure gradients acting on the fluid since the same instant when the plate starts it´s movement. Taking advantage of the orthogonality and superposition properties solutions were foundfor both considered cases. Considering a time-dependent pressure gradient, it was found a general solution for the Couette flow for a particular time function. It was found that the solution for a time-dependent pressure gradient includes the solutions for a zero pressure gradient and for a constant pressure gradient.

Amazon Rainforest Exchange of Carbon and Subcanopy Air Flow: Manaus LBA Site-A Complex Terrain Condition

On the moderately complex terrain covered by dense tropical Amazon Rainforest (Reserva Biologica do Cuieiras-ZF2-02 degrees 36'17.1 '' S, 60 degrees 12'24.4 '' W), subcanopy horizontal and vertical gradients of the air temperature, CO2 concentration and wind field were measured for the dry and wet periods in 2006. We tested the hypothesis that horizontal drainage flow over this study area is significant and can affect the interpretation of the high carbon uptake rates reported by previous works at this site. A similar experimental design as the one by Tota et al. (2008) was used with a network of wind, air temperature, and CO2 sensors above and below the forest canopy. A persistent and systematic subcanopy nighttime upslope (positive buoyancy) and daytime downslope (negative buoyancy) flow pattern on a moderately inclined slope (12%) was observed. The microcirculations observed above the canopy (38 m) over the sloping area during nighttime presents a downward motion indicating vertical convergence and correspondent horizontal divergence toward the valley area. During the daytime an inverse pattern was observed. The microcirculations above the canopy were driven mainly by buoyancy balancing the pressure gradient forces. In the subcanopy space the microcirculations were also driven by the same physical mechanisms but probably with the stress forcing contribution. The results also indicated that the horizontal and vertical scalar gradients (e. g., CO2) were modulated by these micro-circulations above and below the canopy, suggesting that estimates of advection using previous experimental approaches are not appropriate due to the tridimensional nature of the vertical and horizontal transport locally. This work also indicates that carbon budget from tower-based measurement is not enough to close the system, and one needs to include horizontal and vertical advection transport of CO2 into those estimates.

Volume averaging of multiphase flows with hydrate formation in subsea pipelines

In oil and gas pipeline operations, the gas, oil, and water phases simultaneously move through pipe systems. The mixture cools as it flows through subsea pipelines, and forms a hydrate formation region, where the hydrate crystals start to grow and may eventually block the pipeline. The potential of pipe blockage due to hydrate formation is one of the most significant flow-assurance problems in deep-water subsea operations. Due to the catastrophic safety and economic implications of hydrate blockage, it is important to accurately predict the simultaneous flow of gas, water, and hydrate particles in flowlines. Currently, there are few or no studies that account for the simultaneous effects of hydrate growth and heat transfer on flow characteristics within pipelines. This thesis presents new and more accurate predictive models of multiphase flows in undersea pipelines to describe the simultaneous flow of gas, water, and hydrate particles through a pipeline. A growth rate model for the hydrate phase is presented and then used in the development of a new three-phase model. The conservation equations of mass, momentum, and energy are formulated to describe the physical phenomena of momentum and heat transfer between the fluid and the wall. The governing equations are solved based on an analytical-numerical approach using a Newton-Raphson method for the nonlinear equations. An algorithm was developed in Matlab software to solve the equations from the inlet to the outlet of the pipeline. The developed models are validated against a single-phase model with mixture properties, and the results of comparative studies show close agreement. The new model predicts the volume fraction and velocity of each phase, as well as the mixture pressure and temperature profiles along the length of the pipeline. The results from the hydrate growth model reveal the growth rate and location where the initial hydrates start to form. Finally, to assess the impact of certain parameters on the flow characteristics, parametric studies have been conducted. The results show the effect of a variation in the pipe diameter, mass flow rate, inlet pressure, and inlet temperature on the flow characteristics and hydrate growth rates.

Coupled wellbore-reservoir simulation of thermal effects during multiphase CO2 injection

T2Well-ECO2M is a coupled wellbore reservoir simulator still under development at Lawrence Berkeley National Laboratory (USA) with the ability to deal with a mixture of H2O-CO2-NaCl and includes the simulation of CO2 phase transition and multiphase flow. The code was originally developed for the simulation of CO2 injection into deep saline aquifers and the modelling of enhanced geothermal systems; however, the focus of this research was to modify and test T2Well-ECO2M to simulate CO2 injection into depleted gas reservoirs. To this end, the original code was properly changed in a few parts and a dedicated injection case was developed to study CO2 phase transition inside of a wellbore and the corresponding thermal effects. In the first scenario, the injection case was run applying the fully numerical approach of wellbore to formation heat exchange calculation. Results were analysed in terms of wellbore pressure and temperature vertical profiles, wellhead and bottomhole conditions, and characteristic reservoir displacement fronts. Special attention was given to the thorough analysis of bottomhole temperature as the critical parameter for hydrate formation. Besides the expected direct effect of wellbore temperature changes on reservoir conditions, the simulation results indicated also the effect of CO2 phase change in the near wellbore zone on BH pressure distribution. To test the implemented software changes, in a second scenario, the same injection case was reproduced using the improved semi-analytical time-convolution approach for wellbore to formation heat exchange calculation. The comparison of the two scenarios showed that the simulation of wellbore and reservoir parameters after one year of continuous CO2 injection are in good agreement with the computation time to solve the time-convolution semi-analytical reduced. The new updated T2Well-ECO2M version has shown to be a robust and performing wellbore-reservoir simulator that can be also used to simulate the CO2 injection into depleted gas reservoirs.

Flow and axial dispersion simulation for traveling axisymmetric Taylor vortices

Numerical experiments using a finite difference method were carried out to determine the motion of axisymmetric Taylor vortices for narrow-gap Taylor vortex flow. When a pressure gradient is imposed on the flow the vortices are observed to move with an axial speed of 1.16 +/- 0.005 times the mean axial flow velocity. The method of Brenner was used to calculate the long-time axial spread of material in the flow. For flows where there is no pressure gradient, the axial dispersion scales with the square root of the molecular diffusion, in agreement with the results of Rosen-bluth et al. for high Peclet number dispersion in spatially periodic flows with a roll structure. When a pressure gradient is imposed the dispersion increases by an amount approximately equal to 6.5 x 10(-4) (W) over bar(2)d(2)/D-m, where (W) over bar is the average axial velocity in the annulus, analogous to Taylor dispersion for laminar flow in an empty tube.

Analytical and numerical study of the electro-osmotic annular flow of viscoelastic fluids

In this work we present semi-analytical solutions for the electro-osmotic annular flow of viscoelastic fluids modeled by the Linear and Exponential PTT models. The viscoelastic fluid flows in the axial direction between two concentric cylinders under the combined influences of electrokinetic and pressure forcings. The analysis invokes the Debye-Hückel approximation and includes the limit case of pure electro-osmotic flow. The solution is valid for both no slip and slip velocity at the walls and the chosen slip boundary condition is the linear Navier slip velocity model. The combined effects of fluid rheology, electro-osmotic and pressure gradient forcings on the fluid velocity distribution are also discussed.

Annular flow of viscoelastic fluids: Analytical and numerical solutions

This work provides analytical and numerical solutions for the linear, quadratic and exponential Phan–Thien–Tanner (PTT) viscoelastic models, for axial and helical annular fully-developed flows under no slip and slip boundary conditions, the latter given by the linear and nonlinear Navier slip laws. The rheology of the three PTT model functions is discussed together with the influence of the slip velocity upon the flow velocity and stress fields. For the linear PTT model, full analytical solutions for the inverse problem (unknown velocity) are devised for the linear Navier slip law and two different slip exponents. For the linear PTT model with other values of the slip exponent and for the quadratic PTT model, the polynomial equation for the radial location (β) of the null shear stress must be solved numerically. For both models, the solution of the direct problem is given by an iterative procedure involving three nonlinear equations, one for β, other for the pressure gradient and another for the torque per unit length. For the exponential PTT model we devise a numerical procedure that can easily compute the numerical solution of the pure axial flow problem

Remote control of pulmonary blood flow: a dream comes true.

The indication for pulmonary artery banding is currently limited by several factors. Previous attempts have failed to produce adjustable pulmonary artery banding with reliable external regulation. An implantable, telemetrically controlled, battery-free device (FloWatch) developed by EndoArt SA, a medical company established in Lausanne, Switzerland, for externally adjustable pulmonary artery banding was evaluated on minipigs and proved to be effective for up to 6 months. The first human implant was performed on a girl with complete atrioventricular septal defect with unbalanced ventricles, large patent ductus arteriosus and pulmonary hypertension. At one month of age she underwent closure of the patent ductus arteriosus and FloWatch implantation around the pulmonary artery through conventional left thoracotomy. The surgical procedure was rapid and uneventful. During the entire postoperative period bedside adjustments (narrowing or release of pulmonary artery banding with echocardiographic assessment) were repeatedly required to maintain an adequate pressure gradient. The early clinical results demonstrated the clinical benefits of unlimited external telemetric adjustments. The next step will be a multi-centre clinical trial to confirm the early results and adapt therapeutic strategies to this promising technology.

Velocity-jump instabilities in Hele-Shaw flow of associating polymer solutions

We study fracturelike flow instabilities that arise when water is injected into a Hele-Shaw cell filled with aqueous solutions of associating polymers. We explore various polymer architectures, molecular weights, and solution concentrations. Simultaneous measurements of the finger tip velocity and of the pressure at the injection point allow us to describe the dynamics of the finger in terms of the finger mobility, which relates the velocity to the pressure gradient. The flow discontinuities, characterized by jumps in the finger tip velocity, which are observed in experiments with some of the polymer solutions, can be modeled by using a nonmonotonic dependence between a characteristic shear stress and the shear rate at the tip of the finger. A simple model, which is based on a viscosity function containing both a Newtonian and a non-Newtonian component, and which predicts nonmonotonic regions when the non-Newtonian component of the viscosity dominates, is shown to agree with the experimental data.

Modeling of turbulent flow through radial diffuser

The work considers the modeling of turbulent flow in radial diffuser with axial feeding. Due to its claimed capability to predict flow including features such as separation, curvature and adverse pressure gradient, the RNG k-epsilon model of Orzag et al. (1993) is applied in the present analysis. The governing equations are numerically solved using the finite volume methodology. Experiments were conducted to assess the turbulence model. Numerical results of pressure distribution on the front disk surface for different flow conditions when compared to the experimental data indicated that the RNG k-epsilon model is adequate to predict this class of flow.

Flow in porous media of oil-water systems

Fluid flow behaviour in porous media is a conundrum. Therefore, this research is focused on filtration-volumetric characterisation of fractured-carbonate sediments, coupled with their proper simulation. For this reason, at laboratory rock properties such as pore volume, permeability and porosity are measured, later phase permeabilities and oil recovery in function of flow rate are assessed. Furthermore, the rheological properties of three oils are measured and analysed. Finally based on rock and fluid properties, a model using COMSOL Multiphysics is built in order to compare the experimental and simulated results. The rock analyses show linear relation between flow rate and differential pressure, from which phase permeabilities and pressure gradient are determined, eventually the oil recovery under low and high flow rate is established. In addition, the oils reveal thixotropic properties as well as non-Newtonian behaviour described by Bingham model, consequently Carreau viscosity model for the used oil is given. Given these points, the model for oil and water is built in COMSOL Multiphysics, whereupon successfully the reciprocity between experimental and simulated results is analysed and compared. Finally, a two-phase displacement model is elaborated.

Numerical Simulation of Electroosmotic Flow with Step Change in Zeta Potential

Electroosmotic flow is a convenient mechanism for transporting polar fluid in a microfluidic device. The flow is generated through the application of an external electric field that acts on the free charges that exists in a thin Debye layer at the channel walls. The charge on the wall is due to the chemistry of the solid-fluid interface, and it can vary along the channel, e.g. due to modification of the wall. This investigation focuses on the simulation of the electroosmotic flow (EOF) profile in a cylindrical microchannel with step change in zeta potential. The modified Navier-Stoke equation governing the velocity field and a non-linear two-dimensional Poisson-Boltzmann equation governing the electrical double-layer (EDL) field distribution are solved numerically using finite control-volume method. Continuities of flow rate and electric current are enforced resulting in a non-uniform electrical field and pressure gradient distribution along the channel. The resulting parabolic velocity distribution at the junction of the step change in zeta potential, which is more typical of a pressure-driven velocity flow profile, is obtained.