层状弹性半空间非轴对称动力问题的奇异解

在柱坐标系下，利用关于方位角的Ｆｏｕｒｉｅｒ变换及关于径向的Ｈａｎｋｅｌ变换，将弹性力学基本方程组转化为非齐次的一阶常微分方程组的标准形式．采用求解微分方程组的矩阵法，建立了介质层的传递矩阵．由层间完全接触条件，导出了在任意埋藏源作用下层状弹性半空间频域奇异解，时域奇异解可通过关于频率的Ｆｏｕｒｉｅｒ积分得到．该方法可应用到固体、流体层的情况     

Numerical Analysis for the Bingham—Papanastasiou Fluid Flow Over a Rotating Disk

This study is conducted to investigate the Bingham—Papanastasiou fluid flow driven by a rotating infinite disk. The Bingham—Papanastasiou model is a modification of the Bingham plastic model, which is developed by introducing a continuation parameter to overcome its discontinuity. The von K´armán similarity solution is used to transform the flow equations from ordinary differential equations to a nonlinear system of partial differential equations, which is solved numerically. The effect of the Bingham flow parameters on the radial, tangential, and axial velocities, pressure, and radial and tangential skin friction coefficients is discussed.     

Nonlinear dynamic responses of sandwich functionally graded porous cylindrical shells embedded in elastic media under 1:1 internal resonance

In this article, the nonlinear dynamic responses of sandwich functionally graded(FG) porous cylindrical shell embedded in elastic media are investigated. The shell studied here consists of three layers, of which the outer and inner skins are made of solid metal, while the core is FG porous metal foam. Partial differential equations are derived by utilizing the improved Donnell's nonlinear shell theory and Hamilton's principle. Afterwards, the Galerkin method is used to transform the governing equations into nonlinear ordinary differential equations, and an approximate analytical solution is obtained by using the multiple scales method. The effects of various system parameters,specifically, the radial load, core thickness, foam type, foam coefficient, structure damping,and Winkler-Pasternak foundation parameters on nonlinear internal resonance of the sandwich FG porous thin shells are evaluated.     

Scattering by an anisotropic circle

The scattering by a circle is considered when the outside medium is isotropic and the inside medium is anisotropic (orthotropic). The problem is a scalar one and is phrased as a scattering problem for elastic waves with polarization out of the plane of the circle (SH wave), but the solution is with minor modifications valid also for scattering of electromagnetic waves. The equation inside the circle is first transformed to polar coordinates and it then explicitly contains the azimuthal angle through trigonometric functions. Making an expansion in a trigonometric series in the azimuthal coordinate then gives a coupled system of ordinary differential equations in the radial coordinate that is solved by power series expansions. With the solution inside the circle complete the scattering problem is solved essentially as in the classical case. Some numerical examples are given showing the influence of anisotropy, and it is noted that the effects of anisotropy are generally strong except at low frequencies where the dominating scattering only depends on the mean stiffness and not on the degree of anisotropy.     

Exact Analytic Solutions of the Nonlinear Long-Wave Equations in the Case of Axisymmetric Fluid Vibrations in a Parabolic Rotating Vessel

A class of exact analytic solutions of the system of nonlinear long-wave equations is found. This class corresponds to the axisymmetric vibrations of an ideal incompressible homogeneous fluid in a rotating vessel in the shape of a paraboloid of revolution. The radial velocity of these motions is a linear function, and the azimuthal velocity and free surface displacements are polynomials in the radial coordinate with time-dependent coefficients. The nonlinear vibration frequency is equal to the frequency of the lowest mode of linear axisymmetric standing waves in the parabolic vessel.     

Optimal transient growth in turbulent pipe flow

The optimal transient growth process of perturbations driven by the pressure gradient is studied in a turbulent pipe flow. A new computational method is proposed, based on the projection operators which project the governing equations onto the subspace spanned by the radial vorticity and radial velocity. The method is validated by comparing with the previous studies. Two peaks of the maximum transient growth amplification curve are found at different Reynolds numbers ranging from 20 000 to 250 000. The optimal flow structures are obtained and compared with the experiments and DNS results. The location of the outer peak is at the azimuthal wave number n=1, while the location of the inner peak is varying with the Reynolds number. It is observed that the velocity streaks in the buffer layer with a spacing of 100δ_v are the most amplified flow structures. Finally, we consider the optimal transient growth time and its dependence on the azimuthal wave length. It shows a self-similar behavior for perturbations of different scales in the optimal transient growth process.     

