A STUDY OF PENALTY FORMULATIONS USED IN THE NUMERICAL APPROXIMATION OF A RADIALLY SYMMETRIC ELASTICITY PROBLEM

We consider a class of two-dimensional problems in classical linear elasticity for which material overlapping occurs in the absence of singularities. Of course, material overlapping is not physically realistic, and one possible way to prevent it uses a constrained minimization theory. In this theory, a minimization problem consists of minimizing the total potential energy of a linear elastic body subject to the constraint that the deformation field must be locally invertible. Here, we use an interior and an exterior penalty formulation of the minimization problem together with both a standard finite element method and classical nonlinear programming techniques to compute the minimizers. We compare both formulations by solving a plane problem numerically in the context of the constrained minimization theory. The problem has a closed-form solution, which is used to validate the numerical results. This solution is regular everywhere, including the boundary. In particular, we show numerical results which indicate that, for a fixed finite element mesh, the sequences of numerical solutions obtained with both the interior and the exterior penalty formulations converge to the same limit function as the penalization is enforced. This limit function yields an approximate deformation field to the plane problem that is locally invertible at all points in the domain. As the mesh is refined, this field converges to the exact solution of the plane problem.

The random barrier-crossing problem

This paper addresses the time-variant reliability analysis of structures with random resistance or random system parameters. It deals with the problem of a random load process crossing a random barrier level. The implications of approximating the arrival rate of the first overload by an ensemble-crossing rate are studied. The error involved in this so-called ""ensemble-crossing rate"" approximation is described in terms of load process and barrier distribution parameters, and in terms of the number of load cycles. Existing results are reviewed, and significant improvements involving load process bandwidth, mean-crossing frequency and time are presented. The paper shows that the ensemble-crossing rate approximation can be accurate enough for problems where load process variance is large in comparison to barrier variance, but especially when the number of load cycles is small. This includes important practical applications like random vibration due to impact loadings and earthquake loading. Two application examples are presented, one involving earthquake loading and one involving a frame structure subject to wind and snow loadings. (C) 2007 Elsevier Ltd. All rights reserved.

Galerkin finite element method for incompressible thermofluid flows framed within the bond graph theory

In this paper a bond graph methodology is used to model incompressible fluid flows with viscous and thermal effects. The distinctive characteristic of these flows is the role of pressure, which does not behave as a state variable but as a function that must act in such a way that the resulting velocity field has divergence zero. Velocity and entropy per unit volume are used as independent variables for a single-phase, single-component flow. Time-dependent nodal values and interpolation functions are introduced to represent the flow field, from which nodal vectors of velocity and entropy are defined as state variables. The system for momentum and continuity equations is coincident with the one obtained by using the Galerkin method for the weak formulation of the problem in finite elements. The integral incompressibility constraint is derived based on the integral conservation of mechanical energy. The weak formulation for thermal energy equation is modeled with true bond graph elements in terms of nodal vectors of temperature and entropy rates, resulting a Petrov-Galerkin method. The resulting bond graph shows the coupling between mechanical and thermal energy domains through the viscous dissipation term. All kind of boundary conditions are handled consistently and can be represented as generalized effort or flow sources. A procedure for causality assignment is derived for the resulting graph, satisfying the Second principle of Thermodynamics. (C) 2007 Elsevier B.V. All rights reserved.

Assessing the evolution of sustainability reporting in the mining sector

