scispace - formally typeset
Search or ask a question

Showing papers on "Open-channel flow published in 2010"


Journal ArticleDOI
TL;DR: In this article, a method is introduced that is able to extract dynamic information from flow fields that are either generated by a (direct) numerical simulation or visualized/measured in a physical experiment.
Abstract: The description of coherent features of fluid flow is essential to our understanding of fluid-dynamical and transport processes. A method is introduced that is able to extract dynamic information from flow fields that are either generated by a (direct) numerical simulation or visualized/measured in a physical experiment. The extracted dynamic modes, which can be interpreted as a generalization of global stability modes, can be used to describe the underlying physical mechanisms captured in the data sequence or to project large-scale problems onto a dynamical system of significantly fewer degrees of freedom. The concentration on subdomains of the flow field where relevant dynamics is expected allows the dissection of a complex flow into regions of localized instability phenomena and further illustrates the flexibility of the method, as does the description of the dynamics within a spatial framework. Demonstrations of the method are presented consisting of a plane channel flow, flow over a two-dimensional cavity, wake flow behind a flexible membrane and a jet passing between two cylinders.

4,150 citations


Journal ArticleDOI
TL;DR: In this paper, the authors describe the flow and deposition observed in and around a finite patch of vegetation located at the wall of a channel, and demonstrate that the deposition pattern is nearly uniform across the patch width and along the length of the patch.

243 citations


Journal ArticleDOI
TL;DR: In this paper, a direct numerical simulation of channel flow with rough walls comprising staggered arrays of cubes having various plan area densities is presented and discussed, and it is shown that the surface drag is predominantly form drag, which is greatest at an area coverage around 15%, the height variation of the axial pressure force across the obstacles weakens significantly as the area coverage decreases, but is always largest near the top of the obstacles.
Abstract: Computations of channel flow with rough walls comprising staggered arrays of cubes having various plan area densities are presented and discussed. The cube height h is 12.5% of the channel half-depth and Reynolds numbers (u? h/?) are typically around 700 – well into the fully rough regime. A direct numerical simulation technique, using an immersed boundary method for the obstacles, was employed with typically 35 million cells. It is shown that the surface drag is predominantly form drag, which is greatest at an area coverage around 15%. The height variation of the axial pressure force across the obstacles weakens significantly as the area coverage decreases, but is always largest near the top of the obstacles. Mean flow velocity and pressure data allow precise determination of the zero-plane displacement (defined as the height at which the axial surface drag force acts) and this leads to noticeably better fits to the log-law region than can be obtained by using the zero-plane displacement merely as a fitting parameter. There are consequent implications for the value of von K´arm´ an’s constant. As the effective roughness of the surface increases, it is also shown that there are significant changes to the structure of the turbulence field around the bottom boundary of the inertial sublayer. In distinct contrast to twodimensional roughness (longitudinal or transverse bars), increasing the area density of this three-dimensional roughness leads to a monotonic decrease in normalized vertical stress around the top of the roughness elements. Normalized turbulence stresses in the outer part of the flows are nonetheless very similar to those in smooth-wall flows.

204 citations


Journal ArticleDOI
TL;DR: In this article, the dispersion of rigid, highly elongated fibers in a turbulent channel flow is investigated, where fibers are treated as prolate ellipsoidal particles which move according to their inertia and to hydrodynamic drag.
Abstract: In this paper, the dispersion of rigid, highly elongated fibers in a turbulent channel flow is investigated. Fibers are treated as prolate ellipsoidal particles which move according to their inertia and to hydrodynamic drag and rotate according to hydrodynamic torques. The orientational behavior of fibers is examined together with their preferential distribution, near-wall accumulation, and wall deposition: all these phenomena are interpreted in connection with turbulence dynamics near the wall. In this work a wide range of fiber classes, characterized by different elongation (quantified by the fiber aspect ratio, λ) and different inertia (quantified by a suitably defined fiber response time, τp) is considered. A parametric study in the (λ,τp)-space confirms that, in the vicinity of the wall, fibers tend to align with the mean streamwise flow direction. However, this aligned configuration is unstable, particularly for higher inertia of the fiber, and can be maintained for rather short times before fibers ...

