Transient free convection flow of a viscoelastic fluid over vertical surface

In this paper, the viscoelsatic boundary layer flow and the heat transfer near a vertical isothermal impermeable surface and in a quiescent fluid are examined. The gov-erning equations are formulated and solved numerically using MackCormak’s technique. The results show excellent agreement with previously published results by a compari-sion. Representative results for the velocity and temperature profiles, boundary layer thicknesses, Nusselt numbers, and local skin friction coefficients are shown graphically for different values of viscoelsatic parameters. In general, it is found that the velocities increase inside the hydrodynamic boundary layers and the temperatures decrease inside the thermal boundary layers for the viscoelsatic fluid as compared with the Newtonian fluid due to favorable tensile stresses. Consequently, the coefficients of friction and heat transfer enhance for higher viscoelsatic parameters.     

Non-uniform slot suction/injection into mixed convection boundary layer flow over vertical cone

The aim of this work is to study the effect of non-uniform single and double slot suction/injection into a steady mixed convection boundary layer flow over a vertical cone, while the axis of the cone is inline with the flow. The governing boundary layer equations are transformed into a non-dimensional form by a group of non-similar trans- formations. The resulting coupled non-linear partial differential equations are solved nu- merically by employing the quasi-linearization technique and an implicit finite-difference scheme. Numerical computations are performed for different values of the dimensionless parameters to display the velocity and temperature profiles graphically. Also, numerical results are presented for the skin friction and heat transfer coefficients. Results indicate that the skin friction and heat transfer coefficients increase with non-uniform slot suction, but the effect of non-uniform slot injection is just opposite.     

Heat Conduction Problem for Tightly Compressed Adjacent Layers of Plane Plates with Gradually Varying Heat Conduction Coefficients

In dealing with the heat conduction problem of the tightly compressedadjacent plane plates with different constant thermal conduction coefficients,we conventionally solve it layer by layer requiring to fulfill two require-ments of equal temperature and equal heat flux at every two directly con-tacting intermediate surfaces.In this paper,we consider only such a case that the adjacent plane plat-es possess a gradually varying thermal conduction coefficients in orderto offer a chance to solve such problem once for all instead of pursuitingthe old tedious and time-consuming procedure,i.e.,by solving the problemlayer by layer.     

ASYMPTOTIC NON-STABILITY AND BLOW-UP AT BOUNDARY FOR SOLUTIONS OF A FILTRATION EQUATION

For a class of nonlinear filtration equation with nonlinear second-third boundary value condition, it is shown that a priori boundary of the solution can be estimated and controlled by initial data and integral on the boundary of the region. The priori estimate of the solutions was established by iterative method. By using this estimate the solutions may blow-up on the boundary of the region and thus it may have asymptotic non-stability.     

MHD effect of mixed convection boundary-layer flow of Powell-Eyring fluid past nonlinear stretching surface

Sufficient conditions are found for the existence of similar solutions of the mixed convection flow of a Powell-Eyring fluid over a nonlinear stretching permeable sur- face in the presence of magnetic field. To achieve this, one parameter linear group trans- formation is applied. The governing momentum and energy equations are transformed to nonlinear ordinary differential equations by use of a similarity transformation. These equations are solved by the homotopy analysis method （HAM） to obtain the approximate solutions. The effects of magnetic field, suction, and buoyancy on the Powell-Eyring fluid flow with heat transfer inside the boundary layer are analyzed. The effects of the non- Newtonian fluid （Powell-Eyring model） parameters ε and δon the skin friction and local heat transfer coefficients for the cases of aiding and opposite flows are investigated and discussed. It is observed that the momentum boundary layer thickness increases and the thermal boundary layer thickness decreases with the increase in ε whereas the momentum boundary layer thickness decreases and thermal boundary layer thickness increases with the increase in δ for both the aiding and opposing mixed convection flows.     

ON THE HEAT CONDUCTION PROBLEM AND A DESIGN PRINCIPLE FOR PLASMA ARC FURNACE