Exact Solution of the Heat Transfer Problem for a Rotating Disk under Uniform Jet Impingement

An exact solution of the heat transfer problem for a uniform air stream impinging on a rotating disk is found. By introducing self-similar radial velocity and temperature profiles, the problem is reduced to a system of ordinary differential equations which are solved numerically. The Nusselt numbers are calculated for Prandtl numbers equal to 1 and 0.71 and various ratios of the free-stream velocity to the disk rotation velocity. The limits of the flow regime in which the heat transfer is determined solely by the impact jet parameters are found. The results are compared with experimental data for the stagnation point.     

Estimates of azimuthal numbers associated with elementary elliptic cylinder wave functions

The paper deals with issues related to the construction of solutions, 2 π-periodic in the angular variable, of the Mathieu differential equation for the circular elliptic cylinder harmonics, the associated characteristic values, and the azimuthal numbers needed to form the elementary elliptic cylinder wave functions. A superposition of the latter is one possible form for representing the analytic solution of the thermoelastic wave propagation problem in long waveguides with elliptic cross-section contour. The classical Sturm-Liouville problem for the Mathieu equation is reduced to a spectral problem for a linear self-adjoint operator in the Hilbert space of infinite square summable two-sided sequences. An approach is proposed that permits one to derive rather simple algorithms for computing the characteristic values of the angular Mathieu equation with real parameters and the corresponding eigenfunctions. Priority is given to the application of the most symmetric forms and equations that have not yet been used in the theory of the Mathieu equation. These algorithms amount to constructing a matrix diagonalizing an infinite symmetric pentadiagonal matrix. The problem of generalizing the notion of azimuthal number of a wave propagating in a cylindrical waveguide to the case of elliptic geometry is considered. Two-sided mutually refining estimates are constructed for the spectral values of the Mathieu differential operator with periodic and half-periodic (antiperiodic) boundary conditions.     

Instability waves and low-frequency noise radiation in the subsonic chevron jet

Spatial instability waves associated with lowfrequency noise radiation at shallow polar angles in the chevron jet are investigated and are compared to the round counterpart. The Reynolds-averaged Navier–Stokes equations are solved to obtain the mean flow fields, which serve as the baseflow for linear stability analysis. The chevron jet has more complicated instability waves than the round jet, where three types of instability modes are identified in the vicinity of the nozzle, corresponding to radial shear, azimuthal shear,and their integrated effect of the baseflow, respectively. The most unstable frequency of all chevron modes and round modes in both jets decrease as the axial location moves downstream. Besides, the azimuthal shear effect related modes are more unstable than radial shear effect related modes at low frequencies. Compared to a round jet, a chevron jet reduces the growth rate of the most unstable modes at downstream locations. Moreover, linearized Euler equations are employed to obtain the beam pattern of pressure generated by spatially evolving instability waves at a dominant low frequency St = 0.3, and the acoustic efficiencies of these linear wavepackets are evaluated for both jets. It is found that the acoustic efficiency of linear wavepacket is able to be reduced greatly in the chevron jet, compared to the round jet.     

Darcy-Forchheimer flows of copper and silver water nanofluids between two rotating stretchable disks

This investigation describes the nanofluid flow in a non-Darcy porous medium between two stretching and rotating disks. A nanofluid comprises of nanoparticles of silver and copper. Water is used as a base fluid. Heat is being transferred with thermal radiation and the Joule heating. A system of ordinary differential equations is obtained by appropriate transformations. Convergent series solutions are obtained. Effects of various parameters are analyzed for the velocity and temperature. Numerical values of the skin friction coefficient and the Nusselt number are tabulated and examined. It can be seen that the radial velocity is affected in the same manner with both porous and local inertial parameters. A skin friction coefficient depicts the same impact on both disks for both nanofluids with larger stretching parameters.     

Low-frequency flow fluctuations downstream of a sudden enlargement in a circular tube

The low-frequency structure of the turbulent flow downstream of a sudden enlargement in a circular tube, after the flow has reattached to the wall, is studied. In a cross-section of the wide part of the tube beyond the reattachment point the total pressure is synchronously measured using radial and azimuthal rakes equipped with Pitot probes. The values of the coefficient of spatial correlation between the total-pressure fluctuations are considerably different for the azimuthal and radial directions. Under the assumption that the total-pressure loss in a fluid particle is conserved after flow reattachment to the wall, the azimuthal and radial components of the total-pressure fluctuations are estimated on the basis of the measured data.     

