Some numerical experiments on the splitting of Burgers' equation

The combined approach of linearization and splitting up is used for devising new algorithms to solve a one-dimensional Burgers' equation. Two schemes are discussed and the computed solutions are compared with the exact solution. For this problem it is found that the schemes proposed yield excellent numerical results for Reynolds number R ranges from 50 up to 1500. The schemes were also tested for another problem whose R = 10000. In this case a filtering technique is used to overcome the nonlinear instability.     

Viscoelastic flow in a curved channel: A similarity solution for the Oldroyd-B fluid

A similarity solution is used to analyse the flow of the Oldroyd fluid B, which includes the Newtonian and Maxwell fluids, in a curved channel modelled by the narrow annular region between two circular concentric cylinders of large radius. The solution is exact, including inertial forces. It is found that the non-Netonian kinematics are very similar to the Newtonian ones, although some stress components can become very large. At high Reynolds number a boundary layer is developed at the inner cylinder. The structure of this boundary layer is asymptotically analysed for the Newtonian fluid. Non-Newtonian stress boundary layers are also developed at the inner cylinder at large Reynolds numbers.     

MHD turbulence in a channel with spanwise field

We study turbulent channel flow of an electrically conducting liquid with a homogeneous magnetic field imposed in the spanwise direction. The Lorentz force is modelled using the quasistatic approximation. Direct and large–eddy simulations are performed for hydrodynamic Reynolds numbers Re=10000 and Re=20000 and the Hartmann number varying in a wide range. The main effect of the magnetic field is the suppression of turbulent velocity fluctuations and momentum transfer in the wall–normal direction. Comparing the results from direct and large–edddy simulations we show that the dynamic Smagorinsky model accurately reproduces the flow transformation. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Weak solution of incompressible Euler equations with decreasing energy

Weak solution of incompressible Euler equations are L²-vector fields, satisfying integral relations, which express the mass and momentum balance. They are believed to describe the turbulent fluid motion at high Reynolds numbers. We justify this conjecture by constructing a weak solution with decreasing kinetic energy. The construction is based on Generalized Flows, introduced by Y. Brenier.     

Flow separation behind ellipses at Reynolds numbers less than 10

Flow separation behind two-dimensional ellipses with aspect ratios ranging from 0, a flat plate, to 1, a circular cylinder, were investigated for Reynolds numbers less than 10 using both a cellular automata model and a commercial computational fluid dynamics software program. The relationship between the critical aspect ratio for flow separation and Reynolds number was determined to be linear for Reynolds numbers greater than one. At slower velocities, the critical aspect ratio decreases more quickly as the Reynolds number approaches zero. The critical Reynolds numbers estimated for flow separation behind a flat plate and circular cylinder agree with extrapolations from experimental observations. Fluctuations in the values of the stream function for laminar flow behind the ellipses were found at combinations of Reynolds number and aspect ratio near the critical values for separation.     

On the Lagerstrom model for flow at low Reynolds numbers

We show that a model proposed by Lagerstrom for viscous incompressible flow at low Reynolds numbers has a solution. The proof is constructive in that we obtain lower and upper bounds for the solution. The dependence of the initial slope of the solution on a small parameter is also made explicit.     

Simulation of Separating Flows over a Backward-Facing Step

Simulation results are reported for plane two-dimensional viscous incompressible flow in a channel with an abrupt expansion. The mathematical model is provided by the quasi-hydrodynamic equations in the incompressible fluid approximation. The computations are carried out in a range of Reynolds numbers including both laminar and turbulent flow. As the Reynolds number increases, the solution bifurcates and the steady laminar flow changes to time-dependent flow. The computation results are consistent with known experimental data. Turbulence models were not used for large Reynolds number computations.     

Nonstationary mass transfer of a reacting spherical particle in laminar stream of a viscous fluid : PMM vol. 38, n 6, 1974, pp. 1041–1047

The process of the formation of a stationary mass transfer mode for a moving reacting particle is examined. An analytic expression valid for a nonstationary distribution of the concentration of matter in a steady stream of viscous fluid, flowing past a spherical particle, was obtained for the case when at a certain instant a chemical reaction of the first order begins at the surface of the sphere. The problem is solved for small finite Reynolds and Péclet numbers. The solution of the corresponding stationary problem has been obtained in [1]. Paper [2] examined a nonstationary heat transfer of a fluid spherical drop in an inviscid flow with spasmodic change of initial temperature at high Péclet numbers. Paper [3] contains an analysis of the problem of a nonstationary heat transfer of a rigid spherical particle for small Reynolds and Péclet numbers at spasmodic change of temperature of the particle surface. The results obtained in [3] can be used to describe the mass transfer for a moving reacting particle only in the case of a diffusion mode of the chemical reaction.     

