Analytical test functions for free vibration analysis of rotating non-homogeneous Timoshenko beams

In this paper, the governing equations for free vibration of a non-homogeneous rotating Timoshenko beam, having uniform cross-section, is studied using an inverse problem approach, for both cantilever and pinned-free boundary conditions. The bending displacement and the rotation due to bending are assumed to be simple polynomials which satisfy all four boundary conditions. It is found that for certain polynomial variations of the material mass density, elastic modulus and shear modulus, along the length of the beam, the assumed polynomials serve as simple closed form solutions to the coupled second order governing differential equations with variable coefficients. It is found that there are an infinite number of analytical polynomial functions possible for material mass density, shear modulus and elastic modulus distributions, which share the same frequency and mode shape for a particular mode. The derived results are intended to serve as benchmark solutions for testing approximate or numerical methods used for the vibration analysis of rotating non-homogeneous Timoshenko beams.

Differential Games of Mixed Type with Control and Stopping Times

We study a zero sum differential game of mixed type where each player uses both control and stopping times. Under certain conditions we show that the value function for this problem exists and is the unique viscosity solution of the corresponding variational inequalities. We also show the existence of saddle point equilibrium for a special case of differential game.

A characteristic based numerical method with tracking for nonlinear wave equations

Many physical problems can be modeled by scalar, first-order, nonlinear, hyperbolic, partial differential equations (PDEs). The solutions to these PDEs often contain shock and rarefaction waves, where the solution becomes discontinuous or has a discontinuous derivative. One can encounter difficulties using traditional finite difference methods to solve these equations. In this paper, we introduce a numerical method for solving first-order scalar wave equations. The method involves solving ordinary differential equations (ODEs) to advance the solution along the characteristics and to propagate the characteristics in time. Shocks are created when characteristics cross, and the shocks are then propagated by applying analytical jump conditions. New characteristics are inserted in spreading rarefaction fans. New characteristics are also inserted when values on adjacent characteristics lie on opposite sides of an inflection point of a nonconvex flux function, Solutions along characteristics are propagated using a standard fourth-order Runge-Kutta ODE solver. Shocks waves are kept perfectly sharp. In addition, shock locations and velocities are determined without analyzing smeared profiles or taking numerical derivatives. In order to test the numerical method, we study analytically a particular class of nonlinear hyperbolic PDEs, deriving closed form solutions for certain special initial data. We also find bounded, smooth, self-similar solutions using group theoretic methods. The numerical method is validated against these analytical results. In addition, we compare the errors in our method with those using the Lax-Wendroff method for both convex and nonconvex flux functions. Finally, we apply the method to solve a PDE with a convex flux function describing the development of a thin liquid film on a horizontally rotating disk and a PDE with a nonconvex flux function, arising in a problem concerning flow in an underground reservoir.

Unsteady-Flow Of A Viscous-Fluid Between 2 Parallel Disks With A Time-Varying Gap Width And A Magnetic-Field

The unsteady incompressible viscous fluid flow between two parallel infinite disks which are located at a distance h(t*) at time t* has been studied. The upper disk moves towards the lower disk with velocity h'(t*). The lower disk is porous and rotates with angular velocity Omega(t*). A magnetic field B(t*) is applied perpendicular to the two disks. It has been found that the governing Navier-Stokes equations reduce to a set of ordinary differential equations if h(t*), a(t*) and B(t*) vary with time t* in a particular manner, i.e. h(t*) = H(1 - alpha t*)(1/2), Omega(t*) = Omega(0)(1 - alpha t*)(-1), B(t*) = B-0(1 - alpha t*)(-1/2). These ordinary differential equations have been solved numerically using a shooting method. For small Reynolds numbers, analytical solutions have been obtained using a regular perturbation technique. The effects of squeeze Reynolds numbers, Hartmann number and rotation of the disk on the flow pattern, normal force or load and torque have been studied in detail

Unsteady flow over a stretching surface with a magnetic field in a rotating fluid

The effect of the magnetic field on the unsteady flow over a stretching surface in a rotating fluid has been studied. The unsteadiness in the flow field is due to the time-dependent variation of the velocity of the stretching surface and the angular velocity of the rotating fluid. The Navier-Stokes equations and the energy equation governing the flow and the heat transfer admit a self-similar solution if the velocity of the stretching surface and the angular velocity of the rotating fluid vary inversely as a linear function of time. The resulting system of ordinary differential equations is solved numerically using a shooting method. The rotation parameter causes flow reversal in the component of the velocity parallel to the strerching surface and the magnetic field tends to prevent or delay the flow reversal. The surface shear stresses dong the stretching surface and in the rotating direction increase with the rotation parameter, but the surface heat transfer decreases. On the other hand, the magnetic field increases the surface shear stress along the stretching surface, but reduces the surface shear stress in the rotating direction and the surface heat transfer. The effect of the unsteady parameter is more pronounced on the velocity profiles in the rotating direction and temperature profiles.

