Método da continuação aplicado na análise de contingência de linhas de transmissão

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Proposição de uma técnica de parametrização geométrica para o fluxo de carga continuado

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Parametric resonances in a base-excited double pendulum

Two parametrically-induced phenomena are addressed in the context of a double pendulum subject to a vertical base excitation. First, the parametric resonances that cause the stable downward vertical equilibrium to bifurcate into large-amplitude periodic solutions are investigated extensively. Then the stabilization of the unstable upward equilibrium states through the parametric action of the high-frequency base motion is documented in the experiments and in the simulations. It is shown that there is a region in the plane of the excitation frequency and amplitude where all four unstable equilibrium states can be stabilized simultaneously in the double pendulum. The parametric resonances of the two modes of the base-excited double pendulum are studied both theoretically and experimentally. The transition curves (i.e., boundaries of the dynamic instability regions) are constructed asymptotically via the method of multiple scales including higher-order effects. The bifurcations characterizing the transitions from the trivial equilibrium to the periodic solutions are computed by either continuation methods and or by time integration and compared with the theoretical and experimental results.

A robust numerical framework for the analysis of material failure of fibre reinforced soft tissue

In this work, robustness and stability of continuum damage models applied to material failure in soft tissues are addressed. In the implicit damage models equipped with softening, the presence of negative eigenvalues in the tangent elemental matrix degrades the condition number of the global matrix, leading to a reduction of the computational performance of the numerical model. Two strategies have been adapted from literature to improve the aforementioned computational performance degradation: the IMPL-EX integration scheme [Oliver,2006], which renders the elemental matrix contribution definite positive, and arclength-type continuation methods [Carrera,1994], which allow to capture the unstable softening branch in brittle ruptures. The IMPL-EX integration scheme has as a major drawback the need to use small time steps to keep numerical error below an acceptable value. A convergence study, limiting the maximum allowed increment of internal variables in the damage model, is presented. Finally, numerical simulation of failure problems with fibre reinforced materials illustrates the performance of the adopted methodology.

A robust numerical framework to the analysis of material failure of fibre reinforced soft tissue

In this work, robustness and stability of continuum damage models applied to material failure in soft tissues are addressed. In the implicit damage models equipped with softening, the presence of negative eigenvalues in the tangent elemental matrix degrades the condition number of the global matrix, leading to a reduction of the computational performance of the numerical model. Two strategies have been adapted from literature to improve the aforementioned computational performance degradation: the IMPL-EX integration scheme [Oliver,2006], which renders the elemental matrix contribution definite positive, and arclength-type continuation methods [Carrera,1994], which allow to capture the unstable softening branch in brittle ruptures. The IMPL-EX integration scheme has as a major drawback the need to use small time steps to keep numerical error below an acceptable value. A convergence study, limiting the maximum allowed increment of internal variables in the damage model, is presented. Finally, numerical simulation of failure problems with fibre reinforced materials illustrates the performance of the adopted methodology.

Dynamics of coupled modes in two-section semiconductor lasers with saturable absorption

Semiconductor lasers have the potential to address a number of critical applications in advanced telecommunications and signal processing. These include applications that require pulsed output that can be obtained from self-pulsing and mode-locked states of two-section devices with saturable absorption. Many modern applications place stringent performance requirements on the laser source, and a thorough understanding of the physical mechanisms underlying these pulsed modes of operation is therefore highly desirable. In this thesis, we present experimental measurements and numerical simulations of a variety of self-pulsation phenomena in two-section semiconductor lasers with saturable absorption. Our theoretical and numerical results will be based on rate equations for the field intensities and the carrier densities in the two sections of the device, and we establish typical parameter ranges and assess the level of agreement with experiment that can be expected from our models. For each of the physical examples that we consider, our model parameters are consistent with the physical net gain and absorption of the studied devices. Following our introductory chapter, the first system that we consider is a two-section Fabry-Pérot laser. This example serves to introduce our method for obtaining model parameters from the measured material dispersion, and it also allows us to present a detailed discussion of the bifurcation structure that governs the appearance of selfpulsations in two-section devices. In the following two chapters, we present two distinct examples of experimental measurements from dual-mode two-section devices. In each case we have found that single mode self-pulsations evolve into complex coupled dualmode states following a characteristic series of bifurcations. We present optical and mode resolved power spectra as well as a series of characteristic intensity time traces illustrating this progression for each example. Using the results from our study of a twosection Fabry-Pérot device as a guide, we find physically appropriate model parameters that provide qualitative agreement with our experimental results. We highlight the role played by material dispersion and the underlying single mode self-pulsing orbits in determining the observed dynamics, and we use numerical continuation methods to provide a global picture of the governing bifurcation structure. In our concluding chapter we summarise our work, and we discuss how the presented results can inform the development of optimised mode-locked lasers for performance applications in integrated optics.

