A variational analytical approach to time dependent flow of oldroyd B fluid in circular tube

In the present investigation the time dependent flow of an Oldroyd fluid B in a horizontal cylindrical pipe is stuided by the variational analytical approach developed by author. The time dependent problem is mathematically reduced to a partial differential equation of third order. Using the improved variational approach due to Kantorovich the partial differential equation can be reduced to a system of ordinary differential equations for different approximations. The ordinary differential equations are solved by the method of the Laplace transform which is led to an analytical form of the solutions. Project supported by TWAS and Chinese Academy of Sciences and the National Science Foundation of China     

A VARIATIONAL ANALYTICAL APPROACH TO TIME DEPENDENT FLOW OF OLDROYD B FLUID IN CIRCULAR TUBE

ＡＶＡＲＩＡＴＩＯＮＡＬＡＮＡＬＹＴＩＣＡＬＡＰＰＲＯＡＣＨＴＯＴＩＭＥＤＥＰＥＮＤＥＮＴＦＬＯＷＯＦＯＬＤＲＯＹＤＢＦＬＵＩＤＩＮＣＩＲＣＵＬＡＲＴＵＢＥＨａｎｓｈｉ－ｆａｎｓ（韩式方）（Ｒｅｓ．Ｌａｂ．Ｎｏｎ－ＮｅｗｔｏｎｉａｎＦｌｕｉｄＭｅｃｈ...     

Solution of the Rayleigh problem for a powerlaw non-Newtonian conducting fluid via group method

ＩｎｔｒｏｄｕｃｔｉｏｎＦｌｕｉｄｓｔｈａｔｏｂｅｙＮｅｗｔｏｎ’ｓｌａｗｏｆｖｉｓｃｏｓｉｔｙａｒｅｃａｌｌｅｄＮｅｗｔｏｎｉａｎｆｌｕｉｄｓ.Ｎｅｗｔｏｎ’ｓｌａｗｏｆｖｉｓｃｏｓｉｔｙｉｓτ=μｄｕ/ｄｔ,ｗｈｅｒｅτｉｓｔｈｅｓｈｅａｒｓｔｒｅｓｓａｎｄμｉｓｔｈｅｖｉｓｃｏｓｉｔｙ .ＮｏｔａｌｌｆｌｕｉｄｓｆｏｌｌｏｗｔｈｅＮｅｗｔｏｎｉａｎｓｔｒｅｓｓ＿ｓｔｒａｉｎｒｅｌａｔｉｏｎ .Ｓｏｍｅｆｌｕｉｄｓ ,ｓｕｃｈａｓｋｅｔｃｈｕｐ ,ａｒｅｓｈｅａｒｔｈｉｎｎｉｎｇ ,ｔｈａｔｉｓ,ｔｈｅｃ…     

Computer Simulation of Non-Newtonian Flow and Mass Transport Through Coronary Arterial Stenosis

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｍｏｓｔｏｂｖｉｏｕｓｃｈａｒａｃｔｅｒｏｆａｔｈｅｒｏｓｃｌｅｒｏｓｉｓｉｓｔｈｅｌｏｃａｌｉｎｔｉｍａｌａｃｃｕｍｕｌａｔｉｏｎｏｆｌｏｗ＿ｄｅｎｓｉｔｙｌｉｐｏｐｒｏｔｅｉｎ(ＬＤＬ) ,ｔｈｅｉｎｃｒｅａｓｅｏｆｔｈｅｆｉｂｅｒａｌｔｉｓｓｕｅａｎｄｔｈｅｓｔｅｎｏｓｉｓ.Ｏｎｏｎｅｈａｎｄ ,ｉｔｉｓｆｏｕｎｄｔｈａｔｔｈｅｄｉｓｅａｓｅｏｆｔｅｎｏｃｃｕｒｓａｔｔｈｅｃｏｍｐｌｅｘｇｅｏｍｅｔｒｙｒｅｇｉｏｎ ,ｓｕｃｈａｓｂｉｆｕｒｃａｔｉｏｎｚｏｎｅ,ｃｕｒｅ…     

Numerical study of the axially symmetric motion of an incompressible viscous fluid in an annulus between two concentric rotating spheres