Initial-Boundary Value Problem for the Reynolds Lubrication Equation

The nonstationary Reynolds equation of the theory of gas lubrication is considered. The existence and uniqueness of a solution to the initial-boundary value problem for this equation are established in the case of sufficiently smooth data. Estimates for the solution are obtained for large bearing numbers. Bibliography: 8 titles.     

Solvability of the plane problem of thermocapillary convection

For small Reynolds numbers one proves the unique solvability of the problem describing the stationary thermocapillary convection of a viscous incompressible fluid in a closed vessel in the plane. The solution is found in a class of functions belonging to a weighted Hölder space.Translated from Problemy Matematicheskogo Analiza, No. 10, pp. 33–47, 1986.In conclusion, the author expresses her gratitude to her scientific advisor V. A. Solonnikov for his considerable help.     

Integral Equation Formulation of Oseen Flow Problems

An integral equation approach to the problem of Oseen flow pastobstacles is presented and an approximate solution of the appropriateintegral equations is obtained for low Reynolds numbers. Bothobstacles in bounded and unbounded media are examined.     

Continuum models in problems of the hypersonic flow of a rarefied gas around blunt bodiest

All possible continuum (hydrodynamic) models in the case of two-dimensional problems of supersonic and hypersonic flows around blunt bodies in the two-layer model (a viscous shock layer and shock-wave structure) over the whole range of Reynolds numbers, Re, from low values (free molecular and transitional flow conditions) up to high values (flow conditions with a thin leading shock wave, a boundary layer and an external inviscid flow in the shock layer) are obtained from the Navier-Stokes equations using an asymptotic analysis. In the case of low Reynolds numbers, the shock layer is considered but the structure of the shock wave is ignored. Together with the well-known models (a boundary layer, a viscous shock layer, a thin viscous shock layer, parabolized Navier-Stokes equations (the single-layer model) for high, moderate and low Re numbers, respectively), a new hydrodynamic model, which follows from the Navier-Stokes equations and reduces to the solution of the simplified (“local”) Stokes equations in a shock layer with vanishing inertial and pressure forces and boundary conditions on the unspecified free boundary (the shock wave) is found at Reynolds numbers, and a density ratio, k, up to and immediately after the leading shock wave, which tend to zero subject to the condition that (k/Re)^1/2 → 0. Unlike in all the models which have been mentioned above, the solution of the problem of the flow around a body in this model gives the free molecular limit for the coefficients of friction, heat transfer and pressure. In particular, the Newtonian limit for the drag is thereby rigorously obtained from the Navier-Stokes equations. At the same time, the Knudsen number, which is governed by the thickness of the shock layer, which vanishes in this model, tends to zero, that is, the conditions for a continuum treatment are satisfied. The structure of the shock wave can be determined both using continuum as well as kinetic models after obtaining the solution in the viscous shock layer for the weak physicochemical processes in the shock wave structure itself. Otherwise, the problem of the shock wave structure and the equations of the viscous shock layer must be jointly solved. The equations for all the continuum models are written in Dorodnitsyn--Lees boundary layer variables, which enables one, prior to solving the problem, to obtain an approximate estimate of second-order effects in boundary-layer theory as a function of Re and the parameter k and to represent all the aerodynamic and thermal characteristic; in the form of a single dependence on Re over the whole range of its variation from zero to infinity.

An efficient numerical method of global iterations, previously developed for solving viscous shock-layer equations, can be used to solve problems of supersonic and hypersonic flows around the windward side of blunt bodies using a single hydrodynamic model of a viscous shock layer for all Re numbers, subject to the condition that the limit (k/Re)^1/2 → 0 is satisfied in the case of small Re numbers. An aerodynamic and thermal calculation using different hydrodynamic models, corresponding to different ranges of variation Re (different types of flow) can thereby, in fact, be replaced by a single calculation using one model for the whole of the trajectory for the descent (entry) of space vehicles and natural cosmic bodies (meteoroids) into the atmosphere.     

A segregated finite element formulation of Navier-Stokes equations under laminar conditions

