Vorticity–velocity formulation of the 3D Navier–Stokes equations in cylindrical co‐ordinates

A finite difference method is presented for solving the 3D Navier–Stokes equations in vorticity–velocity form. The method involves solving the vorticity transport equations in ‘curl‐form’ along with a set of Cauchy–Riemann type equations for the velocity. The equations are formulated in cylindrical co‐ordinates and discretized using a staggered grid arrangement. The discretized Cauchy–Riemann type equations are overdetermined and their solution is accomplished by employing a conjugate gradient method on the normal equations. The vorticity transport equations are solved in time using a semi‐implicit Crank–Nicolson/Adams–Bashforth scheme combined with a second‐order accurate spatial discretization scheme. Special emphasis is put on the treatment of the polar singularity. Numerical results of axisymmetric as well as non‐axisymmetric flows in a pipe and in a closed cylinder are presented. Comparison with measurements are carried out for the axisymmetric flow cases. Copyright © 2003 John Wiley & Sons, Ltd.     

硬脑膜受压膨出挠度的力学分析

硬脑膜是一种粘弹性材料,为控制硬脑膜在脑压作用下的膨出度,对粘弹性薄膜受压膨出挠度作力学分析。以位移为未知量,从粘弹性材料的分型本构关系出发将Foepple薄膜大挠度理论从弹性推广到粘弹性膜,得到一组非线性积分偏微分方程。先在空间上运用Galerkin方法将积分偏微分方程组化为积分常微分方程组。然后,在时间域上运用数值积分和有限差分将方程离散为非线性代数方程组。本文对四周固定夹紧的圆形、椭圆形和矩形薄膜进行了求解,并将求解结果用于颅底缺损重建膜的膨出量计算,计算值与实验值吻合,为颅底外科提供一个理论分析方法。     

Periodic Traveling Waves for Diffusion Equations with Time Delayed and Non-local Responding Reaction

We develop a singular perturbation technique to study the existence of periodic traveling wave solutions with large wave speed for a class of reaction-diffusion equations with time delay and non-local response. Unlike the classical singular perturbation method, our approach is based on a transformation of the differential equations to integral equations in a Banach space that reduces the singular perturbation problem to a regular perturbation problem. The periodic traveling wave solutions then are obtained by the use of Liapunov-Schmidt method and a generalized implicit function theorem. The general result obtained has been applied to a non-local reaction-diffusion equation derived from an age-structured population model with a logistic type of birth function.     

THE DYNAMIC BUCKLING OF AN ORTHOTROPIC CYLINDRICAL THIN SHELL UNDER TORSION VARYING AS A POWER FUNCTION OF TIME

The subject of this investigation is to study the buckling of orthotropic cylindricalthin shells under torsion,which is a power function of time.The dynamic stability and compati-bility equations are obtained first.These equations are subsequently reduced to a time dependentdifferential equation with variable coefficient by using Galerkin's method.Finally,the critical dy-namic and static loading,the corresponding wave numbers,the dynamic factors,critical time andcritical impulse are found analytically by applying the Ritz type variational method.Using thoseresults,the effects of the variations of the power of time in the torsion load expression,of theloading parameter,the ratio of the Young's moduli and the ratio of the radius to thickness onthe critical parameters are studied numerically.It is observed that these factors have appreciableeffects on the critical parameters of the problem in the heading.     

Hydrodynamics with quadratic pressure. 1. General results

A wide class of solutions of Euler equations with quadratic pressure are described. In Lagrangian coordinates, these solutions linearize exactly momentum equations and are characterized by special initial data: the Jacobian matrix of the initial velocity field has constant algebraic invariants. The equations are integrated using the method of separation of the time and Lagrangian coordinates. Time evolution is defined by elliptic functions. The solutions have a poletype singularity at a finite time. A representation for the velocity vortex is given.     

卷积型GD半解析法及矩形薄板瞬态响应解

卷积型的Gurtin变分原理是目前在数学上唯一能和动力学初值问题完全等价的变分原理,它完全反映了有关初值问题的全部特征。GD法（General Differential Method）是从泰勒展开式出发,推出的一种求解偏微分方程的数值方法,本文系统地介绍了GD法的基本原理,以及权系数的推导。本文通过卷积将矩形薄板原始控...     

Computational Investigation of Wing Flutter

Numerical aeroelastic simulation of a high-aspect-ratio transport type wing model in transonic region is presented. The aeroelastic responses of the wing are extracted by integrating compressible thin-layer Navier-Stokes equations coupled with the equations of motion of the wing structure, in a time dependent manner. The Yee-Harten implicit TVD scheme and the Wilson's θ method are employed to integrate these equations, respectively. Flutter boundaries were found for Mach number range, 0.7 to 0.85 and the results were compared with experimental flutter boundaries. Futhermore, Limit Cycle Oscillations were found and the characteristics of the flutter and limit cycle oscillations are investigated and discussed.     

A tensor method for the derivation of the equations of rigid body dynamics

