Analyzing the influence of compressibility on the rapid pressure–strain rate correlation in turbulent shear flows

The influence of compressibility on the rapid pressure–strain rate tensor is investigated using the Green’s function for the wave equation governing pressure fluctuations in compressible homogeneous shear flow. The solution for the Green’s function is obtained as a combination of parabolic cylinder functions; it is oscillatory with monotonically increasing frequency and decreasing amplitude at large times, and anisotropic in wave-vector space. The Green’s function depends explicitly on the turbulent Mach number M _t , given by the root mean square turbulent velocity fluctuations divided by the speed of sound, and the gradient Mach number M _g , which is the mean shear rate times the transverse integral scale of the turbulence divided by the speed of sound. Assuming a form for the temporal decorrelation of velocity fluctuations brought about by the turbulence, the rapid pressure–strain rate tensor is expressed exactly in terms of the energy (or Reynolds stress) spectrum tensor and the time integral of the Green’s function times a decaying exponential. A model for the energy spectrum tensor linear in Reynolds stress anisotropies and in mean shear is assumed for closure. The expression for the rapid pressure–strain correlation is evaluated using parameters applicable to a mixing layer and a boundary layer. It is found that for the same range of M _t there is a large reduction of the pressure–strain correlation in the mixing layer but not in the boundary layer. Implications for compressible turbulence modeling are also explored.      

Laminar co-rotating Batchelor vortex merging

The dynamics of laminar co-rotating vortex pairs without axial flow have been recently thoroughly studied through theoretical, experimental and numerical studies, which revealed different instabilities contributing to the decay of the vortices. In this paper, the objective is to extend the analysis to the case of co-rotating vortices with axial flow at low Reynolds numbers. A high-order incompressible Navier–Stokes flow solver is used. The momentum equations are spatially discretized on a staggered mesh by finite differences and all derivatives are evaluated with 10th order compact finite difference schemes with RK-4 temporal discretization. The initial condition is a linear superposition of two co-rotating circular Batchelor vortices with q = 1. It is found that there is an initial evolution that resembles the evolution that single q = 1 vortices go through. Azimuthal disturbances grow and result in the appearance of large-scale helical sheets of vorticity. With the development of these instability waves, the axial velocity deficit is weakened. The redistribution of both angular and axial momentum between the core and the surroundings drives the vortex core to a more stable configuration, with a higher q value. After these processes, the evolution is somewhat similar to a pair of co-rotating Lamb–Oseen vortices. A three-dimensional instability develops, with a large band of unstable modes, with the most amplified mode corresponding scaling with the vortex initial separation distance. P. J. S. A. Ferreira de Sousa wishes to acknowledge the support of FCT—SFRH/BD/1129/2000 and SFRH/BPD/21778/2005.     

Possibilities and Limitations of Computer Simulations of Industrial Turbulent Dispersed Multiphase Flows

We give an overview on the usage of computer simulations in industrial turbulent dispersed multiphase flows. We present a few examples of industrial flows: bubble columns and bubbly pipe flows, stirred tanks, cyclones, and a fluid catalytic cracking unit. The fluid catalytic cracking unit is used to illustrate the complexity of the physical phenomena involved, and the possibilities and limitations of the different approaches used: Eulerian–Lagrangian (particle-tracking) and Eulerian–Eulerian (two-fluid). In the first approach, the continuous phase is solved using either RANS simulations (Reynolds-Averaged Navier–Stokes simulations) or DNS/LES (Direct Numerical Simulations/Large-Eddy Simulations), and the individual particles are tracked. In the second approach, the dispersed phase is averaged, leading to two sets equations, which are quite similar to the RANS equations of single-phase flows. The Eulerian–Eulerian approach is the most commonly used in industrial applications, however, it requires a significant amount of modelling. Eulerian–Lagrangian RANS can be simpler to use; in particular in situations involving complex boundary conditions, polydisperse flows and agglomeration/breakup. The key issue for the success of the simulations is to have good models for the complex physics involved. A major weakness is the lack of good models for: the turbulence modification promoted by the particles, the inter-particle interactions, and the near-wall effects. Eulerian–Lagrangian DNS/LES can play an important role as a research tool, in order to get a better physical understanding, and to improve the models used in the RANS simulations (either Eulerian–Eulerian or Eulerian–Lagrangian).     

Numerical simulation of stabilization of the boundary layer on a surface with a porous coating in a supersonic separated flow