Exact solutions for axisymmetric flexural free vibrations of inhomogeneous circular Mindlin plates with variable thickness

Circular plates with radially varying thickness, stiffness, and density are widely used for the structural optimization in engineering. The axisymmetric flexural free vibration of such plates, governed by coupled differential equations with variable coefficients by use of the Mindlin plate theory, is very difficult to be studied analytically.In this paper, a novel analytical method is proposed to reduce such governing equations for circular plates to a pair of uncoupled and easily solvable differential equations of the Sturm-Liouville type. There are two important parameters in the reduced equations.One describes the radial variations of the translational inertia and flexural rigidity with the consideration of the effect of Poisson's ratio. The other reflects the comprehensive effect of the rotatory inertia and shear deformation. The Heun-type equations, recently well-known in physics, are introduced here to solve the flexural free vibration of circular plates analytically, and two basic differential formulae for the local Heun-type functions are discovered for the first time, which will be of great value in enriching the theory of Heun-type differential equations.     

MHD axisymmetric flow of third grade fluid between porous disks with heat transfer

The magnetohydrodynamic(MHD) flow of the third grade fluid between two permeable disks with heat transfer is investigated.The governing partial differential equations are converted into the ordinary differential equations by suitable transformations.The transformed equations are solved by the homotopy analysis method(HAM).The expressions for square residual errors are defined,and the optimal values of convergencecontrol parameters are selected.The dimensionless velocity and temperature fields are examined for various dimensionless parameters.The skin friction coefficient and the Nusselt number are tabulated to analyze the effects of dimensionless parameters.     

极坐标系下的Legendre谱元方法求解Poisson-型方程

r=0处的坐标奇异性是求解极坐标下Poisson-型方程的关键。本文提出一种极坐标系下基于Galerkin变分的Legendre谱元方法用于求解圆形区域内的Poisson-型方程,物理区域的径向和周向划分若干单元,计算单元均采用Legendre多项式展开;圆心所在单元的径向使用LGR(Legendre Gauss Radau)积分点,其他单元径向使用LGL(Legendre Gauss Lobatto)积分点,从而避免了极点处1/r坐标奇异性,周向单元均采用LGL积分点。利用区域分解技术,可以避免节点在极点附近聚集;最后求解了多个Dirichlet或Neumann边界条件下的Poisson-型方程算例。数值结果表明,谱元方法具有很高的精度。     

Three-dimensional free vibration analysis of four-parameter continuous grading fiber reinforced cylindrical panels resting on Pasternak foundations

This research investigates three-dimensional free vibration analysis of four-parameter continuous grading fiber reinforced (CGFR) cylindrical panels resting on Pasternak foundations by using generalized power-law distribution. The functionally graded orthotropic panel is simply supported at the edges, and it is assumed to have an arbitrary variation of matrix volume fraction in the radial direction. A four-parameter power-law distribution presented in literature is proposed. Symmetric and asymmetric volume fraction profiles are presented. Suitable displacement functions that identically satisfy the boundary conditions at the simply supported edges are used to reduce the equilibrium equations to a set of coupled ordinary differential equations with variable coefficients, which are solved by generalized differential quadrature method, and natural frequency is obtained. The fast rate of convergence of the method is demonstrated, and to validate the results, comparisons are made with the available solutions for functionally graded isotropic shells with/without elastic foundations. The effect of the elastic foundation stiffness parameters and various geometrical parameters on the vibration behavior of the CGFR cylindrical panels is investigated. This work mainly contributes to illustrate the influence of the four parameters of power-law distributions on the vibration behavior of functionally graded orthotropic cylindrical panels resting on elastic foundation. This paper is also supposed to present useful results for continuous grading of matrix volume fraction in the thickness direction of a cylindrical panel on elastic foundation and comparison with similar discrete laminated composite cylindrical panel.     

Two-dimensional accessible solitons in PT-symmetric potentials

Two-dimensional parity-time (PT) symmetric potentials are introduced, which allow the existence of spatial solitons in the model of the strongly nonlocal nonlinear Schr?dinger equation. Two-dimensional accessible solitons are found in the form of solutions separating the radial amplitude, given in terms of Laguerre polynomials, a?phase function involving quadratic, linear, and constant phase shifts, and a specific azimuthal modulation function. Shape-preserving solitons are constructed from Laguerre?CGaussian functions containing the self-similar variable and an exponential form of the azimuthal modulation, containing sine and cosine functions, when a suitable PT-symmetric potential is chosen. Interesting soliton profiles and the corresponding PT-symmetric potentials are displayed for different values of the parameters.     

