Phenomenological model of interaction of a plate with a flow

In the present paper, a finite-dimensional phenomenological model of unsteady interaction of a rigid plate with a flow is proposed. It is assumed that the plate performs translational motion across the flow. The internal dynamics of the flow is modeled by the attached second order dynamical system. It is shown that the model allows satisfactory agreement with experimental data. With the developed model an inverse problem of dynamics is examined for the situation where the plate performing uniform translational motion at some moment begins uniform deceleration and finally stops. It is shown that for sufficiently large values of the plate acceleration for some time range the flow does not resist the motion of the plate but “accelerates” it. It is shown also that the equations of motion in the context of the proposed model can be reduced to the integro-differential form, and comparison with the known model of S. M. Belotserkovsky is performed. The structural resemblance of the motion equations for a body in flow in both models is noted. The domain of applicability of the quasi-stationary model is examined. __________ Translated from Fundamentalnaya i Prikladnaya Matematika, Vol. 11, No. 7, pp. 43–62, 2005.     

Assessment of turbulence modeling for gas flow in two-dimensional convergent–divergent rocket nozzle

In the present study, the turbulent gas flow dynamics in a two-dimensional convergent–divergent rocket nozzle is numerically predicted and the associated physical phenomena are investigated for various operating conditions. The nozzle is assumed to have impermeable and adiabatic walls with a flow straightener in the upstream side and is connected to a plenum surrounding the nozzle geometry and extended in the downstream direction. In this integrated component model, the inlet flow is assumed a two-dimensional, steady, compressible, turbulent and subsonic. The physics based mathematical model of the considered flow consists of conservation of mass, momentum and energy equations subject to appropriate boundary conditions as defined by the physical problem stated above. The system of the governing equations with turbulent effects is solved numerically using different turbulence models to demonstrate their numerical accuracy in predicting the characteristics of turbulent gas flow in such complex geometry. The performance of the different turbulence models adopted has been assessed by comparing the obtained results of the static wall pressure and the shock position with the available experimental and numerical data. The dimensionless shear stress at the nozzle wall and the separation point are also computed and the flow field is illustrated. The various implemented turbulence models have shown different behavior of the turbulent characteristics. However, the shear-stress transport (SST) k–ω model exhibits the best overall agreement with the experimental measurements. In general, the proposed numerical procedure applied in the present paper shows good capability in predicting the physical phenomena and the flow characteristics encountered in such kinds of complex turbulent flow.     

Simulation of suspended substance transport on the continental shelf: Horizontal dispersion

Suspended substance dispersion in a water body is simulated in the case when the spread area is considerably larger than the depth of the water body. A model of horizontal dispersion of pollutants is formulated and analyzed. Numerical approaches to the computation of suspended substance dispersion in a water body are discussed. A meshless stochastic numerical algorithm is proposed that combines the advantages of two well-known techniques, namely, the discrete cloud method and the stochastic discrete particle method. The performance the method and its features are demonstrated by comparing numerical results with the exact solution to the model problem of turbulent dispersion of a pollution plume produced by a continuous source of suspension.     

On the nature of turbulence in a problem on the motion of a fluid behind a ledge

In the present paper, we suggest a numerical method for the analysis of the motion of a viscous incompressible fluid under the transition to the turbulent mode for an example of the numerical solution of a three-dimensional space problem on the fluid flow behind a ledge for various values of the Reynolds number. We show that, at the initial stages, the turbulence in the problem is developed via successive bifurcations of generation of a stable cycle, two-dimensional tori, and then three-dimensional tori in the infinite-dimensional phase space of the system.     

The Use of Wall Functions for Simulating the Turbulent Thermal Boundary Layer

Computational Mathematics and Mathematical Physics - An important problem in the numerical simulation of turbulent heat exchange in fluids is accurate prediction of hydrodynamic characteristics of...     

CFD of turbulent transport of particles behind a backward-facing step using a new model—k-ε-Sp

