Compound wall treatment with low‐Re turbulence model

A generalized treatment for the wall boundary conditions relating to turbulent flows is developed that blends the integration to a solid wall with wall functions. The blending function ensures a smooth transition between the viscous and turbulent regions. An improved low Reynolds number k?ε model is coupled with the proposed compound wall treatment to determine the turbulence field. The eddy viscosity formulation maintains the positivity of normal Reynolds stresses and Schwarz' inequality for turbulent shear stresses. The model coefficients/functions preserve the anisotropic characteristics of turbulence. Computations with fine and coarse meshes of a few flow cases yield appreciably good agreement with the direct numerical simulation and experimental data. The method is recommended for computing the complex flows where computational grids cannot satisfy a priori the prerequisites of viscous/turbulence regions. Copyright © 2011 John Wiley & Sons, Ltd.     

Friction relaxation model for fast transient flows application to water hammer in two-phase flow – The WAHA code

The paper deals with the problem of the wall shear stress during rapid transient 1D flows in a piping system caused by water hammers in two-phase flow induced by a fast valve closure. The evolution of the transient wall shear stress is interpreted in terms of two steps. The first step is a sudden and dramatic change of the wall shear stress due to the passage of the pressure wave. The second step is a relaxation process of the shear stress which is modeled from the Extended Irreversible Thermodynamics theory. The friction relaxation model (FRM) presented in the first part of this paper describes both steps of the evolution of the wall shear stress during water hammers. The second part of the paper deals with the application of the FRM model as a closure law in the WAHA code. The main purpose of the WAHA code is to predict various situations relative to single- and two-phase water hammer transients in piping systems. The last part of the paper deals with the simulation of several cases from the UMSICHT databank using the adapted WAHA computer code with the FRM model. The results of these simulations are systematically compared with the experimental data. It is concluded that the new FRM model has a clear effect on water hammer pressure wave damping and on the pressure wave propagation velocity.     

Numerical investigation of fully developed turbulent fluid flow and heat transfer in a square duct

Three-dimensional fully developed turbulent fluid flow and heat transfer in a square duct are numerically investigated with the author's anisotropic low-Reynolds-number k-ε turbulence model. Special attenton has been given to the regions close to the wall and the corner, which are known to influence the characteristics of secondary flow a great deal. Hence, instead of the common wall function approach, the no-slip boundary condition at the wall is directly used. Velocity and temperature profiles are predicted for fully developed turbulent flows with constant wall temperature. The predicted variations of both local wall shear stress and local wall heat flux are shown to be in close agreement with available experimental data. The present paper also presents the budget of turbulent kinetic energy equation and the systematic evaluation for existing wall function forms. The commonly adopted wall function forms that are valid for two-dimensional flows are found to be inadequate for three-dimensional turbulent flows in a square duct.     

Predictions of axisymmetric and two-dimensional impinging turbulent jets

The k − turbulence model and a version of a second-moment closure, modified to include the effect of pressure reflections from a solid surface, have been used as the basis of predictions of the flow that results from the orthogonal impingement of circular and two-dimensional (2-D) jets on a flat surface. Comparison of model predictions has been made with velocity measurements obtained in the stagnation and wall jet regions of the impinging flows. Results, in general, confirm the superiority of the Reynolds stress transport equation model for predicting mean and fluctuating velocities within the latter regions of such flows. In particular, modifications to the second-moment closure to account for the influence of the surface in distorting the fluctuating pressure field away from the wall successfully predict the damping of normal-to-wall velocity fluctuations throughout the impinging flows. In contrast, results derived from the eddy-viscosity-based approach do not, in general, accurately reproduce experimental observations.     

CFD Analysis of a Thermal Plume and the Indoor Air Flow Using k-ε Models with Buoyancy Effects

Turbulence models based on the eddy viscosity concept perform poorly for simulation of non-isothermal flows, which are characterized by strongly non-isotropic features of turbulence. In the present study, a thermal plume is calculated using the WET model and severalk-ε models. Comparison with the experiments reveals that the WET model predicts the flow and temperature fields most accurately. Furthermore, the WET model has been modified to account for the damping effect near the wall and applied to the indoor air flow. The computed flow field agrees well with the measurement, except for higher temperatures appearing near the walls. This revised version was published online in July 2006 with corrections to the Cover Date.     

Three‐dimensional numerical simulation of compound meandering open channel flow by the Reynolds stress model

Turbulent flow in a compound meandering open channel with seminatural cross sections is one of the most complicated turbulent flows as the flow pattern is influenced by the combined action of various forces, such as centrifugal force, pressure, and shear stresses. In this paper, a three‐dimensional (3D) Reynolds stress model (RSM) is adopted to simulate the compound meandering channel flows. Governing equations of the flow are solved numerically with finite‐volume method. The velocity fields, wall shear stresses, and Reynolds stresses are calculated for a range of input conditions. Good agreement between the simulated results and measurements indicates that RSM can successfully predict the complicated flow phenomenon. Copyright © 2008 John Wiley & Sons, Ltd.     

