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 region of the turbulent boundary layer when applied to experimental particle images.
more »
« less
Magnus Force Estimation Using Gauss’s Principle of Least Constraint
The flow around a rotating cylinder is one of the fundamental problems that have piqued the interests of many venerable fluid mechanicians since the time of Rayleigh. The force caused by the rotation of the cylinder has always been considered as an immediate consequence of viscosity, since the potential flow model failed entirely to predict the value of the circulation due to the lack of a Kutta-like condition. On the other hand, Glauert modeled the flow outside the boundary layer of a rotating cylinder as a potential flow with an unknown circulation. He then obtained an approximate solution of Prandtl’s boundary-layer equations and applied the no-slip condition to estimate the circulation in the outer flow. Interestingly, for rapidly rotating cylinders ([Formula: see text]), up to fourth-order in the small parameter [Formula: see text], the obtained circulation is independent of viscosity. In this work, we use Glauert’s model of the outer flow (i.e., a potential flow with an unknown circulation). However, instead of the tedious boundary-layer calculations, we rely on Gauss’s principle of least constraint to obtain the unknown circulation. A perfect match with Glauert’s solution is found. Moreover, our solution, in contrast to Glauert’s, points to the existence of different physics at small rotational speeds. The obtained results, given their perfect matching with Glauert’s solution (relying on the no-slip condition), point to a potential equivalence between the no-slip condition and fluid body forces.
more »
« less
- Award ID(s):
- 2332556
- PAR ID:
- 10539099
- Publisher / Repository:
- AIAA
- Date Published:
- Journal Name:
- AIAA Journal
- Volume:
- 62
- Issue:
- 5
- ISSN:
- 0001-1452
- Page Range / eLocation ID:
- 1962 to 1969
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We obtain regularity results in weighted Sobolev spaces for the solution of the obstacle problem for the integral fractional Laplacian [Formula: see text] in a Lipschitz bounded domain [Formula: see text] satisfying the exterior ball condition. The weight is a power of the distance to the boundary [Formula: see text] of [Formula: see text] that accounts for the singular boundary behavior of the solution for any [Formula: see text]. These bounds then serve us as a guide in the design and analysis of a finite element scheme over graded meshes for any dimension [Formula: see text], which is optimal for [Formula: see text].more » « less
-
null (Ed.)Abstract ‘Horizontal convection’ (HC) is the generic name for the flow resulting from a buoyancy variation imposed along a horizontal boundary of a fluid. We study the effects of rotation on three-dimensional HC numerically in two stages: first, when baroclinic instability is suppressed and, second, when it ensues and baroclinic eddies are formed. We concentrate on changes to the thickness of the near-surface boundary layer, the stratification at depth, the overturning circulation and the flow energetics during each of these stages. Our results show that, for moderate flux Rayleigh numbers ( $$O(1{0}^{11} )$$ ), rapid rotation greatly alters the steady-state solution of HC. When the flow is constrained to be uniform in the transverse direction, rapidly rotating solutions do not support a boundary layer, exhibit weaker overturning circulation and greater stratification at all depths. In this case, diffusion is the dominant mechanism for lateral buoyancy flux and the consequent buildup of available potential energy leads to baroclinically unstable solutions. When these rapidly rotating flows are perturbed, baroclinic instability develops and baroclinic eddies dominate both the lateral and vertical buoyancy fluxes. The resulting statistically steady solution supports a boundary layer, larger values of deep stratification and multiple overturning cells compared with non-rotating HC. A transformed Eulerian-mean approach shows that the residual circulation is dominated by the quasi-geostrophic eddy streamfunction and that the eddy buoyancy flux has a non-negligible interior diabatic component. The kinetic and available potential energies are greater than in the non-rotating case and the mixing efficiency drops from $${\sim }0. 7$$ to $${\sim }0. 17$$ . The eddies play an important role in the formation of the thermal boundary layer and, together with the negatively buoyant plume, help establish deep stratification. These baroclinically active solutions have characteristics of geostrophic turbulence.more » « less
-
null (Ed.)The [Formula: see text]-pinch is a classical steady state for the MHD model, where a confined plasma fluid is separated by vacuum, in the presence of a magnetic field which is generated by a prescribed current along the [Formula: see text] direction. We develop a variational framework to study its stability in the absence of viscosity effect, and demonstrate for the first time that such a [Formula: see text]-pinch is always unstable. Moreover, we discover a sufficient condition such that the eigenvalues can be unbounded, which leads to ill-posedness of the linearized MHD system.more » « less
-
Taylor–Couette (TC) flow, the flow between two independently rotating and co-axial cylinders, is commonly used as a canonical model for shear flows. Unlike plane Couette flow, pinned secondary flows can be found in TC flow. These are known as Taylor rolls and drastically affect the flow behaviour. We study the possibility of modifying these secondary structures using patterns of stress-free and no-slip boundary conditions on the inner cylinder. For this, we perform direct numerical simulations of narrow-gap TC flow with pure inner-cylinder rotation at four different shear Reynolds numbers up to $$Re_s=3\times 10^4$$ . We find that one-dimensional azimuthal patterns do not have a significant effect on the flow topology, and that the resulting torque is a large fraction ( $$\sim$$ 80 %–90 %) of torque in the fully no-slip case. One-dimensional axial patterns decrease the torque more, and for certain pattern frequency disrupt the rolls by interfering with the existing Reynolds stresses that generate secondary structures. For $$Re\geq 10^4$$ , this disruption leads to a smaller torque than what would be expected from simple boundary layer effects and the resulting effective slip length and slip velocity. We find that two-dimensional checkerboard patterns have similar behaviour to azimuthal patterns and do not affect the flow or the torque substantially, but two-dimensional spiral inhomogeneities can move around the pinned secondary flows as they induce persistent axial velocities. We quantify the roll's movement for various angles and the widths of the spiral pattern, and find a non-monotonic behaviour as a function of pattern angle and pattern frequency.more » « less
An official website of the United States government