The k-ε-S_p model, describing two-dimensional gas–solid two-phase turbulent flow, has been developed. In this model, the diffusion flux and slip velocity of solid particles are introduced to represent the particle motion in two-phase flow. Based on this model, the gas–solid two-phase turbulent flow behind a vertical backward-facing step is simulated numerically and the turbulent transport velocities of solid particles with high density behind the step are predicted. The numerical simulation is validated by comparing the results of the numerical calculation with two other two-phase turbulent flow models (k-ε-A_p, k-ε-k_p) by Laslandes and the experimental measurements. This model, not only has the same virtues of predicting the longitudinal transport of the solid particles as the present practical two-phase flow models, but also can predict the lateral transport of the solid particles correctly.     

Fast spectral solutions of the double-gyre problem in a turbulent flow regime

Several semi-analytical models are considered for a double-gyre problem in a turbulent flow regime for which a reference fully numerical eddy-resolving solution is obtained. The semi-analytical models correspond to solving the depth-averaged Navier–Stokes equations using the spectral Galerkin approach. The robustness of the linear and Smagorinsky eddy-viscosity models for turbulent diffusion approximation is investigated. To capture essential properties of the double-gyre configuration, such as the integral kinetic energy, the integral angular momentum, and the jet mean-flow distribution, an improved semi-analytical model is suggested that is inspired by the idea of scale decomposition between the jet and the surrounding flow.     

Approximate explicit analytical expressions of friction factor for flow of Bingham fluids in smooth pipes using Adomian decomposition method

This study concerns with the development of explicit expressions of friction factor for the flow of Bingham fluids in smooth pipes in both laminar and turbulent regimes by using a powerful analytical technique, i.e. Adomian decomposition method (ADM) and one of its variant namely restarted ADM (RADM). The results from these expressions match very well with those obtained by well-known correlations and are found to be valid for all practical ranges. It is hoped that the derived expressions will be useful for the researchers working in the area.     

Application of numerical schemes with singling out the boundary layer for the computation of turbulent flows using eddy-resolving approaches on unstructured grids

The use of eddy-resolving approaches to solving problems on arbitrary unstructured grids is investigated. The applications of such approaches requires the use of low dissipation numerical schemes, which can lead to numerical oscillations of the solution on unstructured grids. Numerical oscillations typically occur in domains with large gradients of velocities, in particular, in the near-wall layer. It is proposed to single out the boundary layer and use a numerical scheme with increased numerical dissipation in it. The algorithm for singling out the boundary layer uses a switching function to change the parameters of the numerical scheme. This algorithm is formulated based on the BCD scheme from the family NVD. Its validity and advantages are investigated using the zonal RANS–LES approach for solving some problems of turbulent flow of incompressible fluids.     

Numerical modelling of aeroelastic behaviour of an airfoil in viscous incompressible flow

In this paper, the problem of the numerical approximation of a two-dimensional incompressible viscous fluid flow interacting with a flexible structure is considered. Due to high Reynolds numbers in the range 10⁴ − 10⁶ the turbulent character of the flow is considered and modelled with the aid of Reynolds equations coupled with the k − ω turbulence model. The structure motion is described by a system of ordinary differential equations for three degrees of freedom: vertical displacement, rotation and rotation of the aileron. The problem is discretized in space by the Galerkin Least-Squares stabilized finite element method and the computational domain is treated with the aid of Arbitrary Lagrangian Eulerian method.     

Integral equation solution for the flow due to the motion of a body of arbitrary shape near a plane interface at small Reynolds number

The problem of determining the slow viscous flow due to an arbitrary motion of a particle of arbitrary shape near a plane interface is formulated exactly as a system of three linear Fredholm integral equations of the first kind, which is shown to have a unique solution. A numerical method based on these integral equations is proposed. In order to test this method valid for arbitrary particle shape, the problem of arbitrary motion of a sphere is worked out and compared with the available analytical solution. This technique can be also extended to low Reynolds number flow due to the motion of a finite number of bodies of arbitrary shape near a plane interface. As an example the case of two equal sized spheres moving parallel and perpendicular to the interface is solved in the limiting case of infinite viscosity ratio.     

