An extended finite element model coupled with level set method for stable pitting corrosion

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), which requires only three inputs, namely the solid metal concentration, saturation concentration of the dissolved metal ions and diffusion coefficient. A combined eXtended Finite Element Model (XFEM) and level set method is developed in this paper. The extended finite element model handles the jump discontinuity in the metal concentrations at the interface, by using discontinuous-derivative enrichment formulation for concentration discontinuity at the interface. This eliminates the requirement of using front conforming mesh and re-meshing after each time step as in conventional finite element method. A numerical technique known as level set method tracks the position of the moving interface and updates it over time. Numerical analysis for pitting corrosion of stainless steel 304 is presented. The above proposed method is validated by comparing the numerical results with experimental results, exact solutions and some other approximate solutions.

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.

A spectral multiscale method for wave propagation analysis: Atomistic-continuum coupled simulation

In this paper, we present a new multiscale method which is capable of coupling atomistic and continuum domains for high frequency wave propagation analysis. The problem of non-physical wave reflection, which occurs due to the change in system description across the interface between two scales, can be satisfactorily overcome by the proposed method. We propose an efficient spectral domain decomposition of the total fine scale displacement along with a potent macroscale equation in the Laplace domain to eliminate the spurious interfacial reflection. We use Laplace transform based spectral finite element method to model the macroscale, which provides the optimum approximations for required dynamic responses of the outer atoms of the simulated microscale region very accurately. This new method shows excellent agreement between the proposed multiscale model and the full molecular dynamics (MD) results. Numerical experiments of wave propagation in a 1D harmonic lattice, a 1D lattice with Lennard-Jones potential, a 2D square Bravais lattice, and a 2D triangular lattice with microcrack demonstrate the accuracy and the robustness of the method. In addition, under certain conditions, this method can simulate complex dynamics of crystalline solids involving different spatial and/or temporal scales with sufficient accuracy and efficiency. (C) 2014 Elsevier B.V. All rights reserved.

Approximate analysis of coupled non-linear non-conservative systems subjected to step-function excitation

The paper deals with the approximate analysis of non-linear non-conservative systems oftwo degrees of freedom subjected to step-function excitation. The method of averaging of Krylov and Bogoliubov is used to arrive at the approximate equations for amplitude and phase. An example of a spring-mass-damper system is presented to illustrate the method and a comparison with numerical results brings out the validity of the approach.

Transient response of coupled non-linear non-conservative systems

This paper deals with an approximate method of analysis of non-linear, non-conservative systems of two degrees of freedom. The approximate equations for amplitude and phase are obtained by a generalized averaging technique based on the ultraspherical polynomial approximation. The method is illustrated by an example of a spring-mass-damper system.

An FETI-preconditioned conjuerate gradient method for large-scale stochastic finite element problems

In the spectral stochastic finite element method for analyzing an uncertain system. the uncertainty is represented by a set of random variables, and a quantity of Interest such as the system response is considered as a function of these random variables Consequently, the underlying Galerkin projection yields a block system of deterministic equations where the blocks are sparse but coupled. The solution of this algebraic system of equations becomes rapidly challenging when the size of the physical system and/or the level of uncertainty is increased This paper addresses this challenge by presenting a preconditioned conjugate gradient method for such block systems where the preconditioning step is based on the dual-primal finite element tearing and interconnecting method equipped with a Krylov subspace reusage technique for accelerating the iterative solution of systems with multiple and repeated right-hand sides. Preliminary performance results on a Linux Cluster suggest that the proposed Solution method is numerically scalable and demonstrate its potential for making the uncertainty quantification Of realistic systems tractable.

Asymptotic analysis for the coupled wavenumbers in an infinite fluid-filled flexible cylindrical shell: The beam mode

