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.

Study of soils by thermal analysis

The study of soils is very important in the geological and geological engineering researches. A study of ten samples of soils was carried out by thermal analysis, and X-Ray Fluorescence Spectrometry to understand soil evolution in Angra dos Reis region, Rio de Janeiro State, Brazil. The sample collection sites were chosen based on geological characteristics, the soil layer thickness, the soil composition pattern, and whether or not it was moved either by erosion or by gravitational shifts. Because of the humid tropical climatic condition, natural soils tend to show great thickness of weathered mantles with formation of saprolites and saprolite soils. Kaolinite is an important secondary mineral which can be formed from many different minerals, like k-mica and k-feldspar and can be weathered to gibbsite. The results from TG/DTG and DTA indicated which soils had more weathering, and the same results were obtained by XRF, when silica/aluminum ratios from samples are compared with thermal analysis results.

Study of the water retention properties of a tropical soil

Some peculiarities of water retention in a tropical lateritic soil of clayey nature are presented and discussed. The typical soil microstructure is shown through thin-layer plates emphasizing soil microaggregation and pore distribution and their repercussion on the soil-water retention curve and on hysteresis. It is shown that the clayey soil has a behavior that to a large extent resembles sandy soil, which is characterized by the relatively high saturated hydraulic conductivity, low air-entry value, and small suction range at which water drainage takes place. The severe weathering processes that originated this soil have produced an altered soil that seems to be homogeneous in terms of physical indices, hydraulic conductivity, and soil-water retention characteristics, up to 4.5 m in depth.

Solution of a heterogeneous modeling of a horizontal-flow anaerobic immobilized biomass (HAIB) reactor by the sequencing method

A modeling study was completed to develop a methodology that combines the sequencing and finite difference methods for the simulation of a heterogeneous model of a tubular reactor applied in the treatment of wastewater. The system included a liquid phase (convection diffusion transport) and a solid phase (diffusion reaction) that was obtained by completing a mass balance in the reactor and in the particle, respectively. The model was solved using a pilot-scale horizontal-flow anaerobic immobilized biomass (HAIB) reactor to treat domestic sewage, with the concentration results compared with the experimental data. A comparison of the behavior of the liquid phase concentration profile and the experimental results indicated that both the numerical methods offer a good description of the behavior of the concentration along the reactor. The advantage of the sequencing method over the finite difference method is that it is easier to apply and requires less computational time to model the dynamic simulation of outlet response of HAIB.

A comparative study of GERIPA ethanol with other fuels

