Effect of surfactants on the instability of a two-layer film flow down an inclined plane

The effect of insoluble surfactants on the instability of a two-layer film flow down an inclined plane is investigated based on the Orr-Sommerfeld boundary value problem. The study, focusing on Stokes flow P. Gao and X.-Y. Lu, ``Effect of surfactants on the inertialess instability of a two-layer film flow,'' J. Fluid Mech. 591, 495-507 (2007)], is further extended by including the inertial effect. The surface mode is recognized along with the interface mode. The initial growth rate corresponding to the interface mode accelerates at sufficiently long-wave regime in the presence of surface surfactant. However, the maximum growth rate corresponding to both interface and surface modes decelerates in the presence of surface surfactant when the upper layer is more viscous than the lower layer. On the other hand, when the upper layer is less viscous than the lower layer, a new interfacial instability develops due to the inertial effect and becomes weaker in the presence of interfacial surfactant. In the limit of negligible surface and interfacial tensions, respectively, two successive peaks of temporal growth rate appear in the long-wave and short-wave regimes when the interface mode is analyzed. However, in the case of the surface mode, only the long-wave peak appears. (C) 2014 AIP Publishing LLC.

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.

A general methodology for calculating mixed mode stress intensity factors and fracture toughness of solder joints with interfacial cracks

Solder joints in electronic packages undergo thermo-mechanical cycling, resulting in nucleation of micro-cracks, especially at the solder/bond-pad interface, which may lead to fracture of the joints. The fracture toughness of a solder joint depends on material properties, process conditions and service history, as well as strain rate and mode-mixity. This paper reports on a methodology for determining the mixed-mode fracture toughness of solder joints with an interfacial starter-crack, using a modified compact mixed mode (CMM) specimen containing an adhesive joint. Expressions for stress intensity factor (K) and strain energy release rate (G) are developed, using a combination of experiments and finite element (FE) analysis. In this methodology, crack length dependent geometry factors to convert for the modified CMM sample are first obtained via the crack-tip opening displacement (CTOD)-based linear extrapolation method to calculate the under far-field mode I and II conditions (f(1a) and f(2a)), (ii) generation of a master-plot to determine a(c), and (iii) computation of K and G to analyze the fracture behavior of joints. The developed methodology was verified using J-integral calculations, and was also used to calculate experimental fracture toughness values of a few lead-free solder-Cu joints. (C) 2014 Elsevier Ltd. All rights reserved.

Two techniques to improve the precision of a demand-driven null-dereference verification approach

The problem addressed in this paper is sound, scalable, demand-driven null-dereference verification for Java programs. Our approach consists conceptually of a base analysis, plus two major extensions for enhanced precision. The base analysis is a dataflow analysis wherein we propagate formulas in the backward direction from a given dereference, and compute a necessary condition at the entry of the program for the dereference to be potentially unsafe. The extensions are motivated by the presence of certain ``difficult'' constructs in real programs, e.g., virtual calls with too many candidate targets, and library method calls, which happen to need excessive analysis time to be analyzed fully. The base analysis is hence configured to skip such a difficult construct when it is encountered by dropping all information that has been tracked so far that could potentially be affected by the construct. Our extensions are essentially more precise ways to account for the effect of these constructs on information that is being tracked, without requiring full analysis of these constructs. The first extension is a novel scheme to transmit formulas along certain kinds of def-use edges, while the second extension is based on using manually constructed backward-direction summary functions of library methods. We have implemented our approach, and applied it on a set of real-life benchmarks. The base analysis is on average able to declare about 84% of dereferences in each benchmark as safe, while the two extensions push this number up to 91%. (C) 2014 Elsevier B.V. All rights reserved.

Influence of external weld flash on the in-plane plane-strain formability of friction stir welded sheets

