Multiscale Expansion Research Articles

A nonlocal Burgers equation in atmospheric dynamical system and its exact solutions**Project supported by the National Natural Science Foundation of China (Grant Nos. 11405110, 11275129, and 11472177) and the Natural Science Foundationof Zhejiang Province of China (Grant No. LY18A050001).

From a two-vortex interaction model in atmospheric and oceanic systems, a nonlocal counterpart with shifted parity and delayed time reversal is derived by a simple AB reduction. To obtain some approximate analytic solutions of this nonlocal system, the multi-scale expansion method is applied to get an AB-Burgers system. Various exact solutions of the AB-Burgers equation, including elliptic periodic waves, kink waves and solitary waves, are obtained and shown graphically. To show the applications of these solutions in describing correlated events, a simple approximate solution for the two-vortex interaction model is given to show two correlated dipole blocking events at two different places. Furthermore, symmetry reduction solutions of the nonlocal AB-Burgers equation are also given by using the standard Lie symmetry method.

Large-Scale Compressibility in Rotating Flows of Astrophysical Plasma in the Shallow Water Approximation

Two systems of magnetohydrodynamic equations in the shallow water approximation are proposed as a basis for studies in the field of plasma astrophysics: the system of equations with a full allowance for the Coriolis force and the system of equations on the β-plane in which the changes of the Coriolis parameter are linear in coordinate. Both systems of equations take into account such fundamentally important phenomena in plasma astrophysics as the compressibility and external magnetic field effects, increasing significantly the potential for applying these equations to study astrophysical objects. Compressibility in plasma astrophysics is shown to change significantly the dispersion laws for magneto-Poincare, magnetostrophic, and magneto-Rossby waves. The same nonlinear interactions as those in the absence of compressibility have been found to be realized in the case of a compressible rotating plasma. Three-wave equations in the weak nonlinearity approximation, in which the interaction coefficients depend on plasma large-scale compressibility and thermodynamic characteristics, have been derived by the method of multiscale expansions. Expressions for the growth rates of the parametric instabilities of three-wave interactions with large-scale compressibility have been derived.

Chemo-mechanical coupling analysis of composites by second-order multiscale asymptotic expansion

A novel second-order multiscale asymptotic expansion is developed in this work to analyze chemo-mechanical properties of composites with periodic microstructure. The chemo-mechanical coupling model which considers mutual interaction between the displacement and concentration fields of composites is introduced at first. Then, the second-order multiscale formulas based on homogenization methods and multiscale asymptotic expansion for evaluating the chemo-mechanical coupling problems are proposed, including the microscale cell functions, homogenized coefficients and homogenized equations. Further, the related numerical algorithms on the basis of the proposed multiscale models are brought forward. Finally, by some representative examples, the efficiency and accuracy of the presented algorithms are verified. The numerical results clearly illustrate that the second-order multiscale methods proposed in this work are effective and valid to predict the chemo-mechanical coupling properties, and can capture the microscale behavior of the composites accurately.

Speech and music classification using spectrogram based statistical descriptors and extreme learning machine

This article proposes a novel feature extraction approach for speech/music classification based on generalized Gaussian distribution descriptors extracted from IIR-CQT spectrogram representation. IIR-CQT spectrogram visual representation provides superior temporal resolution at high frequencies and better spectral resolution for low frequencies compared to the conventional short-time Fourier transform analysis which provides uniform frequency resolution. Multi-level decomposition of the spectrogram image is then performed using the Nonsubsampled Contourlet Transform (NSCT) which a fully shift-invariant, multi-scale, and multi-direction expansion that can preserve the edges of the textural pattern of speech and music. The generalized Gaussian distribution (GGD) parameters are produced using maximum likelihood estimation (MLE) from the NSCT subbands to create the image feature descriptor. Chaos crow search algorithm is employed to chose the most relevant feature sub-set and to discard redundant features and finally the extreme learning machine classifier categorizes input audio segment into speech/music. The experimental results show that the proposed feature descriptor is effective and performs better compared to the existing approaches in the speech/music classification. In addition, mismatched training and testing results are also presented.

An inclined cable excited by a non-ideal massive moving deck: an asymptotic formulation

