Photochemical dynamics of laser-irradiated semiconductor nanostructures

Thesis (Ph.D.)--University of Washington, 2016-08

Properties of a composite material with mixed imperfect contact conditions

We present an analytical solution of a mixed boundary value problem for an unbounded 2D doubly periodic domain which is a model of a composite material with mixed imperfect interface conditions. We find the effective conductivity of the composite material with mixed imperfect interface conditions, and also give numerical analysis of several of their properties such as temperature and flux.

Um problema de dispersão de poluentes em rios e canais por meio do método de separação de variáveis

Neste trabalho obtém-se uma solução analítica para a equação de advecção-difusão aplicada a problemas de dispersão de poluentes em rios e canais. Para tanto, consideram-se os casos unidimensionais e bidimensionais em regime transiente com coeficientes de difusividade e velocidades constantes. A abordagem utilizada para a resolução deste problema é o método de Separação de Variáveis. Os modelos resolvidos foram simulados utilizando o MatLab. Apresentam-se os resultados das simulações numéricas em formato gráfico. Os resultados de algumas simulações numéricas existem na literatura e puderam ser comparados. O modelo proposto mostrou-se coerente em relação aos dados considerados. Para outras simulações não foram encontrados comparativos na literatura, todavia esses problemas governados por equações diferenciais parciais, mesmo lineares, não são de fácil solução analítica. Sendo que, muitas delas representam importantes problemas de matemática e física, com diversas aplicações na engenharia. Dessa forma, é de grande importância a disponibilidade de um maior número de problemas-teste para avaliação de desempenho de formulações numéricas, cada vez mais eficazes, já que soluções analíticas oferecem uma base mais segura para comparação de resultados.

Advanced Adhesion Strength Testing Methods of Thin Film Multilayers in Electronic Packaging Systems

With the continued miniaturization and increasing performance of electronic devices, new technical challenges have arisen. One such issue is delamination occurring at critical interfaces inside the device. This major reliability issue can occur during the manufacturing process or during normal use of the device. Proper evaluation of the adhesion strength of critical interfaces early in the product development cycle can help reduce reliability issues and time-to-market of the product. However, conventional adhesion strength testing is inherently limited in the face of package miniaturization, which brings about further technical challenges to quantify design integrity and reliability. Although there are many different interfaces in today's advanced electronic packages, they can be generalized into two main categories: 1) rigid to rigid connections with a thin flexible polymeric layer in between, or 2) a thin film membrane on a rigid structure. Knowing that every technique has its own advantages and disadvantages, multiple testing methods must be enhanced and developed to be able to accommodate all the interfaces encountered for emerging electronic packaging technologies. For evaluating the adhesion strength of high adhesion strength interfaces in thin multilayer structures a novel adhesion test configuration called “single cantilever adhesion test (SCAT)” is proposed and implemented for an epoxy molding compound (EMC) and photo solder resist (PSR) interface. The test method is then shown to be capable of comparing and selecting the stronger of two potential EMC/PSR material sets. Additionally, a theoretical approach for establishing the applicable testing domain for a four-point bending test method was presented. For evaluating polymeric films on rigid substrates, major testing challenges are encountered for reducing testing scatter and for factoring in the potentially degrading effect of environmental conditioning on the material properties of the film. An advanced blister test with predefined area test method was developed that considers an elasto-plastic analytical solution and implemented for a conformal coating used to prevent tin whisker growth. The advanced blister testing with predefined area test method was then extended by employing a numerical method for evaluating the adhesion strength when the polymer’s film properties are unknown.

Solução semi-analítica para modelagem das condições inicial e de contorno aplicáveis ao ensaio de difusão pura

Dissertação (mestrado)—Universidade de Brasília, Faculdade de Tecnologia, Departamento de Engenharia Civil e Ambiental, 2015.

An hr-adaptive discontinuous Galerkin method for advection-diffusion problems

We propose an adaptive mesh refinement strategy based on exploiting a combination of a pre-processing mesh re-distribution algorithm employing a harmonic mapping technique, and standard (isotropic) mesh subdivision for discontinuous Galerkin approximations of advection-diffusion problems. Numerical experiments indicate that the resulting adaptive strategy can efficiently reduce the computed discretization error by clustering the nodes in the computational mesh where the analytical solution undergoes rapid variation.

Effect of elastic modulus and poisson's ratio on guided wave dispersion using transversely isotropic material modelling

