Nonlinear Dynamics of a Rotating Flexible Link

This paper deals with the study of the nonlinear dynamics of a rotating flexible link modeled as a one dimensional beam, undergoing large deformation and with geometric nonlinearities. The partial differential equation of motion is discretized using a finite element approach to yield four nonlinear, nonautonomous and coupled ordinary differential equations (ODEs). The equations are nondimensionalized using two characteristic velocities-the speed of sound in the material and a velocity associated with the transverse bending vibration of the beam. The method of multiple scales is used to perform a detailed study of the system. A set of four autonomous equations of the first-order are derived considering primary resonances of the external excitation and one-to-one internal resonances between the natural frequencies of the equations. Numerical simulations show that for certain ranges of values of these characteristic velocities, the slow flow equations can exhibit chaotic motions. The numerical simulations and the results are related to a rotating wind turbine blade and the approach can be used for the study of the nonlinear dynamics of a single link flexible manipulator.

Correlation equations for average deposition rate coefficients of nanoparticles in a cylindrical pore

Nanoparticle deposition behavior observed at the Darcy scale represents an average of the processes occurring at the pore scale. Hence, the effect of various pore-scale parameters on nanoparticle deposition can be understood by studying nanoparticle transport at pore scale and upscaling the results to the Darcy scale. In this work, correlation equations for the deposition rate coefficients of nanoparticles in a cylindrical pore are developed as a function of nine pore-scale parameters: the pore radius, nanoparticle radius, mean flow velocity, solution ionic strength, viscosity, temperature, solution dielectric constant, and nanoparticle and collector surface potentials. Based on dominant processes, the pore space is divided into three different regions, namely, bulk, diffusion, and potential regions. Advection-diffusion equations for nanoparticle transport are prescribed for the bulk and diffusion regions, while the interaction between the diffusion and potential regions is included as a boundary condition. This interaction is modeled as a first-order reversible kinetic adsorption. The expressions for the mass transfer rate coefficients between the diffusion and the potential regions are derived in terms of the interaction energy profile. Among other effects, we account for nanoparticle-collector interaction forces on nanoparticle deposition. The resulting equations are solved numerically for a range of values of pore-scale parameters. The nanoparticle concentration profile obtained for the cylindrical pore is averaged over a moving averaging volume within the pore in order to get the 1-D concentration field. The latter is fitted to the 1-D advection-dispersion equation with an equilibrium or kinetic adsorption model to determine the values of the average deposition rate coefficients. In this study, pore-scale simulations are performed for three values of Peclet number, Pe = 0.05, 5, and 50. We find that under unfavorable conditions, the nanoparticle deposition at pore scale is best described by an equilibrium model at low Peclet numbers (Pe = 0.05) and by a kinetic model at high Peclet numbers (Pe = 50). But, at an intermediate Pe (e.g., near Pe = 5), both equilibrium and kinetic models fit the 1-D concentration field. Correlation equations for the pore-averaged nanoparticle deposition rate coefficients under unfavorable conditions are derived by performing a multiple-linear regression analysis between the estimated deposition rate coefficients for a single pore and various pore-scale parameters. The correlation equations, which follow a power law relation with nine pore-scale parameters, are found to be consistent with the column-scale and pore-scale experimental results, and qualitatively agree with the colloid filtration theory. These equations can be incorporated into pore network models to study the effect of pore-scale parameters on nanoparticle deposition at larger length scales such as Darcy scale.

Physical insights into salicylic acid release from poly(anhydrides)

Salicylic acid (SA) based biodegradable polyanhydrides (PAHs) are of great interest for drug delivery in a variety of diseases and disorders owing to the multi-utility of SA. There is a need for the design of SA-based PAHs for tunable drug release, optimized for the treatment of different diseases. In this study, we devised a simple strategy for tuning the release properties and erosion kinetics of a family of PAHs. PAHs incorporating SA were derived from related aliphatic diacids, varying only in the chain length, and prepared by simple melt condensation polymerization. Upon hydrolysis induced erosion, the polymer degrades into cytocompatible products, including the incorporated bioactive SA and diacid. The degradation follows first order kinetics with the rate constant varying by nearly 25 times between the PAH obtained with adipic acid and that with dodecanedioic acid. The release profiles have been tailored from 100% to 50% SA release in 7 days across the different PAHs. The release rate constants of these semi-crystalline, surface eroding PAHs decreased almost linearly with an increase in the diacid chain length, and varied by nearly 40 times between adipic acid and dodecanedioic acid PAH. The degradation products with SA concentration in the range of 30-350 ppm were used to assess cytocompatibility and showed no cytotoxicity to HeLa cells. This particular strategy is expected to (a) enable synthesis of application specific PAHs with tunable erosion and release profiles; (b) encompass a large number of drugs that may be incorporated into the PAH matrix. Such a strategy can potentially be extended to the controlled release of other drugs that may be incorporated into the PAH backbone and has important implications for the rational design of drug eluting bioactive polymers.

