Guidelines for Modeling a 2D Rough Wall Channel Flow

The effects of 2-D roughness elements on the Reynolds stress anisotropy tensor, and of the energy dissipation rate anisotropy tensor of a turbulent channel flow are investigated using data obtained from direct numerical simulations (DNSs). The roughness elements consist of transverse square rods of size , placed on one wall of the channel only. While is kept constant (, is the half-width of the channel), the spacing between the rods is varied from to . The results show that the variation in can dramatically change the structure of the wall region flow. The modeling of the near-wall region needs to reflect the structural changes caused by the variation in . On the basis of the Reynolds stress budgets, attention needs to be given to the turbulent energy and pressure diffusion terms while local isotropy may be a reasonable approximation for the energy dissipation rate, especially over a range of for which the drag is near its maximum.     

Flow Visualization of Superbursts and of the Log-Layer in a DNS at $$ \operatorname{Re} _{\tau } = 950 $$

This paper uses direct numerical simulations (DNS) of turbulent flow in a channel at (Del álamo, Jiménez, Zandonade, Moser J Fluid Mech 500:135–144, 2004) to provide a picture of the turbulent structures making large contributions to the Reynolds shear stress. Considerable work of this type has been done for the viscous wall region at smaller , for which a log-layer does not exist. Recent PIV measurements of turbulent velocity fluctuations in a plane parallel to the direction of flow have emphasized the dominant contribution of large scale structures in the outer flow. This prompted Hanratty and Papavassiliou (The role of wall vortices in producing turbulence. In: Panton, R.L. (ed) Self-sustaining Mechanism of Wall Turbulence. Computational Mechanics Publications, Southampton, pp. 83–108, 1997) to use DNS at to examine these structures in a plane perpendicular to the direction of flow. They identified plumes which extend from the wall to the center of a channel. The data at are used to explore these results further, to examine the structure of the log-layer, and to test present notions about the viscous wall layer.     

A nonlinear kp-εp particle two-scale turbulence model and its application

A particle nonlinear two-scale turbulence model is proposed for simulating the anisotropic turbulent two-phase flow. The particle kinetic energy equation for two-scale fluctuation, particle energy transfer rate equation for large-scale fluctuation, and particle turbulent kinetic energy dissipation rate equation for small-scale fluctuation are derived and closed. This model is used to simulate gas–particle flows in a sudden-expansion chamber. The simulation is compared with the experiment and with those obtained by using another two kinds of tow-phase turbulence model, such as the single-scale two-phase turbulence model and the particle two-scale second-order moment (USM) two-phase turbulence model. It is shown that the present model gives simulation in much better agreement with the experiment than the single-scale two-phase turbulence model does and is almost as good as the particle two-scale USM turbulence model. The project supported by China Postdoctoral Science Foundation (2004036239).     

Simultaneous measurement of turbulent velocity field and surface wave amplitude in the initial stage of an open-channel flow by PIV

Particle image velocimetry (PIV) was employed to measure the two components of the turbulent velocity field in the initial stage of an open-channel flow in a streamwise-wall-normal plane, and the free-surface level was discriminated from the PIV image. The details of this technique was described and demonstrated by showing the instantaneous velocity field together with the free-surface shape, statistics of velocity field, and the wave-turbulence interaction terms. Preliminary experimental results showed that the turbulence intensity of the streamwise velocity fluctuations (u′) decreased, whereas that of the wall-normal velocity (v′) increased near the bottom wall with downstream distance in the initial stage of an open-channel flow; (g is the fluctuation of free-surface level) had a negative value, had a positive value near the free surface, and the surface-wave-affected depth deepened with downstream distance.     

Vortex formation process in gravity-driven starting jets

Experiments were conducted to investigate the vortex formation process in gravity-driven starting jets at three Reynolds numbers ( 3,528 and 4,716, where D is the nozzle diameter, the average discharging velocity and ν the kinematic viscosity). Planar laser induced fluorescence (PLIF) was used to visualize the flow while particle image velocimetry (PIV) was used to quantify the 2D velocity and vorticity fields. Vortex leapfrogging was found to appear at Re = 2,358, while pinch-off for the leading vortex was observed in the other two cases. Conditions for the vortex pinch-off appeared to be different from those found in the piston-driven starting jets. Although the leading vortex attained the maximum circulation and energy level very quickly after the jet was initiated, its detachment from the main stem may not be achieved necessarily. Attempt had been made to elaborate the different flow characteristics for respective cases from the related distribution of circulation and energy level via the measured velocity field by PIV.     

Resolvent Estimates and Maximal Regularity in Weighted L q -spaces of the Stokes Operator in an Infinite Cylinder

