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.     

Initial Algebraic Growth of Small Angular Dependent Disturbances in Pipe Poiseuille Flow

The evolution of small, angular dependent velocity disturbances in laminar pipe flow is studied. In particular, streamwise independent perturbations are considered. To fully describe the flow field, two equations are required, one for the radial and the other for the streamwise velocity perturbation. Whereas the former is homogeneous, the latter has the radial velocity component as a forcing term. First, the normal modes of the system are determined and analytical solutions for eigenfunctions, damping rates, and phase velocities are calculated. As the azimuthal wave number (n) increases, the damping rate increases and the phase velocities decrease. Particularly interesting are results showing that the phase velocities associated with the streamwise eigenfunctions are independent of the radial mode index when n = 1, and when n = 5 the same is obtained for phase velocities associated with the eigenfunctions of the radial component. Then, the initial value problem is treated and the time development of the disturbances is determined. The radial and the azimuthal velocity components always decay but, owing to the forcing, the streamwise component shows an initial algebraic growth, followed by a decay. The kinetic energy density is used to characterize the induced streamwise disturbance. Its dependence on the Reynolds number, the radial mode, and the azimuthal wave number is investigated. With a normalized initial disturbance, n = 1 gives the largest amplification, followed by n = 2 etc. However, for small times, higher values of n are associated with the largest energy density. As n increases, the distribution of the streamwise velocity perturbation becomes more concentrated to the region near the pipe wall.     

Rohrströmung mit schwachem Drall

Summary The steady laminar flow at high Reynolds' numbers in a straight tube of circular cross-section is investigated under the condition that swirl is present. Given any axially symmetrical distribution of the rotational velocities, the flow downstream is calculated in the special case that the rotational velocities are small compared to the axial velocity. The problem is reduced to an eigen-value problem of an ordinary linear differential equation of second order, the theory of which is developed in detail. Numerical values for the five lowest eigen-values and the corresponding eigen-functions are given. As an example, the decay of an initially rigid rotation is considered.

---

---

HerrnFr. A. Willers zum siebzigsten Geburtstag     

Computational study of decaying annular vortex flow using the Rε/k − ε turbulence model

In this article, we present three dimensional CFD study of turbulent vortex flow in an annular passage using OpenFOAM 1.6. The vortex flow is generated by introducing the flow through a tangential entry to the passage. For the analysis presented in this article, turbulence was modeled using the R_ε/k − ε model, in addition, a comparison between such model with the standard k − ε model was conducted and discussed. The main characteristics of the flow such as vortex structure and recirculation zone were investigated. It was found that flow is subjected to Rankine vortex structure with three forced vortex regimes and a free vortex region near to the outer wall. The phenomenon of vortex decay was investigated by depicting the swirl number trend along the axial direction of the flow domain. It was found that the vortex decay is subjected to an exponential decay behavior. New coefficients for the exponential decay correlation were derived based on local values of velocity components in different radial planes.     

Comparison of turbulence models in simulating swirling pipe flows

Swirling flow is a common phenomenon in engineering applications. A numerical study of the swirling flow inside a straight pipe was carried out in the present work with the aid of the commercial CFD code fluent. Two-dimensional simulations were performed, and two turbulence models were used, namely, the RNG k–ε model and the Reynolds stress model. Results at various swirl numbers were obtained and compared with available experimental data to determine if the numerical method is valid when modeling swirling flows. It has been shown that the RNG k–ε model is in better agreement with experimental velocity profiles for low swirl, while the Reynolds stress model becomes more appropriate as the swirl is increased. However, both turbulence models predict an unrealistic decay of the turbulence quantities for the flows considered here, indicating the inadequacy of such models in simulating developing pipe flows with swirl.     

New trends towards lower energy estimates and optimality for nonlinearly damped vibrating systems

We present an approach based on comparison principles for energy and interpolation properties to derive lower energy estimates for nonlinearly either locally damped or boundary damped vibrating systems. We show how the dissipation relation provides strong information on the asymptotic behavior of the energy of solutions. The geometrical situations are either one-dimensional, or radial two-dimensional or three-dimensional for annulus domains. We also consider the case of general domains, but in this case, for solutions with bounded velocities in time and space. In all these cases, the nonlinear damping function is assumed to have arbitrary (strictly sublinear) growth at the origin. We give results for strong solutions and stronger lower estimates for smoother solutions. The results are presented in two forms, either on the side of energy comparison principles, or through time-pointwise lower estimates. Under additional geometric assumptions, we give the resulting lower and upper estimates for four representative examples of damping functions. We further give a “weak” lower estimate (in the sense of a certain lim supt→∞) and an upper estimate of the velocity for smoother solutions in case of general damping functions and for radial, as well as multi-dimensional domains. We also discuss these estimates in the framework of optimality, which is not proved here, and indicate open problems raised by these results.     

Constrained radial symmetry for monotone elliptic quasilinear operators

