A modified Primal-Dual Logarithmic-Barrier Method for solving the Optimal Power Flow problem with discrete and continuous control variables

The aim of solving the Optimal Power Flow problem is to determine the optimal state of an electric power transmission system, that is, the voltage magnitude and phase angles and the tap ratios of the transformers that optimize the performance of a given system, while satisfying its physical and operating constraints. The Optimal Power Flow problem is modeled as a large-scale mixed-discrete nonlinear programming problem. This paper proposes a method for handling the discrete variables of the Optimal Power Flow problem. A penalty function is presented. Due to the inclusion of the penalty function into the objective function, a sequence of nonlinear programming problems with only continuous variables is obtained and the solutions of these problems converge to a solution of the mixed problem. The obtained nonlinear programming problems are solved by a Primal-Dual Logarithmic-Barrier Method. Numerical tests using the IEEE 14, 30, 118 and 300-Bus test systems indicate that the method is efficient. (C) 2012 Elsevier B.V. All rights reserved.

Modeling and simulation of severe slugging with mass transfer effects

A mathematical model and numerical simulations are presented to investigate the dynamics of gas, oil and water flow in a pipeline-riser system. The pipeline is modeled as a lumped parameter system and considers two switchable states: one in which the gas is able to penetrate into the riser and another in which there is a liquid accumulation front, preventing the gas from penetrating the riser. The riser model considers a distributed parameter system, in which movable nodes are used to evaluate local conditions along the subsystem. Mass transfer effects are modeled by using a black oil approximation. The model predicts the liquid penetration length in the pipeline and the liquid level in the riser, so it is possible to determine which type of severe slugging occurs in the system. The method of characteristics is used to simplify the differentiation of the resulting hyperbolic system of equations. The equations are discretized and integrated using an implicit method with a predictor-corrector scheme for the treatment of the nonlinearities. Simulations corresponding to severe slugging conditions are presented and compared to results obtained with OLGA computer code, showing a very good agreement. A description of the types of severe slugging for the three-phase flow of gas, oil and water in a pipeline-riser system with mass transfer effects are presented, as well as a stability map. (C) 2011 Elsevier Ltd. All rights reserved.

Error Upper Bounds for a Computational Method for Nonlinear Boundary and Initial-Value Problems

This work develops a computational approach for boundary and initial-value problems by using operational matrices, in order to run an evolutive process in a Hilbert space. Besides, upper bounds for errors in the solutions and in their derivatives can be estimated providing accuracy measures.

Resting heart rate as a predictor of metabolic dysfunctions in obese children and adolescents

Background: Recent studies have identified that a higher resting heart rate (RHR) is associated with elevated blood pressure, independent of body fatness, age and ethnicity. However, it is still unclear whether RHR can also be applied as a screening for other risk factors, such as hyperglycemia and dyslipidemia. Thus, the purpose of the presented study was to analyze the association between RHR, lipid profile and fasting glucose in obese children and adolescents. Methods: The sample was composed of 180 obese children and adolescents, aged between 7-16 years. Whole-body and segmental body composition were estimated by Dual-energy X-ray absorptiometry. Resting heart rate (RHR) was measured by heart rate monitors. The fasting blood samples were analyzed for serum triglycerides, total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and glucose, using the colorimetric method. Results: Fasting glucose, TC, triglycerides, HDL-C, LDL-C and RHR were similar in both genders. The group of obese subjects with a higher RHR presented, at a lower age, higher triglycerides and TC. There was a significant relationship between RHR, triglycerides and TC. In the multivariate model, triglycerides and TC maintained a significant relationship with RHR independent of age, gender, general and trunk adiposity. The ROC curve indicated that RHR has a high potential for screening elevated total cholesterol and triglycerides as well as dyslipidemia. Conclusion: Elevated RHR has the potential to identify subjects at an increased risk of atherosclerosis development.

Hirota method for oblique solitons in two-dimensional supersonic nonlinear Schrodinger flow

