The existence and stability of spike equilibria in the one-dimensional Gray–Scott model on a finite domain

Results for the existence, stability, and pulse-splitting behavior of spike patterns in the one-dimensional Gray–Scott model on a finite domain in the semi-strong spike interaction regime are summarized. Conditions on the parameters for the existence of competition instabilities, synchronous oscillatory instabilities, or pulse-splitting behavior of spike patterns are given.     

Pulse-Splitting for Some Reaction-Diffusion Systems in One-Space Dimension

Pulse-splitting, or self-replication, behavior is studied for some two-component singularly perturbed reaction-diffusion systems on a one-dimensional spatial domain. For the Gierer-Meinhardt model in the weak interaction regime, characterized by asymptotically small activator and inhibitor diffusivities, a numerical approach is used to verify the key bifurcation and spectral conditions of Ei et al. [Japan. J. Indust. Appl. Math., 18, (2001)] that are believed to be essential for the occurrence of pulse-splitting in a reaction-diffusion system. The pulse-splitting that is observed here is edge-splitting, where only the spikes that are closest to the boundary are able to replicate. For the Gray–Scott model, it is shown numerically that there are two types of pulse-splitting behavior depending on the parameter regime: edge-splitting in the weak interaction regime, and a simultaneous splitting in the semi-strong interaction regime. For the semi-strong spike interaction regime, where only one of the solution components is localized, we construct several model reaction-diffusion systems where all of the pulse-splitting conditions of Ei et al. can be verified analytically, yet no pulse-splitting is observed. These examples suggest that an extra condition, referred to here as the multi-bump transition condition , is also required for pulse-splitting behavior. This condition is in fact satisfied by the Gierer–Meinhardt and Gray–Scott systems in their pulse-splitting parameter regimes.     

The Existence and Stability of Spike Equilibria in the One-Dimensional Gray–Scott Model: The Low Feed-Rate Regime

In a singularly perturbed limit of small diffusivity ɛ of one of the two chemical species, equilibrium spike solutions to the Gray–Scott (GS) model on a bounded one-dimensional domain are constructed asymptotically using the method of matched asymptotic expansions. The equilibria that are constructed are symmetric k -spike patterns where the spikes have equal heights. Two distinguished limits in terms of a dimensionless parameter in the reaction-diffusion system are considered: the low feed-rate regime and the intermediate regime. In the low feed-rate regime, the solution branches of k -spike equilibria are found to have a saddle-node bifurcation structure. The stability properties of these branches of solutions are analyzed with respect to the large eigenvalues λ in the spectrum of the linearization. These eigenvalues, which have the property that  λ= O (1)  as  ɛ→ 0  , govern the stability of the solution on an O (1) time scale. Precise conditions, in terms of the nondimensional parameters, for the stability of symmetric k -spike equilibrium solutions with respect to this class of eigenvalues are obtained. In the low feed-rate regime, it is shown that a large eigenvalue instability leads either to a competition instability, whereby certain spikes in a sequence are annihilated, or to an oscillatory instability (typically synchronous) of the spike amplitudes as a result of a Hopf bifurcation. In the intermediate regime, it is shown that only oscillatory instabilities are possible, and a scaling-law determining the onset of such instabilities is derived. Detailed numerical simulations are performed to confirm the results of the stability theory. It is also shown that there is an equivalence principle between spectral properties of the GS model in the low feed-rate regime and the Gierer–Meinhardt model of morphogenesis. Finally, our results are compared with previous analytical work on the GS model.     

Instability analysis of the split‐step Fourier method on the background of a soliton of the nonlinear Schrödinger equation

We analyze a numerical instability that occurs in the well‐known split‐step Fourier method on the background of a soliton. This instability is found to be very sensitive to small changes of the parameters of both the numerical grid and the soliton, unlike the instability of most finite‐difference schemes. Moreover, the principle of “frozen coefficients,” in which variable coefficients are treated as “locally constant” for the purpose of stability analysis, is strongly violated for the instability of the split‐step method on the soliton background. Our analysis quantitatively explains all these features. It is enabled by the fact that the period of oscillations of the unstable Fourier modes is much smaller than the width of the soliton. Our analysis is different from the von Neumann analysis in that it requires spatially growing or decaying harmonics (not localized near the boundaries) as opposed to purely oscillatory ones. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 641–669, 2012     

Nonuniform Approach to Terminal Velocity for Single Mode Rayleigh-Taylor Instability

