Search

Nitinol exhibits unique (thermo)mechanical properties that make it central to the design of many medical devices. However, nitinol nominally contains 50 atomic percent nickel, which if released in sufficient quantities, can lead to adverse health effects. While nickel release from nitinol devices is typically characterized using in vitro immersion tests, these evaluations require lengthy time periods. We have explored elevated temperature as a potential method to expedite this testing. Nickel release was characterized in nitinol materials with surface oxide thickness ranging from 12 to 1564 nm at four different temperatures from 310 to 360 K. We found that for three of the materials with relatively thin oxide layers, ≤ 87 nm nickel release exhibited Arrhenius behavior over the entire temperature range with activation energies of 80 to 85 kJ/mol. Conversely, the fourth ''black-oxide'' material, with a much thicker, complex oxide layer, was not well characterized by an Arrhenius relationship. Power law release profiles were observed in all four materials; however, the exponent from the thin oxide materials was approximately 1/4 compared with 3/4 for the black-oxide material. To illustrate the potential benefit of using elevated temperature to abbreviate nickel release testing, we demonstrated that a > 50 day 310 K release profile could be accurately recovered by testing for less than 1 week at 340 K. However, because the materials explored in this study were limited, additional testing and mechanistic insight are needed to establish a protective temperature scaling that can be applied to all nitinol medical device components.

Various public health concerns can arise from the unintended leaching of additives and impurities from polymeric medical devices or food packaging, which is directly related to each solute's diffusivity D. Both experimental and simulation methods can be used to quantify D, but slow diffusion at physiologic temperature in glassy polymers can render these approaches impractical. Here, we investigate a simulation approach with the potential to more rapidly calculate D. Specifically, we examine links between dynamic localization, characterized by the Debye-Waller factor, ⟨u2⟩, and D in a variety of polymer/solute systems using atomistic molecular dynamics (MD) simulations. Using short, high-temperature MD simulations to estimate D at physiologic temperature, we find that the relation ln D ∝ 1/⟨u2⟩ quantitatively predicts D for small solutes and produces an upper-bound estimate of D for larger solutes. Upper-bound estimates are useful in certain contexts, and we compare our results with another approach for determining upper bounds, the Piringer model, to show where each method may be useful. Then, we examine a modified relation where the Debye-Waller factor is rescaled by the mode coupling temperature Tc, which can produce better estimates of D if Tc is carefully chosen. Last, we compare our approach with several other models that relate temperature or localized dynamics with diffusivity. Although each of these approaches can be used to model D across wide temperature ranges using one or more adjustable parameters, none of them are truly predictive in glassy polymers. Further developments are needed to predict the optimal values of the adjustable parameters a priori.

Three-dimensional (3D) printing has enabled benchtop fabrication of customized bioengineered constructs with intricate architectures. Various approaches are being explored to enable optimum integration of such constructs into the physiological environment including addition of bioactive fillers. In this work, we incorporated a corticosteroid drug, dexamethasone (Dex), in a low modulus polyester (SC5050) and examined the effect of Dex incorporation on solvent-, initiator-, and monomer-free pneumatic extrusion-based 3D printing of the polymer. Dex-SC5050 interactions were characterized by plotting thermodynamic binary phase diagrams based on the Flory-Huggins theory. The effect of Dex composition on the 3D printability of the SC5050 polyester was examined by rheological characterization and by image analysis of each layer of the 3D printed scaffolds. The drug release and the degradation of the polymer from the 3D printed scaffolds was used to analyze the effect of Dex composition on the performance of the 3D printed scaffolds. We found that Dex was insoluble in SC5050 polyester at relevant 3D printing temperatures and the insoluble drug particles physically reinforced the polymer, increasing the viscosity and the shear modulus of the base polymer. In addition, the reinforcing effect improved the shape fidelity of the printed filaments and the overall quality of the scaffolds. The Dex particles demonstrated a two-phase release, with an initial burst release and a slower sustained release of drug under in vitro conditions. To investigate preliminary host response of the 3D printed SC5050 scaffolds for tissue engineering applications, the printed scaffolds were implanted subcutaneously in Sprague-Dawley rats for 6 weeks and examined for fibrous tissue formation, infiltration of cells, and vascularization into the pores of the scaffolds.

