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.

Nanomechanical and phononic properties of structured soft materials

Significant interest in nanotechnology, is stimulated by the fact that materials exhibit qualitative changes of properties when their dimensions approach ”finite-sizes”. Quantization of electronic, optical and acoustic energies at the nanoscale provides novel functions, with interests spanning from electronics and photonics to biology. The present dissertation involves the application of Brillouin light scattering (BLS) to quantify and utilize material displacementsrnfor probing phononics and elastic properties of structured systems with dimensions comparable to the wavelength of visible light. The interplay of wave propagation with materials exhibiting spatial inhomogeneities at sub-micron length scales provides information not only about elastic properties but also about structural organization at those length scales. In addition the vector nature of q allows, for addressing the directional dependence of thermomechanical properties. To meet this goal, one-dimensional confined nanostructures and a biological system possessing high hierarchical organization were investigated. These applications extend the capabilities of BLS from a characterization tool for thin films to a method for unravelingrnintriguing phononic properties in more complex systems.

Computational prediction of the influence of crosslink distribution on the thermo-mechanical properties of epoxies

Experimental studies on epoxies report that the microstructure consists of highly crosslinked localized regions connected with a dispersed phase of low crosslink density. The various thermo-mechanical properties of epoxies might be affected by the crosslink distribution. But as experiments cannot report the exact number of crosslinked covalent bonds present in the structure, molecular dynamics is thus being used in this work to determine the influence of crosslink distribution on thermo-mechanical properties. Molecular dynamics and molecular mechanics simulations are used to establish wellequilibrated molecular models of EPON 862-DETDA epoxy system with a range of crosslink densities and various crosslink distributions. Crosslink distributions are being varied by forming differently crosslinked localized clusters and then by forming different number of crosslinks interconnecting the clusters. Simulations are subsequently used to predict the volume shrinkage, thermal expansion coefficients, and elastic properties of each of the crosslinked systems. The results indicate that elastic properties increase with increasing levels of overall crosslink density and the thermal expansion coefficient decreases with overall crosslink density, both above and below the glass transition temperature. Elastic moduli and coefficients of linear thermal expansion values were found to be different for systems with same overall crosslink density but having different crosslink distributions, thus indicating an effect of the epoxy nanostructure on physical properties. The values of thermo-mechanical properties for all the crosslinked systems are within the range of values reported in literature.

Predicting the Mechanical Properties of Carbon-Based Materials Using Molecular Dynamics

The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although molecular modeling has been used extensively to predict elastic properties of materials, modeling of more complex phenomenon such as fracture has only recently been possible with the development of new force fields such as ReaxFF, which is used in this work. It is not fully understood what molecular modeling parameters such as thermostat type, thermostat coupling, time step, system size, and strain rate are required for accurate modeling of fracture. Selection of modeling parameters to model fracture can be difficult and non-intuitive compared to modeling elastic properties using traditional force fields, and the errors generated by incorrect parameters may be non-obvious. These molecular modeling parameters are systematically investigated and their effects on the fracture of well-known carbon materials are analyzed. It is determined that for coupling coefficients of 250 fs and greater do not result in substantial differences in the stress-strain response of the materials using any thermostat type. A time step of 0.5 fs of smaller is required for accurate results. Strain rates greater than 2.2 ns-1 are sufficient to obtain repeatable results with slower strain rates for the materials studied. The results of this study indicate that further refinement of the Chenoweth parameter set is required to accurately predict the mechanical response of carbon-based systems. The ReaxFF has been used extensively to model systems in which bond breaking and formation occur. In particular ReaxFF has been used to model reactions of small molecules. Some elastic and fracture properties have been successfully modeled using ReaxFF in materials such as silicon and some metals. However, it is not clear if current parameterizations for ReaxFF are able to accurately reproduce the elastic and fracture properties of carbon materials. The stress-strain response of a new ReaxFF parameterization is compared to the previous parameterization and density functional theory results for well-known carbon materials. The new ReaxFF parameterization makes xv substantial improvements to the predicted mechanical response of carbon materials, and is found to be suitable for modeling the mechanical response of carbon materials. Finally, a new material composed of carbon nanotubes within an amorphous carbon (AC) matrix is modeled using the ReaxFF. Various parameters that may be experimentally controlled are investigated such as nanotube bundling, comparing multi-walled nanotube with single-walled nanotubes, and degree of functionalization of the nanotubes. Elastic and fracture properties are investigated for the composite systems and compared to results of pure-nanotube and pure-AC models. It is found that the arrangement of the nanotubes and degree of crosslinking may substantially affect the properties of the systems, particularly in the transverse directions.