187 citations


Journal ArticleDOI
TL;DR: It is shown that, in a turbulent channel at Re{τ}≈550, large-scale motions can self-sustain even when smaller-scale structures populating the near-wall and logarithmic regions are artificially quenched.
Abstract: Large-scale motions, important in turbulent shear flows, are frequently attributed to the interaction of structures at smaller scales Here we show that, in a turbulent channel at Reτ≈550, large-scale motions can self-sustain even when smaller-scale structures populating the near-wall and logarithmic regions are artificially quenched This large-scale self-sustained mechanism is not active in periodic boxes of width smaller than Lz≈15h or length shorter than Lx≈3h which correspond well to the most energetic large scales observed in the turbulent channel

174 citations


Journal ArticleDOI
TL;DR: In this paper, a large-eddy simulation of turbulent flow in an open channel through staggered arrays of rigid, emergent cylinders, which can be regarded as idealized vegetation, is presented.
Abstract: This paper presents results of several large-eddy simulations (LES) of turbulent flow in an open channel through staggered arrays of rigid, emergent cylinders, which can be regarded as idealized vegetation. In this study, two cylinder Reynolds numbers, RD=1,340 and RD=500, and three vegetation densities are considered. The LES of the lowest density and at RD=1,340 corresponds to a recently completed laboratory experiment, the data of which is used to validate the simulations. Fairly good agreement between calculated and measured first- and second-order statistics along measurement profiles is found, confirming the accuracy of the simulations. The high resolution of the simulations enables an explicit calculation of drag forces, decomposed into pressure and friction drag, that are exerted on the cylinders. The effect of the cylinder Reynolds number and the cylinder density on the drag and hence on the flow resistance is quantified and in agreement with previous experimental studies. Turbulence structures are visualized through instantaneous pressure fluctuations, isosurfaces of the Q-criterion and contours of vertical vorticity in horizontal planes. Analysis of velocity time signals and distributions of drag and lift forces over time reveals that flow and turbulence are more influenced by the vegetation density than by the cylinder Reynolds number.

167 citations



Journal ArticleDOI
TL;DR: In this paper, the shape of the dominant large-scale structure is investigated by conditional averages based on the largescale velocity, determined using a filter width equal to the channel half-height.
Abstract: We investigate statistics of large-scale structures from large-eddy simulation (LES) of turbulent channel flow at friction Reynolds numbers $Re_\tau = 2 {\rm k}$ and $200 {\rm k}$. To properly capture the behaviour of large-scale structures, the channel length is chosen to be 96 times the channel half-height. In agreement with experiments, these large-scale structures are found to give rise to an apparent amplitude modulation of the underlying small-scale fluctuations. This effect is explained in terms of the phase relationship between the large- and small-scale activity. The shape of the dominant large-scale structure is investigated by conditional averages based on the large-scale velocity, determined using a filter width equal to the channel half-height. The conditioned field demonstrates coherence on a scale of several times the filter width, and the small-scale--large-scale relative phase difference increases away from the wall, passing through $\pi/2$ in the overlap region of the mean velocity before approaching $\pi$ further from the wall. We also found that, near the wall, the convection velocity of the large-scales departs slightly, but unequivocally, from the mean velocity.

161 citations


Journal ArticleDOI
TL;DR: In this paper, the shape of the dominant large-scale structure is investigated by conditional averages based on the large scale velocity, determined using a filter width equal to the channel half-height.
Abstract: We investigate statistics of large-scale structures from large-eddy simulation (LES) of turbulent channel flow at friction Reynolds numbers Re_τ = 2K and 200K (where K denotes 1000). In order to capture the behaviour of large-scale structures properly, the channel length is chosen to be 96 times the channel half-height. In agreement with experiments, these large-scale structures are found to give rise to an apparent amplitude modulation of the underlying small-scale fluctuations. This effect is explained in terms of the phase relationship between the large- and small-scale activity. The shape of the dominant large-scale structure is investigated by conditional averages based on the large-scale velocity, determined using a filter width equal to the channel half-height. The conditioned field demonstrates coherence on a scale of several times the filter width, and the small-scale–large-scale relative phase difference increases away from the wall, passing through π/2 in the overlap region of the mean velocity before approaching π further from the wall. We also found that, near the wall, the convection velocity of the large scales departs slightly, but unequivocally, from the mean velocity.

