首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到10条相似文献,搜索用时 218 毫秒
1.
Boundary conforming Delaunay mesh generation   总被引:3,自引:0,他引:3  
A boundary conforming Delaunay mesh is a partitioning of a polyhedral domain into Delaunay simplices such that all boundary simplices satisfy the generalized Gabriel property. It’s dual is a Voronoi partition of the same domain which is preferable for Voronoi-box based finite volume schemes. For arbitrary 2D polygonal regions, such meshes can be generated in optimal time and size. For arbitrary 3D polyhedral domains, however, this problem remains a challenge. The main contribution of this paper is to show that boundary conforming Delaunay meshes for 3D polyhedral domains can be generated efficiently when the smallest input angle of the domain is bounded by arccos 1/3 ≈ 70.53°. In addition, well-shaped tetrahedra and an appropriate mesh size can be obtained. Our new results are achieved by reanalyzing a classical Delaunay refinement algorithm. Note that our theoretical guarantee on the input angle (70.53°) is still too strong for many practical situations. We further discuss variants of the algorithm to relax the input angle restriction and to improve the mesh quality.  相似文献   

2.
In this article we construct and analyze a mixed finite volume method for second‐order nonlinear elliptic problems employing H(div; Ω)‐conforming approximations for the vector variable and completely discontinuous approximations for the scalar variable. The main attractive feature of our method is that, although the vector variable is H(div; Ω)‐conforming, one can eliminate it in a local manner to obtain a discontinuous Galerkin method for the scalar variable. Optimal error estimates will be established for both vector and scalar variables. We also present a fully discrete version of this method that is more convenient for computational purposes. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005  相似文献   

3.
For any 2D triangulation τ, the 1-skeleton mesh of τ is the wireframe mesh defined by the edges of τ, while that for any 3D triangulation τ, the 1-skeleton and the 2-skeleton meshes, respectively, correspond to the wireframe mesh formed by the edges of τ and the “surface” mesh defined by the triangular faces of τ. A skeleton-regular partition of a triangle or a tetrahedra, is a partition that globally applied over each element of a conforming mesh (where the intersection of adjacent elements is a vertex or a common face, or a common edge) produce both a refined conforming mesh and refined and conforming skeleton meshes. Such a partition divides all the edges (and all the faces) of an individual element in the same number of edges (faces). We prove that sequences of meshes constructed by applying a skeleton-regular partition over each element of the preceding mesh have an associated set of difference equations which relate the number of elements, faces, edges and vertices of the nth and (n−1)th meshes. By using these constitutive difference equations we prove that asymptotically the average number of adjacencies over these meshes (number of triangles by node and number of tetrahedra by vertex) is constant when n goes to infinity. We relate these results with the non-degeneracy properties of longest-edge based partitions in 2D and include empirical results which support the conjecture that analogous results hold in 3D.  相似文献   

4.
In this paper, we present pointwise estimates of the streamline diffusion finite element method (SDFEM) for conforming piecewise linears on Shishkin triangular meshes. The method is applied to a model singularly perturbed convection-diffusion problem with characteristic layers. Using a new variant of artificial crosswind diffusion, we prove that uniformly pointwise error bounds away from the layers are of order almost 7/4 (up to a logarithmic factor). In some cases, the convergence order is almost 15/8. Our analysis depends on discrete Green’s functions and sharp estimates of the diffusion and convection parts in the bilinear form. Finally, numerical experiments support our theoretical results.  相似文献   

5.
In this article, we study the superconvergence analysis of conforming bilinear finite element method (FEM) for nonlinear Joule heating equations. Based on the rigorous estimates together with high accuracy analysis of this element, mean value technique and interpolation postprocessing approach, the superclose and superconvergent estimates about the related variables in H1‐norm are derived for semidiscrete and a linearized backward Euler fully discrete schemes, which extends the results of optimal estimates obtained for conforming FEMs in the previous literature. At last, a numerical experiment is performed to verify the theoretical analysis.  相似文献   