Some overdetermined problems associated to monotone elliptic quasilinear operators are investigated. A model operator is the p-Laplacian. Assuming that a solution exists, the domain of our problem is shown to be either a ball centered at the origin or an annulus centered at the origin. In the special case of the Laplace equation, a result of approximate radial symmetry is also obtained. Proofs are based on comparisons with radial solutions.     

Numerical simulation of internal flow fields of swirl coaxial injector in a hot environment

Based on the fractional volume of fluid (VOF), a pure Eulerian model for defining and capturing the gas/liquid interface is developed in this paper. This model can describe gas/liquid interface in high refinement, which is better than the original VOF methodology. To validate the proposed model and the algorithm, the computational code is employed to predict the flow performance in a cylindrical swirl injector under cold-flow condition, and the predicted results agree well with experimental measurements. Furthermore, the proposed model is used to simulate gas-liquid reacting flows inside a gas/liquid coaxial swirl injector operating in a hot environment. The turbulent combustion process is simulated with the k−ε−f−g model. The numerical simulation is carried out under actual operating condition of the coaxial injector. The injector performances, such as liquid film thickness, liquid film injection velocity, spray angle, pressure drop, are obtained based on the detailed information of the internal flow field. The predicted results also show that droplets are shed from the liquid film in the recess cup of the coaxial injector because of the large velocity gradient between the gas and liquid streams, and a burning area, which is characterized by high temperature, is present inside the injector.     

Statistical characteristics and time autocorrelation coefficients of the turbulent swirl flow in pipe

In this paper is presented experimental investigation of the turbulent swirl flow in pipe generated by three axial fans of various geometries. One of these fans (model ZP) generates Rankine swirl, while other two, industrial fans, produce mainly solid body circumferential velocity profile. One-component laser Doppler anemometry (LDA) is employed in this research. Downstream transformations of the non-dimensional time-averaged velocity profiles are analyzed. Distributions of turbulence levels are discussed for all three fans in both measuring sections, as well as statistical moments of higher orders for fluctuations in the axial, radial and circumferential directions. Applied correlation theory revealed turbulence structure and its statistical nature. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the development of a wake vortex in inviscid flow

The evolution of an initial perturbation in an axisymmetric subsonic normal inviscid gas flow through a pipe is directly simulated. The basic (unperturbed) flow has a zero radial velocity component, while its axial velocity component (along the axis of symmetry) increases or decreases linearly with the radius. The perturbation is specified as a swirl (rotation about the axis) with a positive or negative velocity vanishing on the central axis and the lateral surface. Irrespective of its direction, the swirl gives rise to a steady-state vortex carried by the flow. It shape is spherical (contiguous to the rotation axis) or circular (sliding along the impermeable lateral surface).     

Spatial decay bounds for time dependent magnetohydrodynamic geophysical flow

This article is concerned with spatial decay bounds for the time dependent magnetohydrodynamic geophysical flow in an infinite pipe when homogeneous lateral surface boundary conditions are applied. Assuming that the entrance velocity and magnetic field data are small enough and the fluid flow converges to laminar flow as the distance down the pipe tends to infinity, we derive a second order differential inequality that leads to an exponential decay estimate for the “energy” associated with the velocity and magnetic field represented by the difference between the entrance flow and fully developed laminar flow. We also show how to establish the explicit decay bounds for the total energy.     

Burgers turbulence initialized by a regenerative impulse

In this article we look at a one-dimensional infinitesimal particle system governed by the completely inelastic collision rule. Considering uniformly spread mass, we feed the system with initial velocities, so that when time evolves the corresponding velocity field fulfils the inviscid Burgers equation. More precisely, we suppose here that the initial velocities are zero, except for particles located on a stationary regenerative set for which the velocity is some given constant number. We give results of a large deviation type. First, we estimate the probability that a typical particle is located at time 1 at distance at least D from its initial position, when D tends to infinity. Its behaviour is related to the left tail of the gap measure of the regenerative set. We also show the same asymptotics for the tail of the shock interval length distribution. Second, we analyse the event that a given particle stands still at time T as T tends to infinity. The data to which we relate its behaviour are the right tail of the gap measure of the regenerative set. We conclude with some results on the shock structure.     

A model for the corona of a spherical star cluster

A general model for the outer layers of steady spherical star systems is considered. The corona is assumed to be formed by stars flying out of the system along ballistic orbits. Expressions for density runs and velocity dispersions are found. Models of three kinds are studied, namely, systems with isotropic velocity distribution, spheres with purely radial orbits, and models with ellipsoidal velocity distribution. The following isotropic distribution functions were considered as examples: the truncated Maxwellian distribution, King’s model, and the spherical analogue of Perek’s model. The truncated Maxwellian distribution of radial velocities is an example of a radial orbit system. The truncated Schwarzschildian velocity distribution and generalized polytropes introduced by Kuzmin and Veltmann are examples of models with ellipsoidal velocity distribution. Density asymptotics are found for all models provided that the self-gravitation of coronas is neglected.     