In this paper a plane heat conduction problem with variable coefficients of heat conductivity K(T) is analysed with given electric power supplied to the plasma arc. The governing equation for unknown temperature distribution is a nonlinear one with a δ function as its nonhomogeneous term. To make the problem attractable by the method of separation of variables, a set of trans-formation of governing equation is introduced.An explicit simple formula is found for the efficiency of the furnace η,η depends linearly on r_0, the nondimensional distance between the arc and surface of melted material, as well as on another nondimensional quantity Q,we describe the above in detail in this paper. This relationship holds for r_0<0.4 and gives a good guidance for the design of furnace.     

Geometric nonlinear formulation for thermal-rigid-flexible coupling system

This paper develops geometric nonlinear hybrid formulation for flexible multibody system with large deformation considering thermal efect. Diferent from the conventional formulation, the heat flux is the function of the rotational angle and the elastic deformation, therefore, the coupling among the temperature, the large overall motion and the elastic deformation should be taken into account. Firstly,based on nonlinear strain–displacement relationship, variational dynamic equations and heat conduction equations for a flexible beam are derived by using virtual work approach,and then, Lagrange dynamics equations and heat conduction equations of the first kind of the flexible multibody system are obtained by leading into the vectors of Lagrange multiplier associated with kinematic and temperature constraint equations. This formulation is used to simulate the thermal included hub-beam system. Comparison of the response between the coupled system and the uncoupled system has revealed the thermal chattering phenomenon. Then, the key parameters for stability, including the moment of inertia of the central body, the incident angle, the damping ratio and the response time ratio, are analyzed. This formulation is also used to simulate a three-link system applied with heat flux. Comparison of the results obtained by the proposed formulation with those obtained by the approximate nonlinear model and the linear model shows the significance of considering all the nonlinear terms in the strain in case of large deformation. At last, applicability of the approximate nonlinear model and the linear model are clarified in detail.     

HEAT FLOW CONTROL BY TRANSFORMATION METHOD WITH GRID GENERATION METHOD

In the context of the transformation method, we propose a general approach to construct numerically the mapping generated by imposing specific boundary conditions with a targeted function, and the necessary material and heat source spatial distributions are then derived with the help of transformation method. The construction of mapping by grid generation method through solving partial differential equations circumvents the limitation of device geometry, which paves the way for designing more complex heat flow control devices. Two numerical examples are also given to show how to design material properties and heat source in order to control temperature patterns.     

A numerical approach based on the meshless collocation method in elastodynamics

In this paper, a collocation technique with the modified equilibrium on line method （ELM） for imposition of Neumann （natural） boundary conditions is presented for solving the two-dimensional problems of linear elastic body vibrations. In the modified ELM, equilibrium over the lines on the natural boundary is satisfied as Neumann boundary condition equations. In other words, the natural boundary conditions are satisfied naturally by using the weak formulation. The performance of the modified version of the ELM is studied for collocation methods based on two different ways to construct meshless shape functions： moving least squares approximation and radial basis point interpolation. Numerical examples of two-dimensional free and forced vibration analyses show that by using the modified ELM, more stable and accurate results would be obtained in comparison with the direct collocation method.     

Novel wavelet-homotopy Galerkin technique for analysis of lid-driven cavity flow and heat transfer with non-uniform boundary conditions

In this paper, a brand-new wavelet-homotopy Galerkin technique is developed to solve nonlinear ordinary or partial differential equations. Before this investigation,few studies have been done for handling nonlinear problems with non-uniform boundary conditions by means of the wavelet Galerkin technique, especially in the field of fluid mechanics and heat transfer. The lid-driven cavity flow and heat transfer are illustrated as a typical example to verify the validity and correctness of this proposed technique. The cavity is subject to the upper and lower walls' motions in the same or opposite directions.The inclined angle of the square cavity is from 0 to π/2. Four different modes including uniform, linear, exponential, and sinusoidal heating are considered on the top and bottom walls, respectively, while the left and right walls are thermally isolated and stationary.A parametric analysis of heating distribution between upper and lower walls including the amplitude ratio from 0 to 1 and the phase deviation from 0 to 2π is conducted. The governing equations are non-dimensionalized in terms of the stream function-vorticity formulation and the temperature distribution function and then solved analytically subject to various boundary conditions. Comparisons with previous publications are given,showing high efficiency and great feasibility of the proposed technique.     