Assessment of an Eulerian CFD model for prediction of dilute droplet dispersion in a turbulent jet

The turbulent dispersion of non-evaporating droplets in an axisymmetric round jet issuing from a nozzle is investigated both experimentally and theoretically. The experimental data set has a well-defined inlet boundary with low turbulence intensity at the nozzle exit, so that droplet dispersion is not affected by the transport of nozzle-generated fluctuating motion into the jet, and is influenced solely by turbulence in the gas phase produced in the shear layer of the jet. This data set is thus ideal for testing algebraic models of droplet fluctuating motion that assume local equilibrium with the turbulence in the gas phase. Moreover, the droplet flux measurements are sufficiently accurate that conservation of the total volume flow of the droplet phase has been demonstrated. A two-fluid turbulence modelling approach is adopted, which uses the k–ε turbulence model and a simple algebraic model that assumes local equilibrium to predict the fluid and droplet turbulent correlations, respectively. We have shown that the k–ε turbulence model lacks generality for predicting the spread of momentum in jets with and without a potential core. However, in general, the model predicts the radial dispersion of droplets in the considered turbulent jet with reasonable accuracy over a broad range of droplet sizes, once deficiencies in the k–ε turbulence model are taken into account.     

Analysis of some greedy algorithms for the single-sink fixed-charge transportation problem

The single-sink fixed-charge transportation problem (SSFCTP) consists of finding a minimum cost flow from a number of nodes to a single sink. Beside a cost proportional to the amount shipped, the flow cost encompass a fixed charge. The SSFCTP is an important subproblem of the well-known fixed-charge transportation problem. Nevertheless, just a few methods for solving this problem have been proposed in the literature. In this paper, some greedy heuristic solutions methods for the SSFCTP are investigated. It is shown that two greedy approaches for the SSFCTP known from the literature can be arbitrarily bad, whereas an approximation algorithm proposed in the literature for the binary min-knapsack problem has a guaranteed worst case bound if adapted accordingly to the case of the SSFCTP.     

Turbulent Poiseuille flow with near-critical wall transpiration

An incompressible, pressure-driven, fully developed turbulent flow between two parallel walls, with an extra constant transverse velocity component, is considered. A closure condition is formulated, which relates the shear stress with the first and the second derivatives of the longitudinal mean velocity. The closure condition is derived without invoking any special hypotheses on the nature of turbulent motion, only taking advantage of the fact that the flow depends on a finite number of governing parameters. By virtue of the closure condition, the momentum equation is reduced to the boundary-value problem for a second-order differential equation, which is solved by the method of matched asymptotic expansions at high values of the logarithm of the Reynolds number based on the friction velocity. The case of near-critical transpiration, when the shear stress at the injection wall vanishes, is considered. It is shown that the maximum point on the mean velocity profile lies in a thin sublayer near the suction wall in this case. A formula for the position of the maximum point as a function of the transpiration factor is obtained. The mean velocity profiles near the suction wall are calculated. A friction law for Poiseuille flow with near-critical transpiration is found, which makes it possible to describe the relation between the shear stress at the wall, the Reynolds number, and the transpiration velocity by a single function of one variable. Direct numerical simulation of the flow for some transpiration factors is performed. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

填充烧结铜球的T型管中流体混合的大涡模拟

在FLUENT软件平台上,运用大涡模拟湍流模型及Smagorinsky-Lilly亚格子尺度模型,对填充有烧结铜球多孔介质的T型管道内冷热流体混合过程的流动与传热情况进行了数值计算,与未填充多孔介质时混合区域内的平均温度和温度波动、平均速度和速度波动等数据进行了对比,并对温度波动进行了功率谱密度分析．数值结果表明,多孔介质可有效削弱T型通道流体混合区域内的温度和速度波动,有效降低1 Hz至10 Hz频域中的温度波动的功率谱密度．     

