Prediction of the dynamic oscillation threshold in a clarinet model with a linearly increasing blowing pressure

Reed instruments are modeled as self-sustained oscillators driven by the pressure inside the mouth of the musician. A set of nonlinear equations connects the control parameters (mouth pressure, lip force) to the system output, hereby considered as the mouthpiece pressure. Clarinets can then be studied as dynamical systems; their steady behavior being dictated uniquely by the values of the control parameters. Considering the resonator as a lossless straight cylinder is a dramatic yet common simplification that allows for simulations using nonlinear iterative maps. This paper investigates analytically the effect of a linearly increasing blowing pressure on the behavior of this simplified clarinet model. When the control parameter varies, results from the so-called dynamic bifurcation theory are required to properly analyze the system. This study highlights the phenomenon of bifurcation delay and defines a new quantity, the dynamic oscillation threshold. A theoretical estimation of the dynamic oscillation threshold is proposed and compared with numerical simulations.     

Bifurcations and stability of a discrete singular bioeconomic system

The paper analyzes the stability and bifurcations of a discrete singular bioeconomic system in the closed first quadrant $R_{+}^{3}$ . First, applying the Poincaré scheme to a differential-algebraic predator–prey system where the economic interest of harvesting is taken into account, a discrete singular bioeconomic system is proposed. Then, local stability and the existing conditions of the flip bifurcation and Neimark–Sacker bifurcation around the interior equilibria of the proposed model are discussed by using the normal form of the discrete singular bioeconomic system, the center manifold theorem and the bifurcation theory, when choosing the step size δ as the parameter of the bifurcation. Finally, the results are illustrated and the complex dynamical behaviors are exhibited by computer numerical simulations.     

Stability and Neimark-Sacker bifurcation analysis for a discrete single genetic negative feedback autoregulatory system with delay

In this paper, we deal with a discrete single genetic negative feedback autoregulatory system with delay by using Euler method. Choosing the delay $\tau $ as the bifurcation parameter and analyzing the associated characteristic equation corresponding to the unique positive fixed point, it is found that the stability of the positive equilibrium and Neimark-Sacker bifurcation may occur when $\tau $ crosses some critical values. Then the explicit formula which determines the stability, direction, and other properties of bifurcating periodic solution is derived by using the center manifold theorem and normal form theory. Finally, in order to illustrate our theoretical analysis, numerical simulations are also included in the end.     

A Classical, Nonlinear Thermodynamic Theory of Elastic Shells Based on a Single Constitutive Assumption

By integrating along a thickness-like coordinate, exact two-dimensional equations of motion and a Second Law of Thermodynamics (a Clausius-Duhem Inequality) are derived, without approximation, for a shell-like body. The theory derived is called classical because the only stress measures that appear are resultants and couples. Construction of a Virtual Power Identity automatically produces associated extensional and bending strains that are nonlinear in the deformed position $\bar{\mathbf {y}}$ of a reference surface and a rotation tensor  $\hbox{\mathversion {bold}$\mathsf{Q}$}$ . The only approximations come when the Virtual Power Identity is augmented by heating and an internal energy and the resulting expression taken as the First Law of Thermodynamics (Conservation of Energy) for the shell. A Legendre-Fenchel Transformation is introduced to remove possible ill-conditioning when certain approximations are introduced into the strain-energy density.     

Bifurcation and chaotic behavior of a discrete-time Ricardo–Malthus model

The dynamics of a discrete-time Ricardo–Malthus model obtained by numerical discretization is investigated, where the step size δ is regarded as a bifurcation parameter. It is shown that the system undergoes flip bifurcation and Neimark–Sacker bifurcation in the interior of $R^{2}_{+}$ by using the theory of center manifold and normal form. Numerical simulations are presented not only to illustrate our theoretical results, but also to exhibit the system’s complex dynamical behavior, such as the cascade of period-doubling bifurcation in orbits of period 2, 4, 8 16, period-11, 22, 28 orbits, quasiperiodic orbits, and the chaotic sets. These results reveal far richer dynamics of the discrete model compared with the continuous model. The Lyapunov exponents are numerically computed to confirm further the complexity of the dynamical behaviors.     

Three-dimensional flow in a lid-driven cavity with width-to-height ratio of 1.6