Abstract The temporal development of a single mode Rayleigh-Taylor instability consists of three stages:the linear, free fall and terminal velocity regimes. The purpose of this paper is to report on new phenomenaobserved in the approach to terminal velocity. Our numerical study shows an unexpected nonuniform approachto terminal velocity. The nonuniformity applies especially to the spikes, which are fingers of heavy fluid fallinginto the light fluid, but it also applies to the rising bubbles of light fluid. For spikes especially our results callinto question the meaningfulness of a terminal velocity for moderate values of the Atwood number A. After ashort time period of pseudo-terminal plateau, the spike velocity increases to a significantly higher maximum,followed by a decrease. This phenomena appears to be due to a slow evolution in the shape of the spike andbubble. We find a relation between the spike (bubble) acceleration and the tip curvature. In correlation with anincrease in the spike velocity, th     

Pattern formation and control of spatiotemporal chaos in a reaction diffusion prey–predator system supplying additional food

Spatiotemporal dynamics of a predator–prey system in presence of spatial diffusion is investigated in presence of additional food exists for predators. Conditions for stability of Hopf as well as Turing patterns in a spatial domain are determined by making use of the linear stability analysis. Impact of additional food is clear from these conditions. Numerical simulation results are presented in order to validate the analytical findings. Finally numerical simulations are carried out around the steady state under zero flux boundary conditions. With the help of numerical simulations, the different types of spatial patterns (including stationary spatial pattern, oscillatory pattern, and spatiotemporal chaos) are identified in this diffusive predator–prey system in presence of additional food, depending on the quantity, quality of the additional food and the spatial domain and other parameters of the model. The key observation is that spatiotemporal chaos can be controlled supplying suitable additional food to predator. These investigations may be useful to understand complex spatiotemporal dynamics of population dynamical models in presence of additional food.     

Implementation and stability study of phase-locked-loop nonlinear dynamic measurement systems

We study the stability and convergence of a phase-locked-loop applied to a nonlinear system. It has been shown through numerical simulations by previous investigators that nonlinearity gives rise to oscillatory instability. By applying the method of averaging to the nonlinear system, we found that the nonlinear system has the identical criterion for stability as the linear system. However, the stable equilibrium has a shrinking domain of attraction as the nonlinearity increases. We show this by examining the feedback function. Moreover, we propose a nonlinear feedback which has faster convergence rate.     

Some sufficient conditions for the non-negativity preservation property in the discrete heat conduction model

In the course of the numerical approximation of mathematical models there is often a need to solve a system of linear equations with a tridiagonal or a block-tridiagonal matrices. Usually it is efficient to solve these systems using a special algorithm (tridiagonal matrix algorithm or TDMA) which takes advantage of the structure. The main result of this work is to formulate a sufficient condition for the numerical method to preserve the non-negativity for the special algorithm for structured meshes. We show that a different condition can be obtained for such cases where there is no way to fulfill this condition. Moreover, as an example, the numerical solution of the two-dimensional heat conduction equation on a rectangular domain is investigated by applying Dirichlet boundary condition and Neumann boundary condition on different parts of the boundary of the domain. For space discretization, we apply the linear finite element method, and for time discretization, the well-known Θ-method. The theoretical results of the paper are verified by several numerical experiments.     

Numerical Energy Conservation for Multi-Frequency Oscillatory Differential Equations

The long-time near-conservation of the total and oscillatory energies of numerical integrators for Hamiltonian systems with highly oscillatory solutions is studied in this paper. The numerical methods considered are symmetric trigonometric integrators and the St?rmer–Verlet method. Previously obtained results for systems with a single high frequency are extended to the multi-frequency case, and new insight into the long-time behaviour of numerical solutions is gained for resonant frequencies. The results are obtained using modulated multi-frequency Fourier expansions and the Hamiltonian-like structure of the modulation system. A brief discussion of conservation properties in the continuous problem is also included. AMS subject classification (2000) 65L05, 65P10     

一些高振荡积分、高振荡积分方程的高性能计算

电磁、声波散射问题的研究涉及一类数学物理问题, 此类问题具有深刻的理论价值和重要的应用背景, 亟待解决. 高振荡微分、积分方程是刻画这些问题的重要的数学模型, 其数值计算存在许多挑战性研究课题. 本文从积分方程解法角度出发, 综述了求解这类高振荡问题的一些最新进展, 特别是针对广义Fourier 变换、Bessel 变换的高效算法、高振荡核Volterra 积分方程的数值解法作了详细介绍. 这些数值方法共有特点是振荡频率越高算法精度愈高, 且可望为电磁计算的研究提供一些新的高效算法.     

Stress continuity in DEM-FEM multiscale coupling based on the generalized bridging domain method

