求解单侧障碍问题的隐式投影方法

利用不动点原理,得到了求解一类障碍问题的隐式投影算法.采用中心差分格式将障碍问题离散为一个线性互补问题,从而得到了基于投影形式的隐式算法.该方法的每一步迭代只需要求解一个线性方程组.用投影性质很容易证明算法收敛性.给出了具体的算法过程,数值算例结果和理论分析是一致的.     

四阶障碍问题的稳定化混合有限元方法

本文讨论了四阶障碍问题的稳定化混合有限元方法.首先,引入网格依赖范数,通过加罚方法得到了与四阶障碍问题的等价的混合变分形式.随后给出了基于C~0协调有限元空间(W_h,V_h)的混合有限元逼近,例如P_k-P_k三角形有限元.在网格依赖范数下,(W_h,V_h)满足离散的inf-sup条件.最后,我们在不同的假设下,得到了一些误差估计.     

Numerical Approximation of First Eigenvalues for a Heat Conduction Problem

An eigenvalue problem for a differential operator connected with the heat conduction of a non-steady flow is considered. By using an iterative method for computing the eigenvalues of second kind Fredholm operators (see [1]-[11]) we derive approximation of first eigenvalues closer with respect to results appearing in preceding literature.     

A Probabilistic-Numerical Approximation for an Obstacle Problem Arising in Game Theory

We investigate a two-player zero-sum stochastic differential game in which one of the players has more information on the game than his opponent. We show how to construct numerical schemes for the value function of this game, which is given by the solution of a quasilinear partial differential equation with obstacle.     

Error Estimates for the Numerical Approximation of a Semilinear Elliptic Control Problem

We study the numerical approximation of distributed nonlinear optimal control problems governed by semilinear elliptic partial differential equations with pointwise constraints on the control. Our main result are error estimates for optimal controls in the maximum norm. Characterization results are stated for optimal and discretized optimal control. Moreover, the uniform convergence of discretized controls to optimal controls is proven under natural assumptions.     

Numerical Approximation of Solutions of a Nonlinear Inverse Problem Arising in Olfaction Experimentation

Identification of detailed features of neuronal systems is an important challenge in the biosciences today. Cilia are long thin structures that extend from the olfactory receptor neurons into the nasal mucus. Transduction of an odor into an electrical signal occurs in the membranes of the cilia. The cyclic-nucleotide-gated (CNG) channels which reside in the ciliary membrane and are activated by adenosine 3',5'-cyclic monophosphate (cAMP) allow a depolarizing influx of Ca(2+) and Na(+) and are thought to initiate the electrical signal.In this paper, a mathematical model consisting of two nonlinear differential equations and a constrained Fredholm integral equation of the first kind is developed to model experiments involving the diffusion of cAMP into cilia and the resulting electrical activity. The unknowns in the problem are the concentration of cAMP, the membrane potential and, the quantity of most interest in this work, the distribution of CNG channels along the length of a cilium. A simple numerical method is derived that can be used to obtain estimates of the spatial distribution of CNG ion channels along the length of a cilium. Certain computations indicate that this mathematical problem is ill-conditioned.     

Numerical Approximation of the Elastic-Viscoplastic Contact Problem with Non-matching Meshes

Summary.  This work deals with the approximation of a time dependent variational inequality modelling the unilateral contact problem of elastic-viscoplastic bodies in a bidimensional context. The problem is approximated in the space variable with nonconforming finite element methods which allow the handling of nonmatching meshes on the contact zone. Several error estimates are established and the corresponding numerical experiments are achieved. Received Febuary 7, 2001 / Revised version received August 14, 2002 / Published online October 29, 2002 Mathematics Subject Classification (1991): 73T05     

Two-Phase Obstacle Problem

On the basis of the monotonicity formula due to Alt, Caffarelli, and Friedman, the boundedness of the second-order derivatives D ² u of solutions to the equation

Remark on a Problem of Rational Approximation

We show that for any nonincreasing number sequence {a _n} _{n= 0} converging to zero, there exists a continuous 2-periodic function g such that the sequence of its best uniform trigonometric rational approximations {R _n (g,C ₂)} _{n= 0} and the sequence {a _n} _{n= 0} have the same order of decay.     

Finite-Element Approximation of a Plasma Equilibrium Problem