Stability of a supersonic (M _∞ = 5.373) boundary layer with local separation in a compression corner with a passive porous coating partly absorbing flow perturbations is considered by solving two-dimensional Navier-Stokes equations numerically. The second mode of disturbances of a supersonic boundary layer is demonstrated to be the most important one behind the boundary-layer reattachment point. The possibility of effective stabilization of these disturbances behind the reattachment point with the use of porous coatings is confirmed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 2, pp. 39–47, March–April, 2007.     

Grafting of polyamide 6 on a styrene–acrylonitrile maleic anhydride terpolymer: melt rheology at the critical gel state

In this work, we studied the melt rheology of multigraft copolymers with a styrene–acrylonitrile maleic anhydride (SANMA) terpolymer backbone and randomly grafted polyamide 6 (PA 6) chains. The multi-grafted chains were formed by interfacial reactions between the maleic anhydride groups of SANMA and the amino end groups of PA 6 during melt blending. Because of the phase separation of SANMA and PA 6, the grafted SANMA backbones formed nearly circular domains which were embedded in the PA 6 melt with a diameter in the order of 20 to 40 nm. The linear viscoelastic behaviour of PA 6/SANMA blends at a sufficiently large SANMA concentration displayed the characteristics of the critical gel state, i.e. the power relations G′ ∝ G′′ ∝ ω ^0.5. In elongation, the PA 6/SANMA blend at the critical gel state showed a non-linear strain hardening behaviour already at a very small Hencky strain. In contrast to neat PA 6, the elasticity of the PA 6/SANMA blends was strongly pronounced, which was demonstrated by recovery experiments. Rheotens tests agreed with the linear viscoelastic shear oscillations and the measurements using the elongational rheometer RME. Increasing the SANMA concentration led to a larger melt strength and a reduced drawability. The occurrence of the critical gel state can be interpreted by the cooperative motion of molecules which develops between the grafted PA 6 chains of neighbouring micelle-like SANMA domains.     

Suspended load and bed-load transport of particle-laden gravity currents: the role of particle–bed interaction

The development of particle-enriched regions (bed-load) at the base of particle-laden gravity currents has been widely observed, yet the controls and relative partitioning of material into the bed-load is poorly understood. We examine particle-laden gravity currents whose initial mixture (particle and fluid) density is greater than the ambient fluid, but whose interstitial fluid density is less than the ambient fluid (such as occurs in pyroclastic flows produced during volcanic eruptions or when sediment-enriched river discharge enters the ocean, generating hyperpycnal turbidity currents). A multifluid numerical approach is employed to assess suspended load and bed-load transport in particle-laden gravity currents under varying boundary conditions. Particle-laden flows that traverse denser fluid (such as pyroclastic flows crossing water) have leaky boundaries that provide the conceptual framework to study suspended load in isolation from bed-load transport. We develop leaky and saltation boundary conditions to study the influence of flow substrate on the development of bed-load. Flows with saltating boundaries develop particle–enriched basal layers (bed-load) where momentum transfer is primarily a result of particle–particle collisions. The grain size distribution is more homogeneous in the bed-load and the saltation boundaries increase the run-out distance and residence time of particles in the flow by as much as 25% over leaky boundary conditions. Transport over a leaky substrate removes particles that reach the bottom boundary and only the suspended load remains. Particle transport to the boundary is proportional to the settling velocity of particles, and flow dilution results in shear and buoyancy instabilities at the upper interface of these flows. These instabilities entrain ambient fluid, and the continued dilution ultimately results in these currents becoming less dense than the ambient fluid. A unifying concept is energy dissipation due to particle–boundary interaction: leaky boundaries dissipate energy more efficiently at the boundary than their saltating counterparts and have smaller run-out distance.

---

     

Analysis and Modeling of the Turbulent Diffusion of Turbulent Kinetic Energy in Natural Convection