Ethylene oxide (EO) gas is commonly used to sterilize medical devices. Bioavailable residual EO, however, presents a significant toxicity risk to patients. Residual EO is assessed using international standards describing extraction conditions for different medical device applications. We examine a series of polymers and explore different extraction conditions to determine residual EO. Materials were sterilized with EO and exhaustively extracted in water, in one of three organic solvents, or in air using thermal desorption. The EO exhaustively extracted varies significantly and is dictated by two factors: the EO that permeates the material during sterilization; and the effectiveness of the extraction protocol in flushing residual EO from the material. Extracted EO is maximized by a close matches between Hildebrand solubility parameters δpolymer , δEO , and δsolvent . There remain complexities to resolve, however, because maximized EO uptake and detection are accompanied by great variability. These observations may inform protocols for material selection, sterilization, and EO extraction. © 2017 Wiley Periodicals, Inc. J Biomed Mater Res Part B: Appl Biomater, 106B: 2455-2463, 2018.

Liners used in orthopedic devices are often made from ultrahigh molecular weight polyethylene (UHMWPE). A general predictive capability for transport coefficients of small molecules in UHMWPE does not exist, making it difficult to assess properties associated with leaching or uptake of small molecules. To address this gap, we describe here how a form of the Vrentas-Duda free volume model can be used to predict upper-bound diffusion coefficients (D) of arbitrary molecules within UHMWPE on the basis of their size and shape. Within this framework, the free-volume microstructure of UHMWPE is defined by analysis of a curated set of model diffusants. We determined an upper limit on D for vitamin E, a common antioxidant added to UHMWPE, to be 7.1 × 10-12 cm2 s-1 . This means that a liner that contains 0.1 wt % or less Vitamin E and has

Poly (lactic-co-glycolic acid) (PLGA) copolymers have been broadly used in controlled drug release applications. Because these polymers are biodegradable, they provide an attractive option for drug delivery vehicles. There are a variety of material, processing, and physiological factors that impact the degradation rates of PLGA polymers and concurrent drug release kinetics. This work is intended to provide a comprehensive and collective review of the physicochemical and physiological factors that dictate the degradation behavior of PLGA polymers and drug release from contemporary PLGA-based drug-polymer products. In conjunction with the existing experimental results, analytical and numerical theories developed to predict drug release from PLGA-based polymers are summarized and correlated with the experimental observations. © 2016 Wiley Periodicals, Inc. J Biomed Mater Res Part B: Appl Biomater, 105B: 1692-1716, 2017.

An overview of diffuse interface models specific to drug-eluting stent coatings is presented. Microscale heterogeneities, both in the coating and use environment, dictate the performance of these coatings. Using diffuse interface methods, these heterogeneities can be explicitly incorporated into the model equations with relative ease. This enables one to predict the complex microstructures that evolve during coating fabrication and subsequent impact on drug release. Examples are provided that illustrate the wide range of phenomena that can be addressed with diffuse interface models including: crystallization, constrained phase separation, hydrolytic degradation, and heterogeneous binding. Challenges associated with the lack of material property data and numerical solution of the model equations are also highlighted. Finally, in light of these potential drawbacks, the potential to utilize diffuse interface models to help guide product and process development is discussed.

The successful function of drug-eluting devices used in the treatment of atherosclerosis relies on the concentration and retention of the drug in the vessel wall. While drug deposition necessarily depends on the underlying tissue structure, conventional models do not account for the intrinsic structural complexity of arterial tissue and its impact on deposition. By employing only average bulk material properties, the capability to predict the potential for local toxicity or therapeutic failure is limited. To address these limitations, we have developed a model that accounts explicitly for variations in the tissue structure. The approach uses a laminate approximation of the underlying microscopic structure to specify an expression for the continuous spatial dependence of the effective macroscopic material properties. Based on this continuum description, we derive an analytic expression for drug uptake into arterial tissue under typical ex vivo experimental conditions. This expression is used to extract relevant material properties for paclitaxel in bovine arteries based on available literature data. The best fit parameters are then used as the basis for numerical simulations of long-term deposition behavior from a stent with a pure paclitaxel coating. The results of these simulations are quantitatively consistent with previously reported in vivo observations. We also demonstrate that significant errors can arise in both the interpretation of experimental data and the prediction of drug deposition when structural heterogeneities are ignored. Establishing a robust deposition model can ultimately reduce empiricism in the design of drug-eluting devices, providing a facile means to guide the development and refinement of these technologies.

The drug coating process for coated drug-eluting stents (DES) has been identified as a key source of inter- and intra-batch variability in drug elution rates. Quality-by-design (QbD) principles were applied to gain an understanding of the ultrasonic spray coating process of DES. Statistically based design of experiments (DOE) were used to understand the relationship between ultrasonic atomization spray coating parameters and dependent variables such as coating mass ratio, roughness, drug solid state composite microstructure, and elution kinetics. Defect-free DES coatings composed of 70% 85:15 poly(DL-lactide-co-glycolide) and 30% everolimus were fabricated with a constant coating mass. The drug elution profile was characterized by a mathematical model describing biphasic release kinetics. Model coefficients were analyzed as a DOE response. Changes in ultrasonic coating processing conditions resulted in substantial changes in roughness and elution kinetics. Based on the outcome from the DOE study, a design space was defined in terms of the critical coating process parameters resulting in optimum coating roughness and drug elution. This QbD methodology can be useful to enhance the quality of coated DES.

