On the radial oscillations of incompressible, isotropic, elastic and limited elastic thick-walled tubes

The finite amplitude, free radial oscillations of a thick-walled circular cylindrical tube are studied for an arbitrary incompressible, isotropic and homogeneous rubber-like material having limiting molecular chain extensibility. First, based on classical results for hyperelastic tubes, some results for thick-walled Mooney-Rivlin tubes are described graphically in the phase plane. Then the periodicity of the finite amplitude, free oscillations of a general limited elastic, thick-walled tube is studied; and some analytical results for the Gent model are illustrated in several numerical examples. Results for thick-walled Gent tubes are compared with those for corresponding Mooney-Rivlin tubes; and the motion of thin-walled Gent tubes is illustrated in the phase plane. Physical conclusions are presented. The period of small amplitude oscillations of an arbitrary elastic or limited elastic tube is derived from relations obtained by a linearization of a general class of equations of which the tube problem is a special case. Classical results of the linear theory are thereby recovered and compared with results for Mooney-Rivlin and Gent tubes.     

On dual-grid level-set method for contact line modeling during impact of a droplet on hydrophobic and superhydrophobic surfaces

In this paper, a numerical methodology for modeling contact line motion in a dual-grid level-set method (DGLSM) – solved on a uniform grid for interface which is twice that for the flow equations – is presented. A quasi-dynamic contact angle model – based on experimental inputs – is implemented to model the dynamic wetting of a droplet, impacting on a hydrophobic or a superhydrophobic surface. High-speed visualization experiments are also presented for the impact of a water droplet on hydrophobic surfaces, with non-bouncing at smaller and bouncing at larger impact velocity. The experimental results for temporal variation of the droplet shapes, wetted-diameter and maximum height of the droplet match very well with the DGLSM based numerical results. The validation of the numerical results is also presented with already published experimental results, for the non-bouncing on a hydrophobic and bouncing on a superhydrophobic surface, at a constant impact velocity. Finally, a qualitative as well as quantitative performance of the DGLSM as compared to the traditional level set method (LSM) is presented by considering our experimental results. The accuracy of the partially refined DGLSM is close to that of the fine-grid based LSM, at a computation cost which is close to that of the coarse-grid based LSM. The DGLSM is demonstrated as an improved LSM for the computational multi-fluid dynamics (CMFD) simulations involving contact line motion.     

On the elastic-plastic deformation of a sheet containing an edge crack

The solution of the planar tension and bending of an edge-cracked sheet of elastic-plastic material is given when the plastic deformation is represented by a Dugdale model. The analysis assumes conditions of generalized plane stress (for which the model of plastic relaxation is often a suitable one), but the usual transformation of elastic constants may be used to obtain the results also for plane-strain conditions. The method of solution involves the use of a Mellin transform and a Weiner-Hopf technique. Computed results for the size of the plastic zone and the opening at the crack tip are presented, and asymptotic results are obtained for small-scale and large-scale yielding. The results suggest that, when the material is constrained to fracture close to its ultimate tensile stress, the extra severity of a surface flaw compared with a corresponding internal crack is significantly greater than that predicted by linear elastic fracture mechanics.     

Finite amplitude spherically symmetric wave propagation in a prestressed hyperelastic shell

Spherically symmetric finite amplitude wave propagation in a prestressed compressible hyperelastic spherical shell is considered. The prestress results from quasi-static application of internal pressure and a numerical solution for this elastostatic problem is obtained first. Dynamic change of the internal pressure results in the propagation of a spherically symmetric wave. A Godunov type finite difference scheme is proposed for the solution of the wave propagation problem and numerical results, which are valid until the first reflection, are presented for a particular isotropic strain energy function and for the special cases of sudden removal and sudden increase of the internal pressure.     

Crack problem for superconducting strip with finite thickness

The inclined crack problems are considered for a thin strip and a strip with finite thickness in a perpendicular magnetic field. The critical current density is assumed to be a constant. The crack orientation is varied and the effect of crack on the magnetic field distribution is neglected. Based on the analytical results and variational inequality, the field and current distributions are computed for both thin strip and strip with finite thickness cases, respectively. Then, the stress intensity factors at the crack tip are determined using the finite element method for magnetic field loads. The numerical results are presented for different inclined crack angles, magnetization processes and geometry parameters of the strip. The results show that the fracture behavior of the strip with finite thickness is more complicated than that of the thin strip. With the numerical results, we can predict the largest possibility of cracking as the strip is in an external field.     