Buoyant flows often contain regions with unstable and stable thermal stratification from which counter gradient turbulent fluxes are resulting, e.g. fluxes of heat or of any turbulence quantity. Basing on investigations in meteorology an improvement in the standard gradient-diffusion model for turbulent diffusion of turbulent kinetic energy is discussed. The two closure terms of the turbulent diffusion, the velocity-fluctuation triple correlation and the velocity-pressure fluctuation correlation, are investigated based on Direct Numerical Simulation (DNS) data for an internally heated fluid layer and for Rayleigh–Bénard convection. As a result it is decided to extend the standard gradient-diffusion model for the turbulent energy diffusion by modeling its closure terms separately. Coupling of two models leads to an extended RANS model for the turbulent energy diffusion. The involved closure term, the turbulent diffusion of heat flux, is studied based on its transport equation. This results in a buoyancy-extended version of the Daly and Harlow model. The models for all closure terms and for the turbulent energy diffusion are validated with the help of DNS data for internally heated fluid layers with Prandtl number Pr = 7 and for Rayleigh–Bénard convection with Pr = 0.71. It is found that the buoyancy-extended diffusion model which involves also a transport equation for the variance of the vertical velocity fluctuation gives improved turbulent energy diffusion data for the combined case with local stable and unstable stratification and that it allows for the required counter gradient energy flux.     

Modeling jet flows caused by the incidence of a shock wave on a profiled free surface

The problem of the incidence of a shock wave with a front-pressure amplitude of about 30 GPa at the profiled free surface of an aluminum sample is studied. It is shown that in the case of large perturbations (amplitude 1 mm and wavelength 10 mm), jet flows occur on the free surface. The data obtained are described using a kinetic fracture model that takes into account the damage initiation and growth in the material due to tensile stress and shear strain. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 1, pp. 16–23, January–February, 2007.     

Effect of the Source Term on Steady Free Convection Boundary Layer Flows over an Vertical Plate in a Porous Medium. Part I

The Darcy free convection boundary layer flow over a vertical flat plate is considered in the presence of volumetric heat generation/absorption. In the present first part of the paper it is assumed that the heat generation/absorption takes place in a self-consistent way, the source term q ^′′′≡ S of the energy equation being an analytical function of the local temperature difference T − T _∞. In a forthcoming second part, the case of the externally controlled source terms S = S(x,y ) will be considered. It is shown that due to the presence of S, the physical equivalence of the up- and downflows gets in general broken, in the sense that the free convection flow over the upward projecting hot plate (“upflow”) and over its downward projecting cold counterpart (“downflow”) in general become physically distinct. The consequences of this circumstance are examined for different forms of S. Several analytical solutions are given. Some of them describe algebraically decaying boundary layers which can also be recovered as limiting cases of exponentially decayingones. This asymptotic phenomenon is discussed in some detail.     

Effect of Schmidt number on the structure and propagation of density currents

The results of a numerical study of two- and three-dimensional Boussinesq density currents are described. They are aimed at exploring the role of the Schmidt number on the structure and dynamics of density driven currents. Two complementary approaches are used, namely a spectral method and a finite-volume interface capturing method. They allow for the first time to describe density currents in the whole range of Schmidt number 1 ≤ Sc ≤ ∞ and Reynolds number 10² ≤ Re ≤ 10⁴. The present results confirm that the Schmidt number only weakly influences the structure and dynamics of density currents provided the Reynolds number of the flow is large, say of O(10⁴) or more. On the contrary low- to moderate-Re density currents are dependant on Sc as the structure of the mixing region and the front velocities are modified by diffusion effects. The scaling of the characteristic density thickness of the interface has been confirmed to behave as (ScRe)^−1/2. Three-dimensional simulations suggest that the patterns of lobes and clefts are independent of Sc. In contrast the Schmidt number is found to affect dramatically (1) the shape of the current head as a depression is observed at high-Sc, (2) the formation of vortex structures generated by Kelvin–Helmholtz instabilities. A criterion is proposed for the stability of the interface along the body of the current based on the estimate of a bulk Richardson number. This criterion, derived for currents of arbitrary density ratio, is in agreement with present computed results as well as available experimental and numerical data.      

Objective flow classification parameters and their use in general steady flows

Three quantitative flow classification parameters have been studied in the context of Tanner and Huilgol’s suggestion of strong and weak flows. Seen in this context, the different types of streamlines possible for general 3-D flows furnish no indication with respect to the flow strength. This is in total contrast to 2-D flows, where the type of the streamline and the strength of the flow go hand in hand. Astarita’s [J Non-Newton Fluid Mech, 6:69–76, 1979] flow classification parameter takes care of this fact and, if properly generalized, can be applied to more general flows: Two other flow classification parameters also have their basis in homogeneous 2-D flows, but their generalization leads, for general flows, to nonuniqueness and other unacceptable results. For 3-D flows, none of the parameters can quantitatively be used in general, and additional parameters, with their basis outside the 2-D flow regime, seem to be called for.

Prediction of separation flows around a 6：1 prolate spheroid using RANS/LES hybrid approaches

