Averaging Theory of Arbitrary Order for Piecewise Smooth Differential Systems and Its Application

The averaging theory for studying periodic orbits of smooth differential systems has a long history. Whereas the averaging theory for piecewise smooth differential systems appeared only in recent years, where the unperturbed systems are smooth. When the unperturbed systems are only piecewise smooth, there is not an existing averaging theory to study existence of periodic orbits of their perturbed systems. Here we establish such a theory for one dimensional perturbed piecewise smooth periodic differential equations. Then we show how to transform planar perturbed piecewise smooth differential systems to one dimensional piecewise smooth periodic differential equations when the unperturbed planar piecewise smooth differential systems have a family of periodic orbits. Finally as application of our theory we study limit cycle bifurcation of planar piecewise differential systems which are perturbation of a \(\Sigma \)-center.     

Planar Vector Fields with a Given Set of Orbits

We determine all the \({\mathcal{C}^1}\) planar vector fields with a given set of orbits of the form y ? y(x) = 0 satisfying convenient assumptions. The case when these orbits are branches of an algebraic curve is also study. We show that if a quadratic vector field admits a unique irreducible invariant algebraic curve \({g(x, y) = \sum_{j=0}^S a_j(x) y^{S-j}= 0}\) with S branches with respect to the variable y, then the degree of the polynomial g is at most 4S.     

Beltrami Fields with a Nonconstant Proportionality Factor are Rare

We consider the existence of Beltrami fields with a nonconstant proportionality factor f in an open subset U of \({\mathbb{R}^3}\). By reformulating this problem as a constrained evolution equation on a surface, we find an explicit differential equation that f must satisfy whenever there is a nontrivial Beltrami field with this factor. This ensures that there are no nontrivial regular solutions for an open and dense set of factors f in the C^k topology, \({k\geqq 7}\). In particular, there are no nontrivial Beltrami fields whenever f has a regular level set diffeomorphic to the sphere. This provides an explanation of the helical flow paradox of Morgulis et al. (Commun Pure Appl Math 48:571–582, 1995).     

Parameterization of Invariant Manifolds for Periodic Orbits (II): A Posteriori Analysis and Computer Assisted Error Bounds

In this paper we develop mathematically rigorous computer assisted techniques for studying high order Fourier–Taylor parameterizations of local stable/unstable manifolds for hyperbolic periodic orbits of analytic vector fields. We exploit the numerical methods developed in Castelli et al. (SIAM J Appl Dyn Syst 14(1):132–167, 2015) in order to obtain a high order Fourier–Taylor series expansion of the parameterization. The main result of the present work is an a-posteriori theorem which provides mathematically rigorous error bounds. The hypotheses of the theorem are checked with computer assistance. The argument relies on a sequence of preliminary computer assisted proofs where we validate the numerical approximation of the periodic orbit, its stable/unstable normal bundles, and the jets of the manifold to some desired order M. We illustrate our method by implementing validated computations for two dimensional manifolds in the Lorenz equations in \(\mathbb {R}^3\) and a three dimensional manifold of a suspension bridge equation in \(\mathbb {R}^4\).     

Exponential Decay of the Vorticity in the Steady-State Flow of a Viscous Liquid Past a Rotating Body

We construct a Sobolev homeomorphism in dimension \({n \geqq 4,\,f \in W^{1,1}((0, 1)^n,\mathbb{R}^n)}\) such that \({J_f = {\rm det} Df > 0}\) on a set of positive measure and J _f  < 0 on a set of positive measure. It follows that there are no diffeomorphisms (or piecewise affine homeomorphisms) f _k such that \({f_k\to f}\) in \({W^{1,1}_{\rm loc}}\).     

Application of the γ-Reθ Transition Model to Simulations of the Flow Past a Circular Cylinder