Buckling and vibration of a thin tensioned sheet with an elliptical hole

The out-of-plane deflection which can be developed in a thin tensioned sheet containing a central opening is discussed. It is noted that this behavior is of practical importance because it can alter the stress concentration and fatigue properties. The results of an experimental and analytical study of the buckling and vibration behavior of tensioned sheets with an elliptical opening are presented. Data from buckling and vibration experiments for several opening shapes ranging between a crack and a circle are presented. A generalized Galerkin method for solving the governing differential equations is used and the results obtained are presented in terms of buckling loads and natural frequencies. The experimental-buckling-stress estimates are correlated with the analytical results. The data from the vibration experiments indicate that the load-frequency behavior is dependent on the size of the test specimen and the opening shape. The analytical results display a load-frequency behavior which is experimentally observed for a narrow elliptical hole; i.e., for a small value of the ratio of minor axis to major axis. Some features of the analytical and experimental vibration results differ substantially, however. The analysis indicates that the frequency decreases to zero as the applied load approaches the buckling load. The experimental results do not exhibit this behavior. The sources of this difference in behavior are discussed.     

圆柱壳波动振动响应分析

采用Flugge的壳体理论,计算出有阻尼圆柱壳的固有频率并且推导出波动振动响应的模态解,得到了波动共振转速的计算公式。为了提高计算精度,采用高斯迭代法解计算固有频率的超越函数方程;计算了悬臂圆柱壳的固有频率及共波动响应,绘制了振型图,并与NASA试验值作了对比。     

Radiation effects on mixed convection about a cone embedded in a porous medium filled with a nanofluid

The problem of steady, laminar, mixed convection boundary-layer flow over a vertical cone embedded in a porous medium saturated with a nanofluid is studied, in the presence of thermal radiation. The model used for the nanofluid incorporates the effects of Brownian motion and thermophoresis with Rosseland diffusion approximation. The cone surface is maintained at a constant temperature and a constant nanoparticle volume fraction. The resulting governing equations are non-dimensionalized and transformed into a non-similar form and then solved by Keller box method. A comparison is made with the available results in the literature, and our results are in very good agreement with the known results. A parametric study of the physical parameters is made and a representative set of numerical results for the local Nusselt and Sherwood numbers are presented graphically. Also, the salient features of the results are analyzed and discussed.     

Radiation Effects on Mixed Convection over a Wedge Embedded in a Porous Medium Filled with a Nanofluid

The problem of steady, laminar, mixed convection boundary-layer flow over an isothermal vertical wedge embedded in a porous medium saturated with a nanofluid is studied, in the presence of thermal radiation. The model used for the nanofluid incorporates the effects of Brownian motion and thermophoresis with Rosseland diffusion approximation. The wedge surface is maintained at a constant temperature and a constant nanoparticle volume fraction. The resulting governing equations are non-dimensionalized and transformed into a non-similar form and then solved by Keller box method. A comparison is made with the available results in the literature, and our results are in very good agreement with the known results. A parametric study of the physical parameters is made, and a representative set of numerical results for the velocity, temperature, and volume fraction, the local Nusselt and Sherwood numbers are presented graphically. The salient features of the results are analyzed and discussed.     

Griffith crack moving in a piezoelectric strip

Summary  The dynamic problem of an impermeable crack of constant length 2a propagating along a piezoelectric ceramic strip is considered under the action of uniform anti-plane shear stress and uniform electric field. The integral transform technique is employed to reduce the mixed-boundary-value problem to a singular integral equation. For the case of a crack moving in the mid-plane, explicit analytic expressions for the electroelastic field and the field intensity factors are obtained, while for an eccentric crack moving along a piezoelectric strip, numerical results are determined via the Lobatto–Chebyshev collocation method for solving a resulting singular integral equation. The results reveal that the electric-displacement intensity factor is independent of the crack velocity, while other field intensity factors depend on the crack velocity when referred to the moving coordinate system. If the crack velocity vanishes, the present results reduce to those for a stationary crack in a piezoelectric strip. In contrast to the results for a stationary crack, applied stress gives rise to a singular electric field and applied electric field results in a singular stress for a moving crack in a piezoelectric strip. Received 14 August 2001; accepted for publication 24 September 2002 The author is indebted to the AAM Reviewers for their helpful suggestions for improving this paper. The work was supported by the National Natural Science Foundation of China under Grant 70272043.     