The plasma problem studied is: given R⁺ find (, d, u) R ?R ? H¹() such that Let ₁ < ₂ be the first two eigenvalues of the associatedlinear eigenvalue problem: find $$\left(\lambda ,\phi \right)\in\mathrm{R;}\times {\hbox{ H }}_{0}^{1}\left(\Omega \right)$$such that For ₀(0,₂) it is well known that there exists a unique solution(₀, d₀, u₀) to the above problem. We show that the standard continuous piecewise linear Galerkinfinite-element approximatinon $$\left({\lambda }_{0},{\hbox{d }}_{0}^{k},{u}_{0}^{h}\right)$$, for ₀(0,₂), converges atthe optimal rate in the H¹, L², and L norms as h, the mesh length,tends to 0. In addition, we show that dist (, ^h)Ch² ln 1/h,where $${\Gamma }^{\left(h\right)}=\left\{x\in \Omega :{u}_{0}^{\left(h\right)}\left(x\right)=0\right\}$$.Finally we consider a more practical approximation involvingnumerical integration.     

An Obstacle Control Problem with a Source Term

   Abstract. An optimal control problem for an elliptic variational inequality with a source term is considered. The obstacle is the control, and the goal is to keep the solution of the variational inequality close to the desired profile while the H ¹ norm of the obstacle is not too large. The addition of the source term strongly affects the needed compactness result for the existence of a minimizer.     

An Obstacle Control Problem with a Source Term

Abstract. An optimal control problem for an elliptic variational inequality with a source term is considered. The obstacle is the control, and the goal is to keep the solution of the variational inequality close to the desired profile while the H ¹ norm of the obstacle is not too large. The addition of the source term strongly affects the needed compactness result for the existence of a minimizer.     

Approximation of a Two-phase Continuous Casting Stefan Problem

The continuous casting Stefan problem is a mathematical model describing the solidification with convection of a material being cast continuously with a prescribed velocity. We propose a practical piecewise linear finite element scheme motivated by the characteristic finite element method and derive an error estimate for the scheme which is of the same convergence order as that proved for Stefan problem without convection.     

Optimal Control of the Obstacle Problem in a Perforated Domain

We study the problem of optimally controlling the solution of the obstacle problem in a domain perforated by small periodically distributed holes. The solution is controlled by the choice of a perforated obstacle which is to be chosen in such a fashion that the solution is close to a given profile and the obstacle is not too irregular. We prove existence, uniqueness and stability of an optimal obstacle and derive necessary and sufficient conditions for optimality. When the number of holes increase indefinitely we determine the limit of the sequence of optimal obstacles and solutions. This limit depends strongly on the rate at which the size of the holes shrink.     

Approximation of the Posterior Distribution in a Change-Point Problem

We consider a family of models that arise in connection with sharp change in hazard rate corresponding to high initial hazard rate dropping to a more stable or slowly changing rate at an unknown change-point . Although the Bayes estimates are well behaved and are asymptotically efficient, it is difficult to compute them as the posterior distributions are generally very complicated. We obtain a simple first order asymptotic approximation to the posterior distribution of . The accuracy of the approximation is judged through simulation. The approximation performs quite well. Our method is also applied to analyze a real data set.     

Numerical Approximation of Hysteresis Problems

Hysteresis effects are represented by means of continuous or(multivalued) discontinuous Volterra functionals. Constitutiverelations of this type are coupled with parabolic equationsand existence results are given for the corresponding weak formulations.Implicit time discretization and finite elements are used forthe numerical approximation. Several numerical tests show convergenceof the approximate solutions.     

Finite Element Approximation of a Contact Vector Eigenvalue Problem

We consider a nonstandard elliptic eigenvalue problem of second order on a two-component domain consisting of two intervals with a contact point. The interaction between the two domains is expressed through a coupling condition of nonlocal type, more specifically, in integral form. The problem under consideration is first stated in its variational form and next interpreted as a second-order differential eigenvalue problem. The aim is to set up a finite element method for this problem. The error analysis involved is shown to be affected by the nonlocal condition, which requires a suitable modification of the vector Lagrange interpolant on the overall finite element mesh. Nevertheless, we arrive at optimal error estimates. In the last section, an illustrative numerical example is given, which confirms the theoretical results.