Buckling analysis of multi-walled carbon nanotubes: a continuum model accounting for van der Waals interaction

Explicit formulas are derived for the van der Waals (vdW) interaction between any two layers of a multi-walled carbon nanotube (CNT). Based on the derived formulas, an efficient algorithm is established for the buckling analysis of multi-walled CNTs, in which individual tubes are modeled as a continuum cylindrical shell. The explicit expressions are also derived for the buckling of double-walled CNTs. In previous studies by Ru (J. Appl. Phys. 87 (2000b) 7227) and Wang et al. (Int. J. Solids Struct. 40 (2003) 3893), only the vdW interaction between adjacent two layers was considered and the vdW interaction between the other two layers was neglected. The vdW interaction coefficient was treated as a constant that was not dependent on the radii of the tubes. However, the formulas derived herein reveal that the vdW interaction coefficients are dependent on the change of interlayer spacing and the radii of the tubes. With the increase of radii, the coefficients approach constants, and the constants between two adjacent layers are about 10% higher than those reported by Wang et al. (Int. J. Solids. Struct. 40 (2003) 3893). In addition, the numerical results show that the vdW interaction will lead to a higher critical buckling load in multi-walled CNTs. The effect of the tube radius on the critical buckling load of a multi-walled CNT is also examined.     

Curvature effects on axially compressed buckling of a small-diameter double-walled carbon nanotube

The curvature effects of interlayer van der Waals (vdW) forces on axially compressed buckling of a double-walled carbon nanotube (DWNT) of diameter down to 0.7 nm are studied. Unlike most existing models which assume that the interlayer vdW pressure at a point between the inner and outer tubes depends merely on the change of the interlayer spacing at that point, the present model considers the dependence of the interlayer vdW pressure on the change of the curvatures of the inner and outer tubes at that point. A simple expression is derived for the curvature-dependence of the interlayer vdW pressure in which the curvature coefficient is determined. Based on this model, an explicit formula is obtained for the axial buckling strain. It is shown that neglecting the curvature effect alone leads to an under-estimate of the critical buckling strain with a relative error up to −7%, while taking the average radius of two tubes as the representative radius and the curvature effect leads to an over-estimate of the critical buckling strain with a relative error up to 20% when the inner radius downs to 0.35 nm. Therefore, the curvature effects play a significant role in axially compressed buckling problems only for DWNTs of very small radii. In addition, our results show that the effect of the vdW interaction pressure prior to buckling of DWNTs under pure axial stress is small enough and can be negligible whether the vdW interaction curvature effects are neglected or not.     

Buckling of embedded multi-walled carbon nanotubes under combined torsion and axial loading

This paper describes an investigation into elastic buckling of an embedded multi-walled carbon nanotube under combined torsion and axial loading, which takes account of the radial constraint from the surrounding elastic medium and van der Waals force between two adjacent tube walls. Depending on the ratio of radius to thickness, the multi-walled carbon nanotubes discussed here are classified as thin, thick, and nearly solid. Critical buckling load with the corresponding mode is obtained for multi-walled carbon nanotubes under combined torsion and axial loading, with various values of the radius to thickness ratio and surrounded with different elastic media. The study indicates that the buckling mode (m, n) of an embedded multi-walled carbon nanotube under combined torsion and axial loading is unique and it is different from that with axial compression only. New features for the buckling of an embedded multi-walled carbon nanotube under combined torsion and axial loading and the meaningful numerical results are useful in the design of nanodrive device, nanotorsional oscillator and rotational actuators, where multi-walled carbon nanotubes act as basic elements.     

Analytical and numerical techniques to predict carbon nanotubes properties

