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.     

A mesh-free numerical method for solution of the family of Kuramoto-Sivashinsky equations

In this paper, radial basis function (RBFs) based mesh-free method is implemented to find numerical solution of the Kuramoto-Sivashinsky equations. This approach has an edge over traditional methods such as finite-difference and finite element methods because it does not require a mesh to discretize the problem domain, and a set of scattered nodes in the domain of influence provided by initial data is required for the realization of the method. The accuracy of the method is assessed in terms of the error norms L₂,L_∞, number of nodes in the domain of influence, free parameter, dependent parameter RBFs and time step length. Numerical experiments demonstrate accuracy and robustness of the method for solving a class of nonlinear partial differential equations.     

A random Euler scheme for Carathéodory differential equations

We study a random Euler scheme for the approximation of Carathéodory differential equations and give a precise error analysis. In particular, we show that under weak assumptions, this approximation scheme obtains the same rate of convergence as the classical Monte–Carlo method for integration problems.