skip to main content

Title: Front propagation and global bifurcations in a multivariable reaction-diffusion model
We study the existence and stability of propagating fronts in Meinhardt’s multivariable reaction-diffusion model of branching in one spatial dimension. We identify a saddle-node-infinite-period bifurcation of fronts that leads to episodic front propagation in the parameter region below propagation failure and show that this state is stable. Stable constant speed fronts exist only above this parameter value. We use numerical continuation to show that propagation failure is a consequence of the presence of a T-point corresponding to the formation of a heteroclinic cycle in a spatial dynamics description. Additional T-points are identified that are responsible for a large multiplicity of different unstable traveling front-peak states. The results indicate that multivariable models may support new types of behavior that are absent from typical two-variable models but may nevertheless be important in developmental processes such as branching and somitogenesis.  more » « less
Award ID(s):
Author(s) / Creator(s):
Date Published:
Journal Name:
Chaos: An Interdisciplinary Journal of Nonlinear Science
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We study the linear stability properties of spatially localized single- and multi-peak states generated in a subcritical Turing bifurcation in the Meinhardt model of branching. In one spatial dimension, these states are organized in a foliated snaking structure owing to peak–peak repulsion but are shown to be all linearly unstable, with the number of unstable modes increasing with the number of peaks present. Despite this, in two spatial dimensions, direct numerical simulations reveal the presence of stable single- and multi-spot states whose properties depend on the repulsion from nearby spots as well as the shape of the domain and the boundary conditions imposed thereon. Front propagation is shown to trigger the growth of new spots while destabilizing others. The results indicate that multi-variable models may support new types of behavior that are absent from typical two-variable models. 
    more » « less
  2. Abstract

    Wave front propagation with nontrivial bottom topography is studied within the formalism of hyperbolic long wave models. Evolution of nonsmooth initial data is examined, and, in particular, the splitting of singular points and their short time behavior is described. In the opposite limit of longer times, the local analysis of wave fronts is used to estimate the gradient catastrophe formation and how this is influenced by the topography. The limiting cases when the free surface intersects the bottom boundary, belonging to the so‐called “physical” and “nonphysical” vacuum classes, are examined. Solutions expressed by power series in the spatial variable lead to a hierarchy of ordinary differential equations for the time‐dependent series coefficients, which are shown to reveal basic differences between the two vacuum cases: for nonphysical vacuums, the equations of the hierarchy are recursive and linear past the first two pairs, whereas for physical vacuums, the hierarchy is nonrecursive, fully coupled, and nonlinear. The former case may admit solutions that are free of singularities for nonzero time intervals, whereas the latter is shown to develop nonstandard velocity shocks instantaneously. Polynomial bottom topographies simplify the hierarchy, as they contribute only a finite number of inhomogeneous forcing terms to the equations in the recursion relations. However, we show that truncation to finite‐dimensional systems and polynomial solutions is in general only possible for the case of a quadratic bottom profile. In this case, the system's evolution can reduce to, and is completely described by, a low‐dimensional dynamical system for the time‐dependent coefficients. This system encapsulates all the nonlinear properties of the solution for general power series initial data, and, in particular, governs the loss of regularity in finite times at the dry point. For the special case of parabolic bottom topographies, an exact, self‐similar solution class is introduced and studied to illustrate via closed‐form expressions the general results.

    more » « less
  3. Objective: Mortality-trends from alcoholic liver disease (ALD) have recently increased and they differ by various factors in the U.S. However, these trends have only been analyzed using univariate models and in reality they may be influenced by various factors. We thus examined trends in age-standardized mortality from ALD among U.S. adults for 1999-2017, using multivariable piecewise log-linear models. Methods: We collected mortality-data from the Centers for Disease Control and Prevention Wide-ranging Online Data for Epidemiologic Research database, using the Underlying Cause of Death. Results: We identified 296,194 deaths from ALD and 346,386 deaths indirectly attributable to ALD during the period from 1999-2017. The multivariable-adjusted, age-standardized ALD mortality was stable during 1999-2006 (annual percentage change [APC]=-2.24, P=0.24), and increased during 2006-2017 (APC=3.18, P<0.006). Their trends did not differ by sex, race, age or urbanization. Subgroup analyses revealed upward multivariable-adjusted, age-standardized mortality-trends in alcoholic fatty liver (APC=4.64, P<0.001), alcoholic hepatitis (APC=4.38, P<0.001), and alcoholic cirrhosis (APC=5.33, P<0.001), but downward mortality-trends in alcoholic hepatic failure (APC=-1.63, P=0.006) and unspecified ALD (APC=-0.86, P=0.013). Strikingly, non-alcoholic cirrhosis also had an upward multivariable-adjusted, age-standardized mortality-trend (APC=0.69, P=0.046). By contrast, recent mortality-trends were stable for all cause of deaths (APC=-0.39, P=0.379) and downward for malignant neoplasms excluding liver cancer (APC=-2.82, P<0.001), infections (APC=-2.60, P<0.001), cardiovascular disease (APC=-0.69, P=0.044) and respiratory disease (APC=-0.56, P=0.002). The adjusted mortality with ALD as a contributing cause of death also had an upward trend during 2000-2017 (APC=5.47, P<0.001). Strikingly, common comorbidities of ALD, including hepatocellular carcinoma, cerebrovascular and ischemic heart cardiovascular diseases and sepsis, had upward trends during the past 14 to 16 years. Conclusions: ALD had an upward multivariable-adjusted, age-standardized mortality-trend among U.S. adults, without significant differences by sex, race, age or urbanization. Three ALD subtypes (alcoholic fatty liver, alcoholic hepatitis and alcoholic cirrhosis) and non-alcoholic cirrhosis had upward morality-trends, while other ALD subtypes and other causes of death did not. 
    more » « less
  4. Colmcille Caulfield (Ed.)

    When a saturated brine layer is cooled from above, both a convective temperature front as well as a front of sedimenting salt crystals can form. We employ direct numerical simulations to investigate the evolution and interaction of these two density fronts. Depending on the ratio of the temperature front velocity and the crystal settling velocity, which is governed by a dimensionless parameter in the form of a Rayleigh number, we find that either two separate fronts exist for all times, two initially separate fronts combine into a single front after some time or a single front exists at all times. We furthermore propose approximate scaling laws for the propagation of the thermal and crystal fronts in each regime and compare them with the simulation data, with generally good agreement.

    more » « less
  5. Abstract Becker et al. measured the mean free path of Lyman-limit photons in the intergalactic medium (IGM) at z = 6. The short value suggests that absorptions may have played a prominent role in reionization. Here we study physical properties of ionizing photon sinks in the wake of ionization fronts (I-fronts) using radiative hydrodynamic simulations. We quantify the contributions of gaseous structures to the Lyman-limit opacity by tracking the column-density distributions in our simulations. Within Δ t = 10 Myr of I-front passage, we find that self-shielding systems ( N H I > 10 17.2 cm −2 ) are comprised of two distinct populations: (1) overdensity Δ ∼ 50 structures in photoionization equilibrium with the ionizing background, and (2) Δ ≳ 100 density peaks with fully neutral cores. The self-shielding systems contribute more than half of the opacity at these times, but the IGM evolves considerably in Δ t ∼ 100 Myr as structures are flattened by pressure smoothing and photoevaporation. By Δ t = 300 Myr, they contribute ≲10% to the opacity in an average 1 Mpc 3 patch of the universe. The percentage can be a factor of a few larger in overdense patches, where more self-shielding systems survive. We quantify the characteristic masses and sizes of self-shielding structures. Shortly after I-front passage, we find M = 10 4 –10 8 M ⊙ and effective diameters d eff = 1–20 ckpc h −1 . These scales increase as the gas relaxes. The picture herein presented may be different in dark matter models with suppressed small-scale power. 
    more » « less