This paper presents hybrid Reynolds-averaged Navier–Stokes (RANS) and large-eddy-simulation (LES) methods for the separated flows at high angles of attack around a 6:1 prolate spheroid. The RANS/LES hybrid methods studied in this work include the detached eddy simulation (DES) based on Spalart–Allmaras (S–A), Menter’s k–ω shear-stress-transport (SST) and k–ω with weakly nonlinear eddy viscosity formulation (Wilcox–Durbin+, WD+) models and the zonal-RANS/LES methods based on the SST and WD+ models. The switch from RANS near the wall to LES in the core flow region is smooth through the implementation of a flow-dependent blending function for the zonal hybrid method. All the hybrid methods are designed to have a RANS mode for the attached flows and have a LES behavior for the separated flows. The main objective of this paper is to apply the hybrid methods for the high Reynolds number separated flows around prolate spheroid at high-incidences. A fourth-order central scheme with fourth-order artificial viscosity is applied for spatial differencing. The fully implicit lower–upper symmetric-Gauss–Seidel with pseudo time sub-iteration is taken as the temporal differentiation. Comparisons with available measurements are carried out for pressure distribution, skin friction, and profiles of velocity, etc. Reasonable agreement with the experiments, accounting for the effect on grids and fundamental turbulence models, is obtained for the separation flows. The project supported by the National Natural Science Foundation of China (10502030 and 90505005).     

Symmetries and exact solutions of the shallow water equations for a two-dimensional shear flow

This paper considers nonlinear equations describing the propagation of long waves in two-dimensional shear flow of a heavy ideal incompressible fluid with a free boundary. A nine-dimensional group of transformations admitted by the equations of motion is found by symmetry methods. Two-dimensional subgroups are used to find simpler integrodifferential submodels which define classes of exact solutions, some of which are integrated. New steady-state and unsteady rotationally symmetric solutions with a nontrivial velocity distribution along the depth are obtained. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 5, pp. 41–54, September–October, 2008.     

Numerical simulation of the flow in a variable-section channel with pulsed-periodic energy supply

Results of numerical simulations of a quasi-one-dimensional unsteady flow in a channel considered as an element of an air-breathing engine are presented. The influence of parameters of energy supplied in the pulsed-periodic mode (power, pulse frequency, and distribution of energy sources along the channel) on the characteristics of the flow with Mach numbers M ₀ = 2.4–4.0 at the channel entrance is determined. A channel configuration that allows the energy supply distribution to be found from the condition of restriction of the maximum value of the gas temperature is proposed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 1, pp. 3–11, January–February, 2009.     

Analysis of fragmentation of deformation in Ni3Ge single crystals

The strain distribution on two faces ofNi ₃Ge single crystals compressed to a strain ε-14.16% at T=77.673K was studied by the grid method. It is shown that temperature has a significant effect on the strain distribution. Fragmentation of local strain due to shape change in specimens during active loading was established by the method of main components. Tomsk State Architecture-Construction Academy, Tomsk 634003. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 1, pp. 154–159, January–February, 1998.     

Structure of supersonic turbulent flows in the vicinity of inclined backward-facing steps

Supersonic (M_∞ = 2−5) turbulent flows in the vicinity of a two-dimensional backward-facing step with an inclined leeward side are considered by methods of mathematical modeling. The wave structure of the flow with a varied angle of inclination of the leeward side of the step and a varied free-stream Mach number is considered. __________ Translated From Prikladnaya Mekhanika I Tekhnicheskaya Fizika, Vol. 47, No. 6, pp. 48–58, November–December, 2006.     

Numerical analysis of pressure oscillations over axisymmetric spiked blunt bodies at Mach 6.80

The pressure oscillations over a forward facing spike attached to an axisymmetric blunt body are simulated by solving time-dependent compressible Navier–Stokes equations. The governing fluid flow equations are discretized in spatial coordinates employing a finite volume approach which reduces the equations to semidiscretized ordinary differential equations. Temporal integration is performed using the two-stage Runge–Kutta time stepping scheme. A global time step is used to obtain a time-accurate numerical solution. The numerical computation is carried out for a freestream Mach number of 6.80 and for spike length to hemispherical diameter ratios of 0.5, 1.0 and 2.0. The flow features around the spiked blunt body are characterized by a conical shock wave emanating from the spike tip, a region of separated flow in front of the hemispherical cap, and the resulting reattachment shock wave. Comparisons of the numerical results are made with the available experimental results, such as schlieren pictures and the surface pressure distribution along the spiked blunt body. They are found to be in good agreement. Spectral analysis of the computed pressure oscillations are performed employing fast Fourier transforms. The surface pressure oscillations over the spike and phase plots exhibit a behaviour analogous to that of the Van der Pol equation for a self-sustained oscillatory flow. Received 28 February 2001 / Accepted 17 January 2002     

