Numerical solutions of a model for the propagation of a surface-catalysed flame in a tube

Numerical simulations of a surface-catalysed flame in a tubeare performed, corresponding to an experiment where a premixedfuel is fed into a tube whose inner surface is coated with acatalyst. In these experiments, subsequent to ignition, a reactionwave can be seen as a red-hot region which propagates back alongthe tube towards the inlet, and is due to low temperature combustionoccurring only on the inner surface of the tube where the catalystis present. The solutions of a mathematical model for this behaviourshow that initial-value problems do indeed result in such steadilypropagating waves. The numerically obtained wave speeds andsteady solution are compared to a previous large Damköhlernumber (D_a) asymptotic analysis using a simple reaction ratemodel, and agreement is very good even for moderately largevalues of D_a. However, for such Damköhler numbers, thewave speeds are found to be much larger than observed experimentally.Indeed, the simulations show that O(1) values of D_a are requiredto obtain the lower experimental wave speeds. Nevertheless,the wave speeds as a function of flow rate through the tubedo not agree well with the preliminary experimental resultsfor any choice of the parameters. A more realistic, Arrheniusreaction rate model is then considered. The Arrhenius modelpredicts a rapid change in temperature at the wave front, inmuch better agreement with the experiments than for the simplerreaction model.     

Numerical treatment of retarded boundary integral equations by sparse panel clustering

Stefan Sauter Institute for Mathematics, University of Zürich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland ^{We consider the wave equation in a boundary integral formulation.The discretization in time is done by using convolution quadraturetechniques and a Galerkin boundary element method for the spatialdiscretization. In a previous paper, we have introduced a sparseapproximation of the system matrix by cut-off, in order to reducethe storage costs. In this paper, we extend this approach byintroducing a panel clustering method to further reduce thesecosts.}     

Comparison of ODE methods for laminar reacting gas flow simulations

Two‐dimensional transient simulations are presented of the transport phenomena and multispecies, multireaction chemistry in chemical vapor deposition (CVD). The transient simulations are run until steady state, such that the steady state can be validated against the steady state solutions from literature. We compare various time integration methods in terms of efficiency and robustness. Besides stability, which is important due to the stiffness of the problem, preservation of non‐negativity is crucial. It appears that this latter condition on a time integration method is much more restrictive toward the time step size than stability. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

A multidomain integrated‐radial‐basis‐function collocation method for elliptic problems

This article is concerned with the use of integrated radial‐basis‐function networks (IRBFNs) and nonoverlapping domain decompositions (DDs) for numerically solving one‐ and two‐dimensional elliptic problems. A substructuring technique is adopted, where subproblems are discretized by means of one‐dimensional IRBFNs. A distinguishing feature of the present DD technique is that the continuity of the RBF solution across the interfaces is enforced with one order higher than with conventional DD techniques. Several test problems governed by second‐ and fourth‐order differential equations are considered to investigate the accuracy of the proposed technique. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

CONVERGENCE OF AN IMMERSED INTERFACE UPWIND SCHEME FOR LINEAR ADVECTION EQUATIONS WITH PIECEWISE CONSTANT COEFFICIENTS I： L^1-ERROR ESTIMATES

We study the L^l-error estimates for the upwind scheme to the linear advection equations with a piecewise constant coefficients modeling linear waves crossing interfaces. Here the interface condition is immersed into the upwind scheme. We prove that, for initial data with a bounded variation, the numerical solution of the immersed interface upwind scheme converges in L^l-norm to the differential equation with the corresponding interface condition. We derive the one-halfth order L^l-error bounds with explicit coefficients following a technique used in [25]. We also use some inequalities on binomial coefficients proved in a consecutive paper [32].     

Type-II hidden symmetries through weak symmetries for nonlinear partial differential equations

The Type-II hidden symmetries are extra symmetries in addition to the inherited symmetries of the differential equations when the number of independent and dependent variables is reduced by a Lie point symmetry. In [B. Abraham-Shrauner, K.S. Govinder, Provenance of Type II hidden symmetries from nonlinear partial differential equations, J. Nonlinear Math. Phys. 13 (2006) 612-622] Abraham-Shrauner and Govinder have analyzed the provenance of this kind of symmetries and they developed two methods for determining the source of these hidden symmetries. The Lie point symmetries of a model equation and the two-dimensional Burgers' equation and their descendants were used to identify the hidden symmetries. In this paper we analyze the connection between one of their methods and the weak symmetries of the partial differential equation in order to determine the source of these hidden symmetries. We have considered the same models presented in [B. Abraham-Shrauner, K.S. Govinder, Provenance of Type II hidden symmetries from nonlinear partial differential equations, J. Nonlinear Math. Phys. 13 (2006) 612-622], as well as the WDVV equations of associativity in two-dimensional topological field theory which reduces, in the case of three fields, to a single third order equation of Monge-Ampère type. We have also studied a second order linear partial differential equation in which the number of independent variables cannot be reduced by using Lie symmetries, however when is reduced by using nonclassical symmetries the reduced partial differential equation gains Lie symmetries.     

Asymptotic behaviour of the three-dimensional α-Navier-Stokes model with delays

We prove the existence and exponential stability of the stationary solutions for a three-dimensional α-Navier-Stokes model with delays. Instead of working directly with the model, we establish the results for a general abstract delay model, and then we apply them to our particular situation. We use two techniques in the stability analysis: the Lyapunov and Razumikhin approaches.     

Shape reconstruction of an inverse Stokes problem

This paper deals with the shape reconstruction of a viscous incompressible fluid driven by the Stokes flow. For the approximate solution of the ill-posed and nonlinear problem we propose a regularized Newton method. A theoretical foundation for the Newton method is given by establishing the differentiability of the initial boundary value problem with respect to the interior boundary curve in the sense of a domain derivative. The numerical examples show that our theory is useful for practical purpose and the proposed algorithm is feasible.     

Lattice-Boltzmann type relaxation systems and high order relaxation schemes for the incompressible Navier-Stokes equations

A relaxation system based on a Lattice-Boltzmann type discrete velocity model is considered in the low Mach number limit. A third order relaxation scheme is developed working uniformly for all ranges of the mean free path and Mach number. In the incompressible Navier-Stokes limit the scheme reduces to an explicit high order finite difference scheme for the incompressible Navier-Stokes equations based on nonoscillatory upwind discretization. Numerical results and comparisons with other approaches are presented for several test cases in one and two space dimensions.