Based on the finite volume method, the flow past a two-dimensional circular cylinder at a critical Reynolds number (Re = 8.5 × 10⁵) was simulated using the Navier-Stokes equations and the γ-Re_θ transition model coupled with the SST k ? ω turbulence model (hereinafter abbreviated as γ-Re_θ model). Considering the effect of free-stream turbulence intensity decay, the SST k ? ω turbulence model was modified according to the ambient source term method proposed by Spalart and Rumsey, and then the modified SST k ? ω turbulence model is coupled with the γ-Re_θ transition model (hereinafter abbreviated as γ-Re_θ-SR model). The flow past a circular cylinder at different inlet turbulence intensities were simulated by the γ-Re_θ-SR model. At last, the flow past a circular cylinder at subcritical, critical and supercritical Reynolds numbers were each simulated by the γ-Re_θ-SR model, and the three flow states were analyzed. It was found that compared with the SST k ? ω turbulence model, the γ-Re_θ model could simulate the transition of laminar to turbulent, resulting in better consistency with experimental result. Compared with the γ-Re_θ model, for relatively high inlet turbulence intensities, the γ-Re_θ-SR model could better simulate the flow past a circular cylinder; however the improvement almost diminished for relatively low inlet turbulence intensities The γ-Re_θ-SR model could well simulate the flow past a circular cylinder at subcritical, critical and supercritical Reynolds numbers.     

Numerical Simulation of Two Phase Turbulent Flow of Nanofluids in Confined Slot Impinging Jet

In this article, a numerical investigation is performed on flow and heat transfer of confined impinging slot jet, with a mixture of water and Al₂O₃ nanoparticles as the working fluid. Two-dimensional turbulent flow is considered and a constant temperature is applied on the impingement surface. The k ? ω turbulence model is used for the turbulence computations. Two-phase mixture model is implemented to study such a flow field. The governing equations are solved using the finite volume method. In order to consider the effect of obstacle angle on temperature fields in the channel, the numerical simulations were performed for different obstacle angles of 0° ? 60°. Also different geometrical parameters, volume fractions and Reynolds numbers have been considered to study the behavior of the system in terms of stagnation point, average and local Nusselt number and stream function contours. The results showed that the intensity and size of the vortex structures depend on jet- impingement surface distance ratio (H/W) and volume fraction. The maximum Nusselt number occurs at the stagnation point with the highest values at about H/W = 1. Increasing obstacle angle, from 15° to 60°, enhances the heat transfer rate. It was also revealed that the minimum value of average Nusselt number occurs in higher H/W ratios with decreasing the channel length.     

Low order nonconforming mixed finite element method for nonstationary incompressible Navier-Stokes equations

This paper studies a low order mixed finite element method (FEM) for nonstationary incompressible Navier-Stokes equations. The velocity and pressure are approximated by the nonconforming constrained Q ₁ ^rot element and the piecewise constant, respectively. The superconvergent error estimates of the velocity in the broken H ¹-norm and the pressure in the L ²-norm are obtained respectively when the exact solutions are reasonably smooth. A numerical experiment is carried out to confirm the theoretical results.     

Experimental and Theoretical Sensitivity Analysis of Turbulent Flow Past a Square Cylinder

The Siemens SGT-800 3^rd generation DLE burner fitted to an atmospheric combustion rig has been numerically investigated. Pure methane and methane enriched by 80 vol% hydrogen flames have been considered. A URANS (Unsteady Reynolds Averaged Navier-Stokes) approach was used in this study along with the k ? ω SST and the k ? ω SST-SAS models for the turbulence transport. The chemistry is coupled to the turbulent flow simulations by the use of a laminar flamelet library combined with a presumed PDF. The effect of the mesh density in the mixing and the flame region and the effect of the turbulence model and reaction rate model constant are first investigated for the methane/air flame case. The results from the k ? ω SST-SAS along with flamelet libraries are shown to be in excellent agreement with experimental data, whereas the k ? ω SST model is too dissipative and cannot capture the unsteady motion of the flame. The k ? ω SST-SAS model is used for simulation of the 80 vol% hydrogen enriched flame case without further adjusting the model constants. The global features of the hydrogen enrichment are very well captured in the simulations using the SST-SAS model. With the hydrogen enrichment the time averaged flame front location moves upstream towards the burner exit nozzle. The results are consistent with the experimental observations. The model captures the three dominant low frequency unsteady motion observed in the experiments, indicating that the URANS/LES hybrid model indeed is capable of capturing complex, time dependent, features such as an interaction between a PVC and the flame front.     

Turbulent flows in channels with injection. results of large eddy simulation and the two-equation turbulence model

