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.     

Analogue Models of and for Gravity

Condensed matter systems, such as acoustics in flowing fluids, light in moving dielectrics, or quasiparticles in a moving superfluid, can be used to mimic aspects of general relativity. More precisely these systems (and others) provide experimentally accessible models of curved-space quantum field theory. As such they mimic kinematic aspects of general relativity, though typically they do not mimic the dynamics. Although these analogue models are thereby limited in their ability to duplicate all the effects of Einstein gravity they nevertheless are extremely important—they provide black hole analogues (some of which have already been seen experimentally) and lead to tests of basic principles of curved-space quantum field theory. Currently these tests are still in the realm of gedanken-experiments, but there are plausible candidate models that should lead to laboratory experiments in the not too distant future.     

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.     

Radiation bursts from particles in the field of compact, impenetrable, astrophysical objects

The radiation emitted by charged, scalar particles in a Schwarzschild field with maximal acceleration corrections is calculated classically and in the tree approximation of quantum field theory. In both instances the particles emit radiation that has characteristics similar to those of gamma-ray bursters.     

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.     

超强脉冲激光在低密度等离子体中的相对论自导引效应

分析了相对论条件下激光超短脉冲在等离子体中的传输特性 ,在傍轴近似和慢变振幅近似条件下 ,推导了折射率、电子密度、静电场以及电子空腔尺度的表达式。当激光功率超过产生自导引阈值功率时 ,激光束斑沿着传输光轴方向振荡。在有质动力产生的压力非常强时 ,聚焦光束中央部分的电子被全部排开形成电子空腔。给出了电子空腔的尺寸以及在出现电子空腔时的处理方法。在超过形成电子空腔的阈值功率 (Pc≈ 2 .5TW )时 ,空腔的尺度几乎与激光功率无关 ,这意味着电子空腔阻止了激光脉冲的进一步聚焦。     

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.

     

Smart Geodesic Tracing in GRworkbench

We have developed a new tool for numerical work in General Relativity: GRworkbench. We discuss how GRworkbench's implementation of a numerically-amenable analogue to Differential Geometry facilitates the development of robust and chart-independent numerical algorithms. We consider, as an example, geodesic tracing on two charts covering the exterior Schwarzschild space-time.     

关于AOR迭代法的研究

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

Discretization of the solutions to Poisson’s equation

Sharp estimates (in the power scale) are obtained for the discretization error in the solutions to Poisson’s equation whose right-hand side belongs to a Korobov class. Compared to the well-known Korobov estimate, the order is almost doubled and has an ultimate value in the power scale.