Since the 1990s several large companies have been publishing nonfinancial performance reports. Focusing initially on the physical environment, these reports evolved to consider social relations, as well as data on the firm`s economic performance. A few mining companies pioneered this trend, and in the last years some of them incorporated the three dimensions of sustainable development, publishing so-called sustainability reports. This article reviews 31 reports published between 2001 and 2006 by four major mining companies. A set of 62 assessment items organized in six categories (namely context and commitment, management, environmental, social and economic performance, and accessibility and assurance) were selected to guide the review. The items were derived from international literature and recommended best practices, including the Global Reporting Initiative G3 framework. A content analysis was performed using the report as a sampling unit, and using phrases, graphics, or tables containing certain information as data collection units. A basic rating scale (0 or 1) was used for noting the presence or absence of information and a final percentage score was obtained for each report. Results show that there is a clear evolution in report`s comprehensiveness and depth. Categories ""accessibility and assurance"" and ""economic performance"" featured the lowest scores and do not present a clear evolution trend in the period, whereas categories ""context and commitment"" and ""social performance"" presented the best results and regular improvement; the category ""environmental performance,"" despite it not reaching the biggest scores, also featured constant evolution. Description of data measurement techniques, besides more comprehensive third-party verification are the items most in need of improvement.

Modeling the kinetics of the coalescence of water droplets in crude oil emulsions subject to an electric field, with the cellular automata technique

In this study, the concept of cellular automata is applied in an innovative way to simulate the separation of phases in a water/oil emulsion. The velocity of the water droplets is calculated by the balance of forces acting on a pair of droplets in a group, and cellular automata is used to simulate the whole group of droplets. Thus, it is possible to solve the problem stochastically and to show the sequence of collisions of droplets and coalescence phenomena. This methodology enables the calculation of the amount of water that can be separated from the emulsion under different operating conditions, thus enabling the process to be optimized. Comparisons between the results obtained from the developed model and the operational performance of an actual desalting unit are carried out. The accuracy observed shows that the developed model is a good representation of the actual process. (C) 2010 Published by Elsevier Ltd.

Minimizing total tardiness in a stochastic single machine scheduling problem using approximate dynamic programming

This paper addresses the non-preemptive single machine scheduling problem to minimize total tardiness. We are interested in the online version of this problem, where orders arrive at the system at random times. Jobs have to be scheduled without knowledge of what jobs will come afterwards. The processing times and the due dates become known when the order is placed. The order release date occurs only at the beginning of periodic intervals. A customized approximate dynamic programming method is introduced for this problem. The authors also present numerical experiments that assess the reliability of the new approach and show that it performs better than a myopic policy.

Scatter search for a real-life heterogeneous fleet vehicle routing problem with time windows and split deliveries in Brazil

In this paper, we consider a real-life heterogeneous fleet vehicle routing problem with time windows and split deliveries that occurs in a major Brazilian retail group. A single depot attends 519 stores of the group distributed in 11 Brazilian states. To find good solutions to this problem, we propose heuristics as initial solutions and a scatter search (SS) approach. Next, the produced solutions are compared with the routes actually covered by the company. Our results show that the total distribution cost can be reduced significantly when such methods are used. Experimental testing with benchmark instances is used to assess the merit of our proposed procedure. (C) 2008 Published by Elsevier B.V.

A Separation Principle for the Continuous-Time LQ-Problem With Markovian Jump Parameters

In this paper, we devise a separation principle for the finite horizon quadratic optimal control problem of continuous-time Markovian jump linear systems driven by a Wiener process and with partial observations. We assume that the output variable and the jump parameters are available to the controller. It is desired to design a dynamic Markovian jump controller such that the closed loop system minimizes the quadratic functional cost of the system over a finite horizon period of time. As in the case with no jumps, we show that an optimal controller can be obtained from two coupled Riccati differential equations, one associated to the optimal control problem when the state variable is available, and the other one associated to the optimal filtering problem. This is a separation principle for the finite horizon quadratic optimal control problem for continuous-time Markovian jump linear systems. For the case in which the matrices are all time-invariant we analyze the asymptotic behavior of the solution of the derived interconnected Riccati differential equations to the solution of the associated set of coupled algebraic Riccati equations as well as the mean square stabilizing property of this limiting solution. When there is only one mode of operation our results coincide with the traditional ones for the LQG control of continuous-time linear systems.

On the Filtering Problem for Continuous-Time Markov Jump Linear Systems with no Observation of the Markov Chain

