We propose a fractional-step method for the numerical solution of unsteady thermal convection in non-Newtonian fluids with temperature-dependent physical parameters. The proposed method is based on a viscosity-splitting approach, and it consists of four uncoupled steps where the convection and diffusion terms of both velocity and temperature solutions are uncoupled while a viscosity term is kept in the correction step at all times. This fractional-step method maintains the same boundary conditions imposed in the original problem for the corrected velocity solution, and it eliminates all inconsistencies related to boundary conditions for the treatment of the pressure solution. In addition, the method is unconditionally stable, and it allows the temperature to be transported by a non-divergence-free velocity field. In this case, we introduce a methodology to handle the subtle temperature convection term in the error analysis and establish full first-order error estimates for the velocity and the temperature solutions and 1/2-order estimates for the pressure solution in their appropriate norms. Three numerical examples are presented to demonstrate the theoretical results and examine the performance of the proposed method for solving unsteady thermal convection in non-Newtonian fluids. The computational results obtained for the considered examples confirm the convergence, accuracy, and applicability of the proposed time fractional-step method for unsteady thermal convection in non-Newtonian fluids.

We propose a general method for constructing a minimal cover of high-dimensional chaotic attractors by embedded unstable recurrent patterns. By minimal cover we mean a subset of available patterns such that the approximation of chaotic dynamics by a minimal cover with a predefined proximity threshold is as good as the approximation by the full available set. The proximity measure, based on the concept of a directed Hausdorff distance, can be chosen with considerable freedom and adapted to the properties of a given chaotic system. In the context of a spatiotemporally chaotic attractor of the Kuramoto–Sivashinsky system on a periodic domain, we demonstrate that the minimal cover can be faithfully constructed even when the proximity measure is defined within a subspace of dimension much smaller than the dimension of space containing the attractor. We discuss how the minimal cover can be used to provide a reduced description of the attractor structure and the dynamics on it.

In this study, we solve the fourth-order parabolic problem by combining the implicit θ-schemes in time for θ∈[12,1] with the stabilizer free weak Galerkin (SFWG) method. The semi-discrete and full-discrete numerical schemes are proposed. And specifically, the full-discrete scheme is a first-order backward Euler scheme when θ=1, and a second-order Crank–Nicolson scheme for θ=12. Then, we determine the optimal convergence orders of the error in the H2 and L2 norms after analyzing the well-posedness of the schemes. The theoretical findings are validated by numerical experiments.

In this paper we study a nonlinear infection viral propagation model with diffusion, in which, the left boundary is fixed and with homogeneous Dirichlet boundary conditions, while the right boundary is free. We find that the habitat always expands to the half line [0,∞), and that the virus and infected cells always die out when the Basic Reproduction NumberR0≤1, while the virus and infected cells have persistence properties when R0>1. To obtain the persistence properties of virus and infected cells when R0>1, the most work of this paper focuses on the existence and uniqueness of bounded positive equilibrium solutions for subsystems and the existence of positive equilibrium solutions for the entire system.

In this paper, we probe the pth moment exponential stability (p−ES) of stochastic delayed systems subject to event-triggered delayed impulsive control (ETDIC), where the impulsive intensities are assumed to be positive random variables. Based on event-triggered mechanism (ETM) in the sense of expectation, some new sufficient conditions are developed to ensure the stability of the addressed system with Zeno-free behavior. The Lyapunov–Razumikhin method is adopted to handle the time-varying delay in continuous dynamics, and the concept of average random impulsive estimation (ARIE) is introduced to reduce the design requirements of the controller. Especially, the proposed ETDIC strategy not only generates the impulse time sequence according to the predesigned ETM, but also removes the limitations on the size of time delays. Furthermore, the ETM serves as a solution for synchronization problems in complex neural networks. Finally, two examples are given to illustrate the effectiveness of our conclusions.

In the Graph-let based time series analysis, a time series is mapped into a series of graph-lets, representing the local states respectively. The bridges between successive graph-lets are reduced simply to a linkage with an information of occurrence. In the present work, we focus our attention on the bridge series, i.e., preserve the structures of the bridges and reduce the states into nodes. The bridge series can tell us how the system evolves. Technically, the ordinal partition algorithm is adopted to construct the graph-lets and the bridges. Results for the Logistic Map, the Hénon Map, and the Lorenz System show that the statistical properties for transition frequency network for the bridges, e.g., the number of visited bridges and the average out-entropy-degree, have the capability of characterizing chaotic processes, being equivalent with the Lyapunov exponent. What is more, the topological structure can display the details of the contributions of the transitions between the bridges to the statistical properties.