Wall shear stress and void fraction in Poiseuille bubbly flows: Part II: experiments and validity of analytical predictions

In a companion paper, a simple analytical formulation has been established which provides the wall shear stress in laminar bubbly flows for idealised transverse void fraction distributions. In the present paper, this approach is applied to Poiseuille bubbly flows in circular ducts. New measurements of the void fraction profiles and wall friction angular distribution in a pipe are presented for a wide range of flow parameters. Approximating the void profiles by step-functions allows us to evaluate the wall friction with the above mentioned model. Results are shown to agree satisfactorily with measurements. Notably, negative wall shear stress and wall shear stress much higher than their single-phase flow counterpart at the same liquid flow rate are recovered. Therefore, the principal mechanisms responsible for friction modification are captured with this simple model.     

Numerical study of effects of pulsatile amplitude for transitional turbulent pulsatile flow in pipes with ring‐type constrictions

The effects of pulsatile amplitude on sinusoidal transitional turbulent flows through a rigid pipe in the vicinity of a sharp‐edged mechanical ring‐type constriction have been studied numerically. Pulsatile flows were studied for transitional turbulent flow with Reynolds number (Re) of the order of 10⁴, Womersley number (Nw) of the order of 50 with a corresponding Strouhal number (St) of the order of 0.04. The pulsatile flow considered is a sinusoidal flow with dimensionless amplitudes varying from 0.0 to 1.0. Transitional laminar and turbulent flow characteristics in an alternative manner within the pulsatile flow fields were observed and studied numerically. The flow characteristics were studied through the pulsatile contours of streamlines, vorticity, shear stress and isobars. It was observed that fluid accelerations tend to suppress the development of flow disturbances. All the instantaneous maximum values of turbulent kinetic energy, turbulent viscosity, turbulent shear stress are smaller during the acceleration phase when compared with those during deceleration period. Various parametric equations within a pulsatile cycle have also been formulated through numerical experimentations with different pulsatile amplitudes. In the vicinity of constrictions, the empirical relationships were obtained for the instantaneous flow rate (Q), the pressure gradient (dp/dz), the pressure loss (P_loss), the maximum velocity (V_max), the maximum vorticity (ζ_max), the maximum wall vorticity (ζ_w,max), the maximum shear stress (τ_max) and the maximum wall shear stress (τ_w,max). Elliptic relation was observed between flow rate and pressure gradient. Quadratic relations were observed between flow rate and the pressure loss, the maximum values of shear stress, wall shear stress, turbulent kinematic energy and the turbulent viscosity. Linear relationships exist between the instantaneous flow rate and the maximum values of vorticity, wall vorticity and velocity. The time‐average axial pressure gradient and the time average pressure loss across the constriction were observed to increase linearly with the pulsatile amplitude. Copyright © 1999 John Wiley & Sons, Ltd.     

Predicting Buoyant Shear Flows Using Anisotropic Dissipation Rate Models

This paper examines the modeling of two-dimensional homogeneous stratified turbulent shear flows using the Reynolds-stress and Reynolds-heat-flux equations. Several closure models have been investigated; the emphasis is placed on assessing the effect of modeling the dissipation rate tensor in the Reynolds-stress equation. Three different approaches are considered; one is an isotropic approach while the other two are anisotropic approaches. The isotropic approach is based on Kolmogorov's hypothesis and a dissipation rate equation modified to account for vortex stretching. One of the anisotropic approaches is based on an algebraic representation of the dissipation rate tensor, while another relies on solving a modeled transport equation for this tensor. In addition, within the former anisotropic approach, two different algebraic representations are examined; one is a function of the Reynolds-stress anisotropy tensor, and the other is a function of the mean velocity gradients. The performance of these closure models is evaluated against experimental and direct numerical simulation data of pure shear flows, pure buoyant flows and buoyant shear flows. Calculations have been carried out over a range of Richardson numbers (Ri) and two different Prandtl numbers (Pr); thus the effect of Pr on the development of counter-gradient heat flux in a stratified shear flow can be assessed. At low Ri, the isotropic model performs well in the predictions of stratified shear flows; however, its performance deteriorates as Ri increases. At high Ri, the transport equation model for the dissipation rate tensor gives the best result. Furthermore, the results also lend credence to the algebraic dissipation rate model based on the Reynolds stress anisotropy tensor. Finally, it is found that Pr has an effect on the development of counter-gradient heat flux. The calculations show that, under the action of shear, counter-gradient heat flux does not occur even at Ri = 1 in an air flow. This revised version was published online in July 2006 with corrections to the Cover Date.     

