Continuation of travelling‐wave solutions of the Navier–Stokes equations

An efficient way of obtaining travelling waves in a periodic fluid system is described and tested. We search for steady states in a reference frame travelling at the wave phase velocity using a first‐order pseudospectral semi‐implicit time scheme adapted to carry out the Newton's iterations. The method is compared to a standard Newton–Raphson solver and is shown to be highly efficient in performing this task, even when high‐resolution grids are used. This method is well suited to three‐dimensional calculations in cylindrical or spherical geometries. Copyright © 2006 John Wiley & Sons, Ltd.     

THERMAL AND THERMO-ELASTIC-PLASTIC RESPONSE OF CERAMIC-METAL FUNCTIONALLY GRADED MATERIALS—THERMAL SHOCK PROBLEM

The thermal and thermo-elastic-plastic response of newly developed ceramic-metalfunctionally graded materials under a thermal shock load is studied.The materials are heated at the ce-ramic surface with a sudden high-intensity heat flux input,and cooled at the metal surface with aflowing liquid nitrogen.Emphasis is placed on two aspects:(1)the influence of the graded composi-tion of the materials on the temperature and stress response;and(2)the optimum design of the gradedcomposition from a unified viewpoint of the heat insulation property and stress relaxation property.Moreover,a comparison between the thermoelastic stress and the thermo-elastic-plastic stress is alsomade to indicate the plasticity effect.     

叶栅全三维粘性反问题的数值解

本文发展了一种解叶栅全三维粘性反问题的新的数值方法.基于非正交曲线坐标与相应的非正交速度分量下完全守恒型的Navier－Stokes方程,全三维反问题规定叶片表面的无量纲压力分布反求叶型。计算中叶片表面的边界条件采用一种特殊的方式来处理，即一方面强加给定的压力分布条件,另方面叶面的几何位置在迭代过程中又是可移动的，其移动速度将与Navier—Stokes方程在当地的解联系起来，从而形成一种解定常问题的新的不定常过程.试算证明了本文方法的可行性。     

On parabolic equations with a small parameter

A singularly perturbed parabolic equation with a nonlinear right-hand side of a special form is examined. A numerical analytical study of such equations is performed.     

Effect of core-shell micelle formation on the redox properties of phenothiazine-labeled poly(ethyl glycidy ether)-block-poly(ethylene oxide)

Redox properties of phenothiazine-labeled poly(ethyl glycidy ether)-block-poly(ethylene oxide) (PT-EGE_n-b-EO_m) are reversibly changed by core-shell micelle formation. In the temperature range higher than the critical micellization temperature (cmt), the anodic potential of PT group positively shifts and concomitantly its anodic current decrease, or levels off compared to those of the reference polymer PT-EO_m without the thermo-responsive EGE_n segment. The former alteration is caused by incorporation of hydrophobic PT groups into a core of the micelle and the latter by the decrease in the diffusion coefficient of PT groups due to formation of the core-shell micelles. The cmt value and the temperature-dependent alteration in the redox properties strongly depend on the polymer structure, especially the length of thermo-responsive EGE_n segment. The electrochemically determined hydrodynamic radii of the polymer aggregates seem to be overestimated, compared to the values reported for the aggregates of other thermo-responsive polymers with similar molecular weights, implying the presence of electrochemically inactive PT groups in the copolymers having longer thermo-responsive segments.     

Quantitative Comparison of REAPDOR and TRAPDOR Experiments by Numerical Simulations and Determination of H–Al Distances in Zeolites

Quantitative H–Al distances in acid sites of two zeolites with MFI and IFR framework topology were obtained by numerical simulation of ¹H{²⁷Al} rotational echo adiabatic passage double resonance (REAPDOR) experiments. A ²⁷Al offset-dependent data set yields for each resolved ¹H NMR line a corresponding nuclear electric quadrupole coupling constant of the neighboring ²⁷Al site. This information is used for analyzing a second data set for on-resonance irradiation, where the dipolar evolution time (number of rotor cycles) was varied, to yield the ¹H–²⁷Al dipolar coupling constant. Numerical simulations indicate that the REAPDOR method does not depend significantly on the polar angles, defining the orientation of the electric field gradient tensor of ²⁷Al with respect to the Al–H dipolar vector. In contrast, the transfer of populations in double resonance sequence is sensitive to these angles, and it can be thus used to measure them.     

基于频域方法的三维物体分层重现

基于频域方法实现了数字记录和再现三维物体离轴全息图.通过频谱对应关系求得各个面的频谱分布,使用数字滤波方法,成功地消除了零级衍射和共轭像;通过改变再现距离,分别获取了三维物体各个截面的再现像.     

Numerical analysis of the Toda lattice equations

Numerical solutions to three systems of integrable evolutionary equations from the Toda lattice hierarchy are analyzed. These are the classical Toda lattice, the second local dispersive flow, and the second extended dispersive flow. Special attention is given to the properties of soliton solutions. For the equations of the second local flow, two types of solitons interacting in a special manner are found. Solutions corresponding to various initial data are qualitatively outlined.     

Joint detection of a given number of reference fragments in a quasi-periodic sequence and its partition into segments containing series of identical fragments

The problem of joint a posteriori detection of reference fragments in a quasi-periodic sequence and its partition into segments containing series of recurring fragments from the reference tuple is solved. It is assumed that (i) an ordered reference tuple of sequences to be detected is given, (ii) the number of desired fragments is known, (iii) the index of the sequence term corresponding to the beginning of a fragment is a deterministic (not random) value, and (iv) a sequence distorted by an additive uncorrelated Gaussian noise is available for observation. It is established that the problem consists in testing a set of hypotheses about the mean of a random Gaussian vector. The cardinality of the set grows exponentially as the vector dimension (i.e., the sequence length) increases. An efficient a posteriori algorithm producing a maximum-likelihood optimal solution to the problem is substantiated. Time and space complexity bounds related to the parameters of the problem are derived. The results of numerical simulation are presented.     

Theoretical and numerical analysis on multispectral bioluminescence tomography

Recently, molecular imaging has been rapidly developed to studyphysiological and pathological processes in vivo at the cellularand molecular levels. Among molecular imaging modalities, opticalimaging has attracted a major attention for its unique advantages.In this paper, we establish a mathematical framework for multispectralbioluminescence tomography (BLT) that allows simultaneous studiesof multiple optical reporters. We show solution existence, uniquenessand continuous dependence on data as well as the limiting behaviourswhen the regularization parameter approaches zero or when thepenalty parameter approaches infinity. Then, we propose twonumerical schemes for multispectral BLT and derive error estimatesfor the corresponding solutions.