Numerical implementations of an iterative method with boundary condition splitting as applied to the nonstationary stokes problem in the gap between coaxial cylinders

Numerical implementations of a new fast-converging iterative method with boundary condition splitting are constructed for solving the Dirichlet initial-boundary value problem for the nonstationary Stokes system in the gap between two coaxial cylinders. The problem is assumed to be axially symmetric and periodic along the cylinders. The construction is based on finite-difference approximations in time and bilinear finite-element approximations in a cylindrical coordinate system. A numerical study has revealed that the iterative methods constructed have fairly high convergence rates that do not degrade with decreasing viscosity (the error is reduced by approximately 7 times per iteration step). Moreover, the methods are second-order accurate with respect to the mesh size in the max norm for both velocity and pressure.     

On the equilibrium of circular flows under gravity

This paper investigates the stability of an incompressible inviscid rotatory couette flow confined between two corotating and coaxial vertical cylinders under the force field of gravity. The stationary distributions of density and pressure are functions of the radial and axial coordinates with the azimuthal component of velocity being an arbitrary function of radial coordinate. The perturbations concerned are axisymmetric and infinitesimaly small in nature. Sufficient criteria of instability are derived for both vertically stable and unstable density stratifications¹ and this clearly shows that the necessary and sufficient character of Rayleigh's criterion does not hold good in the present circumstances.     

On the stability of jetlike magnetohydrodynamic flows

The linear stability problem is under study for steady axisymmetric translational flows of a density-homogeneous nonviscous incompressible ideal conducting fluid with free surface and “frozen-in” poloidal magnetic field. By the direct Lyapunov method, some sufficient conditions are obtained for the stability of these flows under small long-wave perturbations with the same symmetry. These stability conditions have partial converses; and, for unstable stationary flows, an a priori exponential lower bound is constructed on the growth of small perturbations under consideration, while the increment of the appearing exponent serves as an arbitrary positive parameter. An illustrative analytical example is given of steady flows with superimposed small long-wave axisymmetric perturbations growing in time in accordance with the estimate.     

Singularities on Free Surfaces of Fluid Flows

Isolated singularities on free surfaces of two-dimensional and axially symmetric three-dimensional steady potential flows with gravity are considered. A systematic study is presented, where known solutions are recovered and new ones found. In two dimensions, the singularities found include those described by the Stokes solution with a 120° angle, Craya's flow with a cusp on the free surface, Gurevich's flow with a free surface meeting a rigid plane at 120° angle, and Dagan and Tulin's flow with a horizontal free surface meeting a rigid wall at an angle less than 120°. In three dimensions, the singularities found include those in Garabedian's axially symmetric flow about a conical surface with an approximately 130° angle, flows with axially symmetric cusps, and flows with a horizontal free surface and conical stream surfaces. The Stokes, Gurevich, and Garabedian flows are exact solutions. These are used to generate local solutions, including perturbations of the Stokes solution by Grant and Longuet-Higgins and Fox, perturbations of Gurevich's flow by Vanden-Broeck and Tuck, asymmetric perturbations of Stokes flow and nonaxisymmetric perturbations of Garabedian's flow. A generalization of the Stokes solution to three fluids meeting at a point is also found.     

非牛顿流体非定常旋转流动计算机智能解析理论

计算机符号运算科学是人工智能的前沿方向。计算机软件Macsyma是完成符号运算的有力工具。应用德国Darmstadt大学的计算机软件Macsyma、与数学方法和流变学模型结合,研究了Oldroyd B流体由一类定常状态向另一定常状态转变的非定常流动过程。采用改进的Kantorovich方法和符号运算软件,把该问题的3阶偏微分方程的初、边值问题化为各级近似的2阶常微分方程问题。并给出了1级、2级和3级近似方程的解析形式解答。该研究表明了计算机符号处理解决应用数学和力学问题的潜力,同时指出了由一定常状态向另一定常状态转变的非牛顿流动过程,可以经历无限多途径,这一现象是由于本构方程的非线性性质引起的。     

Linear oscillatory hydromagnetic flow between two coaxial cylinders

Hydromagnetic flow between two coaxial circular cylinders is discussed when the inner cylinder oscillates axially under a radial magnetic field. Exact solution is given for the case of a perfectly conducting fluid. Expressions for velocity, induced magnetic field, current density, electric field, viscous drag and energy transfer are derived and expressed in polar forms so as to facilitate the study of magnitude and phase variations. Current sheets are found to exist on the two boundaries.     

Global subsonic and subsonic-sonic flows through infinitely long axially symmetric nozzles

In this paper, we establish existence of global subsonic and subsonic-sonic flows through infinitely long axially symmetric nozzles by combining variational method, various elliptic estimates and a compensated compactness method. More precisely, it is shown that there exist global subsonic flows in nozzles for incoming mass flux less than a critical value; moreover, uniformly subsonic flows always approach to uniform flows at far fields when nozzle boundaries tend to be flat at far fields, and flow angles for axially symmetric flows are uniformly bounded away from π/2; finally, when the incoming mass flux tends to the critical value, subsonic-sonic flows exist globally in nozzles in the weak sense by using angle estimate in conjunction with a compensated compactness framework.     

