Numerical Approximation of the Spectra of Phan-Thien Tanner Liquids

The linear stability of the linear Phan-Thien Tanner (PTT) fluid model is investigated for plane Poiseuille flow. The PTT model involves parameters that can be used to fit shear and extensional data, which makes it suitable for describing both polymer solutions and melts. The base flow is determined using a Chebyshev-tau method. The linear stability equations are also discretized using Chebyshev approximations to furnish a generalized eigenvalue problem. The spectrum is shown to comprise a continuous part and a discrete part. The theoretical and numerical results are validated for the UCM and Oldroyd-B models, which are special cases of the PTT model, by comparing with results in the literature. It is demonstrated that the linear PTT fluid is stable to infinitesimal disturbances with respect to the range of shear-thinning, extensional and elasticity parameters considered. The computational efficiency and accuracy of the numerical method are also investigated.     

A new class of three-variable orthogonal polynomials and their recurrences relations

A new class of three-variable orthogonai polynomials,defined as eigenfunctions of a second order PDE operator,is studied.These polynomials are orthogonal over a curved tetrahedron region, which can be seen as a mapping from a traditional tetrahedron,and can be taken as an extension of the 2-D Steiner domain.The polynomials can be viewed as Jacobi polynomials on such a domain.Three- term relations are derived explicitly.The number of the individual terms,involved in the recurrences relations,are shown to be independent on the total degree of the polynomials.The numbers now are determined to be five and seven,with respect to two conjugate variables z,(?) and a real variable r, respectively.Three examples are discussed in details,which can be regarded as the analogues of the Chebyshev polynomials of the first and the second kinds,and Legendre polynomials.     

The Use of Multi-Pass Time-of-Flight Mass Analyzers for Nuclear-Decay Spectroscopy of Mass Identified Nuclei

The masses of ions of some keV can be determined in multi-pass time-of-flight mass analyzers [1,2] with high precision. The mass accuracies thus achieved are sufficient to determine the proton and neutron numbers for most short-lived and stable nuclei [3,5]. Recording α- or γ-radiation of the investigated nuclei in delayed coincidence to the ion arrival, one thus can perform nuclear spectroscopy of selected nuclei. This revised version was published online in July 2006 with corrections to the Cover Date.     

High-Order Leap-Frog Based Discontinuous Galerkin Method for the Time-Domain Maxwell Equations on Non-Conforming Simplicial Meshes

A high-order leap-flog based non-dissipative discontinuous Galerkin timedomain method for solving Maxwell's equations is introduced and analyzed. The proposed method combines a centered approximation for the evaluation of fluxes at the interface between neighboring elements, with a Nth-order leap-frog time scheme. Moreover, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary level hanging nodes. The method is proved to be stable under some CFL-like condition on the time step. The convergence of the semi-discrete approximation to Maxwell's equations is established rigorously and bounds on the global divergence error are provided. Numerical experiments with highorder elements show the potential of the method.     

Chebyshev pseudospectral method for computing numerical solution of convection–diffusion equation

A method for computing highly accurate numerical solutions of 1D convection–diffusion equations is proposed. In this method, the equation is first discretized with respect to the spatial variable, transforming the original problem into a set of ordinary differential equations, and then the resulting system is integrated in time by the fourth-order Runge–Kutta method. Spatial discretization is done by using the Chebyshev pseudospectral collocation method. Before describing the method, we review a finite difference-based method by Salkuyeh [D. Khojasteh Salkuyeh, On the finite difference approximation to the convection–diffusion equation, Appl. Math. Comput. 179 (2006) 79–86], and, contrary to the proposal of the author, we show that this method is not suitable for problems involving time dependent boundary conditions, which calls for revision. Stability analysis based on pseudoeigenvalues to determine the maximum time step for the proposed method is also carried out. Superiority of the proposed method over a revised version of Salkuyeh’s method is verified by numerical examples.     

Parameter accuracy analysis of weak-value amplification process in the presence of noise