Using asymptotics, the coupled wavenumbers in an infinite fluid-filled flexible cylindrical shell vibrating in the beam mode (viz. circumferential wave order n = 1) are studied. Initially, the uncoupled wavenumbers of the acoustic fluid and the cylindrical shell structure are discussed. Simple closed form expressions for the structural wavenumbers (longitudinal, torsional and bending) are derived using asymptotic methods for low- and high-frequencies. It is found that at low frequencies the cylinder in the beam mode behaves like a Timoshenko beam. Next, the coupled dispersion equation of the system is rewritten in the form of the uncoupled dispersion equation of the structure and the acoustic fluid, with an added fluid-loading term involving a parameter mu due to the coupling. An asymptotic expansion involving mu is substituted in this equation. Analytical expressions are derived for the coupled wavenumbers (as modifications to the uncoupled wavenumbers) separately for low- and high-frequency ranges and further, within each frequency range, for large and small values of mu. Only the flexural wavenumber, the first rigid duct acoustic cut-on wavenumber and the first pressure-release acoustic cut-on wavenumber are considered. The general trend found is that for small mu, the coupled wavenumbers are close to the in vacuo structural wavenumber and the wavenumbers of the rigid-acoustic duct. With increasing mu, the perturbations increase, until the coupled wavenumbers are better identified as perturbations to the pressure-release wavenumbers. The systematic derivation for the separate cases of small and large mu gives more insight into the physics and helps to continuously track the wavenumber solutions as the fluid-loading parameter is varied from small to large values. Also, it is found that at any frequency where two wavenumbers intersect in the uncoupled analysis, there is no more an intersection in the coupled case, but a gap is created at that frequency. This method of asymptotics is simple to implement using a symbolic computation package (like Maple). (C) 2008 Elsevier Ltd. All rights reserved.

Asymptotic analysis for the coupled wavenumbers in an infinite fluid-filled flexible cylindrical shell: The axisymmetric mode

The coupled wavenumbers of a fluid-filled flexible cylindrical shell vibrating in the axisymmetric mode are studied. The coupled dispersion equation of the system is rewritten in the form of the uncoupled dispersion equation of the structure and the acoustic fluid, with an added fluid-loading term involving a parameter e due to the coupling. Using the smallness of Poisson's ratio (v), a double-asymptotic expansion involving e and v 2 is substituted in this equation. Analytical expressions are derived for the coupled wavenumbers (for large and small values of E). Different asymptotic expansions are used for different frequency ranges with continuous transitions occurring between them. The wavenumber solutions are continuously tracked as e varies from small to large values. A general trend observed is that a given wavenumber branch transits from a rigidwalled solution to a pressure-release solution with increasing E. Also, it is found that at any frequency where two wavenumbers intersect in the uncoupled analysis, there is no more an intersection in the coupled case, but a gap is created at that frequency. Only the axisymmetric mode is considered. However, the method can be extended to the higher order modes.

Analysis of solar radiation pressure induced coupled liberations of a gravity stabilized axisymmetric satellite

This paper presents an analysis of solar radiation pressure induced coupled librations of gravity stabilized cylindrical spacecraft with a special reference to geostationary communication satellites. The Lagrangian approach is used to obtain the corresponding equations of motion. The solar induced torques are assumed to be free of librational angles and are represented by their Fourier expansion. The response and periodic solutions are obtained through linear and nonlinear analyses, using the method of harmonic balance in the latter case. The stability conditions are obtained using Routh-Hurwitz criteria. To establish the ranges of validity the analytic response is compared with the numerical solution. Finally, values of the system parameters are suggested to make the satellite behave as desired. Among these is a possible approach to subdue the solar induced roll resonance. It is felt that the approximate analysis presented here should significantly reduce the computational efforts involved in the design and stability analysis of the systems.

Solitons of Coupled Scalar Field Theories in Two Dimensions

This work offers a method for finding some exact soliton solutions to coupled relativistic scalar field theories in 1+1 dimensions. The method can yield static solutions as well as quasistatic "charged" solutions for a variety of Lagrangians. Explicit solutions are derived as examples. A particularly interesting class of solutions is nontopological without being either charged or time dependent.

Detection and determination of manganese concentration in water using a fiber Bragg grating coupled with nanotechnology

Through this paper we experimentally demonstrate the fabrication of a fiber Bragg grating (FBG) chemical sensor to detect and determine the manganese concentration in water and compare our results with sophisticated spectroscopic methods, such as atomic absorption spectrometry and the inductively coupled plasma method. Here we propose a simple method to develop a thin layer of gold nanoparticles above the etched grating region to enhance the sensitivity of the reflected spectrum of the FBG. By doing so, we achieve a sensitivity of 1.26 nm/parts per million in determining the trace level of Mn in water. Proper reagents are used to detect manganese in water. (C) 2011 Optical Society of America