Exact analytic solutions for the unsteady flow of a non-Newtonian fluid between two cylinders with fractional derivative model

The velocity field and the associated shear stress corresponding to the torsional oscillatory flow of a generalized Maxwell fluid, between two infinite coaxial circular cylinders, are determined by means of the Laplace and Hankel transforms. Initially, the fluid and cylinders are at rest and after some time both cylinders suddenly begin to oscillate around their common axis with different angular frequencies of their velocities. The solutions that have been obtained are presented under integral and series forms in terms of generalized G and R functions. Moreover, these solutions satisfy the governing differential equation and all imposed initial and boundary conditions. The respective solutions for the motion between the cylinders, when one of them is at rest, can be obtained from our general solutions. Furthermore, the corresponding solutions for the similar flow of ordinary Maxwell fluid are also obtained as limiting cases of our general solutions. At the end, flows corresponding to the ordinary Maxwell and generalized Maxwell fluids are shown and compared graphically by plotting velocity profiles at different values of time and some important results are remarked.     

The asymptotics of neutral curve crossing in Taylor–Dean flow

The fluid flow between a pair of coaxial circular cylinders generated by the uniform rotation of the inner cylinder and an azimuthal pressure gradient is susceptible to both Taylor and Dean type instabilities. The flow can be characterised by two parameters: a measure of the relative magnitude of the rotation and pressure effects and a non-dimensional Taylor number. This work considers the small gap, large wavenumber limit for linear perturbations when the onset of the Taylor and Dean instabilities is concurrent. A consistent, matched asymptotic solution is found across the whole annular domain and identifies five regions of interest: two boundary adjustment regions and three internal critical points. Necessary conditions for the Taylor number and wavenumber are found and interpreted with reference to the suggestion of neutral curve kinks occurring at moderate wavelengths. Received: October 21, 2003; revised: November 11, 2004     

Exact periodic traveling water waves with vorticityOndes d'eau avec tourbillons

For the classical inviscid water wave problem under the influence of gravity, described by the Euler equation with a free surface over a flat bottom, we construct periodic traveling waves with vorticity. They are symmetric waves whose profiles are monotone between each crest and trough. We use global bifurcation theory to construct a connected set of such solutions. This set contains flat waves as well as waves that approach flows with stagnation points. To cite this article: A. Constantin, W. Strauss, C. R. Acad. Sci. Paris, Ser. I 335 (2002) 797–800.     

Heat transfer and Couette flow of a chemically reacting non‐linear fluid

In this paper we study the flow and heat transfer in a chemically reacting non‐linear fluid between two long horizontal parallel flat plates that are at different temperatures. The top plate is sheared, whereas the bottom plate is fixed. The fluid is modeled as a generalized power‐law fluid whose viscosity is also assumed to be a function of the concentration. The effects of radiation are neglected. The equations are made dimensionless and the boundary value problem is solved numerically; the velocity and temperature profiles are obtained for various dimensionless numbers. Published in 2009 by John Wiley & Sons, Ltd.     

Localized asymptotic solutions of the navier-stokes equations and laminar wakes in an incompressible fluid

Equations are derived to describe the far-field laminar wake behind a body in incompressible fluid flow with an arbitrary distribution of the free-stream (unperturbed flow) velocity. For certain classes of free-stream flows, analysis of these equations enables various processes in narrow wakes or jets to be described (the interaction of the longitudinal transverse velocity components in a jet, cause it to accelerate or decelerate and conservation of the energy of the wake by distortion of its trajectory regardless of viscous dissipation). In particular, conditions are obtained for the wake growth in spiral flows, analogous to the Rayleigh conditions for the instability of two-dimensionally radially symmetric flows relative to three-dimensional short-wave perturbations.     

On the structure and growth rate of unstable modes to the Rossby‐Haurwitz wave

The normal mode instability study of a steady Rossby‐Haurwitz wave is considered both theoretically and numerically. This wave is exact solution of the nonlinear barotropic vorticity equation describing the dynamics of an ideal fluid on a rotating sphere, as well as the large‐scale barotropic dynamics of the atmosphere. In this connection, the stability of the Rossby‐Haurwitz wave is of considerable mathematical and meteorological interest. The structure of the spectrum of the linearized operator in case of an ideal fluid is studied. A conservation law for perturbations to the Rossby‐Haurwitz wave is obtained and used to get a necessary condition for its exponential instability. The maximum growth rate of unstable modes is estimated. The orthogonality of the amplitude of a non‐neutral or non‐stationary mode to the Rossby‐Haurwitz wave is shown in two different inner products. The analytical results obtained are used to test and discuss the accuracy of a numerical spectral method used for the normal mode stability study of arbitrary flow on a sphere. The comparison of the numerical and theoretical results shows that the numerical instability study method works well in case of such smooth solutions as the zonal flows and Rossby‐Haurwitz waves. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Unsteady helical flows of a generalized Oldroyd-B fluid with fractional derivative