Pseudospectral methods and numerical continuation for the analysis os structured population models.

191 p.

Univariate, bivariate, and multivariate methods in corpus-based lexicography : A study of synonymy

In this dissertation, I present an overall methodological framework for studying linguistic alternations, focusing specifically on lexical variation in denoting a single meaning, that is, synonymy. As the practical example, I employ the synonymous set of the four most common Finnish verbs denoting THINK, namely ajatella, miettiä, pohtia and harkita ‘think, reflect, ponder, consider’. As a continuation to previous work, I describe in considerable detail the extension of statistical methods from dichotomous linguistic settings (e.g., Gries 2003; Bresnan et al. 2007) to polytomous ones, that is, concerning more than two possible alternative outcomes. The applied statistical methods are arranged into a succession of stages with increasing complexity, proceeding from univariate via bivariate to multivariate techniques in the end. As the central multivariate method, I argue for the use of polytomous logistic regression and demonstrate its practical implementation to the studied phenomenon, thus extending the work by Bresnan et al. (2007), who applied simple (binary) logistic regression to a dichotomous structural alternation in English. The results of the various statistical analyses confirm that a wide range of contextual features across different categories are indeed associated with the use and selection of the selected think lexemes; however, a substantial part of these features are not exemplified in current Finnish lexicographical descriptions. The multivariate analysis results indicate that the semantic classifications of syntactic argument types are on the average the most distinctive feature category, followed by overall semantic characterizations of the verb chains, and then syntactic argument types alone, with morphological features pertaining to the verb chain and extra-linguistic features relegated to the last position. In terms of overall performance of the multivariate analysis and modeling, the prediction accuracy seems to reach a ceiling at a Recall rate of roughly two-thirds of the sentences in the research corpus. The analysis of these results suggests a limit to what can be explained and determined within the immediate sentential context and applying the conventional descriptive and analytical apparatus based on currently available linguistic theories and models. The results also support Bresnan’s (2007) and others’ (e.g., Bod et al. 2003) probabilistic view of the relationship between linguistic usage and the underlying linguistic system, in which only a minority of linguistic choices are categorical, given the known context – represented as a feature cluster – that can be analytically grasped and identified. Instead, most contexts exhibit degrees of variation as to their outcomes, resulting in proportionate choices over longer stretches of usage in texts or speech.

Continuation of limit cycles near saddle homoclinic points using splines in phase space

We present a new algorithm for continuation of limit cycles of autonomous systems as a system parameter is varied. The algorithm works in phase space with an ordered set of points on the limit cycle, along with spline interpolation. Currently popular algorithms in bifurcation analysis packages compute time-domain approximations of limit cycles using either shooting or collocation. The present approach seems useful for continuation near saddle homoclinic points, where it encounters a corner while time-domain methods essentially encounter a discontinuity (a relatively short period of rapid variation). Other phase space-based algorithms use rescaled arclength in place of time, but subsequently resemble the time-domain methods. Compared to these, we introduce additional freedom through a variable stretching of arclength based on local curvature, through the use of an auxiliary index-based variable. Several numerical examples are presented. Comparisons with results from the popular package, MATCONT, are favorable close to saddle homoclinic points.

