Optimization heuristic solutions, how good can they be? : With empirical applications in location problems

Combinatorial optimization problems, are one of the most important types of problems in operational research. Heuristic and metaheuristics algorithms are widely applied to find a good solution. However, a common problem is that these algorithms do not guarantee that the solution will coincide with the optimum and, hence, many solutions to real world OR-problems are afflicted with an uncertainty about the quality of the solution. The main aim of this thesis is to investigate the usability of statistical bounds to evaluate the quality of heuristic solutions applied to large combinatorial problems. The contributions of this thesis are both methodological and empirical. From a methodological point of view, the usefulness of statistical bounds on p-median problems is thoroughly investigated. The statistical bounds have good performance in providing informative quality assessment under appropriate parameter settings. Also, they outperform the commonly used Lagrangian bounds. It is demonstrated that the statistical bounds are shown to be comparable with the deterministic bounds in quadratic assignment problems. As to empirical research, environment pollution has become a worldwide problem, and transportation can cause a great amount of pollution. A new method for calculating and comparing the CO2-emissions of online and brick-and-mortar retailing is proposed. It leads to the conclusion that online retailing has significantly lesser CO2-emissions. Another problem is that the Swedish regional division is under revision and the border effect to public service accessibility is concerned of both residents and politicians. After analysis, it is shown that borders hinder the optimal location of public services and consequently the highest achievable economic and social utility may not be attained.

Hamiltonian constraint analysis of vector field theories with spontaneous Lorentz symmetry breaking

Recent investigations of various quantum-gravity theories have revealed a variety of possible mechanisms that lead to Lorentz violation. One of the more elegant of these mechanisms is known as Spontaneous Lorentz Symmetry Breaking (SLSB), where a vector or tensor field acquires a nonzero vacuum expectation value. As a consequence of this symmetry breaking, massless Nambu-Goldstone modes appear with properties similar to the photon in Electromagnetism. This thesis considers the most general class of vector field theories that exhibit spontaneous Lorentz violation-known as bumblebee models-and examines their candidacy as potential alternative explanations of E&M, offering the possibility that Einstein-Maxwell theory could emerge as a result of SLSB rather than of local U(1) gauge invariance. With this aim we employ Dirac's Hamiltonian Constraint Analysis procedure to examine the constraint structures and degrees of freedom inherent in three candidate bumblebee models, each with a different potential function, and compare these results to those of Electromagnetism. We find that none of these models share similar constraint structures to that of E&M, and that the number of degrees of freedom for each model exceeds that of Electromagnetism by at least two, pointing to the potential existence of massive modes or propagating ghost modes in the bumblebee theories.

Otimização Topológica com Termo-elasticidade Utilizando Elementos Tri3

This work proposes a formulation for optimization of 2D-structure layouts submitted to mechanic and thermal shipments and applied an h-adaptive filter process which conduced to computational low spend and high definition structural layouts. The main goal of the formulation is to minimize the structure mass submitted to an effective state of stress of von Mises, with stability and lateral restriction variants. A criterion of global measurement was used for intents a parametric condition of stress fields. To avoid singularity problems was considerate a release on the stress restriction. On the optimization was used a material approach where the homogenized constructive equation was function of the material relative density. The intermediary density effective properties were represented for a SIMP-type artificial model. The problem was simplified by use of the method of finite elements of Galerkin using triangles with linear Lagrangian basis. On the solution of the optimization problem, was applied the augmented Lagrangian Method, that consists on minimum problem sequence solution with box-type restrictions, resolved by a 2nd orderprojection method which uses the method of the quasi-Newton without memory, during the problem process solution. This process reduces computational expends showing be more effective and solid. The results materialize more refined layouts with accurate topologic and shape of structure definitions. On the other hand formulation of mass minimization with global stress criterion provides to modeling ready structural layouts, with violation of the criterion of homogeneous distributed stress

Otimização de forma aplicando B-splines sob critério integral de tensões

