求解二维不可压缩Navier-Stokes方程的混合型微分求积法

微分求积法（ＤＱＭ）能以较少的网格点求得微分方程的高精度数值解,但采用单纯的微分求积法求解二维不可压缩Ｎａｖｉｅｒ＿Ｓｔｏｋｅｓ 方程时,只能对低雷诺数流动获得较好的数值解,当雷诺数较高时会导致数值解不收敛·　为此,提出了一种微分求积法与迎风差分法混合求解二维不可压缩Ｎａｖｉｅｒ＿Ｓｔｏｋｅｓ 方程的预估＿校正数值格式,用伪时间相关算法以较少的网格点获得了较高雷诺数流动的数值解·　作为算例,对１∶１ 和１∶２ 驱动方腔内的流动进行了计算,得到了较好的数值结果·     

考虑到渗透效应的一种血液流动的计算方法

得到了定常情况下,狗二分叉动脉横截面的三维Navier＿Stokes方程的有限元处理方法,并考虑到管壁的渗透影响,数值方法还包括直角坐标和曲线坐标的变换·　详细讨论了渗透性影响下的定常流、分叉流以及切应力情况·　以分支和主干血管的速度比为参量,计算雷诺数为1000情况下管壁切应力,数值结果和先前的实验结果符合得很好·　该文的工作是Sharma等(2001)工作的改进,使计算量更小,能够处理的雷诺数范围更大·     

变截面圆形管道粘性流动的伽辽金-摄动杂交解

采用伽辽金-摄动杂交法来研究壁面是正弦形状的变截面圆形管道的粘性流动,从而避免了摄动小参数的局限性和单纯伽辽金法基函数选取的任意性的困难．讨论了边界和雷诺数对流动的影响,获得流动分离点和附着点的位置,还分析了壁面剪应力和摩擦系数沿轴向的变化情况．在小参数的情况下,计算所获得的结果与摄动解吻合良好．     

关于绕任意机翼非定常流动的一种无条件稳定的欧拉方程解

本文给出了绕二维与三维刚性或弹性振动机翼非定常无粘流动的欧拉方程解。首先利用Jameson的有限体积方法建立了求解欧拉方程的Runge-Kutta方法。为了提高受Runge-Kutta法稳定性限制的时间步长,文中采用了变系数的残值光顺方法。该方法避免了常系数残值光顺引起局部流场损失较大的问题。同时可在保证原计算格式精度要求下,大幅度提高计算时间步长,从而提高了计算效率。文中以二维与三维矩形机翼为例,分别对其在跨音速流场中作则性或弹性振动的非定常气动力进行了计算,研究了不同振动频率对流动产生的影响。部分计算结果与相应实验结果进行了比较。结果证明本方法是可靠的,可以用于求解绕任意运动机翼非定常流动问题。     

粘弹性二阶流体混合层流场拟序结构的数值研究

本文用拟谱方法对随时间发展的二维粘弹性二阶流体混合层流场进行了直接数据值模拟,给出在高雷诺数和低Deborah数下大涡的卷起、配对和合并等过程,通过与相同雷诺数下牛顿流体的比较,揭示了弱粘弹性对混合层中大涡拟序结构演变的影响.     

强迫对流条件下单颗粒水煤浆着火规律的研究

本文对强迫对流条件下单颗粒水煤浆的着火问题进行了实验研究及理论分析。实验表明,在雷诺数较大的情况下,火焰呈轴对称,着火发生于前半球的边界层内。本文建立了二维轴对称层流边界层着火的简化模型,并进行了数值计算。计算结果与实验结果符合较好,可作为在强迫对流条件下水煤浆滴着火温度与着火时间的估算方法。     

定常Stokes问题的边界积分方程法

1.前言 定常Stokes问题本身虽然只反映在小雷诺数情况下不可压缩粘性流体的定常流动,然而却为处理完整的Navier-Stokes方程奠定了基础. Stokes问题一般有两种公式化途径,一是通过流函数,二是利用速度-压力公式.两种公式化途径的区域类型数值方法,如有限差分法及有限单元法,已有不少工作,见[3]和[11].近年来,对这两种公式化途径的边界类型数值方法的研究,也获得一些结果.     

关于绕任意机翼非常流动的一种无条件稳定的欧拉方程解

本给出了绕二维与三维刚性或弹性振动机翼非定常无粘流动的欧拉方程解。首先利用Ｊａｍｅｓｏｎ的有限体积方法建立了求解欧拉方程的Ｒｕｎｇｅ－Ｋｕｔｔａ方法，为了提高受Ｒｕｎｇｅ－Ｋｕｌｔｔａ法稳定性限制的时间步长，中采用了变系数的残值光顺方法，该方法避免了常系数残值光顺引起局部流场损失较大的问题。同时可在保证原计算格式精度要求下，大幅度提高计算时间步长，从而提高了计算效率。中以二维与三维矩形机翼为     

刚性轴对称狭窄血管内压强及其梯度的研究

