Looking at Fokker-Planck dynamics with a noisy instrument

We consider the class of experiments which can be characterized by a Fokker-PIanck dynamics corresponding to the overdamped motion of a state point in a suitable stochastic potential. We assume that the general form of the potential is known (or can be guessed with reasonable accuracy), but that its parameters are to be determined experimentally by measurements made with a noisy instrument. This possible method for determining the potential parameters, which exploits the system's own internal stochastic motion in order to explore rapidly its available parameter space, is substantially more efficient than traditional methods involving time averages of single point measurements, and yet does not appear to have been previously considered. The method could be important when, for example, the experiment must be completed in a limited time owing either to the expense of the experimental materials or to the temporary stationarity of the preparation, situations which are commonly encountered in experimental biochemistry and biology.     

Time-marching numerical schemes for the electric field integral equation on a straight thin wire

We derive and analyse four algorithms for computing the current induced on a thin straight wire by a transient electric field. They all involve solving the thin wire electric field integral equations (EFIEs) and consist of a very accurate differential equations solver together with various schemes to approximate the vector potential integral equation. We carry out a rigorous numerical stability analysis of each of these methods. This has not previously been done for solution schemes for the thin wire EFIEs. Each scheme is shown to be stable and convergent provided the radius of the wire is small enough for the thin wire equations to be a valid model.     

Quantum Thermal Effect of Dirac Particles in a Nonuniformly Rectilinearly Accelerating Black Hole with Electric Charge, Magnetic Charge, and Cosmological Constant

The Hawking radiation of Dirac particles in an arbitrarily rectilinearly accelerating Kinnersley black hole with electromagnetic charge and cosmological constant is investigated by using the generalized tortoise coordinate transformation. Both the location and the temperature of the event horizon depend on the time and the polar angle. The Hawking thermal radiation spectrum of Dirac particles is also derived.     

On the Solvability of the Multidimensional Version of the First Darboux Problem for a Model Second-Order Degenerating Hyperbolic Equation

A multidimensional version of the first Darboux problem is considered for a model second-order degenerating hyperbolic equation. Using the technique of functional spaces with a negative norm, the correct formulation of this problem in the Sobolev weighted space is proved.     

Classical limit of relativistic quantal system with attractive Coulomb interaction

Using the appropriate harmonic oscillator states and reasonable approximations, we construct coherent wavepackets corresponding to the solutions of the Klein-Gordon equation for the attractive potentialV(r)=−k/r, k>0, in two and three space dimensions. We deduce the corresponding classical limit in two dimension by requiring that the expectation value 〈r〉 of the radial variable is large. In the case of three dimensions, besides the condition of large 〈r〉, we make the uncertainty Δr=[〈r ²〉 − 〈r〉²]^1/2 a minimum with respect to certain parameter of the wavepacket. We then investigate the trajectory traversed by the wavepacket in the classical limit. We find that the classical limit of this relativistic quantal problem gives, in the leading order, the same expression for the rate of motion of the perihelion as that given by the solution of the corresponding special relativistic classical dynamical problem. We also briefly discuss some of the subtle aspects of the classical limit of the relativistic quantal system, in general.     

On a fast direct elliptic solver by a modified Fourier method

We describe high order numerical algorithms for the solution of second order elliptic equations in rectangular domains. These algorithms are based on the Fourier method in combination with a subtraction procedure. The singularities at the corner points, arising due to non-smoothness of the boundaries, are treated explicitly using properly constructed singular corner functions. The present algorithm is a generalization of the Fast Poisson Solver developed in our previous paper. This revised version was published online in August 2006 with corrections to the Cover Date.     

Cyclic reduction and FACR methods for piecewise hermite bicubic orthogonal spline collocation

Cyclic reduction and Fourier analysis-cyclic reduction (FACR) methods are presented for the solution of the linear systems which arise when orthogonal spline collocation with piecewise Hermite bicubics is applied to boundary value problems for certain separable partial differential equations on a rectangle. On anN×N uniform partition, the cyclic reduction and Fourier analysis-cyclic reduction methods requireO(N ²log₂ N) andO(N ²log₂log₂ N) arithmetic operations, respectively.     

Trefftz,collocation, and other boundary methods—A comparison

In this article we survey the Trefftz method (TM), the collocation method (CM), and the collocation Trefftz method (CTM). We also review the coupling techniques for the interzonal conditions, which include the indirect Trefftz method, the original Trefftz method, the penalty plus hybrid Trefftz method, and the direct Trefftz method. Other boundary methods are also briefly described. Key issues in these algorithms, including the error analysis, are addressed. New numerical results are reported. Comparisons among TMs and other numerical methods are made. It is concluded that the CTM is the simplest algorithm and provides the most accurate solution with the best numerical stability. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

求解P0函数互补问题的一种微分方程方法

本针对P0函数互补问题，给出了一种微方程方法，并且证明了P0函数互补问题的解是微分方程系统的渐进稳定平衡点。在适当的假设条件下，证明了所给出的算法具有二次收敛速度。几个数值例子表明了该算法的有效性。     

Pointwise and L2 bounds in some nonstandard problems for the heat equation

We consider some initial-boundary value problems for the linear and nonlinear heat equation where the gradient of the solution is prescribed on the boundary. Assuming that a solution exists, we obtain bounds for the solution and its gradient by maximum principle arguments or by means of differential and integral inequalities.