Concurrent multiscale method is a spatial and temporal combination of two different scale models for describing the micro/meso and macro mixed behaviors observed in strain localization, failure and phase transformation processes, etc. Most of the existing coupling schemes use the displacement compatibility conditions to glue different scale models, which leads to displacement continuity and stress discontinuity for the obtained multiscale model. To overcome stress discontinuity, this paper presented a multiscale method based on the generalized bridging domain method for coupling the discrete element (DE) and finite element (FE) models. This coupling scheme adopted displacement and stress mixed compatibility conditions. Displacements that were interpolated from FE nodes were prescribed on the artificial boundary of DE model, while stresses at numerical integration points that were extracted from DE contact forces were applied on the material transition zone of FE model (the coupling domain and the artificial boundary of FE model). In addition, this paper proposed an explicit multiple time-steps integration algorithm and adopted Cundall nonviscous damping for quasi-static problems. DE and FE parameters were calibrated by DE simulations of a biaxial compression test and a deposition process. Numerical examples for a 2D cone penetration test (CPT) show that the proposed multiscale method captures both mesoscopic and macroscopic behaviors such as sand soil particle rearrangement, stress concentration near the cone tip, shear dilation, penetration resistance vibration and particle rotation, etc, during the cone penetration process. The proposed multiscale method is versatile for maintaining stress continuity in coupling different scale models.     

Stabilized approximation to degenerate transport equations via filtering

We analyze a stabilization technique for degenerate transport equations. Of particular interest are coupled parabolic/hyperbolic problems, when the diffusion coefficient is zero in part of the domain. The unstabilized, computed approximations of these problems are highly oscillatory, and several techniques have been proposed and analyzed to mitigate the effects of the sub-grid errors that contribute to the oscillatory behavior. In this paper, we modify a time-relaxation algorithm proposed in [1] and further studied in [10]. Our modification introduces the relaxation operator as a post-processing step. The operator is not time-dependent, so the discrete (relaxation) system need only be factored once. We provide convergence analysis for our algorithm along with numerical results for several model problems.     

Dynamics and Coarsening of Interfaces for the Viscous Cahn—Hilliard Equation in One Spatial Dimension

In one spatial dimension, the metastable dynamics and coarsening process of an n -layer pattern of internal layers is studied for the Cahn–Hilliard equation, the viscous Cahn–Hilliard equation, and the constrained Allen–Cahn equation. These models from the continuum theory of phase transitions provide a caricature of the physical process of the phase separation of a binary alloy. A homotopy parameter is used to encapsulate these three phase separation models into one parameter-dependent model. By studying a differential-algebraic system of ordinary differential equations describing the locations of the internal layers for a metastable pattern for this parameter-dependent model, we are able to provide detailed comparisons between the internal layer dynamics for the three models. Layer collapse events are studied in detail, and the analytical theory is supplemented by numerical results showing the different behaviors for the different models. Finally, an asymptotic-numerical algorithm, based on our asymptotic information of layer collapse events and the conservation of mass condition, is devised to characterize the entire coarsening process for each of these models. Numerical realizations of this algorithm are shown.     

Collective behavior of interacting locally synchronized oscillations in neuronal networks

Local circuits in the cortex and hippocampus are endowed with resonant, oscillatory firing properties which underlie oscillations in various frequency ranges (e.g. gamma range) frequently observed in the local field potentials, and in electroencephalography. Synchronized oscillations are thought to play important roles in information binding in the brain. This paper addresses the collective behavior of interacting locally synchronized oscillations in realistic neural networks. A network of five neurons is proposed in order to produce locally synchronized oscillations. The neuron models are Hindmarsh–Rose type with electrical and/or chemical couplings. We construct large-scale models using networks of such units which capture the essential features of the dynamics of cells and their connectivity patterns. The profile of the spike synchronization is then investigated considering different model parameters such as strength and ratio of excitatory/inhibitory connections. We also show that transmission time-delay might enhance the spike synchrony. The influence of spike-timing-dependence-plasticity is also studies on the spike synchronization.     

Stability of periodic traveling flexural‐gravity waves in two dimensions

In this work, we solve the Euler's equations for periodic waves traveling under a sheet of ice. These waves are referred to as flexural‐gravity waves. We compare and contrast two models for the effect of the ice: a linear model and a nonlinear model. The benefit of this reformulation is that it facilitates the asymptotic analysis. We use it to derive the nonlinear Schrödinger equation that describes the modulational instability of periodic traveling waves. We compare this asymptotic result with the numerical computation of stability using the Fourier–Floquet–Hill method to show they agree qualitatively. We show that different models have different stability regimes for large values of the flexural rigidity parameter. Numerical computations are also used to analyze high‐frequency instabilities in addition to the modulational instability. In the regions examined, these are shown to be the same regardless of the model representing ice.     

