Analysis and optimization of stress wave propagation in two-dimensional granular crystals with defects

Granular crystals are compact periodic assemblies of elastic particles in Hertzian contact whose dynamic response can be tuned from strongly nonlinear to linear by the addition of a static precompression force. This unique feature allows for a wide range of studies that include the investigation of new fundamental nonlinear phenomena in discrete systems such as solitary waves, shock waves, discrete breathers and other defect modes. In the absence of precompression, a particularly interesting property of these systems is their ability to support the formation and propagation of spatially localized soliton-like waves with highly tunable properties. The wealth of parameters one can modify (particle size, geometry and material properties, periodicity of the crystal, presence of a static force, type of excitation, etc.) makes them ideal candidates for the design of new materials for practical applications. This thesis describes several ways to optimally control and tailor the propagation of stress waves in granular crystals through the use of heterogeneities (interstitial defect particles and material heterogeneities) in otherwise perfectly ordered systems. We focus on uncompressed two-dimensional granular crystals with interstitial spherical intruders and composite hexagonal packings and study their dynamic response using a combination of experimental, numerical and analytical techniques. We first investigate the interaction of defect particles with a solitary wave and utilize this fundamental knowledge in the optimal design of novel composite wave guides, shock or vibration absorbers obtained using gradient-based optimization methods.

Stress-induced changes in the Schizosaccharomyces pombe proteome using two-dimensional difference gel electrophoresis, mass spectrometry and a novel integrated robotics platform

Robotic and manual methods have been used to obtain identification of significantly changing proteins regulated when Schizosaccharomyces pombe is exposed to oxidative stress. Differently treated S. pombe cells were lysed, labelled with CyDye and analysed by two-dimensional difference gel electrophoresis. Gel images analysed off-line, using the DeCyder image analysis software [GE Healthcare, Amersham, UK] allowed selection of significantly regulated proteins. Proteins displaying differential expression were excised robotically for manual digestion and identified by matrix-assisted laser desorption/ionisation - mass spectrometry (MALDI-MS). Additionally the same set of proteins displaying differential expression were automatically cut and digested using a prototype robotic platform. Automated MALDI-MS, peak label assignment and database searching were utilised to identify as many proteins as possible. The results achieved by the robotic system were compared to manual methods. The identification of all significantly altered proteins provides an annotated peroxide stress-related proteome that can be used as a base resource against which other stress-induced proteomic changes can be compared.

Stress-induced changes in the Schizosaccharomyces pombe proteome using two-dimensional difference gel electrophoresis, mass spectrometry and a novel integrated robotics platform

Robotic and manual methods have been used to obtain identification of significantly changing proteins regulated when Schizosaccharomyces pombe is exposed to oxidative stress. Differently treated S. pombe cells were lysed, labelled with CyDye (TM) and analysed by two-dimensional difference gel. electrophoresis. Gel images analysed off-line, using the DeCyder (TM) image analysis software [GE Healthcare, Amersham, UK] allowed selection of significantly regulated proteins. Proteins displaying differential expression were excised robotically for manual digestion and identified by matrix-assisted laser desorption/ionisation - mass spectrometry (MALDI-MS). Additionally the same set of proteins displaying differential expression were automatically cut and digested using a prototype robotic platform. Automated MALDI-MS, peak label assignment and database searching were utilised to identify as many proteins as possible. The results achieved by the robotic system were compared to manual methods. The identification of all significantly altered proteins provides an annotated peroxide stress-related proteome that can be used as a base resource against which other stress-induced proteomic changes can be compared.

Comparison of Single-Standing or Connected Implants on Stress Distribution in Bone of Mandibular Overdentures: A Two-Dimensional Finite Element Analysis

This study aimed to compare the influence of single-standing or connected implants on stress distribution in bone of mandibular overdentures by means of two-dimensional finite element analysis. Two finite element models were designed using software (ANSYS) for 2 situations: bar-clip (BC) group-model of an edentulous mandible supporting an overdenture over 2 connected implants with BC system, and o'ring (OR) group-model of an edentulous mandible supporting an overdenture over 2 single-standing implants with OR abutments. Axial loads (100 N) were applied on either central (L1) or lateral (L2) regions of the models. Stress distribution was concentrated mostly in the cortical bone surrounding the implants. When comparing the groups, BC (L1, 52.0 MPa and L2, 74.2 MPa) showed lower first principal stress values on supporting tissue than OR (L1, 78.4 MPa and L2, 76.7 MPa). Connected implants with BC attachment were more favorable on stress distribution over peri-implant-supporting tissue for both loading conditions.

Two-dimensional graphene heterojunctions: The tunable mechanical properties

We report the mechanical properties of different two-dimensional carbon heterojunctions (HJs) made from graphene and various stable graphene allotropes, including α-, β-, γ- and 6612-graphyne (GY), and graphdiyne (GDY). It is found that all HJs exhibit a brittle behaviour except the one with α-GY, which however shows a hardening process due to the formation of triple carbon rings. Such hardening process has greatly deferred the failure of the structure. The yielding of the HJs is usually initiated at the interface between graphene and graphene allotropes, and monoatomic carbon rings are normally formed after yielding. By varying the locations of graphene (either in the middle or at the two ends of the HJs), similar mechanical properties have been obtained, suggesting insignificant impacts from location of graphene allotropes. Whereas, changing the types and percentages of the graphene allotropes, the HJs exhibit vastly different mechanical properties. In general, with the increasing graphene percentage, the yield strain decreases and the effective Young’s modulus increases. Meanwhile, the yield stress appears irrelevant with the graphene percentage. This study provides a fundamental understanding of the tensile properties of the heterojunctions that are crucial for the design and engineering of their mechanical properties, in order to facilitate their emerging future applications in nanoscale devices, such as flexible/stretchable electronics.

