A posteriori error analysis of a quadratic finite volume method for nonlinear elliptic problems

In this article, we construct and analyze a residual-based a posteriori error estimator for a quadratic finite volume method (FVM) for solving nonlinear elliptic partial differential equations with homogeneous Dirichlet boundary conditions. We shall prove that the a posteriori error estimator yields the global upper and local lower bounds for the norm error of the FVM. So that the a posteriori error estimator is equivalent to the true error in a certain sense. Numerical experiments are performed to illustrate the theoretical results.     

Setting local rank constraints by orthogonal projections for image resolution analysis: Application to the determination of a low dose pharmaceutical compound

Raman chemical imaging provides chemical and spatial information about pharmaceutical drug product. By using resolution methods on acquired spectra, the objective is to calculate pure spectra and distribution maps of image compounds. With multivariate curve resolution-alternating least squares, constraints are used to improve the performance of the resolution and to decrease the ambiguity linked to the final solution. Non negativity and spatial local rank constraints have been identified as the most powerful constraints to be used.     

基于科赫曲线分形结构太赫兹带通滤波器

太赫兹滤波器是太赫兹通信、太赫兹成像和太赫兹检测等太赫兹应用系统中不可或缺的功能器件。按照不同的分类方式,滤波器有不同的种类,常见的按照选频功能可分为高通滤波器、低通滤波器、带阻滤波器和带通滤波器。为了实现在太赫兹波段的滤波效果,世界各地的研究人员利用不同的结构、材料和控制方式实现了功能各异的太赫兹滤波器,但是考虑到设计的器件要应用到太赫兹系统中,成本低廉、结构简单、性能优越的太赫兹滤波器一直是研究人员的追求。分形概念自提出以来在很多研究领域都有了快速发展,但是在太赫兹波段的应用还不是很常见,特别是应用于太赫兹功能器件的设计。引入分形中科赫曲线的概念设计并制备了一种新型的太赫兹带通滤波器,该滤波器是在金属薄膜上刻蚀出科赫曲线分形结构,当太赫兹波垂直入射到该滤波器时候实现了在太赫兹波段的窄带滤波。在滤波器的设计过程中,追求理论与实验相结合,首先在电磁仿真软件中建立科赫曲线分形结构滤波器模型进行计算,探究分形结构应用于太赫兹波段进行滤波的可行性,在进行多次计算之后得到优化后的尺寸和结构,然后根据优化后的尺寸加工出科赫曲线分形结构太赫兹滤波器样品,并且将样品放在太赫兹时域光谱系统中进行实验测量,得到实验数据后与仿真结果进行比较。在仿真中利用了时域有限差分法模拟科赫曲线分形结构太赫兹带通滤波器的传输特性,优化后的仿真结果表明:滤波器的谐振频率为0.715 THz,透射系数能够达到0.92,-3 dB带宽为21.9 GHz,将仿真得到的散射参数进行S参数反演得到了太赫兹滤波器样品的电磁参数,这在理论上分析了太赫兹波在谐振点处产生透射增强的原因。利用飞秒激光微加工系统制备了尺寸优化后的科赫曲线分形结构太赫兹带通滤波器样品,然后使用太赫兹时域光谱系统对样品的传输特性进行测试,对实验得到的时域数据进行快速傅里叶变换之后得到频域数据,再把频域数据进行归一化处理后与之前的电磁仿真结果进行对比,发现实验测得的结果与电磁软件仿真得到的结果较为吻合。     

Definition and detection of data-based uniqueness in evaluating bilinear (two-way) chemical measurements

Multivariate curve resolution methods, frequently used in analyzing bilinear data sets, result in ambiguous decomposition in general. Implementing the adequate constraints may lead to reduce the so-called rotational ambiguity drastically, and in the most favorable cases to the unique solution. However, in some special cases, non-negativity constraint as minimal information of the system is a sufficient condition to resolve profiles uniquely. Although, several studies on exploring the uniqueness of the bilinear non-negatively constrained multivariate curve resolution methods have been made in the literature, it has still remained a mysterious question. In 1995, Manne published his profile-based theorems giving the necessary and sufficient conditions of the unique resolution. In this study, a new term, i.e., data-based uniqueness is defined and investigated in details, and a general procedure is suggested for detection of uniquely recovered profile(s) on the basis of data set structure in the abstract space. Close inspection of Borgen plots of these data sets leads to realize the comprehensive information of local rank, and these argumentations furnish a basis for data-based uniqueness theorem. The reported phenomenon and its exploration is a new stage (it can be said fundament) in understanding and describing the bilinear (matrix-type) chemical data in general.     

Thermally induced solid-state transformation of cimetidine. A multi-spectroscopic/chemometrics determination of the kinetics of the process and structural elucidation of one of the products as a stable N3-enamino tautomer

Exposure of cimetidine (CIM) to dry heat (160–180 °C) afforded, upon cooling, a glassy solid containing new and hitherto unknown products. The kinetics of this process was studied by a second order chemometrics-assisted multi-spectroscopic approach. Proton and carbon-13 nuclear magnetic resonance (NMR), as well as ultraviolet and infrared spectroscopic data were jointly used, whereas multivariate curve resolution with alternating least squares (MCR-ALS) was employed as the chemometrics method to extract process information. It was established that drug degradation follows a first order kinetics.     

Efficient Collocation Schemes for Singular Boundary Value Problems

We discuss an error estimation procedure for the global error of collocation schemes applied to solve singular boundary value problems with a singularity of the first kind. This a posteriori estimate of the global error was proposed by Stetter in 1978 and is based on the idea of Defect Correction, originally due to Zadunaisky. Here, we present a new, carefully designed modification of this error estimate which not only results in less computational work but also appears to perform satisfactorily for singular problems. We give a full analytical justification for the asymptotical correctness of the error estimate when it is applied to a general nonlinear regular problem. For the singular case, we are presently only able to provide computational evidence for the full convergence order, the related analysis is still work in progress. This global estimate is the basis for a grid selection routine in which the grid is modified with the aim to equidistribute the global error. This procedure yields meshes suitable for an efficient numerical solution. Most importantly, we observe that the grid is refined in a way reflecting only the behavior of the solution and remains unaffected by the unsmooth direction field close to the singular point.     

A likelihood-based method for haplotype association studies of case-control data with genotyping uncertainty

This paper discusses the associations between traits and haplotypes based on Fl (fluorescent intensity) data sets. We consider a clustering algorithm based on mixtures of t distributions to obtain all possible genotypes of each individual (i.e. "GenoSpec-trum"). We then propose a likelihood-based approach that incorporates the genotyping uncertainty to assessing the associations between traits and haplotypes through a haplotype-based logistic regression model. Simulation studies show that our likelihood-based method can reduce the impact induced by genotyping errors.     

An efficient numerical scheme for precise time integration of a diffusion-dissolution/precipitation chemical system

A numerical scheme based on an operator splitting method and a dense output event location algorithm is proposed to integrate a diffusion-dissolution/precipitation chemical initial-boundary value problem with jumping nonlinearities. The numerical analysis of the scheme is carried out and it is proved to be of order 2 in time. This global order estimate is illustrated numerically on a test case.

     

关于AOR迭代法的研究

本文论证了严格对角占优矩阵之AOR法的误差估计式中的误差估计常数hγ,ω（0≤γ≤ω0）的最小值是h1,1.