We consider in this paper the optimal stationary dynamic linear filtering problem for continuous-time linear systems subject to Markovian jumps in the parameters (LSMJP) and additive noise (Wiener process). It is assumed that only an output of the system is available and therefore the values of the jump parameter are not accessible. It is a well known fact that in this setting the optimal nonlinear filter is infinite dimensional, which makes the linear filtering a natural numerically, treatable choice. The goal is to design a dynamic linear filter such that the closed loop system is mean square stable and minimizes the stationary expected value of the mean square estimation error. It is shown that an explicit analytical solution to this optimal filtering problem is obtained from the stationary solution associated to a certain Riccati equation. It is also shown that the problem can be formulated using a linear matrix inequalities (LMI) approach, which can be extended to consider convex polytopic uncertainties on the parameters of the possible modes of operation of the system and on the transition rate matrix of the Markov process. As far as the authors are aware of this is the first time that this stationary filtering problem (exact and robust versions) for LSMJP with no knowledge of the Markov jump parameters is considered in the literature. Finally, we illustrate the results with an example.

New simple and efficient heuristics for the uncapacitated single allocation hub location problem

Hub-and-spoke networks are widely studied in the area of location theory. They arise in several contexts, including passenger airlines, postal and parcel delivery, and computer and telecommunication networks. Hub location problems usually involve three simultaneous decisions to be made: the optimal number of hub nodes, their locations and the allocation of the non-hub nodes to the hubs. In the uncapacitated single allocation hub location problem (USAHLP) hub nodes have no capacity constraints and non-hub nodes must be assigned to only one hub. In this paper, we propose three variants of a simple and efficient multi-start tabu search heuristic as well as a two-stage integrated tabu search heuristic to solve this problem. With multi-start heuristics, several different initial solutions are constructed and then improved by tabu search, while in the two-stage integrated heuristic tabu search is applied to improve both the locational and allocational part of the problem. Computational experiments using typical benchmark problems (Civil Aeronautics Board (CAB) and Australian Post (AP) data sets) as well as new and modified instances show that our approaches consistently return the optimal or best-known results in very short CPU times, thus allowing the possibility of efficiently solving larger instances of the USAHLP than those found in the literature. We also report the integer optimal solutions for all 80 CAB data set instances and the 12 AP instances up to 100 nodes, as well as for the corresponding new generated AP instances with reduced fixed costs. Published by Elsevier Ltd.

The single machine earliness and tardiness scheduling problem: lower bounds and a branch-and-bound algorithm

This paper addresses the single machine scheduling problem with a common due date aiming to minimize earliness and tardiness penalties. Due to its complexity, most of the previous studies in the literature deal with this problem using heuristics and metaheuristics approaches. With the intention of contributing to the study of this problem, a branch-and-bound algorithm is proposed. Lower bounds and pruning rules that exploit properties of the problem are introduced. The proposed approach is examined through a computational comparative study with 280 problems involving different due date scenarios. In addition, the values of optimal solutions for small problems from a known benchmark are provided.

Landscape unit discrimination for pedological surveys by orbital spectral response

Landscape unit discrimination for pedological surveys by orbital spectral response. The objective of tins study was compare two soil survey methods. The first was performed by methods traditionally used to distinguish landscape units and soil class discrimination. The second was based on soil class distinction through orbital spectral response. In order to establish soil characteristics and their classification, soil samples were collected at two depths in a grid system, with a distance of 500 meters between points. With these samples, physical and chemical analyses were carried out. In the sampling points, the apparent reflectance of the soil, front the orbital image, was determined and, through cluster analysis landscape units were established. In order to evaluate the resemblance reliability between the landscape units established in each method, the Kappa index was used, the value set for the confusion matrix was 0.43, indicating high quality in the comparison, showing that the non-conventional method was as close as the one carried out by photointerpretation.

Natural regeneration of Atlantic Forest species in the understory of Eucalyptus saligna Smith. in a former forest production unit at the Parque das Neblinas, Bertioga, SP

