Asymmetric spike patterns for the one-dimensional Gierer–Meinhardt model: equilibria and stability

2002 ◽  
Vol 13 (3) ◽  
pp. 283-320 ◽  
Author(s):  
M. J. WARD ◽  
J. WEI

Equilibrium solutions to the one-dimensional Gierer–Meinhardt model in the form of sequences of spikes of different heights are constructed asymptotically in the limit of small activator diffusivity ε. For a pattern with k spikes, the construction yields k1 spikes that have a common small amplitude and k2 = k− k1 spikes that have a common large amplitude. A k- spike asymmetric equilibrium solution is obtained from an arbitrary ordering of the small and large spikes on the domain. It is shown that such solutions exist when the inhibitor diffusivity D is less than some critical value Dm that depends upon k1, on k2, and on other parameters associated with the Gierer–Meinhardt model. It is also shown that these asymmetric k-spike solutions bifurcate from the symmetric solution branch sk, for which k spikes have equal height. These asymmetric solutions provide connections between the branch sk and the other symmetric branches sj , for j = 1,…, k− 1. The stability of the asymmetric k-spike patterns with respect to the large O(1) eigenvalues and the small O(ε2) eigenvalues is also analyzed. It is found that the asymmetric patterns are stable with respect to the large O(1) eigenvalues when D > De, where De depends on k1 and k2, on certain parameters in the model, and on the specific ordering of the small and large spikes within a given k-spike sequence. Numerical values for De are obtained from numerical solutions of a matrix eigenvalue problem. Another matrix eigenvalue problem that determines the small eigenvalues is derived. For the examples considered, it is shown that the bifurcating asymmetric branches are all unstable with respect to these small eigenvalues.

2014 ◽  
Vol 2014 ◽  
pp. 1-15 ◽  
Author(s):  
Maryam Sarboland ◽  
Azim Aminataei

The nonlinear Burgers’ equation is a simple form of Navier-Stocks equation. The nonlinear nature of Burgers’ equation has been exploited as a useful prototype differential equation for modeling many phenomena. This paper proposes two meshfree methods for solving the one-dimensional nonlinear nonhomogeneous Burgers’ equation. These methods are based on the multiquadric (MQ) quasi-interpolation operatorℒ𝒲2and direct and indirect radial basis function networks (RBFNs) schemes. In the present schemes, the Taylors series expansion is used to discretize the temporal derivative and the quasi-interpolation is used to approximate the solution function and its spatial derivatives. In order to show the efficiency of the present methods, several experiments are considered. Our numerical solutions are compared with the analytical solutions as well as the results of other numerical schemes. Furthermore, the stability analysis of the methods is surveyed. It can be easily seen that the proposed methods are efficient, robust, and reliable for solving Burgers’ equation.


2021 ◽  
Vol 83 (3) ◽  
Author(s):  
Ginger Egberts ◽  
Fred Vermolen ◽  
Paul van Zuijlen

AbstractTo deal with permanent deformations and residual stresses, we consider a morphoelastic model for the scar formation as the result of wound healing after a skin trauma. Next to the mechanical components such as strain and displacements, the model accounts for biological constituents such as the concentration of signaling molecules, the cellular densities of fibroblasts and myofibroblasts, and the density of collagen. Here we present stability constraints for the one-dimensional counterpart of this morphoelastic model, for both the continuous and (semi-) discrete problem. We show that the truncation error between these eigenvalues associated with the continuous and semi-discrete problem is of order $${{\mathcal {O}}}(h^2)$$ O ( h 2 ) . Next we perform numerical validation to these constraints and provide a biological interpretation of the (in)stability. For the mechanical part of the model, the results show the components reach equilibria in a (non) monotonic way, depending on the value of the viscosity. The results show that the parameters of the chemical part of the model need to meet the stability constraint, depending on the decay rate of the signaling molecules, to avoid unrealistic results.


2021 ◽  
Vol 22 (4) ◽  
pp. 2030
Author(s):  
Hela Ferjani ◽  
Hammouda Chebbi ◽  
Mohammed Fettouhi