Precipitation In Small Systems--II. Mean Field Equations More Effective Than Population Balance

Part I (Manjunath et al., 1994, Chem. Engng Sci. 49, 1451-1463) of this paper showed that the random particle numbers and size distributions in precipitation processes in very small drops obtained by stochastic simulation techniques deviate substantially from the predictions of conventional population balance. The foregoing problem is considered in this paper in terms of a mean field approximation obtained by applying a first-order closure to an unclosed set of mean field equations presented in Part I. The mean field approximation consists of two mutually coupled partial differential equations featuring (i) the probability distribution for residual supersaturation and (ii) the mean number density of particles for each size and supersaturation from which all average properties and fluctuations can be calculated. The mean field equations have been solved by finite difference methods for (i) crystallization and (ii) precipitation of a metal hydroxide both occurring in a single drop of specified initial supersaturation. The results for the average number of particles, average residual supersaturation, the average size distribution, and fluctuations about the average values have been compared with those obtained by stochastic simulation techniques and by population balance. This comparison shows that the mean field predictions are substantially superior to those of population balance as judged by the close proximity of results from the former to those from stochastic simulations. The agreement is excellent for broad initial supersaturations at short times but deteriorates progressively at larger times. For steep initial supersaturation distributions, predictions of the mean field theory are not satisfactory thus calling for higher-order approximations. The merit of the mean field approximation over stochastic simulation lies in its potential to reduce expensive computation times involved in simulation. More effective computational techniques could not only enhance this advantage of the mean field approximation but also make it possible to use higher-order approximations eliminating the constraints under which the stochastic dynamics of the process can be predicted accurately.

Unsteady Natural Convection Flow in a Square Cavity Filled with a Porous Medium Due to Impulsive Change in Wall Temperature

Unsteady natural convection flow in a two- dimensional square cavity filled with a porous material has been studied. The flow is initially steady where the left- hand vertical wall has temperature T-h and the right- hand vertical wall is maintained at temperature T-c ( T-h > T-c) and the horizontal walls are insulated. At time t > 0, the left- hand vertical wall temperature is suddenly raised to (T-h) over bar ((T-h) over bar > T-h) which introduces unsteadiness in the flow field. The partial differential equations governing the unsteady natural convection flow have been solved numerically using a finite control volume method. The computation has been carried out until the final steady state is reached. It is found that the average Nusselt number attains a minimum during the transient period and that the time required to reach the final steady state is longer for low Rayleigh number and shorter for high Rayleigh number.

Stability Of Large Damped Flexible Spacecraft With Stored Angular Momentum

Vibrational stability of large flexible structurally damped spacecraft carrying internal angular momentum and undergoing large rigid body rotations is analysed modeling the systems as elastic continua. Initially, analytical solutions to the motion of rigid gyrostats under torque-free conditions are developed. The solutions to the gyrostats modeled as axisymmetric and triaxial spacecraft carrying three and two constant speed momentum wheels, respectively, with spin axes aligned with body principal axes are shown to be complicated. These represent extensions of solutions for simpler cases existing in the literature. Using these solutions and modal analysis, the vibrational equations are reduced to linear ordinary differential equations. Equations with periodically varying coefficients are analysed applying Floquet theory. Study of a few typical beam- and plate-like spacecraft configurations indicate that the introduction of a single reaction wheel into an axisymmetric satellite does not alter the stability criterion. However, introduction of constant speed rotors deteriorates vibrational stability. Effects of structural damping and vehicle inertia ratio are also studied.

Semi-similar solution of an unsteady compressible three-dimensional stagnation point boundary layer flow with massive blowing

A semi-similar solution of an unsteady laminar compressible three-dimensional stagnation point boundary layer flow with massive blowing has been obtained when the free stream velocity varies arbitrarily with time. The resulting partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme with a quasi-linearization technique in the nodal point region and an implicit finite-difference scheme with a parametric differentiation technique in the saddle point region. The results have been obtained for two particular unsteady free stream velocity distributions: (i) an accelerating stream and (ii) a fluctuating stream. Results show that the skin-friction and heat-transfer parameters respond significantly to the time dependent arbitrary free stream velocity. Velocity and enthalpy profiles approach their free stream values faster as time increases. There is a reverse flow in the y-wise velocity profile, and overshoot in the x-wise velocity and enthalpy profiles in the saddle point region, which increase as injection and wall temperature increase. Location of the dividing streamline increases as injection increases, but as the wall temperature and time increase, it decreases.

Unsteady compressible 3-dimensional boundary-layer flow near an asymmetric stagnation point with mass transfer