A tensor method for the derivation of the equations of rigid body dynamics,based onthe concepts of continuum mechanics,is presented.The formula of time derivative of theinertia tensor with zero corotational rate is used to prove the equivalences of five methods,namely,Lagrange’s equations,Nielsen’s equations,Gibbs-Appell’s equations,Kane’sequations and the generalized momentum type of Kane’s equations.Some differentialidentities on angular velocity and angular acceleration are given.     

Boundary element solution of unsteady magnetohydrodynamic duct flow with differential quadrature time integration scheme

A numerical scheme which is a combination of the dual reciprocity boundary element method (DRBEM) and the differential quadrature method (DQM), is proposed for the solution of unsteady magnetohydrodynamic (MHD) flow problem in a rectangular duct with insulating walls. The coupled MHD equations in velocity and induced magnetic field are transformed first into the decoupled time‐dependent convection–diffusion‐type equations. These equations are solved by using DRBEM which treats the time and the space derivatives as nonhomogeneity and then by using DQM for the resulting system of initial value problems. The resulting linear system of equations is overdetermined due to the imposition of both boundary and initial conditions. Employing the least square method to this system the solution is obtained directly at any time level without the need of step‐by‐step computation with respect to time. Computations have been carried out for moderate values of Hartmann number (M?50) at transient and the steady‐state levels. As M increases boundary layers are formed for both the velocity and the induced magnetic field and the velocity becomes uniform at the centre of the duct. Also, the higher the value of M is the smaller the value of time for reaching steady‐state solution. Copyright © 2005 John Wiley & Sons, Ltd.     

Study of differential control method for solving chaotic solutions of nonlinear dynamic system

In this paper, we focus on the need to solve chaotic solutions of high-dimensional nonlinear dynamic systems of which the analytic solution is difficult to obtain. For this purpose, a Differential Control Method (DCM) is proposed based on the Mechanized Mathematics-Wu Elimination Method (WEM). By sampling, the computer time of the differential operator of the nonlinear differential equation can be substituted by the differential quotient of solving the variable time of the sample. The main emphasis of DCM is placed on substituting the differential quotient of a small neighborhood of the sample time of the computer system for the differential operator of the equations studied. The approximate analytical chaotic solutions of the nonlinear differential equations governing the high-dimensional dynamic system can be obtained by the method proposed. In order to increase the computational efficiency of the method proposed, a thermodynamics modeling method is used to decouple the variable and reduce the dimension of the system studied. The validity of the method proposed for obtaining approximate analytical chaotic solutions of the nonlinear differential equations is illustrated by the example of a turbo-generator system. This work is applied to solving a type of nonlinear system of which the dynamic behaviors can be described by nonlinear differential equations.     

An immersed boundary vector potential-vorticity meshless solver of the incompressible Navier–Stokes equation

We present a strong form meshless solver for numerical solution of the nonstationary, incompressible, viscous Navier–Stokes equations in two (2D) and three dimensions (3D). We solve the flow equations in their stream function-vorticity (in 2D) and vector potential-vorticity (in 3D) formulation, by extending to 3D flows the boundary condition-enforced immersed boundary method, originally introduced in the literature for 2D problems. We use a Cartesian grid, uniform or locally refined, to discretize the spatial domain. We apply an explicit time integration scheme to update the transient vorticity equations, and we solve the Poisson type equation for the stream function or vector potential field using the meshless point collocation method. Spatial derivatives of the unknown field functions are computed using the discretization-corrected particle strength exchange method. We verify the accuracy of the proposed numerical scheme through commonly used benchmark and example problems. Excellent agreement with the data from the literature was achieved. The proposed method was shown to be very efficient, having relatively large critical time steps.     

快速求解叶栅流动的边界元方法

本文利用边界元方法，通过ｋｒｉｃｈｈｏｆｆ变换将描述叶栅流动的控制方程转换成线性方程。并将广义ｋ－Ｊ条件与边界积分方程联立求解，避免了非线性项和叶片出气角的迭代计算。完成了一种快速求解任意迥转面叶栅流场的计算程序，实用表明与其它数值方法及实验结果符合较好，具有快带、简明、实用的特点。     

The relation between CDM instability and Deborah number in differential type rheological equations

Previous studies have argued that rheological equations of the differential type, such as second-order fluid models, are inadequate because they result in unstable solution after cessation of steady shear. If the sign of the viscoelastic coefficient is selected so that the storage modulus is positive, the fluid velocity increases indefinitely and the flow does not decay by viscous dissipation, in contradiction to thermodynamic laws. This study mitigates this problem by demonstrating that the solution of such equations is actually stable at low values of Deborah number De, where these equations are only valid for other reasons. In fact, second order and higher order differential type equations are applicable only if the relaxation time of the fluid is low relative to a characteristic time of the flow. The study shows how to determine the characteristic time and thus clarifies, in practical terms, the limits of the region where differential type equations can be applied.Presented at the Third European Rheology Conference and Golden Jubilee Meeting of the British Society of Rheology, Edinburgh, Sept. 3–7, 1990.     

The DRBEM solution of incompressible MHD flow equations

