Modelo para distribuição de recursos nos municípios brasileiros baseado na lei de diretrizes orçamentárias, análise multicritério e programação linear

The municipal management in any country of the globe requires planning and allocation of resources evenly. In Brazil, the Law of Budgetary Guidelines (LDO) guides municipal managers toward that balance. This research develops a model that seeks to find the balance of the allocation of public resources in Brazilian municipalities, considering the LDO as a parameter. For this using statistical techniques and multicriteria analysis as a first step in order to define allocation strategies, based on the technical aspects arising from the municipal manager. In a second step, presented in linear programming based optimization where the objective function is derived from the preference of the results of the manager and his staff. The statistical representation is presented to support multicriteria development in the definition of replacement rates through time series. The multicriteria analysis was structured by defining the criteria, alternatives and the application of UTASTAR methods to calculate replacement rates. After these initial settings, an application of linear programming was developed to find the optimal allocation of enforcement resources of the municipal budget. Data from the budget of a municipality in southwestern Paraná were studied in the application of the model and analysis of results.

Multivariate orthogonal polynomials and integrable systems

Multivariate orthogonal polynomials in D real dimensions are considered from the perspective of the Cholesky factorization of a moment matrix. The approach allows for the construction of corresponding multivariate orthogonal polynomials, associated second kind functions, Jacobi type matrices and associated three term relations and also Christoffel-Darboux formulae. The multivariate orthogonal polynomials, their second kind functions and the corresponding Christoffel-Darboux kernels are shown to be quasi-determinants as well as Schur complements of bordered truncations of the moment matrix; quasi-tau functions are introduced. It is proven that the second kind functions are multivariate Cauchy transforms of the multivariate orthogonal polynomials. Discrete and continuous deformations of the measure lead to Toda type integrable hierarchy, being the corresponding flows described through Lax and Zakharov-Shabat equations; bilinear equations are found. Varying size matrix nonlinear partial difference and differential equations of the 2D Toda lattice type are shown to be solved by matrix coefficients of the multivariate orthogonal polynomials. The discrete flows, which are shown to be connected with a Gauss-Borel factorization of the Jacobi type matrices and its quasi-determinants, lead to expressions for the multivariate orthogonal polynomials and their second kind functions in terms of shifted quasi-tau matrices, which generalize to the multidimensional realm, those that relate the Baker and adjoint Baker functions to ratios of Miwa shifted tau-functions in the 1D scenario. In this context, the multivariate extension of the elementary Darboux transformation is given in terms of quasi-determinants of matrices built up by the evaluation, at a poised set of nodes lying in an appropriate hyperplane in R^D, of the multivariate orthogonal polynomials. The multivariate Christoffel formula for the iteration of m elementary Darboux transformations is given as a quasi-determinant. It is shown, using congruences in the space of semi-infinite matrices, that the discrete and continuous flows are intimately connected and determine nonlinear partial difference-differential equations that involve only one site in the integrable lattice behaving as a Kadomstev-Petviashvili type system. Finally, a brief discussion of measures with a particular linear isometry invariance and some of its consequences for the corresponding multivariate polynomials is given. In particular, it is shown that the Toda times that preserve the invariance condition lay in a secant variety of the Veronese variety of the fixed point set of the linear isometry.

Recycling Krylov subspaces for efficient large-scale electrical impedance tomography

Electrical impedance tomography (EIT) captures images of internal features of a body. Electrodes are attached to the boundary of the body, low intensity alternating currents are applied, and the resulting electric potentials are measured. Then, based on the measurements, an estimation algorithm obtains the three-dimensional internal admittivity distribution that corresponds to the image. One of the main goals of medical EIT is to achieve high resolution and an accurate result at low computational cost. However, when the finite element method (FEM) is employed and the corresponding mesh is refined to increase resolution and accuracy, the computational cost increases substantially, especially in the estimation of absolute admittivity distributions. Therefore, we consider in this work a fast iterative solver for the forward problem, which was previously reported in the context of structural optimization. We propose several improvements to this solver to increase its performance in the EIT context. The solver is based on the recycling of approximate invariant subspaces, and it is applied to reduce the EIT computation time for a constant and high resolution finite element mesh. In addition, we consider a powerful preconditioner and provide a detailed pseudocode for the improved iterative solver. The numerical results show the effectiveness of our approach: the proposed algorithm is faster than the preconditioned conjugate gradient (CG) algorithm. The results also show that even on a standard PC without parallelization, a high mesh resolution (more than 150,000 degrees of freedom) can be used for image estimation at a relatively low computational cost. (C) 2010 Elsevier B.V. All rights reserved.