An anisotropic elastic-viscoplastic damage model for bone tissue

A new anisotropic elastic-viscoplastic damage constitutive model for bone is proposed using an eccentric elliptical yield criterion and nonlinear isotropic hardening. A micromechanics-based multiscale homogenization scheme proposed by Reisinger et al. is used to obtain the effective elastic properties of lamellar bone. The dissipative process in bone is modeled as viscoplastic deformation coupled to damage. The model is based on an orthotropic ecuntric elliptical criterion in stress space. In order to simplify material identification, an eccentric elliptical isotropic yield surface was defined in strain space, which is transformed to a stress-based criterion by means of the damaged compliance tensor. Viscoplasticity is implemented by means of the continuous Perzyna formulation. Damage is modeled by a scalar function of the accumulated plastic strain D(κ) , reducing all element s of the stiffness matrix. A polynomial flow rule is proposed in order to capture the rate-dependent post-yield behavior of lamellar bone. A numerical algorithm to perform the back projection on the rate-dependent yield surface has been developed and implemented in the commercial finite element solver Abaqus/Standard as a user subroutine UMAT. A consistent tangent operator has been derived and implemented in order to ensure quadratic convergence. Correct implementation of the algorithm, convergence, and accuracy of the tangent operator was tested by means of strain- and stress-based single element tests. A finite element simulation of nano- indentation in lamellar bone was finally performed in order to show the abilities of the newly developed constitutive model.

How fluids bend: the elastic expansion for higher-dimensional black holes

Hydrodynamics can be consistently formulated on surfaces of arbitrary co-dimension in a background space-time, providing the effective theory describing long-wavelength perturbations of black branes. When the co-dimension is non-zero, the system acquires fluid-elastic properties and constitutes what is called a fluid brane. Applying an effective action approach, the most general form of the free energy quadratic in the extrinsic curvature and extrinsic twist potential of stationary fluid brane configurations is constructed to second order in a derivative expansion. This construction generalizes the Helfrich-Canham bending energy for fluid membranes studied in theoretical biology to the case in which the fluid is rotating. It is found that stationary fluid brane configurations are characterized by a set of 3 elastic response coefficients, 3 hydrodynamic response coefficients and 1 spin response coefficient for co-dimension greater than one. Moreover, the elastic degrees of freedom present in the system are coupled to the hydrodynamic degrees of freedom. For co-dimension-1 surfaces we find a 8 independent parameter family of stationary fluid branes. It is further shown that elastic and spin corrections to (non)-extremal brane effective actions can be accounted for by a multipole expansion of the stress-energy tensor, therefore establishing a relation between the different formalisms of Carter, Capovilla-Guven and Vasilic-Vojinovic and between gravity and the effective description of stationary fluid branes. Finally, it is shown that the Young modulus found in the literature for black branes falls into the class predicted by this approach - a relation which is then used to make a proposal for the second order effective action of stationary blackfolds and to find the corrected horizon angular velocity of thin black rings.

Experimental, theoretical and numerical investigation of the nonlinear micromechanical properties of bone