Numerical comparison of conjugate and non-conjugate natural convection for internally heated semi-circular pools

In this numerical investigation, a detailed comparison of the conjugate and non-conjugate natural convection within a semi-cylindrical cavity has been presented. The cavity is assumed to be filled with a fluid containing uniformly distributed internal heating sources. The bottom circular wall of the cavity is taken to be thick with finite conductive properties, while the top wall is considered to be isothermal. The Navier–Stokes and energy equations are solved numerically by using the SIMPLER algorithm. A Rayleigh number range from 3.2×10⁶ to 3.2×10¹¹ has been investigated and the effects of solid-to-fluid conductivity ratios of 1.0, 5.0 and 23.0 have been analysed. The present numerical results for a semi-circular cavity with entirely isothermal walls are compared with known results from the open literature. It was found that these results for the non-conjugate problem are in very good agreement. The present results for a conjugate cavity show a remarkable difference from the non-conjugate analysis. The average Nusselt number for the solid–fluid interface shows a decrease while the top wall average Nu number has increased. It has been concluded that these effects increase for a system with a low solid-to-fluid conductivity ratio. It is evident from the present conjugate results that the assumption of isothermal enclosing walls gives somewhat different results when the walls are thick and the solid-to-fluid conductivity ratio is small.     

Conduction–radiation interaction in 3D irregular enclosures using the finite volume method

Interaction of conduction and radiation in 3D enclosures is carried out with a gray participating media. Application of block structured grid is shown with the finite volume method (FVM). Radiation modeling is performed with the FVM and is coupled with an ‘in-house’ code to solve the set of transport equations. The detailed numerical results are presented for a cubical and a cylindrical enclosure as these results are not available in the literature. The numerical simulation for the cylindrical enclosure is performed using a block-structured ‘O’ grid. Two additional geometries are considered in order to show the applicability of the present work. Results of temperature, radiative heat flux and total heat flux distributions are presented for different optical thicknesses, scattering albedoes, emissivities and conduction–radiation parameters. The 3D results are validated with the available 2D results or results with pure radiation problems as limiting cases.     

MHD boundary layer flow and heat transfer over a stretching sheet with induced magnetic field

In this paper, the problem of steady magnetohydrodynamic boundary layer flow and heat transfer of a viscous and electrically conducting fluid over a stretching sheet is studied. The effect of the induced magnetic field is taken into account. The transformed ordinary differential equations are solved numerically using the finite-difference scheme known as the Keller-box method. Numerical results are obtained for various values of the magnetic parameter, the reciprocal magnetic Prandtl number and the Prandtl number. The effects of these parameters on the flow and heat transfer characteristics are determined and discussed in detail. When the magnetic field is absent, the closed analytical results for the skin friction are compared with the exact numerical results. Also the numerical results for the heat flux from the stretching surface are compared with the results reported by other authors when the magnetic field is absent. It is found that very good agreement exists.     

High-order methods for differential equations with large first-derivative terms

A method is developed to solve elliptic singular perturbation problems. Examples are presented in one and two dimensions for both linear and non-linear problems. In particular, examples are presented for fluid flow problems with boundary layers. In the one-dimensional case an approximating equation is developed using just three points. The method first presented is a fourth-order approximation but is extended to become a higher-order method. Results are included for the fourth-, sixth-, eighth- and tenth-order methods. The results are first compared with results found by Segal in an article about elliptic singular perturbation problems. The elliptic singular perturbation problems are compared with a method by Il'in and also with central and backward difference schemes from Segal's article. There was only one case where the results in Segal's paper were as accurate as the results presented in this paper. However, in this case the method used by Segal did not give accurate values for a second problem presented. The results are also compared with results given by Spalding and by Christie. The method of this paper was also tested on the solution of some non-linear diffusion equations with concentration-dependent diffusion coefficients. The results were superior to results presented by Lee and by Schultz. Finally, the method is extended to several two-dimensional problems. The method developed in this paper is accurate, easy to use and can be generalized to other problems.     