The heat and mass transfer for unsteady laminar compressible boundary-layer flow, which is asymmetric with respect to a 3-dimensional stagnation point (i.e. for a jet incident at an angle on the body), have been studied. It is assumed that the free-stream velocity, wall temperature, and surface mass transfer vary arbitrarily with time and also that the gas has variable properties. The solution in the neighbourhood of the stagnation point has been obtained by series expansion in the longitudinal distance. The resulting partial differential equations have been solved numerically using an implicit finite-difference scheme. The results show that, in contrast with the symmetric flow, the maximum heat transfer does not occur at the stagnation point. The skin-friction and heat-transfer components due to asymmetric flow are only weakly affected by the mass transfer as compared to those components associated with symmetric flow. The variation of the wall temperature with time has a strong effect on the heat transfer component associated with the symmetric part of the flow. The skin friction and heat transfer are strongly affected by the variation of the density-viscosity product across the boundary layer. The skin friction responds more to the fluctuations of the free stream oscillating velocities than the heat transfer. The results have been compared with the available results and they are found to be in excellent agreement.

Incompressible boundary layer for longitudinal flow over a cylinder with an applied magnetic field

The flow, heat and mass transfer problem for a steady laminar incompressible boundary layer flow in an electrically conducting fluid over a longitudinal cylinder with an applied magnetic field has been studied. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. The results are found to be strongly dependent on the magnetic field and dissipation parameter. The effect of the mass transfer is more pronounced on the skin friction than on the heat transfer. The results have been compared with those of the series solution, the asymptotic solution, the Glauert and Lighthill's solution, local similarity, local nonsimilarity and difference-differential methods. Good agreement is found with all of them, except with the results of the local similarity and series solution methods.

On a linear steady-state system of equations in microcontinuum fluid mechanics

Galerkin representations and integral representations are obtained for the linearized system of coupled differential equations governing steady incompressible flow of a micropolar fluid. The special case of 2-dimensional Stokes flows is then examined and further representation formulae as well as asymptotic expressions, are generated for both the microrotation and velocity vectors. With the aid of these formulae, the Stokes Paradox for micropolar fluids is established.

Mass transfer effects on the generalised vortex flow over a stationary surface with or without magnetic field

The effect of injection and suction on the generalised vortex flow of a steady laminar incompressible fluid over a stationary infinite disc with or without magnetic field under boundary-layer approximations has been studied. The coupled nonlinear ordinary differential equations governing the self-similar flow have been numerically solved using the finite-difference scheme. The results indicate that the injection produces a deeper inflow layer and de-stabilises the motion while suction or magnetic field suppresses the inflow layer and produces stability. The effect of decreasingn, the parameter characterising the nature of vortex flow, is similar to that of increasing the injection rate.

Generalized Burgers equations and Euler-Painlev transcendents. I

Initial-value problems for the generalized Burgers equation (GBE) ut+u betaux+lambdaualpha =(delta/2)uxx are discussed for the single hump type of initial data both continuous and discontinuous. The numerical solution is carried to the self-similar ``intermediate asymptotic'' regime when the solution is given analytically by the self-similar form. The nonlinear (transformed) ordinary differential equations (ODE's) describing the self-similar form are generalizations of a class discussed by Euler and Painlevé and quoted by Kamke. These ODE's are new, and it is postulated that they characterize GBE's in the same manner as the Painlev equations categorize the Kortweg-de Vries (KdV) type. A connection problem for some related ODE's satisfying proper asymptotic conditions at x=±[infinity], is solved. The range of amplitude parameter is found for which the solution of the connection problem exists. The other solutions of the above GBE, which display several interesting features such as peaking, breaking, and a long shelf on the left for negative values of the damping coefficient lambda, are also discussed. The results are compared with those holding for the modified KdV equation with damping. Journal of Mathematical Physics is copyrighted by The American Institute of Physics.

Nonsimilar incompressible laminar boundary layers with magnetic field.

The solution of the steady laminar incompressible nonsimilar magneto-hydrodynamic boundary layer flow and heat transfer problem with viscous dissipation for electrically conducting fluids over two-dimensional and axisymmetric bodies with pressure gradient and magnetic field has been presented. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. The computations have been carried out for flow over a cylinder and a sphere. The results indicate that the magnetic field tends to delay or prevent separation. The heat transfer strongly depends on the viscous dissipation parameter. When the dissipation parameter is positive (i.e. when the temperature of the wall is greater than the freestream temperature) and exceeds a certain value, the hot wall ceases to be cooled by the stream of cooler air because the ‘heat cushion’ provided by the frictional heat prevents cooling whereas the effect of the magnetic field is to remove the ‘heat cushion’ so that the wall continues to be cooled. The results are found to be in good agreement with those of the local similarity and local nonsimilarity methods except near the point of separation, but they are in excellent agreement with those of the difference-differential technique even near the point of separation.