Indentation of elastically anisotropic half-spaces by cones and parabolae of revolution

Indentation of ceramic materials with smooth indenters such as parabolae of revolution and spheres can be conducted in the elastic regime to relatively high loads. Ceramic single crystals thus provide excellent calibration media for load-and depth-sensing indentation testing; however, they are generally anisotropic and a complete elastic analysis is cumbersome. This study presents a simplified procedure for the determination of the stiffness of contact for the indentation of an anisotropic half-space by a rigid frictionless parabola of revolution which, to first order, approximates spherical indentation. Using a similar approach, a new procedure is developed for analysing conical indentation of anisotropic elastic media. For both indenter shapes, the contact is found to be elliptical, and equations are determined for the size, shape and orientation of the ellipse and the indentation modulus.

Exponents of Subvarieties of Upper Triangular Matrices over Arbitrary Fields are Integral

Partially supported by grant RFFI 98-01-01020.

Anisotropic characterization of crack growth in the tertiary flow of asphalt mixtures in compression

Asphalt mixtures exhibit primary, secondary, and tertiary stages in sequence during a rutting deterioration. Many field asphalt pavements are still in service even when the asphalt layer is in the tertiary stage, and rehabilitation is not performed until a significant amount of rutting accompanied by numerous macrocracks is observed. The objective of this study was to provide a mechanistic method to model the anisotropic cracking of the asphalt mixtures in compression during the tertiary stage of rutting. Laboratory tests including nondestructive and destructive tests were performed to obtain the viscoelastic and viscofracture properties of the asphalt mixtures. Each of the measured axial and radial total strains in the destructive tests were decomposed into elastic, plastic, viscoelastic, viscoplastic, and viscofracture strains using the pseudostrain method in an extended elastic-viscoelastic correspondence principle. The viscofracture strains are caused by the crack growth, which is primarily signaled by the increase of phase angle in the tertiary flow. The viscofracture properties are characterized using the anisotropic damage densities (i.e., the ratio of the lost area caused by cracks to the original total area in orthogonal directions). Using the decomposed axial and radial viscofracture strains, the axial and radial damage densities were determined by using a dissipated pseudostrain energy balance principle and a geometric analysis of the cracks, respectively. Anisotropic pseudo J-integral Paris' laws in terms of damage densities were used to characterize the evolution of the cracks in compression. The material constants in the Paris' law are determined and found to be highly correlated. These tests, analysis, and modeling were performed on different asphalt mixtures with two binders, two air void contents, and three aging periods. Consistent results were obtained; for instance, a stiffer asphalt mixture is demonstrated to have a higher modulus, a lower phase angle, a greater flow number, and a larger n1 value (exponent of Paris' law). The calculation of the orientation of cracks demonstrates that the asphalt mixture with 4% air voids has a brittle fracture and a splitting crack mode, whereas the asphalt mixture with 7% air voids tends to have a ductile fracture and a diagonal sliding crack mode. Cracks of the asphalt mixtures in compression are inclined to propagate along the direction of the external compressive load. © 2014 American Society of Civil Engineers.

Anisotropic viscoelastic properties of undamaged asphalt mixtures

A test protocol and a data analysis method are developed in this paper on the basis of linear viscoelastic theory to characterize the anisotropic viscoelastic properties of undamaged asphalt mixtures. The test protocol includes three nondestructive tests: (1) uniaxial compressive creep test, (2) indirect tensile creep test, and (3) the uniaxial tensile creep test. All three tests are conducted on asphalt mixture specimens at three temperatures (10, 20, and 30°C) to determine the tensile and compressive properties at each temperature and then to construct the master curve of each property. The determined properties include magnitude and phase angle of the compressive complex modulus in the vertical direction, magnitude and phase angle of the tensile complex modulus, and the magnitude and phase angle of the compressive complex modulus in the horizontal plane. The test results indicate that all tested asphalt mixtures have significantly different tensile properties from compressive properties. The peak value of the master curve of the tensile complex modulus phase angle is within a range from 65 to 85°, whereas the peak value of the compressive moduli phase angle in both directions ranges from 35 to 55°. In addition, the undamaged asphalt mixtures exhibit distinctively anisotropic properties in compression. The magnitude of the compressive modulus in the vertical direction is approximately 1.2 to ̃2 times of the magnitude of the compressive modulus in the horizontal plane. Dynamic modulus tests are performed to verify the results of the proposed test protocol. The test results from the proposed test protocol match well with those from the dynamic tests. © 2012 American Society of Civil Engineers.

Z2-Graded Polynomial Identities for Superalgebras of Block-Triangular Matrices

000 Mathematics Subject Classification: Primary 16R50, Secondary 16W55.

Triangular Models and Asymptotics of Continuous Curves with Bounded and Unbounded Semigroup Generators

2000 Mathematics Subject Classification: Primary 47A48, Secondary 60G12.

Hypoellipticity of Anisotropic Partial Differential Equations

2002 Mathematics Subject Classification: 35S05

Nonlinear Estimates in Anisotropic Gevrey Spaces

2002 Mathematics Subject Classification: 35G20, 47H30

Little G. T. for lp-lattice summing operators

2000 Mathematics Subject Classification: 46B28, 47D15.

Anisotropic opinion dynamics

We consider the process of opinion formation in a society of interacting agents, where there is a set B of socially accepted rules. In this scenario, we observed that agents, represented by simple feed-forward, adaptive neural networks, may have a conservative attitude (mostly in agreement with B) or liberal attitude (mostly in agreement with neighboring agents) depending on how much their opinions are influenced by their peers. The topology of the network representing the interaction of the society's members is determined by a graph, where the agents' properties are defined over the vertexes and the interagent interactions are defined over the bonds. The adaptability of the agents allows us to model the formation of opinions as an online learning process, where agents learn continuously as new information becomes available to the whole society (online learning). Through the application of statistical mechanics techniques we deduced a set of differential equations describing the dynamics of the system. We observed that by slowly varying the average peer influence in such a way that the agents attitude changes from conservative to liberal and back, the average social opinion develops a hysteresis cycle. Such hysteretic behavior disappears when the variance of the social influence distribution is large enough. In all the cases studied, the change from conservative to liberal behavior is characterized by the emergence of conservative clusters, i.e., a closed knitted set of society members that follow a leader who agrees with the social status quo when the rule B is challenged.

Exotic meson decay widths using lattice quantum chromodynamics

Most experiments in particle physics are scattering experiments, the analysis of which leads to masses, scattering phases, decay widths and other properties of one or multi-particle systems. Until the advent of Lattice Quantum Chromodynamics (LQCD) it was difficult to compare experimental results on low energy hadron-hadron scattering processes to the predictions of QCD, the current theory of strong interactions. The reason being, at low energies the QCD coupling constant becomes large and the perturbation expansion for scattering; amplitudes does not converge. To overcome this, one puts the theory onto a lattice, imposes a momentum cutoff, and computes the integral numerically. For particle masses, predictions of LQCD agree with experiment, but the area of decay widths is largely unexplored. ^ LQCD provides ab initio access to unusual hadrons like exotic mesons that are predicted to contain real gluonic structure. To study decays of these type resonances the energy spectra of a two-particle decay state in a finite volume of dimension L can be related to the associated scattering phase shift δ(k) at momentum k through exact formulae derived by Lüscher. Because the spectra can be computed using numerical Monte Carlo techniques, the scattering phases can thus be determined using Lüscher's formulae, and the corresponding decay widths can be found by fitting Breit-Wigner functions. ^ Results of such a decay width calculation for an exotic hybrid( h) meson (JPC = 1-+) are presented for the decay channel h → πa 1. This calculation employed Lüscher's formulae and an approximation of LQCD called the quenched approximation. Energy spectra for the h and πa1 systems were extracted using eigenvalues of a correlation matrix, and the corresponding scattering phase shifts were determined for a discrete set of πa1 momenta. Although the number of phase shift data points was sparse, fits to a Breit-Wigner model were made, resulting in a decay width of about 60 MeV. ^

Lattice Boltzmann method for flow and heat transfer in microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

Effect of water on olivine single crystals plasticity, deformed under high pressure and high temperature

The Earth's upper mantle, mainly composed of olivine, is seismically anisotropic. Seismic anisotropy attenuation has been observed at 220km depth. Karato et al. (1992) attributed this attenuation to a transition between two deformation mechanisms, from dislocation creep above 220km to diffusion creep below 220km, induced by a change in water content. Couvy (2005) and Mainprice et al. (2005) predicted a change in Lattice Preferred Orientation induced by pressure, which comes from a change of slip system, from [100] slip to [001] slip, and is responsible for the seismic anisotropy attenuation. Raterron et al. (2007) ran single crystal deformation experiments under anhydrous conditions and observed that the slip system transition occurs around 8GPa, which corresponds to a depth of 260Km. Experiments were done to quantify the effects of water on olivine single crystals deformed using D-DIA press and synchrotron beam. Deformations were carried out in uniaxial compression along [110]c, [011]c, and [101]c, crystallographic directions, at pressure ranging from 4 to 8GPa and temperature between 1373 and 1473K. Talc sleeves about the annulus of the single crystals were used as source of water in the assembly. Stress and specimen strain rates were calculated by in-situ X-ray diffraction and time resolved imaging, respectively. By direct comparison of single crystals strain rates, we observed that [110]c deforms faster than [011]c below 5GPa. However above 6GPa [011]c deforms faster than [110]c. This revealed that [100](010) is the dominant slip system below 5GPa, and above 6GPa [001](010) becomes dominant. According to our results, the slip system transition, which is induced by pressure, occurs at 6GPa. Water influences the pressure where the switch over occurs, by lowering the transition pressure. The pressure effect on the slip systems activity has been quantified and the hydrolytic weakening has also been estimated for both orientations. Data also shows that temperature affects the slip system activity. The regional variation of the depth for the seismic anisotropy attenuation, which would depend on local hydroxyl content and temperature variations and explains the seismic anisotropy attenuation occurring at about 220Km depth in the mantle, where the pressure is about 6GPa. Deformation of MgO single crystal oriented [100], [110] and [111] were also performed. The results predict a change in the slip system activity at 23GPa, again induced by pressure. This explains the seismic anisotropy observed in the lower mantle.

Lattice boltzmann modeling of fluid flow to determine the permeability of a karst specimen

A limestone sample was scanned using computed tomography (CT) and the hydraulic conductivity of the 3D reconstructed sample was determined using Lattice- Boltzmann methods (LBM) at varying scales. Due to the shape and size of the original sample, it was challenging to obtain a consistent rectilinear test sample. Through visual inspection however, 91 mm and 76 mm samples were digitally cut from the original. The samples had porosities of 58% and 64% and produced hydraulic conductivity values of K= 13.5 m/s and K=34.5 m/s, respectively. Both of these samples were re-sampled to 1/8 and 1/64 of their original size to produce new virtual samples at lower resolutions of 0.542 mm/lu and 1.084 mm/lu, while still representing the same physical dimensions. The hydraulic conductivity tended to increase slightly as the resolution became coarser. In order to determine an REV, the 91 mm sample was also sub-sampled into blocks that were 1/8 and 1/64 the size of the original. The results were consistent with analytical expectations such as those produced by the Kozeny-Carman equation. A definitive REV size was not reached, however, indicating the need for a larger sample. The methods described here demonstrate the ability of LBM to test rock structures and sizes not normally attainable.

Lattice Boltzmann Method for Flow and Heat Transfer in Microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).