A coupled VOF-IBM-enthalpy approach for modeling motion and growth of equiaxed dendrites in a solidifying melt

A coupled methodology for simulating the simultaneous growth and motion of equiaxed dendrites in solidifying melts is presented. The model uses the volume-averaging principles and combines the features of the enthalpy method for modeling growth, immersed boundary method for handling the rigid solid-liquid interfaces, and the volume of fluid method for tracking the advection of the dendrite. The algorithm also performs explicit-implicit coupling between the techniques used. A two-dimensional framework with incompressible and Newtonian fluid is considered. Validation with available literature is performed and dendrite growth in the presence of rotational and buoyancy driven flow fields is studied. It is seen that the flow fields significantly alter the position and morphology of the dendrites. (C) 2012 Elsevier Inc. All rights reserved.

Arbitrary Lagrangian-Eulerian finite-element method for computation of two-phase flows with soluble surfactants

A finite-element scheme based on a coupled arbitrary Lagrangian-Eulerian and Lagrangian approach is developed for the computation of interface flows with soluble surfactants. The numerical scheme is designed to solve the time-dependent Navier-Stokes equations and an evolution equation for the surfactant concentration in the bulk phase, and simultaneously, an evolution equation for the surfactant concentration on the interface. Second-order isoparametric finite elements on moving meshes and second-order isoparametric surface finite elements are used to solve these equations. The interface-resolved moving meshes allow the accurate incorporation of surface forces, Marangoni forces and jumps in the material parameters. The lower-dimensional finite-element meshes for solving the surface evolution equation are part of the interface-resolved moving meshes. The numerical scheme is validated for problems with known analytical solutions. A number of computations to study the influence of the surfactants in 3D-axisymmetric rising bubbles have been performed. The proposed scheme shows excellent conservation of fluid mass and of the total mass of the surfactant. (C) 2012 Elsevier Inc. All rights reserved.

Selective excitation and detection of maximum quantum coherence of a group of scalar coupled protons in chiral molecules: an NMR experiment for enantiodiscrimination

The H-1 NMR spectroscopic discrimination of enantiomers in the solution state and the measurement of enantiomeric composition is most often hindered due to either very small chemical shift differences between the discriminated peaks or severe overlap of transitions from other chemically non-equivalent protons. In addition the use of chiral auxiliaries such as, crown ether and chiral lanthanide shift reagent may often cause enormous line broadening or give little degree of discrimination beyond the crown ether substrate ratio, hampering the discrimination. In circumventing such problems we are proposing the utilization of the difference in the additive values of all the chemical shifts of a scalar coupled spin system. The excitation and detection of appropriate highest quantum coherence yields the measurable difference in the frequencies between two transitions, one pertaining to each enantiomer in the maximum quantum dimension permitting their discrimination and the F-2 cross section at each of these frequencies yields an enantiopure spectrum. The advantage of the utility of the proposed method is demonstrated on several chiral compounds where the conventional one dimensional H-1 NMR spectra fail to differentiate the enantiomers.

A Finite Element Method Based Approach for Evaluating the Sensitivity of Faraday-Type Electromagnetic Flow Meter for Liquid Sodium Flow

Faraday-type electromagnetic flow meters are employed for measuring the flow rate of liquid sodium in fast breeder reactors. The calibration of such flow meters, owing to the required elaborative arrangements is rather difficult. On the other hand, theoretical approach requires solution of two coupled electromagnetic partial differential equation with profile of the flow and applied magnetic field as the inputs. This is also quite involved due to the 3D nature of the problem. Alternatively, Galerkin finite element method based numerical solution is suggested in the literature as an attractive option for the required calibration. Based on the same, a computer code in Matlab platform has been developed in this work with both 20 and 27 node brick elements. The boundary conditions are correctly defined and several intermediate validation exercises are carried out. Finally it is shown that the sensitivities predicted by the code for flow meters of four different dimensions agrees well with the results given by analytical expression, thereby providing strong validation. Sensitivity for higher flow rates, for which analytical approach does not exist, is shown to decrease with increase in flow velocity.