Note on unsteady viscous flow on the outside of an expanding or contracting cylinder

In this paper, the viscous flow on the outside of an expanding or contracting cylinder is studied. The governing Navier-Stokes equations are transformed into a similarity equation, which is solved by a shooting method. The solution is an exact solution to the unsteady Navier-Stokes equations. Results show both trivial and non-trivial solutions. For trivial solutions, there is no axial flow induced during the cylinder expansion or contraction. However, for the non-trivial solutions which only exist for cylinder expansion, an axial flow is generated and its strength increases with the increase in expansion speed.     

The method of lines for the computation of incompressible viscous flows

A common approach to finding numerical solutions of the time-dependent incompressible Navier-Stokes equations is considered within the method of lines framework [1]. For the determination of viscous incompressible flows the stream-function formulation for the fourth-order equation [2, 3], an artificial compressibility method [4], and a modified velocity correction method [5] are designed. Some improved and extended results of numerical simulations obtained by the author in the previous works are presented. Test calculations have been done for various flows inside square, triangular and semicircular cavities with one moving wall, the backward-facing step, double bent channels and for the flow around an aerofoil at large angle of attack. An alternative and practical methodology for resolving the Navier-Stokes equations in arbitrarily complex geometries using Cartesian meshes is proposed. Some of complex geometrical configurations can be decomposed into a set of subdomains. The simplest approach for specifying boundary conditions near curved or irregular boundaries is to transfer all the variables from the boundaries to the nearest grid knots. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Navier－Stokes方程组的一个一致有效渐近解

本文利用多重尺度法[1,2]研究了大雷诺数情况下的平板绕流问题,得到了Navier-Stokes方程的一个一致有效渐近解。     

Dimension splitting method for the three dimensional rotating Navier-Stokes equations

In this paper, we propose a dimensional splitting method for the three dimensional (3D) rotating Navier-Stokes equations. Assume that the domain is a channel bounded by two surfaces and is decomposed by a series of surfaces ■i into several sub-domains, which are called the layers of the flow. Every interface i between two sub-domains shares the same geometry. After establishing a semi-geodesic coordinate (S-coordinate) system based on ■i , Navier-Stoke equations in this coordinate can be expressed as the sum of two operators, of which one is called the membrane operator defined on the tangent space on ■i , another one is called the bending operator taking value in the normal space on ■i . Then the derivatives of velocity with respect to the normal direction of the surface are approximated by the Euler central difference, and an approximate form of Navier-Stokes equations on the surface ■i is obtained, which is called the two-dimensional three-component (2D-3C) Navier-Stokes equations on a two dimensional manifold. Solving these equations by alternate iteration, an approximate solution to the original 3D Navier-Stokes equations is obtained. In addition, the proof of the existence of solutions to 2D-3C Navier-Stokes equations is provided, and some approximate methods for solving 2D-3C Navier-Stokes equations are presented.     

Unconditionally Optimal Error Estimates of the Bilinear-Constant Scheme for Time-Dependent Navier-Stokes Equations

In this paper,the unconditional error estimates are presented for the time-dependent Navier-Stokes equations by the bilinear-constant scheme.The corresponding optimal error estimates for the velocity and the pressure are derived unconditionally,while the previous works require certain time-step restrictions.The analysis is based on an iterated time-discrete system,with which the error function is split into a temporal error and a spatial error.The τ-independent(τ is the time stepsize)error estimate between the numerical solution and the solution of the time-discrete system is proven by a rigorous analysis,which implies that the numerical solution in L^∞-norm is bounded.Thus optimal error estimates can be obtained in a traditional way.Numerical results are provided to confirm the theoretical analysis.     

A locally conservative LDG method for the incompressible Navier-Stokes equations

In this paper a new local discontinuous Galerkin method for the incompressible stationary Navier-Stokes equations is proposed and analyzed. Four important features render this method unique: its stability, its local conservativity, its high-order accuracy, and the exact satisfaction of the incompressibility constraint. Although the method uses completely discontinuous approximations, a globally divergence-free approximate velocity in is obtained by simple, element-by-element post-processing. Optimal error estimates are proven and an iterative procedure used to compute the approximate solution is shown to converge. This procedure is nothing but a discrete version of the classical fixed point iteration used to obtain existence and uniqueness of solutions to the incompressible Navier-Stokes equations by solving a sequence of Oseen problems. Numerical results are shown which verify the theoretical rates of convergence. They also confirm the independence of the number of fixed point iterations with respect to the discretization parameters. Finally, they show that the method works well for a wide range of Reynolds numbers.

     

用SADI方法求解方形空腔中具有高Rayleigh数的非定常自然对流问题

