Experimental investigation of natural convective flows in slowly rotating wide spherical shells

Experimental results are presented on natural convection in a spherical shell of inner and outer radii r ₁ = 14 mm and r ₂ = 35 mm, with the inner sphere cooled and the outer sphere heated. The fluids filling the shell are two different silicon oils having Prandtl numbers 39 and 233. Both spheres are fixed together and can be rotated. In the studied regime, both Coriolis and centrifugal forces become significant. For sufficiently small Rayleigh numbers the resulting flow pattern is axisymmetric and steady, consisting of a plume descending from the south pole of the inner sphere, and returning in the equatorial regions. For greater Rayleigh numbers the flow becomes non-axisymmetric, with azimuthal modes m = 2 to 4 arising. We map out the critical Rayleigh number for the onset of these different modes, and consider how they vary with increasingly rapid overall rotation. Detailed flow measurements are done by converting a standard 2D particle image velocimetry system into a scanning quasi-3D PIV system.     

Numerical investigation of two-phase rotating flow

Finite difference solutions of the two-fluid equations of motion for a particle (droplet)-fluid mixture in a rotating finite axisymmetric cylinder are presented. The numerical method, which can be regarded as an extension of the Harlow & Amsden approach, employs forward time and centred space discretization and treats implicitly the pressure, Coriolis and volume flux terms. The computed flow fields are examined via a detailed comparison to previous analytic approximations, which illuminates both the physical and numerical aspects and the validity of these approximations.     

Multi-symplectic methods for membrane free vibration equation

In this paper,the multi-symplectic formulations of the membrane free vi- bration equation with periodic boundary conditions in Hamilton space are considered. The complex method is introduced and a semi-implicit twenty-seven-points scheme with certain discrete conservation laws—a multi-symplectic conservation law(CLS),a local energy conservation law(ECL)as well as a local momentum conservation law(MCL)—is constructed to discrete the PDEs that are derived from the membrane free vibra- tion equation.The results of the numerical experiments show that the multi-symplectic scheme has excellent long-time numerical behavior.     

An investigation into the spreading of a thin liquid drop under gravity on a slowly rotating disk

The axisymmetric spreading of a thin liquid drop under the influence of gravity and rotation is investigated. The effects of the Coriolis force and surface tension are ignored. The Lie group method is used to analyse the non-linear diffusion-convection equation modelling the spreading of the liquid drop under gravity and rotation. A stationary group invariant solution is obtained. The case when rotation is small is considered next. A straightforward perturbation approach is used to determine the effects of the small rotation on the solution given for spreading under gravity only. Over a short period of time no real difference is observed between the approximate solution and the solution for spreading under gravity only. After a long period of time, the approximate solution tends toward a dewetting solution. We find that the approximate solution is valid only in the interval t∈[0,t∗), where t∗ is the time when dewetting takes place. An approximation to t∗ is obtained.     

Investigation of the effect of the Coriolis force on a thin fluid film on a rotating disk

The effect of the Coriolis force on the evolution of a thin film of Newtonian fluid on a rotating disk is investigated. The thin-film approximation is made in which inertia terms in the Navier–Stokes equation are neglected. This requires that the thickness of the thin film be less than the thickness of the Ekman boundary layer in a rotating fluid of the same kinematic viscosity. A new first-order quasi-linear partial differential equation for the thickness of the thin film, which describes viscous, centrifugal and Coriolis-force effects, is derived. It extends an equation due to Emslie et al. [J. Appl. Phys. 29, 858 (1958)] which was obtained neglecting the Coriolis force. The problem is formulated as a Cauchy initial-value problem. As time increases the surface profile flattens and, if the initial profile is sufficiently negative, it develops a breaking wave. Numerical solutions of the new equation, obtained by integrating along its characteristic curves, are compared with analytical solutions of the equation of Emslie et al. to determine the effect of the Coriolis force on the surface flattening, the wave breaking and the streamlines when inertia terms are neglected.     

New expressions for the radiation force function of spherical targets in stationary and quasi-stationary waves

A new expression for the radiation force function – which is the radiation force per unit energy density and unit cross-sectional surface area – for spheres in a stationary (or standing) and quasi-stationary wave is obtained based on the far-field acoustic scattering field. The radiation force function formulation has been simplified mathematically and improved into a more general form. Numerical results are presented for rigid and elastic spheres, air bubbles in water as well as liquid drops in air to illustrate the theory. It is demonstrated that expressions for the radiation force functions obtained from the far-field derivation approach are equivalent to those obtained from the near-field-based derivation.     

Elastic instabilities for strain-stiffening rubber-like spherical and cylindrical thin shells under inflation