Electrochemical impedance spectroscopy investigation of the electrochemical behaviour of copper coated with artificial patina layers and submitted to wet and dry cycles

Due to rain events historical monuments exposed to the atmosphere are frequently submitted to wet and dry cycles. During drying periods wetness is maintained in some confined regions and the corrosion product layer, generally denominated patinas, builds up and gets thicker. The aim of this study is to use electrochemical impedance spectroscopy (EIS) to investigate the electrochemical behaviour of pure copper coated with two artificial patina layers and submitted either to continuous or to intermittent immersion tests, this latter aiming to simulate wet and dry cycles. The experiments were performed in 0.1 mol dm(-3) NaCl solution and in artificial rainwater containing the most significant pollutants of the city of Sao Paulo. The results of the continuous immersion tests in the NaCl solution have shown that the coated samples behave like a porous electrode with finite pore length. On the other hand, in the intermittent tests a porous electrode response with semi-infinite pore length can be developed. The results were interpreted based on the model of de Levie and a critical comparison with previous interpretations reported in the literature for similar systems is presented. (C) 2011 Elsevier Ltd. All rights reserved.

A two-dimensional analytical solution of groundwater responses to tidal loading in an estuary and ocean

Previous studies on tidal dynamics of coastal aquifers have focussed on the inland propagation of oceanic tides in the cross-shore direction, a configuration that is essentially one-dimensional. Aquifers at natural coasts can also be influenced by tidal waves in nearby estuaries, resulting in a more complex behaviour of head fluctuations in the aquifers. We present an analytical solution to the two-dimensional depth-averaged groundwater flow equation for a semi-infinite aquifer subject to oscillating head conditions at the boundaries. The solution describes the tidal dynamics of a coastal aquifer that is adjacent to a cross-shore estuary. Both the effects of oceanic and estuarine tides on the aquifer are included in the solution. The analytical prediction of the head fluctuations is verified by comparison with numerical solutions computed using a standard finite-difference method. An essential feature of the present analytical solution is the interaction between the cross- and along-shore tidal waves in the aquifer area near the estuary's entry. As the distance from the estuary or coastline increases, the wave interaction is weakened and the aquifer response is reduced, respectively, to the one-dimensional solution for oceanic tides or the solution of Sun (Sun H. A two-dimensional analytical solution of groundwater response to tidal loading in an estuary, Water Resour Res 1997;33:1429-35) for two-dimensional non-interacting tidal waves. (C) 2000 Elsevier Science Ltd. All rights reserved.

Water wave-driven seepage in marine sediments

The action of water waves moving over a porous seabed drives a seepage flux into and out of the marine sediments. The volume of fluid exchange per wave cycle may affect the rate of contaminant transport in the sediments. In this paper, the dynamic response of the seabed to ocean waves is treated analytically on the basis of pore-elastic theory applied to a porous seabed. The seabed is modelled as a semi-infinite, isotropic, homogeneous material. Most previous investigations on the wave-seabed interaction problem have assumed quasi-static conditions within the seabed, although dynamic behaviour often occurs in natural environments. Furthermore, wave pressures used in the previous approaches were obtained from conventional ocean wave theories: which are based on the assumption of an impermeable rigid seabed. By introducing a complex wave number, we derive a new wave dispersion equation, which includes the seabed characteristics (such as soil permeability, shear modulus, etc.). Based on the new closed-form analytical solution, the relative differences of the wave-induced seabed response under dynamic and quasi-static conditions are examined. The effects of wave and soil parameters on the seepage flux per wave cycle are also discussed in detail. (C) 2000 Elsevier Science Ltd. All rights reserved.

The q-deformed supersymmetric t-J model with a boundary

The q-deformed supersymmetric t-J model on a semi-infinite lattice is diagonalized by using the level-one vertex operators of the quantum affine superalgebra U-q[sl(2\1)]. We. give the bosonization of the boundary states. We give an integral expression for the correlation functions of the boundary model, and derive the difference equations which they satisfy.

A simple solution of the laminar dam break wave

The sudden release of a mass of fluid in a channel generates a highly unsteady flow motion, called dam break wave. While industrial fluids exhibit sometimes non-Newtonian behaviours, the viscous fluid flow assumption remains a useful approximation for simplified analyses. In this study, new solutions of laminar dam break wave are proposed for a semi-infinite reservoir based upon the method of characteristics. The solutions yield simple explicit expressions of the wave front location, wave front celerity and instantaneous free-surface profiles that compare favourably with experimental observations. Both horizontal and sloping channel configurations are treated. The simplicity of the equations may allow future extension to more complicated fluid flows.