Let be an infinite cylinder of , n ≥ 3, with a bounded cross-section of C ^1,1-class. We study resolvent estimates and maximal regularity of the Stokes operator in for 1 < q, r < ∞ and for arbitrary Muckenhoupt weights ω ∈ A _r with respect to x′ ∈ Σ. The proofs use an operator-valued Fourier multiplier theorem and techniques of unconditional Schauder decompositions based on the -boundedness of the family of solution operators for a system in Σ parametrized by the phase variable of the one-dimensional partial Fourier transform. Supported by the Gottlieb Daimler- und Karl Benz-Stiftung, grant no. S025/02-10/03.     

Vector potential–vorticity relationship for the Stokes flows: application to the Stokes eigenmodes in 2D/3D closed domain

The unsteady dynamics of the Stokes flows, where , is shown to verify the vector potential–vorticity ( ) correlation , where the field is the pressure-gradient vector potential defined by . This correlation is analyzed for the Stokes eigenmodes, , subjected to no-slip boundary conditions on any two-dimensional (2D) closed contour or three-dimensional (3D) surface. It is established that an asymptotic linear relationship appears, verified in the core part of the domain, between the vector potential and vorticity, , where is a constant offset field, possibly zero.     

Forced convection gaseous slip flow in circular porous micro-channels

Laminar forced convection of gaseous slip flow in a circular micro-channel filled with porous media under local thermal equilibrium condition is studied numerically using the finite difference technique. Hydrodynamically fully developed flow is considered and the Darcy–Brinkman–Forchheimer model is used to model the flow inside the porous domain. The present study reports the effect of several operating parameters (Knudsen number (Kn), Darcy number (Da), Forchhiemer number (Γ), and modified Reynolds number ) on the velocity slip and temperature jump at the wall. Results are given in terms of the velocity distribution, temperature distribution, skin friction , and the Nusselt number (Nu). It is found that the skin friction is increased by (1) decreasing Knudsen number, (2) increasing Darcy number, and (3) decreasing Forchheimer number. Heat transfer is found to (1) decrease as the Knudsen number, or Forchheimer number increase, (2) increase as the Peclet number or Darcy number increase.     

Local Bifurcation Analysis of a Second Gradient Model for Deformations of a Rectangular Slab

In this paper we carry out a derivation of the equilibrium equations of nonlinear elasticity with an added second-gradient term proportional to a small parameter . These equations are given by a fourth order semilinear system of pdes. We discuss different types of possible boundary conditions for these equations. We then specialize the equations to a rectangular slab and study the linearized problem about a homogenous deformation. We show that these equations admit solutions representable as Fourier series in one of the independent variables. Furthermore, we obtain the characteristic equation for the eigenvalues (possible bifurcation points) for the linear problem and derive asymptotic representations for this equation for small . We used these expressions to show that in the limit as the characteristic equation for converges uniformly (in certain regions of the parameter space) to the corresponding characteristic equation for . When the base material () is that of a Blatz–Ko type, we get conditions for the existence of eigenvalues of the linear problem with and small. Our numerical results in this case indicate that the number of bifurcation points is finite when and that this number monotonically increases as . For the problem with we get conditions for the existence of local branches of non-trivial solutions.      

The Steady Motion of a Navier–Stokes Liquid Around a Rigid Body

Let be a body moving by prescribed rigid motion in a Navier–Stokes liquid that fills the whole space and is subject to given boundary conditions and body force. Under the assumptions that, with respect to a frame , attached to , these data are time independent, and that their magnitude is not “too large”, we show the existence of one and only one corresponding steady motion of , with respect to , such that the velocity field, at the generic point x in space, decays like |x|⁻¹. These solutions are “physically reasonable” in the sense of FINN [10]. In particular, they are unique and satisfy the energy equation. Among other things, this result is relevant in engineering applications involving orientation of particles in viscous liquid [14].     

Momentum and heat transfer from an asymmetrically confined circular cylinder in a plane channel

Unsteady momentum and heat transfer from an asymmetrically confined circular cylinder in a plane channel is numerically investigated using FLUENT for the ranges of Reynolds numbers as 10≤Re≤500, of the blockage ratio as 0.1≤β≤0.4, and of the gap ratio as 0.125≤γ≤1 for a constant value of the Prandtl number of 0.744. The transition of the flow from steady to unsteady (characterized by critical Re) is determined as a function of γ and β. The effect of γ on the mean drag and lift coefficients, Strouhal number (St), and Nusselt number (Nu _w ) is studied. Critical Re was found to increase with decreasing γ for all values of β. and St were found to increase with decreasing values of γ for fixed β and Re. The effect of decrease in γ on was found to be negligible for all blockage ratios investigated.     

