In-plane bending of a beam resting on a rigid rough foundation

Summary The bending of a finite-length beam that lies on a rigid, rough, flat foundation and interacts with it in accordance to the dry friction law is considered. Loading by bending moments applied at the ends of the beam is studied in detail. The problem is found to be a self-similar one. For small moments, the central part of the beam remains undeflected, and the problem reduces to the solution of an infinite system of algebraic equations. Large moments deflect the entire length of the beam, and the problem partly loses its self-similarity. In this case, the problem reduces to the solution of a successively decreasing number of ordinary differential equations along with some algebraical equations. The solution for the latter case provides initial conditions for the former one. This permits to obtain a solution for any value of the moment. Received 5 November 1996; accepted for publication 27 January 1997     

Friction and heat transfer in an incompressible fluid near a plane stagnation point

In the neighborhood of a plane stagnation point, the flow and heat transfer of an incompressible fluid are studied. In the inner flow region, the velocity and pressure fields are described by the complete Navier-Stokes equations, and the temperature field is described by the complete energy equation. In the outer flow region, a two-term asymptotic solution of the corresponding equations is obtained. The problem is reduced to the numerical solution of ordinary differential equations. Numerical results are discussed.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 52–65, July–August, 1996.     

Diffusional growth of cloud particles: existence and uniqueness of solutions

Diffusional growth of cloud particles is commonly described by a coupled system of parabolic equations and ordinary differential equations. The Dirichlet boundary condition for the parabolic equation is obtained from the solution of the ordinary differential equations, but this solution itself depends on the solution of the parabolic equations. We first present the governing equations describing diffusional growth of cloud particles. In a second step, we consider a simplified model problem, motivated by the diffusional growth equations. The main difference between the simplified model problem and the diffusional growth equations consists in neglecting the dependence of the domain for the parabolic equations on the solution. For the model problem, we show unique solvability using a fixed point method. Finally, we discuss application of the main result for the model problem to the diffusional growth equations and illustrate these equations with the help of a numerical solution.     

The effect of radiation on transient natural convection past a doubly infinite plate

The present research note is concerned with the transient (short time) simultaneous free convection and radiation analysis of a viscous fluid along a doubly infinite vertical isothermal flat plate. To simplify a very complicated problem, an incompressible flow field is used in the analysis. Generally, the exact numerical solution of this problem is quite lengthy. However, by considering an optically thick radiating gas, expressed by the Rosseland diffusion approximation, the solution is much simpler. Moreover, this case leads to a complete similarity transformation of the governing partial differential equations into a set of ordinary differential equations. An exact numerical solution is obtained of the resulting ordinary differential equations for a Prandtl number equal to 0.733 and for a wide range of involved parameters.     

扇形截面杆扭转问题的差分线法

导出了扇形截面杆扭转问题偏微分方程的差分线法常微分方程组, 并解析求解了该方程组, 得到了扭转应力函数的半解析解, 计算了扭转应力及扭转刚度. 计算过程中, 用追赶法计算 常微分方程组的特解, 用公式计算三对角矩阵的特征值与特征向量, 利用实对阵矩阵的特征 向量相互正交的特性避免矩阵求逆计算, 利用复化梯形公式计算扭转刚度. 整个求解过程在 角度方向离散微分方程和用复化梯形公式进行面积积分时引入了误差, 其他求解过程是精确 的. 计算结果与已有结果进行了对比, 显示了算法的正确性. 该算法对工程中扇形截面扭 转杆的设计有一定的实用价值.     

Quasi-similar solutions for blast waves with internal heat transfer effects

The point explosion problem with internal heat transfer effects taken into account is analysed. The classical inviscid solutions to this problem yield a non-physical phenomenon of infinite temperature and zero density at the center of explosion for all times. With heat transfer fluxes considered, the solution near the center of symmetry is improved and finite values for temperature are obtained. The non-self-similar solution of the problem is based on the quasi-similar approximate technique which reduces the non-linear partial differential equations governing the problem to ordinary differential ones. However, this formulation yields a two-point boundary-value problem. To facilitate the integration, the flow field is first divided into two regions: an outer inviscid region and an inner region where dissipation effects are manifested. This results in two sets of ordinary differential equations expressing the conservation equations in the inner and outer regions which are then solved and matched together to yield the composite solution. Secondly, the problem is then transformed into an initial-value one. Using the results of the composite solution, the governing equations can be integrated directly from the center until the shock front. The structure of the non-self-similar flow fields with internal heat transfer effects is then fully determined for specified values of the heat transfer parameters.     

Jet with variable fluid properties: Free jet and dissipative jet