Asymptotic solutions for the Föppl – von Kármán equations governing deflections of thin axisymmetric annular plates

We are concerned with the deformation of thin, flat annular plates under a force applied orthogonally to the plane of the plate. This mechanical process can be described via a radial formulation of the Föppl – von Kármán equations, a set of nonlinear partial differential equations describing the deflections of thin flat plates. We are able to obtain analytical solutions for the radial Föppl – von Kármán equations with boundary conditions relevant for clamped, loosely clamped, and free inner and outer. This permits us to study the qualitative behavior of the out-of-plane deflections as well as the Airy stress function for a number of cases. Provided that an appropriate non-dimensionalization is taken, we find that the perturbation solutions are surprisingly valid for a wide variety of parameters, and compare favorably with numerical simulations in all cases (rather than just for small parameters). The results demonstrate that the ratio of the inner to outer radius of the annular plate will strongly influence the properties of the solutions, as will the specific boundary data considered. For instance, one may choose to fix the plate in place with a specific set of boundary conditions, in order to minimize the out-of-plane deflections. Other boundary conditions may result in undesirable behaviors.     

Closed-form solution of a shear deformable,extensional ring in contact between two rigid surfaces

Contact of a circular ring with a flat, rigid ground is considered using curved beam theory and analytical methods. Applications include tires, springs, and stiffeners, among others. The governing differential equations are derived using the principle of virtual work and the formulation includes deformations due to bending, transverse shear and circumferential extension. The three associated stiffness quantities, EI, GA and EA, respectively, remain as independent parameters in the differential equations. This allows the special cases such as an inextensible Timoshenko beam (EI and GA) or an extensible Euler beam (EI and EA) to be obtained directly by the appropriate limits. The effect of these three stiffness parameters on the contact pressure solution is studied, which shows how those fundamental parameters can be selected for the purpose of the application. Although the formulation is for small displacement theory, both radial and circumferential distributed loads are considered, which allows the pressure in the deformed state to be vertical rather than radial, which is shown to be important. Closed form expressions for all force and displacement quantities are obtained in terms of the angular location of the edge of contact, which must be determined numerically. Extensibility complicates the analytical expressions within the contact region, and a series solution is proposed in this case. A two-term asymptotic expression for the stiffness of the ring is determined analytically. Finally, all solutions are validated using the commercial finite element software ABAQUS, with attention to non-linear behavior and the range of validity of these solutions.     

Corrigendum of “Similarity analysis in magnetohydrodynamics: Hall effects on free convection flow and mass transfer past a semi-infinite vertical flat plate” [International Journal of Non-Linear Mechanics 38 (2003) 513–520]

The problem of steady, laminar mixed convection heat and mass transfer past a semi-infinite vertical plate in the presence of Hall current has been studied. The governing partial differential equations describing the problem are transformed to a system of non-linear ordinary differential equations with appropriate boundary conditions using Lie's method of infinitesimal transformation groups. The non-linear ordinary differential equations are solved numerically using Chebyshev spectral method. The effects of various parameters on the velocity profiles, temperature and concentration profiles are presented and discussed. This work is an extension and correction for the paper by Megahed et al. [1], published in International Journal of Non-Linear Mechanics 38 (2003) 513–520.     

层合圆柱三维温度场分析的半解析-精细积分法

针对圆柱体的三维温度场分析,提出了一种高效的半解析-精细积分法。将温度场展开为环向坐标的Fourier级数,并对径向坐标进行差分离散,从而把三维热传导方程简化为一系列二阶常微分方程;将这些二阶常微分方程转化为哈密顿体系下的一阶状态方程,并利用两点边值问题的精细积分法求解。由于该方法仅对径向坐标进行差分离散,故相对于传统的数值方法离散规模大幅度减少,不仅提高了计算效率、降低了存贮量,而且缓解了代数方程的病态问题。此外,针对Fourier半解析解,根据热平衡原理推导出了两种材料衔接面的半解析差分方程,从而为求解复合材料层合柱问题打下了基础。算例结果表明,即使对于细长比高达400的圆柱杆件,此方法仍然可以给出精度较高的解答。