本文用三次样条积分计算了在方形空腔中具有高Rayleigh数Ra=107和Ra=2×107的非定常自然对流问题。二维N-S方程和能量方程是在非均匀网格中用两个交替方向的三次样条公式进行计算的。文中简要讨论了过渡流动的主要特征,所得结果与理论予估值[1,2]吻合很好。Ra=107时的稳态结果与近期文献中的结果一致。     

非线性耦合Klein-Gordon方程组和耦合Schr\"{o}dinger-Boussinesq方程组的双行波解

本文提出了一种全新复合$(\frac{G''}{G})$展开方法,运用这种新方法并借助符号计算软件构造了非线性耦合Klein-Gordon方程组和耦合Schr\"{o}dinger-Boussinesq方程组的多种双行波解,包括双双曲正切函数解,双正切函数解,双有理函数解以及它们的混合解. 复合$(\frac{G''}{G})$展开方法不但直接有效地求出了两类非线性偏微分方程的双行波解,而且扩大了解的范围.这种新方法对于研究非线性偏微分方程具有广泛的应用意义.     

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

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

地下水在粘弹性含水层系中不定常渗流

本文研究了地下水在粘弹性含水层系中不定常渗流动态。在前人工作的基础上导出了新的微分-积分方程组。已知的微分方程组是它的特殊情况。新的线性微分-积分方程组描写了弱压缩流体在粘弹性含水层系中流动。用Laplace变换的方法求得了微分-积分方程组的解析解。粘弹性增加了含水层系的非均质性,即具有延迟和补给的性质。数值反演解和解析解符合得较好。它们给出了地下水在这种非均质含水层系中水位变化动态。     

A C0-WEAK GALERKIN FINITE ELEMENT METHOD FOR THE TWO-DIMENSIONAL NAVIER-STOKES EQUATIONS IN STREAM-FUNCTION FORMULATION

We propose and analyze a $C^0$-weak Galerkin (WG) finite element method for the numerical solution of the Navier-Stokes equations governing 2D stationary incompressible flows. Using a stream-function formulation, the system of Navier-Stokes equations is reduced to a single fourth-order nonlinear partial differential equation and the incompressibility constraint is automatically satisfied. The proposed method uses continuous piecewise-polynomial approximations of degree $k+2$ for the stream-function $\psi$ and discontinuous piecewise-polynomial approximations of degree $k+1$ for the trace of $\nabla\psi$ on the interelement boundaries. The existence of a discrete solution is proved by means of a topological degree argument, while the uniqueness is obtained under a data smallness condition. An optimal error estimate is obtained in $L^2$-norm, $H^1$-norm and broken $H^2$-norm. Numerical tests are presented to demonstrate the theoretical results.     

An application of a boundary element method to natural convection

The direct boundary element method is applied to the numerical modelling of thermal fluid flow in a transient state. The Navier-Stokes equations are considered under the Boussinesq approximation and the viscous thermal flow equations are expressed in terms of stream function, vorticity, and temperature in two dimensions. Boundary integral equations are derived using logarithmic potential and time-dependent heat potential as fundamental solutions. Boundary unknowns are discretized by linear boundary elements and flow domains are divided into a series of triangular cells. Charged points are translated upstream in the numerical evaluation of convective terms. Unknown stream function, vorticity, and temperature are staggered in the computational scheme.

Simple iteration is found to converge to the quasi steady-state flow. Boundary solutions for two-dimensional examples at a Reynolds number 100 and Grashoff number 10⁷ are obtained.     

Fujita-Kato Theorem for the Inhomogeneous Incompressible Navier-Stokes Equations with Nonnegative Density

In this paper, we prove the global existence and uniqueness of solutions for the inhomogeneous Navier-Stokes equations with the initial data $(\rho_0,u_0)\in L^∞\times H^s_0$, $s>\frac{1}{2}$ and $||u_0||_{H^s_0}\leq \varepsilon_0$ in bounded domain $\Omega \subset \mathbb{R}^3$, in which the density is assumed to be nonnegative. The regularity of initial data is weaker than the previous $(\rho_0,u_0)\in (W^{1,\gamma}∩L^∞)\times H^1_0$ in [13] and $(\rho_0,u_0)\in L^∞\times H^1_0$ in [7], which constitutes a positive answer to the question raised by Danchin and Mucha in [7]. The methods used in this paper are mainly the classical time weighted energy estimate and Lagrangian approach, and the continuity argument and shift of integrability method are applied to complete our proof.     

Stagnation flows with slip: Exact solutions of the Navier-Stokes Equations