A novel approach for rapid risk assessment of targeted leachables in medical device polymers is proposed and validated. Risk evaluation involves understanding the potential of these additives to migrate out of the polymer, and comparing their exposure to a toxicological threshold value. In this study, we propose that a simple diffusive transport model can be used to provide conservative exposure estimates for phase separated color additives in device polymers. This model has been illustrated using a representative phthalocyanine color additive (manganese phthalocyanine, MnPC) and polymer (PEBAX 2533) system. Sorption experiments of MnPC into PEBAX were conducted in order to experimentally determine the diffusion coefficient, D = (1.6 ± 0.5) × 10−11 cm2/s, and matrix solubility limit, C s = 0.089 wt.%, and model predicted exposure values were validated by extraction experiments. Exposure values for the color additive were compared to a toxicological threshold for a sample risk assessment. Results from this study indicate that a diffusion model-based approach to predict exposure has considerable potential for use as a rapid, screening-level tool to assess the risk of color additives and other small molecule additives in medical device polymers.

An aqueous hyaluronic acid (HA(aq)) pericellular coat, when mediating the tactile aspect of cellular contact inhibition, has three tasks: interface formation, mechanical signal transmission and interface separation. To quantify the interfacial adhesive behavior of HA(aq), we induce simultaneous interface formation and separation between HA(aq) and a model hydrophobic, hysteretic Si-SAM surface. While surface tension γ remains essentially constant, interface formation and separation depend greatly on concentration (5 ≤ C ≤ 30 mg mL(-1)), molecular weight (6 ≤ MW ≤ 2000 kDa) and interfacial velocity (0 ≤ V ≤ 3 mm s(-1)), each of which affect shear elastic and loss moduli G′ and G′′, respectively. Viscoelasticity dictates the mode of interfacial motion: wetting-dewetting, capillary necking, or rolling. Wetting-dewetting is quantified using advancing and receding contact angles θ(A) and θ(R), and the hysteresis between them, yielding data landscapes for each C above the [MW, V] plane. The landscape sizes, shapes, and curvatures disclose the interplay, between surface tension and viscoelasticity, which governs interfacial dynamics. Gel point coordinates modulus G and angular frequency ω appear to predict wetting-dewetting (G75 ω0.2), capillary necking (75 ω0.2G200 ω0.075) or rolling (G200ω0.075). Dominantly dissipative HA(aq) sticks to itself and distorts irreversibly before separating, while dominantly elastic HA(aq) makes contact and separates with only minor, reversible distortion. We propose the dimensionless number (G′V)/(ω(r)γ), varying from 10(-5) to 10(3) in this work, as a tool to predict the mode of interface formation-separation by relating interfacial kinetics with bulk viscoelasticity. Cellular contact inhibition may be thus aided or compromised by physiological or interventional shifts in [C, MW, V], and thus in (G′V)/(ω(r)γ), which affect both mechanotransduction and interfacial dynamics. These observations, understood in terms of physical properties, may be broadened to probe interfacial dynamics of other viscoelastic aqueous biopolymers.

We have devised a novel diffuse interface formulation to model the development of chemical and physical inhomogeneities, i.e. microstructure, during the process of casting drug-eluting coatings. These inhomogeneities, which depend on the coating constituents and manufacturing conditions, can have a profound affect on the rate and extent of drug release, and therefore the ability of coated medical devices to function successfully. By deriving the model equations in a time-dependent reference frame, we find that it is computationally viable to probe a wide, physically relevant range of material and process quantities. To illustrate the application of the model, we have evaluated the impact of manufacturing solvent, coating thickness and evaporation rate on microstructure development. Our results suggest that modifying these process conditions can have a strong and nearly discontinuous effect on coating microstructure, and therefore on drug release. Further, we demonstrate that the model can be applied to processes that involve the incremental application of the coating in layers or passes. This new model formulation, which can also be used to predict the kinetics of drug release, provides a tool to elucidate and quantify the relationships between process variables, microstructure and performance. Establishing these relationships can reduce empiricism in materials selection and process design, providing a facile and efficient means to tailor the underlying microstructure and achieve a desired drug-release behavior.