In a previous work El et al. (2006) [1] exact stable oblique soliton solutions were revealed in two-dimensional nonlinear Schrodinger flow. In this work we show that single soliton solution can be expressed within the Hirota bilinear formalism. An attempt to build two-soliton solutions shows that the system is "close" to integrability provided that the angle between the solitons is small and/or we are in the hypersonic limit. (C) 2012 Elsevier B.V. All rights reserved.

Analysis Methodology for Vortex-Induced Motion (VIM) of a Monocolumn Platform Applying the Hilbert-Huang Transform Method

Vortex-induced motion (VIM) is a highly nonlinear dynamic phenomenon. Usual spectral analysis methods, using the Fourier transform, rely on the hypotheses of linear and stationary dynamics. A method to treat nonstationary signals that emerge from nonlinear systems is denoted Hilbert-Huang transform (HHT) method. The development of an analysis methodology to study the VIM of a monocolumn production, storage, and offloading system using HHT is presented. The purposes of the present methodology are to improve the statistics analysis of VIM. The results showed to be comparable to results obtained from a traditional analysis (mean of the 10% highest peaks) particularly for the motions in the transverse direction, although the difference between the results from the traditional analysis for the motions in the in-line direction showed a difference of around 25%. The results from the HHT analysis are more reliable than the traditional ones, owing to the larger number of points to calculate the statistics characteristics. These results may be used to design risers and mooring lines, as well as to obtain VIM parameters to calibrate numerical predictions. [DOI: 10.1115/1.4003493]

Retrospective Clinical Study of 988 Dual Acid-Etched Implants Placed in Grafted and Native Bone for Single-Tooth Replacement

Purpose: To evaluate the influence of sex, implant characteristics, and bone grafting on the survival rate of dual acid-etched (DAE) implants. Materials and Methods: Patients treated with internal-hex DAE implants for single-tooth replacement in a military dental clinic between January 2005 and December 2010 were included in this study. Clinical data related to implant characteristics, implant location, presence of grafted bone, and implant failures were collected. The primary outcome was implant loss. The survival rate was analyzed using the Kaplan-Meier method. Cox regression modeling was used to determine which factors would predict implant failure. Results: DAE implants were evaluated in a total of 988 patients (80.3% men). Twenty-four (2.4%) implants failed, most were cylindric (54.2%) with regular platforms (70.8%) and were 10 mm long (58.3%). The failure rate was 2.4% for the anterior maxilla, 3.3% for the posterior maxilla, 1.6% for the anterior mandible, and 2.0% for posterior mandible. The cumulative survival rate was 97.6%. The failure rate was 8.8% in implants placed after sinus augmentation, 7.3% in bone block-grafted areas, and 1.6% in native bone. Based on multivariable analysis (Cox regression), sinus augmentation and bone block grafting had a statistically significant effect on implant failure; the hazard ratios were 5.5 and 4.6, respectively. Conclusion: The results revealed that DAE implants had high survival rates, and no influence of sex, location, shape, diameter, or length on failure rates could be observed. However, a significant association was observed between failure and presence of bone graft in the implant area. Int J Oral Maxillofac Implants 2012;27:1243-1248

A generalized finite element method with global-local enrichment functions for confined plasticity problems

The main feature of partition of unity methods such as the generalized or extended finite element method is their ability of utilizing a priori knowledge about the solution of a problem in the form of enrichment functions. However, analytical derivation of enrichment functions with good approximation properties is mostly limited to two-dimensional linear problems. This paper presents a procedure to numerically generate proper enrichment functions for three-dimensional problems with confined plasticity where plastic evolution is gradual. This procedure involves the solution of boundary value problems around local regions exhibiting nonlinear behavior and the enrichment of the global solution space with the local solutions through the partition of unity method framework. This approach can produce accurate nonlinear solutions with a reduced computational cost compared to standard finite element methods since computationally intensive nonlinear iterations can be performed on coarse global meshes after the creation of enrichment functions properly describing localized nonlinear behavior. Several three-dimensional nonlinear problems based on the rate-independent J (2) plasticity theory with isotropic hardening are solved using the proposed procedure to demonstrate its robustness, accuracy and computational efficiency.