This paper is concerned with investigation of the effects of strain-stiffening on the classical limit point instability that is well-known to occur in the inflation of internally pressurized rubber-like spherical thin shells (balloons) and circular cylindrical thin tubes composed of incompressible isotropic non-linearly elastic materials. For a variety of specific strain-energy densities that give rise to strain-stiffening in the stress-stretch response, the inflation pressure versus stretch relations are given explicitly and the non-monotonic character of the inflation curves is examined. While such results are known for constitutive models that exhibit a gradual stiffening (e.g. exponential and power-law models), our primary focus is on materials that undergo severe strain-stiffening in the stress-stretch response. In particular, we consider two phenomenological constitutive models that reflect limiting chain extensibility at the molecular level. It is shown that for materials with sufficiently low extensibility no limit point instability occurs and so stable inflation is then predicted for such materials. Potential applications of the results to the biomechanics of soft tissues are indicated.     

New expression for collision efficiency of spherical nanoparticles in Brownian coagulation

The collision efficiency of dioctyl phthalate nanoparticles in Brownian coag- ulation has been studied. A set of collision equations is solved numerically to find the relationship between the collision efficiency and the particle radius varying in the range of 50 nm to 500 nm in the presence of Stokes resistance, lubrication force, van der Waals force, and elastic deformation force. The calculated results are in agreement with the experimental data qualitatively. The results show that the collision efficiency decreases with the increase of the particle radii from 50 nm to 500 nm. Based on the numerical data, a new expression for collision efficiency is presented.     

The general solution of bending of a spherical thin shallow shell with a circular hole at the center under arbitrary transverse loads

Basing several suggestions appeared in [1], we find out the general solution of the bending of a spherical thin shallow shell with a circular hole at the center. As we know, very few papers had touched upon these problems.     

Two-dimensional integral equation for a thin wind turbine blade rotating in the round tunnel

The paper presents a mathematical treatment of the aerodynamic problem about a thin wind turbine blade rotating in the round tunnel. The radius of the blade is almost the radius of the tunnel. This permits formulation of the boundary condition on the tip vortex line to the simple slip condition over the surface of the tunnel. By applying a standard technique from potential theory, the problem is reduced to a two-dimensional integral equation whose kernel is connected with a special Green's function satisfying the homogeneous Neumann boundary condition on the tunnel surface. This Green's function is constructed in an explicit analytical form.     

旋转机械系统多自由度非线性动力学数值分析

首先将转子系统的动力响应问题归结为2n维未知向量v的一阶非线性动力学方程dv／dt=Hv＋f（v,t）,并给出了求解这一方程的一次近似式法和三次多项式迫近法。在非稳态、非线性油膜力等作用下,以刚性Jeffcott转子与112个自由度的汽轮发电机组低压转子系统为例,用上述求解方法分析了它们的动力响应及非线性动力学特性;其问,还将计算结果与Runge—Kutta法、Newmark法的相应结果进行了比较,并深入讨论了数值稳定性问题。汽轮发电机组的算例表明对一些具有较复杂的非线性右端项,、同时规模又较大的问题,如果采用四阶Runge—Kutta法,才算几步就因数值骤然增大而失控;但若用同样步长的一次近似式,由于它是一种显式的无条件稳定算式,则计算过程迅速且结果合理可靠。     

Neural networks for tip correction of spherical indentation curves from bulk metals and thin metal films

Spherical indentation is widely used to determine a variety of important mechanical properties from small volumes. However, the available nanoindenter tips mostly deviate from the perfect spherical shape making the application of analysis methods developed for perfect spheres uncertain. In this paper, neural network-based methods are presented that are used to correct force-depth curves measured with such indenter tips. Finite element simulations for imperfect and perfect spherical tips with varying material behaviour are used to train the neural networks, which solve the inverse problem of mapping the true tip shape and the measured force-depth curve to one that corresponds to a perfect spherical indenter. Solutions are provided for bulk materials and thin films. The method has been verified experimentally on nanocrystalline nickel and a copper film on a titanium substrate for different spherical tips.     

A model for computational investigation of elasto-plastic normal and tangential contact considering adhesion effect

With the rapid development of Micro-Electro-Mechanical System (MEMS), we enter a field in which the surface effects have dominated many of the micro-scale phenomena, and the adhesive contact is one of the focuses. In this paper, a feasible model for finite element computation is presented via a macroscopic and microscopic combination approach, in which the adhesive forces are simulated by some non-linear spring elements considering the softening stage. Two basic problems concerning the adhesion effect were considered; through specific quantitative analysis, the results show a consistency with the current elastic continuum theories of adhesion and a brief investigation into the effects of adhesion on plastic deformation and tangential contact will be carried out as well. The project supported by the National Natural Science Foundation of China (10172050, 90205022) and Key Grant Project of Chinese MoE (0306)     

薄板几何非线性中的精化元方法及膜闭锁问题

基于放松单元间协调条件的大变形变分原理和全局拉格朗日方法，推导了几何非线性精化三角形薄板单元。对几何刚度矩阵，通过引入特殊的单元位移函数，有效地消除了薄板弯曲问题中伴生的膜闭锁现象。数值结果表明该单元在几何非线性分析中既能消除膜闭锁又具有较高精度。     

