Three-dimensional finite element analysis of composites with coated spherical inclusions

A three-dimensional finite element analysis has been used to determine the internal stresses in a three-phase composite. The stresses have been determined for a variety of interphase properties, the thicknesses of the interphase and the volume fractions of particles. Young's modulus has been calculated from a knowledge of these stresses and the applied deformation. The calculations show that stress distributions in the matrix and the mechanical properties are sensitive to the interphase property in the three-phase composites. The interfacial stresses in the three-dimensional analysis are in agreement with results obtained by an axisymmetric analysis. The predicted bulk modulus in three-dimensional analysis agrees well with the theoretical solution obtained by Qui and Weng, but it presents a great divergence from that in axisymmetric analyses. An investigation indicates that this divergence may be caused by the difference in the unit cell structure between two models. A comparison of the numerically predicted bulk and shear modulus for two-phase composites with the theoretical results indicates that the three-dimensional analysis gives quite satisfactory results.

Methods of multiparameter inversion of seismic data using the acoustic and elastic born approximations

This thesis presents two different forms of the Born approximations for acoustic and elastic wavefields and discusses their application to the inversion of seismic data. The Born approximation is valid for small amplitude heterogeneities superimposed over a slowly varying background. The first method is related to frequency-wavenumber migration methods. It is shown to properly recover two independent acoustic parameters within the bandpass of the source time function of the experiment for contrasts of about 5 percent from data generated using an exact theory for flat interfaces. The independent determination of two parameters is shown to depend on the angle coverage of the medium. For surface data, the impedance profile is well recovered.

The second method explored is mathematically similar to iterative tomographic methods recently introduced in the geophysical literature. Its basis is an integral relation between the scattered wavefield and the medium parameters obtained after applying a far-field approximation to the first-order Born approximation. The Davidon-Fletcher-Powell algorithm is used since it converges faster than the steepest descent method. It consists essentially of successive backprojections of the recorded wavefield, with angular and propagation weighing coefficients for density and bulk modulus. After each backprojection, the forward problem is computed and the residual evaluated. Each backprojection is similar to a before-stack Kirchhoff migration and is therefore readily applicable to seismic data. Several examples of reconstruction for simple point scatterer models are performed. Recovery of the amplitudes of the anomalies are improved with successive iterations. Iterations also improve the sharpness of the images.

The elastic Born approximation, with the addition of a far-field approximation is shown to correspond physically to a sum of WKBJ-asymptotic scattered rays. Four types of scattered rays enter in the sum, corresponding to P-P, P-S, S-P and S-S pairs of incident-scattered rays. Incident rays propagate in the background medium, interacting only once with the scatterers. Scattered rays propagate as if in the background medium, with no interaction with the scatterers. An example of P-wave impedance inversion is performed on a VSP data set consisting of three offsets recorded in two wells.

High pressure states in condensed matter: I. High pressure behavior of the iron-sulphur system with applications to the earth's core. II. Empirical equation of state for organic compounds at high pressures

Part I:

The earth's core is generally accepted to be composed primarily of iron, with an admixture of other elements. Because the outer core is observed not to transmit shear waves at seismic frequencies, it is known to be liquid or primarily liquid. A new equation of state is presented for liquid iron, in the form of parameters for the 4th order Birch-Murnaghan and Mie-Grüneisen equations of state. The parameters were constrained by a set of values for numerous properties compiled from the literature. A detailed theoretical model is used to constrain the P-T behavior of the heat capacity, based on recent advances in the understanding of the interatomic potentials for transition metals. At the reference pressure of 10⁵ Pa and temperature of 1811 K (the normal melting point of Fe), the parameters are: ρ = 7037 kg/m³, K_S0 = 110 GPa, K_S' = 4.53, K_S" = -.0337 GPa-1, and γ = 2.8, with γ α ρ^-1.17. Comparison of the properties predicted by this model with the earth model PREM indicates that the outer core is 8 to 10 % less dense than pure liquid Fe at the same conditions. The inner core is also found to be 3 to 5% less dense than pure liquid Fe, supporting the idea of a partially molten inner core. The density deficit of the outer core implies that the elements dissolved in the liquid Fe are predominantly of lower atomic weight than Fe. Of the candidate light elements favored by researchers, only sulfur readily dissolves into Fe at low pressure, which means that this element was almost certainly concentrated in the core at early times. New melting data are presented for FeS and FeS₂ which indicate that the FeS₂ is the S-hearing liquidus solid phase at inner core pressures. Consideration of the requirement that the inner core boundary be observable by seismological means and the freezing behavior of solutions leads to the possibility that the outer core may contain a significant fraction of solid material. It is found that convection in the outer core is not hindered if the solid particles are entrained in the fluid flow. This model for a core of Fe and S admits temperatures in the range 3450K to 4200K at the top of the core. An all liquid Fe-S outer core would require a temperature of about 4900 K at the top of the core.

