Numerical Investigation on Laminar Flow Due to Sudden Expansion Using Nanofluids

Laminar two-dimensional sudden expansion flow of different nanofluids is studied numerically. The governing equations are solved using stream function-vorticity method. The effect of volume fraction of the nanoparticles and type of nanoparticles on flow behaviour is examined and found significant impact. The flow response to Reynolds number in the presence of nanoparticles is examined. The presence of nanoparticles decreases the flow bifurcation Reynolds number. The size and the reattachment length of the bottom wall recirculation increase with increasing volume fraction and particle density. The effect of volume fraction and density of nanoparticles on friction factor is reported. The bottom wall recirculation strongly respond to the variation in volume faction and type of particles. However, weak response is observed for top wall recirculation.

Analysis of fatigue crack growth in RC beams

In this work, a fatigue crack propagation model developed using dimensional analysis for plain concrete is used in conjunction with the steel closing force to predict the crack growth behavior of reinforced concrete beams. A numerical procedure is followed using the proposed model to compute the fatigue life of RC beams and the dissipated energy in the steel reinforcement due to shake down behavior. Through a sensitivity study, it is found that the structural size is the most sensitive parameter on which the crack growth rate is dependent. Furthermore, the moment carrying capacity of an RC beam is computed as function of crack size by considering the effect of bond slip.

Spatial Wavelet Approach to Local Matrix Crack Detection in Composite Beams with Ply Level Material Uncertainty

Wavelet coefficients based on spatial wavelets are used as damage indicators to identify the damage location as well as the size of the damage in a laminated composite beam with localized matrix cracks. A finite element model of the composite beam is used in conjunction with a matrix crack based damage model to simulate the damaged composite beam structure. The modes of vibration of the beam are analyzed using the wavelet transform in order to identify the location and the extent of the damage by sensing the local perturbations at the damage locations. The location of the damage is identified by a sudden change in spatial distribution of wavelet coefficients. Monte Carlo Simulations (MCS) are used to investigate the effect of ply level uncertainty in composite material properties such as ply longitudinal stiffness, transverse stiffness, shear modulus and Poisson's ratio on damage detection parameter, wavelet coefficient. In this study, numerical simulations are done for single and multiple damage cases. It is observed that spatial wavelets can be used as a reliable damage detection tool for composite beams with localized matrix cracks which can result from low velocity impact damage.

Analysis of fatigue crack growth in reinforced concrete beams

Mechanical behavior of reinforced concrete members is influenced by the action of unknown crack bridging reactions of rebars. Under cyclic loading, due to progressive growth of cracks, this bridging action contributes to the overall strength, stiffness and hysteretic behavior of the member. In this work, fatigue behavior of reinforced concrete beams are studied using a crack propagation law, developed using dimensional analysis for plain concrete with the effect of reinforcement being simulated through constraint exerted on the crack opening. The parameters considered in the model are fracture toughness, crack length, loading ratio and structural size. A numerical procedure is followed to compute fatigue life of RC beams and the dissipated energy in the steel reinforcement due to the shake down phenomenon under cyclic loading. Through a sensitivity study, it is concluded that the structural size is the most sensitive parameter in the fatigue crack propagation phenomenon. Furthermore, the residual moment carrying capacity of an RC member is determined as a function of crack extension by including the bond-slip mechanism.

Energy based approach for the evaluation of damage under partial slip and gross sliding condition

Three possible contact conditions may prevail at a contact interface depending on the magnitude of normal and tangential loads, that is, stick condition, partial slip condition or gross sliding condition. Numerical techniques have been used to evaluate the stress field under partial slip and gross sliding condition. Cattaneo and Mindlin approach has been adapted to model partial slip condition. Shear strain energy density and normalized strain energy release rate have been evaluated at the surface and in the subsurface region. It is apparent from the present study that the shear strain energy density gives a fair prediction for the nucleation of damage, whereas the propagation of the crack is controlled by normalized strain energy release rate. Further, it has been observed that the intensity of damage strongly depends on coefficient of friction and contact conditions prevailing at the contact interface. (C) 2014 Elsevier B.V. All rights reserved.

Crack stability in edge-notched clamped beam specimens: modeling and experiments