基于对Karman-Pohlhausen方法的改进,运用非线性多项式拟合和数值积分,导出了刚性轴对称狭窄管内压强及其梯度的轴向分布,讨论了该分布与雷诺数和狭窄管几何形状之间的关系．结果表明随狭窄度和雷诺数的增加,压强及其梯度在狭窄区域的振荡会迅速加剧,并逐渐导致舒张区出现负压强．尤其在狭窄的轴向区域变宽时,舒张区的血流状态会受到较大影响．在高雷诺数和重度狭窄时,理论计算与过去的实验结果基本一致．     

两同心旋转球间流动问题分歧解的连续算法

1引言本文数值地考察了两个同心旋转球之间的定常轴对称不可压流动．这种流动被称为球面Couette流（SphericalCouetteFlow），简称SCF．SCF对干天体物理，地球物理和工程应用均具重要意义，虽然过往所做的研究甚少（一方面由于分析研究的难度，另方面，所做的实验也少），目前，对其研究的兴趣有增长的趋势．实验发现，SCF在低雷诺数下既是轴对称的，又是关于赤道成反射对称的（Khlebutin［‘］，ZlereP＆SawatskiL‘）．ZiereP＆Sawatski和Wimmer［’］都发现SCF有临界雷诺数Rec，当Re＞Rec，有泰勒旋涡（TaylorVortices）形成…     

Nonlinear simulations of vesicle wrinkling

In this paper, we study nonlinear wrinkling dynamics of a vesicle in an extensional flow. Motivated by the recent experiments and linear theory on wrinkles of a quasi‐spherical membrane, we are interested in examining the linear theory and exploring wrinkling dynamics in a nonlinear regime. We focus on a quasi‐circular vesicle in two dimensions and show that the linear analytical results are qualitatively independent of the number of dimensions. Hence, the two‐dimensional studies can provide insights into the full three‐dimensional problem. We develop a spectral accurate boundary integral method to simulate the nonlinear evolution of surface tension and the nonlinear interactions between flow and membrane morphology. We demonstrate that for a quasi‐circular vesicle, the linear theory well predicts the characteristic wavenumber during the wrinkling dynamics. Nonlinear results of an elongated vesicle show that there exist dumbbell‐like stationary shapes in weak flows. For strong flows, wrinkles with pronounced amplitudes will form during the evolution. As far as the shape transition is concerned, our simulations are able to capture the main features of wrinkles observed in the experiments. Interestingly, numerical results reveal that, in addition to wrinkling, asymmetric rotation can occur for slightly tilted vesicles. The mathematical theory and numerical results are expected to lead to a better understanding of related problems in biology such as cell wrinkling. Copyright © 2013 John Wiley & Sons, Ltd.     

Two dimensional incompressible ideal flow around a thin obstacle tending to a curve

In this work we study the asymptotic behavior of solutions of the incompressible two dimensional Euler equations in the exterior of a single smooth obstacle when the obstacle becomes very thin tending to a curve. We extend results by Iftimie, Lopes Filho and Nussenzveig Lopes, obtained in the context of an obstacle tending to a point, see [D. Iftimie, M.C. Lopes Filho, H.J. Nussenzveig Lopes, Two dimensional incompressible ideal flow around a small obstacle, Comm. Partial Differential Equations 28 (1–2) (2003) 349–379].     

Coupling of compressible and incompressible flow regions using the multiple pressure variables approach

In many cases, multiphase flows are simulated on the basis of the incompressible Navier–Stokes equations. This assumption is valid as long as the density changes in the gas phase can be neglected. Yet, for certain technical applications such as fuel injection, this is no longer the case, and at least the gaseous phase has to be treated as a compressible fluid. In this paper, we consider the coupling of a compressible flow region to an incompressible one based on a splitting of the pressure into a thermodynamic and a hydrodynamic part. The compressible Euler equations are then connected to the Mach number zero limit equations in the other region. These limit equations can be solved analytically in one space dimension that allows to couple them to the solution of a half‐Riemann problem on the compressible side with the help of velocity and pressure jump conditions across the interface. At the interface location, the flux terms for the compressible flow solver are provided by the coupling algorithms. The coupling is demonstrated in a one‐dimensional framework by use of a discontinuous Galerkin scheme for compressible two‐phase flow with a sharp interface tracking via a ghost‐fluid type method. The coupling schemes are applied to two generic test cases. The computational results are compared with those obtained with the fully compressible two‐phase flow solver, where the Mach number zero limit is approached by a weakly compressible fluid. For all cases, we obtain a very good agreement between the coupling approaches and the fully compressible solver. Copyright © 2014 John Wiley & Sons, Ltd.     

Incipient separation over wall irregularities in transonic flow

Incipient separation over wall irregularities in a steady two dimensional flow field of a perfect fluid which has transonic speed characteristics has been investigated considering viscous-inviscid interactions at high Reynolds number. The aim of this work is to investigate dependence of the critical hump height (when a well attached flow over rigid body surface turns into a separated one) on the Karman–Guderley parameter which characterizes of the local flow field. The analysis of the flow field starts with the so-called inspection analysis of the flow properties and then the interaction problem has been constructed using the asymptotic analysis of triple-deck structure of interaction region. Finally, a method based on a semi-direct solution of governing equations of the transonic interaction problem has been used to obtain the numerical solution of the problem.     

