A model of coupled thermal, mechanical, and electrostatic field effects in III-N thin film heterostructures

A one-dimensional coupled multi-physics based model has been developed to accurately compute the effects of electrostatic, mechanical, and thermal field interactions on the electronic energy band structure in group III-nitrides thin film heterostructures. Earlier models reported in published literature assumes electro-mechanical field with uniform temperature thus neglecting self-heating. Also, the effects of diffused interface on the energy band structure were not studied. We include these effects in a self-consistent manner wherein the transport equation is introduced along with the electro-mechanical models, and the lattice structural variation as observed in experiments are introduced at the interface. Due to these effects, the electrostatic potential distribution in the heterostructure is altered. The electron and hole ground state energies decrease by 5% and 9%, respectively, at a relative temperature of 700 K, when compared with the results obtained from the previously reported electro-mechanical model assuming constant and uniform temperature distribution. A diffused interface decreases the ground state energy of electrons and holes by about 11% and 9%, respectively, at a relative temperature of 700 K when compared with the predictions based on uniform temperature based electro-mechanical model. (C) 2013 AIP Publishing LLC.

An equivalent undercooling model to account for flow effect on binary alloy dendrite growth rate

A theoretical analysis is carried out to observe the influence of important flow parameters such as Nusselt number and Sherwood number on the tip speed of an equiaxed dendrite growing in a convecting alloy melt. The effect of thermal and solutal transfer at the interface due to convection is equated to an undercooling of the melt, and an expression is derived for this equivalent undercooling in terms of the flow Nusselt number and Sherwood number. Results for the equivalent undercooling are compared with corresponding numerical values obtained by performing simulations based on the enthalpy method. This method represents a relatively simple procedure to analyze the effects of melt convection on the growth rate of dendrites. (C) 2013 Elsevier Ltd. All rights reserved.

An enthalpy-based model of dendritic growth in a convecting binary alloy melt

Purpose-In the present work, a numerical method, based on the well established enthalpy technique, is developed to simulate the growth of binary alloy equiaxed dendrites in presence of melt convection. The paper aims to discuss these issues. Design/methodology/approach-The principle of volume-averaging is used to formulate the governing equations (mass, momentum, energy and species conservation) which are solved using a coupled explicit-implicit method. The velocity and pressure fields are obtained using a fully implicit finite volume approach whereas the energy and species conservation equations are solved explicitly to obtain the enthalpy and solute concentration fields. As a model problem, simulation of the growth of a single crystal in a two-dimensional cavity filled with an undercooled melt is performed. Findings-Comparison of the simulation results with available solutions obtained using level set method and the phase field method shows good agreement. The effects of melt flow on dendrite growth rate and solute distribution along the solid-liquid interface are studied. A faster growth rate of the upstream dendrite arm in case of binary alloys is observed, which can be attributed to the enhanced heat transfer due to convection as well as lower solute pile-up at the solid-liquid interface. Subsequently, the influence of thermal and solutal Peclet number and undercooling on the dendrite tip velocity is investigated. Originality/value-As the present enthalpy based microscopic solidification model with melt convection is based on a framework similar to popularly used enthalpy models at the macroscopic scale, it lays the foundation to develop effective multiscale solidification.

Early osseointegration of a strontium containing glass ceramic in a rabbit model

The most important property of a bone cement or a bone substitute in load bearing orthopaedic implants is good integration with host bone with reduced bone resorption and increased bone regeneration at the implant interface. Long term implantation of metal-based joint replacements often results in corrosion and particle release, initiating chronic inflammation leading onto osteoporosis of host bone. An alternative solution is the coating of metal implants with hydroxyapatite (HA) or bioglass or the use of bulk bioglass or HA-based composites. In the above perspective, the present study reports the in vivo biocompatibility and bone healing of the strontium (Sr)-stabilized bulk glass ceramics with the nominal composition of 4.5SiO(2)-3Al(2)O(3)-1.5P(2)O(5)-3SrO-2SrF(2) during short term implantation of up to 12 weeks in rabbit animal model. The progression of healing and bone regeneration was qualitatively and quantitatively assessed using fluorescence microscopy, histological analysis and micro-computed tomography. The overall assessment of the present study establishes that the investigated glass ceramic is biocompatible in vivo with regards to local effects after short term implantation in rabbit animal model. Excellent healing was observed, which is comparable to that seen in response to a commercially available implant of HA-based bioglass alone. (C) 2013 Elsevier Ltd. All rights reserved.

An early in-situ stress signature of the AlN-Si pre-growth interface for successful integration of nitrides with (111) Si