The research reported herein involved the study of the transient motion of a system consisting of an incompressible Newtonian fluid in an annulus between two concentric, rotating, rigid spheres. The primary purpose of the research was to study the use of a numerical method for analysing the transient motion that results from the interaction between the fluid in the annulus and the spheres which are started suddenly by the action of prescribed torques. The problems considered in this research included cases where: (a) one or both spheres rotate with prescribed constant angular velocities and (b) one sphere rotates due to the action of an applied constant or impulsive t?orque. In this research the coupled solid and fluid equations were solved numerically by employing the finite difference technique. With the approach adopted in this research, only the derivatives with respect to spatial variables were approximated with the use of the finite difference formulae. The steady state problem was also solved as a separate problem (for verification purposes), and the results were compared with those obtained from the solution of the transient problem. Newton's algorithm was employed to solve the algebraic equations which resulted from the steady state problem, and the Adams fourth-order predictor–corrector method was employed to solve the ordinary differential equations for the transient problem. Results were obtained for the streamfunction, circumferential function, angular velocity of the spheres and viscous torques acting on the spheres as a function of time for various values of the system dimensionless parameters.     

Sensitivity analysis of low Reynolds number channel flow using the finite volume method

An unsteady finite volume‐based fractional step algorithm solved on a staggered grid has been developed for computing design sensitivity parameters in two‐dimensional flows. Verification of the numerical code is performed for the case of low Reynolds number, pressure‐driven flow through a straight channel, which has an exact steady‐state solution to the Navier–Stokes equations. Sensitivity of the flow to the channel height, fluid viscosity, and imposed pressure gradient is considered. Three different numerical techniques for computing the design sensitivity parameters: finite difference, complex‐step differentiation, and sensitivity equation method (SEM), are compared in terms of numerical error (relative to the exact solution), computational expense, and ease of implementation. Results indicate that, of all the three methods, complex step is the most accurate and requires the least computational time. In addition, treatment of the boundary conditions in SEM is addressed, within the framework of the present finite volume approach, with special attention given to parameter dependence in the boundary conditions. Error estimation based on the Grid Convergence Index provides a good indication of the exact error in the SEM solutions. An example of application of the use of sensitivity parameters to estimate the propagation of input uncertainty through the numerical simulation is also provided. Copyright © 2007 John Wiley & Sons, Ltd.     

Flow of a biomagnetic viscoelastic fluid: application to estimation of blood flow in arteries during electromagnetic hyperthermia,a therapeutic procedure for cancer treatment

The paper deals with the theoretical investigation of a fundamental problem of biomaguetic fluid flow through a porous medium subject to a magnetic field by using the principles of biomagnetic fluid dynamics （BFD）. The study pertains to a situation where magnetization of the fluid varies with temperature. The fluid is considered to be non-Newtonian, whose flow is governed by the equation of a second-grade viscoelastic fluid. The walls of the channel are assumed to be stretchable, where the surface velocity is proportional to the longitudinal distance from the origin of coordinates. The problem is first reduced to solving a system of coupled nonlinear differential equations involving seven parameters. Considering blood as a biomagnetic fluid and using the present analysis, an attempt is made to compute some parameters of the blood flow by developing a suitable numerical method and by devising an appropriate finite difference scheme. The computational results are presented in graphical form, and thereby some theoretical predictions are made with respect to the hemodynamical flow of the blood in a hyperthermal state under the action of a magnetic field. The results clearly indicate that the presence of a magnetic dipole bears the potential so as to affect the characteristics of the blood flow in arteries to a significant extent during the therapeutic procedure of electromagnetic hyperthermia. The study will attract the attention of clinicians, to whom the results would be useful in the treatment of cancer patients by the method of electromagnetic hyperthermia.     

ODE conversion techniques and their applications in computational mechanics

In this paper, a number of ordinary differential equation (ODE) conversion techniques for transformation of nonstandard ODE boundary value problems into standard forms are summarised, together with their applications to a variety of boundary value problems in computational solid mechanics, such as eigenvalue problem, geometrical and material nonlinear problem, elastic contact problem and optimal design problems through some simple and representative examples. The advantage of such approach is that various ODE boundary value problems in computational mechanics can be solved effectively in a unified manner by invoking a standard ODE solver. The project is supported by National Natural Science Foundation of China.     

