Unsteady flow simulation of a vertical axis augmented wind turbine: a two-dimensional study

As the integration of vertical axis wind turbines in the built environment is a promising alternative to horizontal axis wind turbines, a 2D computational investigation of an augmented wind turbine is proposed and analysed. In the initial CFD analysis, three parameters are carefully investigated: mesh resolution; turbulence model; and time step size. It appears that the mesh resolution and the turbulence model affect result accuracy; while the time step size examined, for the unsteady nature of the flow, has small impact on the numerical results. In the CFD validation of the open rotor with secondary data, the numerical results are in good agreement in terms of shape. It is, however, observed a discrepancy factor of 2 between numerical and experimental data. Successively, the introduction of an omnidirectional stator around the wind turbine increases the power and torque coefficients by around 30–35% when compared to the open case; but attention needs to be given to the orientation of the stator blades for optimum performance. It is found that the power and torque coefficients of the augmented wind turbine are independent of the incident wind speed considered.

Scattered power from non-thermal, F-region plasma observed by EISCAT—evidence for coherent echoes?

Three rapid, poleward bursts of plasma flow, observed by the U.K.-POLAR EISCAT experiment, are studied in detail. In all three cases the large ion velocities (> 1 kms−1) are shown to drive the ion velocity distribution into a non-Maxwellian form, identified by the characteristic shape of the observed spectra and the fact that analysis of the spectra with the assumption of a Maxwellian distribution leads to excessive rises in apparent ion temperature, and an anticorrelation of apparent electron and ion temperatures. For all three periods the total scattered power is shown to rise with apparent ion temperature by up to 6 dB more than is expected for an isotropic Maxwellian plasma of constant density and by an even larger factor than that expected for non-thermal plasma. The anomalous increases in power are only observed at the lower altitudes (< 300 km). At greater altitudes the rise in power is roughly consistent with that simulated numerically for homogeneous, anisotropic, non-Maxwellian plasma of constant density, viewed using the U.K.-POLAR aspect angle. The spectra at times of anomalously high power are found to be asymmetric, showing an enhancement near the downward Doppler-shifted ion-acoustic frequency. Although it is not possible to eliminate completely rapid plasma density fluctuations as a cause of these power increases, such effects cannot explain the observed spectra and the correlation of power and apparent ion temperature without an unlikely set of coincidences. The observations are made along a beam direction which is as much as 16.5° from orthogonality with the geomagnetic field. Nevertheless, some form of coherent-like echo contamination of the incoherent scatter spectrum is the most satisfactory explanation of these data.

Characterising phosphorus and nitrate inputs to a rural river using high-frequency concentration-flow relationships

The total reactive phosphorus (TRP) and nitrate concentrations of the River Enborne, southern England, were monitored at hourly interval between January 2010 and December 2011. The relationships between these high-frequency nutrient concentration signals and flow were used to infer changes in nutrient source and dynamics through the annual cycle and each individual storm event, by studying hysteresis patterns. TRP concentrations exhibited strong dilution patterns with increasing flow, and predominantly clockwise hysteresis through storm events. Despite the Enborne catchment being relatively rural for southern England, TRP inputs were dominated by constant, non-rain-related inputs from sewage treatment works (STW) for the majority of the year, producing the highest phosphorus concentrations through the spring–summer growing season. At higher river flows, the majority of the TRP load was derived from within-channel remobilisation of phosphorus from the bed sediment, much of which was also derived from STW inputs. Therefore, future phosphorus mitigation measures should focus on STW improvements. Agricultural diffuse TRP inputs were only evident during storms in the May of each year, probably relating to manure application to land. The nitrate concentration–flow relationship produced a series of dilution curves, indicating major inputs from groundwater and to a lesser extent STW. Significant diffuse agricultural inputs with anticlockwise hysteresis trajectories were observed during the first major storms of the winter period. The simultaneous investigation of high-frequency time series data, concentration–flow relationships and hysteresis behaviour through multiple storms for both phosphorus and nitrate offers a simple and innovative approach for providing new insights into nutrient sources and dynamics.

Obstacle avoidance routing scheme through optimal sink movement for home monitoring and mobile robotic consumer devices