The effects of variable specific heat on the stability of hypersonic boundary layer on a flat plate

The gas temperature within hypersonic boundary layer flow is so high that the specific heat of gas is no longer a constant but relates to temperature. How variable specific heat influences on boundary layer flow stability is worth researching. The effect of the variable specific heat on the stability of hypersonic boundary layer flows is studied and compared with the case of constant specific heat based on the linear stability theory. It is found that the variable specific heat indeed has some effects on the neutral curves of both the first-mode and the second-mode waves and on the maximum rate of growth also. Therefore, the relationship between specific heat and temperature should be considered in the study of the stability of the boundary layer.     

Prediction of hypersonic boundary layer transition on sharp wedge flow considering variable specific heat

When the air temperature reaches 600 K or higher, vibration is excited. The specific heat is not a constant but a function of temperature. Under this condition, the transition position of hypersonic sharp wedge boundary layer is predicted by using the improved eN method considering variable specific heat. The transition positions with different Mach numbers of oncoming flow, half wedge angles, and wall conditions are computed condition, the nearer to the Mach number The results show that for the same oncoming flow condition and wall transition positions of hypersonic sharp wedge boundary layer move much leading edge than those of the flat plate. The greater the oncoming flow the closer the transition position to the leading edge.     

A CRITERION FOR THE STABILITY OF MOTION OF NONLINEAR SYSTEM

As to an autonomous nonlinear system, the stability of the equilibrium state in a sufficiently small neighborhood of the equilibrium state can be determined by eigenvalues of the linear part of the nonlinear system provided that the eigenvalues are not in a critical case. Many methods may be used to detect the stability for a linear system. A lot of researches for determining the stability of a nonlinear system are completed by mathematicians and mechanicians but most of them are methods for the special forms of nonlinear systems. Till now, none of these methods can be conveniently applied to all nonlinear systems. The method introduced by this paper gives the necessary and sufficient conditions of the stability of a nonlinear system. The familiar Krasovski's method is a special case of this method     

The numerical method for analysing the transient temperature field and transient phase transformation of metal materials in heat treating

In this paper,the Kirchhoffs transformation is popularized to the nonlinear heat conduction problem which the heat conductivity can be expressd as a multinomial of temperature firstly,the boundary condition of heat conduction problem is determined by analytics.Secondly,the incubation peroid superposition and the linear combination law is employed to simulate the transient phasses transformation in the process of heat treatment of materials.That the begin time of phase transformation,the type of phase transformation and the amount of phase constitution is determined simply.Finally,the three-dimension Dual Reciprocity Boundary Element Method is usedto analysis the total process of various heat treatment of component,the results of numerical calculation of examples show that the method provided in this paper is effectivce.     

Nonlinear dynamical systems and bistability in linearly forced isotropic turbulence

In this paper, we present an extensive study of the linearly forced isotropic turbulence. By using analytical method, we identify two parametric choices, of which they seem to be new as far as our knowledge goes. We prove that the underlying nonlinear dynamical system for linearly forced isotropic turbulence is the general case of a cubic Lienard equation with linear damping. We also discuss a FokkerPlanck approach to this new dynamical system, which is bistable and exhibits two asymmetric and asymptotically stable stationary probability densities.     

传热传质对有抽吸的收缩薄片上的非线性磁流体动力学边界层流动的影响

