New initial-value method for singularly perturbed boundary-value problems

An initial-value method is given for second-order singularly perturbed boundary-value problems with a boundary layer at one endpoint. The idea is to replace the original two-point boundary value problem by two suitable initial-value problems. The method is very easy to use and to implement. Nontrivial text problems are used to show the feasibility of the given method, its versatility, and its performance in solving linear and nonlinear singularly perturbed problems.This work was supported in part by the Consiglio Nazionale delle Ricerche, Contract No. 86.02108.01, and in part by the Ministero della Pubblica Istruzione.     

Numerical differentiation from a viewpoint of regularization theory

In this paper, we discuss the classical ill-posed problem of numerical differentiation, assuming that the smoothness of the function to be differentiated is unknown. Using recent results on adaptive regularization of general ill-posed problems, we propose new rules for the choice of the stepsize in the finite-difference methods, and for the regularization parameter choice in numerical differentiation regularized by the iterated Tikhonov method. These methods are shown to be effective for the differentiation of noisy functions, and the order-optimal convergence results for them are proved.

     

Assessment of thermal expansion coefficient of methylsilsesquioxane thin film

For the methylsilsesquioxane film whose optical birefringence is almost zero, it was recently reported that its vertical thermal expansion coefficient (CTE) was approximately one order of magnitude larger than the lateral CTE. Though the birefringence is not an absolute predictor of anisotropic behavior, the discrepancy in both the CTEs was so remarkable that it was essential to investigate whether the anisotropy was intrinsic property or not. If the effect of Poisson's ratio is considered in the calculation of the vertical CTE and when elastic modulus measured by surface acoustic wave spectroscopy is used in the assessment of the lateral CTE, both the CTEs are coincident with each other. Therefore, it can be concluded that the discrepancy in the CTEs can be attributed to a higher in‐plane polymer chain orientation but it can also arise from the misleadingly assumed modulus and Poisson's ratio. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 3109–3120, 2006     

Atomistic packing models for experimentally investigated swelling states induced by CO2 in glassy polysulfone and poly(ether sulfone)

Atomistic packing models have been created, which help to better understand the experimentally observed swelling behavior of glassy polysulfone and poly (ether sulfone), under CO₂ gas pressures up to 50 bar at 308 K. The experimental characterization includes the measurement of the time‐dependent volume dilation of the polymer samples after a pressure step and the determination of the corresponding gas concentrations by gravimetric gas‐sorption measurements. The models obtained by force‐field‐based molecular mechanics and molecular dynamics methods allow a detailed atomistic analysis of representative swelling states of polymer/gas systems, with respect to the dilation of the matrix. Also, changes of free volume distribution and backbone mobility are accessible. The behavior of gas molecules in unswollen and swollen polymer matrices is characterized in terms of sorption, diffusion, and plasticization. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 1874–1897, 2006     

Approximation of Matrices and a Family of Gander Methods for Polar Decomposition

Two matrix approximation problems are considered: approximation of a rectangular complex matrix by subunitary matrices with respect to unitarily invariant norms and a minimal rank approximation with respect to the spectral norm. A characterization of a subunitary approximant of a square matrix with respect to the Schatten norms, given by Maher, is extended to the case of rectangular matrices and arbitrary unitarily invariant norms. Iterative methods, based on the family of Gander methods and on Higham’s scaled method for polar decomposition of a matrix, are proposed for computing subunitary and minimal rank approximants. Properties of Gander methods are investigated in details. AMS subject classification (2000) 65F30, 15A18     

On disjoint configurations in infinite graphs

For a graph A and a positive integer n, let nA denote the union of n disjoint copies of A; similarly, the union of ?₀ disjoint copies of A is referred to as ?₀A. It is shown that there exist (connected) graphs A and G such that nA is a minor of G for all n??, but ?₀A is not a minor of G. This supplements previous examples showing that analogous statements are true if, instead of minors, isomorphic embeddings or topological minors are considered. The construction of A and G is based on the fact that there exist (infinite) graphs G₁, G₂,… such that G_i is not a minor of G_j for all i ≠ j. In contrast to previous examples concerning isomorphic embeddings and topological minors, the graphs A and G presented here are not locally finite. The following conjecture is suggested: for each locally finite connected graph A and each graph G, if nA is a minor of G for all n ? ?, then ?₀A is a minor of G, too. If true, this would be a far‐reaching generalization of a classical result of R. Halin on families of disjoint one‐way infinite paths in graphs. © 2002 Wiley Periodicals, Inc. J Graph Theory 39: 222–229, 2002; DOI 10.1002/jgt.10016     

