首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
    
This paper considers the numerical simulation of optimal control evolution dam problem by using conjugate gradient method.The paper considers the free boundary value problem related to time dependent fluid flow in a homogeneous earth rectangular dam.The dam is taken to be sufficiently long that the flow is considered to be two dimensional.On the left and right walls of the dam there is a reservoir of fluid at a level dependent on time.This problem can be transformed into a variational inequality on a fixed domain.The numerical techniques we use are based on a linear finite element method to approximate the state equations and a conjugate gradient algorithm to solve the discrete optimal control problem.This algorithm is based on Armijo's rule in the unconstrained optimization theory.The convergence of the discrete optimal solutions to the continuous optimal solutions,and the convergence of the conjugate gradient algorithm are proved.A numerical example is given to determine the location of the minimum surface  相似文献   

2.
1.IntroductionTherehavebeennumerousworkintheanalysisoffiniteelementmethodsfortheunilateralproblem(c.f.[4]andthereferencestherein).ItshouldbementionedthatinF.Scarpiniet.al.[6],I.Hlavaceket.al.[5]andF.Brezziet.al.[1],theconforminglinearelementapproximationtotheunilateralproblemhavebeenconsidered,andthevariouserrorboundshavebeenobtainedinthedifferentassumptionofregularityofsolutionoftheproblem.Inthispaper3weconsiderthreenonconformingfiniteelement(i.e.twoCrouzrixRaviartlinearelementsandWilsone…  相似文献   

3.
THE NONCONFORMING FINITE ELEMENT METHOD FOR SIGNORINI PROBLEM   总被引:1,自引:0,他引:1  
We present the Crouzeix-Raviart linear nonconforming finite element approximation of the variational inequality resulting from Signorini problem. We show if the displacement field is of H2 regularity, then the convergence rate can be improved from O(h3/4) to quasi-optimal O(h|log h|1/4) with respect to the energy norm as that of the continuous linear finite element approximation. If stronger but reasonable regularity is available, the convergence rate can be improved to the optimal O(h) as expected by the linear approximation.  相似文献   

4.
Regular assumption of finite element meshes is a basic condition of most analysis offinite element approximations both for conventional conforming elements and nonconform-ing elements.The aim of this paper is to present a novel approach of dealing with theapproximation of a four-degree nonconforming finite element for the second order ellipticproblems on the anisotropic meshes.The optimal error estimates of energy norm and L~2-norm without the regular assumption or quasi-uniform assumption are obtained based onsome new special features of this element discovered herein.Numerical results are givento demonstrate validity of our theoretical analysis.  相似文献   

5.
For a second-order elliptic boundary value problem, We develop an intergrid transfer operator in multigrid method for the P1-nonconforming finite element method. This intergrid transfer operator needs smaller computation than previous intergrid transfer operators. Multigrid method with this operator converges well.  相似文献   

6.
This paper is devoted to analysis of the nonconforming element approximation to the obstacle problem, and improvement and correction of the results in [11], [12].  相似文献   

7.
Bergouniou  Maïtine  Lenhart  Suzanne 《Positivity》2004,8(3):229-242
We consider an optimal control problem where the state satisfies an obstacle type semilinear variational inequality and the control function is the obstacle. The state is chosen to be close to a desired profile while the obstacle is not too large in H 0 1 (), and H 2-bounded. We prove that an optimal control exists and give necessary optimality conditions, using approximation techniques.  相似文献   

8.
In this paper, we present two nonconforming finite elements for the pure displacement planar elasticity problem. Both of them are locking-free and have two order of convergence. Some numerical results attest the validity of our theoretical analysis.  相似文献   

9.
General optimality conditions are obtained for optimal shape design for systems governed by a class of elliptic variational inequalities. The conditions are established by making use of the necessary conditions for optimal control of systems governed by strongly monotone variational inequalities. These conditions are then applied to an electrochemical machining problem.  相似文献   

10.
    
We introduce a piecewise P2-nonconforming quadrilateral finite element. First, we decompose a convex quadrilateral into the union of four triangles divided by its diagonals. Then the finite element space is defined by the set of all piecewise P2-polynomials that are quadratic in each triangle and continuously differentiable on the quadrilateral. The degrees of freedom (DOFs) are defined by the eight values at the two Gauss points on each of the four edges plus the value at the intersection of the diagonals. Due to the existence of one linear relation among the above DOFs, it turns out the DOFs are eight. Global basis functions are defined in three types: vertex-wise, edge-wise, and element-wise types. The corresponding dimensions are counted for both Dirichlet and Neumann types of elliptic problems. For second-order elliptic problems and the Stokes problem, the local and global interpolation operators are defined. Also error estimates of optimal order are given in both broken energy and L2(Ω) norms. The proposed element is also suitable to solve Stokes equations. The element is applied to approximate each component of velocity fields while the discontinuous P1-nonconforming quadrilateral element is adopted to approximate the pressure. An optimal error estimate in energy norm is derived. Numerical results are shown to confirm the optimality of the presented piecewise P2-nonconforming element on quadrilaterals.https://doi.org/10.1051/m2an/2012044  相似文献   

11.
TRAPEZOIDAL PLATE BENDING ELEMENT WITH DOUBLE SET PARAMETERS   总被引:1,自引:0,他引:1  
Using double set parameter method, a 12-parameter trapezoidal plate bending element is presented. The first set of degrees of freedom, which make the element convergent, are the values at the four vertices and the middle points of the four sides together with the mean values of the outer normal derivatives along four sides. The second set of degree of freedom, which make the number of unknowns in the resulting discrete system small and computation convenient are values and the first derivatives at the four vertices of the element. The convergence of the element is proved.  相似文献   

12.
    
The mathematical model of multilayer viscoelastic system based on the actual structure of asphalt pavement and the properties of its numerical solutions are studied. Firstly, based on the partial differential equations, the variational inequality model of the multilayer viscoelastic system is derived, and the existence and uniqueness of its solutions are further proved. Then, based on the finite element theory, the convergence property and error analysis of the semi-discrete numerical solutions of the variational inequalities are further studied. Finally, the convergence and error analysis of the fully discrete numerical solutions of the variational inequalities are derived by the forward difference method. The conclusion of the above error analysis also confirms the feasibility of studying the mechanical response of asphalt pavement based on variational inequality.  相似文献   

13.
SOME n-RECTANGLE NONCONFORMING ELEMENTS FOR FOURTH ORDER ELLIPTIC EQUATIONS   总被引:1,自引:0,他引:1  
In this paper, three n-rectangle nonconforming elements are proposed with n ≥ 3. They are the extensions of well-known Morley element, Adini element and Bogner-Fox-Schmit element in two spatial dimensions to any higher dimensions respectively. These elements are all proved to be convergent for a model biharmonic equation in n dimensions.  相似文献   

14.
In this paper, the well-known nonconforming Morley element for biharmonic equations in two spatial dimensions is extended to any higher dimensions in a canonical fashion. The general n-dimensional Morley element consists of all quadratic polynomials defined on each n-simplex with degrees of freedom given by the integral average of the normal derivative on each (n-1)-subsimplex and the integral average of the function value on each (n-2)-subsimplex. Explicit expressions of nodal basis functions are also obtained for this element on general n-simplicial grids. Convergence analysis is given for this element when it is applied as a nonconforming finite element discretization for the biharmonic equation. The work was supported by the National Natural Science Foundation of China (10571006). This work was supported in part by NSF DMS-0209497 and NSF DMS-0215392 and the Changjiang Professorship through Peking University.  相似文献   

15.
    
A low order nonconforming mixed finite element method (FEM) is established for the fully coupled non-stationary incompressible magnetohydrodynamics (MHD) problem in a bounded domain in 3D. The lowest order finite elements on tetrahedra or hexahedra are chosen to approximate the pressure, the velocity field and the magnetic field, in which the hydrodynamic unknowns are approximated by inf-sup stable finite element pairs and the magnetic field by $H^1(Omega)$-conforming finite elements, respectively. The existence and uniqueness of the approximate solutions are shown. Optimal order error estimates of $L^2(H^1)$-norm for the velocity field, $L^2(L^2)$-norm for the pressure and the broken $L^2(H^1)$-norm for the magnetic field are derived.  相似文献   

16.
We consider a material that occupies a convex polygonal bounded domain Ω ⊂ ℝn, with regular boundary Γ = Γ1 ∪ Γ2 (with Γ ∩ Γ = ∅︁) with meas (Γ1) = |Γ1| > 0 and |Γ2| > 0. We assume, without loss of generality, that the melting temperature is 0°C. We consider the following steady‐state heat conduction problem in Ω: with α, q, B = Const > 0, and q and α represent the heat flux on Γ2 and the heat transfer coefficient on Γ1, respectively. In a previous article (Tabacman‐ Tarzia, J Diff Eq 77 (1989), 16– 37) sufficient and/or necessary conditions on data α, q, B, Ω, Γ1, Γ2 to obtain a temperature u of nonconstant sign in Ω (that is, a multidimensional steady‐state, two‐phase, Stefan problem) were studied. In this article, we consider a regular triangulation by finite element method of the domain Ω with Lagrange triangles of the type 1, with h > 0 the parameter of the discretization. We study sufficient (and/or necessary) conditions on data α, q, B, Ω, Γ1, and Γ2 to obtain a change of phase (steady‐state, two‐phase, discretized Stefan problem) in corresponding discretized domain, that is, a discrete temperature of nonconstant sign in Ω. Moreover, error bounds as a function of the parameter h, are also obtained. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq. 15: 355–369, 1999  相似文献   

17.
18.
In this paper, using a bubble function, we construct a cuboid element to solve the fourth order elliptic singular perturbation problem in three dimensions. We prove that the nonconforming CO-cuboid element converges in the energy norm uniformly with respect to the perturbation parameter.  相似文献   

19.
We prove an existence theorem for the optimal control of variational inequalities governed by a pseudomonotone operator: the cost is assumed to be quadratic. Then, we give an extension of the theorem to more general costs (assuming the operator to be monotone); we also give a result on a perturbation problem.This work is an extended part of the author's thesis, written under the direction of Professor T. Zolezzi. This research was partially supported by the Consiglio Nazionale delle Ricerche (CNR), Rome, Italy.  相似文献   

20.
In this paper, two new nonconforming hexagonal elements are presented, which are based on the trilinear function space Q1^(3) and are edge-oriented, analogical to the case of the rotated Q1 quadrilateral element. A priori error estimates are given to show that the new elements achieve first-order accuracy in the energy norm and second-order accuracy in the L^2 norm. This theoretical result is confirmed by the numerical tests.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号