Stationary crack tip fields in elastic-plastic solids: an overview of recent numerical simulations

In this paper, an overview of some recent numerical simulations of stationary crack tip fields in elastic-plastic solids is presented. First, asymptotic analyses carried out within the framework of 2D plane strain or plane stress conditions in both pressure insensitive and pressure sensitive plastic solids are reviewed. This is followed by discussion of salient results obtained from recent computational studies. These pertain to 3D characteristics of elastic-plastic near-front fields under mixed mode loading, mechanics of fracture and simulation of near-tip shear banding process of amorphous alloys and influence of crack tip constraint on the structure of near-tip fields in ductile single crystals. These results serve to illustrate several important features associated with stress and strain distributions near the crack tip and provide the foundation for understanding the operative failure mechanisms. The paper concludes by highlighting some of the future prospects for this field of study.

Ruthenium(II)-CO complexes of N-[(2-pyridyl)methyliden]-alpha(or beta)-aminonaphthalene: Synthesis, spectral studies, crystal structure, redox properties and DFT calculation

The characterization and properties of trans-(X)-[RuX2(CO)(2)(alpha/beta-NaiPy)] (1, 2) (alpha-NaiPy (a), beta-NaiPy (b); X = Cl (1), I (2)) are described in this work. The structures are confirmed by single crystal X-ray diffraction studies. Reaction of these compounds with Me3NO in MeCN has isolated monocarbonyl trans-(X)-RuX2(CO)(MeCN)(alpha/beta-NaiPy)] (3, 4). The complexes show intense emission properties. Quantum yields of 1 and 2 (phi= 0.02-0.08) are higher than 3 and 4 (phi = 0.006-0.015). Voltammogram shows higher Ru(III)/Ru(II) (1.3-1.5 V) potential of 1 and 2 than that of 3 and 4 (0.8-0.9 V) that may be due to coordination of two pi-acidic CO groups in former. The electronic spectra and redox properties of the complexes are compared with the results obtained by density functional theory (DFT) and time-dependent density functional theory (TD-DFT) using polarizable continuum model (CPCM).

Numerical Study of Heat Transfer From Pin-Fin Heat Sink Using Steady and Pulsated Impinging Jets

This paper investigates numerically the heat transfer characteristics of confined slot jet impingement on a pin-fin heat sink. A variety of pin-fin heat sinks is investigated, and the resulting enhancement of heat transfer studied. The distribution of heat transfer coefficient on the top surface of the base plate and that along the fin height are examined. Both steady and pulsated jets are studied. It is observed that for a steady jet impingement on a pin-fin heat sink, the effective heat transfer coefficient increases with fin height, leading to a corresponding decrease in base plate temperature for the same heat flux. In the case of pulsated jets, the influence of pulse frequency and the Reynolds number is examined, and their effect on the effective heat transfer coefficient is studied.

A numerical study of projectile impact on thin aluminium plates

This article deals with a simulation-based Study of the impact of projectiles on thin aluminium plates using LS-DYNA by modelling plates with shell elements and projectiles with solid elements. In order to establish the required modelling criterion in terms of element size for aluminium plates, a convergence Study of residual velocity has been carried Out by varying mesh density in the impact zone. Using the preferred material and meshing criteria arrived at here, extremely good prediction of test residual velocities and ballistic limits given by Gupta et al. (2001) for thin aluminium plates has been obtained. The simulation-based pattern of failure with localized bulging and jagged edge of perforation is similar to the perforation with petalling seen in tests. A number Of simulation-based parametric studies have been carried out and results consistent with published test data have been obtained. Despite the robust correlation achieved against published experimental results, it would be prudent to conduct one's own experiments, for a final correlation via the present modelling procedure and analysis with the explicit LS-DYNTA 970 solver. Hence, a sophisticated ballistic impact testing facility and a high-speed camera have been used to conduct additional tests on grade 1100 aluminium plates of 1 mm thickness with projectiles Of four different nose shapes. Finally, using the developed numerical simulation procedure, an excellent correlation of residual velocity and failure modes with the corresponding test results has been obtained.

Structural Investigation of Activated Lattice Oxygen in Ce1-xSnxO2 and Ce1-x-ySnxPdyO2-delta by EXAFS and DFT calculation