The flow in a lid-driven cavity with width-to-height ratio of 1.6 is investigated numerically and experimentally. Experimental investigation use an apparatus with a spanwise length-to-height ratio of $\Uplambda = 10.85.$ Λ = 10.85 . Increasing the Reynolds number, we experimentally find a gradual change from the quasi-two-dimensional basic flow to a three-dimensional flow pattern. The three-dimensional flow has a significant amplitude at considerably low Reynolds numbers. Streak-line photographs and PIV vector maps are presented to illustrate the structure of the finite-amplitude flow pattern. The smooth transition is in contrast to the linear instability predicted by a linear-stability analysis for a cavity with infinite span. LDV measurements confirm the absence of a distinct threshold Reynolds number and indicate an imperfect bifurcation. The deviations between experimental observations and numerical critical Reynolds number for infinite span are explained by conducting three-dimensional simulations for a finite-span geometry. A good agreement between experimental and numerical simulation is obtained. The numerical and experimental data lead to the conjecture of a premature onset of the three-dimensional flow caused by strong secondary flows which are induced by the cavity end walls. Nevertheless, the flow structure in the finite-span cavity carries the same characteristic signatures as the nonlinear flow in the corresponding infinite-length cavity. We conclude that the observed flow can be identified as the continuation of the normal mode C _e ⁴ earlier identified in a linear-stability analysis.     

Periodic solution and heteroclinic bifurcation in a predator–prey system with Allee effect and impulsive harvesting

In this article, we investigate a prey– predator model with Allee effect and state-dependent impulsive harvesting. We obtain the sufficient conditions for the existence and uniqueness of order-1 periodic solution of system (1.2) by means of the geometry theory of semicontinuous dynamic system and the method of successor function. We also obtain that system (1.2) exhibits the phenomenon of heteroclinic bifurcation about parameter $\alpha $ . The methods used in this article are novel and prove the existence of order-1 periodic solution and heteroclinic bifurcation.     

On the effect of non-smooth Coulomb damping on flutter-type self-excitation in a non-gyroscopic circulatory 2-DoF-system

This article deals with self-excited vibrations, attractivity of stationary solutions, and the corresponding bifurcation behavior of two-dimensional differential inclusions of the type $\mathbf{M}\mathbf{q}'' + \mathbf{D}\mathbf{q}' + (\mathbf{K} + \bar{\mu}\mathbf{N})\mathbf{q} \in-\mathbf{R}\operatorname{Sign}(\mathbf{q}')$ . For the smooth case R=0, the equilibrium may become unstable due to non-conservative positional forces stemming from the circulatory matrix N. This type of instability is usually referred to as flutter instability and the loss of stability is related to a Hopf bifurcation of the steady state, which occurs for a critical parameter $\bar{\mu}= \bar{\mu}_{\mathrm{crit}}$ . For R≠0, the steady state is a set of equilibria, which turns out to be attractive for all values of the bifurcation parameter $\bar{\mu}$ . Depending on $\bar{\mu}$ , the basin of attraction of the equilibrium set can be infinite or finite. The transition from an infinite to a finite basin of attraction occurs at the stability threshold $\bar{\mu}_{\mathrm{crit}}$ of the underlying smooth problem. For the finite basin of attraction, its size is proportional to the Coulomb friction and inverse-proportional to $(\bar{\mu}- \bar{\mu}_{\mathrm{crit}})$ . By adding Coulomb damping the notion of steady state stability for the smooth problem is replaced by the question whether the basin of attraction of the steady state is infinite or finite. Simultaneously, the local Hopf-bifurcation is replaced by a global bifurcation. This implies that in the presence of Coulomb damping the occurrence of self-excited vibrations can only be investigated with regard to the perturbation level.     

Numerical and experimental investigations of pseudo-shock systems in a planar nozzle: impact of bypass mass flow due to narrow gaps