This study aimed at characterizing the potential for natural regeneration of native vegetation in the under-story of an earlier Eucalyptus saligna Smith production stand. The study was carried out at the Parque das Neblinas, Bertioga municipality, SP, in a 45 ha third rotation stand; which had been abandoned 15 years ago for natural regeneration to occur. The sampling was done in 24 plots of 20 x 40 m. The sampled area was of 19,200 m(2), with inventory made of 100% of the eucalyptus trees. All regeneration trees with a height >= 1.30 m and DBH >= 5.0 cm were measured, as well as adult individuals with DBH >= 5.0 cm; surveyed in two size classes. 1,417 individuals of E. saligna were measured, with a density of 738,02 individuals/ha and a basal area of 22.69 m(2)/ha. Among 2,763 natural regeneration individuals, 111 species belonged to 66 genera and 34 botanical families. The species represented 43.7% of the tree richness of neighboring native forest fragments. The total estimated density and the basal area were respectively 1,052.6 individuals/ha and 6.4 m(2)/ha of autochthonous trees with DBH >= 5.0 cm (Class 1); while for regeneration there were 3,864.58 individuals/ha, and 2.76 m(2)/ha of individuals with a height >= 1.30 m and DBH <5.0 cm (Class 2). Shannon diversity (H`) was 2.83 and 3.68, respectively, for Classes 1 and 2, and the corrected species richness for a 1000-individual sample (R(1000)) were 75.6 and 87.29 (Fisher`s a index) for the same classes. The majority of the species (34.84%) was typical from the understory of wet tropical forest and had zoochoric fruit dispersal (67.57%). The results indicate that, under these conditions, a eucalyptus forest is able to provide adequate regeneration niches for native vegetation, and may represent a sink habitat for local populations.

Nosocomial outbreak by imipenem-resistant metallo-beta-lactamase-producing Pseudomonas aeruginosa in an adult intensive care unit in a Brazilian teaching hospital

Objective: To describe an outbreak of imipenem-resistant metallo-beta-lactamase-producing Pseudomonas aeruginosa, enzyme type bla, by horizontal transmission in patients admitted to a mixed adult ICU. Methods: A case-control study was carried out, including 47 patients (cases) and 122 patients (control) admitted to the mixed ICU of a university hospital in Minas Gerais. Brazil from November 2003 to July 2005. The infection site, risk factors, mortality, antibiotic susceptibility, metallo-beta-lactamase (MBL) production, enzyme type, and clonal diversity were analyzed, Results: A temporal/spatial relationship was detected in most patients (94%), overall mortality was 55.3%, and pneumonia was the predominant infection (85%). The majority of isolates (95%) were resistant to imipenem and other antibiotics, except for polymyxin, and showed MBL production (76.7%). Only bla SPM-1 (33%) was identified in the 15 specimens analyzed. In addition, 4 clones were identified, with a predominance of clone A (61.5%) and B (23.1%). On multivariate analysis, advanced age, mechanical ventilation, tracheostomy, and previous imipenem use were significant risk factors for imipenem-resistant P. aeruginosa infection. Conclusions: Clonal dissemination of MBL-producing P. aeruginosa strains with a spatial/temporal relationship disclosed problems in the practice of hospital infection control, low adherence to hand hygiene, and empirical antibiotic use. (C) 2008 Elsevier Espana, S.L. All rights reserved.

Aboriginal and Torres Strait Islander Studies Unit Annual Report

Welcome to the 2002 Aboriginal and Torres Strait Islander Studies Unit Annual Report. This report is a brief summary of Unit activities during the 2002 calendar year. The Unit provides personal and academic support for Aboriginal and Torres Strait Islander students and specifically aims to increase the recruitment, retention, academic performance and graduation rates of Indigenous students. The Unit also administers schemes to help Indigenous students gain access to, and receive support in, tertiary studies such as the Alternative Entry scheme and the federally funded Aboriginal Tutorial Assistance Scheme (ATAS). The Unit is also the focus for teaching and research in Aboriginal and Torres Strait Islander Studies at the University of Queensland.