With the advent of nanotechnology, silver nanoparticles increasingly are being used in coatings, especially in medical device applications, to capitalize on their antimicrobial properties. The attractiveness of nanoparticulate silver systems is the expected increased antimicrobial efficacy relative to their bulk counterparts, which may be attributed to an increased silver ion (Ag+) solubility, and hence availability, that arises from capillarity effects in small, nanometer-sized particles. However, a change of the material upon which the antimicrobial nanoparticulate silver is deposited (herein called "substrate") may affect the availability of Ag+ ions and the intended efficacy of the device. We utilize both theory and experiment to determine the effect of substrate on ion release from silver particles in electrochemical environments and find that substrate surface charge, chemical reactivity or affinity of the surface for Ag+ ions, and wettability of the surface all affect availability of Ag+ ions, and hence antimicrobial efficacy. It is also observed that with time of exposure to deionized water, Ag+ ion release increases to a maximum value at 5 min before decreasing to undetectable levels, which is attributed to coarsening of the nanoparticles, and which subsequently reduces the solubility and availability of Ag+ ions. This coarsening phenomenon is also predicted by the theoretical considerations and has been confirmed experimentally by transmission electron microscopy.

In recent years, traditional rigid materials have been replaced with softer elastomers in intraocular lenses to minimize the size of the required surgical incision, thereby reducing patient recuperation time. However, water-filled cavities that may impact visual acuity are found in many of these new implants. We demonstrate that the cavitation observed in vivo can occur due to an osmotic pressure difference between the aqueous solution within the cavity and the external media in which the lens is immersed. By reducing the osmolarity of the external solution from 300 to 0mM, we observe an increase in cavity volume of almost a factor of 30. Further, we have developed a model for cavity growth assuming the controlling factor is diffusion of hydrophilic molecules from the polymer matrix into the cavity. We find that the experimental observations are consistent with the model and suggest that oligomeric species generated during polymerization are responsible for the observed cavitation.

Quantitative homology metrics are proposed for characterizing the thermal–elastic response of polycrystalline materials. Simulations for a calcite-based polycrystal, marble, are used as an illustrative example. The homology metrics are based on topological measurements, such as the number of components and the number of handles of the thermal–elastic response fields for a complex microstructure. These homology metrics are applied to characterize not only the elastic energy density and maximum principal stress response fields in a polycrystal but also the correlated grain-boundary misorientation distributions that influenced the formation of these response fields. It is demonstrated that the topological analysis can quantitatively distinguish between different types of grain-boundary misorientations, as well as between differences in the resulting response fields.

In recent years, controlled release coatings, comprised of drug-polymer composites, have been integrated with medical devices, improving device functionality and performance. However, relationships between material properties, manufacturing environment, composite (micro)structure, and subsequent release kinetics are not well established. We apply a thermodynamically consistent model to probe the influence of drug-polymer chemistry (phobicity), drug loading, and evaporation rate on microstructure development during fabrication. For these structures, we compute release profiles for exposure to polymer-insoluble media and media in which the polymer readily dissolves. We find that with increasing drug-polymer phobicity, structural heterogeneities form at lower loadings and more rapid rates. The heterogeneities remain isolated and compact at low loadings and become interconnected as the drug to polymer ratio approaches 1.0. Release into polymer-insoluble media was dramatically enhanced by heterogeneities, resulting in up to a fourfold increase in drug release. In polymer-soluble media, however, heterogeneities diminished release. Although reductions of only 30% were typically observed, the absolute changes were much larger than observed in polymer-insoluble media. Our results suggest that improved comprehension and quantification of the physico-chemical properties in controlled release systems will enable the microstructure to be tailored to achieve desired responses that are insensitive to manufacturing variations.

Two-dimensional, microstructure-based finite element simulations were used to elucidate the influence of grain orientation configuration on the thermal-elastic response of polycrystalline ceramic materials. The two main constituent minerals of marbles, calcite and dolomite, were considered. The crystallographic configuration of the grains is described in terms of the distribution of grain orientations and grain-boundary misorientations. To probe the influence of grain orientation configuration, we first generated the geometry of a hypothetical microstructure. Next, crystallographic orientations were assigned to each grain in the microstructure such that the grain orientations and grain-boundary misorientations matched predefined distributions. By varying the predefined distributions, we generated 45 unique microstructures covering a wide range of crystallographic configurations. After assigning thermal-elastic properties to each structure, corresponding to either calcite or dolomite, finite-element simulations were performed. The simulations demonstrated that both the orientation and misorientation distributions have a substantial impact on the thermal-elastic response of microstructures with the same material properties: varying the crystallographic configuration results in variations of stored elastic strain energy density from 10 kJ m−3 to 39 kJ m−3 for a 100 °C temperature change. Further, the ratio of the bulk coefficient of thermal expansion in the two principal directions varies from − 0.17 to 1.46. Surprisingly, the grain-boundary misorientations alone had a substantial impact on the thermal-elastic response of the system. By simply rearranging a fixed set of crystallographic orientations to obtain different nearest-neighbor misorientation configurations, we observed variations in strain energy density from 10 kJ m−3 to 39 kJ m−3 and variations in thermal expansion ratio from 0.06 to 0.90. The results suggest that to predict accurately the thermal-elastic response of polycrystalline ceramics, it is critical to consider both the distribution of grain-boundary misorientations as well as the distribution of grain orientations.

