Cyclic reduction and FACR methods for piecewise hermite bicubic orthogonal spline collocation

Cyclic reduction and Fourier analysis-cyclic reduction (FACR) methods are presented for the solution of the linear systems which arise when orthogonal spline collocation with piecewise Hermite bicubics is applied to boundary value problems for certain separable partial differential equations on a rectangle. On anN×N uniform partition, the cyclic reduction and Fourier analysis-cyclic reduction methods requireO(N ²log₂ N) andO(N ²log₂log₂ N) arithmetic operations, respectively.     

Riabouchinsky flow with partially penetrable obstacle

An explicitly solvable Riabouchinsky model with a partially penetrable obstacle is introduced. This model applied to the estimation of the efficiency of free flow turbines allows us to take into account the pressure drop past the lamina.     

Trefftz,collocation, and other boundary methods—A comparison

In this article we survey the Trefftz method (TM), the collocation method (CM), and the collocation Trefftz method (CTM). We also review the coupling techniques for the interzonal conditions, which include the indirect Trefftz method, the original Trefftz method, the penalty plus hybrid Trefftz method, and the direct Trefftz method. Other boundary methods are also briefly described. Key issues in these algorithms, including the error analysis, are addressed. New numerical results are reported. Comparisons among TMs and other numerical methods are made. It is concluded that the CTM is the simplest algorithm and provides the most accurate solution with the best numerical stability. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

A nodal coupling of finite and boundary elements

This article presents and analyzes a simple method for the exterior Laplace equation through the coupling of finite and boundary element methods. The main novelty is the use of a smooth parametric artificial boundary where boundary elements fit without effort together with a straight approximate triangulation in the bounded area, with the coupling done only in nodes. A numerically integrated version of the algorithm is also analyzed. Finally, an isoparametric variant with higher order is proposed. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 555–570, 2003     

Molecular dynamics simulation of entropy and surface tension for grain boundary of α-Fe

The grain boundary is an interface and the surface tension is one of its important thermodynamic properties. In this paper, the surface tension of the ∑9 grain boundary for α-Fe at various temperatures and pressures is calculated by means of Computer Molecular Dynamics (CMD). The results agree satisfactorily with the experimental data. It is shown that the contribution of entropy to surface tension of grain boundary can be ignored. The project supported by the National Natural Science Foundation of China and the Science Foundation of Chinese Academy of Sciences.     

Anisotropic elastodynamics in a half space: An analytic method for polynomial data

An initial boundary value problem for the dynamic system of anisotropic elasticity in a half space is studied in the paper. A novel method of finding an exact solution of this problem for a special polynomial form of initial data and inhomogeneous term of the system is described. On the base of this method the simulation of elastic waves in different homogeneous anisotropic half spaces is implemented.     

HIGH ACCURACY NUMERICAL METHOD OF THIN-FILM PROBLEMS IN MICROMAGNETICS

In this paper, a new high accuracy numerical method for the thin-film problems of micron and submicron size ferromagnetic elements is proposed. For the computation of stray field, we use the finite element method(FEM) by introducing a semi-discrete artificial boundary condition [1, 2]. In our numerical experiments about the domain patterns and their movement, we can see that the results are accordant to that of experiments and other numerical methods. Our method are very convenient to deal with arbitrary shape of thin films such as a polygon with high accuracy.     

On the properties of solutions to the Goursat-Darboux problem with boundary and distributed controls

We consider a control system described by the Goursat-Darboux equation. The system is controlled by distributed and boundary controls. The controls are subject to the constraints given as multivalued mappings with closed, possibly nonconvex, values depending on the phase variable. Alongside the initial constraints, we consider the convexified constraints and the constraints whose values are the extreme points of the convexified constraints. We study the questions of existence of solutions and establish connections between the solutions under various constraints.     

Global existence and blow-up of solutions to a parabolic system with nonlocal sources and boundaries

This paper deals with a semi-linear parabolic system with nonlinear nonlocal sources and nonlocal boundaries.By using super-and sub-solution techniques,we first give the sufficient conditions that the classical solution exists globally and blows up in a finite time respectively,and then give the necessary and sufficient conditions that two components u and v blow up simultaneously.Finally,the uniform blow-up profiles in the interior are presented.     

应变Si1-xGex层中p型杂质电离常温和低温特性

本文采用解析的方法计算了应变Si1-xGex层中p型杂质电离度与Ge组分x、温度T以及掺杂浓度N的关系.发现常温时,在同一Ge组分下,随着掺杂浓度的升高,杂质的电离度的先变小,而后又迅速上升到1.在同一掺杂浓度下,轻掺杂时,杂质的电离度随Ge组分的增加先变大,而后几乎不变;重掺杂时,杂质电离能变为0后,杂质电离度为1.低温下,轻掺杂时,载流子低温冻析效应较为明显,杂质的电离度普遍较小,当掺杂浓度大于Mott转换点时,载流子冻析效应不再明显,电离率迅速上升到1.在同一Ge组分下,随着掺杂浓度的升高,杂质的电离度先变小,后变大,而后又迅速上升到1.在同一掺杂浓度下,轻掺杂时,杂质的电离度随Ge组分的增加变大;重掺杂时,杂质电离能变为0后,杂质电离度为1.