157 citations


Journal ArticleDOI
TL;DR: In this paper, the authors report on the pronounced turbulence modulations and the accompanying drag reduction observed in a two-way coupled simulation of particle-laden channel flow and support the view that drag reduction can be achieved not only by means of polymeric or fiber additives but also with spherical particles.
Abstract: This letter reports on the pronounced turbulence modulations and the accompanying drag reduction observed in a two-way coupled simulation of particle-laden channel flow. The present results support the view that drag reduction can be achieved not only by means of polymeric or fiber additives but also with spherical particles.

141 citations


Journal ArticleDOI
TL;DR: In this article, the vortical structures of near-wall turbulence at moderate Reynolds number are analyzed and compared in datasets obtained from stereoscopic particle image velocimetry (SPIV) in a turbulent boundary layer and from direct numerical simulations (DNS) in turbulent channel flow.
Abstract: The vortical structures of near-wall turbulence at moderate Reynolds number are analyzed and compared in datasets obtained from stereoscopic particle image velocimetry (SPIV) in a turbulent boundary layer and from direct numerical simulations (DNS) in a turbulent channel flow. The SPIV data are acquired in the LTRAC water tunnel at Re=+=820 in a streamwise/wall-normal plane, and in the LML wind tunnel at Re=+=2590 in a streamwise/wall-normal plane and in a plane orthogonal to the mean flow. The DNS data is taken from DelAlamo et al. (Self-similar vortex clusters in the turbulent logarithmic region, J. Fluid Mech. 561 (2006), pp. 329-358), at Re=950. The SPIV database is validated through an analysis of its mean velocity profile and power spectra, which are compared with reference profiles. A common detection algorithm is then applied to both the SPIV and DNS datasets in order to retrieve the streamwise and spanwise vortices. The algorithm employed is based on the 2D swirling strength and on a fit of an Oseen vortex model, which allows to retrieve the vortex characteristics, including its radius, vorticity, and position of the center. At all Reynolds numbers, the near-wall region is found to be the most densely populated region, predominantly with streamwise vortices that are on average smaller and more intense than spanwise vortices. In contrast, the logarithmic region is equally constituted of streamwise and spanwise vortices having equivalent characteristics. Two different scalings were employed to analyze the vortex radius and vorticity: the wall unit scaling and the Kolmogorov scaling. In wall unit scaling, a good universality in Reynolds numbers of the vortices radius and vorticity is observed in the near-wall and logarithmic regions: the vorticity is found to be maximum at the wall, decreasing first rapidly and then increasing slowly with wall-normal distance; the radius is increasing slowly with wall-normal distance in both the regions, except for the streamwise vortices, for which a sharp increase in radius is observed in the near-wall region. However, wall units scaling is found to be deficient in the outer region, where Reynolds number effects are observed. The Kolmogorov scaling appears to be universal both in Reynolds number and wall-normal distance across the investigated three regions, with a mean radius of the order of 8 and a mean vorticity of the order of 1.5-1, but for the SPIV data only. In the DNS dataset, the radius in the Kolmogorov scaling slowly decreases with wall-normal distance. This difference between the SPIV and DNS datasets may be linked to a difference between the boundary layer flow and the channel flow, rather than to the techniques themselves. Finally, the distribution of the vorticity of the detected vortices seems to follow faithfully a log-normal distribution, in good agreement with Kolmogorov's theory (A refinement of previous hypothesis concerning the local structure of turbulence in a viscous incompressible fluid at high Reynolds number, J. Fluid Mech. 13 (1962), pp. 82-85).

Journal ArticleDOI
TL;DR: In this paper, a direct numerical simulation of fully developed turbulent pipe flow is carried out at Reτ≈170 and 500 to investigate the effect of the streamwise periodic length on the convergence of turbulence statistics.
Abstract: In this paper, direct numerical simulation of fully developed turbulent pipe flow is carried out at Reτ≈170 and 500 to investigate the effect of the streamwise periodic length on the convergence of turbulence statistics. Mean flow, turbulence intensities, correlations, and energy spectra were computed. The findings show that in the near-wall region (below the buffer region, r+≤30), the required pipe length for all turbulence statistics to converge needs to be at least a viscous length of O(6300) wall units and should not be scaled with the pipe radius (δ). It was also found for convergence of turbulence statistics at the outer region that the pipe length has to be scaled with pipe radius and a proposed pipe length of 8πδ seems sufficient for the Reynolds numbers considered in this study.

