Dielectric response of graded composites having general power-law-graded cylindrical inclusions

The dielectric response of graded composites having general power-law-graded cylindrical inclusions under a uniform applied electric field is investigated. The dielectric profile of the cylindrical inclusions is modeled by the equation epsilon(i)(r)=c(b+r)(k) (where r is the radius of the cylindrical inclusions and c, b and k are parameters). Analytical solutions for the local electrical potentials are derived in terms of hypergeometric functions and the effective dielectric response of the graded composites is predicted in the dilute limit. Moreover, for a simple power-law dielectric profile epsilon(i)(r) = cr(k) and a linear dielectric profile epsilon(i)(r) = c(b + r), analytical expressions of the electrical potentials and the effective dielectric response are derived exactly from our results by taking the limits b -> 0 and k -> 1, respectively. For a higher concentration of inclusions, the effective dielectric response is estimated by an effective-medium approximation. In addition, we have discussed the effective response of graded cylindrical composites with a more complex dielectric profile of inclusion, epsilon(i)(r)=c(b+r)(k)e(beta r). (c) 2005 American Institute of Physics.

Effective dielectric responses of graded composites of the particles with a general power-law gradation profile

The effective dielectric response of graded spherical composites having general power-law gradient inclusions is investigated under a uniform applied electric field, where the dielectric gradation profile of the spherical inclusions is modeled by the equation epsilon(i) (r) = c(b+r)(k). Analytical solutions of the local electrical potentials are derived in terms of hyper-geometric function and the effective dielectric response of the graded composites is predicted in the dilute limit. From our result, the local potentials of graded spherical composites having both simple power-law dielectric profile epsilon(i)(r) = cr(k) and linear dielectric profile epsilon(i) (r) = c(b+r) are derived exactly by taking the limits b --> 0 and k --> 1, respectively. In the dilute limit, our exact result is used to test the validity of differential effective dipole approximation (DEDA) for estimating the effective response of graded spherical composites, and it is shown that the DEDA is in excellent agreement with exact result. (C) 2005 Elsevier B.V. All rights reserved.

Load-Unload Response Ratio and Accelerating Moment/Energy Release Critical Region Scaling and Earthquake Prediction

The main idea of the Load-Unload Response Ratio (LURR) is that when a system is stable, its response to loading corresponds to its response to unloading, whereas when the system is approaching an unstable state, the response to loading and unloading becomes quite different. High LURR values and observations of Accelerating Moment/Energy Release (AMR/AER) prior to large earthquakes have led different research groups to suggest intermediate-term earthquake prediction is possible and imply that the LURR and AMR/AER observations may have a similar physical origin. To study this possibility, we conducted a retrospective examination of several Australian and Chinese earthquakes with magnitudes ranging from 5.0 to 7.9, including Australia's deadly Newcastle earthquake and the devastating Tangshan earthquake. Both LURR values and best-fit power-law time-to-failure functions were computed using data within a range of distances from the epicenter. Like the best-fit power-law fits in AMR/AER, the LURR value was optimal using data within a certain epicentral distance implying a critical region for LURR. Furthermore, LURR critical region size scales with mainshock magnitude and is similar to the AMR/AER critical region size. These results suggest a common physical origin for both the AMR/AER and LURR observations. Further research may provide clues that yield an understanding of this mechanism and help lead to a solid foundation for intermediate-term earthquake prediction.

Fractals and Scaling in Fracture Induced by Microcrack Coalescence

A numerical model is proposed to simulate fracture induced by the coalescence of numerous microcracks, in which the condition for coalescence between two randomly nucleated microcracks is determined in terms of a load-sharing principle. The results of the simulation show that, as the number density of nucleated microcracks increases, stochastic coalescence first occurs followed by a small fluctuation, and finally a newly nucleated microcrack triggers a cascade coalescence of microcracks resulting in catastrophic failure. The fracture profiles exhibit self-affine fractal characteristics with a universal roughness exponent, but the critical damage threshold is sensitive to details of the model. The spatiotemporal distribution of nucleated microcracks in the vicinity of critical failure follows a power-law behaviour, which implies that the microcrack system may evolve to a critical state.

Similarity consideration of structural bifurcation buckling

