The Reproducing Kernel DMS-FEM: 3D Shape Functions and Applications to Linear Solid Mechanics

We propose a family of 3D versions of a smooth finite element method (Sunilkumar and Roy 2010), wherein the globally smooth shape functions are derivable through the condition of polynomial reproduction with the tetrahedral B-splines (DMS-splines) or tensor-product forms of triangular B-splines and ID NURBS bases acting as the kernel functions. While the domain decomposition is accomplished through tetrahedral or triangular prism elements, an additional requirement here is an appropriate generation of knotclouds around the element vertices or corners. The possibility of sensitive dependence of numerical solutions to the placements of knotclouds is largely arrested by enforcing the condition of polynomial reproduction whilst deriving the shape functions. Nevertheless, given the higher complexity in forming the knotclouds for tetrahedral elements especially when higher demand is placed on the order of continuity of the shape functions across inter-element boundaries, we presently emphasize an exploration of the triangular prism based formulation in the context of several benchmark problems of interest in linear solid mechanics. In the absence of a more rigorous study on the convergence analyses, the numerical exercise, reported herein, helps establish the method as one of remarkable accuracy and robust performance against numerical ill-conditioning (such as locking of different kinds) vis-a-vis the conventional FEM.

Squeezed states, metaplectic group, and operator möbius transformations

We present an analysis, based on the metaplectic group Mp(2), of the recently introduced single-mode inverse creation and annihilation operators and of the associated eigenstates of different two-photon annihilation operators. We motivate and obtain a quantum operator form of the classical Mobius or fractional linear transformation. The subtle relation to the two unitary irreducible representations of Mp(2) is brought out. For problems involving inverse operators the usefulness of the Bargmann analytic function representation of quantum mechanics is demonstrated. Squeezing, bunching, and photon-number distributions of the four families of states that arise in this context are studied both analytically and numerically

Analysis of a long thick orthotropic circular cylindrical shell panel

An exact three-dimensional elasticity solution has been obtained for an infinitely long, thick transversely isotropic circular cylindrical shell panel, simply supported along the longitudinal edges and subjected to a radial patch load. Using a set of three displacement functions, the boundary value problem is reduced to Bessel's differential equation. Numerical results are presented for different thickness to mean radius ratios and semicentral angles of the shell panel. Classical and first-order shear deformation orthotropic shell theories have been examined in comparison with the present elasticity solution.

Gas-liquid nucleation at large metastability: unusual features and a new formalism

Nucleation at large metastability is still largely an unsolved problem, even though it is a problem of tremendous current interest, with wide-ranging practical value, from atmospheric research to materials science. It is now well accepted that the classical nucleation theory (CNT) fails to provide a qualitative picture and gives incorrect quantitative values for such quantities as activation-free energy barrier and supersaturation dependence of nucleation rate, especially at large metastability. In this paper, we present an alternative formalism to treat nucleation at large supersaturation by introducing an extended set of order parameters in terms of the kth largest liquid-like clusters, where k = 1 is the largest cluster in the system, k = 2 is the second largest cluster and so on. At low supersaturation, the size of the largest liquid-like cluster acts as a suitable order parameter. At large supersaturation, the free energy barrier for the largest liquid-like cluster disappears. We identify this supersaturation as the one at the onset of kinetic spinodal. The kinetic spinodal is system-size-dependent. Beyond kinetic spinodal many clusters grow simultaneously and competitively and hence the nucleation and growth become collective. In order to describe collective growth, we need to consider the full set of order parameters. We derive an analytic expression for the free energy of formation of the kth largest cluster. The expression predicts that, at large metastability (beyond kinetic spinodal), the barrier of growth for several largest liquid-like clusters disappears, and all these clusters grow simultaneously. The approach to the critical size occurs by barrierless diffusion in the cluster size space. The expression for the rate of barrier crossing predicts weaker supersaturation dependence than what is predicted by CNT at large metastability. Such a crossover behavior has indeed been observed in recent experiments (but eluded an explanation till now). In order to understand the large numerical discrepancy between simulation predictions and experimental results, we carried out a study of the dependence on the range of intermolecular interactions of both the surface tension of an equilibrium planar gas-liquid interface and the free energy barrier of nucleation. Both are found to depend significantly on the range of interaction for the Lennard-Jones potential, both in two and three dimensions. The value of surface tension and also the free energy difference between the gas and the liquid phase increase significantly and converge only when the range of interaction is extended beyond 6-7 molecular diameters. We find, with the full range of interaction potential, that the surface tension shows only a weak dependence on supersaturation, so the reason for the breakdown of CNT (with simulated values of surface tension and free energy gap) cannot be attributed to the supersaturation dependence of surface tension. This remains an unsettled issue at present because of the use of the value of surface tension obtained at coexistence.