In recent years, ZigBee has been proven to be an excellent solution to create scalable and flexible home automation networks. In a home automation network, consumer devices typically collect data from a home monitoring environment and then transmit the data to an end user through multi-hop communication without the need for any human intervention. However, due to the presence of typical obstacles in a home environment, error-free reception may not be possible, particularly for power constrained devices. A mobile sink based data transmission scheme can be one solution but obstacles create significant complexities for the sink movement path determination process. Therefore, an obstacle avoidance data routing scheme is of vital importance to the design of an efficient home automation system. This paper presents a mobile sink based obstacle avoidance routing scheme for a home monitoring system. The mobile sink collects data by traversing through the obstacle avoidance path. Through ZigBee based hardware implementation and verification, the proposed scheme successfully transmits data through the obstacle avoidance path to improve network performance in terms of life span, energy consumption and reliability. The application of this work can be applied to a wide range of intelligent pervasive consumer products and services including robotic vacuum cleaners and personal security robots1.

Mesh adaptation on the sphere using optimal transport and the numerical solution of a Monge-Ampère type equation

An equation of Monge-Ampère type has, for the first time, been solved numerically on the surface of the sphere in order to generate optimally transported (OT) meshes, equidistributed with respect to a monitor function. Optimal transport generates meshes that keep the same connectivity as the original mesh, making them suitable for r-adaptive simulations, in which the equations of motion can be solved in a moving frame of reference in order to avoid mapping the solution between old and new meshes and to avoid load balancing problems on parallel computers. The semi-implicit solution of the Monge-Ampère type equation involves a new linearisation of the Hessian term, and exponential maps are used to map from old to new meshes on the sphere. The determinant of the Hessian is evaluated as the change in volume between old and new mesh cells, rather than using numerical approximations to the gradients. OT meshes are generated to compare with centroidal Voronoi tesselations on the sphere and are found to have advantages and disadvantages; OT equidistribution is more accurate, the number of iterations to convergence is independent of the mesh size, face skewness is reduced and the connectivity does not change. However anisotropy is higher and the OT meshes are non-orthogonal. It is shown that optimal transport on the sphere leads to meshes that do not tangle. However, tangling can be introduced by numerical errors in calculating the gradient of the mesh potential. Methods for alleviating this problem are explored. Finally, OT meshes are generated using observed precipitation as a monitor function, in order to demonstrate the potential power of the technique.

FORCED CONVECTION TO LAMINAR FLOW OF LIQUID EGG YOLK IN CIRCULAR AND ANNULAR DUCTS

The steady-state heat transfer in laminar flow of liquid egg yolk - an important pseudoplastic fluid food - in circular and concentric annular ducts was experimentally investigated. The average convection heat transfer coefficients, determined by measuring temperatures before and after heating sections with constant temperatures at the tube wall, were used to obtain simple new empirical expressions to estimate the Nusselt numbers for fully established flows at the thermal entrance of the considered geometries. The comparisons with existing correlations for Newtonian and non-Newtonian fluids resulted in excellent agreement. The main contribution of this work is to supply practical and easily applicable correlations, which are, especially for the case of annulus, rather scarce and extensively required in the design of heat transfer operations dealing with similar shear-thinning products. In addition, the experimental results may support existing theoretical analyses.

Fully electrochemical hyphenated flow system for preconcentration, cleanup, stripping, capillary electrophoresis with stacking and contactless conductivity detection of trace heavy metals

Successful coupling of electrochemical preconcentration (EPC) to capillary electrophoresis (CE) with contactless conductivity detection (C(4)D) is reported for the first time. The EPC-CE interface comprises a dual glassy carbon electrode (GCE) block, a spacer and an upper block with flow inlet and outlet, pseudo-reference electrode and a fitting for the CE silica column, consisting of an orifice perpendicular to the surface of a glassy carbon electrode with a bushing inside to ensure a tight press fit. The end of the capillary in contact with the GCE is slant polished, thus defining a reproducible distance from the electrode surface to the column bore. First results with EPC-CE-C(4)D are very promising, as revealed by enrichment factors of two orders of magnitude for Tl, Cu, Pb and Cd ion peak area signals. Detection limits for 10 min deposition time fall around 20 nmol L(-1) with linear calibration curves over a wide range. Besides preconcentration, easy matrix exchange between accumulation and stripping/injection favors procedures like sample cleanup and optimization of pH, ionic strength and complexing power. This was demonstrated for highly saline samples by using a low conductivity buffer for stripping/injection to improve separation and promote field-enhanced sample stacking during electromigration along the capillary. (C) 2010 Elsevier B.V. All rights reserved.

Determination of catechins in green tea infusions by reduced flow micellar electrokinetic chromatography