Turbulent flows in channels with intense distributed injection are modeled using the large eddy method and the two-equation k-? turbulence model. The calculations are carried out for different velocities of injection from the channel walls. For a channel with one-sided injection the results of large eddy simulation are in good agreement with the measured data, whereas the calculations in accordance with the k-? model give a less convex cross-sectional velocity profile and an appreciable error in determining the surface friction coefficient on the impermeable wall and also have certain other shortcomings. In the case of two-sided injection, the results of the calculations by the large eddy method and the k-? model are in good agreement with one another and the data of physical experiments.     

Regularity of Optimal Maps on the Sphere: the Quadratic Cost and the Reflector Antenna

Building on the results of Ma et al. (in Arch. Rational Mech. Anal. 177(2), 151–183 (2005)), and of the author Loeper (in Acta Math., to appear), we study two problems of optimal transportation on the sphere: the first corresponds to the cost function d ²(x, y), where d(·, ·) is the Riemannian distance of the round sphere; the second corresponds to the cost function ?log |x ? y|, known as the reflector antenna problem. We show that in both cases, the cost-sectional curvature is uniformly positive, and establish the geometrical properties so that the results of Loeper (in Acta Math., to appear) and Ma et al. (in Arch. Rational Mech. Anal. 177(2), 151–183 (2005)) can apply: global smooth solutions exist for arbitrary smooth positive data and optimal maps are Hölder continuous under weak assumptions on the data.     

Time Periodic Traveling Curved Fronts in the Periodic Lotka–Volterra Competition–Diffusion System

This paper is concerned with time periodic traveling curved fronts for periodic Lotka–Volterra competition system with diffusion in two dimensional spatial space

