Security constrained optimal active power flow via network model and interior point method

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Fourth-order method for solving the Navier-Stokes equations in a constricting channel

A fourth-order numerical method for solving the Navier-Stokes equations in streamfunction/vorticity formulation on a two-dimensional non-uniform orthogonal grid has been tested on the fluid flow in a constricted symmetric channel. The family of grids is generated algebraically using a conformal transformation followed by a non-uniform stretching of the mesh cells in which the shape of the channel boundary can vary from a smooth constriction to one which one possesses a very sharp but smooth corner. The generality of the grids allows the use of long channels upstream and downstream as well as having a refined grid near the sharp corner. Derivatives in the governing equations are replaced by fourth-order central differences and the vorticity is eliminated, either before or after the discretization, to form a wide difference molecule for the streamfunction. Extra boundary conditions, necessary for wide-molecule methods, are supplied by a procedure proposed by Henshaw et al. The ensuing set of non-linear equations is solved using Newton iteration. Results have been obtained for Reynolds numbers up to 250 for three constrictions, the first being smooth, the second having a moderately sharp corner and the third with a very sharp corner. Estimates of the error incurred show that the results are very accurate and substantially better than those of the corresponding second-order method. The observed order of the method has been shown to be close to four, demonstrating that the method is genuinely fourth-order. © 1977 John Wiley & Sons, Ltd.

Robust control applied to power flow control in single-phase inverter with LCL filter, using droop control and D-stability

This paper proposes a new methodology to control the power flow between a distributed generator (DG) and the electrical power distribution grid. It is used the droop voltage control to manage the active and reactive power. Through this control a sinusoidal voltage reference is generated to be tracked by voltage loop and this loop generates the current reference for the current loop. The proposed control introduces feed-forward states improving the control performance in order to obtain high quality for the current injected to the grid. The controllers were obtained through the linear matrix inequalities (LMI) using the D-stability analysis to allocate the closed-loop controller poles. Therefore, the results show quick transient response with low oscillations. Thus, this paper presents the proposed control technique, the main simulation results and a prototype with 1000VA was developed in the laboratory in order to demonstrate the feasibility of the proposed control. © 2012 IEEE.

Estudo teórico-experimental em via de determinação de lei de atrito em escoamentos de fluidos hiperconcentrados

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A numerical study of temperature effects on hot wire measurements inside turbulent channel flows

A way to investigate turbulence is through experiments where hot wire measurements are performed. Analysis of the in turbulence of a temperature gradient on hot wire measurements is the aim of this thesis work. Actually - to author's knowledge - this investigation is the first attempt to document, understand and ultimately correct the effect of temperature gradients on turbulence statistics. However a numerical approach is used since instantaneous temperature and streamwise velocity fields are required to evaluate this effect. A channel flow simulation at Re_tau = 180 is analyzed to make a first evaluation of the amount of error introduced by temperature gradient inside the domain. Hot wire data field is obtained processing the numerical flow field through the application of a proper version of the King's law, which connect voltage, velocity and temperature. A drift in mean streamwise velocity profile and rms is observed when temperature correction is performed by means of centerline temperature. A correct mean velocity pro�le is achieved correcting temperature through its mean value at each wall normal position, but a not negligible error is still present into rms. The key point to correct properly the sensed velocity from the hot wire is the knowledge of the instantaneous temperature field. For this purpose three correction methods are proposed. At the end a numerical simulation at Re_tau =590 is also evaluated in order to confirm the results discussed earlier.

Does the beta-blocker nebivolol increase coronary flow reserve?

