An extended finite-element model coupled with level set method for analysis of growth of corrosion pits in metallic structures

Mass balance between metal and electrolytic solution, separated by a moving interface, in stable pit growth results in a set of governing equations which are solved for concentration field and interface position (pit boundary evolution). The interface experiences a jump discontinuity in metal concentration. The extended finite-element model (XFEM) handles this jump discontinuity by using discontinuous-derivative enrichment formulation, eliminating the requirement of using front conforming mesh and re-meshing after each time step as in the conventional finite-element method. However, prior interface location is required so as to solve the governing equations for concentration field for which a numerical technique, the level set method, is used for tracking the interface explicitly and updating it over time. The level set method is chosen as it is independent of shape and location of the interface. Thus, a combined XFEM and level set method is developed in this paper. Numerical analysis for pitting corrosion of stainless steel 304 is presented. The above proposed model is validated by comparing the numerical results with experimental results, exact solutions and some other approximate solutions. An empirical model for pitting potential is also derived based on the finite-element results. Studies show that pitting profile depends on factors such as ion concentration, solution pH and temperature to a large extent. Studying the individual and combined effects of these factors on pitting potential is worth knowing, as pitting potential directly influences corrosion rate.

Finite Element Analysis of Stress Evolution in Al-Si Alloy

A 2D multi-particle model is carried out to understand the effect of microstructural variations and loading conditions on the stress evolution in Al-Si alloy under compression. A total of six parameters are varied to create 26 idealized microstructures: particle size, shape, orientation, matrix temper, strain rate, and temperature. The effect of these parameters is investigated to understand the fracture of Si particles and the yielding of Al matrix. The Si particles are modeled as a linear elastic solid and the Al matrix is modeled as an elasto-plastic solid. The results of the study demonstrate that the increase in particle size decreases the yield strength of the alloy. The particles with high aspect ratio and oriented at 0A degrees and 90A degrees to the loading axis show higher stress values. This implies that the particle shape and orientation are dominant factors in controlling particle fracture. The heat treatment of the alloy is found to increase the stress levels of both particles and matrix. Stress calculations also show that higher particle fracture and matrix yielding is expected at higher strain rate deformation. Particle fracture decreases with increase in temperature and the Al matrix plays an important role in controlling the properties of the alloy at higher temperatures. Further, this strain rate and temperature dependence is more pronounced in the heat-treated microstructure. These predictions are consistent with the experimentally observed Si particle fracture in real microstructure.

Generation of Temperature Dependent Transversely Isotropic Properties for Zigzag and Armchair Single-Walled Carbon Nanotubes

Finite element analysis has been carried out to obtain temperature dependent transversely isotropic properties of the single-walled carbon nanotubes (SWCNTs). Finite element models of SWCNTs are generated by specifying the C-C bond rigidities. The five independent transversely isotropic properties for different chiralities are evaluated using the stress fields of thick-walled cylinders and the elastic deformations of SWCNTs subjected to pure extension, internal pressure and pure torsion loads. Empirical relations are provided for the five independent elastic constants useful to armchair and zigzag SWCNTs.

Semi-analytical finite element analysis of a strip footing on an elastic reinforced soil

A plane strain elastic interaction analysis of a strip footing resting on a reinforced soil bed has been made by using a combined analytical and finite element method (FEM). In this approach the stiffness matrix for the footing has been obtained using the FEM, For the reinforced soil bed (halfplane) the stiffness matrix has been obtained using an analytical solution. For the latter, the reinforced zone has been idealised as (i) an equivalent orthotropic infinite strip (composite approach) and (ii) a multilayered system (discrete approach). In the analysis, the interface between the strip footing and reinforced halfplane has been assumed as (i) frictionless and (ii) fully bonded. The contact pressure distribution and the settlement reduction have been given for different depths of footing and scheme of reinforcement in soil. The load-deformation behaviour of the reinforced soil obtained using the above modelling has been compared with some available analytical and model test results. The equivalent orthotropic approach proposed in this paper is easy to program and is shown to predict the reinforcing effects reasonably well.

Unsteady flow and heat transfer between rotating coaxial disks

A study is made on the flow and heat transfer of a viscous fluid confined between two parallel disks. The disks are allowed to rotate with different time dependent angular velocities, and the upper disk is made to approach the lower one with a constant speed. Numerical solutions of the governing parabolic partial differential equations are obtained through a fourth-order accurate compact finite difference scheme. The normal forces and torques that the fluid exerts on the rotating surfaces are obtained at different nondimensional times for different values of the rate of squeezing and disk angular velocities. The temperature distribution and heat transfer are also investigated in the present analysis.