Many structural bifurcation buckling problems exhibit a scaling or power law property. Dimensional analysis is used to analyze the general scaling property. The concept of a new dimensionless number, the response number-Rn, suggested by the present author for the dynamic plastic response and failure of beams, plates and so on, subjected to large dynamic loading, is generalized in this paper to study the elastic, plastic, dynamic elastic as well as dynamic plastic buckling problems of columns, plates as well as shells. Structural bifurcation buckling can be considered when Rn(n) reaches a critical value.

Calculations of longitudinal and transverse velocity structure functions using a vortex model of isotropic turbulence

The longitudinal structure function (LSF) and the transverse structure function (TSF) in isotropic turbulence are calculated using a vortex model. The vortex model is composed of the Rankine and Burgers vortices which have the exponential distributions in the vortex Reynolds number and vortex radii. This model exhibits a power law in the inertial range and satisfies the minimal condition of isotropy that the second-order exponent of the LSF in the inertial range is equal to that of the TSF. Also observed are differences between longitudinal and transverse structure functions caused by intermittency. These differences are related to their scaling differences which have been previously observed in experiments and numerical simulations.

Finite Element Solutions for Plane Strain Mode I Crack with Strain Grading Effect

In this paper, a new phenomenological theory with strain gradient effects is proposed to account for the size dependence of plastic deformation at micro- and submicro-length scales. The theory fits within the framework of general couple stress theory and three rotational degrees of freedom omega(i) are introduced in addition to the conventional three translational degrees of freedom mu(i). omega(i) is called micro-rotation and is the sum of material rotation plus the particles' relative rotation. While the new theory is used to analyze the crack tip field or the indentation problems, the stretch gradient is considered through a new hardening law. The key features of the theory are that the rotation gradient influences the material character through the interaction between the Cauchy stresses and the couple stresses; the term of stretch gradient is represented as an internal variable to increase the tangent modulus. In fact the present new strain gradient theory is the combination of the strain gradient theory proposed by Chen and Wang (Int. J. Plast., in press) and the hardening law given by Chen and Wang (Acta Mater. 48 (2000a) 3997). In this paper we focus on the finite element method to investigate material fracture for an elastic-power law hardening solid. With remotely imposed classical K fields, the full field solutions are obtained numerically. It is found that the size of the strain gradient dominance zone is characterized by the intrinsic material length l(1). Outside the strain gradient dominance zone, the computed stress field tends to be a classical plasticity field and then K field. The singularity of stresses ahead of the crack tip is higher than that of the classical field and tends to the square root singularity, which has important consequences for crack growth in materials by decohesion at the atomic scale. (C) 2002 Elsevier Science Ltd. All rights reserved.

What is indentation hardness?

Using dimensional analysis and finite element calculations, we derive simple scaling relationships for loading and unloading curve, contact depth, and hardness. The relationship between hardness and the basic mechanical properties of solids, such as Young's modulus, initial yield strength, and work-hardening exponent, is then obtained. The conditions for 'piling-up' and 'sinking-in' of surface profiles during indentation are determined. A method for estimating contact depth from initial unloading slope is examined. The work done during indentation is also studied. A relationship between the ratio of hardness to elastic modulus and the ratio of irreversible work to total work is discovered. This relationship offers a new method for obtaining hardness and elastic modulus. Finally, a scaling theory for indentation in power-law creep solids using self-similar indenters is developed. A connection between creep and 'indentation size effect' is established.

Scaling, Dimensional Analysis, and Indentation Measurements

We provide an overview of the basic concepts of scaling and dimensional analysis, followed by a review of some of the recent work on applying these concepts to modeling instrumented indentation measurements. Specifically, we examine conical and pyramidal indentation in elastic-plastic solids with power-law work-hardening, in power-law creep solids, and in linear viscoelastic materials. We show that the scaling approach to indentation modeling provides new insights into several basic questions in instrumented indentation, including, what information is contained in the indentation load-displacement curves? How does hardness depend on the mechanical properties and indenter geometry? What are the factors determining piling-up and sinking-in of surface profiles around indents? Can stress-strain relationships be obtained from indentation load-displacement curves? How to measure time dependent mechanical properties from indentation? How to detect or confirm indentation size effects? The scaling approach also helps organize knowledge and provides a framework for bridging micro- and macroscales. We hope that this review will accomplish two purposes: (1) introducing the basic concepts of scaling and dimensional analysis to materials scientists and engineers, and (2) providing a better understanding of instrumented indentation measurements.

Mode I crack tip field with strain gradient effects

