skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Unsteady Incompressible Thermal Laminar Boundary Layers Subject to Streamwise Pressure Gradients
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
Award ID(s):
2314303 1847241
PAR ID:
10500008
Author(s) / Creator(s):
Publisher / Repository:
AIP
Date Published:
Journal Name:
Proceedings of the 21st International Conference of Numerical Analysis and Applied Mathematics
Format(s):
Medium: X
Location:
Crete, Greece
Sponsoring Org:
National Science Foundation
More Like this
  1. A transient stability flow analysis is performed using the unsteady laminar boundary layer equations. The flow dynamics are studied via the Navier–Stokes equations. In the case of external spatially developing flow, the differential equations are reduced via Prandtl or boundary-layer assumptions, consisting of continuity and momentum conservation equations. Prescription of streamwise pressure gradients (decelerating and accelerating flows) is carried out by an impulsively started Falkner–Skan (FS) or wedge-flow similarity flow solution in the case of flat plate or a Blasius solution for particular zero-pressure gradient case. The obtained mean streamwise velocity and its derivatives from FS flows are then inserted into the well-known Orr–Sommerfeld equation of small disturbances at different dimensionless times (τ). Finally, the corresponding eigenvalues are dynamically computed for temporal stability analysis. A finite difference algorithm is effectively applied to solve the Orr–Sommerfeld equations. It is observed that flow acceleration or favorable pressure gradients (FPGs) lead to a significantly shorter transient period before reaching steady-state conditions, as the developed shear layer is notably thinner compared to cases with adverse pressure gradients (APGs). During the transient phase (i.e., for τ<1), the majority of the flow modifications are confined to the innermost 20–25% of the boundary layer, in proximity to the wall. In the context of temporal flow stability, the magnitude of the pressure gradient is pivotal in determining the streamwise extent of the Tollmien–Schlichting (TS) waves. In highly accelerated laminar flows, these waves experience considerable elongation. Conversely, under the influence of a strong adverse pressure gradient, the characteristic streamwise length of the smallest unstable wavelength, which is necessary for destabilization via TS waves, is significantly reduced. Furthermore, flows subjected to acceleration (β > 0) exhibit a higher propensity to transition towards a more stable state during the initial transient phase. For instance, the time response required to reach the steady-state critical Reynolds number was approximately 1τ for β = 0.18 (FPG) and τ = 6.8 for β = −0.18 (APG). 
    more » « less
  2. null (Ed.)
    An analytical approach has been developed to model the rapid term contribution to the unsteady surface pressure fluctuations in wall jet turbulent boundary layer flows. The formulation is based on solving Poisson’s equation for the turbulent wall pressure by integrating the source terms (Kraichnan, 1956). The inputs for the model are obtained from 2D time-resolved Particle Image Velocimetry measurements performed in a wall jet flow. The wall normal turbulence wavenumber two-point cross-spectra is determined using an extension of the von Kármán homogeneous turbulence spectrum. The model is applied to compare and understand the baseline flow in the wall jet and to study the attenuation in surface pressure fluctuations by unidirectional canopies (Gonzales et al, 2019). Different lengthscale formulations are tested and we observe that the wall jet flow boundary layer contributes to the surface pressure fluctuations from two distinct regions. The high frequency spectrum is captured well. However, the low frequency range of the spectrum is not entirely captured. This is because we have used PIV data only up to a height of 2.3𝜹, whereas the largest turbulent lengthscales in the wall jet are on the order of 𝒚𝟏/𝟐≈𝟔𝜹. Using the flow data obtained from PIV and Pitot probe measurements, the model predicts a reduction in the surface pressure due to canopy at low frequencies. 
    more » « less
  3. Abstract Ocean wave-powered reverse osmosis (RO) desalination is an emerging field of study that combines the utilization of ocean energy and RO desalination techniques. However, due to the significant fluctuations in pressure and flow rate within the hydraulic system, an accurate transient model is necessary to estimate its performance accurately and effectively. This paper presents a two-dimensional transient model based on the pressure-correction algorithm to simulate the channel flow with porous walls and time-dependent inlet conditions. The coupled pressure, velocity, and salt concentration problem is solved iteratively by decoupling each term and updating them separately. The model is validated by comparing the results with analytical film theory which estimates the formation of the concentration polarization layer under constant inlet conditions. The performance of the RO systems, especially the concentration polarization phenomenon at the member surface, is investigated using different input conditions, including constant flow condition and sinusoidal flow condition. The salt concentration and permeate flux at the membrane boundary are studied to understand the effect of the dynamic inputs. Results show that the system can reach a higher maximum wall concentration and higher average recovery ratio in sinusoidal signal compared with the constant input. The model’s adaptability to different flow regimes, from steady to sinusoidal, underscores its potential as a valuable tool in optimizing RO desalination powered by ocean wave energy. 
    more » « less
  4. An incoming canonical spatially developing turbulent boundary layer (SDTBL) over a 2-D curved hill is numerically investigated via the Reynolds-averaged Navier–Stokes (RANS) equations plus two eddy-viscosity models: the K−ω SST (henceforth SST) and the Spalart–Allmaras (henceforth SA) turbulence models. A spatially evolving thermal boundary layer has also been included, assuming temperature as a passive scalar (Pr = 0.71) and a turbulent Prandtl number, Prt, of 0.90 for wall-normal turbulent heat flux modeling. The complex flow with a combined strong adverse/favorable streamline curvature-driven pressure gradient caused by concave/convex surface curvatures has been replicated from wind-tunnel experiments from the literature, and the measured velocity and pressure fields have been used for validation purposes (the thermal field was not experimentally measured). Furthermore, direct numerical simulation (DNS) databases from the literature were also employed for the incoming turbulent flow assessment. Concerning first-order statistics, the SA model demonstrated a better agreement with experiments where the turbulent boundary layer remained attached, for instance, in Cp, Cf, and Us predictions. Conversely, the SST model has shown a slightly better match with experiments over the flow separation zone (in terms of Cp and Cf) and in Us profiles just upstream of the bubble. The Reynolds analogy, based on the St/(Cf/2) ratio, holds in zero-pressure gradient (ZPG) zones; however, it is significantly deteriorated by the presence of streamline curvature-driven pressure gradient, particularly due to concave wall curvature or adverse-pressure gradient (APG). In terms of second-order statistics, the SST model has better captured the positively correlated characteristics of u′ and v′ or positive Reynolds shear stresses ( > 0) inside the recirculating zone. Very strong APG induced outer secondary peaks in and turbulence production as well as an evident negative slope on the constant shear layer. 
    more » « less
  5. ‘Ground effect’ refers to the enhanced performance enjoyed by fliers or swimmers operating close to the ground. We derive a number of exact solutions for this phenomenon, thereby elucidating the underlying physical mechanisms involved in ground effect. Unlike previous analytic studies, our solutions are not restricted to particular parameter regimes, such as ‘weak’ or ‘extreme’ ground effect, and do not even require thin aerofoil theory. Moreover, the solutions are valid for a hitherto intractable range of flow phenomena, including point vortices, uniform and straining flows, unsteady motions of the wing, and the Kutta condition. We model the ground effect as the potential flow past a wing inclined above a flat wall. The solution of the model requires two steps: firstly, a coordinate transformation between the physical domain and a concentric annulus; and secondly, the solution of the potential flow problem inside the annulus. We show that both steps can be solved by introducing a new special function which is straightforward to compute. Moreover, the ensuing solutions are simple to express and offer new insight into the mathematical structure of ground effect. In order to identify the missing physics in our potential flow model, we compare our solutions against new experimental data. The experiments show that boundary layer separation on the wing and wall occurs at small angles of attack, and we suggest ways in which our model could be extended to account for these effects. 
    more » « less