A nonlinear elasticity phantom containing spherical inclusions

The strain image contrast of some in vivo breast lesions changes with increasing applied load. This change is attributed to differences in the nonlinear elastic properties of the constituent tissues suggesting some potential to help classify breast diseases by their nonlinear elastic properties. A phantom with inclusions and long-term stability is desired to serve as a test bed for nonlinear elasticity imaging method development, testing, etc. This study reports a phantom designed to investigate nonlinear elastic properties with ultrasound elastographic techniques. The phantom contains four spherical inclusions and was manufactured from a mixture of gelatin, agar and oil. The phantom background and each of the inclusions have distinct Young's modulus and nonlinear mechanical behavior. This phantom was subjected to large deformations (up to 20%) while scanning with ultrasound, and changes in strain image contrast and contrast-to-noise ratio between inclusion and background, as a function of applied deformation, were investigated. The changes in contrast over a large deformation range predicted by the finite element analysis (FEA) were consistent with those experimentally observed. Therefore, the paper reports a procedure for making phantoms with predictable nonlinear behavior, based on independent measurements of the constituent materials, and shows that the resulting strain images (e. g., strain contrast) agree with that predicted with nonlinear FEA.

Augmented Lagrangian method with nonmonotone penalty parameters for constrained optimization

At each outer iteration of standard Augmented Lagrangian methods one tries to solve a box-constrained optimization problem with some prescribed tolerance. In the continuous world, using exact arithmetic, this subproblem is always solvable. Therefore, the possibility of finishing the subproblem resolution without satisfying the theoretical stopping conditions is not contemplated in usual convergence theories. However, in practice, one might not be able to solve the subproblem up to the required precision. This may be due to different reasons. One of them is that the presence of an excessively large penalty parameter could impair the performance of the box-constraint optimization solver. In this paper a practical strategy for decreasing the penalty parameter in situations like the one mentioned above is proposed. More generally, the different decisions that may be taken when, in practice, one is not able to solve the Augmented Lagrangian subproblem will be discussed. As a result, an improved Augmented Lagrangian method is presented, which takes into account numerical difficulties in a satisfactory way, preserving suitable convergence theory. Numerical experiments are presented involving all the CUTEr collection test problems.

Transverse pseudo-nonlinear effects measured in solid-state laser materials using a sensitive time-resolved technique

We present a detailed study of the Baryscan technique, a new efficient alternative to the widespread Z-scan technique which has been demonstrated [Opt. Lett. 36:8, 2011] to reach among the highest sensitivity levels. This method is based upon the measurement of optical nonlinearities by means of beam centroid displacements with a position sensitive detector and is able to deal with any kind of lensing effect. This technique is applied here to measure pump-induced electronic refractive index changes (population lens), which can be discriminated from parasitic thermal effects by using a time-resolved Baryscan experiment. This method is validated by evaluating the polarizability variation at the origin of the population lens observed in the reference Cr3+:GSGG laser material.

A Validated Chromatographic Method for the Determination of Flavonoids in Copaifera langsdorffii by HPLC

Hydroethanolic extracts of C. langsdorffii leaves have therapeutic potential. This work reports a validated chromatographic method for the quantification of polar compounds in the hydroethanolic extract of C. langsdorffii leaves. A reliable HPLC method was developed using two monolithic columns linked in series (100 x 4.6 mm - C-18), with nonlinear gradient elution, and UV detection set at 257 nm. A procedure for the extraction of flavonols was also developed, which involved the use of 70% aqueous ethanol and the addition of benzophenone as the internal standard. The developed method led to a good detection response as the values for linearity were between 10.3 and 1000 mu g/mL, and those for recovery between 84.2 and 111.1%. The detection limit ranged from 0.02 to 1.70 mu g/mL and the quantitation limit from 0.07 to 5.1 mu g/mL, with a maximum RSD of 5.24%. Five compounds, rutin, quercetin-3-O-alpha-L-rhamnopyranoside, kaempferol-3-O-alpha-L-rhamnopyranoside, quercetin and kaempferol, were quantified. This method could, therefore, be used for the quality control of hydroethanolic extracts of Copaifera leaves and their cosmetic and pharmaceutical products.

