On the uniqueness of large deflections of a uniform cantilever beam under a tip-concentrated rotational load

The problem of a uniform cantilever beam under a tip-concentrated load, which rotates in relation with the tip-rotation of the beam, is studied in this paper. The formulation of the problem results in non-linear ordinary differential equations amenable to numerical integration. A relation is obtained for the applied tip-concentrated load in terms of the tip-angle of the beam. When the tip-concentrated load acts always normal to the undeformed axis of the beam (the rotation parameter, β=0) there is a possibility of obtaining non-unique solution for the applied load. This phenomenon is also observed for other rotation parameters less than unity. When the tip-concentrated load is acting normal to the deformed axis of the beam (β=1), many load parameters are obtained for a tip-angle with different deformed configurations of the beam. However, each load parameter corresponds to a tip-angle, which confirms the uniqueness on the solution of non-linear differential equations.     

Solution of a mixed dynamic problem on the displacements given on the boundary of a half-plane lying on the surface of an isotropic homogeneous elastic half-space

We consider the problem on the motion of an isotropic elastic body occupying the half-space z ≥ 0 on whose boundary, along the half-plane x ≥ 0, the horizontal components of displacement are given, while the remaining part of the boundary is stress-free. We seek the solution by the method of integral Laplace transforms with respect to time t and Fourier transforms with respect to the coordinates x, y; the problem is reduced to a system of Wiener-Hopf equations, which can be solved by the methods of singular-integral equations and circulants. We invert the integral transforms and reduce the solution to the Smirnov-Sobolev form. We calculate the tangential stress intensity coefficients near the boundary z = 0, x = 0, |y| < ∞ of the half-plane. The circulant method for solving the Wiener-Hopf system was proposed in [1]. A static problem similar to that considered in the present paper was solved earlier. The Hilbert problem was reduced to a system of Fredholm integral equations in [2]. In the present paper, we solve the above problem by reducing the solution to quadratures and a quasiregular system of Fredholm integral equations. We give a numerical solution of the Fredholm equations and calculate the integrals for the tangential stress intensity coefficients.     

Stability of dual solutions in stagnation-point flow and heat transfer over a porous shrinking sheet with thermal radiation

An analysis is carried out to study the steady two-dimensional stagnation-point flow and heat transfer of an incompressible viscous fluid over a porous shrinking sheet in the presence of thermal radiation. A set of similarity transformations reduce the boundary layer equations to a set of non-linear ordinary differential equations which are solved numerically using fourth order Runge-Kutta method with shooting technique. The analysis of the result obtained shows that as the porosity parameter β increases, the range of region of existence of similarity solution increases. It is also observed that multiple solutions exist for a certain range of the ratio of the shrinking velocity to the free stream velocity (i.e., α) which again depends on β. We then discuss the stability of the unsteady solutions about each steady solution, showing that one steady state solution corresponds to a stable solution whereas the other corresponds to an unstable solution. The stable solution corresponds to the physically relevant solution. Further we obtain numerical results for each solution, which enable us to discuss the features of the respective solutions.     

Application of a non-linear local analysis method for the problem of mixed convection instability

We consider the problem of laminar mixed convection flow between parallel, vertical and uniformly heated plates where the governing dimensionless parameters are the Prandtl, Rayleigh and Reynolds numbers. Using the method based on the centre manifold theorem which was derived from the general theory of dynamical systems, we reduce a three-dimensional simplified model of ordinary differential amplitude equations emanating from the original Navier-Stokes system of the problem in the vicinity of a trivial stationary solution. We have found that when the forcing parameter, the Rayleigh number, increases beyond the critical value Ra_s, the stationary solution is a pitchfork bifurcation point of the system.     

New exact solution of Euler’s equations (rigid body dynamics) in the case of rotation over the fixed point

A new exact solution of Euler’s equations (rigid body dynamics) is presented here. All the components of angular velocity of rigid body for such a solution differ from both the cases of symmetric rigid rotor (which has two equal moments of inertia: Lagrange’s or Kovalevskaya’s case), and from the Euler’s case when all the applied torques are zero, or from other well-known particular cases. The key features are the next: the center of mass of rigid body is assumed to be located at meridional plane along the main principal axis of inertia of rigid body, besides, the principal moments of inertia are assumed to satisfy to a simple algebraic equality. Also, there is a restriction at choosing of initial conditions. Such a solution is also proved to satisfy to Euler–Poinsot equations, including invariants of motion and additional Euler’s invariant (square of the vector of angular momentum is a constant). So, such a solution is a generalization of Euler’s case.     