Structure and Dynamics of Octamethylcyclotetrasiloxane Confined between Mica Surfaces

Using a molecular model for octamethylcydotetrasiloxane (OMCTS), molecular dynamics simulations are carried out to probe the phase state of OMCTS confined between two mica surfaces in equilibrium With a reservoir. Molecular dynamics simulations are carried out for elevations ranging from 5 to 35 K above the melting point for the OMCTS model used in this study. The Helmholtz free energy is, computed for a specific confinement using the :two-phase thermodynamic (2PT) method. Analysis of the in-plane pair correlation functions did not reveal signatures of freezing even under an extreme confinement of two layers. OMCTS is found to orient with a wide distribution of orientations with respect to the mica surface, with a distinct preference for the surface parallel configuration in the contact layers. The self-intermediate scattering function is found to decay with increasing relaxation times as the surface separation is decreased, and the two-step relaxation in the scattering function, a signature of glassy dynamics, distinctly evolves as the temperature is lowered. However, even at 5 K above the melting point, we did not observe a freezing transition and the self-intermediate scattering functions relax within 200 ps for the seven-layered confined system. The self diffusivity and relaxation times obtained from the Kohlrausch-Williams-Watts stretched exponential fits to the late alpha-relaxation exhibit power law scalings with the packing fraction as predicted by mode coupling theory. A distinct discontinuity in the Helmholtz free energy, potential energy, and a sharp change in the local bond order parameter, Q(4), was observed at 230 K for a five-layered system upon cooling, indicative of a first-order transition. A freezing point depression of about 30 K was observed for this five-layered confined system, and at the lower temperatures, contact layers were found to be disordered with long-range order present only in the inner layers. These dynamical signatures indicate that confined OMCTS undergoes a slowdown akin to a fluid approaching a glass transition upon increasing confinement, and freezing under confinement would require substantial subcooling below the bulk melting point of OMCTS.

A peridynamic theory for linear elastic shells

A state-based peridynamic formulation for linear elastic shells is presented. The emphasis is on introducing, possibly for the first time, a general surface based peridynamic model to represent the deformation characteristics of structures that have one geometric dimension much smaller than the other two. A new notion of curved bonds is exploited to cater for force transfer between the peridynamic particles describing the shell. Starting with the three dimensional force and deformation states, appropriate surface based force, moment and several deformation states are arrived at. Upon application on the curved bonds, such states yield the necessary force and deformation vectors governing the motion of the shell. By incorporating a shear correction factor, the formulation also accommodates analysis of shells that have higher thickness. In order to attain this, a consistent second order approximation to the complementary energy density is considered and incorporated in peridynamics via constitutive correspondence. Unlike the uncoupled constitution for thin shells, a consequence of a first order approximation, constitutive relations for thick shells are fully coupled in that surface wryness influences the in-plane stress resultants and surface strain the moments. Our proposal on the peridynamic shell theory is numerically assessed against simulations on static deformation of spherical and cylindrical shells, that of flat plates and quasi-static fracture propagation in a cylindrical shell. (C) 2016 Elsevier Ltd. All rights reserved.

Effect of W addition on the electrical switching of VO2 thin films

