Levenberg-Marquardt application to two-phase nonlinear parameter estimation for finned-tube coil evaporators

A procedure for calculation of refrigerant mass flow rate is implemented in the distributed numerical model to simulate the flow in finned-tube coil dry-expansion evaporators, usually found in refrigeration and air-conditioning systems. Two-phase refrigerant flow inside the tubes is assumed to be one-dimensional, unsteady, and homogeneous. In themodel the effects of refrigerant pressure drop and the moisture condensation from the air flowing over the external surface of the tubes are considered. The results obtained are the distributions of refrigerant velocity, temperature and void fraction, tube-wall temperature, air temperature, and absolute humidity. The finite volume method is used to discretize the governing equations. Additionally, given the operation conditions and the geometric parameters, the model allows the calculation of the refrigerant mass flow rate. The value of mass flow rate is computed using the process of parameter estimation with the minimization method of Levenberg-Marquardt minimization. In order to validate the developed model, the obtained results using HFC-134a as a refrigerant are compared with available data from the literature.

Measurement of energy flow at large pseudorapidities in pp collisions at root s=0.9 and 7 TeV

The energy flow, dE/d eta, is studied at large pseudorapidities in proton-proton collisions at the LHC, for centre-of-mass energies of 0.9 and 7 TeV. The measurements are made using the CMS detector in the pseudorapidity range 3:15 < vertical bar eta vertical bar < 4.9, for both minimum-bias events and events with at least two high-momentum jets. The data are compared to various pp Monte Carlo event generators whose theoretical models and input parameter values are sensitive to the energy-flow measurements. Inclusion of multiple-parton interactions in the Monte Carlo event generators is found to improve the description of the energy-flow measurements.

Continuation load flow using automatically determined branch megawatt losses as parameters

The conventional Newton's method is considered to be inadequate for the computation of the maximum loading point (MLP) of power systems since: (i) it encounters difficulties in the vicinity of the MLP: and (ii) the load flow Jacobian matrix becomes singular at the MLP. It is well known that continuation methods are powerful and useful tools that are able to trace the solution PV curve without experiencing such diffculties. However, continuation methods require a parameterisation so that a modified, well conditioned set of load flow equations is obtained. In particular, the Jacobian matrix associated with this modified set of equations should not be singular at the MLP. The authors propose that the actual power losses in transmission branches (lines and transformers) are used to parameterise the approach. Specific procedures for the automatic determination of the most appropriate parameter (branch) are proposed. Such procedures include the utilisation of fast voltage-stability indices. Simulation results are presented to show that the proposed method is able to trace the whole solution PV curve very efficiently.

Nonlinear parameter estimation in laminar forced convection within a circular sector tube

In this article we examine an inverse heat convection problem of estimating unknown parameters of a parameterized variable boundary heat flux. The physical problem is a hydrodynamically developed, thermally developing, three-dimensional steady state laminar flow of a Newtonian fluid inside a circular sector duct, insulated in the flat walls and subject to unknown wall heat flux at the curved wall. Results are presented for polynomial and sinusoidal trial functions, and the unknown parameters as well as surface heat fluxes are determined. Depending on the nature of the flow, on the position of experimental points the inverse problem sometimes could not be solved. Therefore, an identification condition is defined to specify a condition under which the inverse problem can be solved. Once the parameters have been computed it is possible to obtain the statistical significance of the inverse problem solution. Therefore, approximate confidence bounds based on standard statistical linear procedure, for the estimated parameters, are analyzed and presented.

Fuzzy sets in distributed traffic control

The task of controlling urban traffic requires flexibility, adaptability and handling uncertain information spread through the intersection network. The use of fuzzy sets concepts convey these characteristics to improve system performance. This paper reviews a distributed traffic control system built upon a fuzzy distributed architecture previously developed by the authors. The emphasis of the paper is on the application of the system to control part of Campinas downtown area. Simulation experiments considering several traffic scenarios were performed to verify the capabilities of the system in controlling a set of coupled intersections. The performance of the proposed system is compared with conventional traffic control strategies under the same scenarios. The results obtained show that the distributed traffic control system outperforms conventional systems as far as average queues, average delay and maximum delay measures are concerned.