This work proposes a computational methodology to solve problems of optimization in structural design. The application develops, implements and integrates methods for structural analysis, geometric modeling, design sensitivity analysis and optimization. So, the optimum design problem is particularized for plane stress case, with the objective to minimize the structural mass subject to a stress criterion. Notice that, these constraints must be evaluated at a series of discrete points, whose distribution should be dense enough in order to minimize the chance of any significant constraint violation between specified points. Therefore, the local stress constraints are transformed into a global stress measure reducing the computational cost in deriving the optimal shape design. The problem is approximated by Finite Element Method using Lagrangian triangular elements with six nodes, and use a automatic mesh generation with a mesh quality criterion of geometric element. The geometric modeling, i.e., the contour is defined by parametric curves of type B-splines, these curves hold suitable characteristics to implement the Shape Optimization Method, that uses the key points like design variables to determine the solution of minimum problem. A reliable tool for design sensitivity analysis is a prerequisite for performing interactive structural design, synthesis and optimization. General expressions for design sensitivity analysis are derived with respect to key points of B-splines. The method of design sensitivity analysis used is the adjoin approach and the analytical method. The formulation of the optimization problem applies the Augmented Lagrangian Method, which convert an optimization problem constrained problem in an unconstrained. The solution of the Augmented Lagrangian function is achieved by determining the analysis of sensitivity. Therefore, the optimization problem reduces to the solution of a sequence of problems with lateral limits constraints, which is solved by the Memoryless Quasi-Newton Method It is demonstrated by several examples that this new approach of analytical design sensitivity analysis of integrated shape design optimization with a global stress criterion purpose is computationally efficient

Otimização topológica de estruturas termoelásticas tridimensionais

This work presents an optimization technique based on structural topology optimization methods, TOM, designed to solve problems of thermoelasticity 3D. The presented approach is based on the adjoint method of sensitivity analysis unified design and is intended to loosely coupled thermomechanical problems. The technique makes use of analytical expressions of sensitivities, enabling a reduction in the computational cost through the use of a coupled field adjoint equation, defined in terms the of temperature and displacement fields. The TOM used is based on the material aproach. Thus, to make the domain is composed of a continuous distribution of material, enabling the use of classical models in nonlinear programming optimization problem, the microstructure is considered as a porous medium and its constitutive equation is a function only of the homogenized relative density of the material. In this approach, the actual properties of materials with intermediate densities are penalized based on an artificial microstructure model based on the SIMP (Solid Isotropic Material with Penalty). To circumvent problems chessboard and reduce dependence on layout in relation to the final optimal initial mesh, caused by problems of numerical instability, restrictions on components of the gradient of relative densities were applied. The optimization problem is solved by applying the augmented Lagrangian method, the solution being obtained by applying the finite element method of Galerkin, the process of approximation using the finite element Tetra4. This element has the ability to interpolate both the relative density and the displacement components and temperature. As for the definition of the problem, the heat load is assumed in steady state, i.e., the effects of conduction and convection of heat does not vary with time. The mechanical load is assumed static and distributed

Simulação numérica da interação entre uma nuvem incidente de vórtices e uma aerofólio através do método dos vórtices

The study of aerodynamic loading variations has many engineering applications, including helicopter rotor blades, wind turbines and turbo machinery. This work uses a Vortex Method to make a lagrangian description of the a twodimensional airfoil/ incident wake vortex interaction. The flow is incompressible, newtonian, homogeneus and the Reynolds Number is 5x105 .The airfoil is a NACA 0018 placed a angle of attack of the 0° and 5°simulates with the Painel Method with a constant density vorticity panels and a generation poit is near the painel. The protector layer is created does not permit vortex inside the body. The vortex Lamb convection is realized with the Euler Method (first order) and Adans-Bashforth (second order). The Random Walk Method is used to simulate the diffusion. The circular wake has 366 vortex all over positive or negative vorticity located at different heights with respect to the airfoil chord. The Lift was calculated based in the algorithm created by Ricci (2002). This simulation uses a ready algorithm vatidated with single body does not have a incident wake. The results are compared with a experimental work The comparasion concludes that the experimental results has a good agrement with this papper

Um estudo na teoria do movimento brownianno com viscosidade variável