This paper considers the inhomogeneous biaxial nematic liquid crystal flow in a smooth bounded domain Ω⊂R2, where the velocity u and the orthogonal unit vector fields (m,n) admit the Dirichlet and Neumann boundary condition, respectively. By applying piecewise estimate and continuity method, we get the global existence of strong solutions, provided that the basic energy is suitably small. Our result may be regarded as an extension and improvement of Gong-Lin (2022) and Li-Liu-Zhong (2017) to the Neumann boundary condition, where the initial vacuum is allowed. Some new techniques are developed in order to deal with integral estimates caused by the boundary condition, and more complicated model.

The research on stability analysis and control design for random nonlinear systems have been greatly popularized in recent ten years, but almost no literature focuses on the fractional-order case. This paper explores the stability problem for a class of random Caputo fractional-order nonlinear systems. As a prerequisite, under the globally and the locally Lipschitz conditions, it is shown that such systems have a global unique solution with the aid of the generalized Gronwall inequality and a Picard iterative technique. By resorting to Laplace transformation and Lyapunov stability theory, some feasible conditions are established such that the considered fractional-order nonlinear systems are respectively Mittag-Leffler noise-to-state stable, Mittag-Leffler globally asymptotically stable. Then, a tracking control strategy is established for a class of random Caputo fractional-order strict-feedback systems. The feasibility analysis is addressed according to the established stability criteria. Finally, a power system and a mass–spring-damper system modeled by the random fractional-order method are employed to demonstrate the efficiency of the established analysis approach. More critically, the deficiency in the existing literatures is covered up by the current work and a set of new theories and methods in studying random Caputo fractional-order nonlinear systems is built up.

The dynamic model for perturbative longitudinal vibration of microresonators subjected to the parallel-plate electrostatic force, which can be converted into a cubic oscillator with nonlinear polynomials, is established in this manuscript. The orbits and global dynamical behaviors of the cubic oscillator at full state are studied both analytically and numerically. The expressions of homoclinic orbits and subharmonic orbits are obtained analytically by solving the Hamilton system. The scenarios of phase portraits and equilibria are given. With the Melnikov method, the critical value of chaos arising from homoclinic intersections is derived analytically. The investigation yields intriguing dynamical phenomena, including the controllable frequencies that regulate the system without inducing chaos. The conditions for the occurrence of subharmonic bifurcations of integer order are presented with the subharmonic Melnikov method. Besides, the results indicate that the system does not undergo fractional order subharmonic bifurcation and it can reach a chaotic state through a finite number of integer order subharmonic bifurcations. On the basis of theoretical analysis, some numerical simulations including time histories, phase portraits, bifurcation diagrams, Poincaré cross-sections, Lyapunov exponential spectrums and basins of attractor are given, which are consistent with theoretical results.

A heterodimensional cycle is formed by the intersection of stable and unstable manifolds of two saddle periodic orbits that have unstable manifolds of different dimensions: connecting orbits exist from one periodic orbit to the other, and vice versa. The difference in dimensions of the invariant manifolds can only be achieved in vector fields of dimension at least four. At least one of the connecting orbits of the heterodimensional cycle will necessarily be structurally unstable, meaning that is does not persist under small perturbations. Nevertheless, the theory states that the existence of a heterodimensional cycle is generally a robust phenomenon: any sufficiently close vector field (in the C1-topology) also has a heterodimensional cycle.We investigate a particular four-dimensional vector field that is known to have a heterodimensional cycle. We continue this cycle as a codimension-one invariant set in a two-parameter plane. Our investigations make extensive use of advanced numerical methods that prove to be an important tool for uncovering the dynamics and providing insight into the underlying geometric structure. We study changes in the family of connecting orbits as two parameters vary and Floquet multipliers of the periodic orbits in the heterodimensional cycle change. In particular the Floquet multipliers of one of the periodic orbits change from real positive to real negative prior to a period-doubling bifurcation. We then focus on the transitions that occur near this period-doubling bifurcation and find that it generates new families of heterodimensional cycles with different geometric properties. Our careful numerical study suggests that further two-parameter continuation of the ‘period-doubled heterodimensional cycles’ gives rise to an abundance of heterodimensional cycles of different types in the limit of a period-doubling cascade.Our results for this particular example vector field make a contribution to the emerging bifurcation theory of heterodimensional cycles. In particular, the bifurcation scenario we present can be viewed as a specific mechanism behind so-called stabilisation of a heterodimensional cycle via the embedding of one of its constituent periodic orbits into a more complex invariant set.