In this paper, we derive a posteriori error estimates for the finiteelementapproximation of quadratic optimal control problem governed by linear parabolic equation. We obtain a posteriori error estimates for both th...
详细信息
In this paper, we derive a posteriori error estimates for the finiteelementapproximation of quadratic optimal control problem governed by linear parabolic equation. We obtain a posteriori error estimates for both the state and the control approximation. Such estimates, which are apparently not available in the literature, are an important step towards developing reliable adaptive finiteelementapproximation schemes for the control problem.
A finite-element approximation of a fourth-order differential equation is given. In the direct implementation, a system of nonlinear equations will be obtained and also a full size matrix will be introduced when Newto...
详细信息
A finite-element approximation of a fourth-order differential equation is given. In the direct implementation, a system of nonlinear equations will be obtained and also a full size matrix will be introduced when Newton's method is adopted to solve the system. This difficulty can be avoided by the use of an iterative scheme which is shown to converge to the solution.
In this paper we consider the numerical approximation of phase field and sharp interface models for diffusion-induced grain boundary motion. The phase field model consists of a double-obstacle Allen-Cahn equation with...
详细信息
In this paper we consider the numerical approximation of phase field and sharp interface models for diffusion-induced grain boundary motion. The phase field model consists of a double-obstacle Allen-Cahn equation with a forcing obtained from the solution of a degenerate diffusion equation. On the other hand the sharp interface model consists of forced mean curvature flow coupled to a diffusion equation holding on the interface itself. Formal asymptotics yield the sharp interface model as the limit of the phase field equations as the width of the associated diffuse interface tends to zero. A finite-element approximation of the phase field model is presented and is shown to be convergent to a weak solution. Numerical simulations of both models are described and compared. It is shown that the two models are consistent.
In non-convex optimisation problems, in particular in non-convex variational problems, there usually does not exist any classical solution but only generalised solutions which involve Young measures. In this paper, fi...
详细信息
In non-convex optimisation problems, in particular in non-convex variational problems, there usually does not exist any classical solution but only generalised solutions which involve Young measures. In this paper, first a suitable relaxation and approximation theory is developed together with optimality conditions, and then an adaptive scheme is proposed for the efficient numerical treatment. The Young measures solving the approximate problems are usually composed only from a few atoms. This is the main argument our effective active-set type algorithm is based on. The support of those atoms is estimated from the Weierstrass maximum principle which involves a Hamiltonian whose good guess is obtained by a multilevel technique. Numerical experiments are performed in a one-dimensional variational problem and support efficiency of the algorithm.
We present a convergence analysis of an algorithm for the numerical computation of the rank-one convex envelope of a function f : M-mxn --> R. A rate of convergence for the scheme is established, and numerical expe...
详细信息
We present a convergence analysis of an algorithm for the numerical computation of the rank-one convex envelope of a function f : M-mxn --> R. A rate of convergence for the scheme is established, and numerical experiments are presented to illustrate the analytical results and applications of the algorithm.
This paper deals with optimal control problems associated with the 2-D Boussinesq equations. The controls considered may be of either the distributed or the Neumann type. These problems are first put into an appropria...
详细信息
This paper deals with optimal control problems associated with the 2-D Boussinesq equations. The controls considered may be of either the distributed or the Neumann type. These problems are first put into an appropriate mathematical formulation. Then the existence of optimal solutions is proved. The use of Lagrange multiplier techniques is justified and an optimality system of equations is derived. (C) 200 Academic Press.
We investigate a general control problem for a class of nonlinear parabolic evolution equations. Applications are related to solid-solid and solid-liquid phase transitions. We prove compactness of the solution operato...
详细信息
We investigate a general control problem for a class of nonlinear parabolic evolution equations. Applications are related to solid-solid and solid-liquid phase transitions. We prove compactness of the solution operator, existence of optimal controls and show convergence of the finite-dimensional approximate control problem to the original one.
We derive quantitative a posteriori estimates for the error caused by replacing an obstacle problem with a linear problem. The error bound depends on the solution of the linear problem, but is independent of the solut...
详细信息
We derive quantitative a posteriori estimates for the error caused by replacing an obstacle problem with a linear problem. The error bound depends on the solution of the linear problem, but is independent of the solution of the obstacle problem. We then give a discrete version of the a posteriori error estimates, which is used in solving a finite-element system of the obstacle problem. A detailed analysis for a one-dimensional example is given, showing the effectiveness of our error estimates.
An algorithm for thermophysical property calculations based on the finiteelementapproximation in the thermodynamic plane has been suggested. As an example of using the algorithm, errors in obtaining the compressibili...
详细信息
An algorithm for thermophysical property calculations based on the finiteelementapproximation in the thermodynamic plane has been suggested. As an example of using the algorithm, errors in obtaining the compressibility of ethane and its derivatives with respect to reduced density and temperature are discussed.
A scheme for the reaction of hemoglobin with ligands is described, which postulates the functional heterogeneity of the chains, considers all possible combinations of the distribution of the ligand on the four chains ...
详细信息
A scheme for the reaction of hemoglobin with ligands is described, which postulates the functional heterogeneity of the chains, considers all possible combinations of the distribution of the ligand on the four chains of hemoglobin, and does not require simplifying assumptions about the hemoglobin reactivity. Ten tetrameric species are considered, together with 16 reactions between these species, each with an “on” and an “off” rate constant. The dissociation of hemoglobin tetramers into dimers is also considered, with four “on” and four “off” rate constants for the reactions between dimers, and ten equilibrium constants for the reactions between tetramers and dimers. Moreover, some side reactions, such as the “trapping” of ligands by a hemoglobin competitor, are included. A FORTRAN program, suitable for microcomputers, is described for handling this scheme, with some examples showing its advantages.
暂无评论