A fractional differential constitutive model for dynamic stress intensity factors of an anti-plane crack in viscoelastic materials

Fractional differential constitutive relationships are introduced to depict the history of dynamic stress inten- sity factors （DSIFs） for a semi-infinite crack in infinite viscoelastic material subjected to anti-plane shear impact load. The basic equations which govern the anti-plane deformation behavior are converted to a fractional wave-like equation. By utilizing Laplace and Fourier integral transforms, the fractional wave-like equation is cast into an ordinary differential equation （ODE）. The unknown function in the solution of ODE is obtained by applying Fourier transform directly to the boundary conditions of fractional wave-like equation in Laplace domain instead of solving dual integral equations. Analytical solutions of DSIFs in Laplace domain are derived by Wiener-Hopf technique and the numerical solutions of DSIFs in time domain are obtained by Talbot algorithm. The effects of four parameters α, β, b1, b2 of the fractional dif- ferential constitutive model on DSIFs are discussed. The numerical results show that the present fractional differential constitutive model can well describe the behavior of DSIFs of anti-plane fracture in viscoelastic materials, and the model is also compatible with solutions of DSIFs of anti-plane fracture in elastic materials.     

Stochastic stability of a viscoelastic column axially loaded by a white noise force

This paper deals with the analysis of stability of a hinged-hinged viscoelastic column subjected to a non-zero mean stochastic axial force. The randomly variable part of this is described by a stationary Gaussian white noise process. The viscosity affects the curvature of the column, for which the classic Euler-Bernoulli's model is adopted. The viscosity is described by the linear Kelvin-Voigt's model. A dynamic stability analysis is performed. Normal modes are introduced in the integro-differential equation of motion so that uncoupled modal equations are retrieved. With reference to the first mode, by using an additional state variable, three Itô’s ODE are obtained, from which the differential equations ruling the response statistical moment evolution are written by means of Itô’s differential rule. The zero solution, that is undeformed straight column, corresponds to zero moments. If the column is perturbed, it is stable when the response moments tend to zero. A necessary and sufficient condition of stability in the moments of order r is that the matrix A_r of the coefficients of the ODE system ruling them has negative real eigenvalues and complex eigenvalues with negative real parts. Because of the linearity of the system the stability of the first two moments is the strongest condition of stability. If the mean axial force μ_P or the white noise intensity wP are increased, there exist critical values μ_Pcr, wPcr for which almost an eigenvalue is positive. The critical mean axial force is found to be inversely proportional to the parameter φ_∞, which measures the amount of viscous deformation. The search for the critical values of wP is made numerically, and several graphs are presented for a representative column.     

Heat transfer analysis of a particle-containing channel flow

This paper describes an analytical model of heat transfer in a two-dimensional, steady, nonreacting particle-containing channel flow. An idealized gas flow of specified uniform velocity between insulated parallel plates is assumed and the nonvaporizing particles are conceptualized as contained within an thin sheet injected at the symmetry plane. Two dimensionless parameters that affect the solution are described. These are the effective gas diffusivityK and the dimensionless particle number densityP. The linear, coupled differential equations governing the energy exchange between the gas and liquid phases are solved by means of the Green's function technique. This procedure yields a Volterra integral-series equation as the solution of the gas-phase energy equation. A series solution of this integral equation is obtained by the method of successive substitutions and terms up to second order are calculated.     

MHD flow of a third-grade fluid due to eccentric rotations of a porous disk and a fluid at infinity

The problem of magnetohydrodynamics (MHD) flow of a conducting, incompressible third-grade fluid due to non-coaxial rotations of a porous disk and a fluid at infinity in the presence of a uniform transverse magnetic field is considered. An exact analysis is carried out to model the governing non-linear partial differential equation. A numerical solution of the third-order non-linear partial differential equation has been obtained. Several graphs and tables have been drawn to show the influence of porosity ε, magnetic parameter N, material parameters α and β on the velocity distribution.     

The state vector methods for space axisymmetric problems in multilayered piezoelectric media

The state vector equations for space axisymmetric problems of transversely isotropic piezoelectric media are established from the basic equations. Using the Hankel transform, the state vector equations are reduced to a system of ordinary differential equations. An analytical solution of the problems in the Hankel transform space is presented in the form of the product of initial state vector and transfer matrix. The transfer matrices are given for the three distinct eigenvalues. Applications of the solutions are discussed. An analytical solution for the transversely isotropic semi-infinite piezoelectric media subjected to concerted point loads on the surface z=0 is presented in the Hankel transform space. Using transfer matrix and the continuity conditions at the layer interfaces, the general solution formulation of N-layered transversely isotropic piezoelectric media is given. A selected set of numerical solutions is presented for a layered semi-infinite piezoelectric solid.     

