Compressibility effects on hydromagnetic surface

The dispersion equation for hydromagnetic surface waves along a plasma-plasma interface has been solved as a function of the compressibility factor c 1/v A1, where c 1 and v A1 are the acoustic and Alfvén wave speed in one of the medium, for general wave propagation direction. Both slow and fast magnetosonic surface waves can exist. The nature and existence of these waves depends on the values of c 1/v A1 and theta, the angle of wave propagation. For low-beta plasmas only fast mode exists. The slow mode does not propagate below a critical value of c 1. When c 1 rarr infin the phase velocity of the slow wave tend to the Alfvén surface wave velocity in the incompressible media and for large theta the phase velocity of the fast wave approaches this value. The phase velocity of the slow wave increases whereas for the fast wave it decreases with increase in the angle theta.

Shear-imposed falling film

The study of a film falling down an inclined plane is revisited in the presence of imposed shear stress. Earlier studies regarding this topic (Smith, J. Fluid Mech., vol. 217, 1990, pp. 469-485; Wei, Phys. Fluids, vol. 17, 2005a, 012103), developed on the basis of a low Reynolds number, are extended up to moderate values of the Reynolds number. The mechanism of the primary instability is provided under the framework of a two-wave structure, which is normally a combination of kinematic and dynamic waves. In general, the primary instability appears when the kinematic wave speed exceeds the speed of dynamic waves. An equality criterion between their speeds yields the neutral stability condition. Similarly, it is revealed that the nonlinear travelling wave solutions also depend on the kinematic and dynamic wave speeds, and an equality criterion between the speeds leads to an analytical expression for the speed of a family of travelling waves as a function of the Froude number. This new analytical result is compared with numerical prediction, and an excellent agreement is achieved. Direct numerical simulations of the low-dimensional model have been performed in order to analyse the spatiotemporal behaviour of nonlinear waves by applying a constant shear stress in the upstream and downstream directions. It is noticed that the presence of imposed shear stress in the upstream (downstream) direction makes the evolution of spatially growing waves weaker (stronger).

Wallner lines, crack velocity and mechanisms of crack nucleation and growth in a brittle bulk metallic glass

Mode I fracture experiments were conducted on brittle bulk metallic glass (BMG) samples and the fracture surface features were analyzed in detail to understand the underlying physical processes. Wollner lines, which result from the interaction between the propagating crack front and shear waves emanating from a secondary source, were observed on the fracture surface and geometric analysis of them indicates that the maximum crack velocity is similar to 800 m s(-1), which corresponds to similar to 0.32 times the shear wave speed. Fractography reveals that the sharp crack nucleation at the notch tip occurs at the mid-section of the specimens with the observation of flat and half-penny-shaped cracks. On this basis, we conclude that the crack initiation in brittle BMGs is stress-controlled and occurs through hydrostatic stress-assisted cavity nucleation ahead of the notch tip. High magnification scanning electron and atomic force microscopies of the dynamic crack growth regions reveal highly organized, nanoscale periodic patterns with a spacing of similar to 79 nm. Juxtaposition of the crack velocity with this spacing suggests that the crack takes similar to 10(-10) s for peak-to-peak propagation. This, and the estimated adiabatic temperature rise ahead of the propagating crack tip that suggests local softening, is utilized to critically discuss possible causes for the nanocorrugation formation. Taylor's fluid meniscus instability is unequivocally ruled out. Then, two other possible mechanisms, viz. (a) crack tip blunting and resharpening through nanovoid nucleation and growth ahead of the crack tip and eventual coalescence, and (b) dynamic oscillation of the crack in a thin slab of softened zone ahead of the crack-tip, are critically discussed. (C) 2014 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Roll Waves in Overland Flow

Roll waves are frequently observed in overland flow, especially in rill flow, which has an important effect on the development of soil erosion. Using one-dimensional St. Venant equations, this paper investigates the dynamics of periodic roll waves based on Dressler’s and Brock’s work. Under the assumption that the average flow depth equals the uniform flow depth, expressions of the roll-wave speed and roll-wave profile were obtained. Testing with the results observed by Brock (1970) for wave properties shows that these expressions can approximately describe the characteristics of periodic permanent roll waves. Numerical solutions of roll waves under specific conditions are found, which show that when a roll wave appears, the shear stress of flow increases, and the soil erosion accelerates.