$$\begin{aligned} {\left\{ \begin{array}{ll} \dfrac{\partial u_{1}}{\partial t}=\Delta u_{1} +u_{1}(x,y,t)\left( r_{1}(t)-a_{1}(t)u_{1}(x,y,t)-b_{1}(t)u_{2}(x,y,t)\right) ,\\ \dfrac{\partial u_{2}}{\partial t}=d\Delta u_{2} +u_{2}(x,y,t)\left( r_{2}(t)-a_{2}(t)u_{1}(x,y,t)-b_{2}(t)u_{2}(x,y,t)\right) , \end{array}\right. } \end{aligned}$$

where \(\Delta \) denotes \(\frac{\partial ^{2}}{\partial x^{2} }+ \frac{\partial ^{2}}{\partial y^{2} }\), \(x,y\in {\mathbb {R}}\) and \(d>0\) is a constant, the functions \(r_i(t),a_i(t)\) and \(b_i(t)\) are T-periodic and Hölder continuous. Under suitable assumptions that the corresponding kinetic system admits two stable periodic solutions (p(t), 0) and (0, q(t)), the existence, uniqueness and stability of one-dimensional traveling wave solution \(\left( \Phi _{1}(x+ct,t),\Phi _{2}(x+ct,t)\right) \) connecting two periodic solutions (p(t), 0) and (0, q(t)) have been established by Bao and Wang ( J Differ Equ 255:2402–2435, 2013) recently. In this paper we continue to investigate two-dimensional traveling wave solutions of the above system under the same assumptions. First, we establish the asymptotic behaviors of one-dimensional traveling wave solutions of the system at infinity. Using these asymptotic behaviors, we then construct appropriate super- and subsolutions and prove the existence and non-existence of two-dimensional time periodic traveling curved fronts. Finally, we show that the time periodic traveling curved front is asymptotically stable.

     

Steady Nearly Incompressible Vector Fields in Two-Dimension: Chain Rule and Renormalization

Given bounded vector field \({b : {\mathbb{R}^{d}} \to {\mathbb{R}^{d}}}\), scalar field \({u : {\mathbb{R}^{d}} \to {\mathbb{R}}}\), and a smooth function \({\beta : {\mathbb{R}} \to {\mathbb{R}}}\), we study the characterization of the distribution \({{\rm div}(\beta(u)b)}\) in terms of div b and div(ub). In the case of BV vector fields b (and under some further assumptions), such characterization was obtained by L. Ambrosio, C. De Lellis and J. Malý, up to an error term which is a measure concentrated on the so-called tangential set of b. We answer some questions posed in their paper concerning the properties of this term. In particular, we construct a nearly incompressible BV vector field b and a bounded function u for which this term is nonzero. For steady nearly incompressible vector fields b (and under some further assumptions), in the case when d = 2, we provide complete characterization of div(\({\beta(u)b}\)) in terms of div b and div(ub). Our approach relies on the structure of level sets of Lipschitz functions on \({{\mathbb{R}^{2}}}\) obtained by G. Alberti, S. Bianchini and G. Crippa. Extending our technique, we obtain new sufficient conditions when any bounded weak solution u of \({\partial_t u + b \cdot \nabla u=0}\) is renormalized, that is when it also solves \({\partial_t \beta(u) + b \cdot \nabla \beta(u)=0}\) for any smooth function \({\beta \colon{\mathbb{R}} \to {\mathbb{R}}}\). As a consequence, we obtain new a uniqueness result for this equation.     

Wall slip and multi-tier yielding in capillary suspensions

Impact of wall slip on the yield stress measurement is examined for capillary suspensions consisting of cocoa powder as the dispersed phase, vegetable oil as the continuous primary fluid, and water as the secondary fluid using smooth and serrated parallel plates. Using dynamic oscillatory measurements, we investigated the yielding behavior of this ternary solid-fluid-fluid system with varying particle volume fraction, ?, from 0.45 to 0.65 and varying water volume fraction, ?_w, from 0.02 to 0.08. Yield stress is defined as the maximum in the elastic stress (G^′γ), which is obtained by plotting the product of elastic modulus (G^′) and strain amplitude (γ) as a function of applied strain amplitude. With serrated plates, which offer minimal slippage, capillary suspensions with ? ≥?0.45 and a fixed ?_w =?0.06 showed a two-step yielding behavior as indicated by two peaks in the plots of elastic stress as a function of strain amplitude. On the other hand with smooth plates, the capillary suspensions showed strong evidence of wall slip as evident by the presence of three distinct peaks and lowered first yield stresses for all ? and ?_w. These results can be interpreted based on the fact that a particle-depleted layer, which is known to be responsible for slip, is present in the vicinity of the smooth surfaces. The slip layer presents itself as an additional “pseudo-microstructure” (characteristic length scale) besides the two microstructures, aqueous bridges and solid particle agglomerates, that may occur in the system. With serrated plates, both the yield stresses (σ₁σ₂) and storage moduli plateau at lower strain (before the first yield point) and at higher strain (before the second yield point) (G\(^{\prime }_{p1}\), G\(^{\prime }_{p2}\)) were found to increase with ? (at a fixed ?_w =?0.06) following power-law dependences. Similarly with increasing ?_w (0.02 – 0.08) at a fixed ? =?0.62, the system behaved as a solid-like material in a jammed state with particles strongly held together as manifested by rapidly increasing σ₁ and σ₂. The usage of smooth surfaces primarily affected σ₁ which was reflected by an approximately 70–90% decrement in the measured σ₁ for all values of ?. By contrast, σ₂ and G\(^{\prime }_{p2}\) were found to be unaffected as shown by close agreement of values obtained using serrated geometry due to vanishing slip layers at higher strain amplitudes.     

On the formation of hidden chaotic attractors and nested invariant tori in the Sprott A system

We consider the well-known Sprott A system, which depends on a single real parameter a and, for \(a=1\), was shown to present a hidden chaotic attractor. We study the formation of hidden chaotic attractors as well as the formation of nested invariant tori in this system, performing a bifurcation analysis by varying the parameter a. We prove that, for \(a=0\), the Sprott A system has a line of equilibria in the z-axis, the phase space is foliated by concentric invariant spheres with two equilibrium points located at the south and north poles, and each one of these spheres is filled by heteroclinic orbits of south pole–north pole type. For \(a\ne 0\), the spheres are no longer invariant algebraic surfaces and the heteroclinic orbits are destroyed. We do a detailed numerical study for \(a>0\) small, showing that small nested invariant tori and a limit set, which encompasses these tori and is the \(\alpha \)- and \(\omega \)-limit set of almost all orbits in the phase space, are formed in a neighborhood of the origin. As the parameter a increases, this limit set evolves into a hidden chaotic attractor, which coexists with the nested invariant tori. In particular, we find hidden chaotic attractors for \(a<1\). Furthermore, we make a global analysis of Sprott A system, including the dynamics at infinity via the Poincaré compactification, showing that for \(a>0\), the only orbit which escapes to infinity is the one contained in the z-axis and all other orbits are either homoclinic to a limit set (or to a hidden chaotic attractor, depending on the value of a), or contained on an invariant torus, depending on the initial condition considered.     

Bi-parametric topology of subharmonics of an asymmetric bubble oscillator at high dissipation rate

The subharmonic topology of a nonlinear, asymmetric bubble oscillator (Keller–Miksis equation) in glycerine is investigated in the parameter space of its external excitation (frequency and pressure amplitude). The bi-parametric investigation revealed that the exoskeleton of the topology can be described as the composition of U-shaped subharmonics of periodic orbits. The fine substructure (higher-order ultra-subharmonic resonances) usually appearing via the well-known period n-tupling phenomenon is completely missing due to the high dissipation rate of the viscous liquid. Moreover, a complex internal structure of the subharmonics has been found, which are composed by interconnected bifurcation blocks (in a zig-zag pattern) each describing the skeleton of a shrimp-shaped domain. The employed numerical techniques are the combination of an in-house initial value problem solver written in C++/CUDA C to harness the high processing power of professional graphics cards, and the boundary value problem solver AUTO to compute periodic orbits directly regardless of their stability.     

Periodic solutions of <Emphasis Type="Italic">n</Emphasis>-dofs autonomous vibro-impact oscillators with one lasting contact phase

In the continuum structural mechanics framework, a unilateral contact condition between two flexible bodies does not generate impulsive contact forces. However, finite-dimensional systems, derived from a finite element semi-discretization in space for instance, and undergoing a unilateral contact condition, require an additional impact law: Unilateral contact occurrences then become impacts of zero duration unless (i) the impact law is purely inelastic, or (ii) the pre-impact velocity is zero. This contribution explores autonomous periodic solutions with one contact phase per period and zero pre-impact velocity [case (ii)], for any n-dof mechanical systems involving linear free-flight dynamics together with a linear unilateral contact constraint. A recent work has shown that such solutions seem to be limits of periodic trajectories with k impacts per period as k increases. Minimal analytic equations governing the existence of such solutions are proposed, and it is proven that, generically, they occur only for discrete values of the period. It is also shown that the graphs of such periodic solutions have two axes of symmetry in time. Results are illustrated on a spring–mass system and on a 4-dof two-dimensional system made of 1D finite elements. Animations of SPPs with up to 30 dofs are provided.     

Boundary Stabilization of Stokes System in Exterior Domains

The problem of stabilizing a solution to the 2D Stokes system defined in the exterior of the bounded domain with smooth boundary is investigated, i.e. for a given initial velocity field and prescribed positive number k > 0 one has to construct a control function defined on the boundary such that the solution stabilizes to zero at the rate of 1/t ^k .     

Two-Scale <Emphasis Type="Italic">Γ</Emphasis>-Convergence of Integral Functionals and its Application to Homogenisation of Nonlinear High-Contrast Periodic Composites

An analytical framework is developed for passing to the homogenisation limit in (not necessarily convex) variational problems for composites whose material properties oscillate with a small period ε and that exhibit high contrast of order \({\varepsilon^{-1}}\) between the constitutive, “stress-strain”, response on different parts of the period cell. The approach of this article is based on the concept of “two-scale Γ-convergence”, which is a kind of “hybrid” of the classical Γ-convergence (De Giorgi and Franzoni in Atti Accad Naz Lincei Rend Cl Sci Fis Mat Natur (8)58:842–850, 1975) and the more recent two-scale convergence (Nguetseng in SIAM J Math Anal 20:608–623, 1989). The present study focuses on a basic high-contrast model, where “soft” inclusions are embedded in a “stiff” matrix. It is shown that the standard Γ-convergence in the L ^p -space fails to yield the correct limit problem as \({\varepsilon \to 0,}\) due to the underlying lack of L ^p -compactness for minimising sequences. Using an appropriate two-scale compactness statement as an alternative starting point, the two-scale Γ-limit of the original family of functionals is determined via a combination of techniques from classical homogenisation, the theory of quasiconvex functions and multiscale analysis. The related result can be thought of as a “non-classical” two-scale extension of the well-known theorem by Müller (Arch Rational Mech Anal 99:189–212, 1987).