Domain decomposition method and numerical analysis of a fluid dynamics problem

A two-dimensional problem obtained by time discretization and linearization of a viscous flow governed by the incompressible Navier-Stokes equations is considered. The original domain is divided into subdomains such that their interface is a smooth (nonclosed, self-avoiding) curve with the ends belonging to the boundary of the domain. A nonconforming finite element method is constructed for the problem, and the convergence rate of the discrete solution of the problem to the exact one is estimated in the L ₂(Ω _h ) norm.     

A numerical method for the dynamics and stability of spiral waves

In this paper we propose a new method of investigating the change of dynamics in reaction-diffusion equations, which is based on approximating the Euclidian norm of state variables along with the introduction of phase space. Our method is simple in implementation and can be applied to study the dynamics of multiple spirals. The method is extended to study the stability of spiral waves by developing an algorithm which is applied to circular and meandering motions.     

A numerical method for determining the localized initial condition in the FitzHugh-Nagumo and Aliev-Panfilov models

The inverse problem for the FitzHugh-Nagumo and Aliev-Panfilov models describing wave propagation in excitable media is considered. The problem lies in determining a localized initial condition from measurements on the external boundary of a plane region. A numerical method for solving the inverse problem is proposed, and the results from a numerical solution of the inverse problem for regions similar to different sections of a heart are presented.     

Convergence and stability of a numerical method for micromagnetics

The convergence and stability of a numerical method, which applies a nonconforming finite element method and an artificial boundary method to a multi-atomic Young measure relaxation model, for micromagnetics are analyzed. By revealing some key properties of the solution sets of both the continuous and discrete problems, we show that our numerical method is stable, and the solution set of the continuous problem is well approximated by those of the discrete problems. The performance of our method is also illustrated by some numerical examples. The research was supported in part by the Major State Basic Research Projects (2005CB321701), NSFC projects (10431050, 10571006, 10528102 and 10871011) and RFDP of China.     

Regularized numerical integration of multibody dynamics with the generalized method

This paper discusses the consistent regularization property of the generalized α method when applied as an integrator to an initial value high index and singular differential-algebraic equation model of a multibody system. The regularization comes from within the discretization itself and the discretization remains consistent over the range of values the regularization parameter may take. The regularization involves increase of the smallest singular values of the ill-conditioned Jacobian of the discretization and is different from Baumgarte and similar techniques which tend to be inconsistent for poor choice of regularization parameter. This regularization also helps where pre-conditioning the Jacobian by scaling is of limited effect, for example, when the scleronomic constraints contain multiple closed loops or singular configuration or when high index path constraints are present. The feed-forward control in Kane’s equation models is additionally considered in the numerical examples to illustrate the effect of regularization. The discretization presented in this work is adopted to the first order DAE system (unlike the original method which is intended for second order systems) for its A-stability and same order of accuracy for positions and velocities.     

A numerical method for a blood-artery interaction problem

We present a numerical method for the resolution of a bidimensional blood flow problem and more generally for a fluid flow surrounded by a time dependent domain. Our approach is based on an ALE formulation which is solved using a Galerkin method with an eigenvectors basis set on the initial fixed domain.     

A stochastic approach for simulating spatially inhomogeneous coagulation dynamics in the gelation regime

We present a stochastic approach for the simulation of coagulation–diffusion dynamics in the gelation regime. The method couples the mass flow algorithm for coagulation processes with a stochastic variant of the diffusion-velocity method in a discretized framework. The simulation of the stochastic processes occurs according to an optimized implementation of the principle of grouping the possible events. A full simulation of a particle system driven by coagulation–diffusion dynamics is performed with a high degree of accuracy. This allows a qualitative and quantitative analysis of the behaviour of the system. The performance of the method becomes more evident especially in the gelation regime, where the computations become usually very time consuming.     

A second-order accurate numerical method for a fractional wave equation

We study a generalized Crank–Nicolson scheme for the time discretization of a fractional wave equation, in combination with a space discretization by linear finite elements. The scheme uses a non-uniform grid in time to compensate for the singular behaviour of the exact solution at t = 0. With appropriate assumptions on the data and assuming that the spatial domain is convex or smooth, we show that the error is of order k ² + h ², where k and h are the parameters for the time and space meshes, respectively.     

Kufarev's method for determining the Schwarz-Christoffel parameters

Summary We present an analysis of a method suggested by Kufarev for the numerical determination of the parameters involved in the Schwarz-Christoffel integral. We test the approach for certain polygonal domains including two with re-entrant corners and one with a cut. Finally, we consider the Motz problem from this point of view.     

On the convergence of a numerical method for optimal control problems

Under suitable restrictions, convergence to the solutions of a class of optimal control problems is proved for a method by Rosen in which the differential equations, constraints, and cost functional are discretized, and the resulting mathematical programming problem is solved approximately by a penalty-function approach.     

A new direct method for the numerical calculation of a swallowtail

A nonlinear operator equation, depending on three real parameters, is given. We designed a new direct method for the calculation of a swallowtail of this equation. It turns out that the special structure of the extended system can be used to construct an efficient algorithm.     

A parameter-uniform numerical method for a Sobolev problem with initial layer

The present study is concerned with the numerical solution, using finite difference method of a one-dimensional initial-boundary value problem for a linear Sobolev or pseudo-parabolic equation with initial jump. In order to obtain an efficient method, to provide good approximations with independence of the perturbation parameter, we have developed a numerical method which combines a finite difference spatial discretization on uniform mesh and the implicit rule on Shishkin mesh(S-mesh) for the time variable. The fully discrete scheme is shown to be convergent of order two in space and of order one expect for a logarithmic factor in time, uniformly in the singular perturbation parameter. Some numerical results confirming the expected behavior of the method are shown.