An analysis of the laminar jet of an incompressible Newtonian fluid emerging from a narrow slot or a circular hole, where the physical properties like viscosity and thermal conductivity depends upon the temperature, is given. Both the cases: the case of In the absence of viscous heat dissipation and the case of In the presence of viscous heat dissipation are considered. The governing partial differential equations of the flow problem are transformed into the ordinary differential equations by group theoretic technique. The Runge–Kutta method is applied to obtained numerical solution of the transformed ordinary differential equations.     

METHOD–OF–LINES SOLUTION OF TIME–DEPENDENT TWO–DIMENSIONAL NAVIER–STOKES EQUATIONS

A novel approach to the development of a code for the solution of the time-dependent two-dimensional Navier–Stokes equations is described. The code involves coupling between the method of lines (MOL) for the solution of partial differential equations and a parabolic algorithm which removes the necessity of iterative solution on pressure and solution of a Poisson-type equation for the pressure. The code is applied to a test problem involving the solution of transient laminar flow in a short pipe for an incompressible Newtonian fluid. Comparisons show that the MOL solutions are in good agreement with the previously reported values. The proposed method described in this paper demonstrates the ease with which the Navier–Stokes equations can be solved in an accurate manner using sophisticated numerical algorithms for the solution of ordinary differential equations (ODEs).     

Hydromagnetic flow and heat transfer adjacent to a stretching vertical sheet with prescribed surface heat flux

The similarity solution for the problem of mixed convection boundary layer flow adjacent to a stretching vertical sheet in an incompressible electrically conducting fluid in the presence of a transverse magnetic field is presented. It is assumed that the sheet is stretched with a power-law velocity and is subjected to a variable surface heat flux. The governing partial differential equations are first transformed into a system of non-linear ordinary differential equations, before being solved numerically by the Keller-box method. The numerical results obtained are then compared with previously reported cases available in the literature as well as the series solution for certain values of parameters, to support their validity. The effects of the governing parameters on the flow field and heat transfer characteristics are obtained and discussed.     

Well posed formulations of holonomic mechanical system dynamics and design sensitivity analysis

Abstract

Manifold theoretic ordinary differential equations of motion for holonomic mechanical systems that depend on problem data, or design variables, are shown to be well posed; i.e., they have a unique solution that depends continuously on problem data. It is proved that these differential equations are equivalent to the d’Alembert variational formulation and the index 3 Lagrange multiplier formulation of differential-algebraic equations of motion, which are also shown to be well posed. These results provide a foundation for dynamic system design sensitivity analysis, which requires differentiability of solutions of the equations of motion with respect to design variables.     

Chemically reactive hydromagnetic flow of a second-grade fluid in a semi-porous channel

An analysis of a second-grade fluid in a semi-porous channel in the presence of a chemical reaction is carried out to study the effects of mass transfer and magnetohydrodynamics. The upper wall of the channel is porous, while the lower wall is impermeable. The basic governing flow equations are transformed into a set of nonlinear ordinary differential equations by means of a similarity transformation. An approximate analytical solution of nonlinear differential equations is constructed by using the homotopy analysis method. The features of the flow and concentration fields are analyzed for various problem parameters. Numerical values of the skin friction coefficient and the rate of mass transfer at the wall are found.     

Homotopy analysis and differential quadrature solution of the problem of free-convective magnetohydrodynamic flow over a stretching sheet with the Hall effect and mass transfer taken into account

This paper presents an analytical solution of the problem of free-convective magnetohydrodynamic flow over a stretched sheet with the Hall effect and mass transfer taken into account. A similarity transform reduces the Navier-Stokes, energy, Ohm law, and mass-transfer equations to a system of nonlinear ordinary differential equations. The governing equations are solved analytically using an analytical method for solving nonlinear problems, namely, the homotopy analysis method. The results are compared with the results of a promising numerical method of differential quadrature developed by the authors. It is shown that there is very good agreement between analytical results and those obtained by the differential quadrature method. The differential quadrature method was validated, and the effects of non-dimensional parameters on the velocity, temperature and concentration profiles were studied.     

Analytical solution of the contact problem for a system of bodies under collective wear

The contact problem is considered for a system of bodies subject to wear on a common base. The deformation properties of the bodies and the base are described by the Winkler model. The problem is reduced to a system of ordinary differential equations and an integral equation of hereditary type with difference kernel. The solution of the problem is constructed with the use of the Laplace transform. The asymptotics of the solution at large times is studied. The continuity conditions for the contact of each of the bodies with the base are derived.     

Motion of a strong shock front in a two-dimensional gas layer with moving internal boundary