Dynamic Fracture Instability Of Tough Bulk Metallic Glass

We report the observations of a clear fractographic evolution from vein pattern, dimple structure, and then to periodic corrugation structure, followed by microbranching pattern, along the crack propagation direction in the dynamic fracture of a tough Zr41.2Ti13.8Cu12.5Ni10Be22.5 (Vit.1) bulk metallic glass (BMGs) under high-velocity plate impact. A model based on fracture surface energy dissipation and void growth is proposed to characterize this fracture pattern transition. We find that once the dynamic crack propagation velocity reaches a critical fraction of Rayleigh wave speed, the crack instability occurs; hence, crack microbranching goes ahead. Furthermore, the correlation between the critical velocity of amorphous materials and their intrinsic strength such as Young's modulus is uncovered. The results may shed new insight into dynamic fracture instability for BMGs. (C) 2008 American Institute of Physics.

Electromechanical influence of crack velocity at bifurcation for poled ferroelectric materials

The piezoelastodynamic field equations are solved to determine the crack velocity at bifurcation for poled ferroelectric materials where the applied electrical field and mechanical stress can be varied. The underlying physical mechanism, however, may not correspond to that assumed in the analytical model. Bifurcation has been related to the occurrence of a pair of maximum circumferential stress oriented symmetrically about the moving crack path. The velocity at which this behavior prevails has been referred to as the limiting crack speed. Unlike the classical approach, bifurcation will be identified with finite distances ahead of a moving crack. Nucleation of microcracks can thus be modelled in a single formulation. This can be accomplished by using the energy density function where fracture initiation is identified with dominance of dilatation in relation to distortion. Poled ferroelectric materials are selected for this study because the microstructure effects for this class of materials can be readily reflected by the elastic, piezoelectic and dielectric permittivity constants at the macroscopic scale. Existing test data could also shed light on the trend of the analytical predictions. Numerical results are thus computed for PZT-4 and compared with those for PZT-6B in an effort to show whether the branching behavior would be affected by the difference in the material microstructures. A range of crack bifurcation speed upsilon(b) is found for different r/a and E/sigma ratios. Here, r and a stand for the radial distance and half crack length, respectively, while E and a for the electric field and mechanical stress. For PZT-6B with upsilon(b) in the range 100-1700 m/s, the bifurcation angles varied from +/-6degrees to +/-39degrees. This corresponds to E/sigma of -0.072 to 0.024 V m/N. At the same distance r/a = 0.1, PZT-4 gives upsilon(b) values of 1100-2100 m/s; bifurcation angles of +/-15degrees to +/-49degrees; and E/sigma of -0.056 to 0.059 V m/N. In general, the bifurcation angles +/-theta(0) are found to decrease with decreasing crack velocity as the distance r/a is increased. Relatively speaking, the speed upsilon(b) and angles +/-theta(0) for PZT-4 are much greater than those for PZT-6B. This may be attributed to the high electromechanical coupling effect of PZT-4. Using upsilon(b)(0) as a base reference, an equality relation upsilon(b)(-) < upsilon(b)(0) < upsilon(b)(+) can be established. The superscripts -, 0 and + refer, respectively, to negative, zero and positive electric field. This is reminiscent of the enhancement and retardation of crack growth behavior due to change in poling direction. Bifurcation characteristics are found to be somewhat erratic when r/a approaches the range 10(-2)-10(-1) where the kinetic energy densities would fluctuate and then rise as the distance from the moving crack is increased. This is an artifact introduced by the far away condition of non-vanishing particle velocity. A finite kinetic energy density prevails at infinity unless it is made to vanish in the boundary value problem. Future works are recommended to further clarify the physical mechanism(s) associated with bifurcation by means of analysis and experiment. Damage at the microscopic level needs to be addressed since it has been known to affect the macrocrack speeds and bifurcation characteristics. (C) 2002 Published by Elsevier Science Ltd.