Approximate fundamental solutions and wave fronts for general anisotropic materials

The time-dependent differential equations of elastodynamics for homogeneous solids with a general structure of anisotropy are considered in the paper. A new method of computation of the fundamental solution for these equations is proposed. This method consists of the following. Applying the Fourier transformation with respect to space variables to these equations, we obtain a system of second order ordinary differential equations whose coefficients depend on Fourier parameters. Using the matrix transformations and properties of the coefficients, the Fourier image of the fundamental solution is computed. Finally, the fundamental solution is calculated by the inverse Fourier transformation to the obtained Fourier transform. The implementation and justification of the suggested method have been made by computational experiments in MATLAB. These experiments confirm the robustness of the suggested method. The visualization of the displacement components in general homogeneous anisotropic solids by modern computer tools allows us to see and evaluate the dependence between the structure of solids and the behavior of the displacement field. Our method allows users to observe the elastic wave propagation, arising from pulse point forces of the form e^mδ(x)δ(t), in monoclinic, triclinic and other anisotropic solids. The visualization of displacement components gives knowledge about the form of fronts of elastic wave propagation in Sodium Thiosulfate with monoclinic and Copper Sulphate Pentahydrate with triclinic structures of anisotropy.     

Steady Motions of Viscoelastic Fluids in Three‐Dimensional Exterior Domains. Existence, Uniqueness and Asymptotic Behaviour

. We study systems of nonlinear partial differential equations governing the steady motion of certain viscoelastic non‐Newtonian fluids around a rigid body ${\cal B}\subset\real^{3}$ . Considering the equations in a suitably decomposed form, we obtain, for small and sufficiently regular data, existence of a unique solution using a fixed‐point argument in an appropriate functional setting. This setting contains also the asymptotic decay of the solution. Our model equations include the second‐grade and the Maxwell fluid.     

Multiple spatial and temporal scales method for numerical simulation of non-classical heat conduction problems: one dimensional case

A multiple spatial and temporal scales method is studied to simulate numerically the phenomenon of non-Fourier heat conduction in periodic heterogeneous materials. The model developed is based on the higher-order homogenization theory with multiple spatial and temporal scales in one dimensional case. The amplified spatial scale and the reduced temporal scale are introduced respectively to account for the fluctuations of non-Fourier heat conduction due to material heterogeneity and non-local effect of the homogenized solution. By separating the governing equations into various scales, the different orders of homogenized non-Fourier heat conduction equations are obtained. The reduced time dependence is thus eliminated and the fourth-order governing differential equations are derived. To avoid the necessity of C¹ continuous finite element implementation, a C⁰ continuous mixed finite element approximation scheme is put forward. Numerical results are shown to demonstrate the efficiency and validity of the proposed method.     

Analysis of a quasistatic contact problem with adhesion and nonlocal friction for viscoelastic materials

A mathematical model is established to describe a contact problem between a deformable body and a foundation. The contact is bilateral and modelled with a nonlocal friction law, in which adhesion is taken into account. Evolution of the bonding field is described by a first-order differential equation. The materials behavior is modelled with a nonlinear viscoelastic constitutive law. A variational formulation of the mechanical problem is derived, and the existence and uniqueness of the weak solution can be proven if the coefficient of friction is sufficiently small. The proof is based on arguments of time-dependent variational inequalities, differential equations, and the Banach fixed-point theorem.     

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.     

Soret and Dufour effects on Double-Diffusive Free Convection From a Corrugated Vertical Surface in a Non-Darcy Porous Medium

Combined heat and mass transfer process by natural convection from a wavy vertical surface immersed in a fluid-saturated semi-infinite porous medium due to Soret and Dufour effects for Forchheimer extended non-Darcy model has been analyzed. A similarity transformation followed by a wavy to flat surface transformation is applied to the governing coupled non-linear partial differential equations, and they are reduced to boundary layer equations. The obtained boundary layer equations are solved by finite difference scheme based on the Keller-Box approach in conjunction with block-tridiagonal solver. Detailed simulations are carried out for a wide range of parameters like Groshof number (Gr*), Lewis number (Le), Buoyancy ratio (B), Wavy wall amplitude (a), Soret number (S _r ), and Dufour number (D _f ). Comparison tables local and average Nusselt (Nu) number, local and average Sherwood (Sh) number plots are presented.     