Unsteady Incompressible Two-Dimensional and Axisymmetric Turbulent Boundary Layer Flows

The unsteady turbulent incompressible boundary-layer flow over two-dimensional and axisymmetric bodies with pressure gradient has been studied. An eddy-viscosity model has been used to model the Reynolds shear stress. The unsteadiness is due to variations in the free stream velocity with time. The nonlinear partial differential equation with three independent variables governing the flow has been solved using Keller's Box method. The results indicate that the free stram velocity distribution exerts strong influence on the boundary-layer characteristics. The point of zero skin friction is found to move upstream as time increases.

A two-dimensional model for crack closure effect in plates under bending

A two-dimensional model is proposed for taking into account the establishment of contact on the compression side of crack faces in plates under bending. An approximate but simple method is developed for evaluating reduction of stress intensity factor due to such ‘crack closure’. Analysis is first carried out permitting interference of the crack faces. Contact forces are then introduced on the crack faces and their magnitudes determined from the consideration that the interference is just eliminated. The method is based partly on finite element analysis and partly on a continuum analysis using Irwin's solution for point loads on the crack line.

On the criteria for reverse transition in a two-dimensional boundary layer flow

Experiments on reverse transition were conducted in two-dimensional accelerated incompressible turbulent boundary layers. Mean velocity profiles, longitudinal velocity fluctuations $\tilde{u}^{\prime}(=(\overline{u^{\prime 2}})^{\frac{1}{2}})$ and the wall-shearing stress (TW) were measured. The mean velocity profiles show that the wall region adjusts itself to laminar conditions earlier than the outer region. During the reverse transition process, increases in the shape parameter (H) are accompanied by a decrease in the skin friction coefficient (Cf). Profiles of turbulent intensity (u’2) exhibit near similarity in the turbulence decay region. The breakdown of the law of the wall is characterized by the parameter \[ \Delta_p (=\nu[dP/dx]/\rho U^{*3}) = - 0.02, \] where U* is the friction velocity. Downstream of this region the decay of $\tilde{u}^{\prime}$ fluctuations occurred when the momentum thickness Reynolds number (R) decreased roughly below 400.

Suspension mechanics: I. Inertial and non-Newtonian migration of neutrally buoyant rigid spheres in two-dimensional unidirectional flows; II. The effect of viscoelasticity on the creeping motion of a train of neutrally buoyant Newtonian drops through a circular tube

The lateral migration of neutrally buoyant rigid spheres in two-dimensional unidirectional flows was studied theoretically. The cases of both inertia-induced migration in a Newtonian fluid and normal stress-induced migration in a second-order fluid were considered. Analytical results for the lateral velocities were obtained, and the equilibrium positions and trajectories of the spheres compared favorably with the experimental data available in the literature. The effective viscosity was obtained for a dilute suspension of spheres which were simultaneously undergoing inertia-induced migration and translational Brownian motion in a plane Poiseuille flow. The migration of spheres suspended in a second-order fluid inside a screw extruder was also considered.

The creeping motion of neutrally buoyant concentrically located Newtonian drops through a circular tube was studied experimentally for drops which have an undeformed radius comparable to that of the tube. Both a Newtonian and a viscoelastic suspending fluid were used in order to determine the influence of viscoelasticity. The extra pressure drop due to the presence of the suspended drops, the shape and velocity of the drops, and the streamlines of the flow were obtained for various viscosity ratios, total flow rates, and drop sizes. The results were compared with existing theoretical and experimental data.

Simulations of Two-Dimensional Foam under Couette Shear

null Sponsorship: Financial support is acknowledged from the University of Wales Aberystwyth Senate Fund, the Ulysses France-Ireland Exchange Scheme and EPSRC (EP/D014956/1, EP/D048397/1, EP/D071127/1).

Orientation, Flow, and Clogging in a Two-Dimensional Hopper: Ellipses vs. Disks

Two-dimensional (2D) hopper flow of disks has been extensively studied. Here, we investigate hopper flow of ellipses with aspect ratio $\alpha = 2$, and we contrast that behavior to the flow of disks. We use a quasi-2D hopper containing photoelastic particles to obtain stress/force information. We simultaneously measure the particle motion and stress. We determine several properties, including discharge rates, jamming probabilities, and the number of particles in clogging arches. For both particle types, the size of the opening, $D$, relative to the size of particles, $\ell$ is an important dimensionless measure. The orientation of the ellipses plays an important role in flow rheology and clogging. The alignment of contacting ellipses enhances the probability of forming stable arches. This study offers insight for applications involving the flow of granular materials consisting of ellipsoidal shapes, and possibly other non-spherical shapes.

Multiscale numerical modelling of crack propagation in two-dimensional metal plate

A novel multiscale model of brittle crack propagation in an Ag plate with macroscopic dimensions has been developed. The model represents crack propagation as stochastic drift-diffusion motion of the crack tip atom through the material, and couples the dynamics across three different length scales. It integrates the nanomechanics of bond rupture at the crack tip with the displacement and stress field equations of continuum based fracture theories. The finite element method is employed to obtain the continuum based displacement and stress fields over the macroscopic plate, and these are then used to drive the crack tip forward at the atomic level using the molecular dynamics simulation method based on many-body interatomic potentials. The linkage from the nanoscopic scale back to the macroscopic scale is established via the Ito stochastic calculus, the stochastic differential equation of which advances the tip to a new position on the macroscopic scale using the crack velocity and diffusion constant obtained on the nanoscale. Well known crack characteristics, such as the roughening transitions of the crack surfaces, crack velocity oscillations, as well as the macroscopic crack trajectories, are obtained.