O método de pontos interiores no planejamento da radioterapia

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

Planejamento de reativos em sistemas de energia elétrica através de um algoritmo de Branch-and-Bound não linear

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

A matheuristic for the selection of beam directions and dose distribution in radiotherapy planning

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

Optimal day-ahead scheduling of integrated urban energy systems

An optimal day-ahead scheduling method (ODSM) for the integrated urban energy system (IUES) is introduced, which considers the reconfigurable capability of an electric distribution network. The hourly topology of a distribution network, a natural gas network, the energy centers including the combined heat and power (CHP) units, different energy conversion devices and demand responsive loads (DRLs), are optimized to minimize the day-ahead operation cost of the IUES. The hourly reconfigurable capability of the electric distribution network utilizing remotely controlled switches (RCSs) is explored and discussed. The operational constraints from the unbalanced three-phase electric distribution network, the natural gas network, and the energy centers are considered. The interactions between the electric distribution network and the natural gas network take place through conversion of energy among different energy vectors in the energy centers. An energy conversion analysis model for the energy center was developed based on the energy hub model. A hybrid optimization method based on genetic algorithm (GA) and a nonlinear interior point method (IPM) is utilized to solve the ODSM model. Numerical studies demonstrate that the proposed ODSM is able to provide the IUES with an effective and economical day-ahead scheduling scheme and reduce the operational cost of the IUES.

An Efficient Precoder Design for Multiuser MIMO Cognitive Radio Networks with Interference Constraints

We consider a linear precoder design for an underlay cognitive radio multiple-input multiple-output broadcast channel, where the secondary system consisting of a secondary base-station (BS) and a group of secondary users (SUs) is allowed to share the same spectrum with the primary system. All the transceivers are equipped with multiple antennas, each of which has its own maximum power constraint. Assuming zero-forcing method to eliminate the multiuser interference, we study the sum rate maximization problem for the secondary system subject to both per-antenna power constraints at the secondary BS and the interference power constraints at the primary users. The problem of interest differs from the ones studied previously that often assumed a sum power constraint and/or single antenna employed at either both the primary and secondary receivers or the primary receivers. To develop an efficient numerical algorithm, we first invoke the rank relaxation method to transform the considered problem into a convex-concave problem based on a downlink-uplink result. We then propose a barrier interior-point method to solve the resulting saddle point problem. In particular, in each iteration of the proposed method we find the Newton step by solving a system of discrete-time Sylvester equations, which help reduce the complexity significantly, compared to the conventional method. Simulation results are provided to demonstrate fast convergence and effectiveness of the proposed algorithm. 

Handgrip strength cut-off for cardiometabolic risk index among Colombian children and adolescents : the FUPRECOL study

Abstract Objective: Evidence shows an association between muscular strength (MS) and health among youth, however low muscular strength cut-points for the detection of high metabolic risk in Latin-American populations are scarce. The aim of this study was two-fold: to explore potential age- and sex-specific thresholds of MS, for optimal cardiometabolic risk categorization among Colombian children and adolescents; and to investigate if cardiometabolic risk differed by MS group by applying the receiver operating characteristic curve (ROC) cut point. Methods: This is a secondary analysis of a cross-sectional study (the FUPRECOL study), published elsewhere. The FUPRECOL study assessments were conducted during the 2014– 2015 school year. MS was estimated by a handle dynamometer on 1,950 children and adolescents from Colombia, using the MS relative to weight (handgrip strength/body mass). A metabolic risk score was computed from the following components: waist circumference, triglycerides, HDL-c, glucose, systolic and diastolic blood pressure. ROC analysis showed a significant discriminatory accuracy of MS in identifying the low/high metabolic risk in children and adolescents and both gender. Results: In children, handgrip strength/body mass level for a low metabolic risk were 0.359 and 0.376 in girls and boys, respectively. In adolescents, these points were 0.440 and 0.447 in girls and boys, respectively. Conclusion: In conclusion, the results suggest a hypothetical MS level relative to weight for having a low metabolic risk, which could be used to identify youths at risk.

Auditoria à qualidade do ar interior de um edifício de serviços - RSECE QAI