Stability of a fracture toughness testing geometry is important to determine the crack trajectory and R-curve behavior of the specimen. Few configurations provide for inherent geometric stability, especially when the specimen being tested is brittle. We propose a new geometrical construction called the single edge notched clamped bend specimen (SENCB), a modified form of three point bending, yielding stable cracking under load control. It is shown to be particularly suitable for small-scale structures which cannot be made free-standing, (e.g., thin films, coatings). The SENCB is elastically clamped at the two ends to its parent material. A notch is inserted at the bottom center and loaded in bending, to fracture. Numerical simulations are carried out through extended finite element method to derive the geometrical factor f(a/W) and for different beam dimensions. Experimental corroborations of the FEM results are carried out on both micro-scale and macro-scale brittle specimens. A plot of vs a/W, is shown to rise initially and fall off, beyond a critical a/W ratio. The difference between conventional SENB and SENCB is highlighted in terms of and FEM simulated stress contours across the beam cross-section. The `s of bulk NiAl and Si determined experimentally are shown to match closely with literature values. Crack stability and R-curve effect is demonstrated in a PtNiAl bond coat sample and compared with predicted crack trajectories from the simulations. The stability of SENCB is shown for a critical range of a/W ratios, proving that it can be used to get controlled crack growth even in brittle samples under load control.

The cold mode: a phenomenological model for the evolution of density perturbations in the intracluster medium

Cool cluster cores are in global thermal equilibrium but are locally thermally unstable. We study a non-linear phenomenological model for the evolution of density perturbations in the intracluster medium (ICM) due to local thermal instability and gravity. We have analysed and extended a model for the evolution of an overdense blob in the ICM. We find two regimes in which the overdense blobs can cool to thermally stable low temperatures. One for large t(cool)/t(ff) (t(cool) is the cooling time and t(ff) is the free-fall time), where a large initial overdensity is required for thermal runaway to occur; this is the regime which was previously analysed in detail. We discover a second regime for t(cool)/t(ff) less than or similar to 1 (in agreement with Cartesian simulations of local thermal instability in an external gravitational field), where runaway cooling happens for arbitrarily small amplitudes. Numerical simulations have shown that cold gas condenses out more easily in a spherical geometry. We extend the analysis to include geometrical compression in weakly stratified atmospheres such as the ICM. With a single parameter, analogous to the mixing length, we are able to reproduce the results from numerical simulations; namely, small density perturbations lead to the condensation of extended cold filaments only if t(cool)/t(ff) less than or similar to 10.

Effect of viscosity contrast on plume formation in density stratified fluids

We perform two and three dimensional numerical simulations of plume formation in density and viscosity stratified fluid systems. We show that the ambient to plume fluid viscosity ratio strongly affects the near wall plume structures (line or sheet plumes) such as plume spacing and shape of plumes. We observe that where mushroom-like plumes are observed for lower viscosity ratios, taller plumes with bulbous heads form for high viscosity ratios. Plume structure and spacing are in good agreement with experimental results. By studying the geometry of the line plumes and the flow in the circulation cells, we discuss the mechanisms of their formation and the dynamics of merging. We show that an increase in the viscosity ratio decreases the total length of line plumes in the planform which indicates a decreased mixing at higher viscosity ratios. (C) 2015 Elsevier Ltd. All rights reserved.

Complex potentials and singular integral equation for curve crack problem in antiplane elasticity

Four types of the fundamental complex potential in antiplane elasticity are introduced: (a) a point dislocation, (b) a concentrated force, (c) a dislocation doublet and (d) a concentrated force doublet. It is proven that if the axis of the concentrated force doublet is perpendicular to the direction of the dislocation doublet, the relevant complex potentials are equivalent. Using the obtained complex potentials, a singular integral equation for the curve crack problem is introduced. Some particular features of the obtained singular integral equation are discussed, and numerical solutions and examples are given.

Numerical Analysis of Rainfall Infiltration in the Slope with a Fracture

With the finite volume method, a 2D numerical model for seepage in unsaturated soil has been established to study the rainfall infiltration in the fractured slope.The result shows that more rain may infiltrate into the slope due to existing fracture and then the pore pressure rises correspondingly. Very probably, it is one of the crucial factors accounting for slope failure. Furthermore a preliminary study has been conducted to investigate the influence of various fracture and rainfall factors such as the depth, width and location of a crack, surface condition, rainfall intensity and duration. Pore pressure and water volumetric content during the transient seepage are carefully examined to reveal the intrinsic mechanism.

Interface crack problems with strain gradient effects

In this paper, the strain gradient theory proposed by Chen and Wang (2001 a, 2002b) is used to analyze an interface crack tip field at micron scales. Numerical results show that at a distance much larger than the dislocation spacing the classical continuum plasticity is applicable; but the stress level with the strain gradient effect is significantly higher than that in classical plasticity immediately ahead of the crack tip. The singularity of stresses in the strain gradient theory is higher than that in HRR field and it slightly exceeds or equals to the square root singularity and has no relation with the material hardening exponents. Several kinds of interface crack fields are calculated and compared. The interface crack tip field between an elastic-plastic material and a rigid substrate is different from that between two elastic-plastic solids. This study provides explanations for the crack growth in materials by decohesion at the atomic scale.

InGaN/GaN LEDs grown on Si(111): Dependence of device performance on threading dislocation density and emission wavelength

We demonstrate the growth of crack-free blue and greenemitting LED structures grown on 2-inch and 6-inch Si(111) substrates by metalorganic vapour phase epitaxy (MOVPE), using AlN nucleation layers and AlGaN buffer layers for stress management. LED device performance and its dependence on threading dislocation (TD) density and emission wavelength were studied. Despite the inherently low light extraction efficiency, an output power of 1.2 mW at 50 mA was measured from a 500 μm square planar device, emitting at 455 nm. The light output decreases dramatically as the emission wavelength increases from 455 nm to 510 nm. For LED devices emitting at similar wavelength, the light output was more than doubled when the TD density was reduced from 5×1 09 cm-2 to 2×109 cm-2. Our results clearly show that high TD density is detrimental to the overall light output, highlighting the need for further TD reduction for structures grown on Si. © 2010 Wiley-VCH Verlag GmbH & Co. KGaA.

Interacting multiple cracks in piezoelectric materials

In this paper, a method is presented to calculate the plane electro-elastic fields in piezoelectric materials with multiple cracks. The cracks may be distributed randomly in locations, orientations and sizes. In the method, each crack is treated as a continuous distributed dislocations with the density function to be determined according to the conditions of external loads and crack surfaces. Some numerical examples are given to show the interacting effect among multiple cracks.

Numerical Study Of Heat Transfer Mechanism In Turbulent Supercritical Co2 Channel Flow

Direct numerical simulation (DNS) of supercritical CO2 turbulent channel flow has been performed to investigate the heat transfer mechanism of supercritical fluid. In the present DNS, full compressible Navier-Stokes equations and Peng-Robison state equation are solved. Due to effects of the mean density variation in the wall normal direction, mean velocity in the cooling region becomes high compared with that in the heating region. The mean width between high-and low-speed streaks near the wall decreases in the cooling region, which means that turbulence in the cooling region is enhanced and lots of fine scale eddies are created due to the local high Reynolds number effects. From the turbulent kinetic energy budget, it is found that compressibility effects related with pressure fluctuation and dilatation of velocity fluctuation can be ignored even for supercritical condition. However, the effect of density fluctuation on turbulent kinetic energy cannot be ignored. In the cooling region, low kinematic viscosity and high thermal conductivity in the low speed streaks modify fine scale structure and turbulent transport of temperature, which results in high Nusselt number in the cooling condition of the supercritical CO2.

Direct Numerical Simulation Of Three-Dimensional Richtmyer-Meshkov Instability

Direct numerical simulation (DNS) is used to study flow characteristics after interaction of a planar shock with a spherical media interface in each side of which the density is different. This interfacial instability is known as the Richtmyer-Meshkov (R-M) instability. The compressible Navier-Stoke equations are discretized with group velocity control (GVC) modified fourth order accurate compact difference scheme. Three-dimensional numerical simulations are performed for R-M instability installed passing a shock through a spherical interface. Based on numerical results the characteristics of 3D R-M instability are analysed. The evaluation for distortion of the interface, the deformation of the incident shock wave and effects of refraction, reflection and diffraction are presented. The effects of the interfacial instability on produced vorticity and mixing is discussed.