This paper describes some aspects of reconstruction of microstructures in three dimensions. A distinction is drawn between tomographic approaches that seek to characterize specific volumes of material, either with or without diffraction, and statistical approaches that focus on particular aspects of microstructure. A specific example of the application of the statistical approach is given for an aerospace aluminum alloy in which the distributions of coarse constituent particles are modeled. Such distributions are useful for modeling fatigue crack initiation and propagation.

We have used electron backscattered diffraction patterns to determine the misorientation of 201 adjacent pairs of grains in a magnesia polycrystal. The width and depth of the thermal grooves formed by these same grain boundaries were also measured by atomic force microscopy (AFM). By simulating the errors associated with the AFM observations and comparing our observations to existing data for magnesia and alumina, we show that, under appropriate experimental conditions, surface dihedral angles, relative grain-boundary energies, and surface diffusivities determined from AFM measurements are consistent with data acquired by more laborious techniques. Correlation of the grain-boundary misorientation and thermal groove geometry leads to the observation that grain boundaries with small misorientations, regardless of the rotation axis, have shallow thermal grooves and relatively low grain-boundary energies. Furthermore, numerous boundaries with relatively large misorientations but shallow thermal grooves correspond to special boundaries near coincident-site-lattice (CSL) misorientations. Finally, the data set indicates that factors other than the boundary misorientation, such as anisotropy of the surface energy and the grain-boundary tangent plane, play a role in determining the groove geometry.

A stereological method is described for estimating the distribution of grain-boundary types in polycrystalline materials on the basis of observations from a single planar section. The grain-boundary distribution is expressed in terms of five macroscopically observable parameters that include: three parameters that describe the lattice misorientation across the boundary and two parameters that describe the orientation of the grain-boundary plane normal. The grain-boundary distribution is derived from measurements of grain orientations and the orientations of the lines formed where grain boundaries intersect the plane of observation. Tests of the method on simulated observations illustrate that the distribution of boundaries in a material with cubic symmetry can be reliably determined with about 10° of resolution from the analysis of 5×104 or more line segments. Furthermore, grain-boundary distributions directly observed from serial sections of a SrTiO3 polycrystal are compared to those resulting from the stereological analysis of a single plane. The comparison shows that the stereological method provides a reasonable estimate of the measured distribution. The differences between the directly observed grain-boundary distribution and that derived from the stereological analysis are consistent with the results from the simulation.

Techniques are described that have been used to create a statistically representative three-dimensional model microstructure for input into computer simulations using the geometric and crystallographic observations from two orthogonal sections through an aluminum polycrystal. Orientation maps collected on the observation planes are used to characterize the sizes, shapes, and orientations of grains. Using a voxel-based tessellation technique, a microstructure is generated with grains whose size and shape are constructed to conform to those measured experimentally. Orientations are then overlaid on the grain structure such that distribution of grain orientations and the nearest-neighbor relationships, specified by the distribution of relative misorientations across grain boundaries, match the experimentally measured distributions. The techniques are applicable to polycrystalline materials with sufficiently compact grain shapes and can also be used to controllably generate a wide variety of hypothetical microstructures for initial states in computer simulations.

Measurements of the grain boundary population as a function of misorientation and boundary plane orientation show that the distribution is inversely correlated to the sum of the energies of the surfaces comprising each boundary. The observed correlation suggests that the difference between the energy of a high-angle grain boundary and the two component surfaces is relatively constant as a function of misorientation. Two exceptions to this correlation were identified: low-misorientation-angle boundaries and the coherent twin boundary, where the (111) planes in the adjoining crystals are parallel to each other, but rotated by 60° around the [111] axis. In these cases, the high degree of coincidence across this interface probably lowers the boundary energy with respect to that of the component surfaces. For all other boundaries, the anisotropy of the population is accurately predicted by the surface energy anisotropy, and in general, boundaries display a preference for {100} orientations, the planes of minimum surface energy.