Molecular dynamics simulation of crack-tip processes in copper

The crack tip processes in copper under mode II loading have been simulated by a molecular dynamics method. The nucleation, emission, dislocation free zone (DFZ) and pile-up of the dislocations are analyzed by using a suitable atom lattice configuration and Finnis & Sinclair potential. The simulated results show that the dislocation emitted always exhibits a dissociated fashion. The stress intensity factor for dislocation nucleation, DFZ and dissociated width of partial dislocations are strongly dependent on the loading rate. The stress distributions are in agreement with the elasticity solution before the dislocation emission, but are not in agreement after the emission. The dislocation can move at subsonic wave speed (less than the shear wave speed) or at transonic speed (greater than the shear wave speed but less than the longitudinal wave speed), but at the longitudinal wave speed the atom lattice breaks down.

Torsional Behaviour of Cello Strings

The behaviour of a bowed string depends, among other things, on the frequency, impedance and internal damping of torsional waves on the string. Very little published information is available about these quantities, especially the torsional damping. Measurements of all relevant torsional properties have been made on cello strings of three different constructions. These show that the torsional modes are harmonically spaced to reasonable accuracy, and that the Q factors are approximately equal for all modes of a given string. These torsional Q factors are roughly an order of magnitude smaller than those of the transverse modes of the same string. The torsional wave speed varies somewhat with the tension in the string, decreasing with higher tension. The damping factors are not significantly influenced by tension. These results have been expressed in terms of a novel "reflection function" [1] suitable for direct incorporation into simulations of the bowing process.

Simulations of thermal upgrading methods for oil shale reservoirs

This paper analyzes reaction and thermal front development in porous reservoirs with reacting flows, such as those encountered in shale oil extraction. A set of dimensionless parameters and a 3D code are developed in order to investigate the important physical and chemical variables of such reservoirs when heated by in situ methods. This contribution builds on a 1D model developed for the precursor study to this work. Theory necessary for this study is presented, namely shale decomposition chemical mechanisms, governing equations for multiphase flow in porous media and necessary closure models. Plotting the ratio of the thermal wave speed to the fluid speed allows one to infer that the reaction wave front ends where this ratio is at a minimum. The reaction front follows the thermal front closely, thus allowing assumptions to be made about the extent of decomposition solely by looking at thermal wave progression. Furthermore, this sensitivity analysis showed that a certain minimum permeability is required in order to ensure the formation of a traveling thermal wave. It was found that by studying the non-dimensional governing parameters of the system one can ascribe characteristic values for these parameters for given initial and boundary conditions. This allows one to roughly predict the performance of a particular method on a particular reservoir given approximate values for initial and boundary conditions. Channelling and flow blockage due to carbon residue buildup impeded each method's performance. Blockage was found to be a result of imbalanced heating. Copyright 2012, Society of Petroleum Engineers.

Dynamic behaviors of mode III interfacial crack under a constant loading rate

In an earlier study on intersonic crack propagation, Gao et al. (J. Mech. Phys. Solids 49: 2113-2132, 2001) described molecular dynamics simulations and continuum analysis of the dynamic behaviors of a mode II dominated crack moving along a weak plane under a constant loading rate. The crack was observed to initiate its motion at a critical time after the onset of loading, at which it is rapidly accelerated to the Rayleigh wave speed and propagates at this speed for a finite time interval until an intersonic daughter crack is nucleated at a peak stress at a finite distance ahead of the original crack tip. The present article aims to analyze this behavior for a mode III crack moving along a bi-material interface subject to a constant loading rate. We begin with a crack in an initially stress-free bi-material subject to a steadily increasing stress. The crack initiates its motion at a critical time governed by the Griffith criterion. After crack initiation, two scenarios of crack propagation are investigated: the first one is that the crack moves at a constant subsonic velocity; the second one is that the crack moves at the lower shear wave speed of the two materials. In the first scenario, the shear stress ahead of the crack tip is singular with exponent -1/2, as expected; in the second scenario, the stress singularity vanishes but a peak stress is found to emerge at a distance ahead of the moving crack tip. In the latter case, a daughter crack supersonic with respect to the softer medium can be expected to emerge ahead of the initial crack once the peak stress reaches the cohesive strength of the interface.

