# Numerical simulation of reaction fronts in dissipative media

Bhanot, Rajeev Prakash (2017) Numerical simulation of reaction fronts in dissipative media. [Thesis (PhD/Research)]  Preview
Text (Whole Thesis)
Bhanot_2017.pdf

## Abstract

Fronts of reaction in certain systems (such as so-called solid flames and detonation fronts) can be simulated by a single-equation phenomenological model of Strunin (1999, 2009). This is a high-order nonlinear partial differential equation describing the shape of the front as a function of spatial coordinates and time. The equation is of active-dissipative type, with 6th-order spatial derivative. For one-dimensional case, the equation was previously solved using the Galerkin method, but only one numerical experiment with limited information on the dynamics was obtained. For two-dimensional case only two numerical ex-
periments were reported so far, in which a low-accuracy infinite difference scheme was used. In this thesis, we use a more recent and sophisticated method, namely the one-dimensional integrated radial basis function networks (1D-IRBFN). The method had been developed by Tran-Cong and May-Duy (2001, 2003) and successfully applied to several problems such as structural analysis, viscoelastic flows and fluid-structure interaction. In contrast to commonly used approaches, where a function of interest is differentiated to give approximate derivatives, leading to a reduction in convergence rate for derivatives (and this reduction increases with derivative order, which magnifes errors), the 1D-IRBFN method uses the integral formulation. It utilizes spectral approximants to represent highest-order derivatives under consideration. They are then integrated analytically to yield approximate expressions for lower-order derivatives and the function itself.

In this thesis the following main results are obtained. A numerical program implementing the 1D-IRBFN method is developed in Matlab to solve the equation of interest. The program is tested by (a) constructing a forced version of the equation, which allows analytical solution, and verifying the numerical solution against the analytical solution; (b) reproducing one-dimensional spinning waves obtained from the model previously. A modified version of the program is successfully applied to similar high-order equations modelling auto-pulses in fluid flows with elastic walls.

We obtained numerically and analyzed a far richer variety of one-dimensional dynamics of the reaction fronts. Two kinds of boundary conditions were used: homogeneous conditions on the edges of the domain, and periodic conditions corresponding to periodicity of the front on a cylinder. The dependence of the dynamics on the size of the domain is explored showing how larger space accommodates multiple spinning waves. We determined the critical domain size (bifurcation point) at which non-trivial settled regimes become possible. We found a regime where the front is shaped as a pair of kinks separated by a relatively short distance. Interestingly, the pair moves in a stable joint formation far from the boundaries. A similar regime for three connected kinks is obtained. We demonstrated that the initial condition determines the direction of motion
of the kinks, but not their size and velocity. This is typical for active-dissipative systems. The settled character of these regimes is demonstrated. We also applied the 1D-IRBFN method to two-dimensional topology corresponding to a solid cylinder. Stable spinning wave solutions are obtained for this case. Statistics for this ePrint Item
Item Type: Thesis (PhD/Research) Live Archive Doctor of Philosophy (PhD) thesis. Historic - Faculty of Health, Engineering and Sciences - School of Agricultural, Computational and Environmental Sciences (1 Jul 2013 - 5 Sep 2019) Historic - Faculty of Health, Engineering and Sciences - School of Agricultural, Computational and Environmental Sciences (1 Jul 2013 - 5 Sep 2019) Strunin, Dmitry; Tran-Cong, Thanh; Ngo-Cong, Duc 31 Jul 2017 04:32 21 Apr 2021 01:29 active dissipative systems; reaction-diffusion systems; nonlinear partial differential equation; nonlinear excitation; finite difference; spinning waves; 1D-RBFN (one dimentional radial basis function network) method 08 Information and Computing Sciences > 0899 Other Information and Computing Sciences > 089999 Information and Computing Sciences not elsewhere classified 46 INFORMATION AND COMPUTING SCIENCES > 4699 Other information and computing sciences > 469999 Other information and computing sciences not elsewhere classified doi:10.26192/5c0096f1c4671 http://eprints.usq.edu.au/id/eprint/32889 Archive Repository Staff Only