Part II.

The abundance of uses for organic compounds in the modern world results in many applications in which these materials are subjected to high pressures. This leads to the desire to be able to describe the behavior of these materials under such conditions. Unfortunately, the number of compounds is much greater than the number of experimental data available for many of the important properties. In the past, one approach that has worked well is the calculation of appropriate properties by summing the contributions from the organic functional groups making up molecules of the compounds in question. A new set of group contributions for the molar volume, volume thermal expansivity, heat capacity, and the Rao function is presented for functional groups containing C, H, and O. This set is, in most cases, limited in application to low molecular liquids. A new technique for the calculation of the pressure derivative of the bulk modulus is also presented. Comparison with data indicates that the presented technique works very well for most low molecular hydrocarbon liquids and somewhat less well for oxygen-bearing compounds. A similar comparison of previous results for polymers indicates that the existing tabulations of group contributions for this class of materials is in need of revision. There is also evidence that the Rao function contributions for polymers and low molecular compounds are somewhat different.

The Equation of State and Petrogenesis of Komatiite

(1) Equation of State of Komatiite

The equation of state (EOS) of a molten komatiite (27 wt% MgO) was detennined in the 5 to 36 GPa pressure range via shock wave compression from 1550°C and 0 bar. Shock wave velocity, U_S, and particle velocity, U_P, in km/s follow the linear relationship U_S = 3.13(±0.03) + 1.47(±0.03) U_P. Based on a calculated density at 1550°C, 0 bar of 2.745±0.005 glee, this U_S-U_P relationship gives the isentropic bulk modulus K_S = 27.0 ± 0.6 GPa, and its first and second isentropic pressure derivatives, K'_S = 4.9 ± 0.1 and K"_S = -0.109 ± 0.003 GPa^-1.

The calculated liquidus compression curve agrees within error with the static compression results of Agee and Walker [1988a] to 6 GPa. We detennine that olivine (FO₉₄) will be neutrally buoyant in komatiitic melt of the composition we studied near 8.2 GPa. Clinopyroxene would also be neutrally buoyant near this pressure. Liquidus garnet-majorite may be less dense than this komatiitic liquid in the 20-24 GPa interval, however pyropic-garnet and perovskite phases are denser than this komatiitic liquid in their respective liquidus pressure intervals to 36 GPa. Liquidus perovskite may be neutrally buoyant near 70 GPa.

At 40 GPa, the density of shock-compressed molten komatiite would be approximately equal to the calculated density of an equivalent mixture of dense solid oxide components. This observation supports the model of Rigden et al. [1989] for compressibilities of liquid oxide components. Using their theoretical EOS for liquid forsterite and fayalite, we calculate the densities of a spectrum of melts from basaltic through peridotitic that are related to the experimentally studied komatiitic liquid by addition or subtraction of olivine. At low pressure, olivine fractionation lowers the density of basic magmas, but above 14 GPa this trend is reversed. All of these basic to ultrabasic liquids are predicted to have similar densities at 14 GPa, and this density is approximately equal to the bulk (PREM) mantle. This suggests that melts derived from a peridotitic mantle may be inhibited from ascending from depths greater than 400 km.