During previous investigations on pseudo-shock systems, we have observed reproducible differences between measurement and simulations for the pressure distribution as well as for size and shape of the pseudo-shock system. A systematic analysis of the deviations leads to the conclusion that small gaps of $\Delta z=O(10^{-4})$  m between quartz glass side walls and metal contour of the test section are responsible for this mismatch. This paper describes a targeted experimental and numerical study of the bypass mass flow within these gaps and its interaction with the main flow. In detail, we analyze how the pressure distribution within the channel as well as the size, shape and oscillation of the pseudo-shock system are affected by the gap size. Numerical simulations are performed to display the flow inside the gaps and to reproduce and explain the experimental results. Numerical and experimental schlieren images of the pseudo-shock system are in good agreement and show that especially the structure of the primary shock is significantly altered by the presence of small gaps. Extensive unsteady flow simulations of the geometry with gaps reveal that the shear layer between subsonic gap flow and supersonic core flow is subject to a Kelvin–Helmholtz instability resulting in small pressure fluctuations. This leads to a shock oscillation with a frequency of $f= O(10^5) \hbox {s}^{-1}$ . The corresponding time scale $\tau $  (s) is 16 times higher than the characteristic time scale $\tau _\delta =\delta /U_\infty $ of the boundary layer given by the ratio of the boundary layer thickness $\delta $ directly ahead of the shock and the undisturbed free stream velocity $U_\infty $ . To assess the reliability of our numerical investigations, the paper includes a grid study as well as an extensive comparison of several RANS turbulence models and their impact on the predicted shape of pseudo-shock systems.     

Dynamics of stochastic non-Newtonian fluids driven by fractional Brownian motion with Hurst parameter H ∈ （1/4,1/2）

A two-dimensional (2D) stochastic incompressible non-Newtonian fluid driven by the genuine cylindrical fractional Brownian motion (FBM) is studied with the Hurst parameter $H \in \left( {\tfrac{1} {4},\tfrac{1} {2}} \right)$ under the Dirichlet boundary condition. The existence and regularity of the stochastic convolution corresponding to the stochastic non-Newtonian fluids are obtained by the estimate on the spectrum of the spatial differential operator and the identity of the infinite double series in the analytic number theory. The existence of the mild solution and the random attractor of a random dynamical system are then obtained for the stochastic non-Newtonian systems with $H \in \left( {\tfrac{1} {2},1} \right)$ without any additional restriction on the parameter H.     

Periodic solution of a pest management Gompertz model with impulsive state feedback control

In this paper, a new model with two state impulses is proposed for pest management. According to different thresholds, an integrated strategy of pest management is considered, that is to say if the density of the pest population reaches the lower threshold \(h_1\) at which pests cause slight damage to the forest, biological control (releasing natural enemy) will be taken to control pests; while if the density of the pest population reaches the higher threshold \(h_2\) at which pests cause serious damage to the forest, both chemical control (spraying pesticide) and biological control (releasing natural enemy) will be taken at the same time. For the model, firstly, we qualitatively analyse its singularity. Then, we investigate the existence of periodic solution by successor functions and Poincaré-Bendixson theorem and the stability of periodic solution by the stability theorem for periodic solutions of impulsive differential equations. Lastly, we use numerical simulations to illustrate our theoretical results.     

Dynamics of stochastic Lorenz–Stenflo system

This paper discusses the Lorenz–Stenflo system under the influence of L \(\acute{\hbox {e}}\) vy noise. We find conditions under which the solution to stochastic Lorenz–Stenflo system is exponentially stable. We then investigate the estimation of the global attractive set and stochastic bifurcation behavior of the stochastic Lorenz–Stenflo system. Results show that the jump noise can make the solution stable, the bounds and bifurcation to undergo change under some conditions. Numerical results show the effectiveness and advantage of our methods.     

Quasilinear and Hessian Type Equations with Exponential Reaction and Measure Data

We prove existence results concerning equations of the type \({-\Delta_pu=P(u)+\mu}\) for p > 1 and F _k [?u] = P(u) + μ with \({1 \leqq k < \frac{N}{2}}\) in a bounded domain Ω or the whole \({\mathbb{R}^N}\) , where μ is a positive Radon measure and \({P(u)\sim e^{au^\beta}}\) with a > 0 and \({\beta \geqq 1}\) . Sufficient conditions for existence are expressed in terms of the fractional maximal potential of μ. Two-sided estimates on the solutions are obtained in terms of some precise Wolff potentials of μ. Necessary conditions are obtained in terms of Orlicz capacities. We also establish existence results for a general Wolff potential equation under the form \({u={\bf W}_{\alpha, p}^R[P(u)]+f}\) in \({\mathbb{R}^N}\) , where \({0 < R \leqq \infty}\) and f is a positive integrable function.     

