Artificial time integration

Many recent algorithmic approaches involve the construction of a differential equation model for computational purposes, typically by introducing an artificial time variable. The actual computational model involves a discretization of the now time-dependent differential system, usually employing forward Euler. The resulting dynamics of such an algorithm is then a discrete dynamics, and it is expected to be “close enough” to the dynamics of the continuous system (which is typically easier to analyze) provided that small – hence many – time steps, or iterations, are taken. Indeed, recent papers in inverse problems and image processing routinely report results requiring thousands of iterations to converge. This makes one wonder if and how the computational modeling process can be improved to better reflect the actual properties sought. In this article we elaborate on several problem instances that illustrate the above observations. Algorithms may often lend themselves to a dual interpretation, in terms of a simply discretized differential equation with artificial time and in terms of a simple optimization algorithm; such a dual interpretation can be advantageous. We show how a broader computational modeling approach may possibly lead to algorithms with improved efficiency. AMS subject classification (2000)  65L05, 65M32, 65N21, 65N22, 65D18     

天光一号预放大器中的放大自发辐射效应

 在三维放大自发辐射（ASE）的模拟计算中，计入了侧壁反射效应以及KrF自发辐射的谱线形状（高斯型、洛伦兹型及平顶型）的影响，并对空间网点的影响进行了考察。经过综合分析，给出了ASE对预放大器放大特性影响的上、下限范围。当计入平顶型频谱分布，壁反射系数为 0 .2时，与不计入ASE的情形对比，天光一号预放大器的提取效率下降为35%。     

Sensitivity of depolarized lidar signals to cloud and aerosol particle properties

Measurements from depolarized lidars provide a promising method to retrieve both cloud and aerosol properties and a versatile complement to passive satellite-based sensors. For lidar observations of clouds and aerosols, multiple scattering plays an important role in the scattering process. Monte Carlo simulations are carried out to investigate the sensitivity of lidar backscattering depolarization to cloud and aerosol properties. Lidar parameters are chosen to be similar to those of the upcoming space-based CALIPSO lidar. Cases are considered that consist of a single cloud or aerosol layer, as well as a case in which cirrus clouds overlay different types of aerosols. It is demonstrated that besides thermodynamic cloud phase, the depolarized lidar signal may provide additional information on ice or aerosol particle shapes. However, our results show little sensitivity to ice or aerosol particle sizes. Additionally, for the case of multiple but overlapping layers involving both clouds and aerosols, the depolarized lidar contains information that can help identify the particle properties of each layer.     

Classical, phenomenological analysis of the kinetics of reactions at the gas-exposed surface of mixed ionic electronic conductors

The kinetics of reactions occurring at the gas-exposed surfaces of charged mixed ionic electronic conductors (MIECs) are examined from theoretical first principles. Analysis based on the classical electrochemical potential-transition state theory model reveals that the nature of the reactions is electrochemical in general. However, the influence of the surface potential on the reaction rate is opposite for adsorption and incorporation reactions. Two-dimensional finite volume models of an MIEC as working electrode in a half-cell configuration are presented. The results for a simple, two-step reduction process show that the effect of the surface potential on the rate of reactions is minimal for incorporation-limited reactions but more influential for adsorption-limited reactions. An erratum to this article is available at .     

用板式电势差计测量电池的电动势和内阻实验的改进

分析了在“用板式电势差计测量电池的电动势和内阻”实验中 ,由于待测干电池的电动势和内阻在实验过程中均不为恒定值而引起的误差 ,依据电化学理论和电阻构成解释了实验现象 ,并给出了实验改进方案 .     

Shape sensitivity of a plane crack front

The 3D‐elasticity model of a solid with a plane crack under the stress‐free boundary conditions at the crack is considered. We investigate variations of a solution and of energy functionals with respect to perturbations of the crack front in the plane. The corresponding expansions at least up to the second‐order terms are obtained. The strong derivatives of the solution are constructed as an iterative solution of the same elasticity problem with specified right‐hand sides. Using the expansion of the potential and surface energy, we consider an approximate quadratic form for local shape optimization of the crack front defined by the Griffith criterion. To specify its properties, a procedure of discrete optimization is proposed, which reduces to a matrix variational inequality. At least for a small load we prove its solvability and find a quasi‐static model of the crack growth depending on the loading parameter. Copyright © 2003 John Wiley & Sons, Ltd.     

Shape Optimization and Fictitious Domain Approach for Solving Free Boundary Problems of Bernoulli Type

This contribution deals with an efficient method for the numerical realization of the exterior and interior Bernoulli free boundary problems. It is based on a shape optimization approach. The state problems are solved by a fictitious domain solver using boundary Lagrange multipliers.     

短脉冲强激光在稀薄等离子体中传播时的电磁类孤立子

 采用二维粒子模拟方法，研究了短脉冲强激光在稀薄等离子体中传播时电磁类孤立子的产生和时空演化过程。通过分析电磁场与等离子体波的非线性能量交换和激光场的频率谱结构等，给出了电磁类孤立子形成的基本物理图像，讨论了等离子体参数对电磁类孤立子形成的影响。模拟结果表明：类孤立子的形成是由于局部电磁波振荡频率减小至等离子体频率引起的，初始等离子体密度越高越容易形成空间局域结构。     

FTIR measurements of CO2 line positions and intensities at high temperature in the 3700-3750 cm spectral region

Positions and intensities for 453 spectral lines in 12 rovibrational bands of ¹²C¹⁶O₂ have been determined between 3700 and 3750 cm⁻¹. At three temperatures (294, 500, and 698 K) eight spectra have been recorded at a pressure around 5 mbar and for an absorption path of about 190 cm⁻¹ using a Bomen DA3 Fourier transform spectrometer (4 × 10⁻³ cm⁻¹ resolution). Some of the measured positions and intensities can be compared with recent experimental results that validate the experimental set-up and the data analysis procedure. The results are also compared with the values listed in the HITRAN 2000 database. If the agreement is generally good, discrepancies are observed for three hot bands.     

Embedding of Dynamical Symmetry Groups U(1, 1) and U(2) of a Free Particle on AdS 2 and S 2 into Parasupersymmetry Algebra

Using two different types of the laddering equations realized simultaneously by the associated Gegenbauer functions, we show that all quantum states corresponding to the motion of a free particle on AdS ₂ and S ² are splitted into infinite direct sums of infinite-and finite-dimensional Hilbert subspaces which represent Lie algebras u(1, 1) and u(2) with infinite- and finite-fold degeneracies, respectively. In addition, it is shown that the representation bases of Lie algebras with rank 1, i.e., gl(2, C), realize the representation of nonunitary parasupersymmetry algebra of arbitrary order. The realization of the representation of parasupersymmetry algebra by the Hilbert subspaces which describe the motion of a free particle on AdS ₂ and S ² with the dynamical symmetry groups U(1, 1) and U(2) are concluded as well.