Timber poles are commonly used for telecommunication and power distribution networks, wharves or jetties, piling or as a substructure of short span bridges. Most of the available techniques currently used for non-destructive testing (NDT) of timber structures are based on one-dimensional wave theory. If it is essential to detect small sized damage, it becomes necessary to consider guided wave (GW) propagation as the behaviour of different propagating modes cannot be represented by one-dimensional approximations. However, due to the orthotropic material properties of timber, the modelling of guided waves can be complex. No analytical solution can be found for plotting dispersion curves for orthotropic thick cylindrical waveguides even though very few literatures can be found on the theory of GW for anisotropic cylindrical waveguide. In addition, purely numerical approaches are available for solving these curves. In this paper, dispersion curves for orthotropic cylinders are computed using the scaled boundary finite element method (SBFEM) and compared with an isotropic material model to indicate the importance of considering timber as an anisotropic material. Moreover, some simplification is made on orthotropic behaviour of timber to make it transversely isotropic due to the fact that, analytical approaches for transversely isotropic cylinder are widely available in the literature. Also, the applicability of considering timber as a transversely isotropic material is discussed. As an orthotropic material, most material testing results of timber found in the literature include 9 elastic constants (three elastic moduli and six Poisson's ratios), hence it is essential to select the appropriate material properties for transversely isotropic material which includes only 5 elastic constants. Therefore, comparison between orthotropic and transversely isotropic material model is also presented in this article to reveal the effect of elastic moduli and Poisson's ratios on dispersion curves. Based on this study, some suggestions are proposed on selecting the parameters from an orthotropic model to transversely isotropic condition.

Machining of coarse grained and ultra fine grained titanium

Machining of titanium is quite difficult and expensive. Heat generated in the process of cutting does not dissipate quickly, which affects tool life. In the last decade ultra fine grained (UFG) titanium has emerged as an option for substitution for more expensive titanium alloys. Extreme grain refinement can be readily performed by severe plastic deformation techniques. Grain refinement of a material achieved in this way was shown to change its mechanical and physical properties. In the present study, the microstructure evolution and the shear band formation in chips of coarse grained and UFG titanium machined to three different depths and three different feeding rates was investigated. A change in thermal characteristics of commercial purity Ti with grain refinement was studied by comparing heating/cooling measurements with an analytical solution of the heat transfer boundary problem. It was demonstrated that an improvement in the machinability can be expected for UFG titanium. © 2012 Springer Science+Business Media, LLC.

Determinations of stress wave velocity in a timber pole using wavelet transform

This paper presents an application of Wavelet Transfonn (WT) for determination of stress wave velocity for Non-destructive Testing of timber utility poles in service. For surface Non-destructive Testing (NDT), the hammer impact, which produces generally broadband frequency excitation, is used to generate stress wave. Moreover, due to practicality the impact location for field testing of a utility pole is on the side of the pole and 1.5 m above ground level. And the geometry of utility pole could not guarantee non-dispersive longitudinal wave. All of these issues have resulted in lack of accuracy and reliability of results from surface NDT in field testing. In recognition of such problem, this research explores methods to reliably calculate desired wave velocity by isolating wave mode and studying dispersive nature of utility pole. Fast Fourier Transfonn (FFT) is firstly conducted to determine the suitable frequency from a stress wave data. Then WT is applied on the wave data mentioned to perfonn time-frequency analysis. Velocity can be detennined by time history data of desired frequency from WT results which will be compared with the available analytical solution for longitudinal wave velocity. The results of the investigation showed that wavelet transfonn analysis can be a reliable signal processing tool for non-destructive testing in tenns of velocity detennination, which in tum also helps to detennine the embedded length of the timber pole.

Strain localisation patterns under equal-channel angular pressing

Instabilities of plastic flow in the form of localised shear bands were experimentally observed to result from equal-channel angular pressing (ECAP) of magnesium alloy AZ31. The appearance of shear bands and their spacing were dependent on velocity of the pressing and applied back-pressure. A generic gradient plasticity theory involving second-order strain gradient terms in a constitutive model was applied to the case of AZ31 deformed by ECAP. Linear stability analysis was applied to the set of equations describing the deformation behaviour in the process zone idealised as a planar shear zone. A full analytical solution providing a dispersion relation between the rate of growth of a perturbation and the wave number was obtained. It was shown that the pattern of incipient localised shear bands exhibits a spectrum of characteristic lengths corresponding to admissible wave numbers. The interval of the spectrum of wave numbers of viable, i.e. growing, perturbations predicted by linear stability analysis was shown to be in good agreement with the experimentally observed spectrum. The effect of back-pressure applied during ECAP was also considered. The predicted displacement of the shear band spectrum towards lower wave numbers, shown to be a result of the decreased shear strain rate in the shear zone, was consistent with the experimentally observed increase of the band spacing with increased back-pressure. A good predictive capability of the general modelling frame used in conjunction with linear stability analysis was thus demonstrated in the instance of the particular alloy system and the specific processing conditions considered.

Fast vision-based catheter 3D reconstruction

Continuum robots offer better maneuverability and inherent compliance and are well-suited for surgical applications as catheters, where gentle interaction with the environment is desired. However, sensing their shape and tip position is a challenge as traditional sensors can not be employed in the way they are in rigid robotic manipulators. In this paper, a high speed vision-based shape sensing algorithm for real-time 3D reconstruction of continuum robots based on the views of two arbitrary positioned cameras is presented. The algorithm is based on the closed-form analytical solution of the reconstruction of quadratic curves in 3D space from two arbitrary perspective projections. High-speed image processing algorithms are developed for the segmentation and feature extraction from the images. The proposed algorithms are experimentally validated for accuracy by measuring the tip position, length and bending and orientation angles for known circular and elliptical catheter shaped tubes. Sensitivity analysis is also carried out to evaluate the robustness of the algorithm. Experimental results demonstrate good accuracy (maximum errors of  ±0.6 mm and  ±0.5 deg), performance (200 Hz), and robustness (maximum absolute error of 1.74 mm, 3.64 deg for the added noises) of the proposed high speed algorithms.

A study of guided wave propagation in timber pole using spectral finite element method

Timber is one of the most widely used structural material all over the world. Round timbers can be seen as a structural component in historical buildings, jetties, short span bridges and also as piles for foundation and poles for electrical and power distribution. To evaluate the current condition of these cylindrical type timber structures, guided wave has a great potential. However, the difficulties associated with the guided wave propagation in timber materials includes orthotropic behaviour of wood, moisture contents, temperature, grain direction, etc. In addition, the effect of fully or partially filled surrounding media, such as soil, water, etc. causes attenuation on the generated stress wave. In order to investigate the effects of these parameters on guided wave propagation, extensive numerical simulation is required to conduct parametric studies. Moreover, due to the presence of multi modes in guided wave propagation, dispersion curves are of great importance. Even though conventional finite element method (FEM) can determine dispersion curves along with wave propagation in time domain, it is highly computationally expensive. Furthermore, incorporating orthotropic behaviour and surrounding media to model a thick cylindrical wave (large diameter cylindrical structures) make conventional FEM inefficient for this purpose. In contrast, spectral finite element method (SFEM) is a semi analytical method to model the guided wave propagation which does not need fine meshes compared to the other methods, such as FEM or finite difference method (FDM). Also, even distribution of mass and stiffness of structures can be obtained with very few elements using SFEM. In this paper, the suitability of SFEM is investigated to model guided wave propagation through an orthotropic cylindrical waveguide with the presence of surrounding soil. Both the frequency domain analysis (dispersion curves) and time domain reconstruction for a multi-mode generated input signal are presented under different loading location. The dispersion curves obtained from SFEM are compared against analytical solution to verify its accuracy. Lastly, different numerical issues to solve for the dispersion curves and time domain results using SFEM are also discussed.

The study of thermal stresses in a single long elastic fiber embedded in an infinite matrix

This work considered the micro-mechanical behavior of a long fiber embedded in an infinite matrix. Using the theory of elasticity, the idea of boundary layer and some simplifying assumptions, an approximate analytical solution was obtained for the normal and shear stresses along the fiber. The analytical solution to the problem was found for the case when the length of the embedded fiber is much greater than its radius, and the Young's modulus of the matrix was much less than that of the fiber. The analytical solution was then compared with a numerical solution based on Finite Element Analysis (FEA) using ANSYS. The numerical results showed the same qualitative behavior of the analytical solution, serving as a validation tool against lack of experimental results. In general this work provides a simple method to determine the thermal stresses along the fiber embedded in a matrix, which is the foundation for a better understanding of the interaction between the fiber and matrix in the case of the classical problem of thermal-stresses.

Molecular Origins of Higher Harmonics in Large-Amplitude Oscillatory Shear Flow: Shear Stress Response

Recent work has focused on deepening our understanding of the molecular origins of the higher harmonics that arise in the shear stress response of polymeric liquids in large-amplitude oscillatory shear flow. For instance, these higher harmonics have been explained by just considering the orientation distribution of rigid dumbbells suspended in a Newtonian solvent. These dumbbells, when in dilute suspension, form the simplest relevant molecular model of polymer viscoelasticity, and this model specifically neglects interactions between the polymer molecules [R.B. Bird et al., J Chem Phys, 140, 074904 (2014)]. In this paper, we explore these interactions by examining the Curtiss-Bird model, a kinetic molecular theory designed specifically to account for the restricted motions that arise when polymer chains are concentrated, thus interacting and specifically, entangled. We begin our comparison using a heretofore ignored explicit analytical solution [Fan and Bird, JNNFM, 15, 341 (1984)]. For concentrated systems, the chain motion transverse to the chain axis is more restricted than along the axis. This anisotropy is described by the link tension coefficient, ε, for which several special cases arise: ε = 0 corresponds to reptation, ε > 1/8 to rod-climbing, 1/2 ≥ ε ≥ 3/4 to reasonable predictions for shear-thinning in steady simple shear flow, and ε = 1 to the dilute solution without hydrodynamic interaction. In this paper, we examine the shapes of the shear stress versus shear rate loops for the special cases ε = (0,1/8, 3/8,1) , and we compare these with those of rigid dumbbell and reptation model predictions.