The flow in a Hele-Shaw cell with a time-increasing gap poses a unique shrinking interface problem. When the upper plate of the cell is lifted perpendicularly at a prescribed speed, the exterior less viscous fluid penetrates the interior more viscous fluid, which generates complex, time-dependent interfacial patterns through the Saffman–Taylor instability. The pattern formation process sensitively depends on the lifting speed and is still not fully understood. For some lifting speeds, such as linear or exponential speed, the instability is transient and the interface eventually shrinks as a circle. However, linear stability analysis suggests there exist shape invariant shrinking patterns if the gap $b(t)$ is increased more rapidly: $b(t)=\left (1-({7}/{2})\tau \mathcal {C} t\right )^{-{2}/{7}}$ , where $\tau$ is the surface tension and $\mathcal {C}$ is a function of the interface perturbation mode $k$ . Here, we use a spectrally accurate boundary integral method together with an efficient time adaptive rescaling scheme, which for the first time makes it possible to explore the nonlinear limiting dynamical behaviour of a vanishing interface. When the gap is increased at a constant rate, our numerical results quantitatively agree with experimental observations (Nase et al. , Phys. Fluids , vol. 23, 2011, 123101). When we use the shape invariant gap $b(t)$ , our nonlinear results reveal the existence of $k$ -fold dominant, one-dimensional, web-like networks, where the fractal dimension is reduced to almost unity at late times. We conclude by constructing a morphology diagram for pattern selection that relates the dominant mode $k$ of the vanishing interface and the control parameter $\mathcal {C}$ .
more »
« less
Stability and optimal decay for a system of 3D anisotropic Boussinesq equations
Abstract This paper focuses on a system of three-dimensional (3D) Boussinesq equations modeling anisotropic buoyancy-driven fluids. The goal here is to solve the stability and large-time behavior problem on perturbations near the hydrostatic balance, a prominent equilibrium in fluid dynamics, atmospherics and astrophysics. Due to the lack of the vertical kinematic dissipation and the horizontal thermal diffusion, this stability problem is difficult. When the spatial domain is Ω = R 2 × T with T = [ − 1 / 2 , 1 / 2 ] being a 1D periodic box, this paper establishes the desired stability for fluids with certain symmetries. The approach here is to distinguish the vertical averages of the velocity and temperature from their corresponding oscillation parts. In addition, the oscillation parts are shown to decay exponentially to zero in time.
more »
« less
- Award ID(s):
- 2104682
- NSF-PAR ID:
- 10322510
- Date Published:
- Journal Name:
- Nonlinearity
- Volume:
- 34
- Issue:
- 8
- ISSN:
- 0951-7715
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Solovjovs, Sergejs (Ed.)In the present paper, we summarize the results of the research devoted to the problem of stability of the fluid flow moving in a channel with flexible walls and interacting with the walls. The walls of the vessel are subject to traveling waves. Experimental data show that the energy of the flowing fluid can be transferred and consumed by the structure (the walls), inducing “traveling wave flutter.” The problem of stability of fluid-structure interaction splits into two parts: (a) stability of fluid flow in the channel with harmonically moving walls and (b) stability of solid structure participating in the energy exchange with the flow. Stability of fluid flow, the main focus of the research, is obtained by solving the initial boundary value problem for the stream function. The main findings of the paper are the following: (i) rigorous formulation of the initial boundary problem for the stream function, ψ x , y , t , induced by the fluid-structure interaction model, which takes into account the axisymmetric pattern of the flow and “no-slip” condition near the channel walls; (ii) application of a double integral transformation (the Fourier transformation and Laplace transformation) to both the equation and boundary and initial conditions, which reduces the original partial differential equation to a parameter-dependent ordinary differential equation; (iii) derivation of the explicit formula for the Fourier transform of the stream function, ψ ˜ k , y , t ; (iv) evaluation of the inverse Fourier transform of ψ ˜ k , y , t and proving that reconstruction of ψ x , y , t can be obtained through a limiting process in the complex k -plane, which allows us to use the Residue theorem and represent the solution in the form of an infinite series of residues. The result of this research is an analytical solution describing blood flowing through a channel with flexible walls that are being perturbed in the form of a traveling wave.more » « less
-
more » « less
In tropical cyclones (TCs), the peak wind speed is typically found near the top of the boundary layer (approximately 0.5–1 km). Recently, it was shown that in a few observed TCs, the wind speed within the eyewall can increase with height within the midtroposphere, resulting in a secondary local maximum at 4–5 km. This study presents additional evidence of such an atypical structure, using dropsonde and Doppler radar observations from Hurricane Patricia (2015). Near peak intensity, Patricia exhibited an absolute wind speed maximum at 5–6-km height, along with a weaker boundary layer maximum. Idealized simulations and a diagnostic boundary layer model are used to investigate the dynamics that result in these atypical wind profiles, which only occur in TCs that are very intense (surface wind speed > 50 m s−1) and/or very small (radius of maximum winds < 20 km). The existence of multiple maxima in wind speed is a consequence of an inertial oscillation that is driven ultimately by surface friction. The vertical oscillation in the radial velocity results in a series of unbalanced tangential wind jets, whose magnitude and structure can manifest as a midlevel wind speed maximum. The wavelength of the inertial oscillation increases with vertical mixing length l∞in a turbulence parameterization, and no midlevel wind speed maximum occurs when l∞is large. Consistent with theory, the wavelength in the simulations scales with (2 K/ I)1/2, where K is the (vertical) turbulent diffusivity, and I2is the inertial stability. This scaling is used to explain why only small and/or strong TCs exhibit midlevel wind speed maxima.
-
Vibrational control is an open loop stabilization technique via the application of highamplitude, high-frequency oscillatory inputs. The averaging theory has been the standard technique for designing vibrational control systems. However, it stipulates too high oscillation frequency that may not be practically feasible. Therefore, although vibrational control is very robust and elegant (stabilization without feedback), it is rarely used in practical applications. The only well-known example is the Kapitza pendulum; an inverted pendulum shose pivot is subject to vertical oscillation. the unstable equilibrium of the inverted pendulum gains asymptotic stability due to the high-frequency oscillation of the pivot. In this paper, we provide a new vibrational control system from Nature; flapping flight dynamics. Flapping flight is a rich dynamical system as a representative model will typically be nonlinear, time-varying, multi-body, multi-time-scale dynamical system. Over the last two decades, using direct averaging, there has been consensus in the flapping flight dynamics community that insects are unstable at the hovering equilibrium due to the lack of pitch stiffness. In this work, we perform higher-order averaging of the time-periodic dynamics of flapping flight to show a vibrational control mechanism due to the oscillation of the driving aerodynamic forces. We also experimentally demonstrate such a phenomenon on a flapping apparatus that has two degrees of freedom: forward translation and pitching motion. It is found that the time-periodic dynamics of the flapping micro-air-vehicle is naturally (without feedback) stabilized beyond a certain threshold. Moreover, if the averaged aerodynamic thrust force is produced by a propeller revolving at a constant speed while maintaining the wings stationary at their mean positions, no stabilization is observed. Hence, it is concluded that the observed stabilization in the flapping system at high frequencies is due to the oscillation of the driving aerodynamic force and, as such, flapping flight indeed enjoys vibrational stabilization.more » « less
-
more » « less
Abstract Very long period (VLP) seismic events (with dominant periods of 15 to 40 s), observed from 2007 to 2018 at the summit of Kīlauea Volcano, Hawai‘i, arise from resonant oscillations in the shallow magma plumbing system. Utilizing an oscillation model developed in the companion paper (Liang et al., 2020), we perform Bayesian inversions on seismic data from four representative VLP events separately for the parameters of the shallow conduit‐reservoir system, exploring both sphere and crack reservoir geometries. Both sphere and crack geometries are preferentially located ∼1–2 km beneath the northeast edge of Halema‘uma‘u crater and produce similar fits to the data. Considering a reasonable range for reservoir storativity, magma density, and density contrast between the top and bottom of the conduit, we favor a spherical reservoir with a radius of 0.8 to 1.2 km and a short conduit of less than a few hundred meters. For this geometry, buoyancy from density stratification in the conduit provides the dominant restoring force for the VLP oscillation. Viscosity is constrained within an order of magnitude for each event (e.g., approximately 2 to 23 Pa s for one event versus 27 to 513 Pa s for another). Changes in VLP period
T and quality factorQ can be explained by changes in viscosity, density stratification, and/or conduit/reservoir geometry. In particular, observed fluctuations inQ over short time intervals (e.g., hours) with minimal changes inT apparently require rapid changes of magma viscosity by over an order of magnitude, assuming geometry remains unchanged, possibly reflecting changes in volatile content, bubble concentration, or conduit flow regime.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1088/1361-6544/ac08e9