Exact similarity solutions of the Navier-Stokes equations are found for stagnation flows towards a plate with slip. The solutions are applicable to the slip regime of rarefied gases.     

The oscillatory channel flow with arbitrary wall injection

In this article, we consider the laminar oscillatory flow in a low aspect ratio channel with porous walls. For small-amplitude pressure oscillations, we derive asymptotic formulations for the flow parameters using three different perturbation approaches. The undisturbed state is represented by an arbitrary mean-flow solution satisfying the Berman equation. For uniform wall injection, symmetric solutions are obtained for the temporal field from both the linearized vorticity and momentum transport equations. Asymptotic solutions that have dissimilar expressions are compared and shown to agree favourably with one another and with numerical experiments. In fact, numerical simulations of both linearly perturbed and nonlinear Navier-Stokes equations are used for validation purposes. As we insist on verifications, the absolute error associated with the total time-dependent velocities is analysed. The order of the cumulative error is established and the formulation based on the two-variable multiple-scale approach is found to be the most general and accurate. The explicit formulations help unveil interesting technical features and vortical structures associated with the oscillatory wave character. A similarity parameter is shown to exist in all formulations regardless of the mean-flow selection.     

Analytical solutions of a coupled nonlinear system arising in a flow between stretching disks

In this paper, we consider the axi-symmetric flow between two infinite stretching disks. By using a similarity transformation, we reduce the governing Navier-Stokes equations to a system of nonlinear ordinary differential equations. We first obtain analytical solutions via a four-term perturbation method for small and large values of the Reynolds number R. Also, we apply the Homotopy Analysis Method (which may be used for all values of R) to obtain analytical solutions. These solutions converge over a larger range of values of the Reynolds number than the perturbation solutions. Our results agree well with the numerical results of Fang and Zhang [22]. Furthermore, we obtain the analytical solutions valid for moderate values of R by use of Homotopy Analysis.     

Effects of Poiseuille flow on peristaltic transport of a particulate suspension

The problem of peristaltic transport induced by sinusoidal waves of a particle-fluid mixture in the presence of a Poiseuille flow, is analysed. The governing equations of motion resulting from the Navier-Stokes equations for both the fluid and particle phases are solved and closed form solutions are obtained for limiting values of Reynolds number, wave number and the Poiseuille flow parameter while the method of Frobenius series solution is used for the general case. It is found that the mean flow is strongly dependent on the Poiseuille flow parameter. The effects of particle concentration in the fluid is well discharged throughout the analysis and the results are compared with the other studies in the literature.     

Finite difference scheme based on proper orthogonal decomposition for the nonstationary Navier-Stokes equations

The proper orthogonal decomposition(POD)and the singular value decomposition(SVD) are used to study the finite difference scheme(FDS)for the nonstationary Navier-Stokes equations. Ensembles of data are compiled from the transient solutions computed from the discrete equation system derived from the FDS for the nonstationary Navier-Stokes equations.The optimal orthogonal bases are reconstructed by the elements of the ensemble with POD and SVD.Combining the above procedures with a Galerkin projection approach yields a new optimizing FDS model with lower dimensions and a high accuracy for the nonstationary Navier-Stokes equations.The errors between POD approximate solutions and FDS solutions are analyzed.It is shown by considering the results obtained for numerical simulations of cavity flows that the error between POD approximate solution and FDS solution is consistent with theoretical results.Moreover,it is also shown that this validates the feasibility and efficiency of POD method.     

Homotopy Analysis Method for Solving(2+1)-dimensional Navier-Stokes Equations with Perturbation Terms

In this paper Homotopy Analysis Method (HAM) is implemented for obtaining approximate solutions of (2+1)-dimensional Navier-Stokes equations with perturbation terms. The initial approximations are obtained using linear systems of the Navier-Stokes equations;by the iterations formula of HAM,the first approxima-tion solutions and the second approximation solutions are successively obtained and Homotopy Perturbation Method(HPM)is also used to solve these equations;finally, approximate solutions by HAM of (2+1)-dimensional Navier-Stokes equations with-out perturbation terms and with perturbation terms are compared. Because of the freedom of choice the auxiliary parameter of HAM,the results demonstrate that the rapid convergence and the high accuracy of the HAM in solving Navier-Stokes equa-tions;due to the effects of perturbation terms,the 3rd-order approximation solutions by HAM and HPM have great fluctuation.     

Bornes sur la densité pour les équations de Navier-Stokes compressibles isentropiques avec conditions aux limites de Dirichlet

We show that the local bounds on the density, obtained in a previous work by the author, for the solutions of compressible, isentropic Navier-Stokes equations with Dirichlet boundary conditions hold in fact up to the boundary.