Recent advances both in experimental instrumentation and computing power have made it possible to interrogate the distribution of internal interfaces in polycrystals and the three dimensional structure of the grain boundary network with an unprecedented level of detail. The purpose of this paper is to review techniques that can be used to study the mesoscopic crystallographic structure of grain boundary networks and to summarize current findings. Recent studies have shown that grain surfaces within dense polycrystals favor the same low energy planes that are found on equilibrium crystal shapes and growth forms of crystals in contact with another phase. In the materials for which comprehensive data exists, the distribution of grain boundaries is inversely correlated to the sum of the energies of the surfaces of the grains on either side of the boundary.

Using measurements of the geometric and crystallographic characteristics of approximately 104 triple junctions in a MgO polycrystal, the relative grain boundary energy has been determined as a function of five macroscopic parameters. The relative energy of a particular grain boundary is inversely correlated with its frequency of occurrence. At all misorientations, grain boundaries with {1 0 0} interface planes have relatively low energies. For low misorientation angle grain boundaries, the results are consistent with the predictions of dislocation models. At high misorientation angles, the population and energy are correlated to the sum of the energies of the free surfaces that comprise the boundary.

A semi-automated method has been used to measure all five macroscopically observable parameters of 4.1×10 6 boundary plane segments making up 5.4 mm 2 of boundary area in a hot-pressed magnesia polycrystal. The observations allow a complete description of the distribution of crystal orientations, grain boundary misorientations, and the crystallographic orientations of grain boundary planes. Among the low misorientation angle grain boundaries, there is a preference for tilt boundaries, especially those with boundary plane normals in the direction. At all fixed misorientations, there is a preference for boundaries with a boundary plane normal in the direction. These boundaries are generally asymmetric and occur at least twice as frequently as the average boundary for each fixed misorientation.

We have developed a technique that allows all five macroscopically observable grain boundary degrees of freedom to be characterized for a statistically significant number of interfaces. Using this technique, we have characterized 5 x 10 μm of grain boundary plane area in a 0.15 mm volume of MgO. The observations demonstrate that there is significant texture in the distribution of boundary planes. Here we compare the observed distribution of grain boundary planes at low misorientation angles (~ 5 °) to the calculated geometrically necessary dislocation content of the same interfaces. Based on the inverse correlation between these two quantities, we conclude that relatively low energy configurations are adopted with the highest frequency.

Microstructure-level residual stresses arise in polycrystalline ceramics during processing as a result of thermal expansion anisotropy and crystallographic disorientation across the grain boundaries. Depending upon the grain size, the magnitude of these stresses can be sufficiently high to cause spontaneous microcracking during the processing of these materials. They are also likely to affect where cracks initiate and propagate under macroscopic loading. The magnitudes of residual stresses in untextured and textured alumina samples were predicted using object oriented finite (OOF) element analysis and experimentally determined grain orientations. The crystallographic orientations were obtained by electron-backscattered diffraction (EBSD). The residual stresses were lower and the stress distributions were narrower in the textured samples compared to those in the untextured samples. Crack initiation and propagation were also simulated using the Griffith fracture criterion. The grain boundary to surface energy ratios required for computations were estimated using AFM groove measurements.

By measuring the geometry and crystallography of the three interfaces that meet at grain boundary thermal grooves, it is possible to determine the anisotropy of the surface free energy. Previously, the surface energy of MgO at 1400°C in air was approximated by a truncated double Fourier series with coefficients that were determined by fitting the observations to Herring's condition for local equilibrium at a triple junction. The purpose of this paper is to describe an alternative analysis of the same data set that is not limited by an assumed functional form of the surface energy. In this case, the space of surface characters is discretized and each orientation is associated with a capillarity vector (according to the Cahn–Hoffmann definition). The set of capillarity vectors that most closely satisfies the condition for local equilibrium at each triple junction is then determined by an iterative method. The relative surface free energies derived from this analysis are more anisotropic than those derived from the series fit and more consistent with the observed faceting of MgO in air at 1400°C. The relative surface energies of the low index planes are γ110/γ100 = 1.07±0.04 and γ111/γ100 = 1.17±0.04.

Geometric and crystallographic data obtained from a well annealed magnesia polycrystal have been used to specify the five macroscopic degrees of freedom for 4665 grain boundaries. The results indicate, that for this sample, the five parameter grain boundary character space is fully occupied. A finite series of symmetrized spherical harmonics has been used to approximate the misorientation dependence of the relative grain boundary energy. Best fit coefficients for this series were determined by assuming that the interfacial tensions at each triple junction are balanced. The grain boundary energy function shows Read-Shockley behavior at small misorientations and a broad minimum near the Σ3 misorientation. Furthermore, misorientations about the ‹100› axis create boundaries with relative energies that are less than those created by misorientations about the ‹110› or ‹111› axes.