Diagonal forms, linear algebraic methods and Ramsey-type problems

This thesis focuses mainly on linear algebraic aspects of combinatorics. Let N_t(H) be an incidence matrix with edges versus all subhypergraphs of a complete hypergraph that are isomorphic to H. Richard M. Wilson and the author find the general formula for the Smith normal form or diagonal form of N_t(H) for all simple graphs H and for a very general class of t-uniform hypergraphs H.

As a continuation, the author determines the formula for diagonal forms of integer matrices obtained from other combinatorial structures, including incidence matrices for subgraphs of a complete bipartite graph and inclusion matrices for multisets.

One major application of diagonal forms is in zero-sum Ramsey theory. For instance, Caro's results in zero-sum Ramsey numbers for graphs and Caro and Yuster's results in zero-sum bipartite Ramsey numbers can be reproduced. These results are further generalized to t-uniform hypergraphs. Other applications include signed bipartite graph designs.

Research results on some other problems are also included in this thesis, such as a Ramsey-type problem on equipartitions, Hartman's conjecture on large sets of designs and a matroid theory problem proposed by Welsh.

Fast numerical methods for mixed, singular Helmholtz boundary value problems and Laplace eigenvalue problems - with applications to antenna design, sloshing, electromagnetic scattering and spectral geometry

This thesis presents a novel class of algorithms for the solution of scattering and eigenvalue problems on general two-dimensional domains under a variety of boundary conditions, including non-smooth domains and certain "Zaremba" boundary conditions - for which Dirichlet and Neumann conditions are specified on various portions of the domain boundary. The theoretical basis of the methods for the Zaremba problems on smooth domains concern detailed information, which is put forth for the first time in this thesis, about the singularity structure of solutions of the Laplace operator under boundary conditions of Zaremba type. The new methods, which are based on use of Green functions and integral equations, incorporate a number of algorithmic innovations, including a fast and robust eigenvalue-search algorithm, use of the Fourier Continuation method for regularization of all smooth-domain Zaremba singularities, and newly derived quadrature rules which give rise to high-order convergence even around singular points for the Zaremba problem. The resulting algorithms enjoy high-order convergence, and they can tackle a variety of elliptic problems under general boundary conditions, including, for example, eigenvalue problems, scattering problems, and, in particular, eigenfunction expansion for time-domain problems in non-separable physical domains with mixed boundary conditions.

Parameterized fast decoupled power flow methods for obtaining the maximum loading point of power systems. Part II. Performance evaluation

The parameterized fast decoupled power flow (PFDPF), versions XB and BX, using either theta or V as a parameter have been proposed by the authors in Part I of this paper. The use of reactive power injection of a selected PVbus (Q(PV)) as the continuation parameter for the computation of the maximum loading point (MLP) was also investigated. In this paper, the proposed versions obtained only with small modifications of the conventional one are used for the computation of the MLP of IEEE test systems (14, 30, 57 and 118 buses). These new versions are compared to each other with the purpose of pointing out their features, as well as the influence of reactive power and transformer tap limits. The results obtained with the new approaches are presented and discussed. The results show that the characteristics of the conventional FDPF method are enhanced and the region of convergence around the singular solution is enlarged. In addition, it is shown that these versions can be switched during the tracing process in order to efficiently determine all the PV curve points with few iterations. A trivial secant predictor, the modified zero-order polynomial, which uses the current solution and a fixed increment in the parameter (V, theta, or mu) as an estimate for the next solution, is used for the predictor step. (C) 2003 Elsevier B.V. All rights reserved.

Continuation fast decoupled power flow with secant predictor

