Close-coupling calculations of elastic and inelastic positron-helium scattering

Elastic and inelastic positron-helium scattering have been investigated in different partial waves at medium energies using the close-coupling approximation with realistic wavefunctions employing the following states: He(1s1s), He(1s2s), He(1s2p), He(1s3s), He(1s3p), Ps(1s), Ps(2s) and Ps(2p). All excitations of the helium atom are in the spin-singlet electronic state. Calculations are reported of cross sections to He(1s1s), He(1s2s), and He(1s2p) transitions for incident positron energies up to 200 eV. These cross sections are in good agreement with experimental results.

Disturbance and repair of solitary waves in blood vessels with aneurysm

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Transmission problem for waves with frictional damping April 22, 2007

In this paper we consider the transmission problem, in one space dimension, for linear dissipative waves with frictional damping. We study the wave propagation in a medium with a component with attrition and another simply elastic. We show that for this type of material, the dissipation produced by the frictional part is strong enough to produce exponential decay of the solution, no matter how small is its size. ©2007 Texas State University.

Elastic Propagation in random media: applications to the imaging of volcano structures

High-frequency seismograms contain features that reflect the random inhomogeneities of the earth. In this work I use an imaging method to locate the high contrast small- scale heterogeneity respect to the background earth medium. This method was first introduced by Nishigami (1991) and than applied to different volcanic and tectonically active areas (Nishigami, 1997, Nishigami, 2000, Nishigami, 2006). The scattering imaging method is applied to two volcanic areas: Campi Flegrei and Mt. Vesuvius. Volcanic and seismological active areas are often characterized by complex velocity structures, due to the presence of rocks with different elastic properties. I introduce some modifications to the original method in order to make it suitable for small and highly complex media. In particular, for very complex media the single scattering approximation assumed by Nishigami (1991) is not applicable as the mean free path becomes short. The multiple scattering or diffusive approximation become closer to the reality. In this thesis, differently from the ordinary Nishigami’s method (Nishigami, 1991), I use the mean of the recorded coda envelope as reference curve and calculate the variations from this average envelope. In this way I implicitly do not assume any particular scattering regime for the "average" scattered radiation, whereas I consider the variations as due to waves that are singularly scattered from the strongest heterogeneities. The imaging method is applied to a relatively small area (20 x 20 km), this choice being justified by the small length of the analyzed codas of the low magnitude earthquakes. I apply the unmodified Nishigami’s method to the volcanic area of Campi Flegrei and compare the results with the other tomographies done in the same area. The scattering images, obtained with frequency waves around 18 Hz, show the presence of high scatterers in correspondence with the submerged caldera rim in the southern part of the Pozzuoli bay. Strong scattering is also found below the Solfatara crater, characterized by the presence of densely fractured, fluid-filled rocks and by a strong thermal anomaly. The modified Nishigami’s technique is applied to the Mt. Vesuvius area. Results show a low scattering area just below the central cone and a high scattering area around it. The high scattering zone seems to be due to the contrast between the high rigidity body located beneath the crater and the low rigidity materials located around it. The central low scattering area overlaps the hydrothermal reservoirs located below the central cone. An interpretation of the results in terms of geological properties of the medium is also supplied, aiming to find a correspondence of the scattering properties and the geological nature of the material. A complementary result reported in this thesis is that the strong heterogeneity of the volcanic medium create a phenomenon called "coda localization". It has been verified that the shape of the seismograms recorded from the stations located at the top of the volcanic edifice of Mt. Vesuvius is different from the shape of the seismograms recorded at the bottom. This behavior is justified by the consideration that the coda energy is not uniformly distributed within a region surrounding the source for great lapse time.

Topological spin waves in the atomic-scale magnetic skyrmion crystal

We study the spin waves of the triangular skyrmion crystal that emerges in a two-dimensional spin lattice model as a result of the competition between Heisenberg exchange, Dzyalonshinkii–Moriya interactions, Zeeman coupling and uniaxial anisotropy. The calculated spin wave bands have a finite Berry curvature that, in some cases, leads to non-zero Chern numbers, making this system topologically distinct from conventional magnonic systems. We compute the edge spin-waves, expected from the bulk-boundary correspondence principle, and show that they are chiral, which makes them immune to elastic backscattering. Our results illustrate how topological phases can occur in self-generated emergent superlattices at the mesoscale.

Sub-wavelength resonators as seismic Rayleigh waves shield

