Stable advection‐reaction‐diffusion with arbitrary anisotropy

Turing first theorized that many biological patterns arise through the processes of reaction and diffusion. Subsequently, reaction‐diffusion systems have been studied in many fields, including computer graphics. We first show that for visual simulation purposes, reaction‐diffusion equations can be made unconditionally stable using a variety of straightforward methods. Second, we propose an anisotropy embedding that significantly expands the space of possible patterns that can be generated. Third, we show that by adding an advection term, the simulation can be coupled to a fluid simulation to produce visually appealing flows. Fourth, we couple fast marching methods to our anisotropy embedding to create a painting interface to the simulation. Unconditional stability is maintained throughout, and our system runs at interactive rates. Finally, we show that on the Cell processor, it is possible to implement reaction‐diffusion on top of an existing fluid solver with no significant performance impact. Copyright © 2007 John Wiley & Sons, Ltd.

[1]  Insung Ihm,et al.  Animation of reactive gaseous fluids through chemical kinetics , 2004, SCA '04.

[2]  J. E. Pearson Complex Patterns in a Simple System , 1993, Science.

[3]  James F. O'Brien,et al.  A semi-Lagrangian contouring method for fluid simulation , 2005, TOGS.

[4]  A. M. Turing,et al.  The chemical basis of morphogenesis , 1952, Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences.

[5]  David S. Ebert,et al.  Texturing and Modeling: A Procedural Approach , 1994 .

[6]  Michael Griebel,et al.  Numerical Simulation in Fluid Dynamics: A Practical Introduction , 1997 .

[7]  J A Sethian,et al.  A fast marching level set method for monotonically advancing fronts. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[8]  William V. Baxter,et al.  A viscous paint model for interactive applications , 2004, SIGGRAPH '04.

[9]  Greg Turk,et al.  Generating textures on arbitrary surfaces using reaction-diffusion , 1991, SIGGRAPH.

[10]  D. Barkley A model for fast computer simulation of waves in excitable media , 1991 .

[11]  James A. Sethian,et al.  The Fast Construction of Extension Velocities in Level Set Methods , 1999 .

[12]  J. Shewchuk An Introduction to the Conjugate Gradient Method Without the Agonizing Pain , 1994 .

[13]  S. Yoshizawa,et al.  An Active Pulse Transmission Line Simulating Nerve Axon , 1962, Proceedings of the IRE.

[14]  J. Tsitsiklis,et al.  Efficient algorithms for globally optimal trajectories , 1994, Proceedings of 1994 33rd IEEE Conference on Decision and Control.

[15]  David H. Laidlaw,et al.  Cellular texture generation , 1995, SIGGRAPH.

[16]  H. Meinhardt Models of biological pattern formation , 1982 .

[17]  L. Trefethen Finite Difference and Spectral Methods for Ordinary and Partial Differential Equations , 1996 .

[18]  J. Verwer,et al.  Numerical solution of time-dependent advection-diffusion-reaction equations , 2003 .

[19]  Andrew Witkin,et al.  Reaction-diffusion textures , 1991, SIGGRAPH.

[20]  Jos Stam,et al.  Stable fluids , 1999, SIGGRAPH.

[21]  Ronald Fedkiw,et al.  Visual simulation of smoke , 2001, SIGGRAPH.