Weissenberg effect and its dependence upon the experimental geometry

The flow of a second-order fluid with a free surface between two coaxially mounted cylinders of finite length, the inner one of which is rotating, is being studied. In the case of slow flow and small shear rates the flow can be divided into a primary flow in the plane perpendicular to the axis of rotation and a secondary flow in the meridional plane. These flow components are numerically calculated and the results are compared with the analytical results for the semi-infinite cylinder approximation. The influence of the finiteness of the cylinders (end effect) upon the free surface deformation is analysed. The numerical results for the secondary flow are compared with results obtained by flow visualisation.     

Forces on oscillating bodies in viscous fluid

This paper describes a method for determining the fluid forces on oscillating bodies in viscous fluid when the corresponding flow problem has been solved using the finite element method. These forces are characterized by the concept of added mass, added damping and added force. Numerical results are obtained for several example body shapes. Comparison is made with exact analytical results and other finite element results for the limiting cases of Stoke's flow and inviscid flow, and good agreement is obtained. The results for finite values of the body amplitude parameter β show the appearance of added force from the steady streaming component of the flow for asymmetric bodies. Results are also obtained for the associated flow where the fluid remote from a fixed body is oscillating.     

范式理论与平均法的等价性分析

分析了范式理论与平均法的等价性。得到的结论是：对含有一对纯虚根的二维非线性系统，使用两种方法得到的结果是等价的，并提供了两个算例来证实其结论的正确性。虽然分析是针对一类二维非线性系统，但其结论同样适合于高维非线性系统。     

Surface cracks in a plate of finite width under extension or bending

In this paper the problem of a finite plate containing collinear surface cracks is considered. The problem is solved by using the line spring model with plane elasticity and Reissner's plate theory. The main purpose of the study is to investigate the effect of interaction between two cracks or between cracks and stress-free plate boundaries on the stress intensity factors and to provide extensive numerical results which may be useful in applications. First, some sample results are obtained and are compared with the existing finite element results. Then the problem is solved for a single (internal) crack, two collinear cracks and two corner cracks for wide range of relative dimensions. Particularly in corner cracks the agreement with the finite element solution is surprisingly very good. The results are obtained for semielliptic and rectangular crack profiles which may, in practice, correspond to two limiting cases of the actual profile of a subcritically growing surface crack.     

Buckling of elastic-plastic cylindrical panel under axial compression

The buckling behaviour is investigated for an axially compressed elastic-plastic cylindrical panel of the type occurring in stiffened shells. The bifurcation stress is determined analytically and an asymptotically exact expansion is obtained for the initial post-bifurcation behaviour in the plastic range. For panels with small initial imperfections the behaviour is analysed asymptotically on the basis of the hypoelastic theory that results from neglecting the effect of elastic unloading. The imperfection-sensitivity of an elastic-plastic panel is also computed numerically by a linear incremental method, and the results are compared with the results of the asymptotic analysis. For a low hardening material the panel is found to be imperfection-sensitive in the whole range of curvatures considered, whereas for a high hardening material the panel is only imperfection-sensitive if the curvature exceeds a certain value.     

Water wave scattering by a submerged circular-arc-shaped plate

The problem of water wave scattering by a thin circular-arc-shaped plate submerged in infinitely deep water is investigated by linear theory. The circular-arc is not necessarily symmetric about the vertical through its center. The problem is formulated in terms of a hypersingular integral equation for a discontinuity of the potential function across the plate. The integral equation is solved approximately using a finite series involving Chebyshev polynomials of the second kind. The unknown constants in the finite series are determined numerically by using the collocation and the Galerkin methods. Both the methods ultimately produce very accurate numerical estimates for the reflection coefficient. The numerical results are depicted graphically against the wave number for a variety of configurations of the arc. Some results are compared with known results available in the literature and good agreement is achieved. The suitability of using a circular-arc-shaped plate as an element of a water wave lens has also been discussed on the basis of the present numerical results.