The EOS of ultrabasic magmas was used to model adiabatic melting in a peridotitic mantle. If komatiites are formed by >15% partial melting of a peridotitic mantle, then komatiites generated by adiabatic melting come from source regions in the lower transition zone (≈500-670 km) or the lower mantle (>670 km). The great depth of incipient melting implied by this model, and the melt density constraint mentioned above, suggest that komatiitic volcanism may be gravitationally hindered. Although komatiitic magmas are thought to separate from their coexisting crystals at a temperature =200°C greater than that for modern MORBs, their ultimate sources are predicted to be diapirs that, if adiabatically decompressed from initially solid mantle, were more than 700°C hotter than the sources of MORBs and derived from great depth.

We considered the evolution of an initially molten mantle, i.e., a magma ocean. Our model considers the thermal structure of the magma ocean, density constraints on crystal segregation, and approximate phase relationships for a nominally chondritic mantle. Crystallization will begin at the core-mantle boundary. Perovskite buoyancy at > 70 GPa may lead to a compositionally stratified lower mantle with iron-enriched mangesiowiistite content increasing with depth. The upper mantle may be depleted in perovskite components. Olivine neutral buoyancy may lead to the formation of a dunite septum in the upper mantle, partitioning the ocean into upper and lower reservoirs, but this septum must be permeable.

(2) Viscosity Measurement with Shock Waves

We have examined in detail the analytical method for measuring shear viscosity from the decay of perturbations on a corrugated shock front The relevance of initial conditions, finite shock amplitude, bulk viscosity, and the sensitivity of the measurements to the shock boundary conditions are discussed. The validity of the viscous perturbation approach is examined by numerically solving the second-order Navier-Stokes equations. These numerical experiments indicate that shock instabilities may occur even when the Kontorovich-D'yakov stability criteria are satisfied. The experimental results for water at 15 GPa are discussed, and it is suggested that the large effective viscosity determined by this method may reflect the existence of ice VII on the Rayleigh path of the Hugoniot This interpretation reconciles the experimental results with estimates and measurements obtained by other means, and is consistent with the relationship of the Hugoniot with the phase diagram for water. Sound waves are generated at 4.8 MHz at in the water experiments at 15 GPa. The existence of anelastic absorption modes near this frequency would also lead to large effective viscosity estimates.

(3) Equation of State of Molybdenum at 1400°C

Shock compression data to 96 GPa for pure molybdenum, initially heated to 1400°C, are presented. Finite strain analysis of the data gives a bulk modulus at 1400°C, K'_S. of 244±2 GPa and its pressure derivative, K'_OS of 4. A fit of shock velocity to particle velocity gives the coefficients of U_S = C_O+S U_P to be C_O = 4.77±0.06 km/s and S = 1.43±0.05. From the zero pressure sound speed, C_O, a bulk modulus of 232±6 GPa is calculated that is consistent with extrapolation of ultrasonic elasticity measurements. The temperature derivative of the bulk modulus at zero pressure, θK_OSθT|_P, is approximately -0.012 GPa/K. A thermodynamic model is used to show that the thermodynamic Grüneisen parameter is proportional to the density and independent of temperature. The Mie-Grüneisen equation of state adequately describes the high temperature behavior of molybdenum under the present range of shock loading conditions.

Brillouin scattering of cluster-assembled carbon films

The low frequency vibrational spectrum of cluster beam deposited carbon films was studied by Brillouin light scattering. In thin films the values of both bulk modulus and shear modulus has been estimated from the shifts of surface phonon peaks. The values found indicate a mainly sp2 coordinated random network with low density. In thick films a bulk longitudinal phonon peak was detected in a spectral range compatible with the value of the index of refraction and of the elastic constants of thin films. High surface roughness, combined with a rather strong bulk central peak, prevented the observation of surface phonon features. © 1998 Elsevier Science Ltd. All rights reserved.

Elastic properties of crystalline-amorphous core-shell silicon nanowires