Aging societies suffer from an increasing incidence of bone fractures. Bone strength depends on the amount of mineral measured by clinical densitometry, but also on the micromechanical properties of the bone hierarchical organization. A good understanding has been reached for elastic properties on several length scales, but up to now there is a lack of reliable postyield data on the lower length scales. In order to be able to describe the behavior of bone at the microscale, an anisotropic elastic-viscoplastic damage model was developed using an eccentric generalized Hill criterion and nonlinear isotropic hardening. The model was implemented as a user subroutine in Abaqus and verified using single element tests. A FE simulation of microindentation in lamellar bone was finally performed show-ing that the new constitutive model can capture the main characteristics of the indentation response of bone. As the generalized Hill criterion is limited to elliptical and cylindrical yield surfaces and the correct shape for bone is not known, a new yield surface was developed that takes any convex quadratic shape. The main advantage is that in the case of material identification the shape of the yield surface does not have to be anticipated but a minimization results in the optimal shape among all convex quadrics. The generality of the formulation was demonstrated by showing its degeneration to classical yield surfaces. Also, existing yield criteria for bone at multiple length scales were converted to the quadric formulation. Then, a computational study to determine the influence of yield surface shape and damage on the in-dentation response of bone using spherical and conical tips was performed. The constitutive model was adapted to the quadric criterion and yield surface shape and critical damage were varied. They were shown to have a major impact on the indentation curves. Their influence on indentation modulus, hardness, their ratio as well as the elastic to total work ratio were found to be very well described by multilinear regressions for both tip shapes. For conical tips, indentation depth was not a significant fac-tor, while for spherical tips damage was insignificant. All inverse methods based on microindentation suffer from a lack of uniqueness of the found material properties in the case of nonlinear material behavior. Therefore, monotonic and cyclic micropillar com-pression tests in a scanning electron microscope allowing a straightforward interpretation comple-mented by microindentation and macroscopic uniaxial compression tests were performed on dry ovine bone to identify modulus, yield stress, plastic deformation, damage accumulation and failure mecha-nisms. While the elastic properties were highly consistent, the postyield deformation and failure mech-anisms differed between the two length scales. A majority of the micropillars showed a ductile behavior with strain hardening until failure by localization in a slip plane, while the macroscopic samples failed in a quasi-brittle fashion with microcracks coalescing into macroscopic failure surfaces. In agreement with a proposed rheological model, these experiments illustrate a transition from a ductile mechanical behavior of bone at the microscale to a quasi-brittle response driven by the growth of preexisting cracks along interfaces or in the vicinity of pores at the macroscale. Subsequently, a study was undertaken to quantify the topological variability of indentations in bone and examine its relationship with mechanical properties. Indentations were performed in dry human and ovine bone in axial and transverse directions and their topography measured by AFM. Statistical shape modeling of the residual imprint allowed to define a mean shape and describe the variability with 21 principal components related to imprint depth, surface curvature and roughness. The indentation profile of bone was highly consistent and free of any pile up. A few of the topological parameters, in particular depth, showed significant correlations to variations in mechanical properties, but the cor-relations were not very strong or consistent. We could thus verify that bone is rather homogeneous in its micromechanical properties and that indentation results are not strongly influenced by small de-viations from the ideal case. As the uniaxial properties measured by micropillar compression are in conflict with the current literature on bone indentation, another dissipative mechanism has to be present. The elastic-viscoplastic damage model was therefore extended to viscoelasticity. The viscoelastic properties were identified from macroscopic experiments, while the quasistatic postelastic properties were extracted from micropillar data. It was found that viscoelasticity governed by macroscale properties has very little influence on the indentation curve and results in a clear underestimation of the creep deformation. Adding viscoplasticity leads to increased creep, but hardness is still highly overestimated. It was possible to obtain a reasonable fit with experimental indentation curves for both Berkovich and spherical indenta-tion when abandoning the assumption of shear strength being governed by an isotropy condition. These results remain to be verified by independent tests probing the micromechanical strength prop-erties in tension and shear. In conclusion, in this thesis several tools were developed to describe the complex behavior of bone on the microscale and experiments were performed to identify its material properties. Micropillar com-pression highlighted a size effect in bone due to the presence of preexisting cracks and pores or inter-faces like cement lines. It was possible to get a reasonable fit between experimental indentation curves using different tips and simulations using the constitutive model and uniaxial properties measured by micropillar compression. Additional experimental work is necessary to identify the exact nature of the size effect and the mechanical role of interfaces in bone. Deciphering the micromechanical behavior of lamellar bone and its evolution with age, disease and treatment and its failure mechanisms on several length scales will help preventing fractures in the elderly in the future.

Measuring the mechanical properties of plant cells by combining micro-indentation with osmotic treatments

Growth in plants results from the interaction between genetic and signalling networks and the mechanical properties of cells and tissues. There has been a recent resurgence in research directed at understanding the mechanical aspects of growth, and their feedback on genetic regulation. This has been driven in part by the development of new micro-indentation techniques to measure the mechanical properties of plant cells in vivo. However, the interpretation of indentation experiments remains a challenge, since the force measures results from a combination of turgor pressure, cell wall stiffness, and cell and indenter geometry. In order to interpret the measurements, an accurate mechanical model of the experiment is required. Here, we used a plant cell system with a simple geometry, Nicotiana tabacum Bright Yellow-2 (BY-2) cells, to examine the sensitivity of micro-indentation to a variety of mechanical and experimental parameters. Using a finite-element mechanical model, we found that, for indentations of a few microns on turgid cells, the measurements were mostly sensitive to turgor pressure and the radius of the cell, and not to the exact indenter shape or elastic properties of the cell wall. By complementing indentation experiments with osmotic experiments to measure the elastic strain in turgid cells, we could fit the model to both turgor pressure and cell wall elasticity. This allowed us to interpret apparent stiffness values in terms of meaningful physical parameters that are relevant for morphogenesis.