The integration of Metal Organic Chemical Vapor Deposition (MOCVD) grown group III-A nitride device stacks on Si (111) substrates is critically dependent on the quality of the first AlN buffer layer grown. A Si surface that is both oxide-free and smooth is a primary requirement for nucleating such layers. A single parameter, the AlN layer growth stress, is shown to be an early (within 50 nm), clear (<0.5 GPa versus > 1GPa), and fail-safe indicator of the pre-growth surface, and the AlN quality required for successful epitaxy. Grain coalescence model for stress generation is used to correlate growth stress, the AlN-Si interface, and crystal quality. (C) 2013 AIP Publishing LLC.

Length-scale dependent mechanical properties of Al-Cu eutectic alloy: Molecular dynamics based model and its experimental verification

This paper attempts to gain an understanding of the effect of lamellar length scale on the mechanical properties of two-phase metal-intermetallic eutectic structure. We first develop a molecular dynamics model for the in-situ grown eutectic interface followed by a model of deformation of Al-Al2Cu lamellar eutectic. Leveraging the insights obtained from the simulation on the behaviour of dislocations at different length scales of the eutectic, we present and explain the experimental results on Al-Al2Cu eutectic with various different lamellar spacing. The physics behind the mechanism is further quantified with help of atomic level energy model for different length scale as well as different strain. An atomic level energy partitioning of the lamellae and the interface regions reveals that the energy of the lamellae core are accumulated more due to dislocations irrespective of the length-scale. Whereas the energy of the interface is accumulated more due to dislocations when the length-scale is smaller, but the trend is reversed when the length-scale is large beyond a critical size of about 80 nm. (C) 2014 Author(s).

In vitro osteogenic cell proliferation, mineralization, and in vivo osseointegration of injection molded high-density polyethylene-based hybrid composites in rabbit animal model

The present work reports the biocompatibility property of injection molded HDPE-HA-Al2O3 hybrid composites. In vitro cytocompatibility results reveal that osteogenic cell viability and bone mineralization are favorably supported in a statistically significant manner on HDPE-20% HA-20% Al2O3 composite, in comparison to HDPE-40 wt.% HA or HDPE-40 wt.% Al2O3. The difference in cytocompatibility property is explained in terms of difference in substrate wettability/surface energy and importantly, both the cell proliferation at 7 days or bone mineralization at 21 days on HDPE-20% HA-20% Al2O3 composite are either comparable or better than sintered HA. The progressive healing of cylindrical femoral bone defects in rabbit animal model was assessed by implantation experiments over 1, 4 and 12 weeks. Based on the histological analysis as well as histomorphometrical evaluation, a better efficacy of HDPE-20% HA-20% Al2O3 over high-density polyethylene (HDPE) for bone regeneration and neobone formation at host bone-implant interface was established. Taken together, the present study unequivocally establishes that despite the presence of 20% Al2O3, HDPE-based hybrid composites are as biocompatible as HA in vitro or better than HDPE in vivo.

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.

Nonlinear preconditioning for efficient and accurate interface capturing in simulation of multicomponent compressible flows

Single fluid schemes that rely on an interface function for phase identification in multicomponent compressible flows are widely used to study hydrodynamic flow phenomena in several diverse applications. Simulations based on standard numerical implementation of these schemes suffer from an artificial increase in the width of the interface function owing to the numerical dissipation introduced by an upwind discretization of the governing equations. In addition, monotonicity requirements which ensure that the sharp interface function remains bounded at all times necessitate use of low-order accurate discretization strategies. This results in a significant reduction in accuracy along with a loss of intricate flow features. In this paper we develop a nonlinear transformation based interface capturing method which achieves superior accuracy without compromising the simplicity, computational efficiency and robustness of the original flow solver. A nonlinear map from the signed distance function to the sigmoid type interface function is used to effectively couple a standard single fluid shock and interface capturing scheme with a high-order accurate constrained level set reinitialization method in a way that allows for oscillation-free transport of the sharp material interface. Imposition of a maximum principle, which ensures that the multidimensional preconditioned interface capturing method does not produce new maxima or minima even in the extreme events of interface merger or breakup, allows for an explicit determination of the interface thickness in terms of the grid spacing. A narrow band method is formulated in order to localize computations pertinent to the preconditioned interface capturing method. Numerical tests in one dimension reveal a significant improvement in accuracy and convergence; in stark contrast to the conventional scheme, the proposed method retains its accuracy and convergence characteristics in a shifted reference frame. Results from the test cases in two dimensions show that the nonlinear transformation based interface capturing method outperforms both the conventional method and an interface capturing method without nonlinear transformation in resolving intricate flow features such as sheet jetting in the shock-induced cavity collapse. The ability of the proposed method in accounting for the gravitational and surface tension forces besides compressibility is demonstrated through a model fully three-dimensional problem concerning droplet splash and formation of a crownlike feature. (C) 2014 Elsevier Inc. All rights reserved.

