A Study on Unconditionally Stable FDTD Methods for the Modeling of Metamaterials

We assess the performance of three unconditionally stable finite-difference time-domain (FDTD) methods for the modeling of doubly dispersive metamaterials: 1) locally one-dimensional FDTD; 2) locally one-dimensional FDTD with Strang splitting; and (3) alternating direction implicit FDTD. We use both double-negative media and zero-index media as benchmarks.

Robot measuring form errors

This work presents an automated system for the measurement of form errors of mechanical components using an industrial robot. A three-probe error separation technique was employed to allow decoupling between the measured form error and errors introduced by the robotic system. A mathematical model of the measuring system was developed to provide inspection results by means of the solution of a system of linear equations. A new self-calibration procedure, which employs redundant data from several runs, minimizes the influence of probes zero-adjustment on the final result. Experimental tests applied to the measurement of straightness errors of mechanical components were accomplished and demonstrated the effectiveness of the employed methodology. (C) 2007 Elsevier Ltd. All rights reserved.

Non-linear boundary element formulation applied to contact analysis using tangent operator

This work presents a non-linear boundary element formulation applied to analysis of contact problems. The boundary element method (BEM) is known as a robust and accurate numerical technique to handle this type of problem, because the contact among the solids occurs along their boundaries. The proposed non-linear formulation is based on the use of singular or hyper-singular integral equations by BEM, for multi-region contact. When the contact occurs between crack surfaces, the formulation adopted is the dual version of BEM, in which singular and hyper-singular integral equations are defined along the opposite sides of the contact boundaries. The structural non-linear behaviour on the contact is considered using Coulomb`s friction law. The non-linear formulation is based on the tangent operator in which one uses the derivate of the set of algebraic equations to construct the corrections for the non-linear process. This implicit formulation has shown accurate as the classical approach, however, it is faster to compute the solution. Examples of simple and multi-region contact problems are shown to illustrate the applicability of the proposed scheme. (C) 2011 Elsevier Ltd. All rights reserved.

A boundary element formulation for analysis of elastoplastic plates with geometrical nonlinearity

In this paper a new boundary element method formulation for elastoplastic analysis of plates with geometrical nonlinearities is presented. The von Mises criterion with linear isotropic hardening is considered to evaluate the plastic zone. Large deflections are assumed but within the context of small strain. To derive the boundary integral equations the von Karman`s hypothesis is taken into account. An initial stress field is applied to correct the true stresses according to the adopted criterion. Isoparametric linear elements are used to approximate the boundary unknown values while triangular internal cells with linear shape function are adopted to evaluate the domain value influences. The nonlinear system of equations is solved by using an implicit scheme together with the consistent tangent operator derived along the paper. Numerical examples are presented to demonstrate the accuracy and the validity of the proposed formulation.

Non-linear boundary element formulation with tangent operator to analyse crack propagation in quasi-brittle materials

This work deals with analysis of cracked structures using BEM. Two formulations to analyse the crack growth process in quasi-brittle materials are discussed. They are based on the dual formulation of BEM where two different integral equations are employed along the opposite sides of the crack surface. The first presented formulation uses the concept of constant operator, in which the corrections of the nonlinear process are made only by applying appropriate tractions along the crack surfaces. The second presented BEM formulation to analyse crack growth problems is an implicit technique based on the use of a consistent tangent operator. This formulation is accurate, stable and always requires much less iterations to reach the equilibrium within a given load increment in comparison with the classical approach. Comparison examples of classical problem of crack growth are shown to illustrate the performance of the two formulations. (C) 2009 Elsevier Ltd. All rights reserved.

Consolidation of elastic-plastic saturated porous media by the boundary element method