基于黎曼几何的非线性振子动力学分析递推解析算法

基于黎曼几何和变分原理,推导了黎曼流形上非线性耗散动力系统的二阶微分动力学方程,并运用流形收缩的概念将动力学方程离散化,进而建立了相应的递推求解格式。选取3个自治非线性阻尼振子系统,分别采用递推解析算法和龙格库塔法求解微分动力学方程,并比较分析了不同的时间步长下两种算法的计算耗时。结果表明,与龙格库塔法相比,基于黎曼几何的递推算法不仅能得到每一时步的解析表达式,而且计算耗时短,计算效率高。基于黎曼流形的动力学方程递推算法为非线性动力学系统的解析求解提供了新思路。     

THE TRANSIENT ELLIPTIC FLOW OF POWER-LAW FLUID IN FRACTAL POROUS MEDIA

1　ＴｈｅＦｌｏｗＭｏｄｅｌｏｆＰｏｗｅｒ＿ＬａｗＦｌｕｉｄｉｎＲａｄｉｃａｌＦｒａｃｔａｌＲｅｓｅｒｖｏｉｒＴｈｅｔｒａｎｓｉｅｎｔｆｌｏｗｏｆｐｏｗｅｒ＿ｌａｗｆｌｕｉｄｉｎｒａｄｉｃａｌｆｒａｃｔａｌｒｅｓｅｒｖｏｉｒｉｓｓｔｕｄｉｅｄｉｎＲｅｆ.[1 ] ,ａｎｄａｎａｌｙｔｉｃａｌｓｏｌｕｔｉｏｎｏｆＬａｐｌａｃｅｓｐａｃｅｉｓｄｅｒｉｖｅｄ .ＩｎＲｅｆ.[2 ] ,ｔｈｅｔｒａｎｓｉｅｎｔｅｌｌｉｐｔｉｃａｌｆｌｏｗｉｓｒｅｓｅａｒｃｈｅｄｏｎｍｏｄｅｌｏｆｅｘｐａｎｄｉｎｇｒｅｃｔａｎｇｌｅ .Ｔ…     

COMPARISON OF ANALYTICAL SOLUTIONS FOR CMSMPR CRYSTALLIZER WITH QMOM POPULATION BALANCE MODELING IN FLUENT

Fluent version 6.2 computational fluid dynamics environment has been enhanced with a population balance capability that operates in conjunction with its multiphase calculations to predict the particle size distribution within the flow field. The population balance is solved by the quadrature method of moments （QMOM）. Fluent＇s prediction capabilities are tested by using a 2-dimensional analogy of a constantly stirred tank reactor with a fluid flow compartment that mixes the fluid quickly and efficiently using wall movement and has a feed stream and a product stream. The results of these Fluent simulations using QMOM population balance solver are compared to steady state analytical solutions for the population balance in a stirred tank where 1） growth, 2） aggregation, and 3） breakage, take place separately and 4） combined nucleation and growth and 5） combined nucleation, growth and aggregation take place. The results of these comparisons show that the moments of the population balance are accurately predicted for nucleation, growth, aggregation and breakage when the flow field is turbulent. With laminar flow the mixing is not ideal and as a result the steady state well mixed solutions are not accurately simulated.     

On the solution of MHD flow over a nonlinear stretching sheet by an efficient semi‐analytical technique

The problem of magneto‐hydrodynamic fluid flow past a nonlinear stretching sheet in the presence of a transverse magnetic field is analyzed. The governing equations are transformed into a nonlinear ordinary differential equation that is solved using a novel spectral homotopy analysis method and the Matlab in‐built numerical solverttbvp4c. The new technique removes some known limitations of the homotopy analysis method and offers a more systematic way of selecting initial approximations and the optimal auxiliary parameter ?. A comparison with the numerical solution confirms the robustness, the computational efficiency, and the accuracy of the technique. Copyright © 2011 John Wiley & Sons, Ltd.     

On the case when steady converging/diverging flow of a non-Newtonian fluid in a round cone permits an exact solution

This communication considers the steady converging/diverging flow of a non-Newtonian viscous power-law fluid in a round cone. The motion is driven by a sink/source of mass at the origin. It is shown that the problem permits exact similarity solution for a particular value (n=4/3) of the fluid index. In this case a complete set of governing equations can be reduced to an ordinary differential equation, which is solved numerically for different values of the main non-dimensional parameters (the cone angle and the dimensionless sink/source intensity).     

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).     