In this work we investigate the stochastic behavior of a large class of systems with variable damping which are described by a time-dependent Lagrangian. Our stochastic approach is based on the Langevin treatment describing the motion of a classical Brownian particle of mass m. Two situations of physical interest are considered. In the first one, we discuss in detail an application of the standard Langevin treatment (white noise) for the variable damping system. In the second one, a more general viewpoint is adopted by assuming a given expression to the so-called collored noise. For both cases, the basic diffententiaql equations are analytically solved and al the quantities physically relevant are explicitly determined. The results depend on an arbitrary q parameter measuring how the behavior of the system departs from the standard brownian particle with constant viscosity. Several types of sthocastic behavior (superdiffusive and subdiffusive) are obteinded when the free pamameter varies continuosly. However, all the results of the conventional Langevin approach with constant damping are recovered in the limit q = 1

Um método primal-dual aplicado na resolução do problema de fluxo de potência ótimo

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

Fluxo de potência ótimo reativo via método da função lagrangiana barreira modificada

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

Loss minimization by the predictor-corrector modified barrier approach

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Multiflavor soldering

In two dimensions the simple addition of two chiral bosons of opposite chiralities does not lead to a full massless scalar field. Similarly, in three dimensions the addition of two Maxwell-Chern-Simons fields of opposite helicities +/- 1 will not produce a parity invariant Maxwell-Proca theory. An interference term between the opposite chiralities (helicities) states is required in order to obtain the expected result. The so-called soldering procedure provides the missing interference Lagrangian in both 2D and 3D cases. In two dimensions such interference term allows to fuse two chiral fermionic determinants into, a non-chiral one. In a recent work we have generalized this procedure by allowing the appearance of an extra parameter which takes two possible values and leads to two different soldered Lagrangians. Here we apply this generalized soldering in a bosonic theory which has appeared in a partial bosonization of the 3D gauged Thirring model with N flavors. The multiplicity of flavors allow new types of solderings and help us to understand the connection between different perturbative approaches to bosonization in 3D. In particular, we obtain an interference term which takes us from a multiflavor Niaxwell-Chern-Simons theory to a pair of self-dual and anti-self-dual theories when we combine together both fermionic determinants of +1/2 and -1/2 helicity fermions. An important role is played by a set of pure non-interacting Chern-Simons fields which amount to a normalization factor in the fermionic determinants and act like spectators in the original theory but play an active role in the soldering procedure. Our results suggest that the generalized soldering could be used to provide dual theories in both 2D and 3D cases. (c) 2007 Elsevier B.V. All rights reserved.

Computations of the flow past a still sphere at moderate reynolds numbers using an immersed boundary method

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Termodinâmica de um gás de fótons no contexto de eletrodinâmicas não-lineares

The objective of this dissertation is the development of a general formalism to analyze the thermodynamical properties of a photon gas under the context of nonlinear electrodynamics (NLED). To this end it is obtained, through the systematic analysis of Maxwell s electromagnetism (EM) properties, the general dependence of the Lagrangian that describes this kind of theories. From this Lagrangian and in the background of classical field theory, we derive the general dispersion relation that photons must obey in terms of a background field and the NLED properties. It is important to note that, in order to achieve this result, an aproximation has been made in order to allow the separation of the total electromagnetic field into a strong background electromagnetic field and a perturbation. Once the dispersion relation is in hand, the usual Bose-Einstein statistical procedure is followed through which the thermodynamical properties, energy density and pressure relations are obtained. An important result of this work is the fact that equation of state remains identical to the one obtained under EM. Then, two examples are made where the thermodynamic properties are explicitly derived in the context of two NLED, Born-Infelds and a quadratic approximation. The choice of the first one is due to the vast appearance in literature and, the second one, because it is a first order approximation of a large class of NLED. Ultimately, both are chosen because of their simplicity. Finally, the results are compared to EM and interpreted, suggesting possible tests to verify the internal consistency of NLED and motivating further developement into the formalism s quantum case

A marker-and-cell approach to free surface 2-D multiphase flows

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Co-orbital satellites of Saturn: congenital formation

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