Turbulence in a skewed three‐dimensional wall‐bounded shear flow: effect of mean vorticity on structure modification

Mean‐flow three‐dimensionalities affect both the turbulence level and the coherent flow structures in wall‐bounded shear flows. A tailor‐made flow configuration was designed to enable a thorough investigation of moderately and severely skewed channel flows. A unidirectional shear‐driven plane Couette flow was skewed by means of an imposed spanwise pressure gradient. Three different cases with 8°, 34°and 52°skewing were simulated numerically and the results compared with data from a purely two‐dimensional plane Couette flow. The resulting three‐dimensional flow field became statistically stationary and homogeneous in the streamwise and spanwise directions while the mean velocity vector V and the mean vorticity vector Ω remained parallel with the walls. Mean flow profiles were presented together with all components of the Reynolds stress tensor. The mean shear rate in the core region gradually increased with increasing skewing whereas the velocity fluctuations were enhanced in the spanwise direction and reduced in the streamwise direction. The Reynolds shear stress is known to be closely related to the coherent flow structures in the near‐wall region. The instantaneous and ensemble‐averaged flow structures were turned by the skewed mean flow. We demonstrated for the medium‐skewed case that the coherent structures should be examined in a coordinate system aligned with V to enable a sound interpretation of 3D effects. The conventional symmetry between Case 1 and Case 2 vortices was broken and Case 1 vortices turned out to be stronger than Case 2. This observation is in conflict with the common understanding on the basis of the spanwise (secondary) mean shear rate. A refined model was proposed to interpret the structure modifications in three‐dimensional wall‐flows. What matters is the orientation of the mean vorticity vector Ω relative to the vortex vorticity vector ω _v, that is, the sign of Ω · ω _v. In the present situation, Ω · ω _v > 0 for the Case 1 vortices causing a strengthening relative to the Case 2 vortices. Copyright © 2011 John Wiley & Sons, Ltd.     

Assessment of nonlinear strain measures for extensional and shearing flows of polymer melts

From stress-strain experiments in extensional and shearing flows, nonlinear strain measures and effective damping functions are derived for a polyisobutylene melt. The strain measures determined in planar extensional flow and in simple shear flow coincide. Experimental results are compared with predictions of two molecular theories, the Doi-Edwards model and the molecular stress function approach of Wagner and Schaeffer. Discrepancies between theories and experiment lead to a reconsideration of the classification of extensional flows. The symmetry of the flow field is identified and quantified as an important parameter influencing the strain measure, and a unifying strain measure for general extensional and shearing flows of polymer melts is presented.     

Wall shear stress modified by bubbles in a horizontal channel flow of silicone oil in the transition region

We performed laboratory experiments on bubbly channel flows using silicone oil, which has a low surface tension and clean interface to bubbles, as a test fluid to evaluate the wall shear stress modification for different regimes of bubble migration status. The channel Reynolds numbers of the flow ranged from 1000 to 5000, covering laminar, transition and turbulent flow regimes. The bubble deformation and swarms were classified as packing, film, foam, dispersed, and stretched states based on visualization of bubbles as a bulk void fraction changed. In the dispersed and film states, the wall shear stress reduced by 9% from that in the single-phase condition; by contrast, the wall shear stress increased in the stretched, packing, and foam states. We carried out statistical analysis of the time-series of the wall shear stress in the transition and turbulent-flow regimes. Variations of the PDF of the shear stress and the higher order moments in the statistic indicated that the injection of bubbles generated pseudo-turbulence in the transition regime and suppressed drag-inducing events in the turbulent regime. Bubble images and measurements of shear stress revealed a correlated wave with a time lag, for which we discuss associated to the bubble dynamics and effective viscosity of the bubble mixture in wall proximity.     

An immersed boundary/wall modeling method for RANS simulation of compressible turbulent flows

In this paper, an immersed boundary (IB) method is developed to simulate compressible turbulent flows governed by the Reynolds‐averaged Navier‐Stokes equations. The flow variables at the IB nodes (interior nodes in the immediate vicinity of the solid wall) are evaluated via linear interpolation in the normal direction to close the discrete form of the governing equations. An adaptive wall function and a 2‐layer wall model are introduced to reduce the near‐wall mesh density required by the high resolution of the turbulent boundary layers. The wall shear stress modified by the wall modeling technique and the no‐penetration condition are enforced to evaluate the velocity at an IB node. The pressure and temperature at an IB node are obtained via the local simplified momentum equation and the Crocco‐Busemann relation, respectively. The SST k ? ω and S‐A turbulence models are adopted in the framework of the present IB approach. For the Shear‐Stress Transport (SST) k ? ω model, analytical solutions in near‐wall region are utilized to enforce the boundary conditions of the turbulence equations and evaluate the turbulence variables at an IB node. For the S‐A model, the turbulence variable at an IB node is calculated by using the near‐wall profile of the eddy viscosity. In order to validate the present IB approach, numerical experiments for compressible turbulent flows over stationary and moving bodies have been performed. The predictions show good agreements with the referenced experimental data and numerical results.     