A moving deck is an important (kinematic) excitation source for the inclined cable in cable-stayed bridges. In ideal cases, the deck motion is assumed to be harmonic oscillation and cable’s dynamic effects on the deck are neglected. As a refined version, an inclined cable excited by a massive non-ideal moving deck, i.e., the deck’s oscillation, is slowly modulated by the cable and thus not exactly harmonic is investigated in an asymptotically coupled formulation for understanding cable–deck dynamic interactions. More explicitly, by ordering the deck/cable mass ratio as a large parameter, the coupled system is reduced using asymptotic approximations and multi-scale expansions. After neglecting the reduced model’s nonlinear terms, firstly, cable–deck linear coupled modes are obtained, leading to two different kinds of linear modal dynamics, i.e., the cable-dominated one and the deck-dominated one, whose asymptotic characteristics are also revealed. Then cable’s forced nonlinear vibrations, excited by the deck’s modulated oscillation (i.e., non-ideal moving deck), are fully investigated. Nonlinear frequency responses of the cable–deck coupled system are found, and the dynamic effects on the cable’s periodic and quasi-periodic behaviors, due to cable–deck coupling (characterized by the deck/cable mass ratio), cable’s inclinations, and boundary damping, are also presented.

Solitary waves of the two-dimensional Camassa–Holm—nonlinear Schrödinger equation

In this work, we study solitary waves in a -dimensional variant of the defocusing nonlinear Schrödinger (NLS) equation, the so-called Camassa–Holm-NLS (CH-NLS) equation. We use asymptotic multiscale expansion methods to reduce this model to a Kadomtsev–Petviashvili (KP) equation. The KP model includes both the KP-I and KP-II versions, which possess line and lump soliton solutions. Using KP solitons, we construct approximate solitary wave solutions on top of the stable continuous-wave solution of the original CH-NLS model, which are found to be of both the dark and anti-dark type. We also use direct numerical simulations to investigate the validity of the approximate solutions, study their evolution, as well as their head-on collisions.

Effect of electromagnetic radiation on the dynamics of spatiotemporal patterns in memristor-based neuronal network

Modulational instability, as a mechanism of wave trains and soliton formation in biological system, is explored in the frame work of the new FitzHugh–Nagumo model. This model considered chain networks with memristive synaptic connection between adjacent neurons. This connection replaces the synaptic coupling and neurons bridged for signal exchange. From the physical law of electromagnetic induction, we interpret the traditional current term as magnetic flux variable. Magnetic flux is used to describe time-varying electromagnetic field setup in cells as a result of internal bioelectricity of the nervous system as well as when cells are exposed to external electromagnetic field. We reduced the whole network dynamical equations through multi-scale expansion to obtain a single differential–difference nonlinear equation of Schrodinger type. Linear stability is then performed with emphasis on memristive synaptic coupling. The conditions under which uniform plane waves propagating in the network become stable or unstable under small perturbation are calculated and plotted. Numerical experiments confirm our analytical predictions as the network supports localized mode excitations, spike-like, identified as quasi-periodic patterns, with some features of synchronization. It is confirmed that under strong electromagnetic radiation, the propagating waves encountered turbulent electrical activities, with patterns breakdown into a homogeneous state. This disordered state, collapse and instability of traveling pulse are monitored and analyzed using the sampled time series for membrane potential. It decreases to quiescent state under strong electromagnetic field. This could provide some guidance to understanding some neurodegenerative manifestations linked with high radiation exposure.

Hydrodynamics and two-dimensional dark lump solitons for polariton superfluids.

We study a two-dimensional incoherently pumped exciton-polariton condensate described by an open-dissipative Gross-Pitaevskii equation for the polariton dynamics coupled to a rate equation for the exciton density. Adopting a hydrodynamic approach, we use multiscale expansion methods to derive several models appearing in the context of shallow water waves with viscosity. In particular, we derive a Boussinesq/Benney-Luke-type equation and its far-field expansion in terms of Kadomtsev-Petviashvili-I (KP-I) equations for right- and left-going waves. From the KP-I model, we predict the existence of vorticity-free, weakly (algebraically) localized two-dimensional dark-lump solitons. We find that, in the presence of dissipation, dark lumps exhibit a lifetime three times larger than that of planar dark solitons. Direct numerical simulations show that dark lumps do exist, and their dissipative dynamics is well captured by our analytical approximation. It is also shown that lumplike and vortexlike structures can spontaneously be formed as a result of the transverse "snaking" instability of dark soliton stripes.

Fast and slow resonant triads in the two-layer rotating shallow water equations

