skip to main content


Title: Reduced modelling and global instability of finite-Reynolds-number flow in compliant rectangular channels
Experiments have shown that flow in compliant microchannels can become unstable at a much lower Reynolds number than the corresponding flow in a rigid conduit. Therefore, it has been suggested that the wall's elastic compliance can be exploited towards new modalities of microscale mixing. While previous studies mainly focused on the local instability induced by the fluid–structure interactions (FSIs) in the system, we derive a one-dimensional (1-D) model to study the FSI's effect on the global instability. The proposed 1-D FSI model is tailored to long, shallow rectangular microchannels with a deformable top wall, similar to the experiments. Going beyond the usual lubrication flows analysed in these geometries, we include finite fluid inertia and couple the reduced flow equations to a novel reduced 1-D wall deformation equation. Although a quantitative comparison with previous experiments is difficult, the behaviours of the proposed model show, qualitatively, agreement with the experimental observations, and capture several key effects. Specifically, we find the critical conditions under which the inflated base state of the 1-D FSI model is linearly unstable to infinitesimal perturbations. The critical Reynolds numbers predicted are in agreement with experimental observations. The unstable modes are highly oscillatory, with frequencies close to the natural frequency of the wall, suggesting that the observed instabilities are resonance phenomena. Furthermore, during the start-up from an undeformed initial state, self-sustained oscillations can be triggered by FSI. Our modelling framework can be applied to other microfluidic systems with similar geometric scale separation under different operating conditions.  more » « less
Award ID(s):
1705637 2029540
NSF-PAR ID:
10376754
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
950
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. This study presents the first 3D two-way coupled fluid structure interaction (FSI) simulation of a hybrid anechoic wind tunnel (HAWT) test section with modeling all important effects, such as turbulence, Kevlar wall porosity and deflection, and reveals for the first time the complete 3D flow structure associated with a lifting model placed into a HAWT. The Kevlar deflections are captured using finite element analysis (FEA) with shell elements operated under a membrane condition. Three-dimensional RANS CFD simulations are used to resolve the flow field. Aerodynamic experimental results are available and are compared against the FSI results. Quantitatively, the pressure coefficients on the airfoil are in good agreement with experimental results. The lift coefficient was slightly underpredicted while the drag was overpredicted by the CFD simulations. The flow structure downstream of the airfoil showed good agreement with the experiments, particularly over the wind tunnel walls where the Kevlar windows interact with the flow field. A discrepancy between previous experimental observations and juncture flow-induced vortices at the ends of the airfoil is found to stem from the limited ability of turbulence models. The qualitative behavior of the flow, including airfoil pressures and cross-sectional flow structure is well captured in the CFD. From the structural side, the behavior of the Kevlar windows and the flow developing over them is closely related to the aerodynamic pressure field induced by the airfoil. The Kevlar displacement and the transpiration velocity across the material is dominated by flow blockage effects, generated aerodynamic lift, and the wake of the airfoil. The airfoil wake increases the Kevlar window displacement, which was previously not resolved by two-dimensional panel-method simulations. The static pressure distribution over the Kevlar windows is symmetrical about the tunnel mid-height, confirming a dominantly two-dimensional flow field. 
    more » « less
  3. 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
  4. 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
  5. We develop a theory of fluid--structure interaction (FSI) between an oscillatory Newtonian fluid flow and a compliant conduit. We consider the canonical geometries of a 2D channel with a deformable top wall and an axisymmetric deformable tube. Focusing on the hydrodynamics, we employ a linear relationship between wall displacement and hydrodynamic pressure, which has been shown to be suitable for a leading-order-in-slenderness theory. The slenderness assumption also allows the use of lubrication theory, and the flow rate is related to the pressure gradient (and the tube/wall deformation) via the classical solutions for oscillatory flow in a channel and in a tube (attributed to Womersley). Then, by two-way coupling the oscillatory flow and the wall deformation via the continuity equation, a one-dimensional nonlinear partial differential equation (PDE) governing the instantaneous pressure distribution along the conduit is obtained, without \textit{a priori} assumptions on the magnitude of the oscillation frequency (\textit{i.e.}, at arbitrary Womersley number). We find that the cycle-averaged pressure (for harmonic pressure-controlled conditions) deviates from the expected steady pressure distribution, suggesting the presence of a streaming flow. An analytical perturbative solution for a weakly deformable conduit is obtained to rationalize how FSI induces such streaming. In the case of a compliant tube, the results obtained from the proposed reduced-order PDE and its perturbative solutions are validated against three-dimensional, two-way-coupled direct numerical simulations. We find good agreement between theory and simulations for a range of dimensionless parameters characterizing the oscillatory flow and the FSI, demonstrating the validity of the proposed theory of oscillatory flows in compliant conduits at arbitrary Womersley number. 
    more » « less