Dendrimer building toolkit: Model building and characterization of various dendrimer architectures

We have developed a graphical user interface based dendrimer builder toolkit (DBT) which can be used to generate the dendrimer configuration of desired generation for various dendrimer architectures. The validation of structures generated by this tool was carried out by studying the structural properties of two well known classes of dendrimers: ethylenediamine cored poly(amidoamine) (PAMAM) dendrimer, diaminobutyl cored poly(propylene imine) (PPI) dendrimer. Using full atomistic molecular dynamics (MD) simulation we have calculated the radius of gyration, shape tensor and monomer density distribution for PAMAM and PPI dendrimer at neutral and high pH. A good agreement between the available simulation and experimental (small angle X-ray and neutron scattering; SAXS, SANS) results and calculated radius of gyration was observed. With this validation we have used DBT to build another new class of nitrogen cored poly(propyl ether imine) dendrimer and study it's structural features using all atomistic MD simulation. DBT is a versatile tool and can be easily used to generate other dendrimer structures with different chemistry and topology. The use of general amber force field to describe the intra-molecular interactions allows us to integrate this tool easily with the widely used molecular dynamics software AMBER. This makes our tool a very useful utility which can help to facilitate the study of dendrimer interaction with nucleic acids, protein and lipid bilayer for various biological applications. © 2012 Wiley Periodicals, Inc.

Dendrimer building toolkit: Model building and characterization of various dendrimer architectures

We have developed a graphical user interface based dendrimer builder toolkit (DBT) which can be used to generate the dendrimer configuration of desired generation for various dendrimer architectures. The validation of structures generated by this tool was carried out by studying the structural properties of two well known classes of dendrimers: ethylenediamine cored poly(amidoamine) (PAMAM) dendrimer, diaminobutyl cored poly(propylene imine) (PPI) dendrimer. Using full atomistic molecular dynamics (MD) simulation we have calculated the radius of gyration, shape tensor and monomer density distribution for PAMAM and PPI dendrimer at neutral and high pH. A good agreement between the available simulation and experimental (small angle X-ray and neutron scattering; SAXS, SANS) results and calculated radius of gyration was observed. With this validation we have used DBT to build another new class of nitrogen cored poly(propyl ether imine) dendrimer and study it's structural features using all atomistic MD simulation. DBT is a versatile tool and can be easily used to generate other dendrimer structures with different chemistry and topology. The use of general amber force field to describe the intra-molecular interactions allows us to integrate this tool easily with the widely used molecular dynamics software AMBER. This makes our tool a very useful utility which can help to facilitate the study of dendrimer interaction with nucleic acids, protein and lipid bilayer for various biological applications. (c) 2012 Wiley Periodicals, Inc.

Recapitulating tumour microenvironment in chitosan-gelatin three-dimensional scaffolds: an improved in vitro tumour model

Owing to the reduced co-relationship between conventional flat Petri dish culture (two-dimensional) and the tumour microenvironment, there has been a shift towards three-dimensional culture systems that show an improved analogy to the same. In this work, an extracellular matrix (ECM)-mimicking three-dimensional scaffold based on chitosan and gelatin was fabricated and explored for its potential as a tumour model for lung cancer. It was demonstrated that the chitosan-gelatin (CG) scaffolds supported the formation of tumoroids that were similar to tumours grown in vivo for factors involved in tumour-cell-ECM interaction, invasion and metastasis, and response to anti-cancer drugs. On the other hand, the two-dimensional Petri dish surfaces did not demonstrate gene-expression profiles similar to tumours grown in vivo. Further, the three-dimensional CG scaffolds supported the formation of tumoroids, using other types of cancer cells such as breast, cervix and bone, indicating a possible wider potential for in vitro tumoroid generation. Overall, the results demonstrated that CG scaffolds can be an improved in vitro tool to study cancer progression and drug screening for solid tumours.

On the dynamic contact angle in simulation of impinging droplets with sharp interface methods

Effects of dynamic contact angle models on the flow dynamics of an impinging droplet in sharp interface simulations are presented in this article. In the considered finite element scheme, the free surface is tracked using the arbitrary Lagrangian-Eulerian approach. The contact angle is incorporated into the model by replacing the curvature with the Laplace-Beltrami operator and integration by parts. Further, the Navier-slip with friction boundary condition is used to avoid stress singularities at the contact line. Our study demonstrates that the contact angle models have almost no influence on the flow dynamics of the non-wetting droplets. In computations of the wetting and partially wetting droplets, different contact angle models induce different flow dynamics, especially during recoiling. It is shown that a large value for the slip number has to be used in computations of the wetting and partially wetting droplets in order to reduce the effects of the contact angle models. Among all models, the equilibrium model is simple and easy to implement. Further, the equilibrium model also incorporates the contact angle hysteresis. Thus, the equilibrium contact angle model is preferred in sharp interface numerical schemes.

Interface capturing simulations of acoustically driven bubble dynamics in and near tissue

Tissue injury during therapeutic ultrasound or lithotripsy is thought, in cases, to be due to the action of cavitation bubbles. Assessing this and mitigating it is challenging since bubble dynamics in the complex confinement of tissues or in small blood vessels are challenging to predict. Simulations tools require specialized algorithms to simultaneously represent strong acoustic waves and shocks, topologically complex liquid‐vapor phase boundaries, and the complex viscoelastic material dynamics of tissue. We discuss advances in a simulation tool for such situations. A single‐mesh Eulerian solver is used to solve the governing equations. Special sharpening terms maintain the liquid‐vapor interface in face of the finite numerical dissipation included in the scheme to accurately capture shocks. A recent enhancement to this formulation has significantly improved this interface capturing procedure, which is demonstrated for simulation of the Rayleigh collapse of a bubble. The solver also transports elastic stresses and can thus be used to assess the effects of elastic properties on bubble dynamics. A shock‐induced bubble collapse adjacent to a model elastic tissue is used to demonstrate this and draw some conclusions regarding the injury suppressing role that tissue elasticity might play.

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.