The pressure behavior of Raman frequencies and line widths of crystalline core-amorphous shell silicon nanowires (SiNWs) with two different core-to-shell ratio thicknesses was studied at pressures up to 8 GPa. The obtained isothermal compressibility (bulk modulus) of SiNWs with a core-to-shell ratio of about 1.8 is ∼20% higher (lower) than reported values for bulk Si. For SiNWs with smaller core-to-shell ratios, a plastic deformation of the shell was observed together with a strain relaxation. A significant increase in the full width at half-maximum of the Raman LTO-peak due to phonon decay was used to determine the critical pressure at which LTO-phonons decay into LO + TA phonons. Our results reveal that this critical pressure in strained core-shell SiNWs (∼4 GPa) is different from the reported value for bulk Si (∼7 GPa), whereas no change is observed for relaxed core-shell SiNWs. © 2013 American Chemical Society.

First-principles prediction of the magnetism of 3d transition-metal-doped Rocksalt MgO

Using first-principles band structure methods, we have systematically studied the electronic structures, magnetic stabilities, and half-metal properties of 3d transition-metal (TM) doped Rocksalt MgO compounds TMMg3O4 (TM = V, Cr, Mn, Fe, Co, and Ni). The calculations reveal that only CrMg3O4 has a ferromagnetic stability among the six compounds, which is explained by double-exchange mechanism. The magnetic stability is affected by the doping concentration of TM if the top valance band is composed of partially occupied t(2g) states. In addition, CrMg3O4 is a half-metallic ferromagnet. The origins of half-metallic and ferromagnetic properties are explored. The Curie temperature (T-c) of CrMg3O4 is 182 K. And it is hard for CrMg3O4 to deform due to the large bulk modulus and shear modulus, so it is a promising spintronic material. Our calculations provide the first available information on the magnetic properties of 3d TM-doped MgO.

Mechanical and chemical bonding properties of ground state BeH2

The crystal structure, mechanical properties and electronic structure of ground state BeH2 are calculated employing the first-principles methods based on the density functional theory. Our calculated structural parameters at equilibrium volume are well consistent with experimental results. Elastic constants, which well obey the mechanical stability criteria, are firstly theoretically acquired. The bulk modulus B, Shear modulus G, Young's modulus E, and Poisson's ratio upsilon are deduced from the elastic constants. The bonding nature in BeH2 is fully interpreted by combining characteristics in band structure, density of states, and charge distribution. The ionicity in the Be-H bond is mainly featured by charge transfer from Be 2s to H 1s atomic orbitals while its covalency is dominated by the hybridization of H 1s and Be 2p states. The Bader analysis of BeH2 and MgH2 are performed to describe the ionic/covalent character quantitatively and we find that about 1.61 (1.6) electrons transfer from each Be (Mg) atom to H atoms.

First-principles study of ground-state properties and high pressure behavior of ThO2

The mechanical properties, electronic structure and phonon dispersion of ground state ThO2 as well as the structure behavior up to 240 GPa are studied using first-principles density-functional theory. Our calculated elastic constants indicate that both the ground-state fluorite structure and high pressure cotunnite structure of ThO2 are mechanically stable. The bulk modulus, shear modulus, and Young's modulus of cotunnite ThO2 are all smaller by approximately 25% compared with those of fluorite ThO2. The Poisson's ratios of both structures are approximately equal to 0.3 and the hardness of fluorite ThO2 is 22.4 GPa. The electronic structure and bonding nature of fluorite ThO2 are fully analyzed, and show that the Th-O bond displays a mixed ionic/covalent character. The phase transition from the fluorite to cotunnite structure is calculated to occur at the pressure of 26.5 GPa, consistent with recent experimental measurement by ldiri et al. [1]. For the cotunnite phase it is further predicted that an isostructural transition takes place in the pressure region of 80-130 GPa.

Elasticity, band structure, and piezoelectricity of BexZn1-xO alloys

Lattice constants, elasticity, band structure and piezoelectricity of hexagonal wideband gap BexZn1-xO ternary alloys are calculatedusing firstprinciples methods. The alloys' lattice constants obey Vegard's law well. As Be concentration increases, the bulk modulus and Young's modulus of the alloys increase, whereas the piezoelectricity decreases. We predict that BexZn1-xO/GaN/substrate (x = 0.022) multilayer structure can be suitable for high-frequency surface acoustic wave device applications. Our calculated results are in good agreement with experimental data and other theoretical calculations. (c) 2008 Elsevier B.V. All rights reserved.

