Boundary layer equations and stretching sheet solutions for the modified second grade fluid

A modified second grade non-Newtonian fluid model is considered. The model is a combination of power-law and second grade fluids in which the fluid may exhibit normal stresses, shear thinning or shear thickening behaviors. The equations of motion are derived for two dimensional incompressible flows. The boundary layer equations are derived from the equations. Symmetries of the boundary layer equations are calculated using Lie Group theory. For a special power law index of m = −1, the principal Lie algebra extends. Using one of the symmetries, the partial differential system is transferred to an ordinary differential system. The ordinary differential equations are numerically integrated for the stretching sheet boundary conditions. Effects of power-law index and second grade coefficient on the boundary layers are shown and solutions are contrasted with the usual second grade fluid solutions. The shear stress on the boundary is also calculated.     

MHD 3D free convective flow of nanofluid over an exponentially stretching sheet with chemical reaction

This paper deals with a problem where the effect of variable magnetic field and chemical reaction on free convective flow of an electrically conducting incompressible water based nanofluid over an exponentially stretching sheet has been investigated. In the present study, Buongiorno model associated with Brownian motion and thermophoretic diffusion is employed to describe the heat transfer enhancement of nanofluids. Some suitable similarity transformations reduced the governing boundary layer non-linear partial differential equations into a set of ordinary non-linear differential equations. The transformed equations are then solved numerically using fourth order Runga-Kutta method along with Shooting technique. The major outcomes of the present study is that the magnetic field impedes the fluid motion while thermal as well as mass buoyancy forces accelerate it, the thermophoretic diffusion enhances dimensionless fluid temperature as well as concentration leading to thicker thermal and concentration boundary layers. On the other hand, concentration exponent, Brownian motion parameter and chemical reaction parameter exhibit reverse trend on temperature and concentration. In addition, the presence of magnetic field under the influence of thermal as well as mass buoyancies supports to reduce the rate of heat transfer as well as wall shear stress while the first order chemical reaction develops a thinner concentration boundary layer.     

Dynamics of the dragonfly wings raised by blood circulation

We investigate how the blood flow in the veins in the flapping wings of a dragonfly affects their dynamic response. An idealized model of an elastic tube conveying fluid and rotating around a fixed axis is adopted in this study, based on which governing partial differential equations of motion are obtained by invoking the extended Hamilton’s principle. Separation of variables techniques and assumed modes method are employed to solve the resulting equations, and the stabilization analysis is performed to assess the stability of the system. In particular, the coupling effects of tube rotation, deformation, and the movement of the fluid inside are evaluated under different flow rates and rotation speeds. This demonstrates that if the blood in the dragonfly wings flows from humeral angle distally to the wing apex, a stabilization effect can be obtained, and the higher the blood flow rate is, the faster the system will be stabilized. Contrary cases are also studied for further validation of the model.     

Coupled heat and mass transfer during nonisothermal absorption by falling droplet with internal circulation

We considered mass and heat transfer during nonisothermal absorption of a gas by a falling droplet with internal circulation. Gas phase is assumed to be free of inert admixtures and mass transfer is liquid phase controlled. Mass flux is directed from a gaseous phase to a droplet, and the interfacial shear stress causes a fluid flow inside the droplet. Droplet deformation under the influence of interface shear stress is neglected. Absorbate accumulation and temperature increase in the bulk of liquid phase are taken into account. The problem is solved in the approximations of a thin concentration and temperature boundary layers in the liquid phase. The thermodynamic parameters of the system are assumed constant. The system of transient partial parabolic differential equations of convective diffusion and energy balance with time-dependent boundary conditions is solved by combining the similarity transformation method with Duhamel's theorem, and the solution is obtained in a form of Volterra integral equation of the second kind which is solved numerically. Theoretical results are compared with available experimental data for water vapor absorption by falling droplets of aqueous solution of LiBr.     

Coupled longitudinal-transverse behaviour of a geometrically imperfect microbeam

Based on the modified couple stress theory, the coupled longitudinal-transverse nonlinear behaviour of an imperfect microbeam is investigated numerically. The equations governing the longitudinal and transverse motions are obtained using Hamilton’s principle for the system with an initial geometric imperfection. The Galerkin scheme is employed to discretize the two partial differential equations of motion, yielding a set of second-order nonlinear ordinary differential equations with coupled terms. This set is cast into new set of first-order nonlinear ordinary differential equations and solved by means of the pseudo-arclength continuation technique. The nonlinear resonant response of the system along with bifurcations are presented via frequency–response curves. Moreover, the effect of different system parameter on the frequency–response curves is highlighted.     