Substitution of Sn4+ ion in CeO2 creates activated oxygen in Ce0.8Sn0.2O2 leading to higher oxygen storage capacity compared to Ce0.8Zr0.2O2. With Pd ion substitution in Ce0.8Sn0.2O2,activation of oxygen is further enhanced as observed from the H-2/TPR study. Both EXAFS analysis and DFT calculation reveal that in the solid solution Ceexhibits 4 + 4 coordination, Sri exhibits 4 + 2 + 2 coordination and Pd has 4 + 3 coordination. While the oxygen in the First four coordination with short M-O bonds are strongly held in the lattice, the oxygens in the second and higher coordinations with long M-O bonds are weakly bound, and they are the activated oxygen ill the lattice. Bond valence analysis shows that oxygen with valencies as low its 1.65 are created by the Sn and Pd ion Substitution. Another interesting observation is that H-2/TPR experiment of Ce1-xSnxO2 shows a broad peak starting from 200 to 500 degrees C, while the same reduction is achieved in a single step at similar to 110 degrees C in presence Pd2+ on. Substitution of Pd2+ ion thus facilitates synergistic reduction of the catalyst at lower temperature. We have shown that simultaneous reduction of the Ce4+ and Sr4+ ions by Pd-0 is the synergistic interaction leading to high oxygen storage capacity at low temperature.

A numerical study of crack tip constraint in ductile single crystals

In this work, the effect of crack tip constraint on near-tip stress and deformation fields in a ductile FCC single crystal is studied under mode I, plane strain conditions. To this end, modified boundary layer simulations within crystal plasticity framework are performed, neglecting elastic anisotropy. The first and second terms of the isotropic elastic crack tip field, which are governed by the stress intensity factor K and T-stress, are prescribed as remote boundary conditions and solutions pertaining to different levels of T-stress are generated. It is found that the near-tip deformation field, especially, the development of kink or slip shear bands, is sensitive to the constraint level. The stress distribution and the size and shape of the plastic zone near the crack tip are also strongly influenced by the level of T-stress, with progressive loss of crack tip constraint occurring as T-stress becomes more negative. A family of near-tip fields is obtained which are characterized by two terms (such as K and T or J and a constraint parameter Q) as in isotropic plastic solids.

A numerical study of the effect of loading conditions on tubular hydroforming

In this work, the mechanics of tubular hydroforming under various types of loading conditions is investigated. The main objective is to contrast the effects of prescribing fluid pressure or volume flow rate, in conjunction with axial displacement, on the stress and strain histories experienced by the tube and the process of bulging. To this end, axisymmetric finite element simulations of free hydroforming (without external die contact) of aluminium alloy tubes are carried out. Hill’s normally anisotropic yield theory along with material properties determined in a previous experimental study [A. Kulkarni, P. Biswas, R. Narasimhan, A. Luo, T. Stoughton, R. Mishra, A.K. Sachdev, An experimental and numerical study of necking initiation in aluminium alloy tubes during hydroforming, Int. J. Mech. Sci. 46 (2004) 1727–1746] are employed in the computations. It is found that while prescribed fluid pressure leads to highly non-proportional strain paths, specified fluid volume flow rate may result in almost proportional ones for the predominant portion of loading. The peak pressure increases with axial compression for the former, while the reverse trend applies under the latter. The implication of these results on failure by localized necking of the tube wall is addressed in a subsequent investigation.

Ru(II)-CO complexes of thioarylazoimidazoles: Syntheses, structures, spectroscopy and DFT calculation

The complexes, cis-(CO)-trans-(Cl)-[Ru(SRaaiNR)(CO)(2)Cl-2] (2) and trans-(Cl)-[Ru(SRaaiNR)(CO)Cl-2] (3) (SRaaiNR = 1-alkyl-2-{(o-thioalkyl)phenylazo}imidazoles; R = Me (1a) and Et (1b)) have been synthesized and characterized. The structural confirmation is achieved by single crystal X-ray structure determinations. The complexes show Ru(III)/Ru(II) couple and ligand reductions. Electronic structure and spectral properties of the complexes have been explained with the DFT and TDDFT calculation. (C) 2009 Elsevier B.V. All rights reserved.

