The present paper reports on wind tunnel tests and analyses carried out to investigate the effect of the static angle of attack on the aerodynamic stability of a twin-box bridge deck section. It is found that the critical wind speed for onset of flutter increases with increasing positive static angles (nose-up) and that this effect relates mainly to a decrease in the loss of aerodynamic stiffness. A simplified flutter analysis, linking the slope of the static moment coefficient to the increase of flutter stability for increasing positive angles of attack. It is concluded that it is desirable to design twin-box bridge deck sections to have a positive moment coefficient at zero angle of attack and a positive, decreasing moment slope for increasing nose-up angles. With these requirements fulfilled, the present study show that the critical wind speeds for onset of flutter increase with increasing angles, and ensures that the elastically supported deck will always meet the wind at ever increasing angles for increasing wind speeds.
When an evaporating water droplet is deposited on a thermally conductive substrate, the minimum temperature will be at the apex due to evaporative cooling. Consequently, density and surface tension gradients emerge within the droplet and at the droplet-gas interface, giving rise to competing flows from, respectively, the apex towards the contact line (thermal-buoyancy-driven flow) and the other way around (thermal Marangoni flow). In small droplets with a diameter below the capillary length, the thermal Marangoni effects are expected to dominate over thermal buoyancy ("thermal Rayleigh") effects. However, contrary to these theoretical predictions, our experiments mostly show a dominant circulation from the apex towards the contact line, indicating a prevailing of thermal Rayleigh convection. Furthermore, our experiments often show an unexpected asymmetric flow that persisted for several minutes. We hypothesise that a tiny amount of contaminants, commonly encountered in experiments with water/air interfaces, act as surfactants and counteract the thermal surface tension gradients at the interface and thereby promote the dominance of Rayleigh convection. Our finite element numerical simulations demonstrate that, under our specified experimental conditions, a mere 0.5% reduction in the static surface tension caused by surfactants leads to a reversal in the flow direction, compared to the theoretical prediction without contaminants. Additionally, we investigate the linear stability of the axisymmetric solutions, revealing that the presence of surfactants also affects the axial symmetry of the flow.
This research explores the development and application of the High-Order Dynamic Integration Method for solving integro-differential equations, with a specific focus on turbulent fluid dynamics. Traditional numerical methods, such as the Finite Difference Method and the Finite Volume Method, have been widely employed in fluid dynamics but struggle to accurately capture the complexities of turbulence, particularly in high Reynolds number regimes. These methods often require significant computational resources and are prone to errors in nonlinear dynamic systems. The High-Order Dynamic Integration Method addresses these challenges by integrating higher-order interpolation techniques with dynamic adaptation strategies, significantly enhancing accuracy and computational efficiency. Through rigorous numerical analysis, this method demonstrates superior performance over the Finite Difference Method and the Finite Volume Method in handling the nonlinear behaviors characteristic of turbulent flows. Furthermore, the High-Order Dynamic Integration Method achieves this without a substantial increase in computational cost, making it a highly efficient tool for simulations in computational fluid dynamics. The research validates the capabilities of the High-Order Dynamic Integration Method through a series of benchmark tests and case studies. Results indicate a marked improvement in both accuracy and stability, particularly in simulations of high-Reynolds-number flows, where traditional methods often falter. This innovative approach offers a robust and efficient alternative for solving complex fluid dynamics problems, contributing to advances in the field of numerical methods and computational fluid dynamics.
Previous studies on peristalsis, the pumping of fluid along a channel by wave-like displacements of the channel walls, have shown that the elastic properties of the channel and the peristaltic wave shape can influence the flow rate. Motivated by the oscillatory flow of cerebrospinal fluid along compliant perivascular spaces, we consider a prescribed wave motion of a single boundary which pumps fluid at small Reynolds number near an elastic half space. We investigate the relationship between flow rate and elastic deformation as a function of the fluid and solid properties. We consider transverse and longitudinal motion of the driving peristaltic wave. We find that a transverse peristaltic wave produces net forward flow and induces elastic motion in which all material points oscillate uniformly. Conversely, a longitudinal peristaltic wave produces a net backward flow and drives elastic motion which is nonuniform in the elastic solid. We use dimensional values relevant to the flow of cerebrospinal fluid and find agreement with in vivo velocity data.
Recent studies have shown that capillary waves generated by bursting of an oil drop at the water-air interface produces a daughter droplet inside the bath while part of it floats above it. Successive bursting events produce next generations of daughter droplets, gradually diminishing in size until the entire volume of oil rests atop the water-air interface. In this work, we demonstrate two different ways to modulate this process by modifying the constitution of the drop. Firstly, we introduce hydrophilic clay particles inside the parent oil drop and show that it arrests the cascade of daughter droplet generation preventing it from floating over the water-air interface. Secondly, we show that bursting behavior can be modified by a compound water-oil-air interface made of a film of oil with finite thickness and design a regime map which displays each of these outcomes. We underpin both of these demonstrations by theoretical arguments providing criteria to predict outcomes resulting therein. Lastly, all our scenarios have a direct relation to control of oil-water separation and stability of emulsified solutions in a wide variety of applications which include drug delivery, enhanced oil recovery, oil spills and food processing where a dispersed oil phase tries to separate from a continuous phase.
This research explores the modulation of Rayleigh-Plateau Instability (RPI) by adjusting the orientation angle and eccentricity of a wire within a nozzle. We demonstrate that both the angle and eccentricity significantly influence the RPI regimes. They both also influence characteristics, such as velocity, wavelength, and droplet volume. Notably, when wires are both angled and eccentric, the effect of angle prevails. Our approach includes an empirical scaling analysis, comparing gravity, surface tension, and viscosity forces on a single bead, yielding a unified empirical viscous force law, enhancing understanding of RP regime dynamics. This new framework enriches our understanding of the forces at play in RPI and provides practical insights into the manipulation of fluid dynamics in industrial applications.
Lubricant viscoelasticity arises due to a finite polymer relaxation time ($\lambda$) and can provide beneficial effects. In applications, such as bearings, gears, biological joints, etc., where the height-to-length ratio is small ($H_0 / \ell_x$) and the shear due to the wall velocity ($U_0$) is high, a simplified two-dimensional computational analysis across the channel length and height reveals a finite increase in the load carrying capacity of the film purely due to polymer elasticity. In channels with a finite length-to-width ratio, $a$, the spanwise effects can be significant, but the resulting mathematical model is computationally intensive. In this work, we propose simpler reduced-order models, namely via a (i) first-order perturbation in the Deborah number ($\lambda U_0 / H_0$), and the (ii) viscoelastic Reynolds approach extended from \textitAhmed, H., \& Biancofiore, L. (2021). A new approach for modeling viscoelastic thin film lubrication. Journal of Non-Newtonian Fluid Mechanics, 292, 104524. We predict the variation in the net vertical force exerted on the channel walls (for a fixed film height) versus increasing viscoelasticity, and the channel aspect ratio. The models predict an increase in the net force, which is zero for the Newtonian case, versus both the Deborah number and the channel aspect ratio. Interestingly, for a fixed $De$, this force varies strongly between the two limiting cases (i) $a << 1$; an infinitely wide, and (ii) $a >> 1$; an infinitely short channel channel, implying a change in the polymers response. Furthermore, we observe a different trend (i) for a spanwise varying channel, in which a peak is observed between the two limits, and (ii) for a spanwise uniform channel, where the largest load value is for $a << 1$. When $a$ is O($1$), the viscoelastic response varies strongly and spanwise effects cannot be ignored.
The heterogeneous distribution of water-repellent materials at the soil surface causes a phenomenon known as fractional wettability. This condition frequently triggers destabilization of the wetting front during water infiltration, resulting in the formation of fingered bypass flow. However, few analytical tools exist to understand and model this behavior. Moreover, existing infiltration models fail to fit certain infiltration curves that exist in experimental data. For these reasons, we introduce a novel infiltration model to simulate water infiltration under fractional wettable conditions. We conceptualize the soil surface as a composite of two distinct portions: a water-repellent fraction, where hydrophobic effects impede water infiltration, and a wettable fraction, where capillarity and gravity are the dominant forces controlling the process. The new model was validated using a dataset comprising infiltration data from 60 field measurements. Additionally, validation was performed using 660 analytically generated infiltration curves from six synthetic soils with varying textures. This innovative approach enabled us to account for the combined influence of these two fractions and to enhance the interpretation of infiltration curves with mixed shapes, which other common methods are unable to reproduce.
This study presents an efficient algebraic scheme known as MULES for sharp interface advection, verified against various schemes including first-order upwind, second-order central, van Leer flux limiter, and Geometric Volume-of-Fluid (VOF). Two problems involving a droplet in a two-dimensional (2D) vortex and a stationary droplet were examined. The model assessed the effects of the Interface Compression (IC) coefficient, ranging from 0 to 2, analyzing parameters such as Interface Advection Error (IAE) and Mass Conservation Error (MCE). Results indicated that increasing IC values enhanced interface tracking accuracy but introduced non-physical instabilities at higher values, compromising mass conservation. Specifically, the IAE decreased from 4.8% to 3.95% as IC increased from 0 to 2, showing a favorable effect until IC surpassed 1.4, where IAE fluctuated around 4%. Conversely, the MCE rose steeply from 0% to 23.19%, driven by parasitic currents and numerical instabilities. Additionally, MULES and van Leer flux limiter schemes evaluated volume fraction smoothing effects. Initial filtering reduced Dimensionless Pressure Difference (DPD) and Capillary Number (Ca), stabilizing the solution, but excessive filtering reintroduced numerical errors and instabilities. With one filtering step, DPD reduced by 0.23 and Ca dropped significantly by 73.31%, improving solution stability. However, further filtering increased DPD and Ca, reflecting the reintroduction of numerical errors. The maximum velocity of parasitic flow around the droplet initially decreased by almost 75% but increased by 30.92% with excessive filtering. IAE increased from 0.7 to 0.9 with initial filtering, then decreased to 0.63 with additional steps, indicating improved solver performance on smoother interfaces.
The study of Large-Eddy Simulations (LES) in turbulent flows continues to be a critical area of research, particularly in understanding the behavior of small-scale turbulence structures and their impact on resolved scales. In this study, we focus on the LES of turbulent flows, particularly the one-dimensional Stochastic Burgers Equation (SBE), using fully conservative higher-order schemes. The interaction between spatial discretization and SubGrid-Scale (SGS) modeling is explored rigorously by validating these schemes against analytical solutions for both the Linear Advection-Diffusion (LAD) equation and the Non-Linear Burgers (NLB) equation under laminar conditions. This ensures robustness before applying the approach to LES of stochastic turbulence. The study investigates how second-order and fourth-order discretization schemes influence the dynamic coefficients of various SGS models, including Constant Smagorinsky (CS), Dynamic Smagorinsky (DS), Dynamic Wong-Lilly (DWL), 1.5-order Turbulent Kinetic Energy Deardorff (TKED), Equilibrium Heinz (EH), and Dynamic Heinz (DH) models. The second-order scheme was found to amplify fluctuations in dynamic SGS coefficients due to its higher numerical dissipation, contrasting with the more stable behavior observed with the fourth order scheme, which better captures resolved scales and results in smaller dynamic coefficients. Despite inherent differences in SGS models and discretization schemes, the final velocity distributions in one dimensional turbulence simulations were remarkably consistent, suggesting a limited influence of SGS modeling on large-scale structures in simplified turbulence scenarios. However, notable variations in resolved-scale kinetic energy were uncovered, emphasizing the importance of accurately capturing small-scale turbulence structures for precise energy dissipation predictions in LES.
This paper introduces the Parsimonious Dynamic Mode Decomposition (parsDMD), a novel algorithm designed to automatically select an optimally sparse subset of dynamic modes for both spatiotemporal and purely temporal data. By incorporating time-delay embedding and leveraging Orthogonal Matching Pursuit (OMP), parsDMD ensures robustness against noise and effectively handles complex, nonlinear dynamics. The algorithm is validated on a diverse range of datasets, including standing wave signals, identifying hidden dynamics, fluid dynamics simulations (flow past a cylinder and transonic buffet), and atmospheric sea-surface temperature (SST) data. ParsDMD addresses a significant limitation of the traditional sparsity-promoting DMD (spDMD), which requires manual tuning of sparsity parameters through a rigorous trial-and-error process to balance between single-mode and all-mode solutions. In contrast, parsDMD autonomously determines the optimally sparse subset of modes without user intervention, while maintaining minimal computational complexity. Comparative analyses demonstrate that parsDMD consistently outperforms spDMD by providing more accurate mode identification and effective reconstruction in noisy environments. These advantages render parsDMD an effective tool for real-time diagnostics, forecasting, and reduced-order model construction across various disciplines.
An implicit coupling framework between hypersonic nonequilibrium flows and material thermal response is proposed for the numerical simulation of ablative thermal protection materials during its flight trajectory. Charring ablative materials, when subjected to aerodynamic heating from hypersonic flows, undergo complex processes such as ablation and pyrolysis, involving heterogeneous and homogeneous chemical reactions. These multi-physical phenomena are simulated by a multicomponent material thermal response (MTR) solver that takes into account the complexity of component of pyrolysis gases. The species concentrations are calculated to improve the accuracy of transport and thermophysical parameters of pyrolysis gases. The MTR solver implements implicit time integration on finite difference discretization form to achieve higher efficiency. The numerical solutions of hypersonic flows and material thermal response are coupled through a gas-surface interaction interface based on surface mass and energy balance on the ablating surface. The coupled simulation employs the dual time-step technique, which introduces pseudo time step to improve temporal accuracy. The explicit coupling mechanism updates the interfacial quantities at physical time steps, which achieves higher computational efficiency, but introduces time discretization errors and numerical oscillations of interfacial quantities. In contrast, the implicit coupling mechanism updates the interfacial quantities at pseudo time steps, which reduces the temporal discretization error and suppresses numerical oscillations, but is less efficient. In addition, a simplified ablation boundary based on steady-state ablation assumption or radiation-equilibrium assumption is proposed to approximate solid heat conduction without coupling the MTR solution, providing quasi-steady flow solutions in the presence of ablation.
Complex fluid flows are important in many real-life problems. For an in-depth understanding, new and more elaborate methods of flow description are necessary. Often experimental and numerical data are accumulated in large quantities however only simple flow properties are extracted. Here we show an approach that allows quantitative extraction of the initial dynamics for magnetic micro-convection. It analyzes the initial time dynamics of the Fourier coefficients of the 2D flow snapshots. We demonstrate its effectiveness in both numeric and experimental data. In addition, the method is refined to account for noise that can be present in real-world data. This method can serve as a tool to directly compare mathematical models with experiments.
Braneshop, and
the US National Science Foundation.
SciRate