Propagation of weak perturbations in cracked porous media

A three-velocity, three-pressure mathematical model is proposed which enables one to study wave processes in the case of a double porosity, deformable, fluid-saturated medium. This model takes account of the differences in the velocities and pressures in pore systems of different characteristic scales of the pores, fluid exchange between these pore systems and the unsteady forces due to interphase interactions. It is established that a single transverse and three longitudinal waves: one deformation wave and two filtration waves, propagate in such a medium. The existence of two filtration waves is associated with the two different characteristic scales of the pores and the difference in the velocities and pressures of the fluid in these pore systems. The filtration waves decay considerably more rapidly than the deformation and transverse waves. The velocities of the deformation and transverse waves are mainly determined by the elastic moduli of the skeleton. The velocity and decay of the first filtration wave depend strongly on the intensity of the interphase interaction force while the velocity of the second filtration wave depends strongly on the rate of mass exchange between the pores and the cracks. The rate of decay of the second filtration wave is significantly higher than that of the first filtration wave.     

Injectivity of rotation invariant windowed Radon transforms

We consider rotation invariant windowed Radon transforms that integrate a function over hyperplanes by using a radial weight (called window). T. Quinto proved their injectivity for square integrable functions of compact support. This cannot be extended in general. Actually, when the Laplace transform of the window has a zero with positive real part δ, the windowed Radon transform is not injective on functions with a Gaussian decay at infinity, depending on δ. Nevertheless, we give conditions on the window that imply injectivity of the windowed Radon transform on functions with a more rapid decay than any Gaussian function.     

Local radial point interpolation method for the fully developed magnetohydrodynamic flow

In this paper, a local radial point interpolation method (LRPIM) is presented to obtain the numerical solutions of the coupled equations in velocity and magnetic field for the fully developed magnetohydrodynamic (MHD) flow through a straight duct of rectangular section with arbitrary wall conductivity and orientation of applied magnetic field. Local weak forms are developed using weighted residual method locally for the governing equations of fully developed MHD flow. The shape functions from LRPIM possess the delta function property. Therefore, essential boundary conditions can be applied as easily as that in the finite-element method. The implementation procedure of LRPIM method is node based, and it doesn’t need any “mesh” or “element”. Computations have been carried out for different Hartmann numbers, wall conductivities and orientations of applied magnetic field.     

川渝裂缝性地层自动压井环空多相压力波速特性研究北大核心CSCD

考虑虚拟质量力、环空沿程压力、气液相间阻力、气体滑脱、环空空隙率等因素,基于小扰动理论,提出了裂缝性地层自动压井环空多相压力波速数学模型,结合半显式差分方法,以彭州PZ-5-3D井(垂深5827 m)为实例,对模型编程求解.结果表明:裂缝性地层出气具有段塞流特点,随空隙率增大,压力波速呈现先减小后增大趋势;空隙率在0%至16%区间,压力波速以液弹为主,压力波速呈急剧下降趋势;空隙率在16%至40%区间,压力波速趋于平缓恒定值;空隙率在42%至100%区间,压力波速呈现增大趋势,压力波速以气弹为主;随环空井深减小,环空空隙率减小,压力波速整体呈现减小趋势;随压井循环排气井口回压增大,压力波速整体呈现增大趋势;环空空隙率在0%至13%区间内,气体滑脱速度对压力波速影响不大;环空空隙率在13%至85%区间内,随气体滑脱速度增大,压力波速呈现减小趋势;节流阀调阀时间间隔与井底压力响应时间具有跟随性,随井底压力响应时间增大,调阀时间间隔增大.     

Uniqueness of smooth radial balanced metrics on the disc

We show that the usual Poincaré metric is the only radial balanced metric on the disc with not too wild boundary behaviour. Additionally, we identify explicitly all radial metrics with such boundary behaviour which satisfy the balanced condition as far as germs at the boundary are concerned. Related results for the annulus and the punctured disc are also established.     

Recent Developments in Extended Thermodynamics of Dense and Rarefied Polyatomic Gases

We summarize the recent results and current open problems in extended thermodynamics (ET) of both dense and rarefied polyatomic gases. (i) We review, in particular, extended thermodynamics with 14 independent fields (ET14), that is, the mass density, the velocity, the temperature, the shear stress, the dynamic pressure, and the heat flux. (ii) We explain that, in the case of rarefied polyatomic gases, molecular extended thermodynamics with 14 independent fields (MET14) basing on the kinetic moment theory with the maximum entropy principle can be developed. ET14 and MET14 are fully consistent with each other. (iii) We show that the ET13 theory of rarefied monatomic gases is derived from the ET14 theory as a singular limit. (iv) We discuss briefly some typical applications of the ET14 theory. (v) We study the simple case of ET theory with 6 independent fields (ET6). (vi) The METn theories (n>14) are presented briefly. We analyze, in particular, the dependence of the characteristic velocities for increasing number of moments.