Numerical Simulation of Groundwater Flow with Gradually Increasing Heterogeneity due to Clogging

Well injection replenishes depleting water levels in a well field. Observation well water levels some distance away from the injection well are the indicators of the success of a well injection program. Simulation of the observation well response, located a few tens of meters from the injection well, is likely to be affected by the effects of nonhomogeneous medium, inclined initial water table, and aquifer clogging. Existing algorithms, such as the U.S. Geological Survey groundwater flow software MODFLOW, are capable of handling the first two conditions, whereas time-dependent clogging effects are yet to be introduced in the groundwater flow models. Elsewhere, aquifer clogging is extensively researched in theory of filtration; scope for its application in a well field is a potential research problem. In the present paper, coupling of one such filtration theory to MODFLOW is introduced. Simulation of clogging effects during “Hansol” well recharge in the parts of western India is found to be encouraging.

Numerical studies on channel formation and growth during solidification: Effect of process parameters

In the present work, solidification of a hyper-eutectic ammonium chloride solution in a bottom-cooled cavity (i.e. with stable thermal gradient) is numerically studied. A Rayleigh number based criterion is developed, which determines the conditions favorable for freckles formation. This criterion, when expressed in terms of physical properties and process parameters, yields the condition for plume formation as a function of concentration, liquid fraction, permeability, growth rate of a mushy layer and thermophysical properties. Subsequently, numerical simulations are performed for cases with initial and boundary conditions favoring freckle formation. The effects of parameters, such as cooling rate and initial concentration, on the formation and growth of freckles are investigated. It was found that a high cooling rate produced larger and more defined channels which are retained for a longer durations. Similarly, a lower initial concentration of solute resulted in fewer but more pronounced channels. The number and size of channels are also found to be related to the mushy zone thickness. The trends predicted with regard to the variation of number of channels with time under different process conditions are in accordance with the experimental observations reported in the literature.

Numerical procedure for second order non-linear ordinary differential equations and application to heat transfer problem

In this paper, we have first given a numerical procedure for the solution of second order non-linear ordinary differential equations of the type y″ = f (x;y, y′) with given initial conditions. The method is based on geometrical interpretation of the equation, which suggests a simple geometrical construction of the integral curve. We then translate this geometrical method to the numerical procedure adaptable to desk calculators and digital computers. We have studied the efficacy of this method with the help of an illustrative example with known exact solution. We have also compared it with Runge-Kutta method. We have then applied this method to a physical problem, namely, the study of the temperature distribution in a semi-infinite solid homogeneous medium for temperature-dependent conductivity coefficient.

Numerical simulations of crack tip fields in polycrystalline plastic solids

In this paper, modes I and II crack tip fields in polycrystalline plastic solids are studied under plane strain, small scale yielding conditions. Two different initial textures of an Al–Mg alloy, viz., continuous cast AA5754 sheets in the recrystallized and cold rolled conditions, are considered. The former is nearly-isotropic, while the latter displays distinct anisotropy. Finite element simulations are performed by employing crystal plasticity constitutive equations along with a Taylor-type homogenization as well as by using the Hill quadratic yield theory. It is found that significant texture evolution occurs close to the notch tip which profoundly influences the stress and plastic strain distributions. Also, the cold rolling texture gives rise to higher magnitude of plastic strain near the tip.

Numerical simulations of crack tip fields in polycrystalline plastic solids

In this paper, modes I and II crack tip fields in polycrystalline plastic solids are studied under plane strain, small scale yielding conditions. Two different initial textures of an Al-Mg alloy, viz.,continuous cast AA5754 sheets in the recrystallized and cold rolled conditions, are considered. The former is nearly-isotropic, while the latter displays distinct anisotropy. Finite element simulations are performed by employing crystal plasticity constitutive equations along with a Taylor-type homogenization as well as by using the Hill quadratic yield theory. It is found that significant texture evolution occurs close to the notch tip which profoundly influences the stress and plastic strain distributions. Also, the cold rolling texture gives rise to higher magnitude of plastic strain near the tip. (C) 2010 Elsevier Ltd. All rights reserved.

A fast algorithm for macromolecular packing calculation

An algorithm to improve the computation time of packing calculations for macromolecules is presented. This is achieved by reducing the three-dimensional search to a small set of two-dimensional searches.