Vanadium Oxide has been a frontrunner in the field of oxide electronics because of its metal-insulator transition (MIT). The interplay of different structures of VO2 has played a crucial role in deciding the magnitude of the first order MIT. Substitution doping has been found to introduce different polymorphs of VO2. Hence the role of substitution doping in stabilizing the competing phases of VO2 in the thin film form remains underexplored. Consequently there have been reports both discounting and approving such a stabilization of competing phases in VO2. It is reported in the literature that the bandwidth of the hysteresis and transition temperature of VO2 can be tuned by substitutional doping of VO2 with W. In this work, we have adopted a novel technique called, Ultrasonic Nebulized Spray Pyrolysis of Aqueous Combustion Mixture (UNSPACM) to deposit VO2 and W- doped VO2 as thin films. XRD and Raman spectroscopy were used to investigate the role of tungsten on the structure of VO2 thin films. Morphology of the thin films was found to be consisting of globular and porous nanoparticles of size similar to 20nm. Transition temperature decreased with the addition of W. We found that for 2.0 at % W doping in VO2, the transition temperature has reduced from 68 degrees C to 25 degrees C. It is noted that W-doping in the process of reducing the transition temperature, alters the local structure and also increases room temperature carrier concentration. (c) 2016 Author(s).

Perturbational Finite Volume Method For The Solution of 2-D Navier-Stokes Equations on Unstructured and Structured Colocated Meshes

Based on the first-order upwind and second-order central type of finite volume( UFV and CFV) scheme, upwind and central type of perturbation finite volume ( UPFV and CPFV) schemes of the Navier-Stokes equations were developed. In PFV method, the mass fluxes of across the cell faces of the control volume (CV) were expanded into power series of the grid spacing and the coefficients of the power series were determined by means of the conservation equation itself. The UPFV and CPFV scheme respectively uses the same nodes and expressions as those of the normal first-order upwind and second-order central scheme, which is apt to programming. The results of numerical experiments about the flow in a lid-driven cavity and the problem of transport of a scalar quantity in a known velocity field show that compared to the first-order UFV and second-order CFV schemes, upwind PFV scheme is higher accuracy and resolution, especially better robustness. The numerical computation to flow in a lid-driven cavity shows that the under-relaxation factor can be arbitrarily selected ranging from 0.3 to 0. 8 and convergence perform excellent with Reynolds number variation from 102 to 104.

An analysis of collective damage for short fatigue cracks based on equilibrium of crack numerical density

Collective damage of short fatigue cracks was analyzed in the light of equilibrium of crack numerical density. With the estimation of crack growth rate and crack nucleation rate, the solution of the equilibrium equation was studied to reveal the distinct feature of saturation distribution for crack numerical density. The critical time that characterized the transition of short and long-crack regimes was estimated, in which the influences of grain size and grain-boundary obstacle effect were investigated. Furthermore, the total number of cracks and the first order of damage moment were discussed.

Plastic constitutive behavior of short-fiber/particle reinforced composites

A cylindrical cell model based on continuum theory for plastic constitutive behavior of short-fiber/particle reinforced composites is proposed. The composite is idealized as uniformly distributed periodic arrays of aligned cells, and each cell consists of a cylindrical inclusion surrounded by a plastically deforming matrix. In the analysis, the non-uniform deformation field of the cell is decomposed into the sum of the first order approximate field and the trial additional deformation field. The precise deformation field are determined based on the minimum strain energy principle. Systematic calculation results are presented for the influence of reinforcement volume fraction and shape on the overall mechanical behavior of the composites. The results are in good agreement with the existing finite element analyses and the experimental results. This paper attempts to stimulate the work to get the analytical constitutive relation of short-fiber/particle reinforced composites.

A Lagrangian for von Karman Equations of Large Deflection Problem of Thin Circular Plate

By the semi-inverse method proposed by He, a Lagrangian is established for the large deflection problem of thin circular plate. Ritz method is used to obtain an approximate analytical solution of the problem. First order approximate solution is obtained, which is similar to those in open literature. By Mathematica a more accurate solution can be deduced.

Kinetics of the chemical looping oxidation of H2 by a co-precipitated mixture of CuO and Al2O3