Measurement of the geometry of triple junctions between grain boundaries in polycrystalline materials generates large sets of dihedral angles from which maps of the grain boundary energy may be extracted. A preliminary analysis has been performed for a sample of magnesia based on a three-parameter description of grain boundaries. An extended form of orientation imaging microscopy (OIM) was used to measure both triple junction geometry via image analysis in the SEM and local grain orientation via electron back scatter diffraction. Serial sectioning with registry of both in-plane images and successive sections characterizes triple junction tangents from which true dihedral angles are calculated. We apply Herring's relation at each triple junction, based on the assumption of local equilibrium at the junction. By limiting grain boundary character to a (three parameter) specification of misorientation for the preliminary analysis, we can neglect the torque terms and apply the sine law to the three boundaries. This provides two independent relations per triple junction between grain boundary energies and dihedral angles. Discretizing the misorientation and employing multiscale statistical analysis on large data sets allows (relative) grain boundary energy as a function of boundary character to be extracted from triple junction geometry. A similar analysis of thermal grooves allows the anisotropy of the surface energy to be measured in MgO.

The purpose of the present work is to model corrosion of implant materials (such as nitinol or stainless steel) in simulated environments both under static and dynamic conditions. The present model seeks to imitate the potentiodynamic testing protocol defined as in ASTM F2129, which is used to evaluate the corrosion susceptibility of small implant devices. The model is currently limited to one dimension and hence the effect of physical features, such as defects and crevices, on corrosion is not the focus of this study.

The multiplicity of distinct grain boundary configurations in polycrystals has made it difficult to determine the relative frequency with which each configuration is adopted. As a result, the physiochemical properties of each boundary and the influence of the distribution of boundaries on macroscopic materials properties are not well understood. Using a semiautomated system, we have measured all five macroscopically observable degrees of freedom of 4.1 × 106 boundary plane segments making up 5.2 × 106μm2 of grain boundary interface area in a magnesia polycrystal. Our observations demonstrate that not all grain boundary configurations occur with the same frequency and that the relative free energies of the different interfacial configurations influence the population distribution. Furthermore, the results indicate that relative grain boundary energies can be estimated based on the free surface energies.

We show here that the boundaries of individual grains in dense polycrystals prefer certain crystallographic habit planes, almost as if they were independent of the neighboring crystals. In MgO, SrTiO3, MgAl2O4, TiO2, and aluminum, the specific habit planes within the polycrystal correspond to the same planes that dominate the external growth forms and equilibrium shapes of isolated crystals of the same phase. The observations decrease the apparent complexity of interfacial networks and suggest that the mechanisms of solid-state grain growth may be analogous to conventional crystal growth. The results also indicate that a model for grain-boundary energy and structure based on grain surface relationships is more appropriate than the widely accepted models based on lattice orientation relationships.

A critical metrology issue for pharmaceutical industries is the application of analytical techniques for the characterization of drug delivery systems to address interrelationships between processing, structure, and drug release. In this study, cast coatings were formed from solutions of poly(styrene- b -isobutylene- b -styrene) (SIBS) and tetracycline in tetrahydrofuran (THF). These coatings were characterized by several imaging modalities, including time- of-flight secondary ion mass spectrometry (TOF-SIMS) for chemical imaging and analysis, atomic force microscopy (AFM) for determination of surface structure and morphology, and laser scanning confocal microscopy (LSCM), which was used to characterize the three-dimensional structure beneath the surface. The results showed phase separation between the drug and copolymer regions. The size of the tetracycline phase in the polymer matrix ranged from hundreds of nanometers to tens of microns, depending on coating composition. The mass of drug released was not found to be proportional to drug loading, because the size and spatial distribution of the drug phase varied with drug loading and solvent evaporation rate, which in turn affected the amount of drug released.

Drug delivery from biodegradable polymer coatings is becoming an important area of research for applications including medical devices. In this work, we probe the impact of polymer chemistry and solvent evaporation rate on drug morphology and the subsequent elution from biodegradable polymer systems. Two different formulations of poly(lactic-glycolic-acid) (PLGA) polymers are used in combination with tetracycline (TC) to fabricate coatings using two different solvent evaporation rates. We have conducted elution studies to quantify the release behavior of tetracycline (TC) at the first 2-3 days of dissolution. Significantly we have correlated the drug release kinetics with the drug microstructure of TC/PLGA coating using atomic force microscopy (AFM) and laser scanning confocal microscopy (LSCM). These results suggest that polymer chemistry affects the rate of water absorption and drug dissolution which in turn alters the rate of TC release during the early stages of drug elution.