A sulfated-beta-cyclodextrin (s-beta-CD) modified reduced flow micellar electrokinetic chromatography (RF-MEKC) method was developed and validated for the determination of catechins in green tea. The optimal electrolyte consisted of 0.2% triethylamine, 50 mmol/L SDS and 0.8% s-beta-CD (pH = 2.9), allowing baseline separation of five catechins in 4 min. The samples and standards were injected at 0.6 psi for 5 s under constant voltage of -30 kV. Sample preparation simply involved extraction of 2 g of tea with 200 mL water at 95 C under constant stirring for 5 min. The method demonstrated excellent performance, with limits of detection (LOD) and quantification (LOQ) of 0.02-0.1 and 0.1-0.5 mu g/mL, respectively, and recovery percentages of 94-101%. The method was applied to six samples of Brazilian green tea infusions. Epigallocatechin gallate (23.4-112.4 mu g/mL) was the major component, followed by epigallocatechin (18.4-78.9 mu g/mL), epicatechin gallate (5.6-29.6 mu g/mL), epicatechin (4.6-14.5 mu g/mL) and catechin (3.2-8.2 mu g/mL). (C) 2011 Elsevier Ltd. All rights reserved.

Pt film electrodes prepared by the Pechini method for electrochemical decolourisation of Reactive Orange 16

Electrochemical decolourisation of Reactive Orange 16 was carried out in an electrochemical flow-cell, using as working electrodes a Pt thin film deposited on a Ti substrate (Pt/Ti) prepared by the Pechini method and a pure platinum (Pt) foil. Using the Pt/Ti electrodes better results for dye decolourisation were obtained under milder conditions than those used for pure Pt. For the Pt electrode, colour removal of 93 % (lambda = 493 nm) was obtained after 60 min, at 2.2 V vs. RHE, using 0.017 mol L(-1) NaCl + 0.5 mol L(-1) H(2)SO(4) solution. For the Pt/Ti electrode there was better colour removal, 98%, than for the Pt electrode. Moreover, we used 0.017 mol L(-1) NaCl solution and the applied potential was 1.8 V. Under this condition after 15 min of electrolysis, more than 80% of colour was removed. The rate reaction constant, assuming a first order reaction, was 0.024 min(-1) and 0.069 min(-1), for Pt and Pt/Ti electrodes, respectively.

Optimal V. O2max-to-mass ratio for predicting 15 km performance among elite male cross-country skiers

The aim of this study was 1) to validate the 0.5 body-mass exponent for maximal oxygen uptake (V. O2max) as the optimal predictor of performance in a 15 km classical-technique skiing competition among elite male cross-country skiers and 2) to evaluate the influence of distance covered on the body-mass exponent for V. O2max among elite male skiers. Twenty-four elite male skiers (age: 21.4±3.3 years [mean ± standard deviation]) completed an incremental treadmill roller-skiing test to determine their V. O2max. Performance data were collected from a 15 km classicaltechnique cross-country skiing competition performed on a 5 km course. Power-function modeling (ie, an allometric scaling approach) was used to establish the optimal body-mass exponent for V . O2max to predict the skiing performance. The optimal power-function models were found to be race speed = 8.83⋅(V . O2max m-0.53) 0.66 and lap speed = 5.89⋅(V . O2max m-(0.49+0.018lap)) 0.43e0.010age, which explained 69% and 81% of the variance in skiing speed, respectively. All the variables contributed to the models. Based on the validation results, it may be recommended that V. O2max divided by the square root of body mass (mL⋅min−1 ⋅kg−0.5) should be used when elite male skiers’ performance capability in 15 km classical-technique races is evaluated. Moreover, the body-mass exponent for V . O2max was demonstrated to be influenced by the distance covered, indicating that heavier skiers have a more pronounced positive pacing profile (ie, race speed gradually decreasing throughout the race) compared to that of lighter skiers.

Contributions to the Theory of Optimal Tests

This paper considers tests which maximize the weighted average power (WAP). The focus is on determining WAP tests subject to an uncountable number of equalities and/or inequalities. The unifying theory allows us to obtain tests with correct size, similar tests, and unbiased tests, among others. A WAP test may be randomized and its characterization is not always possible. We show how to approximate the power of the optimal test by sequences of nonrandomized tests. Two alternative approximations are considered. The rst approach considers a sequence of similar tests for an increasing number of boundary conditions. This discretization allows us to implement the WAP tests in practice. The second method nds a sequence of tests which approximate the WAP test uniformly. This approximation allows us to show that WAP similar tests are admissible. The theoretical framework is readily applicable to several econometric models, including the important class of the curved-exponential family. In this paper, we consider the instrumental variable model with heteroskedastic and autocorrelated errors (HAC-IV) and the nearly integrated regressor model. In both models, we nd WAP similar and (locally) unbiased tests which dominate other available tests.