Explicit Effective Constants for an Inhomogeneous Porothermoelastic Medium

An arbitrary anisotropic micro-inhomogeneous (composite) poroelastic medium is considered, containing a random set of ellipsoidal inhomogeneities with different poroelastic characteristics. The properties of these constituents are described by the linear porothermoelastic theory of Biot. One of the self-consistent schemes named effective field method is used to develop explicit expressions for the effective porothermoelastic constants (tensor of the frame elastic compliances , tensor of the generalized Skempton’s coefficients , tensor of thermal expansion coefficients , Biot’s constants , and the heat capacity at constant stress for the static porothermoelastic theory. It is shown that for two components composite porous material these expressions are interconnected and can be expressed only via the components of tensor . Some special cases are considered for the isotropic main material (matrix).     

Numerical simulation of vortex ring formation in the presence of background flow with implications for squid propulsion

Numerical simulations are used to study laminar vortex ring formation under the influence of background flow. The numerical setup includes a round-headed axisymmetric body with an opening at the posterior end from which a column of fluid is pushed out by a piston. The piston motion is explicitly included into the simulations by using a deforming mesh. A well-developed wake flow behind the body together with a finite-thickness boundary layer outside the opening is taken as the initial flow condition. As the jet is initiated, different vortex evolution behavior is observed depending on the combination of background flow velocity to mean piston velocity () ratio and piston stroke to opening diameter () ratio. For low background flow () with a short jet (), a leading vortex ring pinches off from the generating jet, with an increased formation number. For intermediate background flow () with a short jet (), a leading vortex ring also pinches off but with a reduced formation number. For intermediate background flow () with a long jet (), no vortex ring pinch-off is observed. For high background flow () with both a short () and a long () jet, the leading vortex structure is highly deformed with no single central axis of fluid rotation (when viewed in cross-section) as would be expected for a roll-up vortex ring. For , the vortex structure becomes isolated as the trailing jet is destroyed by the opposite-signed vorticity of the background flow. For , the vortex structure never pinches off from the trailing jet. The underlying mechanism is the interaction between the vorticity layer of the jet and the opposite-signed vorticity layer from the initial wake. This interaction depends on both and . A comparison is also made between the thrust generated by long, continuous jets and jet events constructed from a periodic series of short pulses having the same total mass flux. Force calculations suggest that long, continuous jets maximize thrust generation for a given amount of energy expended in creating the jet flow. The implications of the numerical results are discussed as they pertain to adult squid propulsion, which have been observed to generate long jets without a prominent leading vortex ring. PACS 02.60.Cb, 47.32.cf, 47.32.cb, 47.20.Ft, 47.63.M-     

Investigation of aeroacoustic noise generation by simultaneous particle image velocimetry and microphone measurements

A correlation technique is tested, which enables the identification of flow structures that are involved in sound generation processes. At first, the method is applied to the problem of induced noise from flow over a cylinder. The velocity field around a circular cylinder is measured by particle image velocimetry (PIV), while the radiated sound is recorded with a microphone. Both measurements are conducted in a synchronized manner so as to enable the calculation of the cross-correlation between velocity or vorticity fluctuations and the acoustic pressure. The therewith obtained coefficient matrix provides time- and space-resolved information about the statistical dependency between flow structures and the acoustic pressure. Furthermore, a proper orthogonal decomposition (POD) is applied to the velocity field. Then the correlation between dominating modes and the acoustic pressure is computed to identify which modes are mainly involved in the sound generation. Finally, the developed method is applied to the more applied problem of the flow-field inside a leading-edge slat-cove. The results show that, in this case, the signal-to-noise ratio is too low to allow an identification of noise-relevant flow structures, as opposed to the case of the cylinder wake flow, where 5,000 PIV recordings were sufficient to identify the flow structures, which are involved in the noise-generation process. A maximum in spatial distribution of the cross-correlation coefficient is observed 1.6 diameters downstream of the cylinder; its value decreases as one moves further downstream. In this area of maximal correlation, a rapid acceleration of the released vortices takes place. The cross-correlation coefficient fluctuates over time in a sine-type oscillation with maximum values of about and show a periodic behavior with a phase shift of π/2 with respect to each other. These regular oscillations can be explained by coherent periodic structures in the flow-field. These structures generate a sound field with the same periodicity, which is perceived as tone. Hence, the correlation between the velocity fluctuations and the acoustic pressure show oscillations identical to those of the input signals. A filtering of uncorrelated noise can be observed; this being caused by the averaging process during the cross-correlation calculation. The correlation with the eigenmodes of a POD gives correlation coefficients, which are no larger than the correlation with a local near-field quantity.     

Oscillating Solutions and Large ω-limit Sets in a Degenerate Parabolic Equation

The paper deals with positive solutions of the initial-boundary value problem for with zero Dirichlet data in a smoothly bounded domain . Here is positive on (0,∞) with f(0) = 0, and λ₁ is exactly the first Dirichlet eigenvalue of −Δ in Ω. In this setting, (*) may possess oscillating solutions in presence of a sufficiently strong degeneracy. More precisely, writing , it is shown that if then there exist global classical solutions of (*) satisfying and . Under the additional structural assumption , s > 0, this result can be sharpened: If then (*) has a global solution with its ω-limit set being the ordered arc that consists of all nonnegative multiples of the principal Laplacian eigenfunction. On the other hand, under the above additional assumption the opposite condition ensures that all solutions of (*) will stabilize to a single equilibrium.      

Large amplitude oscillatory elongation flow

A filament stretching rheometer (FSR) was used for measuring the elongation flow with a large amplitude oscillative elongation imposed upon the flow. The large amplitude oscillation imposed upon the elongational flow as a function of the time t was defined as where ε is the Hencky strain, is a constant elongational rate for the base elongational flow, Λ the strain amplitude (Λ ≥ 0), and Ω the strain frequency. A narrow molecular mass distribution linear polystyrene with a molecular weight of 145 kg/mol was subjected to the oscillative flow. The onset of the steady periodic regime is reached at the same Hencky strain as the onset of the steady elongational viscosity ( Λ = 0). The integral molecular stress function formulation within the ‘interchain pressure’ concept agrees qualitatively with the experiments.     

A New Low-Reynolds-Number One-Equation Model of Turbulence

In this study, we propose a new Low-Reynolds-Number (LRN)one-equation model, which is derived from an LRN two-equation(k-ε) model. The derivation of the transport equation, in principle, is based on the assumption that the turbulent structure parameter remains constant. However, the relation for the turbulent structure parameter a ₁(=|− |/k) is modified to account for near-wall turbulence. As a result, the present one-equation model contains a term which takes the near-wall limiting behavior explicitly into account. Thus, the present model provides the correct wall-limiting behavior of turbulence in the vicinity of the wall and can be applied to the analysis of heat transfer. The validity of the present model is tested in channel flows, boundary layer flows with and without pressure gradient, plane wall jet, and flow with separation and reattachment. The calculated results showed good agreement with the direct numerical simulation (DNS) and experimental data. This revised version was published online in July 2006 with corrections to the Cover Date.     

Effects of wall permeability on turbulence

In order to understand the effects of the wall permeability on turbulence near a porous wall, flow field measurements are carried out for turbulent flows in a channel with a porous bottom wall by a two-component particle image velocimetry (PIV) system. The porous media used are three kinds of foamed ceramics which have almost the same porosity (0.8) but different permeability. It is confirmed that the flow becomes more turbulent over the porous wall and tends to be turbulent even at the bulk Reynolds number of Re_b=1300 in the most permeable wall case tested. Corresponding to laminar to turbulent transition, the magnitude of the slip velocity on the porous wall is found to increase drastically in a narrow range of the Reynolds number. To discuss the effects of the wall roughness and the wall permeability, detailed discussions are made of zero-plane displacement and equivalent wall roughness for porous media. The results clearly indicate that the turbulence is induced by not only the wall roughness but the wall permeability. The measurements have also revealed that as Re_b or the wall permeability increases, the wall normal fluctuating velocity near the porous wall is enhanced due to the effects of the wall permeability. This leads to the increase of the turbulent shear stress resulting in higher friction factors of turbulence over porous walls.     

Large Eddy Simulation of Premixed Turbulent Combustion Using Ξ Flame Surface Wrinkling Model

One commonly-used method for deriving the RANS equations for multicomponent flow is the technique of conditional averaging. In this paper the concept is extended to LES, by introducing the operations of conditional filtering and surface filtering. Properties of the filtered indicator function are investigated mathematically and computationally. These techniques are then used to derive conditionally filtered versions of the Navier–Stokes equations which are appropriate for simulating multicomponent flow in LES. Transport equations for the favre-averaged indicator function and the unresolved interface properties (the wrinkling and the surface area per unit volume) are also derived. Since the paper is directed towards modelling premixed combustion in the flamelet regime, closure of the equations is achieved by introducing physical models based on the picture of the flame as a wrinkled surface separating burnt and unburnt components of the fluid. This leads to a set of models for premixed turbulent combustion of varying complexity. The results of applying one of this set of models to propagation of a spherical flame in isotropic homogeneous turbulence are analysed. This revised version was published online in July 2006 with corrections to the Cover Date.