This paper presents a domain boundary element formulation for inelastic saturated porous media with rate-independent behavior for the solid skeleton. The formulation is then applied to elastic-plastic behavior for the solid. Biot`s consolidation theory, extended to include irreversible phenomena is considered and the direct boundary element technique is used for the numerical solution after time discretization by the implicit Euler backward algorithm. The associated nonlinear algebraic problem is solved by the Newton-Raphson procedure whereby the loading/unloading conditions are fully taken into account and the consistent tangent operator defined. Only domain nodes (nodes defined inside the domain) are used to represent all domain values and the corresponding integrals are computed by using an accurate sub-elementation scheme. The developments are illustrated through the Drucker-Prager elastic-plastic model for the solid skeleton and various examples are analyzed with the proposed algorithms. (c) 2008 Elsevier B.V. All rights reserved.

Leak detection by inverse transient analysis in an experimental PVC pipe system

Leakage reduction in water supply systems and distribution networks has been an increasingly important issue in the water industry since leaks and ruptures result in major physical and economic losses. Hydraulic transient solvers can be used in the system operational diagnosis, namely for leak detection purposes, due to their capability to describe the dynamic behaviour of the systems and to provide substantial amounts of data. In this research work, the association of hydraulic transient analysis with an optimisation model, through inverse transient analysis (ITA), has been used for leak detection and its location in an experimental facility containing PVC pipes. Observed transient pressure data have been used for testing ITA. A key factor for the success of the leak detection technique used is the accurate calibration of the transient solver, namely adequate boundary conditions and the description of energy dissipation effects since PVC pipes are characterised by a viscoelastic mechanical response. Results have shown that leaks were located with an accuracy between 4-15% of the total length of the pipeline, depending on the discretisation of the system model.

Void fraction measurement and signal analysis from multiple-electrode impedance sensors

Void fraction sensors are important instruments not only for monitoring two-phase flow, but for furnishing an important parameter for obtaining flow map pattern and two-phase flow heat transfer coefficient as well. This work presents the experimental results obtained with the analysis of two axially spaced multiple-electrode impedance sensors tested in an upward air-water two-phase flow in a vertical tube for void fraction measurements. An electronic circuit was developed for signal generation and post-treatment of each sensor signal. By phase shifting the electrodes supplying the signal, it was possible to establish a rotating electric field sweeping across the test section. The fundamental principle of using a multiple-electrode configuration is based on reducing signal sensitivity to the non-uniform cross-section void fraction distribution problem. Static calibration curves were obtained for both sensors, and dynamic signal analyses for bubbly, slug, and turbulent churn flows were carried out. Flow parameters such as Taylor bubble velocity and length were obtained by using cross-correlation techniques. As an application of the void fraction tested, vertical flow pattern identification could be established by using the probability density function technique for void fractions ranging from 0% to nearly 70%.

Computer-Assisted Numerical Analysis for Oxygen and Carbon Dioxide Mass Transfer in Blood Oxygenators

A two-dimensional numeric simulator is developed to predict the nonlinear, convective-reactive, oxygen mass exchange in a cross-flow hollow fiber blood oxygenator. The numeric simulator also calculates the carbon dioxide mass exchange, as hemoglobin affinity to oxygen is affected by the local pH value, which depends mostly on the local carbon dioxide content in blood. Blood pH calculation inside the oxygenator is made by the simultaneous solution of an equation that takes into account the blood buffering capacity and the classical Henderson-Hasselbach equation. The modeling of the mass transfer conductance in the blood comprises a global factor, which is a function of the Reynolds number, and a local factor, which takes into account the amount of oxygen reacted to hemoglobin. The simulator is calibrated against experimental data for an in-line fiber bundle. The results are: (i) the calibration process allows the precise determination of the mass transfer conductance for both oxygen and carbon dioxide; (ii) very alkaline pH values occur in the blood path at the gas inlet side of the fiber bundle; (iii) the parametric analysis of the effect of the blood base excess (BE) shows that V(CO2) is similar in the case of blood metabolic alkalosis, metabolic acidosis, or normal BE, for a similar blood inlet P(CO2), although the condition of metabolic alkalosis is the worst case, as the pH in the vicinity of the gas inlet is the most alkaline; (iv) the parametric analysis of the effect of the gas flow to blood flow ratio (Q(G)/Q(B)) shows that V(CO2) variation with the gas flow is almost linear up to Q(G)/Q(B) = 2.0. V(O2) is not affected by the gas flow as it was observed that by increasing the gas flow up to eight times, the V(O2) grows only 1%. The mass exchange of carbon dioxide uses the full length of the hollow-fiber only if Q(G)/Q(B) > 2.0, as it was observed that only in this condition does the local variation of pH and blood P(CO2) comprise the whole fiber bundle.

Avaliação do desempenho ambiental de obras de recuperação de rodovias

This paper features a tool to carry out environmental performance evaluation in highway rehabilitation works, as a component of environmental supervision activities. The procedure involves (i) evidence gathering by conducting technical inspections and reviewing environmental compliance reports, (ii) ranking nonconformities according to a proposed weighting framework, and (iii) calculation of an environmental conformity index. For the sake of testing and calibration, the procedure was applied to five road segments that were submitted to rehabilitation works in Sao Paulo State and the results were treated qualitatively.

Estudo de variabilidade e otimização do circuito de moagem SAG da Usina do Sossego

Sossego was the first Vale SAG mill operation to process copper-gold ore. It is located in the State of Para, southeastern Amazon region of Brazil. In the first three years of continuous operation, Vale investigated different alternatives for improving the circuit`s performance by investigating operating conditions, mainly focusing on the SAG mill. It was decided to further assess the performance of the comminution circuit as a function of ore characteristics. A comprehensive ore characterization program was then conducted, together with the calibration of mathematical models on the basis of surveys carried out at the industrial circuit. The simulator was then used to predict the throughput associated to each ore type, as well as to establish the optimized circuit configuration and tailored operating conditions. This paper describes in detail the main aspects of optimizing the industrial circuit performance, as well as the successful method for predicting the production as a function of ore characteristics and circuit configuration.

A procedure to evaluate environmental rehabilitation in limestone quarries