Low-temperature electrical conductivity of LaNi1-xFexO3

We report the electrical conductivity between 2 and 300 K for LaNi1-xFexO3 across the composition-controlled metal-insulator (m-i) transition. Using a method first suggested by Mobius, we identify the critical concentration x(c) to be 0.3 for the m-i transition. The negative temperature coefficient of resistivity observed at low temperatures in the metallic phase follows a temperature dependence characteristic of disorder effects. The semiconducting compositions (x greater than or equal to 0.3) do not show a simple activation energy but exhibit variable-range hopping at high temperatures confirming that the m-i transition in this system is driven by increasing disorder effects.

Unsteady Compressible Boundary-Layer Flow On The Stagnation Line Of An Infinite Swept Cylinder

Numerical solutions of flow and heat transfer process on the unsteady flow of a compressible viscous fluid with variable gas properties in the vicinity of the stagnation line of an infinite swept cylinder are presented. Results are given for the case where the unsteady temperature field is produced by (i) a sudden change in the wall temperature (enthalpy) as the impulsive motion is started and (ii) a sudden change in the free-stream velocity. Solutions for the simultaneous development of the thermal and momentum boundary layers are obtained by using quasilinearization technique with an implicit finite difference scheme. Attention is given to the transient phenomenon from the initial flow to the final steady-state distribution. Results are presented for the skin friction and heat transfer coefficients as well as for the velocity and enthalpy profiles. The effects of wail enthalpy parameter, sweep parameter, fluid properties and transpiration cooling on the heat transfer and skin friction are considered.

High temperature deformation and densification in alumina-yttria composites

Pinning by second phase particles offers a potent means for limiting grain growth and enhancing superplasticity in alumina-based ceramics. In the present study, a colloidal technique was used to produce green bodies of alumina-yttria composites; at elevated temperatures, the yttria particles react with alumina to produce YAG particles. The densification and high temperature deformation characteristics of alumina-YAG composites were studied using conventional free sintering and sinter-forging, which involves the application of a compressive stress without any lateral constraints. It is shown that the YAG particles retard both densification and grain growth. The experiments indicate also that the presence of YAG particles does not significantly alter the stress exponent for creep deformation.

Finite volume evolution Galerkin method for hyperbolic conservation laws with spatially varying flux functions