The main aim of the present work is to analyze the influence of external weld flash on the formability of friction stir welding sheets through in-plane plane-strain formability tests. The load-extension behavior and forming limit strains are measured to quantify the formability. The influence of friction stir welding parameters on the height of weld flash was also studied. The base materials used for welding trials are AA6061T6 and AA5052H32 alloy sheets of 2.1-mm thickness. It is observed that the influence of external weld flash on the maximum load and total extension for all the friction stir welding conditions is negligible. The effect of weld flash on the limiting major strain is also insignificant. But the presence of weld flash has changed the limiting minor strain, more toward plane-strain condition, indicating the change in strain-path toward plane-strain. This is due to the strain taken by weld flash, along with the major strain, minor strain, and thickness strain in the friction stir welding sheet plane because of constancy of volume. The formation of weld flash and its height are affected synergistically by the axial force and temperature development during friction stir welding. The higher the axial force and temperature, the higher the flash height.

Crack propagation in graphene

The crack initiation and growth mechanisms in an 2D graphene lattice structure are studied based on molecular dynamics simulations. Crack growth in an initial edge crack model in the arm-chair and the zig-zag lattice configurations of graphene are considered. Influence of the time steps on the post yielding behaviour of graphene is studied. Based on the results, a time step of 0.1 fs is recommended for consistent and accurate simulation of crack propagation. Effect of temperature on the crack propagation in graphene is also studied, considering adiabatic and isothermal conditions. Total energy and stress fields are analyzed. A systematic study of the bond stretching and bond reorientation phenomena is performed, which shows that the crack propagates after significant bond elongation and rotation in graphene. Variation of the crack speed with the change in crack length is estimated. (C) 2015 AIP Publishing LLC.

Effects of instabilities and coherent structures on the performance of a thermocline based thermal energy storage

We investigate the transient dynamics of disturbances inside a thermocline based molten salt thermal energy storage (TES). Numerical simulations were conducted with four inlet flow configurations. The disturbances introduced at the inlet grow via Rayleigh Taylor instability. The formed vortical motions inside the tank propagate downstream and destroy the thermocline. The vortex-thermocline interaction upsets the stratification inside the TES. The disturbance growth rate, penetration length and vortex Reynolds number are measured. The growth of penetration length prior to the vortex-thermocline interaction is quadratic. The vortex Reynolds number of the eddy which causes thermocline breakdown increases with increase in Atwood number. The impingement of vortex on thermocline is studied. (C) 2015 Elsevier Ltd. All rights reserved.

Morphological and functional characterisation of the burrow systems of six earthworm species (Lumbricidae)

Earthworm burrow systems are generally described based on postulated behaviours associated with the three ecological types. In this study, we used X-ray tomography to obtain 3D information on the burrowing behaviour of six very common anecic (Aporrectodea nocturna and Lumbricus terrestris) and endogeic (Aporrectodea rosea, Allolobophora chlorotica, Aporrectodea caliginosa, Aporrectodea icterica) earthworm species, introduced into repacked soil cores for 6 weeks. A simple water infiltration test, the Beerkan method, was also used to assess some functional properties of these burrow systems. Endogeic worms make larger burrow systems, which are more highly branched, less continuous and of smaller diameter, than those of anecic worms. Among the anecic species, L. terrestris burrow systems are shorter (9.2 vs 21.2 m) with a higher number (14.5 vs 23.5) of less branched burrows (12.2 vs 20.2 branches m(-1)), which are also wider (7.78 vs 5.16 mm) than those of A. nocturna. In comparison, the burrow systems made by endogeic species appeared similar to each other. However, A. rosea burrows were short and narrow, whereas A. icterica had a longer burrow system (15.7 m), more intense bioturbation intensity (refilled macropores or soil lateral compaction around them) and thus a greater number of burrows. Regarding water infiltration, anecic burrow systems were far more efficient due to open burrows linking the top and bottom of the cores. For endogeic species, we observed a linear relationship between burrow length and the water infiltration rate (R (2) = 0.49, p < 0.01). Overall, the three main characteristics significantly influencing water infiltration were burrow length, burrow number and bioturbation volume. This last characteristic highlighted the effect of burrow refilling by casts.

DOES THE CHEMOTHERMAL INSTABILITY HAVE ANY ROLE IN THE FRAGMENTATION OF PRIMORDIAL GAS