Three-dimensional Couette flow of a dusty fluid with heat transfer

This work is focused on the mathematical modeling of three-dimensional Couette flow and heat transfer of a dusty fluid between two infinite horizontal parallel porous flat plates. The problem is formulated using a continuum two-phase model and the resulting equations are solved analytically. The lower plate is stationary while the upper plate is undergoing uniform motion in its plane. These plates are, respectively, subjected to transverse exponential injection and its corresponding removal by constant suction. Due to this type of injection velocity, the flow becomes three dimensional. The closed-form expressions for velocity and temperature fields of both the fluid and dust phases are obtained by solving the governing partial differential equations using the perturbation method. A selective set of graphical results is presented and discussed to show interesting features of the problem.     

Numerical simulation of generalized newtonian blood flow past a couple of irregular arterial stenoses

This paper looked at the numerical investigations of the generalized Newtonian blood flow through a couple of irregular arterial stenoses. The flow is treated to be axisymmetric, with an outline of the stenoses obtained from a three dimensional casting of a mild stenosed artery, so that the flow effectively becomes two‐dimensional. The Marker and Cell (MAC) method is developed for the governing unsteady generalized Newtonian equations in staggered grid for viscous incompressible flow in the cylindrical polar co‐ordinates system. The derived pressure‐Poisson equation was solved using Successive‐Over‐Relaxation (S.O.R.) method and the pressure‐velocity correction formulae have been derived. Computations are performed for the pressure drop, the wall shear stress distribution and the separation region. The presented computations show that in comparison to the corresponding Newtonian model the generalized Newtonian fluid experiences higher pressure drop, lower peak wall shear stress and smaller separation region. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 960–981, 2011     

Advanced continuum modelling of gas-particle flows beyond the hydrodynamic limit

The accurate prediction of dilute gas-particle flows using Euler–Euler models is challenging because particle–particle collisions are usually not dominant in such flows. In other words, in dilute flows the particle Knudsen number is not small enough to justify a Chapman–Enskog expansion about the collision-dominated near-equilibrium limit. Moreover, due to the fluid drag and inelastic collisions, the granular temperature in gas-particle flows is often small compared to the mean particle kinetic energy, implying that the particle-phase Mach number can be very large. In analogy to rarefied gas flows, it is thus not surprising that two-fluid models fail for gas-particle flows with moderate Knudsen and Mach numbers. In this work, a third-order quadrature-based moment method, valid for arbitrary Knudsen number, coupled with a fluid solver has been applied to simulate dilute gas-particle flow in a vertical channel with particle-phase volume fractions between 0.0001 and 0.01. In order to isolate the instabilities that arise due to fluid-particle coupling, a fluid mass flow rate that ensures that turbulence would not develop in a single phase flow (Re = 1380) is employed. Results are compared with the predictions of a two-fluid model with standard kinetic theory based closures for the particle phase. The effect of the particle-phase volume fraction on flow instabilities leading to particle segregation is investigated, and differences with respect to the two-fluid model predictions are examined. The influence of the discretization on the solution of both models is investigated using three different grid resolutions. Radial profiles of phase velocities and particle concentration are shown for the case with an average particle volume fraction of 0.01, showing the flow is in the core-annular regime.     

Analytical solutions of laminar swirl decay in a straight pipe

In this work, the laminar swirl flow in a straight pipe is revisited and solved analytically by using prescribed axial flow velocity profiles. Based on two axial velocity profiles, namely a slug flow and a developed parabolic velocity profiles, the swirl velocity equation is solved by the separation of variable technique for a rather general inlet swirl velocity distribution, which includes a forced vortex in the core and a free vortex near the wall. The solutions are expressed by the Bessel function for the slug flow and by the generalized Laguerre function for the developed parabolic velocity. Numerical examples are calculated and plotted for different combinations of influential parameters. The effects of the Reynolds number, the pipe axial distance, and the inlet swirl profiles on the swirl velocity distribution and the swirl decay are analyzed. The current results offer analytical equations to estimate the decay rate and the outlet swirl intensity and velocity distribution for the design of swirl flow devices.     

Numerical simulation of laminar flow past a circular cylinder

The present paper focuses on the analysis of two- and three-dimensional flow past a circular cylinder in different laminar flow regimes. In this simulation, an implicit pressure-based finite volume method is used for time-accurate computation of incompressible flow using second order accurate convective flux discretisation schemes. The computation results are validated against measurement data for mean surface pressure, skin friction coefficients, the size and strength of the recirculating wake for the steady flow regime and also for the Strouhal frequency of vortex shedding and the mean and RMS amplitude of the fluctuating aerodynamic coefficients for the unsteady periodic flow regime. The complex three dimensional flow structure of the cylinder wake is also reasonably captured by the present prediction procedure.