Journal ArticleDOI
TL;DR: In this article, it was shown that the position and number of focusing positions depend on the number of particles per unit length along the channel, which is a function of both the channel cross-section and the particle volume fraction.
Abstract: Inertial focusing in a pressure-driven flow refers to the positioning of particles transverse to the mean flow direction that occurs as a consequence of a finite particle Reynolds number. In channels with rectangular cross-sections, and for a range of channel aspect ratios and particle confinement, experimental results are presented to show that both the location and the number of focusing positions depend on the number of particles per unit length along the channel. This axial number density is a function of both the channel cross-section and the particle volume fraction. These results are rationalized using simulations of the particle-laden flow to show the manner in which hydrodynamic interactions set the preferred locations in these confined flows. A criterion is presented for the occurrence of a stepwise transition from one to two or more trains of particles.

Journal ArticleDOI
TL;DR: In this paper, the authors proposed a nonlinear model for meander hydrodynamics without curvature restrictions, which reveals the mechanisms that drive the velocity redistribution in meander bends and their dependence on the river's roughness, the flow depth H, the radius of curvature R, the width B, and bathymetric variations.
Abstract: Despite the rapid evolution of computational power, simulation of meander dynamics by means of reduced and computationally less expensive models remains practically relevant for investigation of large?scale and long?term processes, probabilistic predictions, or rapid assessments. Existing meander models are invariantly based on the assumptions of mild curvature and slow curvature variations and fail to explain processes in the high?curvature range. This article proposes a nonlinear model for meander hydrodynamics without curvature restrictions. It provides the distribution of the main flow, the magnitude of the secondary flow, the direction of the bed shear stress, and the curvature?induced additional energy losses. It encompasses existing mild curvature models, remains valid for straight flow, and agrees satisfactorily with experimental data from laboratory experiments under conditions that are more demanding than sharp natural river bends. The proposed model reveals the mechanisms that drive the velocity redistribution in meander bends and their dependence on the river's roughness Cf, the flow depth H, the radius of curvature R, the width B, and bathymetric variations. It identifies Cf?1H/R as the major control parameter for meander hydrodynamics in general and the relative curvature R/B for sharp curvature effects. Both parameters are small in mildly curved bends but O(1) in sharply curved bends, resulting in significant differences in the flow dynamics. Streamwise curvature variations are negligible in mildly curved bends, but they are the major mechanisms for velocity redistribution in sharp bends. Nonlinear feedback between the main and secondary flow also plays a dominant role in sharp bends: it increases energy losses and reduces the secondary flow, the transverse bed slope, and the velocity redistribution.

Journal ArticleDOI
TL;DR: A conservative, second-order accurate immersed interface method for representing incompressible fluid flows over complex three dimensional solid obstacles on a staggered Cartesian grid, suitable for Large-Eddy Simulations of high-Reynolds number flows.

