Mixed finite elements for electromagnetic analysis

The occurrence of spurious solutions is a well-known limitation of the standard nodal finite element method when applied to electromagnetic problems. The two commonly used remedies that are used to address this problem are (i) The addition of a penalty term with the penalty factor based on the local dielectric constant, and which reduces to a Helmholtz form on homogeneous domains (regularized formulation); (ii) A formulation based on a vector and a scalar potential. Both these strategies have some shortcomings. The penalty method does not completely get rid of the spurious modes, and both methods are incapable of predicting singular eigenvalues in non-convex domains. Some non-zero spurious eigenvalues are also predicted by these methods on non-convex domains. In this work, we develop mixed finite element formulations which predict the eigenfrequencies (including their multiplicities) accurately, even for nonconvex domains. The main feature of the proposed mixed finite element formulation is that no ad-hoc terms are added to the formulation as in the penalty formulation, and the improvement is achieved purely by an appropriate choice of finite element spaces for the different variables. We show that the formulation works even for inhomogeneous domains where `double noding' is used to enforce the appropriate continuity requirements at an interface. For two-dimensional problems, the shape of the domain can be arbitrary, while for the three-dimensional ones, with our current formulation, only regular domains (which can be nonconvex) can be modeled. Since eigenfrequencies are modeled accurately, these elements also yield accurate results for driven problems. (C) 2014 Elsevier Ltd. All rights reserved.

Effect of extrusion ratio on the microstructure, texture and mechanical properties of (Mg/AZ91)(m)-SiCp composite

Cast Mg/SiCp and AZ91/SiCp composites were successfully hot extruded vis-a-vis cast and unreinforced Mg and AZ91 alloy up to low (R=15:1) and high (R=54:1) extrusion ratios at 350 degrees C. Significant matrix grain refinement was noticed after extrusion due to dynamic recrystallization; the degree of refinement being relatively higher for the two composites. The AZ91 based materials (AZ91 and AZ91/SiCp) exhibited comparatively finer grain size both in cast condition and after extrusion due to strong pinning effect from alloying elements as well as Mg17Al12 intermetallic phase. Compositional analyses eliminated the possibility of any interfacial reaction between matrix (Mg/AZ91) and second phase reinforcement (SiCp) in case of the composites. Texture evolution shows the formation of < 10 (1) over bar0 >parallel to ED texture fibre for all the materials after extrusion irrespective of SiCp addition or alloying which is primarily due to the deformation of the matrix phase. Micro-hardness did not significantly increased on extrusion in comparison to the respective cast materials for both composites and unreinforced alloys. Dynamic mechanical analysis, however, confirmed that the damping properties were affected by the extrusion ratio and to a lesser extent, due to the presence of second phase at room temperature as well as at higher temperature (300 degrees C). (C) 2014 Elsevier B.V. All rights reserved.

Hybrid Elements for Modelling Squeeze Film Effects Coupled with Structural Interactions in Vibratory MEMS Devices

We present a hybrid finite element based methodology to solve the coupled fluid structure problem of squeeze film effects in vibratory MEMS devices, such as gyroscopes, RF switches, and 2D resonators. The aforementioned devices often have a thin plate like structure vibrating normally to a fixed substrate, and are generally not perfectly vacuum packed. This results in a thin air film being trapped between the vibrating plate and the fixed substrate which behaves like a squeeze film offering both stiffness and damping. For accurate modelling of such devices the squeeze film effects must be incorporated. Extensive literature is available on squeeze film modelling, however only a few studies address the coupled fluid elasticity problem. The majority of the studies that account for the plate elasticity coupled with the fluid equation, either use approximate mode shapes for the plate or use iterative solution strategies. In an earlier work we presented a single step coupled methodology using only one type of displacement based element to solve the coupled problem. The displacement based finite element models suffer from locking issues when it comes to modelling very thin structures with the lateral dimensions much larger than the plate thickness as is typical in MEMS devices with squeeze film effects. In this work we present another coupled formulation where we have used hybrid elements to model the structural domain. The numerical results show a huge improvement in convergence and accuracy with coarse hybrid mesh as compared to displacement based formulations. We further compare our numerical results with experimental data from literature and find them to be in good accordance.

Estimation of mechanical properties of single wall carbon nanotubes using molecular mechanics approach

Molecular mechanics based finite element analysis is adopted in the current work to evaluate the mechanical properties of Zigzag, Armchair and Chiral Single wall Carbon Nanotubes (SWCNT) of different diameters and chiralities. Three different types of atomic bonds, that is Carbon Carbon covalent bond and two types of Carbon Carbon van der Waals bonds are considered in the carbon nanotube system. The stiffness values of these bonds are calculated using the molecular potentials, namely Morse potential function and Lennard-Jones interaction potential function respectively and these stiffness's are assigned to spring elements in the finite element model of the CNT. The geometry of CNT is built using a macro that is developed for the finite element analysis software. The finite element model of the CNT is constructed, appropriate boundary conditions are applied and the behavior of mechanical properties of CNT is studied.

Microstructure and mechanical properties at TiCp/Ni-alloy interfaces in laser-synthesized coatings

Titanium carbide particle (TiCp) reinforced Ni alloy composite coatings were synthesized by laser cladding using a cw 3 kW CO2 laser. Two kinds of coatings were present in terms of TiCp origins, i.e. undissolved and in situ reacted TiCp, respectively. The former came from the TiCp pre-coated on the sample, whereas the latter from in situ reaction between titanium and graphite in the molten pool during laser irradiation. Conventional and high-resolution transmission electron microscope observations showed the epitaxial growth of TiC, the precipitation of CrB, and the chemical reaction between Ti and B elements around phase interfaces of undissolved TiCp. The hardness, H, and elastic modulus, E, were measured by nanoindentation of the matrix near the TiCp interface. For undissolved TiCp, the loading curve revealed pop-in phenomena caused by the plastic deformation of the crack formation or debounding of TiCp from the matrix. As for in situ generated TiCp, no pop-in mark appears. On the other hand, in situ reacted TiCp led to much higher hardness and modulus than that in the case of undissolved TiCp. The coating reinforced by in situ generated TiCp displayed the highest impact wear resistance at both low and high impact conditions, as compared with coatings with undissolved TiCp and without TiCp. The impact wear resistance of the coating reinforced by undissolved TiCp increases at a low impact work but decreases at a high impact work, as compared with the single Ni alloy coating. The degree of wear for the composite coating depends primarily on the debonding removal of TiCp.

Determination Ofinterfacial Mechanical Parameters For An Al/Epoxy/Al2O3 System Byusing Peel Test Simulations

Peel test measurements and simulations of the interfacial mechanical parameters for the Al/Epoxy/Al2O3 system are performed in the present investigation. A series of Al film thicknesses between 20 and 250 microns and three peel angles of 90, 135 and 180 degrees are considered. Two types of epoxy adhesives are adopted to obtain both strong and weak interface adhesions. A finite element model with cohesive zone elements is used to identify the interfacial parameters and simulate the peel test process. By simulating and recording normal stress near the crack tip, the separation strength is obtained. Furthermore, the cohesive energy is identified by comparing the simulated steady-state peel force and the experimental result. It is found from the research that both the cohesive energy and the separation strength can be taken as the intrinsic interfacial parameters which are dependent on the thickness of the adhesive layer and independent of the film thickness and peel angle.

I. Quantum mechanics of one-dimensional two-particle models. II. A quantum mechanical treatment of inelastic collisions

Part I

Solutions of Schrödinger’s equation for system of two particles bound in various stationary one-dimensional potential wells and repelling each other with a Coulomb force are obtained by the method of finite differences. The general properties of such systems are worked out in detail for the case of two electrons in an infinite square well. For small well widths (1-10 a.u.) the energy levels lie above those of the noninteresting particle model by as much as a factor of 4, although excitation energies are only half again as great. The analytical form of the solutions is obtained and it is shown that every eigenstate is doubly degenerate due to the “pathological” nature of the one-dimensional Coulomb potential. This degeneracy is verified numerically by the finite-difference method. The properties of the square-well system are compared with those of the free-electron and hard-sphere models; perturbation and variational treatments are also carried out using the hard-sphere Hamiltonian as a zeroth-order approximation. The lowest several finite-difference eigenvalues converge from below with decreasing mesh size to energies below those of the “best” linear variational function consisting of hard-sphere eigenfunctions. The finite-difference solutions in general yield expectation values and matrix elements as accurate as those obtained using the “best” variational function.

The system of two electrons in a parabolic well is also treated by finite differences. In this system it is possible to separate the center-of-mass motion and hence to effect a considerable numerical simplification. It is shown that the pathological one-dimensional Coulomb potential gives rise to doubly degenerate eigenstates for the parabolic well in exactly the same manner as for the infinite square well.

Part II

A general method of treating inelastic collisions quantum mechanically is developed and applied to several one-dimensional models. The formalism is first developed for nonreactive “vibrational” excitations of a bound system by an incident free particle. It is then extended to treat simple exchange reactions of the form A + BC →AB + C. The method consists essentially of finding a set of linearly independent solutions of the Schrödinger equation such that each solution of the set satisfies a distinct, yet arbitrary boundary condition specified in the asymptotic region. These linearly independent solutions are then combined to form a total scattering wavefunction having the correct asymptotic form. The method of finite differences is used to determine the linearly independent functions.

The theory is applied to the impulsive collision of a free particle with a particle bound in (1) an infinite square well and (2) a parabolic well. Calculated transition probabilities agree well with previously obtained values.

Several models for the exchange reaction involving three identical particles are also treated: (1) infinite-square-well potential surface, in which all three particles interact as hard spheres and each two-particle subsystem (i.e. BC and AB) is bound by an attractive infinite-square-well potential; (2) truncated parabolic potential surface, in which the two-particle subsystems are bound by a harmonic oscillator potential which becomes infinite for interparticle separations greater than a certain value; (3) parabolic (untruncated) surface. Although there are no published values with which to compare our reaction probabilities, several independent checks on internal consistency indicate that the results are reliable.

Ultrastructural mechanical and material characterization of fossilized bone

Bone plays a key role in the paleontological and archeological records and can provide insight into the biology, ecology and the environment of ancient vertebrates. Examination of bone at the tissue level reveals a definitive relationship between nanomechanical properties and the local organic content, mineral content, and microstructural organization. However, it is unclear as to how these properties change following fossilization, or diagenesis, where the organic phase is rapidly removed and the remaining mineral phase is reinforced by the deposition of apatites, calcites, and other minerals. While the process of diagenesis is poorly understood, its outcome clearly results in the potential for dramatic alteration of the mechanical response of biological tissues. In this study, fossilized specimens of mammalian long bones, collected from Colorado and Wyoming, were studied for mechanical variations. Nanoindentation performed in both longitudinal and transverse directions revealed preservation of bone's natural anisotropy as transverse modulus values were consistently smaller than longitudinal values. Additionally, modulus values of fossilized bone from 35.0 to 89.1 GPa increased linearly with logarithm of the sample's age. Future studies will aim to clarify what mechanical and material elements of bone are retained during diagenesis as bone becomes part of the geologic milieu. © 2007 Materials Research Society.

An approach to functional synthesis of solutions in mechanical conceptual design. Part III: Spatial configuration

These three papers describe an approach to the synthesis of solutions to a class of mechanical design problems; these involve transmission and transformation of mechanical forces and motion, and can be described by a set of inputs and outputs. The approach involves (1) identifying a set of primary functional elements and rules of combining them, and (2) developing appropriate representations and reasoning procedures for synthesising solution concepts using these elements and their combination rules; these synthesis procedures can produce an exhaustive set of solution concepts, in terms of their topological as well as spatial configurations, to a given design problem. This paper (Part III) describes a constraint propagation procedure which, using a knowledge base of spatial information about a set of primary functional elements, can produce possible spatial configurations of solution concepts generated in Part II.

A coupled hydro-mechanical analysis for prediction of hydraulic fracture propagation in saturated porous media using EFG mesh-less method

The details of the Element Free Galerkin (EFG) method are presented with the method being applied to a study on hydraulic fracturing initiation and propagation process in a saturated porous medium using coupled hydro-mechanical numerical modelling. In this EFG method, interpolation (approximation) is based on nodes without using elements and hence an arbitrary discrete fracture path can be modelled.The numerical approach is based upon solving two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Displacement increment and pore water pressure increment are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system of equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on the penalty method. In order to model discrete fractures, the so-called diffraction method is used.Examples are presented and the results are compared to some closed-form solutions and FEM approximations in order to demonstrate the validity of the developed model and its capabilities. The model is able to take the anisotropy and inhomogeneity of the material into account. The applicability of the model is examined by simulating hydraulic fracture initiation and propagation process from a borehole by injection of fluid. The maximum tensile strength criterion and Mohr-Coulomb shear criterion are used for modelling tensile and shear fracture, respectively. The model successfully simulates the leak-off of fluid from the fracture into the surrounding material. The results indicate the importance of pore fluid pressure in the initiation and propagation pattern of fracture in saturated soils. © 2013 Elsevier Ltd.

Quantum-mechanical calculations of vibrationally resolved cross sections for non-dissociative charge transfer of O3+ with H-2

Charge transfer due to collisions of ground state O3+ (2s(2)2p P-2) ions with molecular hydrogen is investigated using the quantum-mechanical molecular-orbital close-coupling (MOCC) method, and electronic and vibrational state-selective cross sections along with the corresponding differential cross sections are calculated for projectile energies of 100, 500, 1000 and 5000 eV/u at the orientation angles of 25 degrees,45 degrees and 89 degrees. The adiabatic potentials and radial coupling matrix elements utilized in the QMOCC calculations were obtained with the spin-coupled valence-bond approach. The infinite order sudden approximation (IOSA) and the vibrational sudden approximation (VSA) are utilized to deal with the rotation of H-2 and the coupling between the electron and the vibration of H-2. It is found that the distribution of vibrationally resolved cross sections with the vibrational quantum number upsilon' of H-2(+) (upsilon') varies with the increment of the projectile energy; and the electronic and vibrational stateselective differential cross sections show similar behaviors: there is a highest platform within a very small scattering angle, beyond which the differential cross sections decrease as the scattering angle increases and lots of oscillating structures appear, where the scattering angle of the first structure decreases as E-P(-1/2) with the increment of the projectile energy E-P; and the structure and amplitude of the differential cross sections are sensitive to the orientation of molecule H-2, which provides a possibility to identify the orientations of molecule H-2 by the vibrational state-selective differential scattering processes.

Effects of cerium on the microstructure and mechanical properties of Mg–20Zn–8Al alloy

Mg-20Zn-8Al-xCe(x=0-2 wt.%) alloys were prepared by metal mould casting method, the effects of Ce on the microstructure and mechanical properties of the alloys were investigated. The results showed that the dendrite as well as gram size were refined by the addition of Ce, and the best refinement was obtained in 1.39% Ce containing alloy.

Microstructures, mechanical properties and corrosion behavior of high-pressure die-cast Mg-4Al-0.4Mn-xPr (x=1, 2, 4, 6) alloys

Mg-4Al-0.4Mn-xPr (x = 1, 2, 4 and 6 wt.%) magnesium alloys were prepared successfully by the high-pressure die-casting technique. The microstructures, mechanical properties, corrosion behavior as well as strengthening mechanism were investigated. The die-cast alloys were mainly composed of small equiaxed dendrites and the matrix. The fine rigid skin region was related to the high cooling rate and the aggregation of alloying elements, such as Pr. With the Pr content increasing, the alpha-Mg grain sizes were reduced gradually and the amounts of the Al2Pr phase and All, Pr-3 phase which mainly concentrated along the grain boundaries were increased and the relative volume ratio of above two phases was changed. Considering the performance-price ratio, the Pr content added around 4 wt.% was suitable to obtain the optimal mechanical properties which can keep well until 200 degrees C as well as good corrosion resistance. The outstanding mechanical properties were mainly attributed to the rigid casting surface layer, grain refinement, grain boundary strengthening obtained by an amount of precipitates as well as solid solution strengthening.

Effects of cerium on the microstructure and mechanical properties of Mg-20Zn-8Al alloy

Mg-20Zn-8Al-xCe(x=0-2 wt.%) alloys were prepared by metal mould casting method, the effects of Ce on the microstructure and mechanical properties of the alloys were investigated. The results showed that the dendrite as well as gram size were refined by the addition of Ce, and the best refinement was obtained in 1.39% Ce containing alloy. The main phases in the as cast alloys were alpha-Mg and tau-Mg-32 (Al, Zn)(49), and Al4Ce phase was found in the alloys contained more than 1.39% Ce. The addition of Ce improved the mechanical properties of the alloys. The strengthening mechanism was attributed to grain refinement and compound reinforced.

Structure and mechanical properties of Mg-Zn-RE system alloys

The lightest density of Mg has stimulated renewed interest in Mg based alloys for applications in the automotive, aerospace and communications industries. However, Mg in the pure form has relatively low strength, limited ductility and is susceptible to corrosion. Great efforts have been made to improve the mechanical properties of Mg alloys. Alloying Mg with other elements is one of the most important methods. An important class of Mg alloys is the Mg-Zn-RE system (RE = rare earth elements). In recent few decades, a series of new Mg-Zn-RE system alloys have been obtained, and detailed the structure and mechanical properties of the alloys. In this paper, the structure and mechanical properties of the Mg-Zn-RE alloys have been summarized. It showed that these alloys have high strength and they are prospected to be widely used in the future.