A plane strain mode I crack tip field with strain gradient effects is investigated. A new strain gradient theory is used. An elastic-power law hardening strain gradient material is considered and two hardening laws, i.e. a separation law and an integration Law are used respectively. As for the material with the separation law hardening, the angular distributions of stresses are consistent with the HRR field, which differs from the stress results([19]); the angular distributions of couple stresses are the same as the couple stress results([19]). For the material with the integration law hardening, the stress field and the couple stress field can not exist simultaneously, which is the same as the conclusion([19]), but for the stress dominated field, the angular distributions of stresses are consistent with the HRR field; for the couple stress dominated field, the angular distributions of couple stresses are consistent with those in Ref. [19]. However, the increase in stresses is not observed in strain gradient plasticity because the present theory is based on the rotation gradient of the deformation only, while the crack tip field of mode I is dominated by the tension gradient, which will be shown in another paper.

Avalanche behavior and statistical properties in a microcrack coalescence process

Fracture owing to the coalescence of numerous microcracks can be described by a simple statistical model, where a coalescence event stochastically occurs as the number density of nucleated microcracks increases. Both numerical simulation and statistical analysis reveal that a microcrack coalescence process may display avalanche behavior and that the final failure is catastrophic. The cumulative distribution of coalescence events in the vicinity of critical fracture follows a power law and the fracture profile has self-affine fractal characteristic. Some macromechanical quantities may be traced back and extracted from the mesoscopic process based on the statistical analysis of coalescence events.

Sediment incipience in turbulence generated in a square tank by a vertically oscillating grid

The failure of hydraulic structures in many estuaries and coastal regions around the world has been attributed to sediment transport and local scour. The sediment incipience in homogenous turbulence generated by oscillating grid is studied in this paper. The turbulent flow is measured by particle tracer velocimetry (PTV) technique. The integral length scale and time scale of turbulence are obtained. The turbulent flow near the wall is measured by local optical magnification. The sediment incipience is described by static theory. The relationship of probability of sediment incipience and the turbulent kinetic energy were obtained experimentally and theoretically. The distribution of the turbulent kinetic energy near the wall is found to obey the power law and the turbulent energy is further identified as the dynamic mechanism of sediment incipience.

Multi-scale analysis of AFM tip and surface interactions

Thoroughly understanding AFM tip-surface interactions is crucial for many experimental studies and applications. It is important to realize that despite its simple appearance, the system of tip and sample surface involves multiscale interactions. In fact, the system is governed by a combination of molecular force (like the van der Waals force), its macroscopic representations (such as surface force) and gravitational force (a macroscopic force). Hence, in the system, various length scales are operative, from sub-nanoscale (at the molecular level) to the macroscopic scale. By integrating molecular forces into continuum equations, we performed a multiscale analysis and revealed the nonlocality effect between a tip and a rough solid surface and the mechanism governing liquid surface deformation and jumping. The results have several significant implications for practical applications. For instance, nonlocality may affect the measurement accuracy of surface morphology. At the critical state of liquid surface jump, the ratio of the gap between a tip and a liquid dome (delta) over the dome height (y(o)) is approximately (n-4) (for a large tip), which depends on the power law exponent n of the molecular interaction energy. These findings demonstrate that the multiscale analysis is not only useful but also necessary in the understanding of practical phenomena involving molecular forces. (c) 2007 Elsevier Ltd. All rights reserved.

Direct dynamical test for deterministic chaos and optimal embedding of a chaotic time-series

We propose here a local exponential divergence plot which is capable of providing an alternative means of characterizing a complex time series. The suggested plot defines a time-dependent exponent and a ''plus'' exponent. Based on their changes with the embedding dimension and delay time, a criterion for estimating simultaneously the minimal acceptable embedding dimension, the proper delay time, and the largest Lyapunov exponent has been obtained. When redefining the time-dependent exponent LAMBDA(k) curves on a series of shells, we have found that whether a linear envelope to the LAMBDA(k) curves exists can serve as a direct dynamical method of distinguishing chaos from noise.

Direct dynamical test for deterministic chaos

We present a direct and dynamical method to distinguish low-dimensional deterministic chaos from noise. We define a series of time-dependent curves which are closely related to the largest Lyapunov exponent. For a chaotic time series, there exists an envelope to the time-dependent curves, while for a white noise or a noise with the same power spectrum as that of a chaotic time series, the envelope cannot be defined. When a noise is added to a chaotic time series, the envelope is eventually destroyed with the increasing of the amplitude of the noise.