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.

     

混合边界条件下的三维定常 旋转Navier-Stokes方程的原始变量有限元计算

本文利用原始变量有限元法求解混合边界条件下的三维定常旋转Navier-Stokes方程,证明了离散问题解的存在唯一性,得到了有限元解的最优误差估计.给出了求解原始变量有限元逼近解的简单迭代算法,并证明了算法的收敛性.针对三维情况下计算资源的限制,采用压缩的行存储格式存储刚度矩阵的非零元素,并利用不完全的LU分解作预处理的GMRES方法求解线性方程组.最后分析了简单迭代和牛顿迭代的优劣对比,数值算例表明在同样精度下简单迭代更节约计算时间.     

Some exact solutions of the ideal MHD equations through symmetry reduction method

We use the symmetry reduction method based on Lie group theory to obtain some exact solutions, the so-called invariant solutions, of the ideal magnetohydrodynamic equations in (3+1) dimensions. In particular, these equations are invariant under a Galilean-similitude Lie algebra for which the classification by conjugacy classes of r-dimensional subalgebras (1?r?4) was already known. We restrict our study to the three-dimensional Galilean-similitude subalgebras that give us systems composed of ordinary differential equations. Here, some examples of these solutions are presented with a brief physical interpretation.