Seismic properties of volcanic rocks from DSDP Hole 62-465A

Hess Rise, a major structural feature in the northern Pacific, is one of several oceanic plateaus which apparently possess anomalous seismic structures (Sutton et al., 1971). Hence, Laboratory measurements of compressional- and shear-wave velocities in rocks from oceanic plateau regions are of considerable interest. Several questions come to mind: (1) Are compressional - wave velocities of volcanic rocks from oceanic plateaus similar to basalts of equivalent age from normal oceanic crust? (2) Do velocity-density relations for plateau rocks fit the well-established trends for Layer 2 basalts? (3) How do Poisson's ratios, determined from compressional- and shear-wave velocities, of oceanic plateau rocks compare with those of normal sea-floor basalts? To answer these questions, we have selected for velocity measurements five volcanic samples from different depths from Hole 465A, on southern Hess Rise. It is particularly significant that all the rocks are vesicular and have been highly altered, which significantly affects their elastic properties.

Compressional and shear wave velocities and elastic constants at DSDP Leg 59 Holes

Prior to the Deep Sea Drilling Project the composition of the oceanic crust could only be inferred from seismic-refraction and gravity data and the recovery of a wide variety of dredged rocks. Through the success of the Deep Sea Drilling Project, it is now clear that the top of oceanic Layer 2 usually consists of basalt. Several laboratory studies (e.g., Fox et al., 1972; Christensen and Shaw, 1970; Hyndman and Drury, 1976) have demonstrated that the seismic velocities of oceanic basalt are similar to velocities reported from refraction studies of Layer 2 and that the variability in Layer 2 velocities has many causes, the most important being fracturing and sea-floor alteration produced by the interaction of basalt and sea water (Christensen and Salisbury, 1973). To date, most reported measurements of velocities in oceanic basalts are from samples obtained from the main ocean basins. With the exception of an earlier study of velocities and related elastic properties of a suite of rocks from DSDP Sites 292, 293, 294, and 296 located in the Philippine Sea (Christensen et al., 1975; Fountain et al., 1975), elastic properties have not been determined for oceanic rocks from marginal basins. In this chapter compressional- and shear-wave velocities and elastic constants are reported at elevated confining pressures for basalt and volcanic breccias from Holes 447A, 448, and 448A.

Index properties, consolidation characteristics and sediment rebound of ODP Sites 167-1020 and 167-1021 sediments

Drilling during Leg 167 at the California margin was scheduled to recover continuous sedimentary sections. Multiple advanced piston core (APC) holes drilled at different depth offsets provided core overlap in successive APCs. Correlation of high-resolution laboratory physical properties data from adjacent APC holes was used to compile composite depth sections for each site. The composite depth sections were used to confirm continuous recovery and enable high-resolution sampling. The meters composite depth (mcd) scale differs from the shipboard meters below seafloor (mbsf) scale because of (1) core expansion following recovery (MacKillop et al., 1995, doi:10.2973/odp.proc.sr.138.118.1995), (2) coring gaps, and (3) stretching/compression of sediment during coring (Lyle, Koizumi, Richter, et al., 1997, doi:10.2973/odp.proc.ir.167.1997). Moran (1997, doi:10.2973/odp.proc.sr.154.132.1997) calculated that sediment expansion accounted for 90%-95% of the Leg 154 depth offset between shipboard mbsf and the mcd scales. Terzaghi's one-dimensional theory of consolidation (Terzaghi, 1943) describes the response of sediments to stress loading and release. Mechanical loading in marine environments is provided by the buoyant weight of the overlying sediments. The load increases with depth below seabed, resulting in sediment volume reduction as water is "squeezed" out of the voids in the sediment. Stress release during core recovery results in expansion of the sediment and volume increase as water returns to the sediment. The sediment expansion or rebound defines the elastic properties of the sediment. In this study we examine the elastic deformation properties of sediments recovered from Sites 1020 and 1021. These results are used to (1) correct the laboratory index properties measurements to in situ values and (2) determine the contribution of sediment rebound to the depth offset between the mbsf and mcd scales.