Analysing diffusion and flow-driven instability using semidefinite programming

Yutaka Hori, Hiroki Miyazako

Research output: Contribution to journalArticle

Abstract

Diffusion and flow-driven instability, or transport-driven instability, is one of the central mechanisms to generate inhomogeneous gradient of concentrations in spatially distributed chemical systems. However, verifying the transport-driven instability of reaction-diffusion-advection systems requires checking the Jacobian eigenvalues of infinitely many Fourier modes, which is computationally intractable. To overcome this limitation, this paper proposes mathematical optimization algorithms that determine the stability/instability of reaction-diffusion-advection systems by finite steps of algebraic calculations. Specifically, the stability/instability analysis of Fourier modes is formulated as a sum-of-squares optimization program, which is a class of convex optimization whose solvers are widely available as software packages. The optimization program is further extended for facile computation of the destabilizing spatial modes. This extension allows for predicting and designing the shape of the concentration gradient without simulating the governing equations. The streamlined analysis process of self-organized pattern formation is demonstrated with a simple illustrative reaction model with diffusion and advection.

Original languageEnglish
Article number20180586
JournalJournal of the Royal Society Interface
Volume16
Issue number150
DOIs
Publication statusPublished - 2019 Jan 1

Fingerprint

Advection
Computer Communication Networks
Fourier Analysis
Software
Convex optimization
Software packages

Keywords

  • Convex optimization
  • Reaction-diffusion-advection model
  • Self-organized pattern formation
  • Semidefinite programming
  • Transport-driven instability

ASJC Scopus subject areas

  • Biotechnology
  • Biophysics
  • Bioengineering
  • Biomaterials
  • Biochemistry
  • Biomedical Engineering

Cite this

Analysing diffusion and flow-driven instability using semidefinite programming. / Hori, Yutaka; Miyazako, Hiroki.

In: Journal of the Royal Society Interface, Vol. 16, No. 150, 20180586, 01.01.2019.

Research output: Contribution to journalArticle

@article{c1a269ea1a4b4260bd0301b2341c9c56,
title = "Analysing diffusion and flow-driven instability using semidefinite programming",
abstract = "Diffusion and flow-driven instability, or transport-driven instability, is one of the central mechanisms to generate inhomogeneous gradient of concentrations in spatially distributed chemical systems. However, verifying the transport-driven instability of reaction-diffusion-advection systems requires checking the Jacobian eigenvalues of infinitely many Fourier modes, which is computationally intractable. To overcome this limitation, this paper proposes mathematical optimization algorithms that determine the stability/instability of reaction-diffusion-advection systems by finite steps of algebraic calculations. Specifically, the stability/instability analysis of Fourier modes is formulated as a sum-of-squares optimization program, which is a class of convex optimization whose solvers are widely available as software packages. The optimization program is further extended for facile computation of the destabilizing spatial modes. This extension allows for predicting and designing the shape of the concentration gradient without simulating the governing equations. The streamlined analysis process of self-organized pattern formation is demonstrated with a simple illustrative reaction model with diffusion and advection.",
keywords = "Convex optimization, Reaction-diffusion-advection model, Self-organized pattern formation, Semidefinite programming, Transport-driven instability",
author = "Yutaka Hori and Hiroki Miyazako",
year = "2019",
month = "1",
day = "1",
doi = "10.1098/rsif.2018.0586",
language = "English",
volume = "16",
journal = "Journal of the Royal Society Interface",
issn = "1742-5689",
publisher = "Royal Society of London",
number = "150",

}

TY - JOUR

T1 - Analysing diffusion and flow-driven instability using semidefinite programming

AU - Hori, Yutaka

AU - Miyazako, Hiroki

PY - 2019/1/1

Y1 - 2019/1/1

N2 - Diffusion and flow-driven instability, or transport-driven instability, is one of the central mechanisms to generate inhomogeneous gradient of concentrations in spatially distributed chemical systems. However, verifying the transport-driven instability of reaction-diffusion-advection systems requires checking the Jacobian eigenvalues of infinitely many Fourier modes, which is computationally intractable. To overcome this limitation, this paper proposes mathematical optimization algorithms that determine the stability/instability of reaction-diffusion-advection systems by finite steps of algebraic calculations. Specifically, the stability/instability analysis of Fourier modes is formulated as a sum-of-squares optimization program, which is a class of convex optimization whose solvers are widely available as software packages. The optimization program is further extended for facile computation of the destabilizing spatial modes. This extension allows for predicting and designing the shape of the concentration gradient without simulating the governing equations. The streamlined analysis process of self-organized pattern formation is demonstrated with a simple illustrative reaction model with diffusion and advection.

AB - Diffusion and flow-driven instability, or transport-driven instability, is one of the central mechanisms to generate inhomogeneous gradient of concentrations in spatially distributed chemical systems. However, verifying the transport-driven instability of reaction-diffusion-advection systems requires checking the Jacobian eigenvalues of infinitely many Fourier modes, which is computationally intractable. To overcome this limitation, this paper proposes mathematical optimization algorithms that determine the stability/instability of reaction-diffusion-advection systems by finite steps of algebraic calculations. Specifically, the stability/instability analysis of Fourier modes is formulated as a sum-of-squares optimization program, which is a class of convex optimization whose solvers are widely available as software packages. The optimization program is further extended for facile computation of the destabilizing spatial modes. This extension allows for predicting and designing the shape of the concentration gradient without simulating the governing equations. The streamlined analysis process of self-organized pattern formation is demonstrated with a simple illustrative reaction model with diffusion and advection.

KW - Convex optimization

KW - Reaction-diffusion-advection model

KW - Self-organized pattern formation

KW - Semidefinite programming

KW - Transport-driven instability

UR - http://www.scopus.com/inward/record.url?scp=85061330150&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85061330150&partnerID=8YFLogxK

U2 - 10.1098/rsif.2018.0586

DO - 10.1098/rsif.2018.0586

M3 - Article

VL - 16

JO - Journal of the Royal Society Interface

JF - Journal of the Royal Society Interface

SN - 1742-5689

IS - 150

M1 - 20180586

ER -