INTRODUCTION: Nebivolol, a highly selective beta1-adrenergic receptor-blocker, increases basal and stimulated endothelial nitric oxide (NO)-release. It is unknown, whether coronary perfusion is improved by the increase in NO availability. Therefore, we sought to evaluate the effect of nebivolol on coronary flow reserve (CFR) and collateral flow. METHODS: Doppler-flow wire derived coronary flow velocity measurements were obtained in ten controls and eight patients with coronary artery disease (CAD) at rest and after intracoronary nebivolol. CFR was defined as maximal flow during adenosine-induced hyperemia divided by resting flow. In the CAD group, collateral flow was determined after dilatation of a flow-limiting coronary stenosis. Collateral flow index (CFI) was defined as the ratio of flow velocity during balloon inflation divided by resting flow. RESULTS: CFR at rest was 3.0+/-0.6 in controls and 2.1+/-0.4 in CAD patients. After intracoronary doses of 0.1, 0.25, and 0.5 mg nebivolol, CFR increased to 3.4+/-0.7, 3.9+/-0.9, and 4.0+/-0.1 (p<0.01) in controls, and to 2.3+/-0.7, 2.6+/-0.9, and 2.6+/-0.5 (p<0.05) in CAD patients. CFI decreased significantly with intracoronary nebivolol and correlated to changes in heart rate (r=0.75, p<0.001) and rate-pressure product (r=0.59, p=0.001). DISCUSSION: Intracoronary nebivolol is associated with a significant increase in CFR due to reduction in resting flow (controls), or due to an increase in maximal coronary flow (CAD patients). CFI decreased with nebivolol parallel to the reduction in myocardial oxygen consumption.

Reynolds-stress turbulence model in the KIVA code for engine simulation

A Reynolds-Stress Turbulence Model has been incorporated with success into the KIVA code, a computational fluid dynamics hydrocode for three-dimensional simulation of fluid flow in engines. The newly implemented Reynolds-stress turbulence model greatly improves the robustness of KIVA, which in its original version has only eddy-viscosity turbulence models. Validation of the Reynolds-stress turbulence model is accomplished by conducting pipe-flow and channel-flow simulations, and comparing the computed results with experimental and direct numerical simulation data. Flows in engines of various geometry and operating conditions are calculated using the model, to study the complex flow fields as well as confirm the model’s validity. Results show that the Reynolds-stress turbulence model is able to resolve flow details such as swirl and recirculation bubbles. The model is proven to be an appropriate choice for engine simulations, with consistency and robustness, while requiring relatively low computational effort.

Membrane-bound state of the colicin E1 channel domain as an extended two-dimensional helical array

Atomic level structures have been determined for the soluble forms of several colicins and toxins, but the structural changes that occur after membrane binding have not been well characterized. Changes occurring in the transition from the soluble to membrane-bound state of the C-terminal 190-residue channel polypeptide of colicin E1 (P190) bound to anionic membranes are described. In the membrane-bound state, the α-helical content increases from 60–64% to 80–90%, with a concomitant increase in the average length of the helical segments from 12 to 16 or 17 residues, close to the length required to span the membrane bilayer in the open channel state. The average distance between helical segments is increased and interhelix interactions are weakened, as shown by a major loss of tertiary structure interactions, decreased efficiency of fluorescence resonance energy transfer from an energy donor on helix V of P190 to an acceptor on helix IX, and decreased resonance energy transfer at higher temperatures, not observed in soluble P190, implying freedom of motion of helical segments. Weaker interactions are also shown by a calorimetric thermal transition of low cooperativity, and the extended nature of the helical array is shown by a 3- to 4-fold increase in the average area subtended per molecule to 4,200 Å2 on the membrane surface. The latter, with analysis of the heat capacity changes, implies the absence of a developed hydrophobic core in the membrane-bound P190. The membrane interfacial layer thus serves to promote formation of a highly helical extended two-dimensional flexible net. The properties of the membrane-bound state of the colicin channel domain (i.e., hydrophobic anchor, lengthened and loosely coupled α-helices, and close association with the membrane interfacial layer) are plausible structural features for the state that is a prerequisite for voltage gating, formation of transmembrane helices, and channel opening.

An N-methyl-d-aspartate receptor channel blocker with neuroprotective activity

Excitotoxicity, resulting from sustained activation of glutamate receptors of the N-methyl-d-aspartate (NMDA) subtype, is considered to play a causative role in the etiology of ischemic stroke and several neurodegenerative diseases. The NMDA receptor is therefore a target for the development of neuroprotective agents. Here, we identify an N-benzylated triamine (denoted as NBTA) as a highly selective and potent NMDA-receptor channel blocker selected by screening a reduced dipeptidomimetic synthetic combinatorial library. NBTA blocks recombinant NMDA receptors expressed in Xenopus laevis oocytes with a mean IC50 of 80 nM; in contrast, it does not block GluR1, a glutamate receptor of the non-NMDA subtype. The blocking activity of NBTA on NMDA receptors exhibits the characteristics of an open-channel blocker: (i) no competition with agonists, (ii) voltage dependence, and (iii) use dependence. Significantly, NBTA protects rodent hippocampal neurons from NMDA receptor, but not kainate receptor-mediated excitotoxic cell death, in agreement with its selective action on the corresponding recombinant receptors. Mutagenesis data indicate that the N site, a key asparagine on the M2 transmembrane segment of the NR1 subunit, is the main determinant of the blocker action. The results highlight the potential of this compound as a neuroprotectant.