A least-squares front-tracking finite element method analysis of phase change with natural convection

This paper focuses on the numerical modelling of phase-change processes with natural convection. In particular, two-dimensional solidification and melting problems are studied for pure metals using an energy preserving deforming finite element model. The transient Navier–Stokes equations for incompressible fluid flow are solved simultaneously with the transient heat flow equations and the Stefan condition. A least-squares variational finite element method formulation is implemented for both the heat flow and fluid flow equations. The Boussinesq approximation is used to generate the bulk fluid motion in the melt. The mesh motion and mesh generation schemes are performed dynamically using a transfinite mapping. The consistent penalty method is used for modelling incompressibility. The effect of natural convection on the solid/liquid interface motion, the solidification rate and the temperature gradients is found to be important. The proposed method does not possess some of the false diffusion problems associated with the standard Galerkin formulations and it is shown to produce accurate numerical solutions for convection dominated phase-change problems.     

The nonlinear vibration of an initially stressed laminated plate

Nonlinear partial differential equations of motion for a laminated plate in a general state of non-uniform initial stress are presented in various plate theories. This study uses Lo’s displacement field to derive the governing equations. The higher-order terms in Lo’s theory can be disregarded, to obtain the equations of simpler forms and even other theories for laminated plate. These nonlinear partial equations are transformed to ordinary nonlinear differential equations using the Galerkin method. The Runge–Kutta method is used to obtain the ratio of nonlinear frequency to linear frequency. The numerical solutions of an initially stressed laminate plate based on various plate theories obtained by the Galerkin and Runge–Kutta method are presented herein. Using these equations with various theories, the nonlinear vibration behavior of laminated plate is studied. The results show that apparent discrepancies exist among the various displacement fields, which indicates the transverse shear strain, normal strain and initial stress state have great effect on the vibration behavior of laminate plate under nonlinear vibration.     

严重段塞流引起的海洋立管振动响应

基于改进的严重段塞流瞬态数学模型和平面刚架理论,建立严重段塞流海洋立管耦合振动数学模型,对数学模型进行求解,研究了严重段塞流引起的海洋立管振动响应。数值模拟过程中,采用欧拉法计算严重段塞流流动参数,利用Galerkin法对立管结构动力学方程进行有限元离散,Newmark-β法求解离散方程。为了提高计算效率和精度,采用了变时间步长逐步积分方案。将数值模拟结果与实验数据进行对比,验证了数学模型和数值方法的准确性。对严重段塞流引起的立管系统位移响应、内力和支承力变化进行了深入分析。结果表明:立管系统的振动响应、内力变化规律与严重段塞流的周期性特征密切相关;弹性基础可以极大降低管道结构的振动幅度及内力,尤其是下倾管的弯曲内力;上升管的高频振动幅度较大,轴力和弯曲内力变化也较大。严重段塞流引起海洋立管振动响应的分析对海洋立管设计及其支承防护具有重要的指导意义。     

Influence of temperature dependent viscosity on the MHD-channel flow of dusty fluid with heat transfer

Summary This paper studies the effect of variable viscosity on the transient flow of dusty fluid with heat transfer. The fluid is acted upon by a constant pressure gradient, and an external uniform magnetic field is applied perpendieular to the plates. The governing nonlinear partial differential equations are solved numerically, and some important effects for the variable viscosity and the uniform magnetic field on the transient flow and heat transfer of both the fluid and dust particles are indicated.     

复合材料悬臂板系统的多脉冲同宿轨道

针对复合材料悬臂板系统的非线性动力学行为进行了分析。模型考虑高阶横向剪切效应、几何大变形和横向阻尼的影响,基于Reddy的高阶剪切变形理von Karman的大变形理论,利用Hamilton原理,Galerkin离散和多尺度法得到系统横向位移的平均方程。应用广义Melnikov方法研究了复合材料悬臂板的非线性混沌动力学行为。分析了在共振带附近,复合材料悬臂板系统存在的Shilnikov类型多脉冲跳跃同宿轨道。最后结合数值模拟,进一步揭示系统存在多脉冲跳跃现象。     