Journal ArticleDOI
TL;DR: In this paper, an analysis of the average flow pattern in the local coordinate system defined by the eigenvectors of the strain rate tensor is presented, showing that the directions of maximum and minimum strain in a shear layer are in the plane normal to the vorticity vector, this vector aligns with the remaining strain direction, i.e. the intermediate eigenvector of the tensor.
Abstract: Two aspects of small-scale turbulence are currently regarded universal, as they have been reported for a wide variety of turbulent flows. Firstly, the vorticity vector has been found to display a preferential alignment with the eigenvector corresponding to the intermediate eigenvalue of the strain rate tensor; and secondly, the joint probability density function (p.d.f.) of the second and third invariant of the velocity gradient tensor, Q and R, has a characteristic teardrop shape. This paper provides an explanation for these universal aspects in terms of a spatial organization of coherent structures, which is based on an evaluation of the average flow pattern in the local coordinate system defined by the eigenvectors of the strain rate tensor. The approach contrasts with previous investigations, which have relied on assumed model flows. The present average flow patterns have been calculated for existing experimental (particle image velocimetry) or numerical (direct numerical simulation) datasets of a turbulent boundary layer (TBL), a turbulent channel flow and for homogeneous isotropic turbulence. All results show a shear-layer structure consisting of aligned vortical motions, separating two larger-scale regions of relatively uniform flow. Because the directions of maximum and minimum strain in a shear layer are in the plane normal to the vorticity vector, this vector aligns with the remaining strain direction, i.e. the intermediate eigenvector of the strain rate tensor. Further, the QR joint p.d.f. for these average flow patterns reveals a shape reminiscent of the teardrop, as seen in many turbulent flows. The above-mentioned organization of the small-scale motions is not only found in the average patterns, but is also frequently observed in the instantaneous velocity fields of the different turbulent flows. It may, therefore, be considered relevant and universal.

Journal ArticleDOI
TL;DR: An unstructured hybrid mesh numerical method is developed to simulate open channel flows and offers a framework to unify many mesh topologies into a single formulation.
Abstract: An unstructured hybrid mesh numerical method is developed to simulate open channel flows. The method is applicable to arbitrarily shaped mesh cells and offers a framework to unify many mesh topologies into a single formulation. A finite-volume discretization is applied to the two-dimensional depth-averaged equations such that mass conservation is satisfied both locally and globally. An automatic wetting-drying procedure is incorporated in conjunction with a segregated solution procedure that chooses the water surface elevation as the main variable. The method is applicable to both steady and unsteady flows and covers the entire flow range: subcritical, transcritical, and supercritical. The proposed numerical method is well suited to natural river flows with a combination of main channels, side channels, bars, floodplains, and in-stream structures. Technical details of the method are presented, verification studies are performed using a number of simple flows, and a practical natural river is modeled to il...

Journal ArticleDOI
TL;DR: Numerical simulations indicate that the present formulation is second order accurate and the difference of adopting different local known energy distribution functions is, as expected, negligible, which are consistent with the results from the derived discrete macroscopic energy equation.
Abstract: Consistent 2D and 3D thermal boundary conditions for thermal lattice Boltzmann simulations are proposed. The boundary unknown energy distribution functions are made functions of known energy distribution functions and correctors, where the correctors at the boundary nodes are obtained directly from the definition of internal energy density. This boundary condition can be easily implemented on the wall and corner boundary using the same formulation. The discrete macroscopic energy equation is also derived for a steady and fully developed channel flow to assess the effect of the boundary condition on the solutions, where the resulting second order accurate central difference equation predicts continuous energy distribution across the boundary, provided the boundary unknown energy distribution functions satisfy the macroscopic energy level. Four different local known energy distribution functions are experimented with to assess both this observation and the applicability of the present formulation, and are scrutinized by calculating the 2D thermal Poiseuille flow, thermal Couette flow, thermal Couette flow with wall injection, natural convection in a square cavity, and 3D thermal Poiseuille flow in a square duct. Numerical simulations indicate that the present formulation is second order accurate and the difference of adopting different local known energy distribution functions is, as expected, negligible, which are consistent with the results from the derived discrete macroscopic energy equation.

Journal ArticleDOI
TL;DR: A two-layer approach for depth-limited open-channel flow with submerged vegetation is described in this article, where a momentum balance is applied to each layer and expressions for the mean velocities are proposed.
Abstract: A two-layer approach for depth-limited open-channel flow with submerged vegetation is described. A momentum balance is applied to each layer and expressions for the mean velocities are proposed. The velocity is assumed to be uniform in the vegetation layer and logarithmic in the upper layer. The proposed relationship successfully predicts the mean velocity distribution when compared with the measured data. Using the velocity formula, the layer-averaged mean velocities in the upper layer and over the entire layer are derived. An expression for the roughness coefficient increased by vegetation is also presented, performing better for the roughness coefficient than other formulas. Another relationship is proposed for predicting the distribution of suspended sediment in depth-limited flow with submerged vegetation by using an eddy-viscosity profile. The predicted profiles moderately agree with the measured data. Comparisons with simulated data from the Reynolds-averaged Navier–Stokes equations with the k–ϵ mo...