We examine the problem of planar one-dimensional motion of a strong shock wave with moving internal boundary in which the initial position of the front, its intensity, the mass of the gas involved in the motion, and the energy contained in this gas are known. The problem is not self-similar and its exact solution, which involves working with partial differential equations, presents serious difficulties. In the following we determine the law of shock-front motion in this problem via the method of [1], which makes it possible to find a system of ordinary differential equations for the problem. The method is based on an initial specification of the power-law coupling between the dimensionless Lagrangian and Eulerian variables and replacement of the energy equation by this coupling and the energy integral. The solution is sought in the first approximation.     

On heat and mass transfer in the unsteady squeezing flow between parallel plates

This paper reports the heat and mass transfer characteristics in a viscous fluid which is squeezed between parallel plates. The governing partial differential equations for unsteady two-dimensional flow with heat and mass transfer of a viscous fluid are reduced to ordinary differential equations by similarity transformations. Homotopy analysis method (HAM) is employed to construct the series solution of the problem. Physical interpretation to various embedding parameters is assigned through graphs for temperature and concentration profiles and tables for skin friction coefficient, local Nusselt number and local Sherwood number.     

Heat transfer in a radiating, nonsteady, three-dimensional stagnation flow

A similar solution has been obtained to the problem of simultaneous radiation and convection for nonsteady stagnation point flow over a three-dimensional blunt body with both boundary layer suction and injection. The diffusion approximation is used to characterize the radiative heat flux. The three-dimensional, time-dependent equations of motion and the energy equation have been transformed into a set of ordinary differential equations by the similarity transformation and the resulting ordinary differential equations have been solved numerically. The effects of accelerating and decelerating flow, the three-dimensional geometry, injection and suction, hot and cold wall conditions, and the conduction-to-radiation parameter on the temperature distribution within the flow have been investigated.     

On the cone and plate deformation of a rubber-like material: Instability

In this report it is shown that the place boundary-value problem, for a small deformation superimposed on the large cone-and-plate deformation of a Mooney-Rivlin material, has no unique solution at some displacement angles of the cone. The kinematics of the small deformation are chosen such that the problem is reduced to a set of ordinary differential equations. With respect to this restricted class of kinematics the cone-and-plate deformation is unstable.     

The use of Sinc‐collocation method for solving Falkner–Skan boundary‐layer equation

The MHD Falkner–Skan equation arises in the study of laminar boundary layers exhibiting similarity on the semi‐infinite domain. The proposed approach is equipped by the orthogonal Sinc functions that have perfect properties. This method solves the problem on the semi‐infinite domain without truncating it to a finite domain and transforming domain of the problem to a finite domain. In addition, the governing partial differential equations are transformed into a system of ordinary differential equations using similarity variables, and then they are solved numerically by the Sinc‐collocation method. It is shown that the Sinc‐collocation method converges to the solution at an exponential rate. Copyright © 2010 John Wiley & Sons, Ltd.     

Bending analyses of 1D orthorhombic quasicrystal plates

The meshless Petrov–Galerkin method (MLPG) is applied to plate bending analysis in 1D orthorhombic quasicrystals (QCs) under static and transient dynamic loads. The Bak and elasto-hydrodynamic models are applied for phason governing equation in the elastodynamic case. The phason displacement for the orthorhombic QC in the first-order shear deformation plate theory depends only on the in-plane coordinates on the mean plate surface. Nodal points are randomly distributed over the mean surface of the considered plate. Each node is the center of a circle surrounding this node. The coupled governing partial differential equations are satisfied in a weak-form on small fictitious subdomains. The spatial variations of the phonon and phason displacements are approximated by the moving least-squares (MLS) scheme. After performing the spatial MLS approximation, a system of ordinary differential equations (ODEs) for nodal unknowns is obtained. The system of the ODEs of the second order is solved by the Houbolt finite-difference scheme. Our numerical examples demonstrate clearly the effect of the coupling parameter on both static and dynamic phonon/phason deflections.     

Exact piezothermoelastic solution of simply-supported orthotropic circular cylindrical panel in cylindrical bending

Summary This work presents an exact piezothermoelastic solution of infinitely long, simply supported, cylindrically orthotropic, piezoelectric, radially polarised, circular cylindrical shell panel in cylindrical bending under thermal and electrostatic excitation. The general solution of the governing differential equations is obtained by separation of variables. The displacements, electric potential and temperature are expanded in appropriate Fourier series in the circumferential coordinate to satisfy the boundary conditions at the simply-supported longitudinal edges. The governing equations reduce to Euler-Cauchy type of ordinary differential equations. Their general solution involves six constants for each Fourier component. These are solved from the algebraic equations obtained by satisfying the boundary conditions at the lateral surfaces. The solution of the inverse problem of inferring the applied temperature field from the given measured distribution of electrical potential difference between the lateral surfaces of the shell has also been presented. Numerical results are presented for typical thermal and electrostatic loadings for various values of radius to thickness ratio.