Finite element simulations of notch tip fields in magnesium single crystals

Recent experiments using three point bend specimens of Mg single crystals have revealed that tensile twins of {10 (1) over bar2}-type form profusely near a notch tip and enhance the fracture toughness through large plastic dissipation. In this work, 3D finite element simulations of these experiments are carried out using a crystal plasticity framework which includes slip and twinning to gain insights on the mechanics of fracture. The predicted load-displacement curves, slip and tensile twinning activities from finite element analysis corroborate well with the experimental observations. The numerical results are used to explore the 3D nature of the crack tip stress, plastic slip and twin volume fraction distributions near the notch root. The occurrence of tensile twinning is rationalized from the variation of normal stress ahead of the notch tip. Further, deflection of the crack path at twin-twin intersections observed in the experiments is examined from an energy standpoint by modeling discrete twins close to the notch root.

Large deformation dynamic finite element analysis of delaminated composite plates using contact-impact conditions

A new C-0 composite plate finite element based on Reddy's third order theory is used for large deformation dynamic analysis of delaminated composite plates. The inter-laminar contact is modeled with an augmented Lagrangian approach. Numerical results show that the widely used ``unconditionally stable'' beta-Newmark method presents instability problems in the transient simulation of delaminated composite plate structures with large deformation. To overcome this instability issue, an energy and momentum conserving composite implicit time integration scheme presented by Bathe and Baig is used. It is found that a proper selection of the penalty parameter is very crucial in the contact simulation. (C) 2014 Elsevier Ltd. All rights reserved.

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.

Wave propagation analysis in adhesively bonded composite joints using the wavelet spectral finite element method

A wavelet spectral finite element (WSFE) model is developed for studying transient dynamics and wave propagation in adhesively bonded composite joints. The adherands are formulated as shear deformable beams using the first order shear deformation theory (FSDT) to obtain accurate results for high frequency wave propagation. Equations of motion governing wave motion in the bonded beams are derived using Hamilton's principle. The adhesive layer is modeled as a line of continuously distributed tension/compression and shear springs. Daubechies compactly supported wavelet scaling functions are used to transform the governing partial differential equations from time domain to frequency domain. The dynamic stiffness matrix is derived under the spectral finite element framework relating the nodal forces and displacements in the transformed frequency domain. Time domain results for wave propagation in a lap joint are validated with conventional finite element simulations using Abaqus. Frequency domain spectrum and dispersion relation results are presented and discussed. The developed WSFE model yields efficient and accurate analysis of wave propagation in adhesively-bonded composite joints. (C) 2014 Elsevier Ltd. All rights reserved.

A Smooth Discretization Bridging Finite Element and Mesh-free Methods Using Polynomial Reproducing Simplex Splines