Biological experimental demonstration of bifurcations from bursting to spiking predicted by theoretical models

A series of bifurcations from period-1 bursting to period-1 spiking in a complex (or simple) process were observed with increasing extra-cellular potassium concentration during biological experiments on different neural pacemakers. This complex process is composed of three parts: period-adding sequences of burstings, chaotic bursting to chaotic spiking, and an inverse period-doubling bifurcation of spiking patterns. Six cases of bifurcations with complex processes distinguished by period-adding sequences with stochastic or chaotic burstings that can reach different bursting patterns, and three cases of bifurcations with simple processes, without the transition from chaotic bursting to chaotic spiking, were identified. It reveals the structures closely matching those simulated in a two-dimensional parameter space of the Hindmarsh–Rose model, by increasing one parameter \(I\) and fixing another parameter \(r\) at different values. The experimental bifurcations also resembled those simulated in a physiologically based model, the Chay model. The experimental observations not only reveal the nonlinear dynamics of the firing patterns of neural pacemakers but also provide experimental evidence of the existence of bifurcations from bursting to spiking simulated in the theoretical models.     

Multiple Blow-Up Phenomena for the Sinh-Poisson Equation

We consider the sinh-Poisson equation $$(P) _ \lambda - \Delta{u} = \lambda \, {\rm sinh} \, u \quad {\rm in} \, \Omega, \quad u = 0 \quad {\rm on} \, \partial\Omega$$ , where Ω is a smooth bounded domain in ${\mathbb{R}^2}$ and λ is a small positive parameter. If ${0 \in \Omega}$ and Ω is symmetric with respect to the origin, for any integer k if λ is small enough, we construct a family of solutions to (P) _λ , which blows up at the origin, whose positive mass is 4πk(k?1) and negative mass is 4πk(k + 1). This gives a complete answer to an open problem formulated by Jost et al. (Calc Var PDE 31(2):263–276, 2008).     

Mixed Convection Boundary-Layer Flow Along a Vertical Cylinder Embedded in a Porous Medium Filled by a Nanofluid

The steady mixed convection boundary-layer flow on a vertical circular cylinder embedded in a porous medium filled by a nanofluid is studied for both cases of a heated and a cooled cylinder. The governing system of partial differential equations is reduced to ordinary differential equations by assuming that the surface temperature of the cylinder and the velocity of the external (inviscid) flow vary linearly with the axial distance x measured from the leading edge. Solutions of the resulting ordinary differential equations for the flow and heat transfer characteristics are evaluated numerically for various values of the governing parameters, namely the nanoparticle volume fraction ${\phi}$ , the mixed convection or buoyancy parameter ?? and the curvature parameter ??. Results are presented for the specific case of copper nanoparticles. A critical value ?? _c of ?? with ?? _c <?0 is found, with the values of | ?? _c| increasing as the curvature parameter ?? or nanoparticle volume fraction ${\phi}$ is increased. Dual solutions are seen for all values of ?? >??? _c for both aiding, ?? >?0 and opposing, ?? <?0, flows. Asymptotic solutions are also determined for both the free convection limit ${(\lambda \gg 1)}$ and for large curvature parameter ${(\gamma \gg 1)}$ .     

Periodic solutions induced by an upright position of small oscillations of a sleeping symmetrical gyrostat

The aim of this paper is to provide sufficient conditions for the existence of periodic solutions emerging from an upright position of small oscillations of a sleeping symmetrical gyrostat with equations of motion being α and β parameters satisfying Δ=α ²?4β>0 and $\beta-\frac{\alpha^{2}}{2}\pm \frac{\alpha \sqrt{\varDelta }}{2}<0$ , ε a small parameter and, F ₁ and F ₂ smooth periodic maps in the variable t in resonance p:q with some of the periodic solutions of the system for ε=0, where p and q are positive integers relatively prime. The main tool used is the averaging theory.     

Universal Moduli of Continuity for Solutions to Fully Nonlinear Elliptic Equations