Molecular dynamics simulation of fractures using an N-body potential

Molecular dynamics has been employed to model the fracture of a twodimensional triangular atomic lattice. The N-body Sutton-Chen potential developed for fcc metals and its extended version (Rafii-Tabar and Sutton) for fcc random binary alloys were used for the interatomic interactions. It is shown that at low temperatures cleavage fractures can occur in both an elemental metal and an alloy. At elevated temperatures the nucleation of dislocations is shown to cause a brittle-to-ductile transition. For the brittle crack propagation in the elemental metal, crack propagation speeds have been computed for different stress rates, and a crack instability found to exist as the speed reaches a critical value of about 32% of the Rayleigh wave speed. For the random alloy, we find that the dislocation movement can be affected by the distorted lattice.

Molecular dynamics simulation of fractures using an N-body potential

Molecular dynamics has been employed to model the fracture of a two dimensional triangular atomic lattice. The N-body Sutton-Chen potential developed for fcc metals and its extended version (Rafii-Tabar and Sutton) for fcc random binary alloys were used for the interatomic interactions. It is shown that at low temperatures cleavage fractures can occur in both an elemental metal and an alloy. At elevated temperatures the nucleation of dislocations is shown to cause a brittle-to-ductile transition. For the brittle crack propagation in the elemental metal, crack propagation speeds have been computed for different stress rates, and a crack instability found to exist as the speed reaches a critical value of about 32% of the Rayleigh wave speed. For the random alloy, we find that the dislocation movement can be affected by the distorted lattice.

Dust-ion-acoustic supersolitons in dusty plasmas with nonthermal electrons

Supersolitons are a recent addition to the literature on large-amplitude solitary waves in multispecies plasmas. They are distinguished from the usual solitons by their associated electric field profiles which are inherently distinct from traditional bipolar structures. In this paper, dust-ion-acoustic modes in a dusty plasma with stationary negative dust, cold fluid protons, and nonthermal electrons are investigated through a Sagdeev pseudopotential approach to see where supersolitons fit between ranges of ordinary solitons and double layers, as supersolitons always have finite amplitudes. They therefore cannot be described by reductive perturbation treatments, which rely on a weak amplitude assumption. A systematic methodology and discussion is given to distinguish the existence domains in solitary wave speed and amplitude for the different solitons, supersolitons and double layers, in terms of compositional parameters for the plasma model under consideration. © 2013 American Physical Society.

Some numerical experiments on the effect of the variation of static stability in two-layer quasi-geostrophic models

A method to solve a quasi-geostrophic two-layer model including the variation of static stability is presented. The divergent part of the wind is incorporated by means of an iterative procedure. The procedure is rather fast and the time of computation is only 60–70% longer than for the usual two-layer model. The method of solution is justified by the conservation of the difference between the gross static stability and the kinetic energy. To eliminate the side-boundary conditions the experiments have been performed on a zonal channel model. The investigation falls mainly into three parts: The first part (section 5) contains a discussion of the significance of some physically inconsistent approximations. It is shown that physical inconsistencies are rather serious and for these inconsistent models which were studied the total kinetic energy increased faster than the gross static stability. In the next part (section 6) we are studying the effect of a Jacobian difference operator which conserves the total kinetic energy. The use of this operator in two-layer models will give a slight improvement but probably does not have any practical use in short periodic forecasts. It is also shown that the energy-conservative operator will change the wave-speed in an erroneous way if the wave-number or the grid-length is large in the meridional direction. In the final part (section 7) we investigate the behaviour of baroclinic waves for some different initial states and for two energy-consistent models, one with constant and one with variable static stability. According to the linear theory the waves adjust rather rapidly in such a way that the temperature wave will lag behind the pressure wave independent of the initial configuration. Thus, both models give rise to a baroclinic development even if the initial state is quasi-barotropic. The effect of the variation of static stability is very small, qualitative differences in the development are only observed during the first 12 hours. For an amplifying wave we will get a stabilization over the troughs and an instabilization over the ridges.

Traveling waves in the Lethargic Crab Disease

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)