Numerical modeling of the non-isothermal liquid droplet impact on a hot solid substrate

The heat transfer from a solid phase to an impinging non-isothermal liquid droplet is studied numerically. A new approach based on an arbitrary Lagrangian-Eulerian (ALE) finite element method for solving the incompressible Navier Stokes equations in the liquid and the energy equation within the solid and the liquid is presented. The novelty of the method consists in using the ALE-formulation also in the solid phase to guarantee matching grids along the liquid solid interface. Moreover, a new technique is developed to compute the heat flux without differentiating the numerical solution. The free surface and the liquid solid interface of the droplet are represented by a moving mesh which can handle jumps in the material parameter and a temperature dependent surface tension. Further, the application of the Laplace-Beltrami operator technique for the curvature approximation allows a natural inclusion of the contact angle. Numerical simulation for varying Reynold, Weber, Peclet and Biot numbers are performed to demonstrate the capabilities of the new approach. (C) 2014 Elsevier Ltd. All rights reserved.

An improved free surface modeling for incompressible SPH

In Incompressible Smooth Particle Hydrodynamics (ISPH), a pressure Poisson equation (PPE) is solved to obtain a divergence free velocity field. When free surfaces are simulated using this method a Dirichlet boundary condition for pressure at the free surface has to be applied. In existing ISPH methods this is achieved by identifying free surface particles using heuristically chosen threshold of a parameter such as kernel sum, density or divergence of the position, and explicitly setting their pressure values. This often leads to clumping of particles near the free surface and spraying off of surface particles during splashes. Moreover, surface pressure gradients in flows where surface tension is important are not captured well using this approach. We propose a more accurate semi-analytical approach to impose Dirichlet boundary conditions on the free surface. We show the efficacy of the proposed algorithm by using test cases of elongation of a droplet and dam break. We perform two dimensional simulations of water entry and validate the proposed algorithm with experimental results. Further, a three dimensional simulation of droplet splash is shown to compare well with the Volume-of-Fluid simulations. (C) 2014 Elsevier Ltd. All rights reserved.

Analysis of cracked and un-cracked semicircular rings under symmetric loading

Edge cracked specimens have been widely utilized for fracture testing. Edge cracked semicircular disk (ECSD) specimen has now been well characterized with regard to its form factor and weight function. This paper presents a modified semicircular ring version of this specimen to enhance the form factor in general while retaining other desirable features. The efficacy of the modified design is proved by combining theory of elasticity solutions with finite element results to arrive at the optimum design geometry. New insights emerging from this work are used to theoretically re-examine the arch-tension and the four-point bend specimens. (C) 2014 Elsevier Ltd. All rights reserved.

STM verification of the reduction of the Young's modulus of CdS nanoparticles at smaller sizes