In this paper, two different approaches for modeling the behaviour of carbon nanotubes are presented. The first method models carbon nanotubes as an inhomogeneous cylindrical network shell using the asymptotic homogenization method. Explicit formulae are derived representing Young’s and shear moduli of single-walled nanotubes in terms of pertinent material and geometric parameters. As an example, assuming certain values for these parameters, the Young’s modulus was found to be 1.71 TPa, while the shear modulus was 0.32 TPa. The second method is based on finite element models. The inter-atomic interactions due to covalent and non-covalent bonds are replaced by beam and spring elements, respectively, in the structural model. Correlations between classical molecular mechanics and structural mechanics are used to effectively model the physics governing the nanotubes. Finite element models are developed for single-, double- and multi-walled carbon nanotubes. The deformations from the finite element simulations are subsequently used to predict the elastic and shear moduli of the nanotubes. The variation of mechanical properties with tube diameter is investigated for both zig-zag and armchair configurations. Furthermore, the dependence of mechanical properties on the number of nanotubules in multi-walled structures is also examined. Based on the finite element model, the value for the elastic modulus varied from 0.9 to 1.05 TPa for single and 1.32 to 1.58 TPa for double/multi-walled nanotubes. The shear modulus was found to vary from 0.14 to 0.47 TPa for single-walled nanotubes and 0.37 to 0.62 for double/multi-walled nanotubes.     

Effects of layer number and initial pressure on continuum-based buckling analysis of multi-walled carbon nanotubes accounting for van der Waals interaction

The structural instability of multi-walled carbon nanotubes (MWCNTs) has captured extensive attention due to the unique characteristic of extremely thin hollow cylinder structure. The previous studies usually focus on the buckling behavior without considering the effects of the wall number and initial pressure. In this paper, the axial buckling behavior of MWCNTs with the length-to-outermost radius ratio less than 20 is investigated within the framework of the Donnell shell theory. The governing equations for the infinitesimal buckling of MWCNTs are established, accounting for the van der Waals (vdW) interaction between layers. The effects of the wall number, initial pressure prior to buckling, and aspect ratio on the critical buckling mode, buckling load, and buckling strain are discussed, respectively. Specially, the four-walled and twenty-walled CNTs are studied in detail, indicating the fact that the buckling instability may occur in other layers besides the outermost layer. The obtained results extend the buckling analysis of the continuum-based model, and provide theoretical support for the application of CNTs.     

Experimental study on the compressive behaviors of CFRP tube-encased concrete columns

Nine square concrete columns including 6 CFRP/ECCs and 3 concrete columns are prepared, which have cross-section of 200 mm × 200 mm and height of 600 mm. The CFRP tubes with fibers oriented at hoop direction were manufactured to have 3 or 5 layers of CFRP with 10 mm, 20 mm, or 40 mm rounding corner radii at vertical edges. A 100 mm overlap in the direction of fibers was provided to ensure proper bond. Uniaxial compression tests were conducted to investigate the compressive behavior. It is evident that the CFRP tube confinement can improve the behavior of concrete core, in terms of axial compressive strength or axial deformability. Test results show that the stress-strain behavior of CFRP/ECCs vary with different confinement parameters, such as the number of confinement layers and the rounding corner radius.     

Theoretical stress analysis of intersecting cylindrical shells subjected to external loads transmitted through branch pipes

A thin shell theoretical solution of two normally intersecting cylindrical shells subjected to thrust-out force and three kinds of moments transmitted through branch pipes is presented in this paper. The solutions of modified Morley equation, which can be applicable up to ρ₀ = d/D ⩽ 0.8 and λ = d/(DT)^1/2 ⩽ 8 and the order of accuracy is raised to O(T/D), for the four loading cases are given. The accurate continuity conditions of generalized forces and displacements at the intersecting curve of two cylindrical shells for the four loading cases and the condition of the uniqueness of displacements are derived in this paper. The presented results are verified by experimental and numerical results successfully. They are in agreement with WRC Bulletin 297 when d/D is small.     

Carbon nanotube glycol nanofluids: Photo-thermal properties,thermal conductivities and rheological behavior

The efficiency and effectiveness of solar energy capture and storage are to a large extent functions of the heat transfer and storage capacity of the medium used. This paper investigates the potential of using carbon nanotube (CNT)-glycol nanosuspension as such a medium, prepared by freeze drying-ultrasonic dispersing after oxidation treatment with HNO₃. The influences of the mass fraction of CNTs glycol nanofluids and temperatures on photo-thermal properties, thermal conductivities and rheological behavior were investigated. The results show that CNTs with oxidation treatment exhibited good dispersing performance. Strong optical absorption of the CNTs glycol nanofluids was detected in the range of 200–2500 nm. At room temperature, 18% enhancement was found in the photo-thermal conversion efficiency of the 0.5% mass fraction CNTs glycol nanofluids in comparison to the basic fluids, without significant increase in viscosity. At 55 °C, CNTs glycol nanofluids with 4.0% mass fraction exhibited much lower viscosity and 25.4% higher thermal conductivity in comparison to that of pure glycol at room temperature.     