Multifractal analysis, a method to investigate the morphology of materials

Using the multifractal formalism, we discuss the results obtained to characterized the morphology of polymer alloys and granular discontinuous metallic thin films. In the first case we have found a correlation between the multifractality and the mechanical properties of the alloys. In the second case, we have found that it is possible to measure the differences between the morphology of thin films induced by a growth process on a subtrate and that of percolation clusters of the classical theory of percolation.

Prediction of nonlocal scaling parameter for armchair and zigzag single-walled carbon nanotubes based on molecular structural mechanics, nonlocal elasticity and wave propagation

Atomic vibration in the Carbon Nanotubes (CNTs) gives rise to non-local interactions. In this paper, an expression for the non-local scaling parameter is derived as a function of the geometric and electronic properties of the rolled graphene sheet in single-walled CNTs. A self-consistent method is developed for the linearization of the problem of ultrasonic wave propagation in CNTs. We show that (i) the general three-dimensional elastic problem leads to a single non-local scaling parameter (e(0)), (ii) e(0) is almost constant irrespective of chirality of CNT in the case of longitudinal wave propagation, (iii) e(0) is a linear function of diameter of CNT for the case of torsional mode of wave propagation, (iv) e(0) in the case of coupled longitudinal-torsional modes of wave propagation, is a function which exponentially converges to that of axial mode at large diameters and to torsional mode at smaller diameters. These results are valid in the long-wavelength limit. (C) 2011 Elsevier Ltd. All rights reserved.

Critical buckling temperature of single-walled carbon nanotubes embedded in a one-parameter elastic medium based on nonlocal continuum mechanics

In this paper, the critical budding temperature of single-walled carbon nanotubes (SWCNTs), which are embedded in one-parameter elastic medium (Winkler foundation) is estimated under the umbrella of continuum mechanics theory. Nonlocal continuum theory is incorporated into Timoshenko beam model and the governing differential equations of motion are derived. An explicit expression for the non-dimensional critical buckling temperature is also derived in this work. The effect of the nonlocal small scale coefficient, the Winkler foundation parameter and the ratio of the length to the diameter on the critical buckling temperature is investigated in detail. It can be observed that the effects of nonlocal small scale parameter and the Winkler foundation parameter are significant and should be considered for thermal analysis of SWCNTs. The results presented in this paper can provide useful guidance for the study and design of the next generation of nanodevices that make use of the thermal buckling properties of embedded single-walled carbon nanotubes. (C) 2011 Elsevier B.V. All rights reserved.

Characterization of estuarine marine clays for coastal reclamation in Pusan, Korea

Large scale reclamation works in coastal areas of the Nakdong River plain are at various stages of progress, since early 1990's on in-situ soft marine clay deposits. These deposits are of the order of 30 to 40 m thick. A realistic rapid characterization of soft ground would ensure success of any reclamation work in this area. In order to cope with the work carried out with different agencies, it is desirable to evolve a systematic methodology. In this study, engineering properties of clays at three coastal areas, Gadukdo, Noksan and Shinho, have been generated. The analysis of data has been done within the framework of classical developments in soil mechanics. Analysis has also been made by making use of the recent developments in dealing with soft clays. The dominant factors, namely, stress, time, and environment influencing the response of clay to loading are identified.