O objectivo principal deste trabalho é a realização de uma auditoria, à qualidade do ar interior (QAI), a um edifício de serviços – COCIGA, SA, tendo como base o Regulamento dos Sistemas Energéticos de Climatização dos Edifícios (RSECE). A auditoria QAI implica a medição de vários parâmetros físicos, químicos, microbiológicos e também a inspecção aos componentes do sistema de climatização com a finalidade de averiguar o seu estado de limpeza e manutenção. Assim, foram seleccionados 3 espaços, para a realização de amostragens designados por Comercial - Produtos, AVAC e Mezaninne das oficinas, nos quais foi efectuada a medição de diversos parâmetros, de acordo com as imposições do RSECE, utilizando medidores portáteis ou recorrendo a métodos analíticos. Relativamente aos parâmetros físicos, registaram-se valores de temperatura, para os três espaços estudados, entre os 21 e os 24 ºC e valores médios de humidade relativa de cerca de 50 %. Outro parâmetro medido, e de grande importância para garantir o conforto dos ocupantes, foi a velocidade do ar nos postos de trabalho. De acordo com o RSECE este valor não deve ser superior a 0,2 m/s, o que se verificou em todos os pontos medidos. O último parâmetro físico medido foi a concentração de partículas (PM10) tendo-se obtido valores de cerca de 23 μg/m3ar, valor bastante inferior ao máximo permitido pelo RSECE (150 μg/m3ar). Também no que diz respeito aos parâmetros químicos, ou seja, CO2, CO, formaldeído e ozono, não se verificaram valores superiores aos regulamentares. No caso do CO2, o valor máximo encontrado, nestes três espaços, foi de 745 ppm na Mezaninne das Oficinas e para o CO, na zona AVAC com uma concentração de 0,73 ppm. A medição do formaldeído registou valores perto dos 45 μg/m3ar e o ozono apenas foi detectado, em concentração muito reduzida, na zona Comercial – Produtos. Por fim, as concentrações de bactérias e fungos, de acordo com o RSECE, não devem ultrapassar as 500 UFC/m3ar (parâmetros microbiológicos). Em qualquer dos espaços, os valores medidos foram inferiores ao máximo legal, não ultrapassando as 50 UFC/m3ar. Da avaliação do projecto AVAC, e através da medição dos caudais de insuflação/ extracção em cada zona, concluiu-se que os seus valores não estão de acordo com os valores do projecto inicial que poderá ser imputada a uma insuficiência no funcionamento do sistema detectada na altura das medições. No que diz respeito ao estado de limpeza do sistema AVAC, apenas foi possível inspeccionar as unidades de tratamento de ar, tendo-se constatado que se encontram em boas condições. Ou seja, do ponto de vista do RSECE, e referindo-nos apenas à vertente da Qualidade do Ar Interior, o edifício em causa, cumpre todos os limites impostos para as concentrações de poluentes mas, apresenta algumas deficiências no que respeita aos caudais de ar novo insuflados em cada espaço.

How do current term structure model behave beyond the last liquid point?: A comparison of the DNS and Smith-Wilson methods

A Work Project, presented as part of the requirements for the Award of a Master’s Double Degree in Finance from Maastricht University and NOVA – School of Business and Economics

Stabilized Petrov-Galerkin methods for the convection-diffusion-reaction and the Helmholtz equations

We present two new stabilized high-resolution numerical methods for the convection–diffusion–reaction (CDR) and the Helmholtz equations respectively. The work embarks upon a priori analysis of some consistency recovery procedures for some stabilization methods belonging to the Petrov–Galerkin framework. It was found that the use of some standard practices (e.g. M-Matrices theory) for the design of essentially non-oscillatory numerical methods is not feasible when consistency recovery methods are employed. Hence, with respect to convective stabilization, such recovery methods are not preferred. Next, we present the design of a high-resolution Petrov–Galerkin (HRPG) method for the 1D CDR problem. The problem is studied from a fresh point of view, including practical implications on the formulation of the maximum principle, M-Matrices theory, monotonicity and total variation diminishing (TVD) finite volume schemes. The current method is next in line to earlier methods that may be viewed as an upwinding plus a discontinuity-capturing operator. Finally, some remarks are made on the extension of the HRPG method to multidimensions. Next, we present a new numerical scheme for the Helmholtz equation resulting in quasi-exact solutions. The focus is on the approximation of the solution to the Helmholtz equation in the interior of the domain using compact stencils. Piecewise linear/bilinear polynomial interpolation are considered on a structured mesh/grid. The only a priori requirement is to provide a mesh/grid resolution of at least eight elements per wavelength. No stabilization parameters are involved in the definition of the scheme. The scheme consists of taking the average of the equation stencils obtained by the standard Galerkin finite element method and the classical finite difference method. Dispersion analysis in 1D and 2D illustrate the quasi-exact properties of this scheme. Finally, some remarks are made on the extension of the scheme to unstructured meshes by designing a method within the Petrov–Galerkin framework.

GPS to LRM: Integration of Spatial Point Features with Linear Referencing Methods, 2001