The GERIPA project aimed at generating renewable energy integrated with food production has led to a beneficial option for producing ethanol and electricity. Ethanol has economic, social and environmental potential. Considering just the first one, Brazil consumes 39 billion litres per year-L(D)/yr of diesel oil, 18% of it being imported. The Federal Government has a recovery programme for the soybean agribusiness aimed at soybean biodiesel (SBD) production in which a 10% addition to diesel has been proposed. This 10% involves producing 10.7 million L(SB)/d. Soybean bio-diesel production is not self-sustainable and such proposal could require an annual subsidy of up to US$1.33 billion. Soybean plantations would need about 10 to 12 times more land than is necessary for sugarcane plantations to produce the same equivalent thermal energy (ETE). Sixty-seven GERIPA projects (GP) producing 80,000 litres of ethanol per day (GP80) could be set up with the sum of US$1.33 billion; this would substitute current Brazilian biodiesel demand by 4.28%, adding the some value for each new subsidiary. Considering ETE, ethanol-GP cost would be 37% to 50% below that for a litre of SBD on account of its raw material (RM) and region. The diesel cycle`s thermal efficiency (eta(1)) yield is around 50% and that of the Otto cycle engine eta(1) is around 37%. The cost per km driven (CKD) by substituting SBD for ethanol-GP80 would thus indicate an 18% minimum and 59% maximum cost reduction for vehicle engines.

Modeling and experimental study of nucleate boiling on a vertical array of horizontal plain tubes

An investigation of nucleate boiling on a vertical array of horizontal plain tubes is presented in this paper. Experiments were performed with refrigerant RI 23 at reduced pressures varying from 0.022 to 0.64, tube pitch to diameter ratios of 1.32, 1.53 and 2.00, and heat fluxes from 0.5 to 40 kW/m(2). Brass tubes with external diameters of 19.05 mm and average roughness of 0.12 mu m were used in the experiments. The effect of the tube spacing on the local heat transfer coefficient along the tube array was negligible within the present range of experimental conditions. For partial nucleate boiling, characterized by low heat fluxes, and low reduced pressures, the tube positioning shows a remarkable effect on the heat transfer coefficient. Based on these data, a general correlation for the prediction of the nucleate boiling heat transfer coefficient on a vertical array of horizontal tubes under flooded conditions was proposed. According to this correlation, the ratio between the heat transfer coefficients of a given tube and the lowest tube in the array depends only on the tube row number, the reduced pressure and the heat flux. By using the proposed correlation, most of the experimental heat transfer coefficients obtained in the present study were predicted within +/- 15%. The new correlation compares reasonably well with independent data from the literature. (C) 2008 Elsevier Inc. All rights reserved.

Wavelet-Galerkin method for one-dimensional elastoplasticity and damage problems: Constitutive modeling and computational aspects

This work presents an analysis of the wavelet-Galerkin method for one-dimensional elastoplastic-damage problems. Time-stepping algorithm for non-linear dynamics is presented. Numerical treatment of the constitutive models is developed by the use of return-mapping algorithm. For spacial discretization we can use wavelet-Galerkin method instead of standard finite element method. This approach allows to locate singularities. The discrete formulation developed can be applied to the simulation of one-dimensional problems for elastic-plastic-damage models. (C) 2007 Elsevier Inc. All rights reserved.

Experimental study of a Tollmien-Schlichting wave interacting with a shallow 3D roughness element

The paper is devoted to an experimental study of the effect of a shallow 3D roughness element on the evolution of a 2D Tollmien-Schlichting wave in a Blasius boundary layer. The experiments were carried out under controlled disturbance conditions on an airfoil section which could provide a long run with zero pressure gradient flow. A pneumatically driven slit source was used to introduce the Tollmien-Schilichting wave upstream of the lower branch of the neutral stability curve. A few wavelengths downstream, the T-S wave interacts with a cylindrical roughness element. The height of the roughness was slowly oscillating in time, which allows a continuous measurement of the T-S wave response downstream the roughness. The oscillation frequency was approximately 1500 times lower than the frequency of the studied Tollmien-Schlichting wave and therefore, behaved as a steady roughness with respect to the T-S wave. Hot wire anemometry was used to measure wall normal profiles and spanwise scans close to the maximum of the eigenfunction of the T-S wave. The oscillation of the roughness and the synchronization of all-equipments permitted the use of ensemble average techniques. Two different amplitudes of T-S waves with a non-dimensional frequency of F120E-06 were studied. They show a strong amplification of the disturbances in a small spanwise wave number range. The analysis of the wall normal T-S profiles suggests the growth of oblique modes.

Fuzzy computational control for real Chua circuit

We preserit a computational procedure to control art experimental chaotic system by applying the occasional proportional feedback (OPF) method. The method implementation uses the fuzzy theory to relate the variable correction to the necessary adjustment in the control parameter. As an application We control the chaotic attractors of the Chua circuit. We present file developed circuits and algorithms to implement this control in real time. To simplify the used procedure, we use it low resolution analog to digital converter compensated for a lowpass filter that facilitates similar applications to control other systems. (C) 2007 Elsevier Ltd. All rights reserved.

A Study of Influence of Configurations of the Windings on the Performance of Single Phase Induction Motor with Capacitor Start

The single phase induction motors needs two stator windings to produce rotating magnetic field : one main winding and the other auxiliary winding. The aim of the auxiliary winding is to create the rotating electromagnetic field when the machine is started-up and is afterwards turned off, generally through the centrifugal switch coupled together with the shaft of the machine rotor. The main purpose of this document is to evaluate the influence that the two windings have on the external characteristics of the single phase induction motor. For this purpose, two different kinds of windings were carried out and simulated, with the proposal to obtain some benefits. The main winding and the auxiliary winding were prepared and mounted on a prototype. The simulation was done via software based FEM, to make the extraction and results analysis possible. This results are shown at the end this document.

Technology assessment for power quality mitigation devices - Micro-DVR case study

This work presents a case study on technology assessment for power quality improvement devices. A system compatibility test protocol for power quality mitigation devices was developed in order to evaluate the functionality of three-phase voltage restoration devices. In order to validate this test protocol, the micro-DVR, a reduced power development platform for DVR (dynamic voltage restorer) devices, was tested and the results are discussed based on voltage disturbances standards. (C) 2011 Elsevier B.V. All rights reserved.

Finite element computational modeling of externally bonded CFRP composites flexural behavior in RC beams

This paper focuses on the flexural behavior of RC beams externally strengthened with Carbon Fiber Reinforced Polymers (CFRP) fabric. A non-linear finite element (FE) analysis strategy is proposed to support the beam flexural behavior experimental analysis. A development system (QUEBRA2D/FEMOOP programs) has been used to accomplish the numerical simulation. Appropriate constitutive models for concrete, rebars, CFRP and bond-slip interfaces have been implemented and adjusted to represent the composite system behavior. Interface and truss finite elements have been implemented (discrete and embedded approaches) for the numerical representation of rebars, interfaces and composites.

On the finite element modeling of fatigue crack growth in pressurized cylindrical shells

A methodology for the computational modeling of the fatigue crack growth in pressurized shell structures, based on the finite element method and concepts of Linear Elastic Fracture Mechanics, is presented. This methodology is based on that developed by Potyondy [Potyondy D, Wawrzynek PA, Ingraffea, AR. Discrete crack growth analysis methodology for through crack in pressurized fuselage structures. Int J Numer Methods Eng 1995;38:1633-1644], which consists of using four stress intensity factors, computed from the modified crack integral method, to predict the fatigue propagation life as well as the crack trajectory, which is computed as part of the numerical simulation. Some issues not presented in the study of Potyondy are investigated herein such as the influence of the crack increment size and the number of nodes per element (4 or 9 nodes) on the simulation results by means of a fatigue crack propagation simulation of a Boeing 737 airplane fuselage. The results of this simulation are compared with experimental results and those obtained by Potyondy [1]. (C) 2008 Elsevier Ltd. All rights reserved.

Computational evaluation of flexural toughness of FRC and fracture properties of plain concrete

In this paper, a computational tool concerning the computation of flexural and fracture toughness of cement based composites is presented. Firstly, RILEM`s (Reunion Internationale des Laboratoires d`Essais de Materiaux) recommendations related to the analysis of FRC in three-point bend tests are discussed in their relevant aspects regarding the computational implementations. The determination of other mechanical properties such as the Young modulus has been added to the program. Taking this into account, a new formulation based on displacements is used. In the second part of the paper, the determination of fracture properties of concrete, such as the fracture energy, G(F) , and the fracture toughness, K-IC(S), is discussed regarding the computational strategies used in the implementations. Several features whereby anterior data can be reanalyzed, obtained from other standards and recommendations, have been incorporated into the program, therefore allowing comparative studies and back analysis activities.