The collapse of the primordial gas in the density regime similar to 10(8)-10(10) cm(-3) is controlled by the three-body H-2 formation process, in which the gas can cool faster than free-fall time-a condition proposed as the chemothermal instability. We investigate how the heating and cooling rates are affected during the rapid transformation of atomic to molecular hydrogen. With a detailed study of the heating and cooling balance in a 3D simulation of Pop III collapse, we follow the chemical and thermal evolution of the primordial gas in two dark matter minihalos. The inclusion of sink particles in modified Gadget-2 smoothed particle hydrodynamics code allows us to investigate the long-term evolution of the disk that fragments into several clumps. We find that the sum of all the cooling rates is less than the total heating rate after including the contribution from the compressional heating (pdV). The increasing cooling rate during the rapid increase of the molecular fraction is offset by the unavoidable heating due to gas contraction. We conclude that fragmentation occurs because H-2 cooling, the heating due to H-2 formation and compressional heating together set a density and temperature structure in the disk that favors fragmentation, not the chemothermal instability.

Estimation of activation energy for electroporation and pore growth rate in liquid crystalline and gel phases of lipid bilayers using molecular dynamics simulations

Molecular dynamics simulations of electroporation in POPC and DPPC lipid bilayers have been carried out at different temperatures ranging from 230 K to 350 K for varying electric fields. The dynamics of pore formation, including threshold field, pore initiation time, pore growth rate, and pore closure rate after the field is switched off, was studied in both the gel and liquid crystalline (L-alpha) phases of the bilayers. Using an Arrhenius model of pore initiation kinetics, the activation energy for pore opening was estimated to be 25.6 kJ mol(-1) and 32.6 kJ mol(-1) in the L-alpha phase of POPC and DPPC lipids respectively at a field strength of 0.32 V nm(-1). The activation energy decreases to 24.2 kJ mol(-1) and 23.7 kJ mol(-1) respectively at a higher field strength of 1.1 V nm(-1). At temperatures below the melting point, the activation energy in the gel phase of POPC and DPPC increases to 28.8 kJ mol(-1) and 34.4 kJ mol(-1) respectively at the same field of 1.1 V nm(-1). The pore closing time was found to be higher in the gel than in the L-alpha phase. The pore growth rate increases linearly with temperature and quadratically with field, consistent with viscosity limited growth models.

Sum Secrecy Rate in Multicarrier DF Relay Beamforming

In this paper, we study sum secrecy rate in multicarrier decode-and-forward relay beamforming. We obtain the optimal source power and relay weights on each subcarrier which maximize the sum secrecy rate. For a given total power on a given subcarrier k, P-0(k), we reformulate the optimization problem by relaxing the rank-1 constraint on the complex positive semidefinite relay weight matrix, and solve using semidefinite programming. We analytically prove that the solution to the relaxed optimization problem is indeed rank 1. We show that the subcarrier secrecy rate, R-s (P-0(k)), is a concave function in total power P-0(k) if R-s (P-0(k)) > 0 for any P-0(k) > 0. Numerical results show that the sum secrecy rate with optimal power allocation across subcarriers is more than the sum secrecy rate with equal power allocation. We also propose a low complexity suboptimal power allocation scheme which outperforms equal power allocation scheme.

Study of H/D exchange rates to derive the strength of intramolecular hydrogen bonds in halo substituted organic building blocks: An NMR spectroscopic investigation

Rates of hydrogen/deuterium (H/D) exchange determined by H-1 NMR spectroscopy are utilized to derive the strength of hydrogen bonds and to monitor the electronic effects in the site-specific halogen substituted benzamides and anilines. The theoretical fitting of the time dependent variation of the integral areas of H-1 NMR resonances to the first order decay function permitted the determination of HID exchange rate constants (k) and their precise half-lives (t(1/2)) with high degree of reproducibility. The comparative study also permitted the unambiguous determination of relative strength of hydrogen bonds and the contribution from electronic effects on the HID exchange rate. (C) 2015 Elsevier B.V. All rights reserved.

Role of Si modification on the compressive flow behavior of Al-Si based alloy