Remarks of some problems for rectangular thin plates with free edges on elastic foundations

For the bending,stability and vibrations of rectangular thin plates with free edges onelastic foundations,in this paper we give a flexural function which exactly satisfies not onlyall the boundary conditions on free edges but also the conditions at free corner points.Applying energy variation principle,we give equations defining parameters in flexuralfunction,stability equation,frequency equation,and general formulae of minimum criticalforce and minimum eigenfrequency as well.     

Multi‐objective optimization of a rotating cooling channel with staggered pin‐fins for heat transfer augmentation

A rotating channel with staggered pin‐fins is formulated numerically and optimized for performance (heat transfer/required pumping power) using a Kriging meta‐model and hybrid multi‐objective evolutionary algorithm. Two design variables related to cooling channel height, pin diameter, and spacing between the pins are selected for optimization, and two‐objective functions related to the heat transfer and friction loss are employed. A design of experiment is performed, and 20 designs are generated by Latin hypercube sampling. The objective function values are evaluated using a Reynolds‐averaged Navier–Stokes solver, and a Kriging model is constructed to obtain a Pareto‐optimal front through a multi‐objective evolutionary algorithm. Rotation in a cooling channel with staggered pin‐fins induces Coriolis force that causes a heat transfer discrepancy between the trailing (pressure) and leading (suction) surfaces, with a higher Nusselt number on the trailing surface. The tradeoff between the two competing objective functions is determined, and the distribution of the Pareto‐optimal solutions in the design space is discussed through k‐means clustering. In the optimal designs, with a decrease in spacing between the pins, heat transfer is enhanced whereas friction loss is increased. Copyright © 2011 John Wiley & Sons, Ltd.     

Experimental investigation of a maximum entropy assumption for acceleration terms within a poly‐disperse moment framework

Moment transport methods are being developed to model poly‐dispersed multiphase flows by transporting statistical moments of the particle size–velocity joint probability density function (JPDF). A common feature of these methods is the requirement to reproduce or approximate the form of the JPDF from the transported moments for calculation of body force terms and other source terms. This paper examines the application of a maximum entropy technique against phase Doppler anemometry data sets from an electrostatically charged kerosene spray and also an automotive pressure swirl atomizer. An assessment of which moments are required to reproduce the JPDFs using a maximum entropy assumption to a sufficient level of accuracy is made. It is found that it is possible to reproduce the JPDFs to a high level of accuracy using a large number of moments; however, this incurs large computational overheads. If the moments to be transported are chosen on the basis of physical reasoning (such as the relationship between size and velocity due to drag) it is possible to reduce the number of moments to those which would be conserved via balance equations. This permits an approximation to the JPDF commensurate with the closure level of the moment transport method and thus the closure model method is naturally scalable with the degree of information from available conservation equations. Copyright © 2008 John Wiley & Sons, Ltd.     

商用飞机撞击核电站屏蔽厂房数值模拟

为研究大型商用飞机撞击核电站屏蔽厂房的毁伤特性,建立了Boeing 767商用飞机和双钢板混凝土屏蔽厂房的有限元模型,模拟两者的相互撞击作用过程,得到飞机不同部位(机身、发动机)对屏蔽厂房的撞击力时程曲线,对每一个作用部分给出简化撞击力曲线和作用面积,确定了撞击力分布形式。结果分析表明:飞机轴向网格尺寸对撞击力影响较大;屏蔽厂房被撞击部位变形明显,其他区域变形较小;撞击速度对撞击作用时间影响较小,而对结构响应位移影响很大,撞击力合力随着撞击速度的降低迅速下降。     

Multiphysics numerical investigation of photothermal self-driving characteristics of SiO2@Au Janus microparticles for drug delivery

The photothermal self-driving process of Janus microparticles has wide application prospects in the fields of biomedicine. Since silica and gold have good biocompatibility and high photothermal conversion efficiency, the SiO₂@Au Janus microparticles are widely used as drug carriers. Based on the multiphysics coupling method, the photothermal self-driving process of SiO₂@Au Janus microparticles was investigated, wherein the substrate was SiO₂ particles and one side of the particles was coated with gold film. Under a continuous wave laser with irradiation of 20 W/cm², the distance covered by the Janus particles was increased by increasing the thickness of the gold film and reducing the size of the SiO₂ particles; the self-driving characteristics of the Janus particles were controlled substantially by increasing the intensity of the incident laser. Based on the simulation results, the thermophoretic motion and Brownian motion of particles can be measured by comparing the absolute values of the thermophoretic force impulse, Brownian force impulse, and drag force impulse. The Brownian force acting on Janus microparticles under low laser power cannot be ignored. Furthermore, the minimum laser power required for Janus particles to overcome Brownian motion was calculated. The results can effectively guide the design of Janus particles in biomedicine and systematically analyze the mechanism of particle thermophoretic motion during drug delivery.