We theoretically introduce the statistical uncertainty of photon number and phase error to discuss the precision of parameters to be measured based on weak measurements. When the photon counting scheme is used, we discuss the relative accuracy of the system in the presence of phase error by using the orthogonal and nonorthogonal pre-and postselected states, respectively. When using the measurement scheme of pointer shift, we discuss the measurement accuracy in the presence of phase error, pointer resolution, and statistical uncertainty. These results give a guide way to get the smallest relative precision and deepen our understanding about weak measurement.     

Optimal Control of a Tethered Subsatellite of Three Degrees of Freedom

The paper presents the optimal control of the deployment and retrieval processes of a tethered subsatellite system of three degrees of freedom, which takes not only the in-plane motion, but also the out-of-plane motions, into account. After the statement of the optimal control problem of the tethered subsatellite system based on the dynamic equation of the system, with the control cost and the state constraints included, the paper introduces the quasilinearization and the truncated Chebyshev series to approximate the state variables of the system such that the original problem of constrained nonlinear optimal control is simplified into a set of linear quadratic programming problems which can be easily solved. The case studies in the paper not only support the new method, but also show that the controlled trajectories of the deployment process and the retrieval process are geometrically symmetric to each other with respect to the local vertical axis, and that the subsatellite always undergoes a slow, damped oscillation when it is in the beginning of a deployment process or at the end of a retrieval process.     

基于两类切比雪夫多项式的循环矩阵行列式的计算

利用多项式因式分解的逆变换,结合循环矩阵和切比雪夫多项式的特殊结构,首先研究第三类和第四类切比雪夫多项式的通项公式,并给出第三类、第四类切比雪夫多项式的关于行首加r尾r右循环矩阵和行尾加r首r左循环矩阵的行列式的显式表达式,最后给出算法实施步骤.     

A study of the precision and accuracy of peak quantification in comprehensive two-dimensional liquid chromatography in time

Simulated chromatographic data were used to determine the precision and accuracy in the estimation of peak volumes (i.e., peak sizes) in comprehensive two-dimensional liquid chromatography in time (LC × LC). Peak volumes were determined both by summing the areas in the second dimension chromatograms and by fitting the second dimension areas to a Gaussian peak. The Gaussian method is better at predicting the peak volume than the moments method provided there are at least three second dimension injections above the limit of detection (LOD). However, when only two of the second dimension signals are substantially above baseline, the accuracy and precision of the Gaussian fit method become quite poor because the results from the fitting algorithm become indeterminate. Based on simulations in which the modulation ratio (M_R = 4¹σ/t_s) and sampling phase (?) were varied, we conclude for well-resolved peaks that the optimum precision in peak volumes in 2D separations will be obtained when the M_R is between two and five, such that there are typically four to ten second dimension peaks recorded over the eight σ width of the first dimension peak. This sampling rate is similar to that suggested by the Murphy–Schure–Foley criterion. This provides an RSD of approximately 2% for the signal-to-noise ratio used in the present simulations. The precision of the peak volume of experimental data was also assessed, and RSD values were in the range of 4–5%. We conclude that the poorer precision found in the LC × LC experimental data as compared to LC may be due to experimental imprecision in sampling the effluent from the first dimension column.     

标准镜的高精度柔性支撑镜框的优化设计

为了实现干涉仪标准镜中光学元件的高精度定位,设计了一种柔性支撑镜框,研究了该结构的力学模型、结构参数、定位精度和透镜变形。首先,根据材料力学原理将柔性镜框等效为一个弹簧系统;根据力学方程和几何关系,建立了透镜中心位置与柔性结构的挠度之间的二元方程。然后,分析了安装位置、温度、结构参数对透镜位置以及作用力的影响。最后,应用有限元仿真分析了所设计结构的力学性能,并进行对比验证。结果表明,数值仿真分析的结果与有限元仿真结果基本相同,柔性镜框的柔性结构厚度最优值为1.5 mm。该设计方案完全满足干涉仪标准镜对镜框在定位精度、稳定性方面的要求。