Developing structure of two-phase flow in a large diameter pipe at low liquid flow rate

In order to develop the interfacial area transport equation for the interfacial transfer terms in the two-fluid model, accurate data sets on axial development of local parameters such as void fraction, interfacial area concentration, interfacial gas velocity and Sauter mean diameter are indispensable to verify the modeled source and sink terms in the interfacial area transport equation. From this point of view, local measurements of both group 1 spherical/distorted bubbles and group 2 cap/slug bubbles in vertical upward air–water two-phase flow in a large diameter pipe with 200 mm in inner diameter and 26 m in height were performed at three axial locations of z/D = 41.5, 82.8 and 113 as well as 11 radial locations from r/R = 0–0.95 by using four-sensor probe method. Here, z, r, D and R are the axial distance from the inlet, radial distance from the pipe center, pipe diameter and pipe radius, respectively. The liquid flow rate and the void fraction ranged from 0.0505 m/s to 0.312 m/s and from 1.98% to 32.6%, respectively in the present experiment. The flow condition covered extensive region of bubbly flow, cap turbulent flow as well as their transition. The extensive analysis on the radial profiles of local flow parameters and their axial developments demonstrate the development of interfacial structures along the flow direction due to the bubble coalescence and breakup and the gas expansion. The significant decrease in void faction and interfacial area concentration and the increase in Sauter mean diameter and interfacial velocity were observed when the gradual flow regime transition occurred. Finally, the net change in the interfacial area concentration due to the bubble coalescence and breakup was quantitatively investigated in the present paper to reflect the true transfer mechanisms in observed two-phase flows.     

Adaptive and anisotropic mesh strategy for thin shell problems. Case of inhibited parabolic shells

The aim of this paper is to show the reliability of an adaptive and anisotropic mesh procedure for thin shell problems. We consider singular perturbation problems only for parabolic shells whose behavior is described by the Koiter model. The corresponding system of equations, which depends on the relative thickness ε of the shell, is elliptic except at the limit for ε = 0 where it is parabolic. In a first part of this paper, we study theoretically the phenomena of internal layers appearing during the singular perturbation process, when the loading is somewhat singular. These layers have very different structures either they are along or across the asymptotic lines of the middle surface of the shell. In a second part, numerical computations are performed using a finite element software coupled with an adaptive anisotropic mesh generator. This technique enables to approach accurately the singularities and the layers predicted by the theory especially for very small values of the thickness. The efficiency of such a procedure in comparison with uniform meshes is put in a prominent position.     

Experimental study of single bubble motion in a liquid metal column exposed to a DC magnetic field

The motion of single Argon bubbles rising in the eutectic alloy GaInSn under the influence of a DC longitudinal magnetic field (parallel to the direction of bubble motion) was examined. The magnetic field strength was varied up to 0.3 T corresponding to a magnetic interaction parameter N (which measures the ratio of electromagnetic forces to inertial forces) slightly greater than 1. The liquid metal was at rest in a cylindrical container. Bubble and liquid velocities were measured using ultrasound Doppler velocimetry (UDV). The measured bubble terminal velocity showed oscillations indicating a zigzag movement of ellipsoidal bubbles. For small bubbles (d_e ⩽ 4.6 mm) an increase of the drag coefficient with increasing magnetic interaction parameter N was observed, whereas for larger bubbles (d_e ⩾ 5.4 mm) the application of the magnetic field reduces the drag coefficient. The measurements revealed a distinct electromagnetic damping of the bubble induced liquid velocity leading to more rectilinear bubble trajectories when the magnetic field is applied. Moreover, significant modifications of the bubble wake structure were observed. Raising of the magnetic field strength caused an enlargement of the eddies in the wake. The Strouhal number decreases with increasing magnetic interaction parameter N.     