On a difficulty in the formulation of initial and boundary conditions for eigenfunction expansion solutions for the start-up of fluid flow

Most mathematics and engineering textbooks describe the process of “subtracting off” the steady state of a linear parabolic partial differential equation as a technique for obtaining a boundary-value problem with homogeneous boundary conditions that can be solved by separation of variables (i.e., eigenfunction expansions). While this method produces the correct solution for the start-up of the flow of, e.g., a Newtonian fluid between parallel plates, it can lead to erroneous solutions to the corresponding problem for a class of non-Newtonian fluids. We show that the reason for this is the non-rigorous enforcement of the start-up condition in the textbook approach, which leads to a violation of the principle of causality. Nevertheless, these boundary-value problems can be solved correctly using eigenfunction expansions, and we present the formulation that makes this possible (in essence, an application of Duhamel's principle). The solutions obtained by this new approach are shown to agree identically with those obtained by using the Laplace transform in time only, a technique that enforces the proper start-up condition implicitly (hence, the same error cannot be committed).     

Coupling of the methods of successive approximations and undetermined coefficients for the prediction of the thermal behaviour of uniform circumferential fins

This study addresses a distinct, unsophisticated computational procedure for solving approximately, but analytically, the one-dimensional heat equation for circumferential fins of uniform thickness with constant properties. This differential equation with variable coefficients, called the modified Bessel equation of zero order, is subject to a prescribed temperature at the base and zero heat rejection at the tip. Approximate temperature distributions and companion heat transfer rates of excellent quality have been obtained by adequately blending a polynomial curve fit, the method of successive approximations and the method of undetermined coefficients. Detailed error distributions are also presented for real uniform circumferential fins using the exact solution by modified Bessel functions as the baseline case. The calculations of analytic character were carried out with a symbolic algebra software, Maple V, on a personal computer. Received on 23 October 1997     

On poro-hyperelastic shear

The paper examines the problem of the shear of a porous hyperelastic material, the pore space of which is saturated with an incompressible fluid. Poro-hyperelasticity provides a suitable approach for modelling the mechanical behaviour of highly deformable materials in engineering applications and particularly soft tissues encountered in biomechanical applications. Unlike with the infinitesimal theory of poroelasticity, the application of pure shear generates pore fluid pressures that dissipate with time as fluid migrates either from or into the pore space due to the generated fluid pressure gradients. The analytical results provide benchmark problems that can be used to examine the accuracy of computational approaches.     

CHAOTIC TRANSIENTS IN A CURVED FLUID CONVEYING TUBE

I. INTRODUCTIONChaotic motion is a kind of reciprocal non-periodic motion caused by a deterministic system. Itis very sensitive to the initial conditions, apparently random and incapable of long-term prediction.Chaotic transient is such a motion that the system will undergo a final steady state, which is one ofseveral or much more possible steady motions of the system. But this final steady state is very sensitiveto the initial conditions of the system. As it is e?ected by many uncertain…     

Unsteady MHD mixed convection flow over an impulsively stretched permeable vertical surface in a quiescent fluid

The unsteady mixed convection flow of an incompressible laminar electrically conducting fluid over an impulsively stretched permeable vertical surface in an unbounded quiescent fluid in the presence of a transverse magnetic field has been investigated. At the same time, the surface temperature is suddenly increased from the surrounding fluid temperature or a constant heat flux is suddenly imposed on the surface. The problem is formulated in such a way that for small time it is governed by Rayleigh type of equation and for large time by Crane type of equation. The non-linear coupled parabolic partial differential equations governing the unsteady mixed convection flow under boundary layer approximations have been solved analytically by using the homotopy analysis method as well as numerically by an implicit finite difference scheme. The local skin friction coefficient and the local Nusselt number are found to decrease rapidly with time in a small time interval and they tend to steady-state values for t*≥5. They also increase with the buoyancy force and suction, but decrease with injection rate. The local skin friction coefficient increases with the magnetic field, but the local Nusselt number decreases. There is a smooth transition from the unsteady state to the steady state.