An official website of the United States government
Motivated by the need for accurate determination of wall shear stress from profile measurements in turbulent boundary layer flows, the total shear stress balance is analysed and reformulated using several well-established semi-empirical relations. The analysis highlights the significant effect that small pressure gradients can have on parameters deduced from data even in nominally zero pressure gradient boundary layers. Using the comprehensive shear stress balance together with the log-law equation, it is shown that friction velocity, roughness length and zero-plane displacement can be determined with only velocity and turbulent shear stress profile measurements at a single streamwise location for nominally zero pressure gradient turbulent boundary layers. Application of the proposed analysis to turbulent smooth- and rough-wall experimental data shows that the friction velocity is determined with accuracy comparable to force balances (approximately 1 %–4 %). Additionally, application to boundary layer data from previous studies provides clear evidence that the often cited discrepancy between directly measured friction velocities (e.g. using force balances) and those derived from traditional total shear stress methods is likely due to the small favourable pressure gradient imposed by a fixed cross-section facility. The proposed comprehensive shear stress analysis can account for these small pressure gradients and allows more accurate boundary layer wall shear stress or friction velocity determination using commonly available mean velocity and shear stress profile data from a single streamwise location.
more »
« less
Volumetric lattice Boltzmann method for wall stresses of image-based pulsatile flows
Abstract Image-based computational fluid dynamics (CFD) has become a new capability for determining wall stresses of pulsatile flows. However, a computational platform that directly connects image information to pulsatile wall stresses is lacking. Prevailing methods rely on manual crafting of a hodgepodge of multidisciplinary software packages, which is usually laborious and error-prone. We present a new computational platform, to compute wall stresses in image-based pulsatile flows using the volumetric lattice Boltzmann method (VLBM). The novelty includes: (1) a unique image processing to extract flow domain and local wall normality, (2) a seamless connection between image extraction and VLBM, (3) an en-route calculation of strain-rate tensor, and (4) GPU acceleration (not included here). We first generalize the streaming operation in the VLBM and then conduct application studies to demonstrate its reliability and applicability. A benchmark study is for laminar and turbulent pulsatile flows in an image-based pipe (Reynolds number: 10 to 5000). The computed pulsatile velocity and shear stress are in good agreements with Womersley's analytical solutions for laminar pulsatile flows and concurrent laboratory measurements for turbulent pulsatile flows. An application study is to quantify the pulsatile hemodynamics in image-based human vertebral and carotid arteries including velocity vector, pressure, and wall-shear stress. The computed velocity vector fields are in reasonably well agreement with MRA (magnetic resonance angiography) measured ones. This computational platform is good for image-based CFD with medical applications and pore-scale porous media flows in various natural and engineering systems.
more »
« less
- Award ID(s):
- 1803845
- PAR ID:
- 10381938
- Date Published:
- Journal Name:
- Scientific Reports
- Volume:
- 12
- Issue:
- 1
- ISSN:
- 2045-2322
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The transition from laminar to turbulent flow is of great interest since it is one of the most difficult and unsolved problems in fluids engineering. The transition processes are significantly important because the transition has a huge impact on almost all systems that come in contact with a fluid flow by altering the mixing, transport, and drag properties of fluids even in simple pipe and channel flows. Generally, in most transportation systems, the transition to turbulence causes a significant increase in drag force, energy consumption, and, therefore, operating cost. Thus, understanding the underlying mechanisms of the laminar-to-turbulent transition can be a major benefit in many ways, especially economically. There have been substantial previous studies that focused on testing the stability of laminar flow and finding the critical amplitudes of disturbances necessary to trigger the transition in various wall-bounded systems, including circular pipes and square ducts. However, there is still no fundamental theory of transition to predict the onset of turbulence. In this study, we perform direct numerical simulations (DNS) of the transition flows from laminar to turbulence in a channel flow. Specifically, the effects of different magnitudes of perturbations on the onset of turbulence are investigated. The perturbation magnitudes vary from 0.001 (0.1%) to 0.05 (5%) of a typical turbulent velocity field, and the Reynolds number is from 5,000 to 40,000. Most importantly, the transition behavior in this study was found to be in good agreement with other reported studies performed for fluid flow in pipes and ducts. With the DNS results, a finite amplitude stability curve was obtained. The critical magnitude of perturbation required to cause transition was observed to be inversely proportional to the Reynolds number for the magnitude from 0.01 to 0.05. We also investigated the temporal behavior of the transition process, and it was found that the transition time or the time required to begin the transition process is inversely correlated with the Reynolds number only for the magnitude from 0.02 to 0.05, while different temporal behavior occurs for smaller perturbation magnitudes. In addition to the transition time, the transition dynamics were investigated by observing the time series of wall shear stress. At the onset of transition, the shear stress experiences an overshoot, then decreases toward sustained turbulence. As expected, the average values of the wall shear stress in turbulent flow increase with the Reynolds number. The change in the wall shear stress from laminar to overshoot was, of course, found to increase with the Reynolds number. More interestingly was the observed change in wall shear stress from the overshoot to turbulence. The change in magnitude appears to be almost insensitive to the Reynolds number and the perturbation magnitude. Because the change in wall shear stress is directly proportional to the pumping power, these observations could be extremely useful when determining the required pumping power in certain flow conditions. Furthermore, the stability curve and wall shear stress changes can be considered robust features for future applications, and ultimately interpreted as evidence of progress toward solving the unresolved fluids engineering problem.more » « less
-
Inlet and outlet boundary conditions (BCs) play an important role in newly emerged image-based computational hemodynamics for blood flows in human arteries anatomically extracted from medical images. We developed physiological inlet and outlet BCs based on patients’ medical data and integrated them into the volumetric lattice Boltzmann method. The inlet BC is a pulsatile paraboloidal velocity profile, which fits the real arterial shape, constructed from the Doppler velocity waveform. The BC of each outlet is a pulsatile pressure calculated from the three-element Windkessel model, in which three physiological parameters are tuned by the corresponding Doppler velocity waveform. Both velocity and pressure BCs are introduced into the lattice Boltzmann equations through Guo’s non-equilibrium extrapolation scheme. Meanwhile, we performed uncertainty quantification for the impact of uncertainties on the computation results. An application study was conducted for six human aortorenal arterial systems. The computed pressure waveforms have good agreement with the medical measurement data. A systematic uncertainty quantification analysis demonstrates the reliability of the computed pressure with associated uncertainties in the Windkessel model. With the developed physiological BCs, the image-based computation hemodynamics is expected to provide a computation potential for the noninvasive evaluation of hemodynamic abnormalities in diseased human vessels.more » « less
-
null (Ed.)Abstract Coherent structures are critical for controlling turbulent boundary layers due to their roles in momentum and heat transfer in the flow. Turbulent coherent structures can be detected by measuring wall shear stresses that are footprints of coherent structures. In this study, wall shear stress fluctuations were measured simultaneously in a zero pressure gradient turbulent boundary layer using two house-made wall shear stress probes aligned in the spanwise direction. The wall shear stress probe consisted of two hot-wires on the wall aligned in a V-shaped configuration for measuring streamwise and spanwise shear stresses, and their performance was validated in comparison with a direct numerical simulation result. Relationships between measured wall shear stress fluctuations and streamwise velocity fluctuations were analyzed using conditional sampling techniques. The peak detection method and the variable-interval time-averaging (VITA) method showed that quasi-streamwise vortices were inclined toward the streamwise direction. When events were simultaneously detected by the two probes, stronger fluctuations in streamwise velocity were detected, which suggests that stronger coherent structures were detected. In contrast to the former two methods, the hibernating event detection method detects events with lower wall shear stress fluctuations. The ensemble-averaged mean velocity profile of hibernating events was shifted upward compared to the law of the wall, which suggests low drag status of the coherent structures related with hibernating events. These methods suggest significant correlations between wall shear stress fluctuations and coherent structures, which could motivate flow control strategies to fully exploit these correlations.more » « less
-
Abstract. In this paper, a three-dimensional two-phase flow solver, SedFoam-2.0, is presented for sediment transport applications. The solver is extended from twoPhaseEulerFoam available in the 2.1.0 release of the open-source CFD (computational fluid dynamics) toolbox OpenFOAM. In this approach the sediment phase is modeled as a continuum, and constitutive laws have to be prescribed for the sediment stresses. In the proposed solver, two different intergranular stress models are implemented: the kinetic theory of granular flows and the dense granular flow rheology μ(I). For the fluid stress, laminar or turbulent flow regimes can be simulated and three different turbulence models are available for sediment transport: a simple mixing length model (one-dimensional configuration only), a k − ε, and a k − ω model. The numerical implementation is demonstrated on four test cases: sedimentation of suspended particles, laminar bed load, sheet flow, and scour at an apron. These test cases illustrate the capabilities of SedFoam-2.0 to deal with complex turbulent sediment transport problems with different combinations of intergranular stress and turbulence models.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1038/s41598-022-05269-w