Experimental investigation of a developing two-phase bubbly flow in horizontal pipe

Experimental results for various water and air superficial velocities in developing adiabatic horizontal two-phase pipe flow are presented. Flow pattern maps derived from videos exhibit a new boundary line in intermittent regime. This transition from water dominant to water–gas coordinated regimes corresponds to a new transition criterion C_T = 2, derived from a generalized representation with the dimensionless coordinates of Taitel and Dukler.Velocity, turbulent kinetic energy and dissipation rate, void fraction and bubble size radial profiles measured at 40 pipe diameters for J_L = 4.42 m/s by hot film velocimetry and optical probes confirm this transition: the gas influence is not continuous but strongly increases beyond J_G = 0.06 m/s. The maximum dissipation rate, derived from spectra, is increased in two-phase flow by a factor 5 with respect to the single phase case.The axial evolution of the bubble intercept length histograms also reveal the flow organization in horizontal layers, driven by buoyancy effects. Bubble coalescence is attested by a maximum bubble intercept evolving from 2.5 to 4.5 mm along the pipe. Turbulence generated by the bubbles is also manifest by the 4-fold increase of the maximum turbulent dissipation rate along the pipe.     

VIBRATION OF FLUID-FILLED MULTI-WALLED CARBON NANOTUBES SEEN VIA NONLOCAL ELASTICITY THEORY

Vibration characteristics of fluid-filled multi-walled carbon nanotubes axe studied by using nonlocal elastic Fliigge shell model. Vibration governing equations of an N-layer carbon nanotube are formulated by considering the scale effect. In the numerical simulations, the effects of different theories, small-scale, variation of wavenumber, the innermost radius and length of double- walled and triple-walled carbon nanotubes are considered. Vibrational frequencies decrease with an increase of scale coefficient, the innermost radius, length of nanotube and effects of wall number are negligible. The results show that the cut-off frequencies can be influenced by the wall number of nanotubes.     

CFD based mini- vs. micro-system delineation in elongated bubble flow regime

Delineation of mini- and micro-scale channels with respect to two-phase flow has been the subject of many research papers. There is no consensus on when the small channel can be characterized as a mini-channel or micro-channel. The idea proposed by this paper is to use the normalized bubble nose radius, liquid film thickness top over bottom ratio, and bubble shape contour, which are found under normal gravity conditions in slug flow through a horizontal adiabatic channel, as the delineation criteria. The input parameters are bubble nose radius and bubble nose velocity as the characteristic length scale and characteristic velocity scale respectively. 3D numerical simulation with ANSYS FLUENT was used to obtain the necessary data. Following CFD practice, a mesh independence study and a numerical model validation against published experimental data were both conducted. Analysis of the numerical simulation results showed that channels with D ⩽ 100 μm can be characterized as a micro-system, while channels with D ⩾ 400 μm belong to mini-systems. The region 200 μm ⩽ D ⩽ 300 μm represents a transition from the micro-scale to mini-scale.     

LES of turbulent flow in a concentric annulus with rotating outer wall

Fully-developed turbulent flow in a concentric annulus, r₁/r₂ = 0.5, Re_h = 12,500, with the outer wall rotating at a range of rotation rates N = U_θ,wall/U_b from 0.5 up to 4 is studied by large-eddy simulations. The focus is on the effects of moderate to very high rotation rates on the mean flow, turbulence statistics and eddy structure. For N up to ∼2, an increase in the rotation rate dampens progressively the turbulence near the rotating outer wall, while affecting only mildly the inner-wall region. At higher rotation rates this trend is reversed: for N = 2.8 close to the inner wall turbulence is dramatically reduced while the outer wall region remains turbulent with discernible helical vortices as the dominant turbulent structure. The turbulence parameters and eddy structures differ significantly for N = 2 and 2.8. This switch is attributed to the centrifuged turbulence (generated near the inner wall) prevailing over the axial inertial force as well as over the counteracting laminarizing effects of the rotating outer wall. At still higher rotation, N = 4, the flow gets laminarized but with distinct spiralling vortices akin to the Taylor–Couette rolls found between the two counter-rotating cylinders without axial flow, which is the limiting case when N approaches to infinity. The ratio of the centrifugal to axial inertial forces, Ta/Re² ∝ N² (where Ta is the Taylor number) is considered as a possible criterion for defining the conditions for the above regime change.     

