Depletion of nonlinearity in magnetohydrodynamic turbulence: Insights from analysis and simulations

It is shown how suitably scaled, order-m moments, D-m(+/-), of the Elsasser vorticity fields in three-dimensional magnetohydrodynamics (MHD) can be used to identify three possible regimes for solutions of the MHD equations with magnetic Prandtl number P-M = 1. These vorticity fields are defined by omega(+/-) = curl z(+/-) = omega +/- j, where z(+/-) are Elsasser variables, and where omega and j are, respectively, the fluid vorticity and current density. This study follows recent developments in the study of three-dimensional Navier-Stokes fluid turbulence Gibbon et al., Nonlinearity 27, 2605 (2014)]. Our mathematical results are then compared with those from a variety of direct numerical simulations, which demonstrate that all solutions that have been investigated remain in only one of these regimes which has depleted nonlinearity. The exponents q(+/-) that characterize the inertial range power-law dependencies of the z(+/-) energy spectra, epsilon(+/-)(k), are then examined, and bounds are obtained. Comments are also made on (a) the generalization of our results to the case P-M not equal 1 and (b) the relation between D-m(+/-) and the order-m moments of gradients of magnetohydrodynamic fields, which are used to characterize intermittency in turbulent flows.

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.

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.

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.

Creep response of sandwich beams with a metallic foam core

The creep response of metallic foam sandwich beams in 3-point bend is investigated numerically for the case of a metallic foam core and two steel faces. The face sheets are treated as elastic, while the foam core is modeled by a viscoplastic extension of the Deshpande-Fleck yield surface. This power-law creeping constitutive law has been implemented within the commercial finite element code ABAQUS. It is found that the beams creep by a variety of competing mechanisms, depending upon the choice of material properties and the geometric parameters. A failure map is constructed and effect of rate dependence on the load-deflection curves is quantified, and compared against the available experimental data.

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.

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.

Experimental laws of cratering for hypervelocity impacts of spherical projectiles into thick target

It is shown in this paper that the laws of cratering in a thick target under hypervelocity impact by a spherical projectile can be approximately expressed by the so-called iso-deviation law and a 2/3 power law. Moreover, hypervelocity impact should be characterized by the isotropic expansion of a crater. In the special case, when the projectile and target are of the same material, the laws mentioned above reduce to the result of a semi-spherical crater and the energy criterion. Generally speaking, a semi-spherical crater and the energy criterion are both approximations, which only take projectile density and target strength into account, and can be used for a rough estimation on the order of magnitude. The inconsistency in various fitted power laws in the literature was also clarified and explained in the paper.

Crack growth by grain boundary cavitation in the transient and extensive creep regimes

Crack growth due to cavity growth and coalescence along grain boundaries is analyzed under transient and extensive creep conditions in a compact tension specimen. Account is taken of the finite geometry changes accompanying crack tip blunting. The material is characterized as an elastic-power law creeping solid with an additional contribution to the creep rate arising from a given density of cavitating grain boundary facets. All voids are assumed present from the outset and distributed on a given density of cavitating grain boundary facets. The evolution of the stress fields with crack growth under three load histories is described in some detail for a relatively ductile material. The full-field plane strain finite element calculations show the competing effects of stress relaxation due to constrained creep, diffusion and crack tip blunting. and of stress increase due to the instantaneous elastic response to crack growth. At very high crack growth rates the Hui-Riedel fields dominate the crack tip region. However. the high growth rates are not sustained for any length of time in the compact tension geometry analyzed. The region of dominance of the Hui-Riedel field shrinks rapidly so that the near-tip fields are controlled by the HRR-type field shortly after the onset of crack growth. Crack growth rates under various conditions of loading and spanning the range of times from small scale creep to extensive creep are obtained. We show that there is a strong similarity between crack growth history and the behaviour of the C(t) and C(t) parameters. so that crack growth rates correlate rather well with C(t) and C(t). A relatively brittle material is also considered that has a very different near-tip stress field and crack growth history.

Elastic-plastic finite element analyses of short cracks

Stress and strain distributions and crack opening displacement characteristics of short cracks have been studied in single edge notch bend and centre cracked panel specimens using elastic–plastic finite element analyses incorporating both a non strain hardening and a power law hardening behaviour. J contour integral solutions to describe stress strain conditions at crack tips for short cracks differ from those for long cracks. The analyses show that (i) short cracks can propagate at stress levels lower than those required for long cracks and (ii) a two-parameter description of crack tip fields is necessary for crack propagation.