Journal ArticleDOI
TL;DR: In this article, an analysis of published data on a gas-liquid two-phase flow in capillaries of various shapes is given, which makes it possible to systematize the collected body of information.
Abstract: Capillary hydrodynamics has three considerable distinctions from macrosystems: first, there is an increase in the ratio of the surface area of the phases to the volume that they occupy; second, a flow is characterized by small Reynolds numbers at which viscous forces predominate over inertial forces; and third, the microroughness and wettability of the wall of the channel exert a considerable influence on the flow pattern. In view of these differences, the correlations used for tubes with a larger diameter cannot be used to calculate the boundaries of the transitions between different flow regimes in microchannels. In the present review, an analysis of published data on a gas-liquid two-phase flow in capillaries of various shapes is given, which makes it possible to systematize the collected body of information. The specific features of the geometry of a mixer and an inlet section, the hydraulic diameter of a capillary, and the surface tension of a liquid exert the strongest influence on the position of the boundaries of two-phase flow regimes. Under conditions of the constant geometry of the mixer, the best agreement in the position of the boundaries of the transitions between different hydrodynamic regimes in capillaries is observed during the construction of maps of the regimes with the use of the Weber numbers for a gas and a liquid as coordinate axes.

Journal ArticleDOI
TL;DR: Direct measurements of the wavenumber-frequency spectrum of wall pressure fluctuations beneath a turbulent plane channel flow have been performed in an anechoic wind tunnel, showing that for aerodynamically induced car interior noise, both contributions to the surface pressure fluctuations on car windows have to be taken into account.
Abstract: Direct measurements of the wavenumber-frequency spectrum of wall pressure fluctuations beneath a turbulent plane channel flow have been performed in an anechoic wind tunnel. A rotative array has been designed that allows the measurement of a complete map, 63×63 measuring points, of cross-power spectral densities over a large area. An original post-processing has been developed to separate the acoustic and the aerodynamic exciting loadings by transforming space-frequency data into wavenumber-frequency spectra. The acoustic part has also been estimated from a simple Corcos-like model including the contribution of a diffuse sound field. The measured acoustic contribution to the surface pressure fluctuations is 5% of the measured aerodynamic surface pressure fluctuations for a velocity and boundary layer thickness relevant for automotive interior noise applications. This shows that for aerodynamically induced car interior noise, both contributions to the surface pressure fluctuations on car windows have to be taken into account.

Journal ArticleDOI
TL;DR: The LES WALE model, which can be used without wall functions or global damping functions, is investigated within the lattice Boltzmann framework, which produces an efficient and fast scheme due to its algebraic character.
Abstract: Turbulence models which can perform the transition from laminar flow to fully developed turbulent flow are of key importance in industrial applications. A promising approach is the LES WALE model, which can be used without wall functions or global damping functions. The model produces an efficient and fast scheme due to its algebraic character. Additionally, its prediction of the transition from laminar to turbulent regimes has shown promising results. In this work, the LES WALE model is investigated within the lattice Boltzmann framework. For validation purposes, various test cases are presented. First, a channel flow at a Reynolds number of 6876 is investigated. Secondly, the flow around a wall-mounted cube at various Reynolds numbers is determined. The flow regime varies from laminar, to transitional, to fully turbulent conditions at a Reynolds number of 40,000 with respect to the cube height.

Journal ArticleDOI
TL;DR: In this paper, a quasi-counter flow parallel-plates total heat exchanger is constructed and investigated, and a detailed mathematical modeling is conducted and the model is experimentally verified, and the fluid flow, heat and mass transport equations in the entry regions are solved directly.

Journal ArticleDOI
TL;DR: In this article, the Homogenous Relaxation Model (HRM) is used to study thermal non-equilibrium, two-phase flows with flash-boiling and condensation.