Laminar,transitional and turbulent annular flow of drag-reducing polymer solutions

Mean and rms axial velocity-profile data obtained using laser Doppler anemometry are presented together with pressure-drop data for the flow through a concentric annulus (radius ratio κ = 0.506) of a Newtonian (a glycerine–water mixture) and non-Newtonian fluids—a semi-rigid shear-thinning polymer (a xanthan gum) and a polymer known to exhibit a yield stress (carbopol). A wider range of Reynolds numbers for the transitional flow regime is observed for the more shear-thinning fluids. In marked contrast to the Newtonian fluid, the higher shear stress on the inner wall compared to the outer wall does not lead to earlier transition for the non-Newtonian fluids where more turbulent activity is observed in the outer wall region. The mean axial velocity profiles show a slight shift (～5%) of the location of the maximum velocity towards the outer pipe wall within the transitional regime only for the Newtonian fluid.     

Neutrality in the case of N-phase elliptical inclusions with internal uniform hydrostatic stresses

A novel method is proposed to design neutral N-phase (N ? 3) elliptical inclusions with internal uniform hydrostatic stresses. We focus on the study of the internal and external stress states of an N-phase elliptical inclusion which is bonded to an infinite matrix through (N ? 2) interphase layers. The interfaces of the N-phase elliptical inclusion are (N ? 1) confocal ellipses. The design of the resulting overall composite material consists of four stages: (i) an inner perfectly bonded interphase/inclusion interface which is necessary to make the internal uniform stress state hydrostatic; (ii) outer imperfect interphase layers properly designed to make the coated inclusion harmonic (i.e., the uniform mean stress of the original field within the matrix is unperturbed); (iii) the aspect ratio of the elliptic inclusion uniquely chosen for a given material and thickness parameters to make the resulting coated inclusion neutral (i.e., the prescribed uniform stress field in the matrix remains undisturbed); and finally (iv) the derivation of a simple condition relating the remote uniform stresses and the thickness parameters of the (N ? 2) interphase layers for given material parameters which lead to internal uniform hydrostatic stresses. We note that another interesting feature of the present results is that the mean stress is found to be constant within each interphase layer, and the hoop stress in the innermost interphase layer is uniform along the entire interphase/inclusion interface.     

Continuum interpretation of virial stress in molecular simulations

The equivalence of the virial stress and Cauchy stress is reviewed using both theoretical arguments and numerical simulations. Using thermoelasticity problems as examples, we numerically demonstrate that virial stress is equivalent to the continuum Cauchy stress. Neglecting the velocity terms in the definition of virial stress as many authors have recently suggested, can cause significant errors in interpreting MD simulation results at elevated temperatures (T > 0 K).     

Nonaxisymmetric interaction of a spherical radiator in a fluid-filled permeable borehole

The Biot theory of poroelasticity along with the proper cylindrical/spherical wave-field transformations are used to investigate general (nonaxisymmetric) harmonic radiation from a spherical surface vibrating at the center of a fluid-filled circular cylindrical cavity embedded within a fluid-saturated porous elastic formation. This configuration, which is a realistic idealization of an acoustic logging tool suspended in a fluid-filled borehole, is of practical importance with a multitude of possible applications in seismic engineering and geophysics. The analytical results are illustrated with numerical examples in which the spherical source suspended at the center of a water-filled borehole embedded within water-saturated soils of distinct frame properties (i.e., soft or stiff soils), is excited in vibrational modes of various orders. The basic acoustic and elastic field quantities such as the resistive/reactive components of the modal acoustic radiation impedance load as well as the radial displacement and stress components induced within the surrounding formation for a pulsating (n = 0), an oscillating (n = 1), and a quadrupole-like (n = 2) spherical source are evaluated and discussed for representative values of the parameters characterizing the system. Special attention is paid to the effects of source excitation frequency, size, surface velocity profile, and internal impedance as well as soil type on the modal impedance values and the displacement/stress amplitudes. Limiting cases are considered and fair agreements with well-known solutions are obtained.