Optimal power flow via interior-exterior method

This paper presents a new approach to the resolution of the Optimal Power Flow problem. In this approach the inequality constraints are treated by the Modified Barrier and Primal-Dual Logarithmic Barrier methods. The inequality constraints are transformed into equalities by introducing positive auxiliary variables, which are perturbed by the barrier parameter. A Lagrangian function is associated with the modified problem. The first-order necessary conditions are applied to the Lagrangian, generating a nonlinear system which is solved by Newton's method. The perturbation of the auxiliary variables results in an expansion of the feasible set of the original problem, allowing the limits of the inequality constraints to be reached. Numerical tests on the Brazilian CESP and South-Southeast systems and a comparative test indicated that the new approach efficiently resolves of the Optimal Power Flow problem. © 2007 IEEE.

Voltage regulator modeling for the three-phase power flow in distribution networks

In this paper a three-phase power flow for electrical distribution systems considering different models of voltage regulators is presented. A voltage regulator (VR) is an equipment that maintains the voltage level in a predefined value in a distribution line in spite of the load variations within its nominal power. Three different types of connections are analyzed: 1) wye-connected regulators, 2) open delta-connected regulators and 3) closed delta-connected regulators. To calculate the power flow, the three-phase backward/forward sweep algorithm is used. The methodology is tested on the IEEE 34 bus distribution system. ©2008 IEEE.

A Numerical Study of the FENE-CR Model Applied to a Jet Flow Problem

The FENE-CR model is investigated through a numerical algorithm to simulate the time-dependent moving free surface flow produced by a jet impinging on a flat surface. The objective is to demonstrate that by increasing the extensibility parameter L, the numerical solutions converge to the solutions obtained with the Oldroyd-B model. The governing equations are solved by an established free surface flow solver based on the finite difference and marker-and-cell methods. Numerical predictions of the extensional viscosity obtained with several values of the parameter L are presented. The results show that if the extensibility parameter L is sufficiently large then the extensional viscosities obtained with the FENE-CR model approximate the corresponding Oldroyd-B viscosity. Moreover, the flow from a jet impinging on a flat surface is simulated with various values of the extensibility parameter L and the fluid flow visualizations display convergence to the Oldroyd-B jet flow results.

Optimal reactive power flow via the modified barrier Lagrangian function approach

A new approach called the Modified Barrier Lagrangian Function (MBLF) to solve the Optimal Reactive Power Flow problem is presented. In this approach, the inequality constraints are treated by the Modified Barrier Function (MBF) method, which has a finite convergence property: i.e. the optimal solution in the MBF method can actually be in the bound of the feasible set. Hence, the inequality constraints can be precisely equal to zero. Another property of the MBF method is that the barrier parameter does not need to be driven to zero to attain the solution. Therefore, the conditioning of the involved Hessian matrix is greatly enhanced. In order to show this, a comparative analysis of the numeric conditioning of the Hessian matrix of the MBLF approach, by the decomposition in singular values, is carried out. The feasibility of the proposed approach is also demonstrated with comparative tests to Interior Point Method (IPM) using various IEEE test systems and two networks derived from Brazilian generation/transmission system. The results show that the MBLF method is computationally more attractive than the IPM in terms of speed, number of iterations and numerical conditioning. (C) 2011 Elsevier B.V. All rights reserved.

Saturated flow boiling heat transfer and critical heat flux in small horizontal flattened tubes