Similarity and dimensional analysis in the plane problem of the propagation of a vertical hydraulic fracture crack in an elastic medium

The problem of the growth of a vertical hydraulic fracture crack in an unbounded elastic medium under the pressure produced by a viscous incompressible fluid is studied qualitatively and by numerical methods. The fluid motion is described in the approximation of lubrication theory. Near the crack tip a fluid-free domain may exist. To find the crack length, Irwin’s fracture criterion is used. The symmetry groups of the equations describing the hydraulic fracture process are studied for all physically meaningful cases of the degeneration of the problem with respect to the control parameters. The condition of symmetry of the system of equations under the group of scaling and time-shift transformations enables the self-similar variables and the form of the time dependence of the quantities involved in the problem to be found. It is established that at non-zero rock pressure the well-known solution of Spence and Sharp is an asymptotic form of the initial-value problem, whereas the solution of Zheltov and Khristianovich is a limiting self-similar solution of the problem. The problem of the formation of a hydraulic fracture crack taking into account initial data is solved using numerical methods, and the problem of arriving at asymptotic mode is investigated. It is shown that the solution has a self-similar asymptotic form for any initial conditions, and the convergence of the exact solutions to the asymptotic forms is non-uniform in space and time.     

The dynamics of pyramidal bodies within the framework of the local interaction model

Two-dimensional inertial motion of pyramidal bodies in a medium is investigated, on the assumption that the force exerted by the medium on their surface is described by the local interaction model. Assuming unseparated flow around the bodies and small perturbations applied at the initial time to the parameters of rectilinear motion, an analytical solution is constructed of the problem of the two-dimensional motion of slender bodies with bases whose contour is a rhombus or a star consisting of four symmetrical cycles. It is shown that the solution provides the basis for a complete parameterc analysis of the dynamics of the body and for evaluating the forces and torques experienced by the body along its trajectory. A criterion for the stability of the body is found, using which, knowing the velocity, mass and position of the body's centre of gravity, one can determine the form of the perturbed motion of the pyramidal body. It is shown that the body shape is one of the most important factors affecting the stability of motion, and that, of all bodies with the same shape and position of the centre of mass, those with the least mass have the largest reserve of stability. The analytical results are confirmed by numerical solution of the Cauchy problem for the system of equations of motion obtained without the simplifying assumptions.     

On the nature of turbulence in Rayleigh-Benard convection

In the present paper, we suggest a numerical method for the analysis of the motion of a viscous incompressible fluid under the passage to the turbulent mode in Rayleigh-Benard convection. We show that one possible scenario of the onset of turbulence in the problem is given by a subharmonic cascade of bifurcations of stable two-dimensional tori.     

Some mathematical models related to the Zhukovskii problem of the flow around a sheet pile

The seepage under a Zhukovskii sheet pile through a layer of soil underlain by a highly permeable pressurized horizon is considered. The left semi-infinite part of the roof of this horizon is simulated by an impermeable foundation. The flow when the velocity on the edges of the sheet pile is equal to infinity and, on the two water permeable parts of the boundary of the domain of motion, the flow rate takes extremal values, is investigated. The limiting cases, associated with the absence of both a backwater and an impermeable inclusion, are mentioned. The problem of seepage from a foundation pit formed by two Zhukovskii sheet piles is solved within the limits of a flow with a highly permeable pressurized stratum lying below. In the case when there is no infiltration onto the free surface, a solution of the well-known Vedernikov problem is obtained. A contact scheme, arising when there are no such indicated critical points, is considered; it is described outside the scope of the constraints imposed on the unknown conforming mapping parameters ensuring the realization of the basic mathematical model. Solutions are given for two schemes of motion in a semi-inverse formulation. The classical Zhukovskii problem is the limiting case of one of them. The special features of such models are mentioned. The Polubarinova-Kochina method is used to study all the above-mentioned flows. This method enables exact analytical representations of the elements of the motion to be obtained. The results of numerical calculations and an analysis of the effect of all the physical factors on the seepage characteristics are presented.