The new organic–inorganic compound (C6H9N2)2BiCl5 (I) has been grown by the solvent evaporation method. The one-dimensional (1D) structure of the allylimidazolium chlorobismuthate (I) has been determined by single crystal X-ray diffraction. It crystallizes in the centrosymmetric space group C2/c and consists of 1-allylimidazolium cations and (1D) chains of the anion BiCl52−, built up of corner-sharing [BiCl63−] octahedra which are interconnected by means of hydrogen bonding contacts N/C–H⋯Cl. The intermolecular interactions were quantified using Hirshfeld surface analysis and the enrichment ratio established that the most important role in the stability of the crystal structure was provided by hydrogen bonding and H···H interactions. The highest value of E was calculated for the contact N⋯C (6.87) followed by C⋯C (2.85) and Bi⋯Cl (2.43). These contacts were favored and made the main contribution to the crystal packing. The vibrational modes were identified and assigned by infrared and Raman spectroscopy. The optical band gap (Eg = 3.26 eV) was calculated from the diffuse reflectance spectrum and showed that we can consider the material as a semiconductor. The density functional theory (DFT) has been used to determine the calculated gap, which was about 3.73 eV, and to explain the electronic structure of the title compound, its optical properties, and the stability of the organic part by the calculation of HOMO and LUMO energy and the Fukui indices.


2003 ◽  
Vol 14 (08) ◽  
pp. 1087-1105 ◽  
Author(s):  
ZHONGCHENG WANG ◽  
YONGMING DAI

A new twelfth-order four-step formula containing fourth derivatives for the numerical integration of the one-dimensional Schrödinger equation has been developed. It was found that by adding multi-derivative terms, the stability of a linear multi-step method can be improved and the interval of periodicity of this new method is larger than that of the Numerov's method. The numerical test shows that the new method is superior to the previous lower orders in both accuracy and efficiency and it is specially applied to the problem when an increasing accuracy is requested.


2010 ◽  
Vol 24 (02) ◽  
pp. 183-193
Author(s):  
HAI-YONG DING ◽  
HONG-XIANG YANG ◽  
YE-PENG SUN ◽  
LI-LI ZHU

By considering a new four-by-four matrix eigenvalue problem, a hierarchy of Lax integrable evolution equations with four potentials is derived. The Hamiltonian structures of the resulting hierarchy are established by means of the generalized trace identity. The Liouville integrability for the hierarchy of the resulting Hamiltonian equations is presented.


1999 ◽  
Author(s):  
Alexander V. Kasharin ◽  
Jens O. M. Karlsson

Abstract The process of diffusion-limited cell dehydration is modeled for a planar system by writing the one-dimensional diffusion-equation for a cell with moving, semipermeable boundaries. For the simplifying case of isothermal dehydration with constant diffusivity, an approximate analytical solution is obtained by linearizing the governing partial differential equations. The general problem must be solved numerically. The Forward Time Center Space (FTCS) and Crank-Nicholson differencing schemes are implemented, and evaluated by comparison with the analytical solution. Putative stability criteria for the two algorithms are proposed based on numerical experiments, and the Crank-Nicholson method is shown to be accurate for a mesh with as few as six nodes.


Author(s):  
Nariman Ashrafi

The nonlinear stability and bifurcation of the one-dimensional channel (Poiseuille) flow is examined for a Johnson-Segalman fluid. The velocity and stress are represented by orthonormal functions in the transverse direction to the flow. The flow field is obtained from the conservation and constitutive equations using the Galerkin projection method. Both inertia and normal stress effects are included. The stability picture is dramatically influenced by the viscosity ratio. The range of shear rate or Weissenberg number for which the base flow is unstable increases from zero as the fluid deviates from the Newtonian limit as decreases. Typically, two turning points are observed near the critical Weissenberg numbers. The transient response is heavily influenced by the level of inertia. It is found that the flow responds oscillatorily. When the Reynolds number is small, and monotonically at large Reynolds number when elastic effects are dominated by inertia.


SIAM Review ◽  
1980 ◽  
Vol 22 (1) ◽  
pp. 99-100
Author(s):  
T. Sekiguchi ◽  
N. Kimura

Sign in / Sign up

Export Citation Format

Share Document