This work is concerned with Magnetohydrodynamic viscous flow due to a shrinking sheet in the presence of suction. The cases of two dimensional and axisymmetric shrinking are discussed. The governing boundary layer equations are written into a dimensionless form by similarity transformations. The transformed coupled nonlinear ordinary differential equations are numerically solved by using an advanced numeric technique. Favorability comparisons with previously published work are presented. Numerical results for the dimensionless velocity, temperature and concentration profiles as well as for the skin friction, heat and mass transfer and deposition rate are obtained and displayed graphically for pertinent parameters to show interesting aspects of the solution.     

THE METHOD OF COMPOSITE EXPANSIONS APPLIED TO BOUNDARY LAYER PROBLEMS IN SYMMETRIC BENDING OF THE SPHERICAL SHELLS

In this paper,the method of composite expansions which was proposed by W.Z.Chien(1948)is extended to investigate two-parameter boundary layer problems.For the problems of symmetric deformations of the spherical shells under the action of uniformly distribution load q,its nonlinear equilibrium equations can be written as follows:where εand δare undetermined parameters.If δ=1 and εis a small parameter,the above-mentioned problem is called first boundary layer problem;if εis a small parameter,and δis a small parameter,too,the above-mentioned problem is called second boundary layer problem.For the above-mentioned problems,however,we assume that the constants ε,δand p satisfy the following equation:In defining this condition by using the extended method of composite expansions,we find the asymptotic solution of the above-mentioned problems with the clamped boundary conditions.     

FIRST INTEGRAL AND STABILITY OF MOTION IN THE CRITICAL CASES

In this paper. we indicate that after the Liapunov function by using linear combinationof mechanical first integral was suggested by Chetayev in 1946.He and his students solvedstability of conservative system by means of this method. But he had trouble to solve theproblems by means of cut and try.Moreover,the condition of stability is imperfect.Solution by this method is limited for problems of purely imaginary roots.The cases of zeroroots have not been considered Condition of stability secured is more strict.This paper suggests that the differential equation can be transformed into standardform by method of cancellation of cyclic coordinates(method of lowering degree of order),and condition of stability can be determined by energy integral.By this method not only thecomputation is clear and concise.But also zero roots can be considered.Therefore theproblems of two cyclic coordinates can be transformed into second-order system,and we getnew conclusion of the condition of stability simply .As for problem     

Studies on nonlinear stability of three-dimensional H-type disturbance

The three-dimensional H-type nonlinear evolution process for the problem of boundary layer stability is studied by using a newly developed method called parabolic stability equations (PSE). The key initial conditions for sub-harmonic disturbances are obtained by means of the secondaryinstability theory. The initial solutions of two-dimensional harmonic waves are expressed in Landau expansions. The numerical techniques developed in this paper, including the higher order spectrum method and the more effective algebraic mapping for dealing with the problem of an infinite region, increase the numerical accuracy and the rate of convergence greatly. With the predictor-corrector approach in the marching procedure, the normalization, which is very important for PSE method, is satisfied and the stability of the numerical calculation can be assured. The effects of different pressure gradients, including the favorable and adverse pressure gradients of the basic flow, on the “H-type“ evolution are studied in detail. The results of the three-dimensional nonlinear “H-type“ evolution are given accurately and show good agreement with the data of the experiment and the results of the DNS from the curves of the amplitude variation, disturbance velocity profile and the evolution of velocity.     

Effects of nonlinear velocity slip and temperature jump on pseudo-plastic power-law fluid over moving permeable surface in presence of magnetic field

Flow and heat transfer of a pseudo-plastic power-law fluid over a stretching permeable surface with the magnetic effect is investigated. In the boundary conditions,the nonlinear temperature jump and the velocity slip are considered. Semi-similarity equations are obtained and solved by bvp4c with MATLAB. The problem can be considered as an extension of the previous work done by Mahmoud(Mahmoud, M. A. A. Slip velocity effect on a non-Newtonian power-law fluid over a moving permeable surface with heat generation. Mathematical and Computer Modelling, 54, 1228–1237(2011)). Efforts are made to discuss the effects of the power-law number, slip velocity, and temperature jump on the dimensionless velocity and temperature distribution.