Role of interface curvature on stress distribution under indentation for ZrN/Zr multilayer coating

Contact damage in curved interface nano-layeredmetal/nitride (150 (ZrN)/10 (Zr) nm) multilayer is investigated in order to understand the role of interface morphology on contact damage under indentation. A finite element method (FEM) model was formulated with different wavelengths of 1000 nm, 500 nm, 250 nm and common height of 50 nm, which gives insight on the effect of different curvature on stress field generated under indentation. Elastic-plastic properties were assigned to the metal layer and substrate while the nitride layer was assigned perfectly elastic properties. Curved interface multilayers show delamination along the metal/nitride interface and vertical cracks emanating from the ends of the delamination. FEM revealed the presence of tensile stress normal to the interface even under the contact, along with tensile radial stresses, both present at the valley part of the curve, which leads to vertical cracks associated with interfacial delamination. Stress enhancement was seen to be relatively insensitive to curvature. (C) 2014 Elsevier B.V. All rights reserved.

Colored polydimethylsiloxane micropillar arrays for high throughput measurements of forces applied by genetic model organisms

Measuring forces applied by multi-cellular organisms is valuable in investigating biomechanics of their locomotion. Several technologies have been developed to measure such forces, for example, strain gauges, micro-machined sensors, and calibrated cantilevers. We introduce an innovative combination of techniques as a high throughput screening tool to assess forces applied by multiple genetic model organisms. First, we fabricated colored Polydimethylsiloxane (PDMS) micropillars where the color enhances contrast making it easier to detect and track pillar displacement driven by the organism. Second, we developed a semiautomated graphical user interface to analyze the images for pillar displacement, thus reducing the analysis time for each animal to minutes. The addition of color reduced the Young's modulus of PDMS. Therefore, the dye-PDMS composite was characterized using Yeoh's hyperelastic model and the pillars were calibrated using a silicon based force sensor. We used our device to measure forces exerted by wild type and mutant Caenorhabditis elegans moving on an agarose surface. Wild type C. elegans exert an average force of similar to 1 mu N on an individual pillar and a total average force of similar to 7.68 mu N. We show that the middle of C. elegans exerts more force than its extremities. We find that C. elegans mutants with defective body wall muscles apply significantly lower force on individual pillars, while mutants defective in sensing externally applied mechanical forces still apply the same average force per pillar compared to wild type animals. Average forces applied per pillar are independent of the length, diameter, or cuticle stiffness of the animal. We also used the device to measure, for the first time, forces applied by Drosophila melanogaster larvae. Peristaltic waves occurred at 0.4Hz applying an average force of similar to 1.58 mu N on a single pillar. Our colored microfluidic device along with its displacement tracking software allows us to measure forces applied by multiple model organisms that crawl or slither to travel through their environment. (C) 2015 AIP Publishing LLC.

Effects of Interface Treatment on the Fatigue Behavior of Shape Memory Alloy Reinforced Polymer Composites

Interfacial properties of Shape Memory Alloy (SMA) reinforced polymer matrix composites can be enhanced by improving the interfacial bonding. This paper focuses on studying the interfacial stresses developed in the SMA-epoxy interface due to various laser shot penning conditions. Fiber-pull test-setup is designed to understand the role of mechanical bias stress cycling and thermal actuation cycling. Phase transformation is tracked over mechanical and thermal fatigue cycles. A micromechanics based model developed earlier based on shear lag in SMA and energy based consistent homogenization is extended here to incorporate the stress-temperature phase diagram parameters for modeling fatigue.

Effects of Patterned Substrate on Thin Films Simulated by Family Model

Patterned substrate growth has been a subject of much interest. In this work, characteristics of some statistical properties of a film grown on triangular and vicinal substrates using the Family model are studied. Substrate size and tilt angle are varied. It is found that the interface width and the correlation function increase as the roughness of the pattern is increased. The new scaling exponents are calculated and anomalous scaling is obtained. The transient persistence probability does not show a power law relation when the initial surface is sufficiently rough. The initial rough surface also causes multifractal behavior in the model.

A micropolar cohesive damage model for delamination of composites

A micropolar cohesive damage model for delamination of composites is proposed. The main idea is to embed micropolarity, which brings an additional layer of kinematics through the micro-rotation degrees of freedom within a continuum model to account for the micro-structural effects during delamination. The resulting cohesive model, describing the modified traction separation law, includes micro-rotational jumps in addition to displacement jumps across the interface. The incorporation of micro-rotation requires the model to be supplemented with physically relevant material length scale parameters, whose effects during delamination of modes I and II are brought forth using numerical simulations appropriately supported by experimental evidences. (C) 2015 Elsevier Ltd. All rights reserved.