Elasticity, band-gap bowing, and polarization of AlxGa1-xN alloys

Elastic constants, the bulk modulus, Young's modulus, band-gap bowing coefficients, spontaneous and piezoelectric polarizations, and piezoelectric coefficients of hexagonal AlxGa1-xN ternary alloys are calculated using first-principles methods. The fully relaxed structures and the structures subjected to homogeneous biaxial and uniaxial tension are investigated. We show that the biaxial tension in the plane perpendicular to the c axis and the uniaxial tension along the c axis all reduce the bulk modulus, whereas they reduce and enhance Young's modulus, respectively. We find that the biaxial and uniaxial tension can enhance the bowing coefficients. We also find that the biaxial tension can enhance the total polarization, while the uniaxial tension will suppress the total polarization. (C) 2008 American Institute of Physics.

Structural phase transformations of ZnS nanocrystalline under high pressure

In-situ energy dispersive x-ray diffraction on ZnS nanocrystalline was carried out under high pressure by using a diamond anvil cell. Phase transition of wurtzite of 10 nm ZnS to rocksalt occurred at 16.0 GPa, which was higher than that of the bulk materials. The structures of ZnS nanocrystalline at different pressures were built by using materials studio and the bulk modulus, and the pressure derivative of ZnS nanocrystalline were derived by fitting the equation of Birch-Murnaghan. The resulting modulus was higher than that of the corresponding bulk material, which indicates that the nanomaterial has higher hardness than its bulk materials.

An investigation on the pressure-induced phase transition of nanocrystalline ZnS

An in situ energy dispersive x-ray diffraction study on nanocrystalline ZnS was carried out under high pressure up to 30.8 GPa by using a diamond anvil cell. The phase transition from the wurtzite to the zinc-blende structure occurred at 11.5 GPa, and another obvious transition to a new phase with rock-salt structure also appeared at 16.0 GPa-which was higher than the value for the bulk material. The bulk modulus and the pressure derivative of nanocrystalline ZnS were derived by fitting the Birch-Murnaghan equation. The resulting modulus was higher than that of the corresponding bulk material, indicating that the nanomaterial has higher hardness than the bulk material.

Pressure-induced phase transition of nanocrystalline zinc sulfide

In situ energy dispersive X-ray diffraction measurements on nanocrystalline zinc sulfide have been performed by using diamond anvil cell with synchrotron radiation. There is a phase transition which the ultimate structure is rocksalt when the pressure is up to 16.0GPa. Comparing the structure of body materials, the pressure of the phase transition of nano zinc sulfide is high. We fit the: Birch-Murnaghan equation of state and obtained its ambient pressure bulk modulus and its pressure derivative. The bulk modulus of nanocrystalline zinc sulfide is higher than that of body materials, it indicate that the rigidity of nanocrystalline zinc sulfide is high.

Electronic, structural, and optical properties of crystalline yttria

The electronic structure of crystalline Y2O3 is investigated by first-principles calculations within the local-density approximation (LDA) of the density-functional theory. Results are presented for the band structure, the total density of states (DOS), the atom-and orbital-resolved partial DOS. effective charges, bond order, and charge-density distributions. Partial covalent character in the Y-O bonding is shown, and the nonequivalency of the two Y sites is demonstrated. The calculated electronic structure is compared with a variety of available experimental data. The total energy of the crystal is calculated as a function of crystal volume. A bulk modulus B of 183 Gpa and a pressure coefficient B' of 4.01 are obtained, which are in good agreement with compression data. An LDA band gap of 4.54 eV at Gamma is obtained which increases with pressure at a rate of dE(g)/dP = 0.012 eV/Gpa at the equilibrium volume. Also investigated are the optical properties of Y2O3 up to a photon energy of 20 eV. The calculated complex dielectric function and electron-energy-loss function are in good agreement with experimental data. A static dielectric constant of epsilon(O)= 3.20 is obtained. It is also found that the bottom of the conduction band consists of a single band, and direct optical transition at Gamma between the top of the valence band and the bottom of the conduction band may be symmetry forbidden.