In this paper, we examine triad resonances in a rotating shallow water system when there are two free interfaces. This allows for an examination in a relatively simple model of the interplay between baroclinic and barotropic dynamics in a context where there is also a geostrophic mode. In contrast to the much-studied one-layer rotating shallow water system, we find that as well as the usual slow geostrophic mode, there are now two fast waves, a barotropic mode and a baroclinic mode. This feature permits triad resonances to occur between three fast waves, with a mixture of barotropic and baroclinic modes, an aspect that cannot occur in the one-layer system. There are now also two branches of the slow geostrophic mode, with a repeated branch of the dispersion relation. The consequences are explored in a derivation of the full set of triad interaction equations, using a multiscale asymptotic expansion based on a small-amplitude parameter. The derived nonlinear interaction coefficients are confirmed using energy and enstrophy conservation. These triad interaction equations are explored, with an emphasis on the parameter regime with small Rossby and Froude numbers.

Constructive Renormalization of the 2-Dimensional Grosse–Wulkenhaar Model

We study a quartic matrix model with partition function $Z=\int d\ M\exp{\rm Tr}\ (-\Delta M^2-\frac{\lambda}{4}M^4)$. The integral is over the space of Hermitian $(\Lambda+1)\times(\Lambda+1)$ matrices, the matrix $\Delta$, which is not a multiple of the identity matrix, encodes the dynamics and $\lambda>0$ is a scalar coupling constant. We proved that the logarithm of the partition function is the Borel sum of the perturbation series, hence is a well defined analytic function of the coupling constant in certain analytic domain of $\lambda$, by using the multi-scale loop vertex expansions. All the non-planar graphs generated in the perturbation expansions have been taken care of on the same footing as the planar ones. This model is derived from the self-dual $\phi^4$ theory on the 2 dimensional Moyal space, also called the 2 dimensional Grosse-Wulkenhaar model. This would also be the first fully constructed matrix model which is non-trivial and not solvable.

Multiscale numerical algorithms for elastic wave equations with rapidly oscillating coefficients

This paper reports a multiscale analysis and numerical algorithms for the elastic wave equations with rapidly oscillating coefficients. We mainly focus on the first-order and the second-order multiscale asymptotic expansions for the wave equations, which is proved to enjoy an explicit convergence rate. In our method, the homogenized equations are discretized by the finite element method in space and a symplectic geometric scheme in time. The multiscale solutions are then obtained efficiently by the standard multisclae asymptotic expansion framework. Several numerical simulations are carried out to validate the predicted convergence results.

A three-scale homogenization algorithm for coupled conduction-radiation problems in porous materials with multiple configurations

This work establishes a novel three-scale homogenization algorithm to predict heat transfer performance of porous materials with multiple periodic configurations. The heterogeneities of porous structures are considered by periodic distributions of unit cells on the microscale and mesoscale. A new micro-meso-macro formula based on homogenization methods and multiscale asymptotic expansions is given at first. Two types of unit cell solutions in microscale and mesoscale are obtained by solving the distinct multiscale cell functions. Also, two kinds of homogenization coefficients are calculated by up-scaling procedure, and the homogenization equations are defined on global structure. Further, the temperature and heat flux fields are established as three-scale approximate solutions by assembling the various local cell solutions and homogenization solutions. Then, the associated finite element algorithm based on the three-scale homogenization methods is proposed in detail. Finally, some numerical examples are reported to validate the methods. They illustrate that the three-scale homogenization methods presented in this work are effective and accurate for calculating the heat transfer performance of the porous materials with multiple configurations.

Self-organized pattern dynamics of somitogenesis model in embryos

Somitogenesis, the sequential formation of a periodic pattern along the anteroposterior axis of vertebrate embryos, is one of the most obvious examples of the segmental patterning processes that take place during embryogenesis and also one of the major unresolved events in developmental biology. In this paper, we investigate the effect of diffusion on pattern formation use a modified two dimensional model which can be used to explain somitogenesis during embryonic development. This model is suitable for exploring a design space of somitogenesis and can explain many aspects of somitogenesis that previous models cannot. In the present paper, by analyzing the local linear stability of the equation, we acquired the conditions of Hopf bifurcation and Turing bifurcation. In addition, the amplitude equation near the Turing bifurcation point is obtained by using the methods of multi-scale expansion and symmetry analysis. By analyzing the stability of the amplitude equation, we know that there are various complex phenomena, including Spot pattern, mixture of spot–stripe patterns and labyrinthine. Finally, numerical simulation are given to verify the correctness of our theoretical results. Somitogenesis occupies an important position in the process of biological development, and as a pattern process can be used to investigate many aspects of embryogenesis. Therefore, our study helps greatly to cell differentiation, gene expression and embryonic development. What is more, it is of great significance for the diagnosis and treatment of human diseases to study the related knowledge of model biology.