A procedure to evaluate mine rehabilitation practices during the operational phase was developed and validated. It is based on a comparison of actually observed or documented practices with internationally recommended best practices (BP). A set of 150 BP statements was derived from international guides in order to establish the benchmark. The statements are arranged in six rehabilitation programs under three categories: (1) planning (2) operational and (3) management, corresponding to the adoption of the plan-do-check-act management systems model to mine rehabilitation. The procedure consists of (i) performing technical inspections guided by a series of field forms containing BP statements; (ii) classifying evidences in five categories; and (iii) calculating conformity indexes and levels. For testing and calibration purposes, the procedure was applied to nine limestone quarries and conformity indexes were calculated for the rehabilitation programs in each quarry. Most quarries featured poor planning practices, operational practices reached high conformity levels in 50% of the cases and management practices scored moderate conformity. Despite all quarries being ISO 14001 certified, their management systems pay low attention to issues pertaining to land rehabilitation and biodiversity. The best results were achieved by a quarry whose expansion was recently submitted to the environmental impact assessment process, suggesting that public scrutiny may play a positive role in enhancing rehabilitation practices. Conformity indexes and levels can be used to chart the evolution of rehabilitation practices at regular intervals, to establish corporate goals and for communication with stakeholders. (C) 2010 Elsevier Ltd. All rights reserved.

MEASUREMENT OF ELASTIC PROPERTIES OF MATERIALS BY THE ULTRASONIC THROUGH-TRANSMISSION TECHNIQUE

The elastic mechanical behavior of elastic materials is modeled by a pair of independent constants (Young`s modulus and Poisson`s coefficient). A precise measurement for both constants is necessary in some applications, such as the quality control of mechanical elements and standard materials used for the calibration of some equipment. Ultrasonic techniques have been used because wave velocity depends on the elastic properties of the propagation medium. The ultrasonic test shows better repeatability and accuracy than the tensile and indentation test. In this work, the theoretical and experimental aspects related to the ultrasonic through-transmission technique for the characterization of elastic solids is presented. Furthermore, an amorphous material and some polycrystalline materials were tested. Results have shown an excellent repeatability and numerical errors that are less than 3% in high-purity samples.

Evaluation of therapeutic ultrasound equipments performance

The therapeutic ultrasound (US) is one of the resources mostly used by physiotherapists; however the use of uncalibrated equipments results in inefficient or even harmful therapies to the patient. In this direction, the objective of this study was to evaluate the performance and the procedures of utilization and maintenance of US in use in clinics and Physical-therapy offices. A questionnaire with questions related to the procedures applied in service during the use of therapeutic ultrasound was applied to physiotherapists. The performance of 31 equipments of 6 different brands and 13 different models was evaluated according to the IEC 61689 norm. The parameters measured were: acoustic power; effective radiating area (AER); non-uniformity ratio of the beam (RBN); maximum effective intensity; acoustic frequency of operation, modulation factor and wave form on pulsate mode. As for the questionnaires, it was evident that the professionals are not concerned about the calibration of the equipment. The results demonstrated that only 32.3% of the equipments were in accordance with the norms for the variables power and effective radiation area. The frequency analysis indicated that 20% of the 3 MHz transducers and 12.5% of the 1 MHz contemplated the norms. In the pulsate mode, 12.7% presented relation rest/duration inside allowed limits. A great variation of the ultrasonic field was observed on the obtained images, which presented beams not centered, sometimes with bifurcation of its apex. The results allow concluding that, although used in therapeutic sessions with the population, none of the equipments presents all the analyzed variables inside technical norms. (C) 2010 Elsevier B. V. All rights reserved.

Three-dimensional analysis of fracture, corrosion and wear surfaces

A brief look at the history of fractography has shown a recent trend in the quantification of topographic parameters through the use of three-dimensional reconstruction techniques, which associate SEM stereoscopy and stereophotogrammetry software, allowing the calculation of the elevation measurement at numerous points of the topography due to the parallax that takes place during the tilting of the sample along the microscope eucentric plane. Several investigators have used reconstruction techniques to correlate some fractographic parameters, such as fractal dimension and fractured to projected area ratio, to the mechanical properties of materials, such as fracture toughness and tensile strength. So far, the search for a clear relationship between the fracture topography and mechanical properties has provided ambiguous results. The present work applied a surface metrology software to reconstruct three-dimensionally fracture surfaces (transgranular cleavage, intergranular and dimple fracture), corrosion pits and tribo-surfaces in order to explore the potential of this stereophotogrammetry technique. The existence of a variation in the calculated topographic parameters with the conditions of SEM image acquisition reinforces the importance of both good image acquisition and accurate calibration methods in order to validate this 3D reconstruction technique in metrological terms. Preliminary results did not indicate the existence of a clear relationship between either the true to project area ratio and CVN absorbed energy or the fractal dimension and CVN absorbed energy. It is likely that each fracture mechanism presents a proper relationship between the fractographic parameters and mechanical properties. (C) 2009 Elsevier Ltd. All rights reserved.