This paper presents a dual reciprocity boundary element method (DRBEM) formulation coupled with an implicit backward difference time integration scheme for the solution of the incompressible magnetohydrodynamic (MHD) flow equations. The governing equations are the coupled system of Navier‐Stokes equations and Maxwell's equations of electromagnetics through Ohm's law. We are concerned with a stream function‐vorticity‐magnetic induction‐current density formulation of the full MHD equations in 2D. The stream function and magnetic induction equations which are poisson‐type, are solved by using DRBEM with the fundamental solution of Laplace equation. In the DRBEM solution of the time‐dependent vorticity and current density equations all the terms apart from the Laplace term are treated as nonhomogeneities. The time derivatives are approximated by an implicit backward difference whereas the convective terms are approximated by radial basis functions. The applications are given for the MHD flow, in a square cavity and in a backward‐facing step. The numerical results for the square cavity problem in the presence of a magnetic field are visualized for several values of Reynolds, Hartmann and magnetic Reynolds numbers. The effect of each parameter is analyzed with the graphs presented in terms of stream function, vorticity, current density and magnetic induction contours. Then, we provide the solution of the step flow problem in terms of velocity field, vorticity, current density and magnetic field for increasing values of Hartmann number. Copyright © 2010 John Wiley & Sons, Ltd.     

Exact solutions for axisymmetric flexural free vibrations of inhomogeneous circular Mindlin plates with variable thickness

Circular plates with radially varying thickness, stiffness, and density are widely used for the structural optimization in engineering. The axisymmetric flexural free vibration of such plates, governed by coupled differential equations with variable coefficients by use of the Mindlin plate theory, is very difficult to be studied analytically.In this paper, a novel analytical method is proposed to reduce such governing equations for circular plates to a pair of uncoupled and easily solvable differential equations of the Sturm-Liouville type. There are two important parameters in the reduced equations.One describes the radial variations of the translational inertia and flexural rigidity with the consideration of the effect of Poisson's ratio. The other reflects the comprehensive effect of the rotatory inertia and shear deformation. The Heun-type equations, recently well-known in physics, are introduced here to solve the flexural free vibration of circular plates analytically, and two basic differential formulae for the local Heun-type functions are discovered for the first time, which will be of great value in enriching the theory of Heun-type differential equations.     

Solution of the discretized incompressible Navier-Stokes equations with the GMRES method

We describe some experiences using interative solution methods of GMRES type to solve the discretized Navier-Stokes equations. The discretization combined with a pressure correction scheme leads to two different systems of equations: the momentum equations and the pressure equation. It appears that a fast solution method for the pressure equation is obtained by applying the recently proposed GMRESR method, or GMRES combined with a MILU preconditioner. The diagonally scaled momentum equations are solved by GMRES(m), a restarted version of GMRES.     

A control-volume based finite element method for three-dimensional incompressible turbulent fluid flow,heat transfer,and related phenomena

A control-volume based finite element method of equal-order type for three-dimensional incompressible turbulent fluid flow, heat transfer, and related phenomena is presented. The discretization equations are based mainly on the physics of the phenomena under consideration, more than on mathematical arguments. Special emphasis is devoted to the discretization of the convective terms and the continuity equation, and to the treatment of the boundary conditions imposed by the use of a high Reynolds k-?, type turbulence model. The pressure-velocity coupling in the fluid flow calculation is made from a derivative of the original SIMPLER method, without pressure correction. The discretized equations are solved in a sequential, rather than a coupled, form with significant advantage in the required computer time and storage. The method is an extension of a former version proposed by us for two-dimensional, laminar problems, and is here successfully applied to the following situations: three-dimensional deflected turbulent jet, and flows in 90° and 45° junctions of ducts with rectangular cross sections. The calculated results are in very good agreement with the experimental and numerical (obtained with the well established finite difference method) data available in the literature.     

不可压粘流N－S方程的边界积分解法

对原变量的Ｎ－Ｓ方程进行一阶时间离散，采用共轭梯度法解除压强－速度的耦合．对所得的一系列Ｌａｐｌａｃｅ方程、Ｐｏｓｓｉｏｎ方程和Ｈｅｌｍｈｏｔｚ方程均进行边界积分法求解，首次得到了粘性Ｎ－Ｓ方程的边界积分表示式．圆柱的定常、非定常尾迹计算结果表明了本文方法的有效性．     

NEW THEORY FOR EQUATIONS OF NON-FUGHSIAN TYPE REPRESENTATION THEOREM OF TREE SERIES SOLUTION (Ⅰ)

In the analytic theory of differential equations the exact explicit analytic solution has not been obtained for equations of the non-Fuchsian type (Poincare's problem). The new theory proposed in this paper for the first time affords a general method of finding exact analytic expres-sion for irregular integrals.By discarding the assumption of formal solution of classical theory,our method consists in deriving a cor-respondence relation from the equation itself and providing the analytic structure of irregular integrals naturally by the residue theorem. Irregular integrals are made up of three parts: noncontracted part,represented by ordinary recursion series,all-and semi-contracted part by the so-called tree series. Tree series solutions belong to analytic function of the new kind with recursion series as the special case only.