A conditional probability analysis of cystic fibrosis transmembrane conductance regulator gating indicates that ATP has multiple effects during the gating cycle

ATP-binding cassette (ABC) transporters bind and hydrolyze ATP. In the cystic fibrosis transmembrane conductance regulator Cl− channel, this interaction with ATP generates a gating cycle between a closed (C) and two open (O1 and O2) conformations. To understand better how ATP controls channel activity, we examined gating transitions from the C to the O1 and O2 states and from these open states to the C conformation. We made three main observations. First, we found that the channel can open into either the O1 or O2 state, that the frequency of transitions to both states was increased by ATP concentration, and that ATP increased the relative proportion of openings into O1 vs. O2. These results indicate that ATP can interact with the closed state to open the channel in at least two ways, which may involve binding to nucleotide-binding domains (NBDs) NBD1 and NBD2. Second, ATP prolonged the burst duration and altered the way in which the channel closed. These data suggest that ATP also interacts with the open channel. Third, the channel showed runs of specific types of open–closed transitions. This finding suggests a mechanism with more than one cycle of gating transitions. These data suggest models to explain how ATP influences conformational transitions in cystic fibrosis transmembrane conductance regulator and perhaps other ABC transporters.

Mutations within the selectivity filter of the NMDA receptor-channel influence voltage dependent block by 5-hydroxytryptamine

Background and purpose: Voltage-dependent block by Mg2+ is a cardinal feature of NMDA receptors which acts as a coincidence detector to prevent the receptor from over-activation. Inhibition of NMDA receptor currents by 5-hydroxytryptamine (5-HT) indicated that 5-HT, similar to Mg2+, binds within the membrane electric field. In the present study, we assessed whether point mutations of critical asparagine residues located within the selectivity filter of NR1 and NR2A subunits of NMDA receptor-channel affect voltage-dependent block by 5-HT. Experimental approach: The mode of action of 5-HT and Mg2+ on wild-type and mutated NMDA receptor-channels expressed in Xenopus oocytes was investigated using the two-electrode voltage clamp recording technique. Key results: The mutation within the NR1 subunit NR1(N0S or N0Q) strongly reduced the voltage dependent block by 5-HT and increased the IC50. The corresponding mutations within the NR2 subunits NR2A(N0Q or N + 1Q) reduced the block by 5-HT to a lesser extent. This is in contrast to the block produced by external Mg2+ where a substitution at the NR2A(N0) and NR2A(N + 1) sites but not at the NR1(N0) site significantly reduced Mg2+ block. Conclusion and implications: The block of NMDA receptor-channels by 5-HT depends on the NR1-subunit asparagine residue and to a lesser extent on the NR2A-subunit asparagine residues. These data suggest that the interaction of 5-HT with functionally important residues in a narrow constriction of the pore of the NMDA receptor-channel provides a significant barrier to ionic fluxes through the open channel due to energetic factors governed by chemical properties of the binding site and the electric field.

Direct numerical simulation of turbulent flow and mass transfer around a rotating circular cylinder

Turbulent flow around a rotating circular cylinder has numerous applications including wall shear stress and mass-transfer measurement related to the corrosion studies. It is also of interest in the context of flow over convex surfaces where standard turbulence models perform poorly. The main purpose of this paper is to elucidate the basic turbulence mechanism around a rotating cylinder at low Reynolds numbers to provide a better understanding of flow fundamentals. Direct numerical simulation (DNS) has been performed in a reference frame rotating at constant angular velocity with the cylinder. The governing equations are discretized by using a finite-volume method. As for fully developed channel, pipe, and boundary layer flows, a laminar sublayer, buffer layer, and logarithmic outer region were observed. The level of mean velocity is lower in the buffer and outer regions but the logarithmic region still has a slope equal to the inverse of the von Karman constant. Instantaneous flow visualization revealed that the turbulence length scale typically decreases as the Reynolds number increases. Wavelet analysis provided some insight into the dependence of structural characteristics on wave number. The budget of the turbulent kinetic energy was computed and found to be similar to that in plane channel flow as well as in pipe and zero pressure gradient boundary layer flows. Coriolis effects show as an equivalent production for the azimuthal and radial velocity fluctuations leading to their ratio being lowered relative to similar nonrotating boundary layer flows.

