The turbulent channel flow database is produced from a direct numerical simulation (DNS) of wall bounded flow with periodic boundary conditions in the longitudinal and transverse directions, and no-slip conditions at the top and bottom walls. In the simulation, the Navier-Stokes equations are solved using a wall {normal, velocity {vorticity formulation. Solutions to the governing equations are provided using a Fourier-Galerkin pseudo-spectral method for the longitudinal and transverse directions and seventh-order Basis-splines (B-splines) collocation method in the wall normal direction. De-aliasing is performed using the 3/2-rule [3]. Temporal integration is performed using a low-storage, third-order Runge-Kutta method. Initially, the flow is driven using a constant volume flux control (imposing a bulk channel mean velocity of U = 1) until stationary conditions are reached. Then the control is changed to a constant applied mean pressure gradient forcing term equivalent to the shear stress resulting from the prior steps. Additional iterations are then performed to further achieve statistical stationarity before outputting fields.
more »
« less
Transitional Boundary Layer Data Set
Direct numerical simulation (DNS) of a transitional boundary layer over a plate with an elliptical leading edge. Navier-Stokes was discretized on a curvilinear grid and solved using a finite volume DNS code. A fractional-step algorithm was adopted, and the spatial discretization was a staggered volume-flux formulation. The viscous terms were integrated in time implicitly using the Crank-Nicolson and the advections terms were treated explicitly using the Adams-Bashforth. Pressure was treated using implicit Euler in the δp-form. The pressure equation was Fourier transformed in the span, and the resulting Helmholtz equation was solved for every spanwise wavenumber using two-dimensional multi-grid. After the simulation has reached a statistical stationary state, 4701 frames of data, which includes the 3 components of the velocity vector and the pressure, are generated and written in files that can be accessed directly by the database (FileDB system). Since the grid is staggered, data at the wall are not stored in the database. However, JHTDB provides values in the region between the wall and the first grid point, y∈[0, 0.0036], using 4th-order Lagrange polynomial inter- and extrapolation. The y-locations of the grid points in the vertical direction can be downloaded from this text file.
more »
« less
- Award ID(s):
- 2103874
- NSF-PAR ID:
- 10423316
- Publisher / Repository:
- Johns Hopkins Turbulence Databases
- Date Published:
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The data is from a direct numerical simulation (DNS) of homogeneous buoyancy driven turbulence on a 1024-cubed periodic grid. (See README-HBDT.pdf linked document for equations and details.) The simulation was performed with the variable-density version of the petascale CFDNS code. The database covers both the buoyancy driven increase in turbulence intensity as well as the buoyancy mediated turbulence decay.more » « less
-
more » « less
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 the
-
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
-
This study is part of the preliminary experimental investigations designed to assess the feasibility of using supercritical carbon dioxide (sCO2) in the vicinity of its critical point for thermal management applications. In the present study, sCO2 was used as a working fluid in a diffusion bonded 316/316L stainless steel test section having staggered micro pin fin array flow passages of hydraulic diameter 679 µm (0.679 mm). The test section was subjected to a single wall non-uniform heat flux boundary condition and was operated in a horizontal orientation. The primary objective was to characterize the heat transfer performance of sCO2 as it flows through the staggered pin fin array for experimental conditions that span its critical and pseudocritical point. Data analysis methods employing 2-D and 3-D heat transfer models of the test section were used to calculate the average heat transfer coefficients for a given set of experimental conditions. Experiments were conducted by varying the inlet temperature (18 ≤ T_in ≤ 50 °C) and for fixed mass flux (300 kg m-2 s-1), heat flux (40 W cm-2), and reduced pressure (1.1). Experimental data were also compared against the predictions of a correlation proposed for single phase flows in microchannel staggered diamond pin arrays. The correlation predicted the data within 4.3 % when the ratio, T_Bulk/T_PC exceeded 1. It was also found that the enhancement in the heat transfer, a result of employing staggered pin array flow geometry instead of microchannels, carries a commensurate penalty in pressure drop.more » « less