This paper presents experimental results for flow boiling heat transfer coefficient and critical heat flux (CHF) in small flattened tubes. The tested flattened tubes have the same equivalent internal diameter of 2.2 mm, but different aspect height/width ratios (H/W) of 1/4, 1/2, 2 and 4. The experimental data were compared against results for circular tubes using R134a and R245fa as working fluids at a nominal saturation temperature of 31 degrees C. For mass velocities higher than 200 kg/m(2)s, the flattened and circular tubes presented similar heat transfer coefficients. Such a behavior is related to the fact that stratification effects are negligible under conditions of higher mass velocities. Heat transfer correlations from the literature, usually developed using only circular-channel experimental data, predicted the flattened tube results for mass velocities higher than 200 kg/m(2)s with mean absolute error lower than 20% using the equivalent diameter to account for the geometry effect. Similarly, the critical heat flux results were found to be independent of the tube aspect ratio when the same equivalent length was kept. Equivalent length is a new parameter which takes into account the channel heat transfer area. The CHF correlations for round tubes predicted the flattened tube data relatively well when using the equivalent diameter and length. Furthermore, a new proposed CHF correlation predicted the present flattened tube data with a mean absolute error of 5%. (C) 2012 Elsevier Ltd. All rights reserved.

Analysis of a new cross flow heat exchanger flow arrangement - Extension to several rows

The present paper presents a theoretical analysis of a cross flow heat exchanger with a new flow arrangement comprehending several tube rows. The thermal performance of the proposed flow arrangement is compared with the thermal performance of a typical counter cross flow arrangement that is used in chemical, refrigeration, automotive and air conditioning industries. The thermal performance comparison has been performed in terms of the following parameters: heat exchanger effectiveness and efficiency, dimensionless entropy generation, entransy dissipation number, and dimensionless local temperature differences. It is also shown that the uniformity of the temperature difference field leads to a higher thermal performance of the heat exchanger. In the present case this is accomplished thorough a different organization of the in-tube fluid circuits in the heat exchanger. The relation between the recently introduced "entransy dissipation number" and the conventional thermal effectiveness has been obtained in terms of the "number of transfer units". A case study has been solved to quantitatively to obtain the temperature difference distribution over two rows units involving the proposed arrangement and the counter cross flow one. It has been shown that the proposed arrangement presents better thermal performance regardless the comparison parameter. (C) 2012 Elsevier Masson SAS. All rights reserved.

Deviation from quark number scaling of the anisotropy parameter v(2) of pions, kaons, and protons in Au+Au collisions at root s(NN)=200 GeV

Measurements of the anisotropy parameter v(2) of identified hadrons (pions, kaons, and protons) as a function of centrality, transverse momentum p(T), and transverse kinetic energy KET at midrapidity (vertical bar eta vertical bar < 0.35) in Au + Au collisions at root s(N N) = 200 GeV are presented. Pions and protons are identified up to p(T) = 6 GeV/c, and kaons up to p(T) = 4 GeV/c, by combining information from time-of-flight and aerogel Cerenkov detectors in the PHENIX Experiment. The scaling of v(2) with the number of valence quarks (n(q)) has been studied in different centrality bins as a function of transverse momentum and transverse kinetic energy. A deviation from previously observed quark-number scaling is observed at large values of KET/n(q) in noncentral Au + Au collisions (20-60%), but this scaling remains valid in central collisions (0-10%).

Echocardiographic and hemodynamic determinants of right coronary artery flow reserve and phasic flow pattern in advanced non-ischemic cardiomyopathy