Numerical modelling of unsteady mixed flows

This thesis concerns mixed flows (which are characterized by the simultaneous occurrence of free-surface and pressurized flow in sewers, tunnels, culverts or under bridges), and contributes to the improvement of the existing numerical tools for modelling these phenomena. The classic Preissmann slot approach is selected due to its simplicity and capability of predicting results comparable to those of a more recent and complex two-equation model, as shown here with reference to a laboratory test case. In order to enhance the computational efficiency, a local time stepping strategy is implemented in a shock-capturing Godunov-type finite volume numerical scheme for the integration of the de Saint-Venant equations. The results of different numerical tests show that local time stepping reduces run time significantly (between −29% and −85% CPU time for the test cases considered) compared to the conventional global time stepping, especially when only a small region of the flow field is surcharged, while solution accuracy and mass conservation are not impaired. The second part of this thesis is devoted to the modelling of the hydraulic effects of potentially pressurized structures, such as bridges and culverts, inserted in open channel domains. To this aim, a two-dimensional mixed flow model is developed first. The classic conservative formulation of the 2D shallow water equations for free-surface flow is adapted by assuming that two fictitious vertical slots, normally intersecting, are added on the ceiling of each integration element. Numerical results show that this schematization is suitable for the prediction of 2D flooding phenomena in which the pressurization of crossing structures can be expected. Given that the Preissmann model does not allow for the possibility of bridge overtopping, a one-dimensional model is also presented in this thesis to handle this particular condition. The flows below and above the deck are considered as parallel, and linked to the upstream and downstream reaches of the channel by introducing suitable internal boundary conditions. The comparison with experimental data and with the results of HEC-RAS simulations shows that the proposed model can be a useful and effective tool for predicting overtopping and backwater effects induced by the presence of bridges and culverts.

Estimation of Mass Transfer Velocity Based on Measured Turbulence Parameters

The aim of this study is to quantify the mass transfer velocity using turbulence parameters from simultaneous measurements of oxygen concentration fields and velocity fields. The surface divergence model was considered in more detail, using data obtained for the lower range of beta (surface divergence). It is shown that the existing models that use the divergence concept furnish good predictions for the transfer velocity also for low values of beta, in the range of this study. Additionally, traditional conceptual models, such as the film model, the penetration-renewal model, and the large eddy model, were tested using the simultaneous information of concentration and velocity fields. It is shown that the film and the surface divergence models predicted the mass transfer velocity for all the range of the equipment Reynolds number used here. The velocity measurements showed viscosity effects close to the surface, which indicates that the surface was contaminated with some surfactant. Considering the results, this contamination can be considered slight for the mass transfer predictions. (C) 2009 American Institute of Chemical Engineers AIChE J, 56: 2005-2017; 2010

Large-conductance calcium-activated potassium channels in neonatal rat intracardiac ganglion neurons

The properties of single Ca2+-activated K+ (BK) channels in neonatal rat intracardiac neurons were investigated using the patch-clamp recording technique. In symmetrical 140 mM K+, the single-channel slope conductance was linear in the voltage range -60/+60 mV. and was 207+/-19 pS. Na+ ions were not measurably permeant through the open channel. Channel activity increased with the cytoplasmic free Ca2+ concentration ([Ca2+],) with a Hill plot giving a half-saturating [Ca2+] (K-0.5) of 1.35 muM and slope of congruent to3. The BK channel was inhibited reversibly by external tetraethylammonium (TEA) ions, charybdotoxin, and quinine and was resistant to block by 4-aminopyridine and apamin. Ionomycin (1-10 muM) increased BK channel activity in the cell-attached recording configuration. The resting activity was consistent with a [Ca2+](i)