The boundedness of penalty parameters in an augmented Lagrangian method with constrained subproblems

Augmented Lagrangian methods are effective tools for solving large-scale nonlinear programming problems. At each outer iteration, a minimization subproblem with simple constraints, whose objective function depends on updated Lagrange multipliers and penalty parameters, is approximately solved. When the penalty parameter becomes very large, solving the subproblem becomes difficult; therefore, the effectiveness of this approach is associated with the boundedness of the penalty parameters. In this paper, it is proved that under more natural assumptions than the ones employed until now, penalty parameters are bounded. For proving the new boundedness result, the original algorithm has been slightly modified. Numerical consequences of the modifications are discussed and computational experiments are presented.

Analysis of nonlinear dynamics of vocal folds using high-speed video observation and biomechanical modeling

Primary voice production occurs in the larynx through vibrational movements carried out by vocal folds. However, many problems can affect this complex system resulting in voice disorders. In this context, time-frequency-shape analysis based on embedding phase space plots and nonlinear dynamics methods have been used to evaluate the vocal fold dynamics during phonation. For this purpose, the present work used high-speed video to record the vocal fold movements of three subjects and extract the glottal area time series using an image segmentation algorithm. This signal is used for an optimization method which combines genetic algorithms and a quasi-Newton method to optimize the parameters of a biomechanical model of vocal folds based on lumped elements (masses, springs and dampers). After optimization, this model is capable of simulating the dynamics of recorded vocal folds and their glottal pulse. Bifurcation diagrams and phase space analysis were used to evaluate the behavior of this deterministic system in different circumstances. The results showed that this methodology can be used to extract some physiological parameters of vocal folds and reproduce some complex behaviors of these structures contributing to the scientific and clinical evaluation of voice production. (C) 2010 Elsevier Inc. All rights reserved.

A new approach for modeling and control of nonlinear systems via norm-bounded linear differential inclusions

A systematic approach to model nonlinear systems using norm-bounded linear differential inclusions (NLDIs) is proposed in this paper. The resulting NLDI model is suitable for the application of linear control design techniques and, therefore, it is possible to fulfill certain specifications for the underlying nonlinear system, within an operating region of interest in the state-space, using a linear controller designed for this NLDI model. Hence, a procedure to design a dynamic output feedback controller for the NLDI model is also proposed in this paper. One of the main contributions of the proposed modeling and control approach is the use of the mean-value theorem to represent the nonlinear system by a linear parameter-varying model, which is then mapped into a polytopic linear differential inclusion (PLDI) within the region of interest. To avoid the combinatorial problem that is inherent of polytopic models for medium- and large-sized systems, the PLDI is transformed into an NLDI, and the whole process is carried out ensuring that all trajectories of the underlying nonlinear system are also trajectories of the resulting NLDI within the operating region of interest. Furthermore, it is also possible to choose a particular structure for the NLDI parameters to reduce the conservatism in the representation of the nonlinear system by the NLDI model, and this feature is also one important contribution of this paper. Once the NLDI representation of the nonlinear system is obtained, the paper proposes the application of a linear control design method to this representation. The design is based on quadratic Lyapunov functions and formulated as search problem over a set of bilinear matrix inequalities (BMIs), which is solved using a two-step separation procedure that maps the BMIs into a set of corresponding linear matrix inequalities. Two numerical examples are given to demonstrate the effectiveness of the proposed approach.