This work explores simple relations that follow from the momentum-integral equation absent a pressure gradient. The resulting expressions enable us to relate the boundary-layer characteristics of a velocity profile, u(y), to an assumed flow function and its wall derivative relative to the wall-normal coordinate, y. Consequently, disturbance, displacement, and momentum thicknesses, as well as skin friction and drag coefficients, which are typically evaluated and tabulated in classical monographs, can be readily determined for a given profile, F(ξ) = u/U. Here, ξ = y/δ denotes the boundary-layer coordinate. These expressions are then employed to provide a rational explanation for the 1921 Pohlhausen polynomial paradox, namely, the reason why a quartic representation of the velocity leads to less accurate predictions of the disturbance, displacement, and momentum thicknesses than using cubic or quadratic polynomials. Not only do we identify the factors underlying this behavior but also we proceed to outline a procedure to overcome its manifestation at any order. This enables us to derive optimal piecewise approximations that do not suffer from the particular limitations affecting Pohlhausen’s F = 2ξ − 2ξ3 + ξ4. For example, our alternative profile, F = (5ξ − 3ξ3 + ξ4)/3, leads to an order-of-magnitude improvement in precision when incorporated into the Kármán–Pohlhausen approach in both viscous and thermal analyses. Then, noting the significance of the Blasius constant, s¯≈1.630 398, this approach is extended to construct a set of uniformly valid solutions, including F=1−exp[−s¯ξ(1+12s¯ξ+ξ2)], which continues to hold beyond the boundary-layer edge as y → ∞. Given its substantially reduced error, the latter is shown, through comparisons to other models, to be practically equivalent to the Blasius solution.
more »
« less
On the Kármán momentum-integral approach and the Pohlhausen paradox: Extension to a cylinder in crossflow with a potential farfield motion
In this work, the Kármán–Pohlhausen (KP) momentum-integral approach based on optimized fourth-order (MX4) polynomial approximations of the velocity and temperature profiles is applied to a classical benchmark problem, namely, that of a cylinder in crossflow with a variable pressure gradient. This enables us to extract closed-form expressions for both hydrodynamic and thermal boundary-layer parameters and then compare the newly found solutions to their counterparts obtained using Pohlhausen's cubic (KP3) and quartic (KP4) polynomials. As usual, the farfield around the cylinder is modeled using potential flow theory and the momentum-integral analysis is paired with Walz's empirical expression for the momentum thickness, which is based on a wide collection of experiments. This procedure permits retrieving explicit relations for the pressure-sensitive KP3, KP4, and MX4 velocity profiles across the boundary layer; one also obtains accurate approximations for the pressure distribution around the cylinder as well as an improved prediction of the separation point, namely, to within 0.87% of the actual location. In this process, refined estimates are produced for several characteristic parameters whose distributions are found to be in favorable agreement with experimental measurements and numerical simulations. These include the disturbance, momentum, and displacement thicknesses as well as the skin friction, pressure, and total drag coefficients. Finally, the thermal analysis is undertaken using both isothermal and isoflux boundary conditions. For each of these cases, closed-form analytical solutions are obtained for the local Nusselt number distribution around the cylinder, and these distributions are found to exhibit noticeably reduced errors relative to their classical values.
more »
« less
- Award ID(s):
- 1761675
- PAR ID:
- 10340984
- Date Published:
- Journal Name:
- Physics of Fluids
- Volume:
- 34
- Issue:
- 6
- ISSN:
- 1070-6631
- Page Range / eLocation ID:
- 063107
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
This paper focuses on the laminar boundary layer startup process (momentum and thermal) in incompressible flows. The unsteady boundary layer equations can be solved via similarity analysis by normalizing the stream-wise (x), wall-normal (y) and time (t) coordinates by a variable η and τ, respectively. The resulting ODEs are solved by a finite difference explicit algorithm. This can be done for two cases: flat plate flow where the change in pressure are zero (Blasius solution) and wedge or Falkner-Skan flow where the changes in pressure can be favorable (FPG) or adverse (APG). In addition, transient passive scalar transport is examined by setting several Prandtl numbers in the governing equation at two different wall thermal conditions: isothermal and isoflux. Numerical solutions for the transient evolution of the momentum and thermal boundary layer profiles are compared with analytical approximations for both small times (unsteady flow) and large (steady-state flow) times.more » « less
-
Shear sheltering is defined as the effect of the mean flow velocity profile in a boundary layer on the turbulence caused by an imposed gust. It has been studied extensively in applications involving boundary layer transition, where the primary concern is flow instabilities that are enhanced by turbulence in the flow outside the boundary layer. In aeroacoustic applications turbulent boundary layers interacting with blade trailing edges or roughness elements are an important source of sound, and the effect of shear sheltering on these noise sources has not been studied in detail. Since the surface pressure spectrum below the boundary layer is the primary driver of trailing edge and roughness noise, we will consider the effect that shear sheltering has on the surface pressure spectrum below a boundary layer. We will model the incoming turbulence as vortex sheets at specified heights above the surface and show, using classical boundary layer profiles and approximations to numerical results, how the mean flow velocity can be manipulated to alter the surface pressure spectrum and hence the radiated trailing edge noise.more » « less
-
Recent studies reveal that at large friction Reynolds number delta^+ the outer, inertially-dominated region of the turbulent boundary layer is composed of large scale zones of uniform momentum segregated by narrow fissures of concentrated vorticity. Experiments show that, when scaled by the boundary layer thickness, the fissure thickness is O(1/sqrt{delta^+}), while the dimensional jump in streamwise velocity across each fissure scales in proportion to the friction velocity u_tau. A simple model that exploits these essential elements of the turbulent boundary layer structure at large delta^+ is developed. First, a master wall-normal profile of streamwise velocity is constructed by placing a discrete number of fissures across the boundary layer. The number of fissures and their wall-normal locations follow scalings informed by analysis of the mean momentum equation. The fissures are then randomly displaced in the wall-normal direction, exchanging momentum as they move, to create an instantaneous velocity profile. This process is repeated to generate ensembles of streamwise velocity profiles from which statistical moments are computed. The modelled statistical moments are shown to agree remarkably well with those acquired from direct numerical simulations of turbulent channel flow at large delta^+. In particular, the model robustly reproduces the empirically observed sub-Gaussian behaviour for the skewness and kurtosis profiles over a large range of input parameters.more » « less
-
In this paper, we report a novel symmetry-breaking phenomenon that occurs in turbulent convective flow in periodic porous media in the intermediate porosity flow regime for values of porosities in between 0.8 and 0.9. Large eddy simulation is used to numerically simulate the momentum and thermal transport inside the porous medium at the microscale level. The phenomenon is observed to occur for periodically repeating porous media consisting of an in-line arrangement of circular cylinder solid obstacles, such as typically found in heat exchangers. The transition from symmetric to asymmetric flow occurs in between the pore scale Reynolds numbers of 37 (laminar) and 100 (turbulent), and asymmetric flow patterns are reported for Reynolds numbers up to 1,000. A Hopf bifurcation resulting in unsteady oscillatory laminar flow marks the origin of a secondary flow instability arising from the interaction of the shear layers around the solid obstacle. In turbulent flow, stochastic phase difference in the vortex wake oscillations caused by the secondary flow instability results in asymmetrical velocity and temperature distributions in the pore space. Consequently, high and low velocity flow channels are formed in the pore space that leads to the asymmetrical velocity and pressure distributions. At the macroscale level, symmetry-breaking results in residual transverse drag force components acting on the solid obstacle surfaces. The vortex wake oscillations caused by the secondary flow instability promote attached flow on the solid obstacle surface, which improves the surface averaged heat flux from the solid obstacle surface.more » « less
An official website of the United States government