Freezing of classical fluids in a quenched random potential

The phase diagram of a hard-sphere fluid in the presence of a random pinning potential is studied analytically and numerically. In the analytic work, replicas are introduced for averaging over the quenched disorder, and the hypernetted chain approximation is used to calculate density correlations in the replicated liquid. The freezing transition of the liquid into a nearly crystalline state is studied using a density-functional approach, and the liquid to glass transition is studied using a phenomenological replica symmetry breaking approach. In the numerical work, local minima of a discretized version of the Ramakrishnan-Yussouff free-energy functional are located and the phase diagram in the density-disorder plane is obtained from an analysis of the relative stability of these minima. Both approaches lead to similar results for the phase diagram. The first-order liquid to crystalline solid transition is found to change to a continuous liquid to glass transition as the strength of the disorder is increased above a threshold value.

Simulation of the infrared spectra of thioacetamide by the extended molecular mechanics method

The infrared spectrum of the matrix-isolated species of thioacetamide has been simulated using the extended molecular mechanics method. The equilibrium structure, vibrational frequencies, dipole moment and infrared absorption intensities of thioacetamide have been calculated in good agreement with the experiment. The vibrational frequencies and infrared absorption intensities for the isotopic molecules (CH2CSNH2)-C-13, (CH3CSNH2)-N-15 and (CH2CSND2)-C-13 have also been calculated consistent with the experiment. The infrared spectra of the matrix isolated species of N- and C- deuterated isotopomers of thioacetamide, CH3CSND2 and CD3CSNH2 have also been simulated in satisfactory agreement with the experimental spectra.

Simulation of the infrared spectra of N-methylformamide and N-methylacetamide by the extended molecular mechanics method

The infrared spectra of the matrix isolated species of N-methylformamide (NMF) and N-methylacetamide (NMA) and their N-deuterated molecules have been simulated by the extended molecular mechanics method using an empirical force field which includes charges and charge fluxes as coulombic potential parameters. The structural parameters and dipole. moments of NMF and NMA have. also been computed in satisfactory agreement with the experiment. Good agreement between experimental and calculated vibrational frequencies and infrared absorption band intensities for NMF and NMA and their deuterated molecules has been obtained. The vibrational assignments of NMF and NMA are-discussed taking also into account the infrared absorption intensities.

Axial wave propagation in coupled nanorod system with nonlocal small scale effects

This article deals with the axial wave propagation properties of a coupled nanorod system with consideration of small scale effects. The nonlocal elasticity theory has been incorporated into classical rod/bar model to capture unique features of the coupled nanorods under the umbrella of continuum mechanics theory. Nonlocal rod model is developed for coupled nanorods. The strong effect of the nonlocal scale has been obtained which leads to substantially different wave behavior of nanorods from those of macroscopic rods. Explicit expressions are derived for wavenumber, cut-off frequency and escape frequency of nanorods. The analysis shows that the wave characteristics of nanorods are highly over estimated by the classical rod model, which ignores the effect of small-length scale. The studies also shows that the nonlocal scale parameter introduces certain band gap region in axial or longitudinal wave mode, where no wave propagation occurs. This is manifested in the spectrum cures as the region, where the wavenumber tends to infinite or wave speed tends to zero. The effect of the coupled spring stiffness is also capture in the present analysis. It has been also shown that the cut-off frequency increases as the stiffness of the coupled spring increases and also the coupled spring stiffness has no effect on escape frequency of the axial wave mode in the nanorod. This cut-off frequency is also independent of the nonlocal small scale parameter. The present study may bring in helpful insights while investigating multiple-nanorod-system-models for future nano-optomechanical systems applications. The results can also provide useful guidance for the study and design of the next generation of nanodevices that make use of the wave propagation properties of coupled single-walled carbon nanotubes or coupled nanorods. (C) 2011 Elsevier Ltd. All rights reserved.