We explore the thesis that tall structures can be protected by means of seismic metamaterials. Seismic metamaterials can be built as some elements are created over soil layer with different shapes, dimensions, patterns and from different materials. Resonances in these elements are acting as locally resonant metamaterials for Rayleigh surface waves in the geophysics context. Analytically we proved that if we put infinite chain of SDOF resonator over the soil layer as an elastic, homogeneous and isotropic material, vertical component of Rayleigh wave, longitudinal resonance of oscillators will couple with each other, they would create a Rayleigh bandgap frequency, and wave will experience attenuation before it reaches the structure. As it is impossible to use infinite chain of resonators over soil layer, we considered finite number of resonators throughout our simulations. Analytical work is interpreted using finite element simulations that demonstrates the observed attenuation is due to bandgaps when oscillators are arranged at sub-wavelength scale with respect to the incident Rayleigh wave. For wavelength less than 5 meters, the resulting bandgaps are remarkably large and strongly attenuating when impedance of oscillators matches impedance of soil. Since longitudinal resonance of SDOF resonator are proportional to its length inversely, a formed array of resonators that attenuates Rayleigh waves at frequency ≤10 Hz could be designed starting from vertical pillars coupled to the ground. Optimum number of vertical pillars and their interval spacing called effective area of resonators are investigated. For 10 pillars with effective area of 1 meter and resonance frequency of 4.9 Hz, bandgap frequency causes attenuation and a sinusoidal impulsive force illustrate wave steering down phenomena. Simulation results proved analytical findings of this work.

Small strain elastic behaviour of unsaturated soil investigated by bender/extender element testing

The aim of this project was to investigate very small strain elastic behaviour of soils under unsaturated conditions, using bender/extender element (BEE) testing. The behaviour of soils at very small strains has been widely studied under saturated conditions, whereas much less work has been performed on very small strain behaviour under unsaturated conditions. A suction-controlled double wall triaxial apparatus for unsaturated soil testing was modified to incorporate three pairs of BEEs transmitting both shear and compression waves with vertical and horizontal directions of wave transmission and wave polarisation. Various different techniques for measuring wave travel time were investigated in both the time domain and the frequency domain and it was concluded that, at least for the current experimental testing programme, peak-to-first-peak in the time domain was the most reliable technique for determining wave travel time. An experimental test programme was performed on samples of compacted speswhite kaolin clay. Two different forms of compaction were employed (i.e. isotropic and anisotropic). Compacted kaolin soil samples were subjected to constant suction loading and unloading stages at three different values of suction, covering both unsaturated conditions (s= 50kPa and s= 300kPa) and saturated conditions (s=0). Loading and unloading stages were performed at three different values of stress ratio (η=0, η=1 and η=-1 ). In some tests a wetting-drying cycle was performed before or within the loading stage, with the wetting-drying cycles including both wetting-induced swelling and wetting-induced collapse compression. BEE tests were performed at regular intervals throughout all test stages, to measure shear wave velocity Vs and compression wave velocity Vp and hence to determine values of shear modulus G and constrained modulus M. The experimental test programme was designed to investigate how very small strain shear modulus G and constrained modulus M varied with unsaturated state variables, including how anisotropy of these parameters developed either with stress state (stress-induced anisotropy) or with previous straining (strain-induced anisotropy). A new expression has been proposed for the very small strain shear modulus G of an isotropic soil under saturated and unsaturated conditions. This expression relates the variation of G to only mean Bishop’s stress p* and specific volume v, and it converges to a well-established expression for saturated soils as degree of saturation approaches 1. The proposed expression for G is able to predict the variation of G under saturated and unsaturated conditions at least as well as existing expressions from the literature and it is considerably simpler (employing fewer state variables and fewer soil constants). In addition, unlike existing expressions from the literature, the values of soil constants in the proposed new expression can be determined from a saturated test. It appeared that, in the current project at least, any strain-induced anisotropy of very small strain elastic behaviour was relatively modest, with the possible exception of loading in triaxial extension. It was therefore difficult to draw any firm conclusion about evolution of strain-induced anisotropy and whether it depended upon the same aspects of soil fabric as evolution of anisotropy of large strain plastic behaviour. Stress-induced anisotropy of very small strain elastic behaviour was apparent in the experimental test programme. An attempt was made to extend the proposed expression for G to include the effect of stress-induced anisotropy. Interpretation of the experimental results indicated that the value of shear modulus was affected by the values of all three principal Bishop’s stresses (in the direction of wave transmission, the direction of wave polarisation and the third mutually perpendicular direction). However, prediction of stress-induced anisotropy was only partially successful, and it was concluded that the effect of Lode angle was also significant.

Dynamics of Granular Crystals with Elastic-Plastic Contacts