Abstract Background In patients with advanced non-ischemic cardiomyopathy (NIC), right-sided cardiac disturbances has prognostic implications. Right coronary artery (RCA) flow pattern and flow reserve (CFR) are not well known in this setting. The purpose of this study was to assess, in human advanced NIC, the RCA phasic flow pattern and CFR, also under right-sided cardiac disturbances, and compare with left coronary circulation. As well as to investigate any correlation between the cardiac structural, mechanical and hemodynamic parameters with RCA phasic flow pattern or CFR. Methods Twenty four patients with dilated severe NIC were evaluated non-invasively, even by echocardiography, and also by cardiac catheterization, inclusive with Swan-Ganz catheter. Intracoronary Doppler (Flowire) data was obtained in RCA and left anterior descendent coronary artery (LAD) before and after adenosine. Resting RCA phasic pattern (diastolic/systolic) was compared between subgroups with and without pulmonary hypertension, and with and without right ventricular (RV) dysfunction; and also with LAD. RCA-CFR was compared with LAD, as well as in those subgroups. Pearson's correlation analysis was accomplished among echocardiographic (including LV fractional shortening, mass index, end systolic wall stress) more hemodynamic parameters with RCA phasic flow pattern or RCA-CFR. Results LV fractional shortening and end diastolic diameter were 15.3 ± 3.5 % and 69.4 ± 12.2 mm. Resting RCA phasic pattern had no difference comparing subgroups with vs. without pulmonary hypertension (1.45 vs. 1.29, p = NS) either with vs. without RV dysfunction (1.47 vs. 1.23, p = NS); RCA vs. LAD was 1.35 vs. 2.85 (p < 0.001). It had no significant correlation among any cardiac mechanical or hemodynamic parameter with RCA-CFR or RCA flow pattern. RCA-CFR had no difference compared with LAD (3.38 vs. 3.34, p = NS), as well as in pulmonary hypertension (3.09 vs. 3.10, p = NS) either in RV dysfunction (3.06 vs. 3.22, p = NS) subgroups. Conclusion In patients with chronic advanced NIC, RCA phasic flow pattern has a mild diastolic predominance, less marked than in LAD, with no effects from pulmonary artery hypertension or RV dysfunction. There is no significant correlation between any cardiac mechanical-structural or hemodynamic parameter with RCA-CFR or RCA phasic flow pattern. RCA flow reserve is still similar to LAD, independently of those right-sided cardiac disturbances.

Reliable estimation of dense optical flow fields with large displacements

[EN] In this paper we show that a classic optical flow technique by Nagel and Enkelmann can be regarded as an early anisotropic diffusion method with a diffusion tensor. We introduce three improvements into the model formulation that avoid inconsistencies caused by centering the brightness term and the smoothness term in different images use a linear scale-space focusing strategy from coarse to fine scales for avoiding convergence to physically irrelevant local minima, and create an energy functional that is invariant under linear brightness changes.  Applying a gradient descent method to the resulting energy functional leads to a system of diffusion-reaction equations. We prove that this system has a unique solution under realistic assumptions on the initial data, and we present an efficient linear implicit numerical scheme in detail. Our method creates flow fields with 100% density over the entire image domain, it is robust under a large range of parameter variations, and it can recover displacement fields that are far beyond the typical one-pixel limits which are characteristic for many differential methods for determining optical flow. We show that it performs better than the classic optical flow methods with 100%  density that are evaluated by Barron et al. (1994). Our software is available from the Internet.

Discontinuity preserving in TV-L1 optical flow methods

[EN] We analyze the discontinuity preserving problem in TV-L1 optical flow methods. This type of methods typically creates rounded effects at flow boundaries, which usually do not coincide with object contours. A simple strategy to overcome this problem consists in inhibiting the diffusion at high image gradients. In this work, we first introduce a general framework for TV regularizers in optical flow and relate it with some standard approaches. Our survey takes into account several methods that use decreasing functions for mitigating the diffusion at image contours. Consequently, this kind of strategies may produce instabilities in the estimation of the optical flows. Hence, we study the problem of instabilities and show that it actually arises from an ill-posed formulation. From this study, it is possible to come across with different schemes to solve this problem. One of these consists in separating the pure TV process from the mitigating strategy. This has been used in another work and we demonstrate here that it has a good performance. Furthermore, we propose two alternatives to avoid the instability problems: (i) we study a fully automatic approach that solves the problem based on the information of the whole image; (ii) we derive a semi-automatic approach that takes into account the image gradients in a close neighborhood adapting the parameter in each position. In the experimental results, we present a detailed study and comparison between the different alternatives. These methods provide very good results, especially for sequences with a few dominant gradients. Additionally, a surprising effect of these approaches is that they can cope with occlusions. This can be easily achieved by using strong regularizations and high penalizations at image contours.