Highly efficient NURBS-based isogeometric analysis for coupled nonlinear diffusion-reaction equations with and without advection

Ilham Asmouh,Alexander Ostermann

2024-04-05

Abstract:Nonlinear diffusion-reaction systems model a multitude of physical phenomena. A common situation is biological development modeling where such systems have been widely used to study spatiotemporal phenomena in cell biology. Systems of coupled diffusion-reaction equations are usually subject to some complicated features directly related to their multiphysics nature. Moreover, the presence of advection is source of numerical instabilities, in general, and adds another challenge to these systems. In this study, we propose a NURBS-based isogeometric analysis (IgA) combined with a second-order Strang operator splitting to deal with the multiphysics nature of the problem. The advection part is treated in a semi-Lagrangian framework and the resulting diffusion-reaction equations are then solved using an efficient time-stepping algorithm based on operator splitting. The accuracy of the method is studied by means of a advection-diffusion-reaction system with analytical solution. To further examine the performance of the new method on complex geometries, the well-known Schnakenberg-Turing problem is considered with and without advection. Finally, a Gray-Scott system on a circular domain is also presented. The results obtained demonstrate the efficiency of our new algorithm to accurately reproduce the solution in the presence of complex patterns on complex geometries. Moreover, the new method clarifies the effect of geometry on Turing patterns.

Numerical Analysis

What problem does this paper attempt to address?

This paper focuses on the efficient solution of coupled nonlinear diffusion-reaction equation systems, which are widely used in physical phenomena such as biological development modeling. The paper proposes a geometric integration analysis method based on Non-Uniform Rational B-Splines (NURBS), combined with second-order Strang operator splitting technique to handle the multi-physical properties of the problem. In the case of advection, this method can effectively handle numerical instability and solve complex geometric shape problems.

Specifically, the equations involved in the study are a set of nonlinear diffusion-reaction systems, which describe the change of concentrations of activators and inhibitors over time and space. In the paper, the authors use a semi-Lagrangian framework to deal with the transport part, and then solve the diffusion-reaction equation using a time-stepping algorithm based on operator splitting. The method they used does not require linearization, can accurately reproduce solutions in complex patterns, and reveals the influence of geometric shape on Turing patterns.

In addition, the paper also discusses various numerical methods used in past research, such as finite volume methods, finite element methods, and spectral element methods, but these methods are mostly limited to simple geometric shapes or one-dimensional problems. The advantage of the new method is its ability to handle complex geometric shapes and maintain stability and efficiency when considering transport effects.

Finally, the paper verifies the accuracy of the new method through several test problems, including transport-diffusion-reaction systems with analytical solutions, Schnakenberg-Turing model, and Gray-Scott model. The results show that the new algorithm can accurately simulate solutions on complex geometric shapes and analyze the effects of transport on Turing pattern formation.

In conclusion, this paper aims to improve the accuracy and efficiency of simulating coupled nonlinear diffusion-reaction systems in complex environments through an innovative numerical method.