We study the behavior of granular crystals subjected to impact loading that creates plastic deformation at the contacts between constituent particles. Granular crystals are highly periodic arrangements of spherical particles, arranged into densely packed structures resembling crystals. This special class of granular materials has been shown to have unique dynamics with suggested applications in impact protection. However, previous work has focused on very low amplitude impacts where every contact point can be described using the Hertzian contact law, valid only for purely elastic deformation. In this thesis, we extend previous investigation of the dynamics of granular crystals to significantly higher impact energies more suitable for the majority of applications. Additionally, we demonstrate new properties specific to elastic-plastic granular crystals and discuss their potential applications as well. We first develop a new contact law to describe the interaction between particles for large amplitude compression of elastic-plastic spherical particles including a formulation for strain-rate dependent plasticity. We numerically and experimentally demonstrate the applicability of this contact law to a variety of materials typically used in granular crystals. We then extend our investigation to one-dimensional chains of elastic-plastic particles, including chains of alternating dissimilar materials. We show that, using the new elastic-plastic contact law, we can predict the speed at which impact waves with plastic dissipation propagate based on the material properties of the constituent particles. Finally, we experimentally and numerically investigate the dynamics of two-dimensional and three-dimensional granular crystals with elastic-plastic contacts. We first show that the predicted wave speeds for 1D granular crystals can be extended to 2D and 3D materials. We then investigate the behavior of waves propagating across oblique interfaces of dissimilar particles. We show that the character of the refracted wave can be predicted using an analog to Snell's law for elastic-plastic granular crystals and ultimately show how it can be used to design impact guiding "lenses" for mitigation applications.

Resonant metamaterials for the control of Rayleigh waves

Since their emergence, locally resonant metamaterials have found several applications for the control of surface waves, from micrometer-sized electronic devices to meter-sized seismic barriers. The interaction between Rayleigh-type surface waves and resonant metamaterials has been investigated through the realization of locally resonant metasurfaces, thin elastic interfaces constituted by a cluster of resonant inclusions or oscillators embedded near the surface of an elastic waveguide. When such resonant metasurfaces are embedded in an elastic homogeneous half-space, they can filter out the propagation of Rayleigh waves, creating low-frequency bandgaps at selected frequencies. In the civil engineering context, heavy resonating masses are needed to extend the bandgap frequency width of locally resonant devices, a requirement that limits their practical implementations. In this dissertation, the wave attenuation capabilities of locally resonant metasurfaces have been enriched by proposing (i) tunable metasurfaces to open large frequency bandgaps with small effective inertia, and by developing (ii) an analytical framework aimed at studying the propagation of Rayleigh waves propagation in deep resonant waveguides. In more detail, inertial amplified resonators are exploited to design advanced metasurfaces with a prescribed static and a tunable dynamic response. The modular design of the tunable metasurfaces allows to shift and enlarge low-frequency spectral bandgaps without modifying the total inertia of the metasurface. Besides, an original dispersion law is derived to study the dispersive properties of Rayleigh waves propagating in thick resonant layers made of sub-wavelength resonators. Accordingly, a deep resonant wave barrier of mechanical resonators embedded inside the soil is designed to impede the propagation of seismic surface waves. Numerical models are developed to confirm the analytical dispersion predictions of the tunable metasurface and resonant layer. Finally, a medium-size scale resonant wave barrier is designed according to the soil stratigraphy of a real geophysical scenario to attenuate ground-borne vibration.

Analytical methods for modeling elastic metasurfaces

This dissertation aims at developing advanced analytical tools able to model surface waves propagating in elastic metasurfaces. In particular, four different objectives are defined and pursued throughout this work to enrich the description of the metasurface dynamics. First, a theoretical framework is developed to describe the dispersion properties of a seismic metasurface composed of discrete resonators placed on a porous medium considering part of it fully saturated. Such a model combines classical elasticity theory, Biot’s poroelasticity and an effective medium approach to describe the metasurface dynamics and its coupling with the poroelastic substrate. Second, an exact formulation based on the multiple scattering theory is developed to extend the two-dimensional classical Lamb’s problem to the case of an elastic half-space coupled to an arbitrary number of discrete surface resonators. To this purpose, the incident wavefield generated by a harmonic source and the scattered field generated by each resonator are calculated. The substrate wavefield is then obtained as solutions of the coupled problem due to the interference of the incident field and the multiple scattered fields of the oscillators. Third, the above discussed formulation is extended to three-dimensional contexts. The purpose here is to investigate the dynamic behavior and the topological properties of quasiperiodic elastic metasurfaces. Finally, the multiple scattering formulation is extended to model flexural metasurfaces, i.e., an array of thin plates. To this end, the resonant plates are modeled by means of their equivalent impedance, derived by exploiting the Kirchhoff plate theory. The proposed formulation permits the treatment of a general flexural metasurface, with no limitation on the number of plates and the configuration taken into account. Overall, the proposed analytical tools could pave the way for a better understanding of metasurface dynamics and their implementation in engineered devices.

Analysis of the tip roundness effects on the micro- and macroindentation response of elastic-plastic materials