This paper deals with the unsteady helical flows of a generalized Oldroyd-B fluid between two infinite coaxial cylinders and within an infinite cylinder. The fractional calculus approach is used in the constitutive relationship of fluid model. Exact analytical solutions are obtained with the help of integral transforms (Laplace transform, Weber transform and finite Hankel transform). The corresponding solutions for generalized second grade and Maxwell fluids as well as those for the Newtonian and ordinary Oldroyd-B fluids are also given in limiting cases. Finally, the influence of model parameters on the velocity field is also analyzed by graphical illustrations.     

Reducing of birefringence of [111]-cut crystal rod using side-pumping and suitable crystal rotation

The beam quality and output power of high power solid-state lasers is influenced by birefringence. Inhomogeneous distribution of the thermal field inside the laser crystal rod occurs due to non-uniform absorption of the pump light inside the crystal and a heat sink only at boundaries. Due to the photoelastic effect, this distribution leads to inhomogeneous thermal strains and birefringence inside the rod. Plane stress and plane strain assumptions for an axially symmetric pumped crystal have been used formerly for analytical models for calculating the birefringence. This model leads in case of an [111]-cut crystal to an axially symmetric birefringence pattern. However, the shear strains in the axial-radial plane are neglected in this former models using plane stress and plane strain assumptions. This shear strains are taken into account by full 3D numerical calculations. A threefold symmetry pattern due to the anisotropic behaviour of the photoelastic tensor, which is contradictory to the ideal use of a radial or azimuthal polarized beam, is shown by results of the birefringence simulation. A laser rod pumped at three sides with threefold symmetry is analysed in order to reduce the effect of birefringence. In this case the absorption is not axially symmetric anymore. Within the crystal in regions where pumping is stronger, the pump light absorption and consequently the temperature, the strains and birefringence are higher. The degree of three-fold symmetry of birefringence will be reduced, if the region having a low birefringence due to the photoelastic effect is more strongly pumped than the rest of domain. This means the birefringence is affected by the rotation of crystal around its [111]-axis. By an optimal rotation with respect to the edges of the crystal, smallest birefringence can be obtained. For generating radial or azimuthal polarizations, the output beam of this laser device is therefore more suitable. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A large data existence result for the stationary Boltzmann equation in a cylindrical geometry

AnL ¹-existence theorem is proved for the nonlinear stationary Boltzmann equation with hard forces and no small velocity truncation—only the Grad angular cut-off-in a setting between two coaxial rotating cylinders when the indata are given on the cylinders.     

The linear stability of inviscid shear flow of a vibrationally excited diatomic gas

The problem of the linear stability of plane-parallel shear flows of a vibrationally excited compressible diatomic gas is investigated using a two-temperature gas dynamics model. The necessary and sufficient conditions for stability of the flows considered are obtained using the energy integrals of the corresponding linearized system for the perturbations. It is proved that thermal relaxation produces an additional dissipation factor, which enhances the flow stability. A region of eigenvalues of unstable perturbations is distinguished in the upper complex half-plane. Numerical calculations of the eigenvalues and eigenfunctions of the unstable inviscid modes are carried out. The dependence on the Mach number of the carrier stream, the vibrational relaxation time τ and the degree of non-equilibrium of the vibrational mode is analysed. The most unstable modes with maximum growth rate are obtained. It is shown that in the limit there is a continuous transition to well-known results for an ideal fluid as the Mach number and τ approach zero and for an ideal gas when τ → 0.     

Approximations of steady periodic water waves in flows with constant vorticity

We provide high-order approximations to periodic travelling wave profiles and to the velocity field and the pressure beneath the waves, in flows with constant vorticity over a flat bed.     

Pulsatile flow of Herschel–Bulkley fluid through catheterized arteries – A mathematical model

The pulsatile flow of blood through catheterized artery has been studied in this paper by modeling blood as Herschel–Bulkley fluid and the catheter and artery as rigid coaxial circular cylinders. The Herschel–Bulkley fluid has two parameters, the yield stress θ and the power index n. Perturbation method is used to solve the resulting quasi-steady nonlinear coupled implicit system of differential equations. The effects of catheterization and non-Newtonian nature of blood on yield plane locations, velocity, flow rate, wall shear stress and longitudinal impedance of the artery are discussed. The existence of two yield plane locations is investigated and their dependence on yield stress θ, amplitude A, and time t are analyzed. The width of the plug core region increases with increasing value of yield stress at any time. The velocity and flow rate decrease, whereas wall shear stress and longitudinal impedance increase for increasing value of yield stress with other parameters held fixed. On the other hand, the velocity, flow rate and wall shear stress decrease but resistance to flow increases as the catheter radius ratio (ratio of catheter radius to vessel radius) increases with other parameters fixed. The results for power law fluid, Newtonian fluid and Bingham fluid are obtained as special cases from this model.