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
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
- PAR ID:
- 10376754
- 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
-
-
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
-
This experimental study explores the physical mechanisms by which a transverse jet’s upstream shear layer can transition from being a convective instability to an absolute/global instability as the jet-to-cross-flow momentum flux ratio $$J$$ is reduced. As first proposed in computational studies by Iyer & Mahesh ( J. Fluid Mech. , vol. 790, 2016, pp. 275–307), the upstream shear layer just beyond the jet injection may be analogous to a local counter-current shear layer, which is known for a planar geometry to become absolutely unstable at a large enough counter-current shear layer velocity ratio, $$R_{1}$$ . The present study explores this analogy for a range of transverse jet momentum flux ratios and jet-to-cross-flow density ratios $$S$$ , for jets containing differing species concentrations (nitrogen, helium and acetone vapour) at several different jet Reynolds numbers. These studies make use of experimental data extracted from stereo particle image velocimetry as well as simultaneous stereo particle image velocimetry and acetone planar laser-induced fluorescence imaging. They provide experimental evidence for the relevance of the counter-current shear layer analogy to upstream shear layer instability transition in a nozzle-generated transverse jet.more » « less
-
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
-
Abstract Monami is the synchronous waving of a submerged seagrass bed in response to unidirectional fluid flow. Here we develop a multiphase model for the dynamical instabilities and flow-driven collective motions of buoyant, deformable seagrass. We show that the impedance to flow due to the seagrass results in an unstable velocity shear layer at the canopy interface, leading to a periodic array of vortices that propagate downstream. Our simplified model, configured for unidirectional flow in a channel, provides a better understanding of the interaction between these vortices and the seagrass bed. Each passing vortex locally weakens the along-stream velocity at the canopy top, reducing the drag and allowing the deformed grass to straighten up just beneath it. This causes the grass to oscillate periodically even in the absence of water waves. Crucially, the maximal grass deflection is out of phase with the vortices. A phase diagram for the onset of instability shows its dependence on the fluid Reynolds number and an effective buoyancy parameter. Less buoyant grass is more easily deformed by the flow and forms a weaker shear layer, with smaller vortices and less material exchange across the canopy top. While higher Reynolds number leads to stronger vortices and larger waving amplitudes of the seagrass, waving amplitude is maximized at intermediate grass buoyancy. All together, our theory and computations develop an updated schematic of the instability mechanism consistent with experimental observations.more » « less