skip to main content


This content will become publicly available on September 1, 2024

Title: Sensitivity Analysis of Direct Numerical Simulation of a Spatially Developing Turbulent Mixing Layer to the Domain Dimensions
Abstract

Understanding spatial development of a turbulent mixing layer is essential for many engineering applications. However, the flow development is difficult to replicate in physical or numerical experiments. For this reason, the most attractive method for the mixing layer analysis is the direct numerical simulation (DNS), with the most control over the simulation inputs and free from modeling assumptions. On the other hand, the DNS cost often prevents conducting the sensitivity analysis of the simulation results to variations in the numerical procedure and thus, separating numerical and physical effects. In this paper, effects of the computational domain dimensions on statistics collected from DNS of a spatially developing incompressible turbulent mixing layer are analyzed with the focus on determining the domain dimensions suitable for studying the flow asymptotic state. In the simulations, the mixing layer develops between two coflowing laminar boundary layers formed on two sides of a sharp-ended splitter plate of a finite thickness with characteristics close to those of the untripped boundary layers in the experiments by Bell and Mehta, AIAA J., 28(12), 2034 (1990). The simulations were conducted using the spectral-element code Nek5000.

 
more » « less
Award ID(s):
1944568
NSF-PAR ID:
10483349
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
ASME
Date Published:
Journal Name:
Journal of Verification, Validation and Uncertainty Quantification
Volume:
8
Issue:
3
ISSN:
2377-2158
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Populated urban areas along many coastal regions are vulnerable to landfalling tropical cyclones (TCs). To the detriment of surface parameterizations in mesoscale models, the complexities of turbulence at high TC wind speeds in urban canopies are presently poorly understood. Thus, this study explores the impacts of urban morphology on TC-strength winds and boundary layer turbulence in landfalling TCs. To better quantify how urban structures interact with TC winds, large-eddy simulations (LESs) are conducted with the Cloud Model 1 (CM1). This implementation of CM1 includes immersed boundary conditions (IBCs) to represent buildings and eddy recycling to maintain realistic turbulent flow perturbations. Within the IBCs, an idealized coastal city with varying scales is introduced. TC winds impinge perpendicularly to the urbanized coastline. Numerical experiments show that buildings generate distinct, intricate flow patterns that vary significantly as the city structure is varied. Urban IBCs produce much stronger turbulent kinetic energy than is produced by conventional surface parameterizations. Strong effective eddy viscosity due to resolved eddy mixing is displayed in the wake of buildings within the urban canopy, while deep and enhanced effective eddy viscosity is present downstream. Such effects are not seen in a comparison LES using a simple surface parameterization with high roughness values. Wind tunneling effects in streamwise canyons enhance pedestrian-level winds well beyond what is possible without buildings. In the arena of regional mesoscale modeling, this type of LES framework with IBCs can be used to improve parameters in surface and boundary layer schemes to more accurately represent the drag coefficient and the eddy viscosity in landfalling TC boundary layers.

    Significance Statement

    This is among the first large-eddy simulation model studies to examine the impacts of tropical cyclone–like winds around explicitly resolved buildings. This work is a step forward in bridging the gap between engineering studies that use computational fluid dynamics models or laboratory experiments for flow through cities and mesoscale model simulations of landfalling tropical cyclones that use surface parameterizations specialized for urban land use.

     
    more » « less
  2. Abstract

    High fidelity near-wall velocity measurements in wall bounded fluid flows continue to pose a challenge and the resulting limitations on available experimental data cloud our understanding of the near-wall velocity behavior in turbulent boundary layers. One of the challenges is the spatial averaging and limited spatial resolution inherent to cross-correlation-based particle image velocimetry (PIV) methods. To circumvent this difficulty, we implement an explicit no-slip boundary condition in a wavelet-based optical flow velocimetry (wOFV) method. It is found that the no-slip boundary condition on the velocity field can be implemented in wOFV by transforming the constraint to the wavelet domain through a series of algebraic linear transformations, which are formulated in terms of the known wavelet filter matrices, and then satisfying the resulting constraint on the wavelet coefficients using constrained optimization for the optical flow functional minimization. The developed method is then used to study the classical problem of a turbulent channel flow using synthetic data from a direct numerical simulation (DNS) and experimental particle image data from a zero pressure gradient, high Reynolds number turbulent boundary layer. The results obtained by successfully implementing the no-slip boundary condition are compared to velocity measurements from wOFV without the no-slip condition and to a commercial PIV code, using the velocity from the DNS as ground truth. It is found that wOFV with the no-slip condition successfully resolves the near-wall profile with enhanced accuracy compared to the other velocimetry methods, as well as other derived quantities such as wall shear and turbulent intensity, without sacrificing accuracy away from the wall, leading to state of the art measurements in they+<1region of the turbulent boundary layer when applied to experimental particle images.

     
    more » « less
  3. The interaction of a turbulent, spatially developing crossflow with a transverse jet possesses several engineering and technological applications such as film cooling of turbine blades, exhaust plumes, thrust vector control, fuel injection, etc. Direct Numerical Simulation (DNS) of a jet in a crossflow under different streamwise pressure gradients (zero and favorable pressure gradient) is carried out. The purpose is to study the physics behind the transport phenomena and coherent structure dynamics in turbulent crossflow jets at different streamwise pressure gradients and low/high-velocity ratios (0.5 and 1). The temperature was regarded as a passive scalar with a molecular Prandtl number of 0.71. The analysis is performed by prescribing accurate turbulent information (instantaneous velocity and temperature) at the inlet of a computational domain. The upward motion of low-momentum fluid created by the “legs” of the counter-rotating vortex pair (CVP) encounters the downward inviscid flow coming from outside of the turbulent boundary layer, inducing a stagnation point and a shear layer. This layer is characterized by high levels of turbulent mixing, turbulence production, turbulent kinetic energy (TKE) and thermal fluctuations. The formation and development of the above-mentioned shear layer are more evident at higher velocity ratios. 
    more » « less
  4. Advisor: Dr. Guillermo Araya (Ed.)
    The present study provides fundamental knowledge on an issue in fluid dynamics that is not well understood: flow separation and its association with heat and contaminant transport. In the separated region, a swirling motion increases the fluid drag force on the object. Very often, this is undesirable because it can seriously reduce the performance of engineered devices such as aircraft and turbines. Furthermore, Computational Fluid Dynamics (CFD) has gained ground due to its relatively low cost, high accuracy, and versatility. The principal aim of this study is to numerically elucidate the details behind momentum and passive scalar transport phenomena during turbulent boundary layer separation resulting from a wall-curvature-driven pressure gradient. With Open- FOAM CFD software, the numerical discretization of Reynolds-Averaged Navier-Stokes and passive scalar transport equations will be described in two-dimensional domains via the assessment of two popular turbulence models (i.e., the Spalart-Allmaras and the K-w SST model). The computational domain reproduces a wind tunnel geometry from previously performed experiments by Baskaran et al. (JFM, vol. 182 and 232 “A turbulent flow over a curved hill.” Part 1 and Part 2). Only the velocity and pressure distribution were measured there, which will be used for validation purposes in the present study. A second aim in the present work is the scientific visualization of turbulent events and coherent structures via the ParaView toolkit and Unity game engine. Thus, fully immersive visualization approaches will be used via virtual reality (VR) and augmented reality (AR) technologies. A Virtual Wind Tunnel (VWT), developed for the VR approach, emulates the presence in a wind tunnel laboratory and has already employed fluid flow visualization from an existing numerical database with high temporal/spatial resolution, i.e., Direct Numeric Simulation (DNS). In terms of AR, a FlowVisXR app for smartphones and HoloLens has been developed for portability. It allows the user to see virtual 3D objects (i.e., turbulent coherent structures) invoked into the physical world using the device as the lens. 
    more » « less
  5. Abstract

    This study evaluates the methods of identifying the heightziof the top of the convective boundary layer (CBL) during winter (December and January) over the Great Lakes and nearby land areas using observations taken by the University of Wyoming King Air research aircraft during the Lake-Induced Convection Experiment (1997/98) and Ontario Winter Lake-effect Systems (2013/14) field campaigns. Since CBLs facilitate vertical mixing near the surface, the most direct measurement ofziis that above which the vertical velocity turbulent fluctuations are weak or absent. Thus, we usezifrom the turbulence method as the “reference value” to whichzifrom other methods, based on bulk Richardson number (Rib), liquid water content, and vertical gradients of potential temperature, relative humidity, and water vapor mixing ratio, are compared. The potential temperature gradient method using a threshold value of 0.015 K m−1for soundings over land and 0.011 K m−1for soundings over lake provided the estimates ofzithat are most consistent with the turbulence method. The Ribthreshold-based method, commonly used in numerical simulation studies, underestimatedzi. Analyzing the methods’ performance on the averaging windowzavgwe recommend usingzavg= 20 or 50 m forziestimations for lake-effect boundary layers. The present dataset consists of both cloudy and cloud-free boundary layers, some having decoupled boundary layers above the inversion top. Because cases of decoupled boundary layers appear to be formed by nearby synoptic storms, we recommend use of the more general term, elevated mixed layers.

    Significance Statement

    The depthziof the convective atmospheric boundary layer (CBL) strongly influences precipitation rates during lake-effect snowstorms (LES). However, variousziapproximation methods produce significantly different results. This study utilizes extensive concurrently collected observations by project aircraft during two LES field studies [Lake-Induced Convection Experiment (Lake-ICE) and OWLeS] to assess howzifrom common estimation methods compare with “reference”ziderived from turbulent fluctuations, a direct measure of CBL mixing. For soundings taken both over land and lake; with cloudy or cloud-free conditions, potential temperature gradient (PTG) methods provided the best agreement with the referencezi. A method commonly employed in numerical simulations performed relatively poorly. Interestingly, the PTG method worked equally well for “coupled” and elevated decoupled CBLs, commonly associated with nearby cyclones.

     
    more » « less