We demonstrate the first STM evaluation of the Young's modulus (E) of nanoparticles (NPs) of different sizes. The sample deformation induced by tip-sample interaction has been determined using current-distance (I-Z) spectroscopy. As a result of tip-sample interaction, and the induced surface deformations, the I-z curves deviates from pure exponential dependence. Normally, in order to analyze the deformation quantitatively, the tip radius must be known. We show, that this necessity is eliminated by measuring the deformation on a substrate with a known Young's modulus (Au(111)) and estimating the tip radius, and afterwards, using the same tip (with a known radius) to measure the (unknown) Young's modulus of another sample (nanoparticles of CdS). The Young's modulus values found for 3 NP's samples of average diameters of 3.7, 6 and 7.5 nm, were E similar to 73%, 78% and 88% of the bulk value, respectively. These results are in a good agreement with the theoretically predicted reduction of the Young's modulus due to the changes in hydrostatic stresses which resulted from surface tension in nanoparticles with different sizes. Our calculation using third order elastic constants gives a reduction of E which scales linearly with 1/r (r is the NP's radius). This demonstrates the applicability of scanning tunneling spectroscopy for local mechanical characterization of nanoobjects. The method does not include a direct measurement of the tip-sample force but is rather based on the study of the relative elastic response. (C) 2014 Elsevier B.V. All rights reserved.

Bearing capacity of circular foundations reinforced with geogrid sheets

The ultimate bearing capacity of a circular footing, placed over a soil mass which is reinforced with horizontal layers of circular reinforcement sheets, has been determined by using the upper bound theorem of the limit analysis in conjunction with finite elements and linear optimization. For performing the analysis, three different soil media have been separately considered, namely, (i) fully granular, (ii) cohesive frictional, and (iii) fully cohesive with an additional provision to account for an increase of cohesion with depth. The reinforcement sheets are assumed to be structurally strong to resist axial tension but without having any resistance to bending; such an approximation usually holds good for geogrid sheets. The shear failure between the reinforcement sheet and adjoining soil mass has been considered. The increase in the magnitudes of the bearing capacity factors (N-c and N-gamma) with an inclusion of the reinforcement has been computed in terms of the efficiency factors eta(c) and eta(gamma). The results have been obtained (i) for different values of phi in case of fully granular (c=0) and c-phi soils, and (ii) for different rates (m) at which the cohesion increases with depth for a purely cohesive soil (phi=0 degrees). The critical positions and corresponding optimum diameter of the reinforcement sheets, for achieving the maximum bearing capacity, have also been established. The increase in the bearing capacity with an employment of the reinforcement increases continuously with an increase in phi. The improvement in the bearing capacity becomes quite extensive for two layers of the reinforcements as compared to the single layer of the reinforcement. The results obtained from the study are found to compare well with the available theoretical and experimental data reported in literature. (C) 2014 The Japanese Geotechnical Society. Production and hosting by Elsevier B.V. All rights reserved.

Bearing Capacity of Circular Footings on Reinforced Soils

A method is presented for determining the ultimate bearing capacity of a circular footing reinforced with a horizontal circular sheet of reinforcement placed over granular and cohesive-frictional soils. It was assumed that the reinforcement sheet could bear axial tension but not the bending moment. The analysis was performed based on the lower-bound theorem of the limit analysis in combination with finite elements and linear optimization. The present research is an extension of recent work with strip foundations reinforced with different layers of reinforcement. To incorporate the effect of the reinforcement, the efficiency factors eta(gamma) and eta(c), which need to be multiplied by the bearing capacity factors N-gamma and N-c, were established. Results were obtained for different values of the soil internal friction angle (phi). The optimal positions of the reinforcements, which would lead to a maximum improvement in the bearing capacity, were also determined. The variations of the axial tensile force in the reinforcement sheet at different radial distances from the center were also studied. The results of the analysis were compared with those available from literature. (C) 2014 American Society of Civil Engineers.

Wave propagation analysis in adhesively bonded composite joints using the wavelet spectral finite element method

A wavelet spectral finite element (WSFE) model is developed for studying transient dynamics and wave propagation in adhesively bonded composite joints. The adherands are formulated as shear deformable beams using the first order shear deformation theory (FSDT) to obtain accurate results for high frequency wave propagation. Equations of motion governing wave motion in the bonded beams are derived using Hamilton's principle. The adhesive layer is modeled as a line of continuously distributed tension/compression and shear springs. Daubechies compactly supported wavelet scaling functions are used to transform the governing partial differential equations from time domain to frequency domain. The dynamic stiffness matrix is derived under the spectral finite element framework relating the nodal forces and displacements in the transformed frequency domain. Time domain results for wave propagation in a lap joint are validated with conventional finite element simulations using Abaqus. Frequency domain spectrum and dispersion relation results are presented and discussed. The developed WSFE model yields efficient and accurate analysis of wave propagation in adhesively-bonded composite joints. (C) 2014 Elsevier Ltd. All rights reserved.

The contribution of grain boundary sliding in tensile deformation of an ultrafine-grained aluminum alloy having high strength and high ductility

An as-cast Al-7 % Si alloy was processed by high-pressure torsion (HPT) for up to 10 turns at temperatures of 298 or 445 K. The HPT-processed samples had ultrafine-grained structures and they were tested in tension at room temperature at various strain rates in the range from 1.0 x 10(-4) to 1.0 x 10(-2) s(-1). The contributions of grain boundary sliding (GBS) to the total strain were measured directly using atomic force microscopy. Samples simultaneously showing both high strength and high ductility contained the highest fractions of high-angle grain boundaries (HAGB) and exhibited the highest contributions from GBS, whereas samples showing high strength but low ductility gave negligible values for the sliding contributions. It is concluded that high strength and high ductility require both an ultrafine grain size and a high fraction of HAGB.

Non-uniform beams and stiff strings isospectral to axially loaded uniform beams and piano strings

In this paper, we derive analytical expressions for mass and stiffness functions of transversely vibrating clamped-clamped non-uniform beams under no axial loads, which are isospectral to a given uniform axially loaded beam. Examples of such axially loaded beams are beam columns (compressive axial load) and piano strings (tensile axial load). The Barcilon-Gottlieb transformation is invoked to transform the non-uniform beam equation into the axially loaded uniform beam equation. The coupled ODEs involved in this transformation are solved for two specific cases (pq (z) = k (0) and q = q (0)), and analytical solutions for mass and stiffness are obtained. Examples of beams having a rectangular cross section are shown as a practical application of the analysis. Some non-uniform beams are found whose frequencies are known exactly since uniform axially loaded beams with clamped ends have closed-form solutions. In addition, we show that the tension required in a stiff piano string with hinged ends can be adjusted by changing the mass and stiffness functions of a stiff string, retaining its natural frequencies.

Water-responsive carbon nanotubes for selective detection of toxic gases

Ammonia plays an important role in our daily lives and hence its quantitative and qualitative sensing has become necessary. Bulk structure of carbon nanotubes (CNTs) has been employed to detect the gas concentration of 10 ppm. Hydrophobic CNTs were turned to hydrophilic via the application of a ramp electric field that allowed confinement of a controlled amount of water inside CNT microstructure. These samples were then also used to detect different gases. A comparative study has been performed for sensing three reducing gases, namely, ammonia, sulphur-di-oxide, and hydrogen sulphide to elaborate the selectivity of the sensor. A considerable structural bending in the bulk CNT was observed on evaporation of the confined water, which can be accounted to the zipping of individual nanotubes. However, the rate of the stress induced on these bulk microstructures increased on the exposure of ammonia due to the change in the surface tension of the confined solvent. A prototype of an alarm system has been developed to illustrate sensing concept, wherein the generated stress in the bulk CNT induces a reversible loss in electrical contact that changes the equivalent resistance of the electrical circuit upon exposure to the gas. (C) 2015 AIP Publishing LLC.

THREE-LAYER FLUID FLOW OVER A SMALL OBSTRUCTION ON THE BOTTOM OF A CHANNEL

Many boundary value problems occur in a natural way while studying fluid flow problems in a channel. The solutions of two such boundary value problems are obtained and analysed in the context of flow problems involving three layers of fluids of different constant densities in a channel, associated with an impermeable bottom that has a small undulation. The top surface of the channel is either bounded by a rigid lid or free to the atmosphere. The fluid in each layer is assumed to be inviscid and incompressible, and the flow is irrotational and two-dimensional. Only waves that are stationary with respect to the bottom profile are considered in this paper. The effect of surface tension is neglected. In the process of obtaining solutions for both the problems, regular perturbation analysis along with a Fourier transform technique is employed to derive the first-order corrections of some important physical quantities. Two types of bottom topography, such as concave and convex, are considered to derive the profiles of the interfaces. We observe that the profiles are oscillatory in nature, representing waves of variable amplitude with distinct wave numbers propagating downstream and with no wave upstream. The observations are presented in tabular and graphical forms.

Airblast spray in crossflow - Structure, trajectory and droplet sizing

This study reports results of an experimental investigation of airblast spray of water and ethanol in crossflow. Laser shadowgraphy and Particle/Droplet Imaging Analysis (PDIA) are used to derive spray trajectory and drop size information while Particle Tracking Velocimetry (PTV) is used to measure droplet velocities. A new phenomenon of spray bifurcation is observed for low Gas to Liquid Ratio (GLR) cases. The reasons for the spatial bifurcation can be attributed to a combination of reasons. These are (a) presence of large ligaments and droplets in the near-nozzle region for low GLRs (b) secondary breakup experienced by ligaments/droplets leading to formation of a large number of small droplets, and (c) the crossflow causing differential dispersion of the small and large droplets. A novel correlation for spray trajectory is proposed incorporating the momentum ratio and liquid surface tension. This correlation is shown to be effective in predicting the non-linear spray trajectory over a large range of conditions for not only water but ethanol and Jet-A also. It is observed that the larger droplets penetrate further into the crossflow, in the direction of injection. Thus, with increase in height of the measurement location from the injection plane, the droplet Sauter Mean Diameter (SMD) is found to increase. Moreover, as the droplets travel downstream in the crossflow direction, the droplet SMD is observed to decrease. The effect of drag is assessed by comparing velocity of different sizes of droplets at various locations. Smaller droplets are entrained into the crossflow at much lower elevations, whereas larger droplets tend to penetrate further into the crossflow. (C) 2015 Elsevier Ltd. All rights reserved.

Mechanical properties of CNT-Bisphenol E cyanate ester-based CFRP nanocomposite developed through VARTM process

This paper reports on the effect of multiwall carbon nanotubes (CNTs) without and with chemical functionalization on the mechanical properties of Bisphenol E cyanate ester resin (BECy) based carbon fibre reinforced plastic (CFRP) laminated composites. BECy with its low viscosity, low moisture uptake and superior mechanical properties is selected for its application in CFRP laminates through the cost-effective Vacuum Assisted Resin Transfer Moulding (VARTM) process. However, unlike CNT-epoxy-CFRP composites, processing and performance issues such as dispersion of CNTs, chemical bonding with resin, functionalization effects, effects on mechanical properties, etc. for BECy-CNT-CFRP composite system are not well reported. The objective of this study is to improve the mechanical properties of BECy resin with small additions of CNTs and functionalized CNTs in CFRP laminates. CNTs and fCNTs are infused into BECy using ultrasonication and standard mixing methods. Improvements in Young's modulus and strength in tension, compression, shear, flexure and interlaminar shear strength are analysed. It is observed that addition of 0.5wt% CNTs effected for maximum mechanical properties of the resin and 1wt% CNTs for the mechanical properties of CNT-CFRP nanocomposite. Further, improvements obtained with fCNTs are marginal. Dispersion behaviour and effect of CNTs/fCNTs in load transfer corroborated with SEM pictures are presented. The enhanced mechanical properties realized in VARTM processing of BECy-CFRP laminate indicate clear advantage of CNT based modification of the process.

Unified Continuum Damage Model for Matrix Cracking in Composite Rotor Blades

This paper deals with modeling of the first damage mode, matrix micro-cracking, in helicopter rotor/wind turbine blades and how this effects the overall cross-sectional stiffness. The helicopter/wind turbine rotor system operates in a highly dynamic and unsteady environment leading to severe vibratory loads present in the system. Repeated exposure to this loading condition can induce damage in the composite rotor blades. These rotor/turbine blades are generally made of fiber-reinforced laminated composites and exhibit various competing modes of damage such as matrix micro-cracking, delamination, and fiber breakage. There is a need to study the behavior of the composite rotor system under various key damage modes in composite materials for developing Structural Health Monitoring (SHM) system. Each blade is modeled as a beam based on geometrically non-linear 3-D elasticity theory. Each blade thus splits into 2-D analyzes of cross-sections and non-linear 1-D analyzes along the beam reference curves. Two different tools are used here for complete 3-D analysis: VABS for 2-D cross-sectional analysis and GEBT for 1-D beam analysis. The physically-based failure models for matrix in compression and tension loading are used in the present work. Matrix cracking is detected using two failure criterion: Matrix Failure in Compression and Matrix Failure in Tension which are based on the recovered field. A strain variable is set which drives the damage variable for matrix cracking and this damage variable is used to estimate the reduced cross-sectional stiffness. The matrix micro-cracking is performed in two different approaches: (i) Element-wise, and (ii) Node-wise. The procedure presented in this paper is implemented in VABS as matrix micro-cracking modeling module. Three examples are presented to investigate the matrix failure model which illustrate the effect of matrix cracking on cross-sectional stiffness by varying the applied cyclic

Interaction of a vortex ring with a single bubble: bubble and vorticity dynamics

The interaction of a single bubble with a single vortex ring in water has been studied experimentally. Measurements of both the bubble dynamics and vorticity dynamics have been done to help understand the two-way coupled problem. The circulation strength of the vortex ring (Gamma) has been systematically varied, while keeping the bubble diameter (D-b) constant, with the bubble volume to vortex core volume ratio (V-R) also kept fixed at about 0.1. The other important parameter in the problem is a Weber number based on the vortex ring strength. (We = 0.87 rho(Gamma/2 pi a)(2)/(sigma/D-b); a = vortex core radius, sigma = surface tension), which is varied over a large range, We = 3-406. The interaction between the bubble and ring for each of the We cases broadly falls into four stages. Stage I is before capture of the bubble by the ring where the bubble is drawn into the low-pressure vortex core, while in stage II the bubble is stretched in the azimuthal direction within the ring and gradually broken up into a number of smaller bubbles. Following this, in stage III the bubble break-up is complete and the resulting smaller bubbles slowly move around the core, and finally in stage IV the bubbles escape. Apart from the effect of the ring on the bubble, the bubble is also shown to significantly affect the vortex ring, especially at low We (We similar to 3). In these low-We cases, the convection speed drops significantly compared to the base case without a bubble, while the core appears to fragment with a resultant large decrease in enstrophy by about 50 %. In the higher-We cases (We > 100), there are some differences in convection speed and enstrophy, but the effects are relatively small. The most dramatic effects of the bubble on the ring are found for thicker core rings at low We (We similar to 3) with the vortex ring almost stopping after interacting with the bubble, and the core fragmenting into two parts. The present idealized experiments exhibit many phenomena also seen in bubbly turbulent flows such as reduction in enstrophy, suppression of structures, enhancement of energy at small scales and reduction in energy at large scales. These similarities suggest that results from the present experiments can be helpful in better understanding interactions of bubbles with eddies in turbulent flows.