6.
Compared with standard Galerkin finite element methods, mixed methods for second‐order elliptic problems give readily available flux approximation, but in general at the expense of having to deal with a more complicated discrete system. This is especially true when conforming elements are involved. Hence it is advantageous to consider a direct method when finding fluxes is just a small part of the overall modeling processes. The purpose of this article is to introduce a direct method combining the standard Galerkin Q1 conforming method with a cheap local flux recovery formula. The approximate flux resides in the lowest order Raviart‐Thomas space and retains local conservation property at the cluster level. A cluster is made up of at most four quadrilaterals. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 104–127, 2004  相似文献   

7.
This paper presents an a posteriori error analysis for the stationary Stokes–Darcy coupled problem approximated by finite element methods on anisotropic meshes in or 3. Korn's inequality for piecewise linear vector fields on anisotropic meshes is established and is applied to non‐conforming finite element method. Then the existence and uniqueness of the approximation solution are deduced for non‐conforming case. With the obtained finite element solutions, the error estimators are constructed and based on the residual of model equations plus the stabilization terms. The lower error bound is proved by means of bubble functions and the corresponding anisotropic inverse inequalities. In order to prove the upper error bound, it is vital that an anisotropic mesh corresponds to the anisotropic function under consideration. To measure this correspondence, a so‐called matching function is defined, and its discussion shows it to be useful tool. With its help, the upper error bound is shown by means of the corresponding anisotropic interpolation estimates and a special Helmholtz decomposition in both media. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

8.
We consider a regularization method for nonlinear complementarity problems with F being a P0-function which replaces the original problem with a sequence of the regularized complementarity problems. In this paper, this sequence of regularized complementarity problems are solved approximately by applying the generalized Newton method for an equivalent augmented system of equations, constructed by the generalized Fischer–Burmeister (FB) NCP-functions φp with p>1. We test the performance of the regularization semismooth Newton method based on the family of NCP-functions through solving all test problems from MCPLIB. Numerical experiments indicate that the method associated with a smaller p, for example p[1.1,2], usually has better numerical performance, and the generalized FB functions φp with p[1.1,2) can be used as the substitutions for the FB function φ2.  相似文献   

9.
In this paper, we derive robust a posteriori error estimates for conforming approximations to a singularly perturbed reaction-diffusion problem on anisotropic meshes, since the solution in general exhibits anisotropic features, e.g., strong boundary or interior layers. Based on the anisotropy of the mesh elements, we improve the a posteriori error estimates developed by Cheddadi et al., which are reliable and efficient on isotropic meshes but fail on anisotropic ones. Without the assumption that the mesh is shape-regular, the resulting mesh-dependent error estimator is shown to be reliable, efficient and robust with respect to the reaction coefficient, as long as the anisotropic mesh sufficiently reflects the anisotropy of the solution. We present our results in the framework of the vertex-centered finite volume method but their nature is general for any conforming one, like the piecewise linear finite element one. Our estimates are based on the usual H(div)-conforming, locally conservative flux reconstruction in the lowest-order Raviart-Thomas space on a dual mesh associated with the original anisotropic simplex one. Numerical experiments in 2D confirm that our estimates are reliable, efficient and robust on anisotropic meshes.  相似文献   

10.
Summary. We consider the approximation of the vibration modes of an elastic plate in contact with a compressible fluid. The plate is modelled by Reissner-Mindlin equations while the fluid is described in terms of displacement variables. This formulation leads to a symmetric eigenvalue problem. Reissner-Mindlin equations are discretized by a mixed method, the equations for the fluid with Raviart-Thomas elements and a non conforming coupling is used on the interface. In order to prove that the method is locking free we consider a family of problems, one for each thickness , and introduce appropriate scalings for the physical parameters so that these problems attain a limit when . We prove that spurious eigenvalues do not arise with this discretization and we obtain optimal order error estimates for the eigenvalues and eigenvectors valid uniformly on the thickness parameter t. Finally we present numerical results confirming the good performance of the method. Received February 4, 1998 / Revised version received May 26, 1999 / Published online June 21, 2000  相似文献   

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

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