The conventional Newton and fast decoupled power flow methods are considered inadequate for obtaining the maximum loading point of power systems due to ill-conditioning problems at and near this critical point. At this point, the Jacobian matrix of the Newton method becomes singular. In addition, it is widely accepted that the P-V and Q-theta decoupling assumptions made for the fast decoupled power flow formulation no longer hold. However, in this paper, it is presented a new fast decoupled power flow that becomes adequate for the computation of the maximum loading point by simply using the reactive power injection of a selected PV bus as a continuation parameter. Besides, fast decoupled methods using V and 0 as parameters and a secant predictor are also presented. These new versions are compared to each other with the purpose of pointing out their features, as well as the influence of reactive power and transformer tap limits. The results obtained for the IEEE systems (14 and 118 buses) show that the characteristics of the conventional method are enhanced and the region of convergence around the singular solution is enlarged.

Initiation and continuation of randomized trials after the publication of a trial stopped early for benefit asking the same study question: STOPIT-3 study design

BACKGROUND Randomized control trials (RCTs) stopped early for benefit (truncated RCTs) are increasingly common and, on average, overestimate the relative magnitude of benefit by approximately 30%. Investigators stop trials early when they consider it is no longer ethical to enroll patients in a control group. The goal of this systematic review is to determine how investigators of ongoing or planned RCTs respond to the publication of a truncated RCT addressing a similar question. METHODS/DESIGN We will conduct systematic reviews to update the searches of 210 truncated RCTs to identify similar trials ongoing at the time of publication, or started subsequently, to the truncated trials ('subsequent RCTs'). Reviewers will determine in duplicate the similarity between the truncated and subsequent trials. We will analyze the epidemiology, distribution, and predictors of subsequent RCTs. We will also contact authors of subsequent trials to determine reasons for beginning, continuing, or prematurely discontinuing their own trials, and the extent to which they rely on the estimates from truncated trials. DISCUSSION To the extent that investigators begin or continue subsequent trials they implicitly disagree with the decision to stop the truncated RCT because of an ethical mandate to administer the experimental treatment. The results of this study will help guide future decisions about when to stop RCTs early for benefit.

Bevacizumab continuation versus no continuation after first-line chemotherapy plus bevacizumab in patients with metastatic colorectal cancer: a randomized phase III non-inferiority trial (SAKK 41/06)

BACKGROUND Chemotherapy plus bevacizumab is a standard option for first-line treatment in metastatic colorectal cancer (mCRC) patients. We assessed whether no continuation is non-inferior to continuation of bevacizumab after completing first-line chemotherapy. PATIENTS AND METHODS In an open-label, phase III multicentre trial, patients with mCRC without disease progression after 4-6 months of standard first-line chemotherapy plus bevacizumab were randomly assigned to continuing bevacizumab at a standard dose or no treatment. CT scans were done every 6 weeks until disease progression. The primary end point was time to progression (TTP). A non-inferiority limit for hazard ratio (HR) of 0.727 was chosen to detect a difference in TTP of 6 weeks or less, with a one-sided significance level of 10% and a statistical power of 85%. RESULTS The intention-to-treat population comprised 262 patients: median follow-up was 36.7 months. The median TTP was 4.1 [95% confidence interval (CI) 3.1-5.4] months for bevacizumab continuation versus 2.9 (95% CI 2.8-3.8) months for no continuation; HR 0.74 (95% CI 0.58-0.96). Non-inferiority could not be demonstrated. The median overall survival was 25.4 months for bevacizumab continuation versus 23.8 months (HR 0.83; 95% CI 0.63-1.1; P = 0.2) for no continuation. Severe adverse events were uncommon in the bevacizumab continuation arm. Costs for bevacizumab continuation were estimated to be ∼30,000 USD per patient. CONCLUSIONS Non-inferiority could not be demonstrated for treatment holidays versus continuing bevacizumab monotheray, after 4-6 months of standard first-line chemotherapy plus bevacizumab. Based on no impact on overall survival and increased treatment costs, bevacizumab as a single agent is of no meaningful therapeutic value. More efficient treatment approaches are needed to maintain control of stabilized disease following induction therapy. CLINICAL TRIAL REGISTRATION ClinicalTrials.gov, number NCT00544700.