Global positioning systems (GPS) offer a cost-effective and efficient method to input and update transportation data. The spatial location of objects provided by GPS is easily integrated into geographic information systems (GIS). The storage, manipulation, and analysis of spatial data are also relatively simple in a GIS. However, many data storage and reporting methods at transportation agencies rely on linear referencing methods (LRMs); consequently, GPS data must be able to link with linear referencing. Unfortunately, the two systems are fundamentally incompatible in the way data are collected, integrated, and manipulated. In order for the spatial data collected using GPS to be integrated into a linear referencing system or shared among LRMs, a number of issues need to be addressed. This report documents and evaluates several of those issues and offers recommendations. In order to evaluate the issues associated with integrating GPS data with a LRM, a pilot study was created. To perform the pilot study, point features, a linear datum, and a spatial representation of a LRM were created for six test roadway segments that were located within the boundaries of the pilot study conducted by the Iowa Department of Transportation linear referencing system project team. Various issues in integrating point features with a LRM or between LRMs are discussed and recommendations provided. The accuracy of the GPS is discussed, including issues such as point features mapping to the wrong segment. Another topic is the loss of spatial information that occurs when a three-dimensional or two-dimensional spatial point feature is converted to a one-dimensional representation on a LRM. Recommendations such as storing point features as spatial objects if necessary or preserving information such as coordinates and elevation are suggested. The lack of spatial accuracy characteristic of most cartography, on which LRM are often based, is another topic discussed. The associated issues include linear and horizontal offset error. The final topic discussed is some of the issues in transferring point feature data between LRMs.

Desenvolupament d'una tecnologia per a la millora de les recomanacions en la fertilització nitrogenada dels cereals d'hivern. (I) Cereals d'hivern a la Catalunya interior

Es va realitzar una sèrie d'assaigs d'adobat nitrogenat en diferents comarques de la Catalunya interior. En el conjunt d'aquests assaigs es varen comprovar tres mètodes diferents que es va considerar que eren prometedors per tal de millorar la fertilització nitrogenada. Els mètodes assajats eren el mètode del balanç de nitrogen, el del nitrogen mineral i el del contingut de nitrats al suc de la base de les tiges (CNSBT). Els sòls on es van realitzar els assaigs no presentaven cap limitació especial per al cultiu del blat i eren profunds, ben drenats, no salins i de textura mitjana; l'única excepció era un assaig sobre sòl moderadament profund. Per tant, i també pel que fa a la fertilitat química, els sòls s'han de considerar d'un potencial productiu mitjàalt. El mètode del balanç de nitrogen s'ha mostrat com a molt prometedor de cara a definir si cal la magnitud de l'adobat de cobertora per a les condicions estudiades. El mètode de nitrogen mineral també ha estat efectiu en aquest sentit, mentre que el del CNSBT s'ha revelat com a no aplicable en les condicions assajades, on en molts casos l'aigua és també factor limitant. Al llarg dels assaigs s'han identificat un seguit de factors que impedeixen ajustar la fertilitat nitrogenada. Entre aquests cal esmentar la mala estimació de la producció objectiu, la dificultat de predir el N disponible a partir dels adobs orgànics, dificultats de mostreig pel nitrogen nítric i l'efecte crític que té l'erràtica disponibilitat d'aigua que complica molt l'estratègia de fertilització nitrogenada a adoptar.

Numerical Methods for Stiff Two-Point Boundary Value Problems

We consider the two-point boundary value problem for stiff systems of ordinary differential equations. For systems that can be transformed to essentially diagonally dominant form with appropriate smoothness conditions, a priori estimates are obtained. Problems with turning points can be treated with this theory, and we discuss this in detail. We give robust difference approximations and present error estimates for these schemes. In particular we give a detailed description of how to transform a general system to essentially diagonally dominant form and then stretch the independent variable so that the system will satisfy the correct smoothness conditions. Numerical examples are presented for both linear and nonlinear problems.

Parameterized fast decoupled power flow methods for obtaining the maximum loading point of power systems. Part I. Mathematical modeling

The conventional Newton and fast decoupled power flow (FDPF) methods have been considered inadequate to obtain the maximum loading point of power systems due to ill-conditioning problems at and near this critical point. It is well known that the PV and Q-theta decoupling assumptions of the fast decoupled power flow formulation no longer hold in the vicinity of the critical point. Moreover, the Jacobian matrix of the Newton method becomes singular at this point. However, the maximum loading point can be efficiently computed through parameterization techniques of continuation methods. In this paper it is shown that by using either theta or V as a parameter, the new fast decoupled power flow versions (XB and BX) become adequate for the computation of the maximum loading point only with a few small modifications. The possible use of reactive power injection in a selected PV bus (Q(PV)) as continuation parameter (mu) for the computation of the maximum loading point is also shown. 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 in predictor step. 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 approach for the IEEE test systems (14, 30, 57 and 118 buses) are presented and discussed in the companion paper. The results show that the characteristics of the conventional method are enhanced and the region of convergence around the singular solution is enlarged. In addition, it is shown that parameters can be switched during the tracing process in order to efficiently determine all the PV curve points with few iterations. (C) 2003 Elsevier B.V. All rights reserved.

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.