Self-consistent modeling of entangled network strands and linear dangling structures in a single-strand mean-field slip-link model

Linear viscoelastic (LVE) measurements as well as non-linear elongation measurements have been performed on stoichiometrically imbalanced polymeric networks to gain insight into the structural influence on the rheological response (Jensen et al., Rheol Acta 49(1):1–13, 2010). In particular, we seek knowledge about the effect of dangling ends and soluble structures. To interpret our recent experimental results, we exploit a molecular model that can predict LVE data and non-linear stress–strain data. The slip-link model has proven to be a robust tool for both LVE and non-linear stress–strain predictions for linear chains (Khaliullin and Schieber, Phys Rev Lett 100(18):188302–188304, 2008, Macromolecules 42(19):7504–7517, 2009; Schieber, J Chem Phys 118(11):5162–5166, 2003), and it is thus used to analyze the experimental results. Initially, we consider a stoichiometrically balanced network, i.e., all strands in the ensemble are attached to the network in both ends. Next we add dangling strands to the network representing the stoichiometric imbalance, or imperfections during curing. By considering monodisperse network strands without dangling ends, we find that the relative low-frequency plateau, G₀/G_N⁰G_0/G_N^0, decreases linearly with the average number of entanglements. The decrease from G_N⁰G_N^0 to G ₀ is a result of monomer fluctuations between entanglements, which is similar to “longitudinal modes” in tube theory. It is found that the slope of G′ is dependent on the fraction of network strands and the structural distribution of the network. The power-law behavior of G ^″ is not yet captured quantitatively by the model, but our results suggest that it is a result of polydisperse dangling and soluble structures.     

Rheological properties of binary associating polymers

Dynamics of associating polymer solutions above the reversible gelation point are studied. Each macromolecule consists of a soluble backbone (B) and a small fraction of specific strongly interacting groups (A or C stickers) attached to B. A mixture of B–A and B–C associating polymers with 1:1 stoichiometric ratio is considered. As a result of AC association, the polymers reversibly gelate above the overlap concentration. It is shown that (1) the network strands are linear complexes (double chains) of B–A and B–C; (2) “diffusion” of the network junction points is characterized by an apparent activation energy, which can be significantly higher than the energy of one AC bond; (3) most importantly, the randomness of sticker distribution along the chain can significantly slow down the network relaxation leading to a markedly non-Maxwellian viscoelastic behavior. The theory elucidates the most essential features of rheological behavior of polysaccharide associating systems (with A = adamantyl moiety, C = β-cyclodextrin, B = either chitosan or hyaluronan) including similar behavior of G ^′ and G ^″ in a wide frequency range, strong temperature dependence of the characteristic frequency ω _x , and an extremely strong effect of added free stickers (fC) on the dynamics. This paper was presented at Annual European Rheology Conference (AERC) held in Hersonisos, Crete, Greece, April 27–29, 2006.     

Numerical simulation of supersonic flow past reentry capsules

The flow fields over ARD (ESA's atmospheric reentry demonstrator), OREX (orbital reentry experiments) and spherically blunted cone-flare reentry configurations are numerically obtained by solving time-dependent, axisymmetric, compressible Navier–Stokes equations for freestream Mach numbers range of 1.2–6.0. The fluid dynamics are discretized in spatial coordinates employing a finite volume approach which reduces the governing equations to semi discretized ordinary differential equations. Temporal integration is performed using the multistage Runge–Kutta time-stepping scheme. A local time step is used to achieve steady-state solution. The numerical simulation is carried out on a structured grid. The flow-field features around the reentry capsule, such as bow shock wave, sonic line, expansion fan and recirculating flow in the base region are obtained. A good agreement is found between the calculated value of aerodynamic drag coefficient of the spherically blunted cone/fare reentry configuration with the experimental data. The effects of geometrical parameters, such as radius of the spherical cap, half cone angle, with sharp shoulder edge and with smooth shoulder edge on the flow-field have been numerically investigated for various reentry configuration which will be useful for optimization of the reentry capsule. PACS 47.11.Df, 47.40.Ki