Dynamics of liquid membranes. II: Adaptive finite difference methods

Two domain-adaptive finite difference methods are presented and applied to study the dynamic response of incompressible, inviscid, axisymmetric liquid membranes subject to imposed sinusoidal pressure oscillations. Both finite difference methods map the time-dependent physical domain whose downstream boundary is unknown onto a fixed computational domain. The location of the unknown time-dependent downstream boundary of the physical domain is determined from the continuity equation and results in an integrodifferential equation which is non-linearly coupled with the partial differential equations which govern the conservation of mass and linear momentum and the radius of the liquid membrane. One of the finite difference methods solves the non-conservative form of the governing equations by means of a block implicit iterative method. This method possesses the property that the Jacobian matrix of the convection fluxes has an eigenvalue of algebraic multiplicity equal to four and of geometric multiplicity equal to one. The second finite difference procedure also uses a block implicit iterative method, but the governing equations are written in conservation law form and contain an axial velocity which is the difference between the physical axial velocity and the grid speed. It is shown that these methods yield almost identical results and are more accurate than the non-adaptive techniques presented in Part I. It is also shown that the actual value of the pressure coefficient determined from linear analyses can be exceeded without affecting the stability and convergence of liquid membranes if the liquid membranes are subjected to sinusoidal pressure variations of sufficiently high frequencies.     

GF(2)上本原多项式的三项倍式的次数研究

研究了二元域上本原多项式的三项倍式的计数问题,通过对三项倍式进行分类,导出了本原多项式三项倍式最低次数的一个上界.利用这一结果,使构造本原多项式的次数最低的三项倍式的计算复杂性降到原来的九分之一.     

HIGH ACCURACY FINITE VOLUME ELEMENT METHOD FOR TWO-POINT BOUNDARY VALUE PROBLEM OF SECOND ORDER ORDINARY DIFFERENTIAL EQUATIONS

1　ＩｎｔｒｏｄｕｃｔｉｏｎＦｉｎｉｔｅｖｏｌｕｍｅｅｌｅｍｅｎｔｍｅｔｈｏｄｏｒＦＶＥＭｕｓｅｓａｖｏｌｕｍｅｉｎｔｅｇｒａｌｆｏｒｍｕｌａｔｉｏｎｏｆｔｈｅｄｉｆ ｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｗｉｔｈａｆｉｎｉｔｅｐａｒｔｉｔｉｏｎｉｎｇｓｅｔｏｆｖｏｌｕｍｅｔｏｄｉｓｃｒｅｔｉｚｅｔｈｅｅｑｕａｔｉｏｎ ,ｔｈｅｎｒｅ ｓｔｒｉｃｔｓｔｈｅａｄｍｉｓｓｉｂｌｅｆｕｎｃｔｉｏｎｓｔｏａｌｉｎｅａｒｆｉｎｉｔｅｅｌｅｍｅｎｔｓｐａｃｅｔｏｄｉｓｃｒｅｔｉｚｅｔｈｅｓｏｌｕｔｉｏｎ ([1 -4 ] ) .…     

NEW ALTERNATING DIRECTION FINITE ELEMENT SCHEME FOR NONLINEAR PARABOLIC EQUATION

A new alternating direction (AD) finite element (FE) scheme for 3-dimensional nonlinear parabolic equation and parabolic integro-differential equation is studied. By using AD, the 3-dimensional problem is reduced to a family of single space variable problems, calculation work is simplified; by using FE, high accuracy is kept; by using various techniques for priori estimate for differential equations such as inductive hypothesis reasoning, the difficulty arising from the nonlinearity is treated. For both FE and ADFE schemes, the convergence properties are rigorously demonstrated, the optimal H1-and L2-norm space estimates and the 0((△t)2) estimate for time variable are obtained.