This paper provides universal, optimal moduli of continuity for viscosity solutions to fully nonlinear elliptic equations F(X, D ² u) =  f(X), based on the weakest and borderline integrability properties of the source function f in different scenarios. The primary result established in this work is a sharp Log-Lipschitz estimate on u based on the L ⁿ norm of f, which corresponds to optimal regularity bounds for the critical threshold case. Optimal C ^1,α regularity estimates are also delivered when ${f\in L^{n+\varepsilon}}$ . The limiting upper borderline case, ${f \in L^\infty}$ , also has transcendental importance to elliptic regularity theory and its applications. In this paper we show, under the convexity assumption on F, that ${u \in C^{1,{\rm Log-Lip}}}$ , provided f has bounded mean oscillation. Once more, such an estimate is optimal. For the lower borderline integrability condition allowed by the theory, we establish interior a priori estimates on the ${C^{0,\frac{n-2\varepsilon}{n-\varepsilon}}}$ norm of u based on the L ^n-ε norm of f, where ? is the Escauriaza universal constant. The exponent ${\frac{n-2\varepsilon}{n-\varepsilon}}$ is optimal. When the source function f lies in L ^q , n > q > n?ε, we also obtain the exact, improved sharp Hölder exponent of continuity.     

Pore-Scale Simulations of Gas Displacing Liquid in a Homogeneous Pore Network Using the Lattice Boltzmann Method

A lattice Boltzmann high-density-ratio model, which uses diffuse interface theory to describe the interfacial dynamics and was proposed originally by Lee and Liu (J Comput Phys 229:8045–8063, 2010), is extended to simulate immiscible multiphase flows in porous media. A wetting boundary treatment is proposed for concave and convex corners. The capability and accuracy of this model is first validated by simulations of equilibrium contact angle, injection of a non-wetting gas into two parallel capillary tubes, and dynamic capillary intrusion. The model is then used to simulate gas displacement of liquid in a homogenous two-dimensional pore network consisting of uniformly spaced square obstructions. The influence of capillary number (Ca), viscosity ratio ( $M$ M ), surface wettability, and Bond number (Bo) is studied systematically. In the drainage displacement, we have identified three different regimes, namely stable displacement, capillary fingering, and viscous fingering, all of which are strongly dependent upon the capillary number, viscosity ratio, and Bond number. Gas saturation generally increases with an increase in capillary number at breakthrough, whereas a slight decrease occurs when Ca is increased from $8.66\times 10^{-4}$ 8.66 × 10 - 4 to $4.33\times 10^{-3}$ 4.33 × 10 - 3 , which is associated with the viscous instability at high Ca. Increasing the viscosity ratio can enhance stability during displacement, leading to an increase in gas saturation. In the two-dimensional phase diagram, our results show that the viscous fingering regime occupies a zone markedly different from those obtained in previous numerical and experimental studies. When the surface wettability is taken into account, the residual liquid blob decreases in size with the affinity of the displacing gas to the solid surface. Increasing Bo can increase the gas saturation, and stable displacement is observed for $Bo>1$ B o > 1 because the applied gravity has a stabilizing influence on the drainage process.     

The critical tube diameter in a two reaction steps detonation: the H2/NO2 mixture case

We present experimental results on the detonability of the H₂/NO₂ mixture whose detonation exhibits a single cellular structure (λ₁) for the lean mixtures and a double cellular structure (fine cells of size λ₁ inside larger cells of size λ₂) for stoichiometric and rich mixtures. Whatever the equivalence ratio ${\phi}$ , the chemical energy is released in two successive exothermic steps of heat of reaction Q ₁ and Q ₂ (Q ₁ + Q ₂ = Q, the total heat release) and characterised (for ${\phi > 1}$ ) by two chemical lengths. The detonability is evaluated on the basis of critical conditions of self-sustained detonations transmission from a cylindrical tube of i.d. d to free space. Results show that for the critical tube diameter relationship d/λ₁ = k, with respect to the equivalence ratio ${\phi}$ ranging from 0.5 to 1.3 at ambient temperature, k is higher than the classical value 13 and its variation is rather complex. Indeed, d/λ₁ increases with ${\phi}$ from 17–18 for ${\phi = 0.5}$ to 45–50 for ${\phi = 1}$ and to 90–100 for ${\phi = 1.3}$ . The highest detonability obtained for ${\phi = 0.6}$ is explained on the basis of the highest relative contribution of the first exothermic step to the total energy Q. We conclude that, as d/λ₁ drops with Q ₂ decreasing, it should tend to 13 with the vanishing second exothermic reaction.