Optimal two-sided tests for instrumental variables regression with heteroskedastic and autocorrelated errors

This paper considers two-sided tests for the parameter of an endogenous variable in an instrumental variable (IV) model with heteroskedastic and autocorrelated errors. We develop the nite-sample theory of weighted-average power (WAP) tests with normal errors and a known long-run variance. We introduce two weights which are invariant to orthogonal transformations of the instruments; e.g., changing the order in which the instruments appear. While tests using the MM1 weight can be severely biased, optimal tests based on the MM2 weight are naturally two-sided when errors are homoskedastic. We propose two boundary conditions that yield two-sided tests whether errors are homoskedastic or not. The locally unbiased (LU) condition is related to the power around the null hypothesis and is a weaker requirement than unbiasedness. The strongly unbiased (SU) condition is more restrictive than LU, but the associated WAP tests are easier to implement. Several tests are SU in nite samples or asymptotically, including tests robust to weak IV (such as the Anderson-Rubin, score, conditional quasi-likelihood ratio, and I. Andrews' (2015) PI-CLC tests) and two-sided tests which are optimal when the sample size is large and instruments are strong. We refer to the WAP-SU tests based on our weights as MM1-SU and MM2-SU tests. Dropping the restrictive assumptions of normality and known variance, the theory is shown to remain valid at the cost of asymptotic approximations. The MM2-SU test is optimal under the strong IV asymptotics, and outperforms other existing tests under the weak IV asymptotics.

Uso de técnicas de otimização baseadas em derivadas como suporte do planejamento operacional de redes de distribuição de energia elétrica

The usual programs for load flow calculation were in general developped aiming the simulation of electric energy transmission, subtransmission and distribution systems. However, the mathematical methods and algorithms used by the formulations were based, in majority, just on the characteristics of the transmittion systems, which were the main concern focus of engineers and researchers. Though, the physical characteristics of these systems are quite different from the distribution ones. In the transmission systems, the voltage levels are high and the lines are generally very long. These aspects contribute the capacitive and inductive effects that appear in the system to have a considerable influence in the values of the interest quantities, reason why they should be taken into consideration. Still in the transmission systems, the loads have a macro nature, as for example, cities, neiborhoods, or big industries. These loads are, generally, practically balanced, what reduces the necessity of utilization of three-phase methodology for the load flow calculation. Distribution systems, on the other hand, present different characteristics: the voltage levels are small in comparison to the transmission ones. This almost annul the capacitive effects of the lines. The loads are, in this case, transformers, in whose secondaries are connected small consumers, in a sort of times, mono-phase ones, so that the probability of finding an unbalanced circuit is high. This way, the utilization of three-phase methodologies assumes an important dimension. Besides, equipments like voltage regulators, that use simultaneously the concepts of phase and line voltage in their functioning, need a three-phase methodology, in order to allow the simulation of their real behavior. For the exposed reasons, initially was developped, in the scope of this work, a method for three-phase load flow calculation in order to simulate the steady-state behaviour of distribution systems. Aiming to achieve this goal, the Power Summation Algorithm was used, as a base for developing the three phase method. This algorithm was already widely tested and approved by researchers and engineers in the simulation of radial electric energy distribution systems, mainly for single-phase representation. By our formulation, lines are modeled in three-phase circuits, considering the magnetic coupling between the phases; but the earth effect is considered through the Carson reduction. It s important to point out that, in spite of the loads being normally connected to the transformer s secondaries, was considered the hypothesis of existence of star or delta loads connected to the primary circuit. To perform the simulation of voltage regulators, a new model was utilized, allowing the simulation of various types of configurations, according to their real functioning. Finally, was considered the possibility of representation of switches with current measuring in various points of the feeder. The loads are adjusted during the iteractive process, in order to match the current in each switch, converging to the measured value specified by the input data. In a second stage of the work, sensibility parameters were derived taking as base the described load flow, with the objective of suporting further optimization processes. This parameters are found by calculating of the partial derivatives of a variable in respect to another, in general, voltages, losses and reactive powers. After describing the calculation of the sensibility parameters, the Gradient Method was presented, using these parameters to optimize an objective function, that will be defined for each type of study. The first one refers to the reduction of technical losses in a medium voltage feeder, through the installation of capacitor banks; the second one refers to the problem of correction of voltage profile, through the instalation of capacitor banks or voltage regulators. In case of the losses reduction will be considered, as objective function, the sum of the losses in all the parts of the system. To the correction of the voltage profile, the objective function will be the sum of the square voltage deviations in each node, in respect to the rated voltage. In the end of the work, results of application of the described methods in some feeders are presented, aiming to give insight about their performance and acuity

Técnicas de simplificação de redes e otimização multiobjetivo para análise de variações de tensão em regime permanente provocadas por parques eólicos integrados ao sistema elétrico

This study presents a description of the development model of a representation of simplified grid applied in hybrid load flow for calculation of the voltage variations in a steady-state caused by the wind farm on power system. Also, it proposes an optimal load-flow able to control power factor on connection bar and to minimize the loss. The analysis process on system, led by the wind producer, it has as base given technician supplied by the grid. So, the propose model to the simplification of the grid that allows the necessity of some knowledge only about the data referring the internal network, that is, the part of the network that interests in the analysis. In this way, it is intended to supply forms for the auxiliary in the systematization of the relations between the sector agents. The model for simplified network proposed identifies the internal network, external network and the buses of boulders from a study of vulnerability of the network, attributing them floating liquid powers attributing slack models. It was opted to apply the presented model in Newton-Raphson and a hybrid load flow, composed by The Gauss-Seidel method Zbarra and Summation Power. Finally, presents the results obtained to a developed computational environment of SCILAB and FORTRAN, with their respective analysis and conclusion, comparing them with the ANAREDE

Uso de técnicas de otimização baseadas em derivadas como suporte do planejamento operacional de redes de distribuição de energia elétrica

The usual programs for load flow calculation were in general developped aiming the simulation of electric energy transmission, subtransmission and distribution systems. However, the mathematical methods and algorithms used by the formulations were based, in majority, just on the characteristics of the transmittion systems, which were the main concern focus of engineers and researchers. Though, the physical characteristics of these systems are quite different from the distribution ones. In the transmission systems, the voltage levels are high and the lines are generally very long. These aspects contribute the capacitive and inductive effects that appear in the system to have a considerable influence in the values of the interest quantities, reason why they should be taken into consideration. Still in the transmission systems, the loads have a macro nature, as for example, cities, neiborhoods, or big industries. These loads are, generally, practically balanced, what reduces the necessity of utilization of three-phase methodology for the load flow calculation. Distribution systems, on the other hand, present different characteristics: the voltage levels are small in comparison to the transmission ones. This almost annul the capacitive effects of the lines. The loads are, in this case, transformers, in whose secondaries are connected small consumers, in a sort of times, mono-phase ones, so that the probability of finding an unbalanced circuit is high. This way, the utilization of three-phase methodologies assumes an important dimension. Besides, equipments like voltage regulators, that use simultaneously the concepts of phase and line voltage in their functioning, need a three-phase methodology, in order to allow the simulation of their real behavior. For the exposed reasons, initially was developped, in the scope of this work, a method for three-phase load flow calculation in order to simulate the steady-state behaviour of distribution systems. Aiming to achieve this goal, the Power Summation Algorithm was used, as a base for developping the three phase method. This algorithm was already widely tested and approved by researchers and engineers in the simulation of radial electric energy distribution systems, mainly for single-phase representation. By our formulation, lines are modeled in three-phase circuits, considering the magnetic coupling between the phases; but the earth effect is considered through the Carson reduction. Its important to point out that, in spite of the loads being normally connected to the transformers secondaries, was considered the hypothesis of existence of star or delta loads connected to the primary circuit. To perform the simulation of voltage regulators, a new model was utilized, allowing the simulation of various types of configurations, according to their real functioning. Finally, was considered the possibility of representation of switches with current measuring in various points of the feeder. The loads are adjusted during the iteractive process, in order to match the current in each switch, converging to the measured value specified by the input data. In a second stage of the work, sensibility parameters were derived taking as base the described load flow, with the objective of suporting further optimization processes. This parameters are found by calculating of the partial derivatives of a variable in respect to another, in general, voltages, losses and reactive powers. After describing the calculation of the sensibility parameters, the Gradient Method was presented, using these parameters to optimize an objective function, that will be defined for each type of study. The first one refers to the reduction of technical losses in a medium voltage feeder, through the installation of capacitor banks; the second one refers to the problem of correction of voltage profile, through the instalation of capacitor banks or voltage regulators. In case of the losses reduction will be considered, as objective function, the sum of the losses in all the parts of the system. To the correction of the voltage profile, the objective function will be the sum of the square voltage deviations in each node, in respect to the rated voltage. In the end of the work, results of application of the described methods in some feeders are presented, aiming to give insight about their performance and acuity