To improve functionality and performance, controlled drug-release coatings comprised of drug and polymer are integrated with traditional medical devices, e.g., drug eluting stents. Depending on manufacturing conditions, these coatings can exhibit complex microstructures. Previously, a thermodynamically consistent model was developed for microstructure evolution in these systems to establish relationships between process variables, microstructure, and the subsequent release kinetics. Calculations based on the model were, in general, consistent with experimental findings. However, because of assumptions regarding the evaporation of solvent during fabrication, the model was unable to capture variations through the coating thickness that are observed experimentally. Here, a straightforward method is introduced to incorporate solvent evaporation explicitly into the model. Calculations are used to probe the impact of solvent evaporation rate and drug loading on the microstructure that forms during manufacturing and subsequent drug release kinetics. The predicted structures and release kinetics are found to be consistent with experimental observations. Further, the calculations demonstrate that solvent evaporation rate can be as critical to device performance as the amount of drug within the coating. For example, changes of a factor of five in the amount of drug released were observed by modifying the rate of solvent evaporation during manufacturing.

A phase eld approach is developed to model wetting and heterogeneous crystal nucleation of an undercooled pure liquid in contact with a sharp wall. We discuss various choices for the boundary condition at the wall and determine the properties of critical nuclei, including their free energy of formation and the contact angle as a function of undercooling. We nd for particular choices of boundary conditions, we may realize either an analog of the classical spherical cap model or decidedly non-classical behavior, where the contact angle decreases from its value taken at the melting point towards complete wetting at a critical undercooling., 4 figures

Hematopoietic stem cell (HSC) therapy can significantly lower instances of infection in chemotherapy patients by accelerating the recovery of white blood cells in the body. However, therapy requires that HSCs be stored at cryogenic temperatures to retain the cells’ ability to proliferate. Currently, cells are stored in polymeric blood bags that are subject to fracture at the extremely low storage temperatures, which leads to cell contamination, thereby reducing their effectiveness. Therefore, we have developed an analytical model to predict the accumulation of stresses that ultimately lead to crack initiation and bag fracture during cryogenic storage. Our model gives explicit relationships between stress state in the container and thermoelastic properties of the container material, container geometry, and environmental factors that include temperature of the system and pressure induced by excess gas evolving from the stored medium. Predictions based on the model are consistent with experimental observations of bag failures that occurred during cryogenic storage applications. Finally, the model can provide guidance in material selection and bag design to fabricate bags that will be less susceptible to fracture.

The grain boundary plane distributions in MgO, SrTiO3, MgAl2O4, and Al are compared at lattice misorientations with a coincident site density of greater than or equal to 1/9. In most situations, the most frequently adopted grain boundary orientation is a habit plane of low index and low surface energy that depends on the particular material. Cases where the most common boundary orientation is a plane of high planar coincident site density instead of a characteristic habit plane are rare. In fact, in most cases, the distributions of grain boundary planes at misorientations with high lattice coincidence are not substantially different from the distributions at other, more general misorientations. The results indicate that a model for grain boundary energy and structure based on grain surface relationships is more appropriate than the widely accepted models based on lattice orientation relationships.

Measurement of the geometry of triple junctions between grain boundaries in polycrystalline materials is used to generate large sets of dihedral angles from which maps of the grain boundary energy are extracted. A preliminary analysis has been performed for samples of magnesia and aluminum based on a three-parameter description of grain boundaries. An extended form of orientation imaging microscopy (OIM) was used to measure both triple junction geometry via image analysis in the SEM and local grain orientation via electron back scatter diffraction. Serial sectioning with registry of both in-plane images and successive sections characterizes triple junction tangents from which true dihedral angles are calculated. If there is local equilibrium at each triple junction, we may apply Herring's relation. By limiting grain boundary character to a (three parameter) specification of misorientation for the preliminary analysis, we can neglect the torque terms and apply the sine law to the three boundaries. This provides two independent relations per triple junction between grain boundary energies and dihedral angles. By discretizing the misorientation and employing multiscale statistical analysis on large data sets, (relative) grain boundary energy as a function of boundary character can be extracted from triple junction geometry. The results are discussed with respect to current understanding of grain boundary structure based on their crystallography. The results suggest that a three parameter characterization of grain boundaries (lattice disorientation) is not an adequate description of boundary character. A full analysis including torque terms and a five parameter boundary description is under development.