This paper describes the application of the Taylor-Galerkin based finite element method to the calculation of incompressible viscous flows. The fundamental concepts and characteristics of the formulations and the associated solution methodology used are described in technical detail. This scheme is second-order accurate for velocity and pressure in time as well as in space. The three test cases chosen for the performance evaluation of this formulation are the standing vortex problem, the lid-driven cavity flow at Reynolds numbers up to 5000, and the sudden expansion problem at Reynolds numbers up to 800. These results are provided to illustrate the quality of the computed transient solutions in two dimensions under laminar conditions.     

Conforming Chebyshev Spectral Collocation Methods for the Solution of Laminar flow in a Constricted Channel

The numerical simulation of steady planar two-dimensional laminarflow of an incompressible fluid through an abruptly contractingchannel using spectral domain decomposition methods is described.The key features of the method are the decomposition of theflow region into a number of rectangular subregions and spectralapproximations that are pointwise C¹ continuous across subregioninterfaces. Spectral approximations to the solution are obtainedfor Reynolds numbers in the range [0, 500]. The size of thesalient corner vortex decreases as the Reynolds number increasesfrom 0 to around 45. As the Reynolds number is increased further,the vortex grows slowly. A vortex is detected downstream ofthe contraction at a Reynolds number of around 175 that continuesto grow as the Reynolds number is increased further.     

Rotation-Symmetric Referencebodys For Energy Efficient Flow Around Bodies

Topology optimization is used to optimize problems of flow around bodies and problems of guided flow. Within the context of research, optimization criteria are developed to increase the energy efficiency of these problems [1] [2] [3] [4] [5]. In order to evaluate the new criteria in respect to the increasing of energy efficiency, reference bodies for different Reynolds numbers in combination with given design space limitations are needed. Therefore, an optimal body at Reynolds number against 0 was analytically determined by Bourot [6]. At higher Reynolds numbers, in the range of laminar and turbulent flows, no analytical solution is known. Accordingly, reference bodies are calculated by CFD calculations at three technical relevant Reynolds numbers (1.000, 32.000, 100.000) in combination with parameter optimization. The cross section of the bodies is described by a parameterized model. To get the optimal body, a parameter optimization based on a “brute force”; algorithm is used to optimize with regard to the friction loss and pressure loss in order to minimize the total loss (c_d-value). The result is an optimal parameter constellation, depending on the Reynolds number. Within the results, it is possible to develop the optimal geometries. The identified characteristics of the flow field around these bodies are used as base for new optimization criteria. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical study of pressure distribution in entrance pipe flow

This article describes the computation of pipe flow in the entrance region. The pressure distribution and flow characteristics, particularly the effect of vorticity in the vicinity of the wall, were analyzed for moderate Reynolds numbers (Re) ranging from 500 to 10,000. It was found, for the first time, that a large pressure gradient in the radial direction exists near the pipe inlet. The pressure gradient is caused by the radial component of the curl of vorticity, which decreases as Re increases. The pressure at the wall is lower than that at the central core for Re ≤$\le$ 5000. This result is beyond the scope of the boundary-layer assumption for pressure, although it applies to flows at high Reynolds numbers.     

A second-order finite element variational multiscale scheme for the fully discrete unsteady Navier–Stokes equations

In this report, we present and study a fully discrete finite element variational multiscale scheme for the unsteady incompressible Navier–Stokes equations where high Reynolds numbers are allowed. The scheme uses conforming finite element pairs for spatial discretization and a three-point difference formula for temporal discretization which is of second-order, where a stabilization term based on two local Gauss integrations is employed to stabilize the numerical scheme. We prove stability of the scheme, derive a priori error estimates for the fully discrete solution, and finally, give some numerical results to verify the theoretical predictions and demonstrate the effectiveness of the proposed numerical scheme.     

On singular limit for the compressible Navier–Stokes system with non-monotone pressure law

We consider a singular limit for the compressible Navier–Stokes system with general non-monotone pressure law in the asymptotic regime of low Mach number and large Reynolds numbers. We show that any dissipative weak solution approaches the solution of incompressible Euler equation both for well-prepared initial data and ill-prepared initial data.     

General solutions of the Stokes' flow equations

The Papkovich-Neuber solution of the Stokes' flow equations for a viscous incompressible flow at very low Reynolds numbers is derived using three alternative approaches. In each case the class of solutions which is valid varies depending on the assumptions made on the domain considered. The technique is illustrated by several examples. The solution is identical to that of incompressible linear elasticity upon appropriate modification.     

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

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