Experiments on the linear instability of flow in a wavy channel

The effects of wall corrugation on the stability of wall-bounded shear flows have been examined experimentally in plane channel flows. One of the channel walls has been modified by introduction of the wavy wall model with the amplitude of 4% of the channel half height and the wave number of 1.02. The experiment is focused on the two-dimensional travelling wave instability and the results are compared with the theory [J.M. Floryan, Two-dimensional instability of flow in a rough channel, Phys. Fluids 17 (2005) 044101 (also: Rept. ESFD-1/2003, Dept. of Mechanical and Materials Engineering, The University of Western Ontario, London, Ontario, Canada, 2003)]. It is shown that the flow is destabilized by the wall corrugation at subcritical Reynolds numbers below 5772, as predicted by the theory. For the present corrugation geometry, the critical Reynolds number is decreased down to about 4000. The spatial growth rates, the disturbance wave numbers and the distribution of disturbance amplitude measured over such wavy wall also agree well with the theoretical results.     

A general one-equation turbulence model for free shear and wall-bounded flows

The purpose of thiswork is to introduce a complete and general one-equation model capable of correctly predicting a wide class of fundamental turbulent flows like boundary layer, wake, jet, and vortical flows. The starting point is the mature and validated two-equation k−ω turbulence model of Wilcox. The newly derived one-equation model has several advantages and yields better predictions than the Spalart-Allmaras model for jet and vortical flows while retaining the same efficiency and quality of the results for near-wall turbulent flows without using a wall distance. The derivation and validation of the new model using findings computed by the Spalart-Allmaras and the k−ω models are presented and discussed for several free shear and wall-bounded flows.     

弯曲动脉的血流动力学数值分析

利用计算流体力学的理论和方法对弯曲动脉中的血流动力学进行数值分析，是研究心血管疾病流体动力学机理的一种行之有效的方法。本文将升主动脉、主动脉弓和降主动脉联系起来作为弯曲动脉几何模型，给出了血液流动的边界条件以及计算条件。根据生理脉动流条件，对狗的弯曲动脉几何模型内发展中的血液流动进行了有限元数值模拟，并利用可视化方法对血液流动的轴向速度、二次流、壁面切应力等计算结果进行了分析。研究结果表明，在弯管内侧壁处，同时存在主流方向和二次流方向的回流，此处容易形成涡流。弯管内侧壁比外侧壁的壁面切应力具有更强的脉动性。     

Optimal control of unsteady compressible viscous flows

The control of complex, unsteady flows is a pacing technology for advances in fluid mechanics. Recently, optimal control theory has become popular as a means of predicting best case controls that can guide the design of practical flow control systems. However, most of the prior work in this area has focused on incompressible flow which precludes many of the important physical flow phenomena that must be controlled in practice including the coupling of fluid dynamics, acoustics, and heat transfer. This paper presents the formulation and numerical solution of a class of optimal boundary control problems governed by the unsteady two‐dimensional compressible Navier–Stokes equations. Fundamental issues including the choice of the control space and the associated regularization term in the objective function, as well as issues in the gradient computation via the adjoint equation method are discussed. Numerical results are presented for a model problem consisting of two counter‐rotating viscous vortices above an infinite wall which, due to the self‐induced velocity field, propagate downward and interact with the wall. The wall boundary control is the temporal and spatial distribution of wall‐normal velocity. Optimal controls for objective functions that target kinetic energy, heat transfer, and wall shear stress are presented along with the influence of control regularization for each case. Copyright © 2002 John Wiley & Sons, Ltd.     

RAS one-equation turbulence model with non-singular eddy-viscosity coefficient

A simplified consistency formulation for P_k/ε (production to dissipation ratio) is devised to obtain a non-singular C_μ (coefficient of eddy-viscosity) in the explicit algebraic Reynolds stress model of Gatski and Speziale. The coefficient C_μ depends non-linearly on both rotational/irrotational strains and is used in the framework of an improved RAS (Rahman–Agarwal–Siikonen) one-equation turbulence model to calculate a few well-documented turbulent flows, yielding predictions in good agreement with the direct numerical simulation and experimental data. The strain-dependent C_μ assists the RAS model in constructing the coefficients and functions such as to benefit complex flows with non-equilibrium turbulence. Comparisons with the Spalart–Allmaras one-equation model and the shear stress transport k-ω model demonstrate that C_μ improves the response of RAS model to non-equilibrium effects.