Constructive tensorial group field theory II: the model

In this paper, we continue our program of non-pertubative constructions of tensorial group field theories (TGFT). We prove analyticity and Borel summability in a suitable domain of the coupling constant of the simplest super-renormalizable TGFT which contains some ultraviolet divergencies, namely the color-symmetric quartic melonic rank-four model with Abelian gauge invariance, nicknamed . We use a multiscale loop vertex expansion. It is an extension of the loop vertex expansion (the basic constructive technique for non-local theories) which is required for theories that involve non-trivial renormalization.

Three-Component Soliton States in Spinor F=1 Bose-Einstein Condensates.

Dilute-gas Bose-Einstein condensates are an exceptionally versatile test bed for the investigation of novel solitonic structures. While matter-wave solitons in one- and two-component systems have been the focus of intense research efforts, an extension to three components has never been attempted in experiments. Here, we experimentally demonstrate the existence of robust dark-bright-bright (DBB) and dark-dark-bright solitons in a multicomponent F=1 condensate. We observe lifetimes on the order of hundreds of milliseconds for these structures. Our theoretical analysis, based on a multiscale expansion method, shows that small-amplitude solitons of these types obey universal long-short wave resonant interaction models, namely, Yajima-Oikawa systems. Our experimental and analytical findings are corroborated by direct numerical simulations highlighting the persistence of, e.g., the DBB soliton states, as well as their robust oscillations in the trap.

Multiscale analysis and computation for coupled conduction, convection and radiation heat transfer problem in porous materials

This paper discusses the multiscale analysis and numerical algorithms for coupled conduction, convection and radiation heat transfer problem in periodic porous materials. First, the multiscale asymptotic expansion of the solution for the coupled problem is presented, and high-order correctors are constructed. Then, error estimates and their proofs will be given on some regularity hypothesis. Finally, the corresponding finite element algorithms based on multiscale method are introduced and some numerical results are given in detail. The numerical tests demonstrate that the developed method is feasible and valid for predicting the heat transfer performance of periodic porous materials, and support the approximate convergence results proposed in this paper.

Rigorous derivation of nonlinear Dirac equations for wave propagation in honeycomb structures

We consider a nonlinear Schrödinger equation in two spatial dimensions subject to a periodic honeycomb lattice potential. Using a multi-scale expansion together with rigorous error estimates, we derive an effective model of nonlinear Dirac type. The latter describes the propagation of slowly modulated, weakly nonlinear waves spectrally localized near a Dirac point.

A Multiscale Perturbation Expansion Approach for Markov State Modeling of Nonstationary Molecular Dynamics

We investigate metastable dynamical systems subject to nonstationary forcing as they appear in molecular dynamics for systems driven by external fields. We show that if the strength of the forcing is inversely proportional to the length of the slow metastable time scales of the unforced system, then the effective behavior of the forced system on slow time scales can be described by a low-dimensional reduced master equation. Our construction is explicit and uses the multiscale perturbation expansion method called two-timing, or method of multiple scales. The reduced master equation---a Markov state model---can be assembled by constructing two equilibrium Markov state models, one for the unforced system and one for a slightly perturbed one.

Asymptotic solution for the Darcy-Brinkman-Boussinesq flow in a pipe with helicoidal shape

We study the fluid flow and heat transfer in a helical pipe filled with a sparsely packed porous medium. Motivated by the engineering applications, pipe?s thickness and the distance between two coils of the helix have the same small order of magnitude, whereas the fluid inside the pipe is assumed to be cooled (or heated) by the exterior medium. After writing the dimensionless Darcy?Brinkman?Boussinesq system in curvilinear coordinates, we employ the multi-scale expansion technique to formally derive the effective model valid for small Brinkman?Darcy number. The obtained asymptotic solution is given in the explicit form which is important with regards to numerical simulations. Comparison with our previous results on the straight-pipe flow is also provided.

Global-in-Time Stability of 2D MHD Boundary Layer in the Prandtl--Hartmann Regime

In this paper, we prove global existence of solutions with analytic regularity to the 2D MHD boundary layer equations in the mixed Prandtl and Hartmann regime derived by formal multi-scale expansion in \cite{GP}. The analysis shows that the combined effect of the magnetic diffusivity and transveral magnetic field on the boundary leads to a linear damping on the tangential velocity field near the boundary. And this damping effect yields the global in time analytic norm estimate in the tangential space variable on the perturbation of the classical steady Hartmann profile.