Coupled size-dependent behavior of shear deformable microplates

The aim of the present study is to investigate the nonlinear motion characteristics of a shear deformable microplate based on the modified couple stress theory. The microplate is modeled via the third-order shear deformation theory retaining in-plane displacements and inertia. Using the Lagrange equations together with an assumed-mode method, five sets of second-order nonlinear ordinary differential equations of motion with coupled terms are obtained. These five sets of equations (two for the in-plane motions, one for the out-of-plane motion, and two for rotations) are transformed into ten sets of first-order nonlinear ordinary differential equations. These resultant equations are then solved by means of a direct time integration technique and the pseudo-arclength continuation method in order to analyze the nonlinear response of the system. Apart from the nonlinear analysis, the linear natural frequencies of the system are obtained using an eigenvalue analysis. Results are shown through frequency–response and force–response curves. Points of interest in the parameter space in the form of time histories, phase-plane portraits, and fast Fourier transforms are also highlighted. Moreover, a comparison is made between the motion characteristics of the system based on the modified couple stress and classical continuum theories.     

Hydromagnetic flow and heat transfer of a second-grade viscoelastic fluid in a channel with oscillatory stretching walls: application to the dynamics of blood flow

The characteristics of flow and heat transfer of a fluid in a channel with oscillatory stretching walls in the presence of an externally applied magnetic field are investigated. The fluid considered is a second-grade viscoelastic electrically conducting fluid. The partial differential equations that govern the flow are solved by developing a suitable numerical technique. The computational results for the velocity, temperature and the wall shear stress are presented graphically. The study reveals that flow reversal takes place near the central line of the channel. This flow reversal can be reduced to a considerable extent by applying a strong external magnetic field. The results are found to be in good agreement with those of earlier investigations.     

Coupled longitudinal-transverse-rotational behaviour of shear deformable microbeams

In this paper, for the first time, the nonlinear motion characteristics of a hinged-hinged third-order shear deformable microbeam are examined, based on the modified couple stress theory and the third-order shear deformation theory. The extensibility of the microbeam is modelled by taking into account the longitudinal displacement. The nonlinear equations governing the longitudinal, transverse, and rotational motions are derived by means of Hamilton's principle in conjunction with the modified couple stress theory (to take into account small-scale effects). The three coupled nonlinear partial differential equations are discretized via the Galerkin method and the resulting set of ordinary differential equations is solved by means of the pseudo-arclength continuation technique and via direct time-integration. The effects of the system parameters on the behaviour of the microbeam are studied. Results are presented in the form of frequency-responses and force-responses. Points of interest in the parameter space are also highlighted in the form of time histories, phase-plane portraits, and fast Fourier transforms (FFTs). Moreover, the similarities and differences in the response of the system obtained via the modified couple stress and classical continuum mechanics theories are discussed.     

Dynamic stress intensity factors due to scattering of harmonic waves by a crack in an infinite medium

An analytical method for calculating dynamic stress intensity factors in the mixed mode (combination of opening and sliding modes) using complex functions theory is presented. The crack is in infinite medium and subjected to the plane harmonic waves. The basis of the method is grounded on solving the two‐dimensional wave equations in the frequency domain and complex plane using mapping technique. In this domain, solution of the resulting partial differential equations is found in the series of the Hankel functions with unknown coefficients. Applying the boundary conditions of the crack, these coefficients are calculated. After solving the wave equations, the stress and displacement fields, also the J‐integrals are obtained. Finally using the J‐integrals, dynamic stress intensity factors are calculated. Numerical results including the values of dynamic stress intensity factors for a crack in an infinite medium subjected to the dilatation and shear harmonic waves are presented.     

Accuracy and security analysis of transient flows in relatively long pipelines

Consideration is given, in this paper, to the accurate numerical solution of transient flows in relatively long pipelines caused by the water hammer phenomenon. The governing equations for such flows are a set of two coupled non-linear hyperbolic partial differential equations where the friction factor is a very important parameter. In these equations, the gravity force is introduced to take into account the pipeline inclination. The mathematical equations are solved in the time domain by the method of characteristics using linear integration. To reproduce correctly the transient flow, very small time increments were used in the constructed computer program. This program permits to get some systematic indications on the evolution and the damping of the pressure head waves due to a fast closing valve at the downstream end of a long pipeline. The study shows that results are different from those of short pipelines and small viscous fluids. To check the validity of the numerical model, computed results have been successfully compared with those found in the relevant literature. These results show that the gravity lift may have an important effect on the maximum pressures, which may become very important near the valve and cause failure of the pipeline especially in presence of defect. The safety factor, computed at equidistant sections of the line, determines the distance between the reservoir and the defect from which the failure may happen. For the considered application, results show that the pipeline is safe near the supply upstream tank.     