Chemical-looping combustion (CLC) has the inherent property of separating CO2 from flue gases. Instead of air, it uses an oxygen-carrier, usually in the form of a metal oxide, to provide oxygen for combustion. When used for the combustion of gaseous fuels, such as natural gas, or synthesis gas from the gasification of coal, the technique gives a stream of CO2 which, on an industrial scale, would be sufficiently pure for geological sequestration. An important issue is the form of the metal oxide, since it must retain its reactivity through many cycles of complete reduction and oxidation. Here, we report on the rates of oxidation of one constituent of synthesis gas, H2, by co-precipitated mixtures of CuO+Al2O3 using a laboratory-scale fluidised bed. To minimise the influence of external mass transfer, and also of errors in the measurement of [H2], particles sized to 355-500μm were used at low [H2], with the temperature ranging from 450 to 900°C. Under such conditions, the reaction was slow enough for meaningful measurements of the intrinsic kinetics to be made. The reaction was found to be first order with respect to H2. Above ∼800°C, the reaction of CuO was fast and conformed to the shrinking core mechanism, proceeding via the intermediate, Cu2O, in: 2CuO+H2→Cu2O+H2O, ΔH1073 K0=- 116.8 kJ/mol; Cu2O+H2→2Cu+H2O, ΔH1073 K0-80.9 kJ/mol. After oxidation of the products Cu and Cu2O back to CuO, the kinetics in subsequent cycles of chemical looping oxidation of H2 could be approximated by those in the first. Interestingly, the carrier was found to react at temperatures as low as 300°C. The influence of the number of cycles of reduction and oxidation is explored. Comparisons are drawn with previous work using reduction by CO. Finally, these results indicate that the kinetics of reaction of the oxygen carrier with gasifier synthesis gases is very much faster than rates of gasification of the original fuel. © 2010 The Institution of Chemical Engineers.

Hydrodynamic interaction between two vertical cylinders in water waves

The hydrodynamic interaction between two vertical cylinders in water waves is investigated based on the linearized potential flow theory. One of the two cylinders is fixed at the bottom while the other is articulated at the bottom and oscillates with small amplitudes in the direction of the incident wave. Both the diffracted wave and the radiation wave are studied in the present paper. A simple analytical expression for the velocity potential on the surface of each cylinder is obtained by means of Graf's addition theorem. The wave-excited forces and moments on the cylinders, the added masses and the radiation damping coefficients of the oscillating cylinder are all expressed explicitly in series form. The coefficients of the series are determined by solving algebraic equations. Several numerical examples are given to illustrate the effects of various parameters, such as the separation distance, the relative size of the cylinders, and the incident angle, on the first-order and steady second-order forces, the added masses and radiation-damping coefficients as well as the response of the oscillating cylinder.

Vibration analysis of laminated shells using a refined shear deformation theory

A previously published refined shear deformation theory is used to analyse free vibration of laminated shells. The theory includes the assumption that the transverse shear strains across any two layers are linearly dependent on each other. The theory has the same dependent variables as first-order shear deformation theory, hut the set of governing differential equations is of twelfth order. No shear correction factors are required. Free vibration of symmetric cross-ply laminated cylindrical shells, symmetric and antisymmetric cross-ply cylindrical panels is calculated. The numerical results are in good agreement with those from three-dimensional elasticity theory.

Modeling, numerical methods, and simulation for particle-fluid two-phase flow problems

In this paper, we study the issues of modeling, numerical methods, and simulation with comparison to experimental data for the particle-fluid two-phase flow problem involving a solid-liquid mixed medium. The physical situation being considered is a pulsed liquid fluidized bed. The mathematical model is based on the assumption of one-dimensional flows, incompressible in both particle and fluid phases, equal particle diameters, and the wall friction force on both phases being ignored. The model consists of a set of coupled differential equations describing the conservation of mass and momentum in both phases with coupling and interaction between the two phases. We demonstrate conditions under which the system is either mathematically well posed or ill posed. We consider the general model with additional physical viscosities and/or additional virtual mass forces, both of which stabilize the system. Two numerical methods, one of them is first-order accurate and the other fifth-order accurate, are used to solve the models. A change of variable technique effectively handles the changing domain and boundary conditions. The numerical methods are demonstrated to be stable and convergent through careful numerical experiments. Simulation results for realistic pulsed liquid fluidized bed are provided and compared with experimental data. (C) 2004 Elsevier Ltd. All rights reserved.

Vibration analysis of laminated plates using a refined shear deformation-theory

A previously published discrete-layer shear deformation theory is used to analyze free vibration of laminated plates. The theory includes the assumption that the transverse shear strains across any two layers are linearly dependent on each other. The theory has the same dependent variables as first order shear deformation theory, but the set of governing differential equations is of twelfth order. No shear correction factors are required. Free vibration of simply supported symmetric and antisymmetric cross-ply plates is calculated. The numerical results are in good agreement with those from three-dimensional elasticity theory.