An analytical solution for conduction-dominated unidirectional solidification of binary mixtures

In this paper, a fully analytical solution technique is established for the solution of unidirectional, conduction-dominated, alloy solidification problems. By devising appropriate averaging techniques for temperature and phase-fraction gradients, governing equations inside the mushy region are made inherently homogeneous. The above formulation enables one to obtain complete analytical solutions for solid, liquid and mushy regions, without resorting to any numerical iterative procedure. Due considerations are given to account for variable properties and different microscopic models of alloy solidification (namely, equilibrium and non-equilibrium models) in the two-phase domain. The results are tested for the problem of solidification of a NH₄Cl–H₂O solution, and compared with those from existing analytical models as well as with the corresponding results from a fully numerical simulation. The effects of different microscopic models on solidification behaviour are illustrated, and transients in temperature and heat flux distribution are also analysed. A good agreement between the present solutions and results from computational simulation is observed.     

Long-wave oscillatory Marangoni instability in a horizontal liquid layer

The long-wave instability in the problem of thermocapillary convection in a horizontal layer with a free deformable boundary and a solid bottom is investigated. The transcendental equation for the main asymptotic term of the spectral parameter is written in explicit form. The main attention is paid to investigating oscillatory instability. For the frequency of neutral oscillations, simple transcendental equations are obtained that contain the Prandtl and Biot numbers. In a number of cases, exact solutions are indicated. Explicit formulae are given for the main asymptotic term of the Marangoni number. In the case of a non-heat-conducting solid wall, the relation between the critical values of the parameters for inverse Prandtl numbers is found. It is shown that, for different Prandtl numbers, the asymptotic values are in good agreement with the numerical values.     

Constitutive modeling of ferromagnetic and ferroelectric behaviors and application to multiferroic composites

In this paper, the constitutive modeling of nonlinear multifield behavior as well as the finite element implementation are presented. Nonlinear material models describing the magneto-ferroelectric or electro-ferromagnetic behaviors are presented. Both physically and phenomenologically motivated constitutive models have been developed for the numerical calculation of principally different nonlinear magnetostrictive behaviors. Further, the nonlinear ferroelectric behavior is based on a physically motivated constitutive model. On this basis, the polarization in the ferroelectric and magnetization in the ferromagnetic and magnetostrictive phases, respectively, are simulated and the resulting effects analyzed. Numerical simulations focus on the calculation of magnetoelectric coupling and on the prediction of local domain orientations going along with the poling process, thus supplying information on favorable electric-magnetic loading sequences. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dark solitons in external potentials

We consider the persistence and stability of dark solitons in the Gross–Pitaevskii (GP) equation with a small decaying potential. We show that families of black solitons with zero speed originate from extremal points of an appropriately defined effective potential and persist for sufficiently small strength of the potential. We prove that families at the maximum points are generally unstable with exactly one real positive eigenvalue, while families at the minimum points are generally unstable with exactly two complex-conjugated eigenvalues with positive real part. This mechanism of destabilization of the black soliton is confirmed in numerical approximations of eigenvalues of the linearized GP equation and full numerical simulations of the nonlinear GP equation. We illustrate the monotonic instability associated with the real eigenvalues and the oscillatory instability associated with the complex eigenvalues and compare the numerical results of evolution of a dark soliton with the predictions of Newton’s particle law for its position.      

A differential quadrature analysis of unsteady open channel flow

A rapid, convergent and accurate differential quadrature method (DQM) is employed for numerical simulation of unsteady open channel flow. To the best of authors’ knowledge, this is the first attempt to use the DQM in open channel hydraulics. The Saint-Venant equations and the related nonhomogenous, time dependent boundary conditions are discretized in spatial and temporal domain by DQ rules. The unknowns in the entire domain are computed by satisfying governing equations, boundary and initial conditions simultaneously. By employing DQM, accurate results can be obtained using dramatically less grid points in spatial and time domain. The stability of DQM solution is not sensitive to choosing time step or Courant number unlike other methods. Although numerical problems such as instability, oscillation and underestimation near critical depth can be seen by using other methods but DQM solution is smooth and accurate in this case. The results are sensitive to grid distribution in time domain. In light of this, Chebyshev–Gauss–Lobatto distribution performance is excellent. To validate the DQM solutions, the obtained results are compared with those of the characteristic method. In conclusion, DQM is a potential powerful method with minimum computational effort for unsteady flow simulation.