Accurate solution for free vibration analysis of functionally graded thick rectangular plates resting on elastic foundation

Vibration analysis of a functionally graded rectangular plate resting on two parameter elastic foundation is presented here. The displacement filed based on the third order shear deformation plate theory is used. By considering the in-plane displacement components of an arbitrary material point on the mid-plane of the plate and using Hamilton’s principle, the governing equations of motion are obtained which are five highly coupled partial differential equations. An analytical approach is employed to decouple these partial differential equations. The decoupled equations of functionally graded rectangular plate resting on elastic foundation are solved analytically for levy type of boundary conditions. The numerical results are presented and discussed for a wide range of plate and foundation parameters. The results show that the Pasternak (shear) elastic foundation drastically changes the natural frequency. It is also observed that in some boundary conditions, the in-plane displacements have significant effects on natural frequency of thick functionally graded plates and they cannot be ignored.     

BEM modeling of surface water wave motion with laminar boundary layers

This study is concerned with numerical modeling of viscous surface wave motion using boundary element method (BEM). The equations of motion for thin boundary layers at the solid surfaces are coupled with the potential flow in the bulk of the fluid, and a mixed BEM-finite difference technique is used to obtain the viscosity-related quantities such as wave damping rate, shear stress, and velocity distribution inside the boundary layer. The technique is presented for standing surface wave motion. An excellent agreement is obtained between the numerical predictions and the previous results. The extension to other free surface problems is straightforward.     

电流变技术在抑制深孔切削颤振中的应用

针对深孔加工中出现的切削颤振问题,利用电流变液的电流变效应,设计一套基于流动和剪切混合工作模式的电流变减振器。通过对深孔切削系统的动力学模型的建立和分析,得到系统的运动微分方程。利用MATLAB/Sinmulink对动力学模型进行仿真分析并进行实验验证,结果表明电流变减振器能够减小颤振的振幅,且在不同电场强度下具有不同的减振效果。因此通过控制电流变减振器的电场强度可以有效地抑制深孔机床中切削颤振的发生。     

Determination of high‐order fields for multianisotropic material antiplane V‐notches and inclusions by the asymptotic expansion technique and an overdeterministic method

In this paper, the singular behavior for anisotropic multimaterial V‐notched plates is investigated under antiplane shear loading condition. Firstly, the elastic governing equations are transformed into eigen ordinary differential equations through introducing the asymptotic expansions of displacements near the notch tip. The stress singularity exponents, including the higher‐order terms, and the corresponding eigen angular functions are then obtained by solving the established equations by using the interpolating matrix method. Thus, using the combination of the results from finite element analyses and the derived asymptotic expansion, an overdeterministic method is employed to calculate the amplitudes of the coefficients in the asymptotic expansions. Finally, the stress and displacement fields in the vicinity of the notch tip, consisting of both singular terms and higher‐order terms, are determined. The effects of material properties and geometry characteristic on the singular behaviour of the notch tip are discussed in detail.     

On the flow of an electrically conducting nonlocal viscous fluid between parallel plates in the presence of a transverse magnetic field in magnetohydrodynamics

A continuum theory is constructed for the flow of an electrically conducting nonlocal viscous fluid between two nonconducting parallel plates. The flow is subject to the influence of a transverse magnetic field. The effects of long range or nonlocal interactions at a material point in the fluid arising from all material points in the rest of the fluid are taken into account by means of a nonlocal influence function. Equations of motion governing the nonlocal viscous flow are derived from localized forms of global balance laws and constitutive equations appropriate for electromagnetically active media. These field equations are analytically solved for the nonlocal velocity and the nonlocal stress fields. The effects of varying the magnetic field strength on the shear stress are investigated. The effects of such variations on the shear stress exerted on the walls of microscopic channels are also determined. Numerical computations are provided for these results.