Représentation d''une particule quantique de spin 1/2 par le produit de Moyal et une super-variété symplectique plate

The aim of this paper is to show that it is possible to represent a quantum particle with a spin with the help of a star product and a flat supermanifold. We extend to the supermanifolds the Schouten bracket, the exterior differentiation, the Lie derivative. These notions allow us to define the Poisson bracket and the symplectic supermanifolds. It is easy to define the Moyal product: the Grassman algebra of the superfunctions becomes a Clifford algebra. We suppose that the observables of a quantum particle with a spin are superfunctions defined on a flat supermanifold. The study of these star algebras proves the existence of a spectral resolution for some particular elements. Consequently the equation of motion for an observable admits stationary solutions. The examples considered show that the spin splits a given state into two other ones.     

Geometric Analysis and Symbol Calculus: Fourier Transform Magnetic Resonance Imaging and Wavelets

Due to its unequalled advantages, the magnetic resonance imaging (MRI) modality has truly revolutionized the diagnosis and evaluation of pathology. Because many morphological anatomic details that may not be visualized by other high tech imaging methods can now be readily shown by diagnostic MRI, it has already become the standard modality by which all other clinical imaging techniques are measured. The unique quantum physical basis of the MRI modality combined with the imaging capabilities of current computer technology has made this imaging modality a target of interdisciplinary interest for clinicians, physicists, biologists, engineers, and mathematicians. Due to the fact that MRI scanners perform corticomorphic processing, this modality is by far more complex than all the other high tech clinical imaging techniques. The purpose of this paper is to outline a phase coherent wavelet approach to Fourier transform MRI. It is based on distributional harmonic analysis on the Heisenberg nilpotent Lie group G and the associated symplectically invariant symbol calculus of pseudodifferential operators. The contour and contrast resolution of MRI scans which is controlled by symplectic filter bank processing gives the noninvasive MRI modality superiority over X-ray computed tomography (CT) in soft tissue differentiation.     

On theoretical and numerical aspects of symplectic Gram–Schmidt-like algorithms

Gram–Schmidt-like orthogonalization process with respect to a given skew-symmetric scalar product is a key step in model reduction methods, structure-preserving, for large sparse Hamiltonian eigenvalue problem. Theoretical as well as numerical aspects of this step do not benefit of enough attention, compared to the one allowed to the classical Gram–Schmidt algorithm and its modified version. The aim of this paper is to revisit the symplectic Gram–Schmidt algorithms, to built some modified versions and to deal with their theoretical and numerical features. AMS subject classification 65F15, 65F50     

关于辛算法稳定性的若干注记

我们讨论辛算法的线性稳定性和非线性稳定性,从动力系统和计算的角度论述了研究辛算法的这两类稳定性问题的重要性,分析总结了相关重要结果.我们给出了解析方法的明确定义,证明了稳定函数是亚纯函数的解析辛方法是绝对线性稳定的.绝对线性稳定的辛方法既有解析方法（如Runge-Kutta辛方法）,也有非解析方法（如基于常数变易公式对线性部分进行指数积分而对非线性部分使用其它数值积分的方法）.我们特别回顾并讨论了R.I.McLachlan,S.K.Gray和S.Blanes,F.Casas,A.Murua等关于分裂算法的线性稳定性结果,如通过选取适当的稳定多项式函数构造具有最优线性稳定性的任意高阶分裂辛算法和高效共轭校正辛算法,这类经优化后的方法应用于诸如高振荡系统和波动方程等线性方程或者线性主导的弱非线性方程具有良好的数值稳定性.我们通过分析辛算法在保持椭圆平衡点的稳定性,能量面的指数长时间慢扩散和KAM不变环面的保持等三个方面阐述了辛算法的非线性稳定性,总结了相关已有结果.最后在向后误差分析基础上,基于一个自由度的非线性振子和同宿轨分析法讨论了辛算法的非线性稳定性,提出了一个新的非线性稳定性概念,目的是为辛算法提供一个实际可用的非线性稳定性判别法.     

A condition analysis of the weighted linear least squares problem using dual norms

In this paper, based on the theory of adjoint operators and dual norms, we define condition numbers for a linear solution function of the weighted linear least squares problem. The explicit expressions of the normwise and componentwise condition numbers derived in this paper can be computed at low cost when the dimension of the linear function is low due to dual operator theory. Moreover, we use the augmented system to perform a componentwise perturbation analysis of the solution and residual of the weighted linear least squares problems. We also propose two efficient condition number estimators. Our numerical experiments demonstrate that our condition numbers give accurate perturbation bounds and can reveal the conditioning of individual components of the solution. Our condition number estimators are accurate as well as efficient.     

计算最优控制辛数值方法

针对最优控制问题（OCP）的辛数值方法研究及应用进行综述。主要涉及内容包括,动力学系统为常微分方程描述的一般无约束、含不等式约束和状态时滞的最优控制问题,微分代数方程描述的一般无约束、含不等式约束和含切换系统的最优控制问题,以及闭环最优控制问题。从间接法和直接法两个求解框架出发,重点介绍本课题组在保辛算法方面的研究工作。在间接法框架下,首先基于生成函数和变分原理,将OCP保辛离散为非线性方程组,再数值求解方程组。在直接法框架下,将OCP保辛离散为有限维的非线性规划问题（NLP）,再数值求解。针对闭环最优控制问题,提出了保辛模型预测控制、滚动时域估计和瞬时最优控制算法。研究表明,保辛算法具有高精度和高效率的特点,在航空航天和机器人等领域有着广泛应用前景和价值。     

An Improved Hybrid Adjoint Method in External Aerodynamics Using Variational Technique for the Boundary Integral Based Optimal Objective Function Gradient

An improved hybrid adjoint method to the viscous, compressible Reynold-Averaged Navier-Stokes Equation (RANS) is developed for the computation of objective function gradient and demonstrated for external aerodynamic design optimization. In this paper, the main idea is to extend the previous coupling of the discrete and continuous adjoint method by the grid-node coordinates variation technique for the computation of the variation in the gradients of flow variables. This approach in combination with the Jacobian matrices of flow fluxes refrained the objective function from field integrals and coordinate transformation matrix. Thus, it opens up the possibility of employing the hybrid adjoint method to evaluate the subsequent objective function gradient analogous to many shape parameters, comprises of only boundary integrals. This avoids the grid regeneration in the geometry for every surface perturbation in a structured and unstructured grid. Hence, this viable technique reduces the overall CPU cost. Moreover, the new hybrid adjoint method has been successfully applied to the computation of accurate sensitivity derivatives. Finally, for the investigation of the presented numerical method, simulations are carried out on NACA0012 airfoil in a transonic regime and its accuracy and effectiveness related to the new gradient equation have been verified with the Finite Difference Method (FDM). The analysis reveals that the presented methodology for the optimization provides the designer with an indispensable CPU-cost effective tool to reshape the complex geometry airfoil surfaces, useful relative to the state-of-the-art, in a less computing time.     

Symplectic Runge-Kutta methods of high order based on W-transformation

In this paper , characterizations of symmetric and symplectic Runge-Kutta methods based on the W-transformation of Hairer and Wanner are presented. Using these characterizations, we construct two families symplectic (symmetric and algebraically stable or algebraically stable) Runge-Kutta methods of high order. Methods constructed in this way and presented in this paper include and extend the known classes of high order implicit Runge-Kutta methods.