We present a generalization of the finite volume evolution Galerkin scheme [M. Lukacova-Medvid'ova,J. Saibertov'a, G. Warnecke, Finite volume evolution Galerkin methods for nonlinear hyperbolic systems, J. Comp. Phys. (2002) 183 533-562; M. Luacova-Medvid'ova, K.W. Morton, G. Warnecke, Finite volume evolution Galerkin (FVEG) methods for hyperbolic problems, SIAM J. Sci. Comput. (2004) 26 1-30] for hyperbolic systems with spatially varying flux functions. Our goal is to develop a genuinely multi-dimensional numerical scheme for wave propagation problems in a heterogeneous media. We illustrate our methodology for acoustic waves in a heterogeneous medium but the results can be generalized to more complex systems. The finite volume evolution Galerkin (FVEG) method is a predictor-corrector method combining the finite volume corrector step with the evolutionary predictor step. In order to evolve fluxes along the cell interfaces we use multi-dimensional approximate evolution operator. The latter is constructed using the theory of bicharacteristics under the assumption of spatially dependent wave speeds. To approximate heterogeneous medium a staggered grid approach is used. Several numerical experiments for wave propagation with continuous as well as discontinuous wave speeds confirm the robustness and reliability of the new FVEG scheme.

Temperature dependent vibrational lifetimes in supercritical fluids near the critical point

Vibrational relaxation measurements on the CO asymmetric stretching mode (similar to 1980 cm(-1)) of tungsten hexacarbonyl (W(CO)(6)) as a function of temperature at constant density in several supercritical solvents in the vicinity of the critical point are presented. In supercritical ethane, at the critical density, there is a region above the critical temperature (Tc) in which the lifetime increases with increasing temperature. When the temperature is raised sufficiently (similar to T-c + 70 degrees C), the lifetime decreases with further increase in temperature. A recent hydrodynamic/thermodynamic theory of vibrational relaxation in supercritical fluids reproduces this behavior semiquantitatively. The temperature dependent data for fixed densities somewhat above and below the critical density is in better agreement with the theory. In fluoroform solvent at the critical density, the vibrational lifetime also initially increases with increasing temperature. However, in supercritical CO2 at the critical density, the temperature dependent vibrational lifetime decreases approximately linearly with temperature beginning almost immediately above T-c. The theory does not reproduce this behavior. A comparison between the absolute lifetimes in the three solvents and the temperature trends is made.

Virtues of marginally metallic oxides

Transition-metal oxides at the metal-insulator boundary, especially those belonging to the perovskite family, exhibit fascinating phenomena such as insulator-metal transitions controlled by composition, high-temperature superconductivity and giant magnetoresistance (GMR), Interestingly, many of these marginally metallic oxides obey the established criteria for metallicity and have a finite density of states at the Fermi;level. The perovskite manganates exhibiting GMR, on the other hand, are unusual in that they possess very high resistivities in the 'metallic' state and show no significant density of states at the Fermi level, Marginal metallicity in oxide systems is a problem of great complexity and contemporary interest and its understanding is of crucial significance to the diverse phenomena exhibited by these materials.

The high temperature mechanical characteristics of superplastic 3 mol% yttria stabilized zirconia

A detailed study was undertaken to characterize the deformation behavior of a superplastic 3 mol% yttria-stabilized tetragonal zirconia (3YTZ) over a wide range of strain rates, temperatures and grain sizes. The experimental data were analyzed in terms of the following equation for high temperature deformation: Image Full-size image ∞ σn d−pexp(−Q/RT), where Image Full-size image is the strain rate, σ is the flow stress, d is the grain size, Q is the activation energy, R is the gas constant, T is the absolute temperature, and n and p are constants termed the stress exponent and the inverse grain size exponent, respectively. The experimental data over a wide range of stresses revealed a transition in stress exponent. Deformation in the low and high stress regions was associated with n not, vert, similar 3 and p not, vert, similar 1, and n not, vert, similar 2 and p not, vert, similar 3, respectively. The transition stress between the two regions decreased with increasing grain size. The activation energy was similar for both regions with a value of not, vert, similar 550 kJ mol−1. Microstructural measurements revealed that grains remained essentially equiaxed after the accumulation of large strains, and very limited concurrent grain growths occurred in most experiments. Assessment of possible rate controlling creep mechanisms and comparison with previous studied indicate that in the n not, vert, similar 2 region, deformation occurs by a grain boundary sliding process whose rate is independent of impurity content. Deformation in the n not, vert, similar 3 region is controlled by an interface reaction that is highly sensitive to impurity content. It is concluded that an increase in impurity content increases yttrium segregation to grain boundaries, which enhances the rate of the interface reaction, thereby decreasing the apparent transition stress between the n not, vert, similar 2 and n not, vert, similar 3 regions. This unified approach incorporating two sequential mechanisms can rationalize many of the apparently dissimilar results that have been reported previously for deformation of 3YTZ.

Electrical Conductivity and Magneto-Resistance Studies on Polyaniline-Pmma Blends at Very Low Temperature

D.C. electrical conductivity of polyaniline (33%,40%) blended with PMMA was measured from 5K to 300mK. The conductivity behaviour is consistent with fluctuation induced tunneling. Magneto-resistance (MR) was measured between 300K and 2K. From 20K to 2K, a large positive MR was observed. At 2K, for low magnetic fields (<1 Tesla), a deviation from the normal H-2 behaviour was observed.

Unsteady free convection flow in the stagnation-point region of a rotating sphere

The unsteady free convection boundary-layer flow in the forward stagnation-point region of a sphere, which is rotating with time-dependent angular velocity in an ambient fluid, has been studied. Both constant wall temperature and constant hear flux conditions have been considered. The non-linear coupled parabolic partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. The skin friction and the heat transfer are enhanced by the buoyancy force. The effect of the buoyancy force is found to be more pronounced for smaller Prandtl numbers than for larger Prandtl numbers. For a given buoyancy force, the heat transfer increases with an increase in Prandtl number, but the skin friction decreases.

Unsteady free convection flow in the stagnation-point region of a three-dimensional body

The unsteady free convection flow in the stagnation-point region of a heated three-dimensional body placed in an ambient fluid is studied under boundary layer approximations. We have considered the case where there is an initial steady state that is perturbed by a step-change in the wall temperature. The non-linear coupled partial differential equations governing the free convection flow are solved numerically using a finite difference scheme. The presented results show the temporal development of the momentum and thermal boundary layer characteristics.