In this work, the effects of indenter tip roundness oil the load-depth indentation curves were analyzed using finite element modeling. The tip roundness level was Studied based on the ratio between tip radius and maximum penetration depth (R/h(max)), which varied from 0.02 to 1. The proportional Curvature constant (C), the exponent of depth during loading (alpha), the initial unloading slope (S), the correction factor (beta), the level of piling-up or sinking-in (h(c)/h(max)), and the ratio h(max)/h(f) are shown to be strongly influenced by the ratio R/h(max). The hardness (H) was found to be independent of R/h(max) in the range studied. The Oliver and Pharr method was successful in following the variation of h(c)/h(max) with the ratio R/h(max) through the variation of S with the ratio R/h(max). However, this work confirmed the differences between the hardness values calculated using the Oliver-Pharr method and those obtained directly from finite element calculations; differences which derive from the error in area calculation that Occurs when given combinations of indented material properties are present. The ratio of plastic work to total work (W(p)/W(t)) was found to be independent of the ratio R/h(max), which demonstrates that the methods for the Calculation of mechanical properties based on the *indentation energy are potentially not Susceptible to errors caused by tip roundness.

Analysis of the effects of conical indentation variables on the indentation response of elastic-plastic materials through factorial design of experiment

In this work, the effects of conical indentation variables on the load-depth indentation curves were analyzed using finite element modeling and dimensional analysis. A factorial design 2(6) was used with the aim of quantifying the effects of the mechanical properties of the indented material and of the indenter geometry. Analysis was based on the input variables Y/E, R/h(max), n, theta, E, and h(max). The dimensional variables E and h(max) were used such that each value of dimensionless Y/E was obtained with two different values of E and each value of dimensionless R/h(max) was obtained with two different h(max) values. A set of dimensionless functions was defined to analyze the effect of the input variables: Pi(1) = P(1)/Eh(2), Pi(2) = h(c)/h, Pi(3) = H/Y, Pi(4) = S/Eh(max), Pi(6) = h(max)/h(f) and Pi(7) = W(P)/W(T). These six functions were found to depend only on the dimensionless variables studied (Y/E, R/h(max), n, theta). Another dimension less function, Pi(5) = beta, was not well defined for most of the dimensionless variables and the only variable that provided a significant effect on beta was theta. However, beta showed a strong dependence on the fraction of the data selected to fit the unloading curve, which means that beta is especially Susceptible to the error in the Calculation of the initial unloading slope.

Using the ultrasound and instrumented indentation techniques to measure the elastic modulus of engineering materials

Currently, the acoustic and nanoindentation techniques are two of the most used techniques for material elastic modulus measurement. In this article fundamental principles and limitations of both techniques are shown and discussed. Last advances in nanoindentation technique are also reviewed. An experimental study in ceramic, metallic, composite and single crystals was also done. Results shown that ultrasonic technique is capable to provide results in agreement with those reported in literature. However, ultrasonic technique does not allow measuring the elastic modulus of some small samples and single crystals. On the other hand, the nanoindentation technique estimates the elastic modulus values in reasonable agreement with those measured by acoustic methods, particularly in amorphous materials, while in some policristaline materials some deviation from expected values was obtained.

Quasiharmonic elastic constants corrected for deviatoric thermal stresses

The quasiharmonic approximation (QHA), in its simplest form also called the statically constrained (SC) QHA, has been shown to be a straightforward method to compute thermoelastic properties of crystals. Recently we showed that for noncubic solids SC-QHA calculations develop deviatoric thermal stresses at high temperatures. Relaxation of these stresses leads to a series of corrections to the free energy that may be taken to any desired order, up to self-consistency. Here we show how to correct the elastic constants obtained using the SC-QHA. We exemplify the procedure by correcting to first order the elastic constants of MgSiO(3) perovskite and MgSiO(3) postperovskite, the major phases of the Earth's lower mantle. We show that this first-order correction is quite satisfactory for obtaining the aggregated elastic averages of these minerals and their velocities in the lower mantle. This type of correction is also shown to be applicable to experimental measurements of elastic constants in situations where deviatoric stresses can develop, such as in diamond-anvil cells.

Reduction of chaotic particle transport driven by drift waves in sheared flows

Investigations of chaotic particle transport by drift waves propagating in the edge plasma of tokamaks with poloidal zonal flow are described. For large aspect ratio tokamaks, the influence of radial electric field profiles on convective cells and transport barriers, created by the nonlinear interaction between the poloidal flow and resonant waves, is investigated. For equilibria with edge shear flow, particle transport is seen to be reduced when the electric field shear is reversed. The transport reduction is attributed to the robust invariant tori that occur in nontwist Hamiltonian systems. This mechanism is proposed as an explanation for the transport reduction in Tokamak Chauffage Alfven Bresilien [R. M. O. Galvao , Plasma Phys. Controlled Fusion 43, 1181 (2001)] for discharges with a biased electrode at the plasma edge.