This work sets forth a `hybrid' discretization scheme utilizing bivariate simplex splines as kernels in a polynomial reproducing scheme constructed over a conventional Finite Element Method (FEM)-like domain discretization based on Delaunay triangulation. Careful construction of the simplex spline knotset ensures the success of the polynomial reproduction procedure at all points in the domain of interest, a significant advancement over its precursor, the DMS-FEM. The shape functions in the proposed method inherit the global continuity (Cp-1) and local supports of the simplex splines of degree p. In the proposed scheme, the triangles comprising the domain discretization also serve as background cells for numerical integration which here are near-aligned to the supports of the shape functions (and their intersections), thus considerably ameliorating an oft-cited source of inaccuracy in the numerical integration of mesh-free (MF) schemes. Numerical experiments show the proposed method requires lower order quadrature rules for accurate evaluation of integrals in the Galerkin weak form. Numerical demonstrations of optimal convergence rates for a few test cases are given and the method is also implemented to compute crack-tip fields in a gradient-enhanced elasticity model.

A spectral element for wave propagation in honeycomb sandwich construction considering core flexibility

Spectral elements are found to be extremely resourceful to study the wave propagation characteristics of structures at high frequencies. Most of the aerospace structures use honeycomb sandwich constructions. The existing spectral elements use single layer theories for a sandwich construction wherein the two face sheets vibrate together and this model is sufficient for low frequency excitations. At high frequencies, the two face sheets vibrate independently. The Extended Higher order SAndwich Plate theory (EHSaPT) is suitable for representing the independent motion of the face sheets. A 1D spectral element based on EHSaPT is developed in this work. The wave number and the wave speed characteristics are obtained using the developed spectral element. It is shown that the developed spectral element is capable of representing independent wave motions of the face sheets. The propagation speeds of a high frequency modulated pulse in the face sheets and the core of a honeycomb sandwich are demonstrated. Responses of a typical honeycomb sandwich beam to high frequency shock loads are obtained using the developed spectral element and the response match very well with the finite element results. It is shown that the developed spectral element is able to represent the flexibility of the core resulting into independent wave motions in the face sheets, for which a finite element method needs huge degrees of freedom. (C) 2015 Elsevier Ltd. All rights reserved.

Multi-transform based spectral element to include first order shear deformation in plates

For obtaining dynamic response of structure to high frequency shock excitation spectral elements have several advantages over conventional methods. At higher frequencies transverse shear and rotary inertia have a predominant role. These are represented by the First order Shear Deformation Theory (FSDT). But not much work is reported on spectral elements with FSDT. This work presents a new spectral element based on the FSDT/Mindlin Plate Theory which is essential for wave propagation analysis of sandwich plates. Multi-transformation method is used to solve the coupled partial differential equations, i.e., Laplace transforms for temporal approximation and wavelet transforms for spatial approximation. The formulation takes into account the axial-flexure and shear coupling. The ability of the element to represent different modes of wave motion is demonstrated. Impact on the derived wave motion characteristics in the absence of the developed spectral element is discussed. The transient response using the formulated element is validated by the results obtained using Finite Element Method (FEM) which needs significant computational effort. Experimental results are provided which confirms the need to having the developed spectral element for the high frequency response of structures. (C) 2015 Elsevier Ltd. All rights reserved.

Environmental monitoring of bacterial contamination and antibiotic resistance patterns of the fecal coliforms isolated from Cauvery River, a major drinking water source in Karnataka, India

The present study focuses prudent elucidation of microbial pollution and antibiotic sensitivity profiling of the fecal coliforms isolated from River Cauvery, a major drinking water source in Karnataka, India. Water samples were collected from ten hotspots during the year 2011-2012. The physiochemical characteristics and microbial count of water samples collected from most of the hotspots exhibited greater biological oxygen demand and bacterial count especially coliforms in comparison with control samples (p <= 0.01). The antibiotic sensitivity testing was performed using 48 antibiotics against the bacterial isolates by disk-diffusion assay. The current study showed that out of 848 bacterial isolates, 93.51 % (n=793) of the isolates were found to be multidrug-resistant to most of the current generation antibiotics. Among the major isolates, 96.46 % (n=273) of the isolates were found to be multidrug-resistant to 30 antibiotics and they were identified to be Escherichia coli by 16S rDNA gene sequencing. Similarly, 93.85 % (n=107), 94.49 % (n=103), and 90.22 % (n=157) of the isolates exhibited multiple drug resistance to 32, 40, and 37 antibiotics, and they were identified to be Enterobacter cloacae, Pseudomonas trivialis, and Shigella sonnei, respectively. The molecular studies suggested the prevalence of blaTEM genes in all the four isolates and dhfr gene in Escherichia coli and Sh. sonnei. Analogously, most of the other Gram-negative bacteria were found to be multidrug-resistant and the Gram-positive bacteria, Staphylococcus spp. isolated from the water samples were found to be methicillin and vancomycin-resistant Staphylococcus aureus. This is probably the first study elucidating the bacterial pollution and antibiotic sensitivity profiling of fecal coliforms isolated from River Cauvery, Karnataka, India.

Use of permeable reactive barrier to mitigate groundwater nitrate contamination from on-site sanitation

Nearly 50% of India's population depends on variants of pit-toilet systems for human waste disposal. Nitrate contamination of groundwater by pit-toilet leachate is a major environmental concern in the country as it sources a major proportion (50-80%) of potable water from aquifers. Therefore, minimizing nitrate contamination of groundwater due to leachate infiltration from pit-toilet systems is essential. Batch and column experiments demonstrated the capability of bentonite-enhanced sand (BES) specimens to reduce nitrate concentrations in synthetic solutions (initial NO3-N concentration = 22.7 mg/L, C/N = 3) by about 85-90% in 10 to 24 hour by a heterotrophic denitrification process. Based on the laboratory results, it is recommended that use of a BES-permeable reactive barrier layer at the base of pit-toilets will facilitate heterotrophic denitrification and mitigate nitrate contamination of the underlying aquifer.

A Reliable Residual Based A Posteriori Error Estimator for a Quadratic Finite Element Method for the Elliptic Obstacle Problem

A residual based a posteriori error estimator is derived for a quadratic finite element method (FEM) for the elliptic obstacle problem. The error estimator involves various residuals consisting of the data of the problem, discrete solution and a Lagrange multiplier related to the obstacle constraint. The choice of the discrete Lagrange multiplier yields an error estimator that is comparable with the error estimator in the case of linear FEM. Further, an a priori error estimate is derived to show that the discrete Lagrange multiplier converges at the same rate as that of the discrete solution of the obstacle problem. The numerical experiments of adaptive FEM show optimal order convergence. This demonstrates that the quadratic FEM for obstacle problem exhibits optimal performance.

Effect of a mesh on boundary layer transitions induced by free-stream turbulence and an isolated roughness element

Streamwise streaks, their lift-up and streak instability are integral to the bypass transition process. An experimental study has been carried out to find the effect of a mesh placed normal to the flow and at different wall-normal locations in the late stages of two transitional flows induced by free-stream turbulence (FST) and an isolated roughness element. The mesh causes an approximately 30% reduction in the free-stream velocity, and mild acceleration, irrespective of its wall-normal location. Interestingly, when located near the wall, the mesh suppresses several transitional events leading to transition delay over a large downstream distance. The transition delay is found to be mainly caused by suppression of the lift-up of the high-shear layer and its distortion, along with modification of the spanwise streaky structure to an orderly one. However, with the mesh well away from the wall, the lifted-up shear layer remains largely unaffected, and the downstream boundary layer velocity profile develops an overshoot which is found to follow a plane mixing layer type profile up to the free stream. Reynolds stresses, and the size and strength of vortices increase in this mixing layer region. This high-intensity disturbance can possibly enhance transition of the accelerated flow far downstream, although a reduction in streamwise turbulence intensity occurs over a short distance downstream of the mesh. However, the shape of the large-scale streamwise structure in the wall-normal plane is found to be more or less the same as that without the mesh.

Multiobjective fuzzy optimization for sustainable groundwater management using particle swarm optimization and analytic element method

Groundwater management involves conflicting objectives as maximization of discharge contradicts the criteria of minimum pumping cost and minimum piping cost. In addition, available data contains uncertainties such as market fluctuations, variations in water levels of wells and variations of ground water policies. A fuzzy model is to be evolved to tackle the uncertainties, and a multiobjective optimization is to be conducted to simultaneously satisfy the contradicting objectives. Towards this end, a multiobjective fuzzy optimization model is evolved. To get at the upper and lower bounds of the individual objectives, particle Swarm optimization (PSO) is adopted. The analytic element method (AEM) is employed to obtain the operating potentio metric head. In this study, a multiobjective fuzzy optimization model considering three conflicting objectives is developed using PSO and AEM methods for obtaining a sustainable groundwater management policy. The developed model is applied to a case study, and it is demonstrated that the compromise solution satisfies all the objectives with adequate levels of satisfaction. Sensitivity analysis is carried out by varying the parameters, and it is shown that the effect of any such variation is quite significant. Copyright (c) 2015 John Wiley & Sons, Ltd.

Atmospheric pressure plasma chemical vapor deposition reactor for 100 mm wafers, optimized for minimum contamination at low gas flow rates

Gas discharge plasmas used for thinfilm deposition by plasma-enhanced chemical vapor deposition (PECVD) must be devoid of contaminants, like dust or active species which disturb the intended chemical reaction. In atmospheric pressure plasma systems employing an inert gas, the main source of such contamination is the residual air inside the system. To enable the construction of an atmospheric pressure plasma (APP) system with minimal contamination, we have carried out fluid dynamic simulation of the APP chamber into which an inert gas is injected at different mass flow rates. On the basis of the simulation results, we have designed and built a simple, scaled APP system, which is capable of holding a 100 mm substrate wafer, so that the presence of air (contamination) in the APP chamber is minimized with as low a flow rate of argon as possible. This is examined systematically by examining optical emission from the plasma as a function of inert gas flow rate. It is found that optical emission from the plasma shows the presence of atmospheric air, if the inlet argon flow rate is lowered below 300 sccm. That there is minimal contamination of the APP reactor built here, was verified by conducting an atmospheric pressure PECVD process under acetylene flow, combined with argon flow at 100 sccm and 500 sccm. The deposition of a polymer coating is confirmed by infrared spectroscopy. X-ray photoelectron spectroscopy shows that the polymer coating contains only 5% of oxygen, which is comparable to the oxygen content in polymer deposits obtained in low-pressure PECVD systems. (C) 2015 AIP Publishing LLC.

Bearing capacity of circular footings over rock mass by using axisymmetric quasi lower bound finite element limit analysis

The ultimate bearing capacity of a circular footing, placed over rock mass, is evaluated by using the lower bound theorem of the limit analysis in conjunction with finite elements and nonlinear optimization. The generalized Hoek-Brown (HB) failure criterion, but by keeping a constant value of the exponent, alpha = 0.5, was used. The failure criterion was smoothened both in the meridian and pi planes. The nonlinear optimization was carried out by employing an interior point method based on the logarithmic barrier function. The results for the obtained bearing capacity were presented in a non-dimensional form for different values of GSI, m(i), sigma(ci)/(gamma b) and q/sigma(ci). Failure patterns were also examined for a few cases. For validating the results, computations were also performed for a strip footing as well. The results obtained from the analysis compare well with the data reported in literature. Since the equilibrium conditions are precisely satisfied only at the centroids of the elements, not everywhere in the domain, the obtained lower bound solution will be approximate not true. (C) 2015 Elsevier Ltd. All rights reserved.

A FRAMEWORK FOR THE ERROR ANALYSIS OF DISCONTINUOUS FINITE ELEMENT METHODS FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS AND APPLICATIONS TO C-0 IP METHODS

In this article, an abstract framework for the error analysis of discontinuous Galerkin methods for control constrained optimal control problems is developed. The analysis establishes the best approximation result from a priori analysis point of view and delivers a reliable and efficient a posteriori error estimator. The results are applicable to a variety of problems just under the minimal regularity possessed by the well-posedness of the problem. Subsequently, the applications of C-0 interior penalty methods for a boundary control problem as well as a distributed control problem governed by the biharmonic equation subject to simply supported boundary conditions are discussed through the abstract analysis. Numerical experiments illustrate the theoretical findings.