The flow characteristics of a near-eutectic heat-treated Al-Si based cast alloy have been examined in compression at strain rates varying from 3 x 10(-4) to 10(2) s(-1) and at three different temperatures, i.e., room temperature (RT), 100 degrees C and 200 degrees C. The dependence of flow behavior on modification is examined by testing the alloy in both the unmodified and modified conditions. Modification has strong influence on strain rate sensitivity (SRS), strength and work hardening behavior of the alloy. The strength of the alloy is found to increase with increase in strain rate for both the conditions. The increase is more rapid above the strain rate of 10(-1) s(-1) for the unmodified alloy at all the temperatures. This rapid increase is observed at 1 s(-1) at RT and 100 degrees C, and at 10(-2) s(-1) at 200 degrees C for the modified alloy. The thermally dependent process of the Al matrix is rate controlling in the unmodified alloy. On the other hand, the thermally dependent process of both Al matrix and Si particles are rate controlling, which is responsible for the higher strain rate sensitivity (SRS) in the modified alloy. The unmodified alloy exhibits a larger work hardening rate than the modified alloy during the initial stages of straining due to fiber loading of unmodified Si particles. However, the hardening rate decreases sharply at higher strains for the unmodified alloy due to a higher rate of Si particle fracture. Thermal softening is observed for both alloys at 200 degrees C due to precipitate coarsening, which leads to a decrease in SRS at higher temperatures. Stress simulations by microstructure based finite element method support the experimentally observed particle and matrix fracture behavior. Negative SRS and serrated flow are observed at lower strain rate regime (3 x 10(-4) to 10(-2) s(-1)) at RT and 100 degrees C, in both alloys. The critical onset strain is found to be lower and the magnitude of serration is found to be higher for the modified alloy, which suggests that, in addition to dynamic strain aging, Si particle size and morphology also play a role in serrated flow. (C) 2015 Elsevier Inc All rights reserved.

Divergence in Patterns of Leaf Growth Polarity Is Associated with the Expression Divergence of miR396

Lateral appendages often show allometric growth with a specific growth polarity along the proximo-distal axis. Studies on leaf growth in model plants have identified a basipetal growth direction with the highest growth rate at the proximal end and progressively lower rates toward the distal end. Although the molecular mechanisms governing such a growth pattern have been studied recently, variation in leaf growth polarity and, therefore, its evolutionary origin remain unknown. By surveying 75 eudicot species, here we report that leaf growth polarity is divergent. Leaf growth in the proximo-distal axis is polar, with more growth arising from either the proximal or the distal end; dispersed with no apparent polarity; or bidirectional, with more growth contributed by the central region and less growth at either end. We further demonstrate that the expression gradient of the miR396-GROWTH-REGULATING FACTOR module strongly correlates with the polarity of leaf growth. Altering the endogenous pattern of miR396 expression in transgenic Arabidopsis thaliana leaves only partially modified the spatial pattern of cell expansion, suggesting that the diverse growth polarities might have evolved via concerted changes in multiple gene regulatory networks.

Effect of Ni content on the diffusion-controlled growth of the product phases in the Cu(Ni)-Sn system

A strong influence of Ni content on the diffusion-controlled growth of the (Cu,Ni)(3)Sn and (Cu,Ni)(6)Sn-5 phases by coupling different Cu(Ni) alloys with Sn in the solid state is reported. The continuous increase in the thickness ratio of (Cu,Ni)(6)Sn-5 to (Cu,Ni)(3)Sn with the Ni content is explained by combined kinetic and thermodynamic arguments as follows: (i) The integrated interdiffusion coefficient does not change for the (Cu,Ni)(3)Sn phase up to 2.5 at.% Ni and decreases drastically for 5 at.% Ni. On the other hand, there is a continuous increase in the integrated interdiffusion coefficient for (Cu,Ni)(6)Sn-5 as a function of increasing Ni content. (ii) With the increase in Ni content, driving forces for the diffusion of components increase for both components in both phases but at different rates. However, the magnitude of these changes alone is not large enough to explain the high difference in the observed growth rate of the product phases because of Ni addition. (iv) Kirkendall marker experiments indicate that the Cu6Sn5 phase grows by diffusion of both Cu and Sn in the binary case. However, when Ni is added, the growth is by diffusion of Sn only. (v) Also, the observed grain refinement in the Cu6Sn5 phase with the addition of Ni suggests that the grain boundary diffusion of Sn may have an important role in the observed changes in the growth rate.