Izergin-Korepin model with a boundary

The Izergin-Korepin model on a semi-infinite lattice is diagonalized by using the level-one vertex operators of the twisted quantum affine algebra U-q[((2))(2)]. We give the bosonization of the vacuum state with zero particle content. Excitation states are given by the action of the vertex operators on the vacuum state. We derive the boundary S-matrix. We give an integral expression of the correlation functions of the boundary model, and derive the difference equations which they satisfy. (C) 2001 Elsevier Science B.V. All rights reserved.

Invariant integrals applied to nematic liquid crystals with small Ericksen number and topological defects

Invariant integrals are derived for nematic liquid crystals and applied to materials with small Ericksen number and topological defects. The nematic material is confined between two infinite plates located at y = -h and y = h (h is an element of R+) with a semi-infinite plate at y = 0 and x < 0. Planar and homeotropic strong anchoring boundary conditions to the director field are assumed at these two infinite and semi-infinite plates, respectively. Thus, a line disclination appears in the system which coincides with the z-axis. Analytical solutions to the director field in the neighbourhood of the singularity are obtained. However, these solutions depend on an arbitrary parameter. The nematic elastic force is thus evaluated from an invariant integral of the energy-momentum tensor around a closed surface which does not contain the singularity. This allows one to determine this parameter which is a function of the nematic cell thickness and the strength of the disclination. Analytical solutions are also deduced for the director field in the whole region using the conformal mapping method. (C) 2013 Elsevier Ltd. All rights reserved.

Critical Velocity obtained using Simplified Models of the Railway Track: Viability and Applicability

Increased demands on the capacity of the railway network gave rise to new issues related to the dynamic response of railway tracks subjected to moving vehicles. Thus, it becomes important to evaluate the applicability of traditionally used simplified models which have a closed form solution. Regarding simplified models, transversal vibrations of a beam on a visco-elastic foundation subjected to a moving load are considered. Governing equations are obtained by Hamilton’s principle. Shear distortion, rotary inertia and effect of axial force are accounted for. The load is introduced as a time varying force moving at a constant velocity. Transversal vibrations induced by the load are solved by the normal-mode analysis. Reflected waves at the extremities of the full beam are avoided by introduction of semi-infinite elements. Firstly, the critical velocity obtained from this model is compared with results of an undamped Euler- Bernoulli formulation with zero axial force. Secondly, a finite element model in ABAQUS is examined. The new contribution lies in the introduction of semi- infinite elements and in the first step to a systematic comparison, which have not been published so fa

The application of Fick's law to evaluate the period of stay of mega-olivines in alkaline lavas

The integration of the differential equation of the second law of Fick applied to the diffusion of chemical elements in a semi-infinite solid made it easier to estimate the time of stay of olivine mega-cristals in contact with the host lava The results of this research show the existence of two groups of olivine. The first remained in contact with the magmatic liquid during 19 to 22 days, while the second remained so during only 5 to 9 days. This distinction is correlative to that based on the qualitative observation.

Semiclassical image potential at a solid surface

A surface dielectric function of a semi-infinite plane-bounded metal is defined in the spirit of the plasmon-pole dielectric function of the bulk. It is modeled in such a way that the surface-plasmon dispersion relation is recovered for small momentum transfer. This function is employed to compute the image potential at all distances outside the surface. Interaction with bulk modes is neglected for simplicity and clarity. The interaction of a massive point charge with a metal surface is also considered in the context of a boson model for surface-plasmon excitation. We present a new definition of the image potential for this case.

Variational Wigner-Kirkwood approach to relativistic mean field theory

The recently developed variational Wigner-Kirkwood approach is extended to the relativistic mean field theory for finite nuclei. A numerical application to the calculation of the surface energy coefficient in semi-infinite nuclear matter is presented. The new method is contrasted with the standard density functional theory and the fully quantal approach.

Solutions of the telegrapher's equation in the presence of traps

Several problems in the theory of photon migration in a turbid medium suggest the utility of calculating solutions of the telegrapher¿s equation in the presence of traps. This paper contains two such solutions for the one-dimensional problem, the first being for a semi-infinite line terminated by a trap, and the second being for a finite line terminated by two traps. Because solutions to the telegrapher¿s equation represent an interpolation between wavelike and diffusive phenomena, they will exhibit discontinuities even in the presence of traps.