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.     

Identifying aspirin polymorphs from combined DFT-based crystal structure prediction and solid-state NMR

A combined experimental and computational approach was used to distinguish between different polymorphs of the pharmaceutical drug aspirin. This method involves the use of ab initio random structure searching (AIRSS), a density functional theory (DFT)-based crystal structure prediction method for the high-accuracy prediction of polymorphic structures, with DFT calculations of nuclear magnetic resonance (NMR) parameters and solid-state NMR experiments at natural abundance. AIRSS was used to predict the crystal structures of form-I and form-II of aspirin. The root-mean-square deviation between experimental and calculated ¹H chemical shifts was used to identify form-I as the polymorph present in the experimental sample, the selection being successful despite the large similarities between the molecular environments in the crystals of the two polymorphs.     

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.     

Second-order random wave solutions for interfacial internal waves in N-layer density-stratified fluid

This paper studies the random internal wave equations describing the density interface displacements and the velocity potentials of N-layer stratified fluid contained between two rigid walls at the top and bottom. The density interface displacements and the velocity potentials were solved to the second-order by an expansion approach used by Longuet-Higgins （1963） and Dean （1979） in the study of random surface waves and by Song （2004） in the study of second- order random wave solutions for internal waves in a two-layer fluid. The obtained results indicate that the first-order solutions are a linear superposition of many wave components with different amplitudes, wave numbers and frequencies, and that the amplitudes of first-order wave components with the same wave numbers and frequencies between the adjacent density interfaces are modulated by each other. They also show that the second-order solutions consist of two parts： the first one is the first-order solutions, and the second one is the solutions of the second-order asymptotic equations, which describe the second-order nonlinear modification and the second-order wave-wave interactions not only among the wave components on same density interfaces but also among the wave components between the adjacent density interfaces. Both the first-order and second-order solutions depend on the density and depth of each layer. It is also deduced that the results of the present work include those derived by Song （2004） for second-order random wave solutions for internal waves in a two-layer fluid as a particular case.     

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.     

On a stochastic nonlinear equation in one-dimensional viscoelasticity

In this paper we discuss an initial-boundary value problem for a stochastic nonlinear equation arising in one-dimensional viscoelasticity. We propose to use a new direct method to obtain a solution. This method is expected to be applicable to a broad class of nonlinear stochastic partial differential equations.

     

Generic replica symmetric field-theory for short range Ising spin glasses

Symmetry considerations and a direct, Hubbard-Stratonovich type, derivation are used to construct a replica field-theory relevant to the study of the spin glass transition of short range models in a magnetic field. A mean-field treatment reveals that two different types of transitions exist, whenever the replica number n is kept larger than zero. The Sherrington-Kirkpatrick critical point in zero magnetic field between the paramagnet and replica magnet (a replica symmetric phase with a nonzero spin glass order parameter) separates from the de Almeida-Thouless line, along which replica symmetry breaking occurs. We argue that for studying the de Almeida-Thouless transition around the upper critical dimension d = 6, it is necessary to use the generic cubic model with all the three bare masses and eight cubic couplings. The critical role n may play is also emphasized. To make perturbative calculations feasible, a new representation of the cubic interaction is introduced. To illustrate the method, we compute the masses in one-loop order. Some technical details and a list of vertex rules are presented to help future renormalisation-group calculations. Received 9 October 2001     

Diffusion from sessile droplets through membranes

We analyze diffusion from a periodic array of hemispherical droplets through a membrane. We find that the multiple sources do not interact strongly, even when the droplets are closely spaced, so that the flux through the membrane appears nearly additive.     

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.

     

仓库容量有限条件下的随机存贮管理问题的评注

给出2005年全国部分高校研究生数学建模竞赛D题的研究背景,对参赛论文作了简要评注,并提供了一种参考解答.