Journal ArticleDOI
TL;DR: In this paper, the authors simulate turbulent channel flow of drag-reducing polymer solutions in minimal flow geometries and find intervals of "hibernating" turbulence that display many features of the universal maximum drag reduction asymptote observed in polymer solutions: weak streamwise vortices, nearly nonexistent streamwise variations, and a mean velocity gradient.
Abstract: Turbulent channel flow of drag-reducing polymer solutions is simulated in minimal flow geometries. Even in the Newtonian limit, we find intervals of "hibernating" turbulence that display many features of the universal maximum drag reduction asymptote observed in polymer solutions: weak streamwise vortices, nearly nonexistent streamwise variations, and a mean velocity gradient that quantitatively matches experiments. As viscoelasticity increases, the frequency of these intervals also increases, while the intervals themselves are unchanged, leading to flows that increasingly resemble maximum drag reduction.

Journal ArticleDOI
TL;DR: In this paper, a wall-resolving large eddy simulation (LES) is used to analyze turbulent flow in a channel with irregular two-dimensional rough surfaces, and the effects of roughness on the turbulent flow are greatly related to the increase of the height of the maximum peaks of the corrugations.
Abstract: Turbulent flow in a channel with irregular two-dimensional rough surfaces is analysed through wall-resolving large eddy simulation (LES). Both walls of the channel are roughened through the superimposition of sinusoidal functions having random amplitude and four different wavelengths. The downward shift of the velocity profile in the log region due to the roughness, known as roughness function, is well captured in the simulations. The spanwise and wall-normal turbulence intensities are found to increase with the roughness height, while the streamwise component decreases. The analysis of the Reynolds stress anisotropy tensor highlights a tendency towards isotropisation, confirmed by the vorticity rms. The analysis of the statistics shows that the effects of the roughness on the turbulent flow are greatly related to the increase of the height of the maximum peaks of the corrugations. Although the inner layer is dramatically affected by the wall irregularities, the outer layer appears not affected by the spe...

Journal ArticleDOI
TL;DR: In this article, it was shown that if a layer consists of partially molten felsic rocks, then conductances in the range 7000-27,000 s will produce flow velocities of the order 1.5 cm/a.

Journal ArticleDOI
TL;DR: In this paper, the authors measured the streamwise component of the turbulence intensity in a fully developed pipe flow at Reynolds numbers up to 145 000 and showed that the magnitude of the near-wall peak is invariant with Reynolds number in location and magnitude.
Abstract: New measurements of the streamwise component of the turbulence intensity in a fully developed pipe flow at Reynolds numbers up to 145 000 indicate that the magnitude of the near-wall peak is invariant with Reynolds number in location and magnitude. The results agree with previous pipe flow data that have sufficient spatial resolution to avoid spatial filtering effects, but stand in contrast to similar results obtained in boundary layers, where the magnitude of the peak displays a prominent Reynolds number dependence, although its position is fixed at the same location as in pipe flow. This indicates that the interaction between the inner and outer regions is different in pipe flows and boundary layers.

Journal ArticleDOI
Abstract: The mutual interaction of laminar-turbulent transition and mean flow evolution is studied in a pressure-induced laminar separation bubble on a flat plate. The flat-plate boundary layer is subjected to a sufficiently strong adverse pressure gradient that a separation bubble develops. Upstream of the bubble a small-amplitude disturbance is introduced which causes transition. Downstream of transition, the mean flow strongly changes and, due to viscous-inviscid interaction, the overall pressure distribution is changed as well. As a consequence, the mean flow also changes upstream of the transition location. The difference in the mean flow between the forced and the unforced flows is denoted the mean flow deformation. Two different effects are caused by the mean flow deformation in the upstream, laminar part: a reduction of the size of the separation region and a stabilization of the flow with respect to small, linear perturbations. By carrying out numerical simulations based on the original base flow and the time-averaged deformed base flow, we are able to distinguish between direct and indirect nonlinear effects. Direct effects are caused by the quadratic nonlinearity of the Navier-Stokes equations, are associated with the generation of higher harmonics and are predominantly local. In contrast, the stabilization of the flow is an indirect effect, because it is independent of the Reynolds stress terms in the laminar region and is solely governed by the non-local alteration of the mean flow via the pressure. © 2010 Cambridge University Press.

Journal ArticleDOI
TL;DR: In this article, numerical simulations were performed for three-dimensional two-phase water/oxygen flow in the flow field plate at the anode side of a PEM electrolysis cell.