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
Thin liquid films in a funnel
We explore flow of a completely wetting fluid in a funnel, with particular focus on contact line instabilities at the fluid front. While the flow in a funnel may be related to a number of other flow configurations as limiting cases, understanding its stability is complicated due to the presence of additional azimuthal curvature, as well as due to convergent flow effects imposed by the geometry. The convergent nature of the flow leads to thickening of the film, therefore influencing its stability properties. In this work, we analyse these stability properties by combining physical experiments, asymptotic modelling, self-similar type of analysis and numerical simulations. We show that an appropriate long-wave-based model, supported by the input from experiments, simulations and linear stability analysis that originates from the flow down an incline plane, provides a basic insight allowing an understanding of the development of contact line instability and emerging length scales.
more »
« less
- Award ID(s):
- 1815613
- NSF-PAR ID:
- 10287148
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 924
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
In nanoporous rocks, potential size/mobility exclusion and fluid-rock interactions in nano-sized pores and pore throats may turn the rock into a semi-permeable membrane, blocking or hindering the passage of certain molecules while allowing other molecules to pass freely. In this work, we conducted several experiments to investigate whether Niobrara samples possess such sieving properties on hydrocarbon molecules. Molecular dynamics simulation of hydrocarbon adsorption was performed to help understand the trends observed in the experiments. The procedure of the experiments includes pumping of liquid binary hydrocarbon mixtures (C10, C17) of known compositions into Niobrara samples, collecting of the effluents from the samples, and analysis of the compositions of the effluents. A specialized experimental setup that uses an in-line filter as a mini-core holder was built for this investigation. Niobrara samples were cored and machined into 0.5-inch diameter and 0.7-inch length mini-cores. Hydrocarbon mixtures were injected into the mini-cores and effluents were collected periodically and analyzed using gas chromatography. To understand the potential effects of hydrocarbon-rock interactions on their transport, molecular dynamics simulations were performed to clarify the adsorption of C10 and C17 molecules on calcite surfaces using all-atom models. Experimental results show that the heavier component (C17) in the injected fluid was noticeably hindered. After the start of the experiment, the fraction of the lighter component (C10) in the produced fluid gradually increased and eventually reached levels that fluctuated within a range above the fraction of C10 in the original fluid; besides, the fraction of C17 increased in the fluid upstream of the sample. Both observations indicate the presence of membrane properties of the sample to this hydrocarbon mixture. Simulation results suggest that, for a calcite surface in equilibrium with a binary mixture of C10 and C17, more C17 molecules adsorb on the carbonate surface than the C10 molecules, providing a mechanism that directly supports the experimental observations. Some experimental observations suggest that size/mobility exclusion should also exist. This experimental study is the first evidence that nanoporous reservoir rocks may possess membrane properties that can filter hydrocarbon molecules. Component separation due to membrane properties has not been considered in any reservoir simulation models. The consequence of this effect and its dependence on the mixture and environmental conditions (surface, pressure, temperature) are worthy of further investigations.more » « less
-
Partially wetting droplets under an airflow can exhibit complex behaviours that arise from the coupling of surface tension, inertia of the external flow and contact-line dynamics. Recent experiments by Hooshanginejad et al. ( J. Fluid Mech. , vol. 901, 2020) revealed that a millimetric partially wetting water droplet under an impinging jet can oscillate in place, split or depin away from the jet, depending on the magnitude (i.e. $5\unicode{x2013}20\ {\rm m}\ {\rm s}^{-1}$ ) and position of the jet. To rationalise the experimental observations, we develop a two-dimensional lubrication model of the droplet that incorporates the external pressure of the impinging high-Reynolds-number jet, in addition to the capillary and hydrostatic pressures of the droplet. Distinct from the previous model by Hooshanginejad et al. ( J. Fluid Mech. , vol. 901, 2020), we simulate the motion of the contact line using precursor film and disjoining pressure, which allows us to capture a wider range of droplet behaviours, including the droplet dislodging to one side. Our simulations exhibit a comparable time-scale of droplet deformations and similar outcomes as the experimental observations. We also obtain the analytical steady-state solutions of the droplet shapes and construct the minimum criteria for splitting and depinning.more » « less
-
Legged robots require knowledge of pose and velocity in order to maintain stability and execute walking paths. Current solutions either rely on vision data, which is susceptible to environmental and lighting conditions, or fusion of kinematic and contact data with measurements from an inertial measurement unit (IMU). In this work, we develop a contact-aided invariant extended Kalman filter (InEKF) using the theory of Lie groups and invariant observer design. This filter combines contact-inertial dynamics with forward kinematic corrections to estimate pose and velocity along with all current contact points. We show that the error dynamics follows a log-linear autonomous differential equation with several important consequences: (a) the observable state variables can be rendered convergent with a domain of attraction that is independent of the system’s trajectory; (b) unlike the standard EKF, neither the linearized error dynamics nor the linearized observation model depend on the current state estimate, which (c) leads to improved convergence properties and (d) a local observability matrix that is consistent with the underlying nonlinear system. Furthermore, we demonstrate how to include IMU biases, add/remove contacts, and formulate both world-centric and robo-centric versions. We compare the convergence of the proposed InEKF with the commonly used quaternion-based extended Kalman filter (EKF) through both simulations and experiments on a Cassie-series bipedal robot. Filter accuracy is analyzed using motion capture, while a LiDAR mapping experiment provides a practical use case. Overall, the developed contact-aided InEKF provides better performance in comparison with the quaternion-based EKF as a result of exploiting symmetries present in system.more » « less
-
Experiments were carried out to observe the flow inside counterflow atomizers over a range of operating conditions and fluid properties. Liquids used were water and propylene glycol, while the gas was either air or helium. Liquid flow rates ranged from 10 ml/min to 40 ml/min, with gas liquid ratio (GLR) ranging from 0.1 to 0.6. The primary experiments used the 7-BM line of the Advanced Photon Source in Argonne National Laboratories with a 2.6 mm atomizer produced from (Poly)Ethyl-Ether-Ketone (PEEK). The X-Ray beam was operated in phase contrast mode, leading to interference patterns near the gas-liquid interface and enabling a qualitative understanding of the flow structure. Complementary optical work applied laser shadowgraphy to a 1 mm orifice atomizer constructed with quartz capillary tubing. A diffuse pulsed Nd:YAG laser backlight captured instantaneous gas-liquid interface positions in the internal flow. With both techniques, two distinct flow behaviors are observed corresponding to low and high GLR values. At low GLR, the inertia of the injected gas is insufficient to penetrate the liquid downflow. The gas stream entering the mixing chamber in the upstream direction is immediately deflected by the denser liquid and enters the discharge tube around a central liquid jet, which is sheared and accelerated by the surrounding gas, leading to breakup. A distinct frequency of jet breakup is observed inside the discharge tube, with the liquid jet oscillating and fragmenting against the walls. The situation at high GLR is quite different, however, as the incoming gas stream asymmetrically penetrates upstream into the mixing chamber, taking the form of a high-speed jet confined along one wall, and displaying a flapping instability as it encounters the liquid flowing downstream. This flapping causes violent mixing, resulting in a highly disturbed interface, along with the generation of liquid ligaments and gas bubbles. This two-phase mixture enters the discharge tube with no liquid jet formation evident for this case. The transition between these two regimes is explored by changing the liquid viscosity and gas molar mass, and weak sensitivity to fluid properties is observed. Further, quantitative image analysis techniques applied to the low and high GLR cases allow extraction of the frequencies of the liquid jet in the discharge tube at low GLR, as well as the flapping mode at high GLR.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2021.626