Free convection from a horizontal moving sheet

In this paper the free convection flow through a thin rigid hot sheet moving horizontally out of a slot is considered. It is found that there is a similarity formulation of the boundary-layer equations so that the problem reduces to solving a system of coupled ordinary differential equations with suitable boundary conditions. This system of equations is solved numerically for various values of the Prandtl number,Pr, namely 0.45≤Pr≤10000. It is found that for the flow under the sheet there is a reverse flow region near the sheet for small values ofPr, whilst in the case of the flow above the sheet there is no reverse flow region for any value ofPr we have investigated. For the flow under the sheet an asymptotic behaviour, which is valid near the minimum value of the Prandtl number for which it is possible to obtain a numerical solution, is proposed.     

Analysis of a benchmark solution for non-Newtonian radial displacement in porous media

We present an analytical formulation useful to interpret the key phenomena involved in non-Newtonian displacement in porous media and an analysis of the results obtained by considering the uncertainty associated with relevant problem parameters. To derive a benchmark solution, we consider the radial dynamics of a moving stable interface in a porous domain saturated by two fluids, displacing and displaced, both non-Newtonian of shear-thinning power-law behavior, assuming the pressure and velocity to be continuous at the interface, and constant initial pressure. The flow law for both fluids is a modified Darcy’s law. Coupling the nonlinear flow law with the continuity equation, and taking into account compressibility effects, yields a set of nonlinear second-order partial differential equations. Considering two fluids with the same flow behavior index n allows transformation of the latter equations via a self-similar variable; further transformation of the equations incorporating the conditions at the interface shows for n<1 the existence of a compression front ahead of the moving interface. Solving the resulting set of nonlinear equations yields the positions of the moving interface and compression front, and the pressure distributions; the latter are derived in closed form for any value of n. A sensitivity analysis of the model responses is conducted both in a deterministic and a stochastic framework. In the latter case, Global Sensitivity Analysis (GSA) of the benchmark analytical model is adopted to study how the effects of uncertainty affecting selected parameters: (a) the fluids flow behavior index, (b) the relative total compressibility and mobility in the displaced and displacing fluid domains, and (c) the domain permeability and porosity, propagate to state variables. The relative influence of input parameters on model outputs is evaluated by means of associated Sobol indices, calculated via the Polynomial Chaos Expansion (PCE) technique. The goodness of the results obtained by the PCE is assessed by comparison against a traditional Monte Carlo (MC) approach.     

Lie-group method of solution for steady two-dimensional boundary-layer stagnation-point flow towards a heated stretching sheet placed in a porous medium

The boundary-layer equations for two-dimensional steady flow of an incompressible, viscous fluid near a stagnation point at a heated stretching sheet placed in a porous medium are considered. We apply Lie-group method for determining symmetry reductions of partial differential equations. Lie-group method starts out with a general infinitesimal group of transformations under which the given partial differential equations are invariant. The determining equations are a set of linear differential equations, the solution of which gives the transformation function or the infinitesimals of the dependent and independent variables. After the group has been determined, a solution to the given partial differential equations may be found from the invariant surface condition such that its solution leads to similarity variables that reduce the number of independent variables of the system. The effect of the velocity parameter λ, which is the ratio of the external free stream velocity to the stretching surface velocity, permeability parameter of the porous medium k ₁, and Prandtl number Pr on the horizontal and transverse velocities, temperature profiles, surface heat flux and the wall shear stress, has been studied.     

Boundary layer flow of Reiner-Philippoff fluids

An analysis is made of the boundary layer flow of Reiner-Philippoff fluids. This work is an extension of a previous analysis by Hansen and Na [A.G. Hansen and T.Y. Na, Similarity solutions of laminar, incompressible boundary layer equations of non-Newtonian fluids. ASME 67-WA/FE-2, presented at the ASME Winter Annual Meeting, November (1967)], where the existence of similar solutions of the boundary layer equations of a class of general non-Newtonian fluids were investigated. It was found that similarity solutions exist only for the case of flow over a 90° wedge and, being similar, the solution of the non-linear boundary layer equations can be reduced to the solution of non-linear ordinary differential equations. In this paper, the more general case of the boundary layer flow of Reiner-Philippoff fluids over other body shapes will be considered. A general formulation is given which makes it possible to solve the boundary layer equations for any body shape by a finite-difference technique. As an example, the classical solution of the boundary layer flow over a flat plate, known as the Blasius solution, will be considered. Numerical results are generated for a series of values of the parameters in the Reiner-Philippoff model.