首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 0 毫秒
1.
    
Anisotropic diffusion is important to many different types of common materials and media. Based on structured Cartesian meshes, we develop a three-dimensional (3D) nonhomogeneous immersed finite-element (IFE) method for the interface problem of anisotropic diffusion, which is characterized by an anisotropic elliptic equation with discontinuous tensor coefficient and nonhomogeneous flux jump. We first construct the 3D linear IFE space for the anisotropic nonhomogeneous jump conditions. Then we present the IFE Galerkin method for the anisotropic elliptic equation. Since this method can efficiently solve interface problems on structured Cartesian meshes, it provides a promising tool to solve the physical models with complex geometries of different materials, hence can serve as an efficient field solver in a simulation on Cartesian meshes for related problems, such as the particle-in-cell simulation. Numerical examples are provided to demonstrate the features of the proposed method.  相似文献   

2.
    
We introduce a new method to triangulate planar, curved domains that transforms a specific collection of triangles in a background mesh to conform to the boundary. In the process, no new vertices are introduced, and connectivities of triangles are left unaltered. The method relies on a novel way of parameterizing an immersed boundary over a collection of nearby edges with its closest point projection. To guarantee its robustness, we require that the domain be C2‐regular, the background mesh be sufficiently refined near the boundary, and that specific angles in triangles near the boundary be strictly acute. The method can render both straight‐edged and curvilinear triangulations for the immersed domain. The latter includes curved triangles that conform exactly to the immersed boundary, and ones constructed with isoparametric mappings to interpolate the boundary at select points. High‐order finite elements constructed over these curved triangles achieve optimal accuracy, which has customarily proven difficult in numerical schemes that adopt nonconforming meshes. Aside from serving as a quick and simple tool for meshing planar curved domains with complex shapes, the method provides significant advantages for simulating problems with moving boundaries and in numerical schemes that require iterating over the geometry of domains. With no conformity requirements, the same background mesh can be adopted to triangulate a large family of domains immersed in it, including ones realized over several updates during the coarse of simulating problems with moving boundaries. We term such a background mesh as a universal mesh for the family of domains it can be used to triangulate. Universal meshes hence facilitate a framework for finite element calculations over evolving domains while using only fixed background meshes. Furthermore, because the evolving geometry can be approximated with any desired order, numerical solutions can be computed with high‐order accuracy. We present demonstrative examples using universal meshes to simulate the interaction of rigid bodies with Stokesian fluids. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

3.
    
A numerical method to approximate partial differential equations on meshes that do not conform to the domain boundaries is introduced. The proposed method is conceptually simple and free of user‐defined parameters. Starting with a conforming finite element mesh, the key ingredient is to switch those elements intersected by the Dirichlet boundary to a discontinuous‐Galerkin approximation and impose the Dirichlet boundary conditions strongly. By virtue of relaxing the continuity constraint at those elements, boundary locking is avoided and optimal‐order convergence is achieved. This is shown through numerical experiments in reaction–diffusion problems. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

4.
    
We discuss the discretization using discontinuous Galerkin (DG) formulation of an elliptic Poisson problem. Two commonly used DG schemes are investigated: the original average flux proposed by Bassi and Rebay (J. Comput. Phys. 1997; 131 :267) and the local discontinuous Galerkin (LDG) (SIAM J. Numer. Anal. 1998; 35 :2440–2463) scheme. In this paper we expand on previous expositions (Discontinuous Galerkin Methods: Theory, Computation and Applications. Springer: Berlin, 2000; 135–146; SIAM J. Sci. Comput. 2002; 24 (2):524–547; Int. J. Numer. Meth. Engng. 2003; 58 (2): 1119–1148) by adopting a matrix based notation with a view to highlighting the steps required in a numerical implementation of the DG method. Through consideration of standard C0‐type expansion bases, as opposed to elementally orthogonal expansions, with the matrix formulation we are able to apply static condensation techniques to improve efficiency of the direct solver when high order expansions are adopted. The use of C0‐type expansions also permits the direct enforcement of Dirichlet boundary conditions through a ‘lifting’ approach where the LDG flux does not require further stabilization. In our construction we also adopt a formulation of the continuous DG fluxes that permits a more general interpretation of their numerical implementation. In particular it allows us to determine the conditions under which the LDG method provides a near local stencil. Finally a study of the conditioning and the size of the null space of the matrix systems resulting from the DG discretization of the elliptic problem is undertaken. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

5.
    
The superconvergent patch recovery technique with constraints (SPR-C) consists in improving the accuracy of the recovered stresses obtained with the original SPR technique by considering known information about the exact solution, like the internal equilibrium equation, the compatibility equation or the Neumann boundary conditions, during the recovery process. In this paper the SPR-C is extended to consider the equilibrium around the contact area when solving contact problems with the Cartesian grid Finite Element Method. In the proposed method, the Finite Element stress fields of both bodies in contact are considered during the recovery process and the equilibrium is enforced by means of the continuity of tractions along the contact surface.  相似文献   

6.
    
Large isoparametric macroelements with closed‐form cardinal global shape functions under the label ‘Coons‐patch macroelements’ (CPM) have been previously proposed and used in conjunction with the finite element method and the boundary element method. This paper continues the research on the performance of CPM in conjunction with the collocation method. In contrast to the previous CPM that was based on a Galerkin/Ritz formulation, no domain integration is now required, a fact that justifies the name ‘integration‐free Coons macroelements’. Therefore, in addition to avoiding mesh generation, and saving human effort, the proposed technique has the additional advantage of further reducing the computer effort. The theory is supported by five test cases concerning Poisson and Laplace problems within 2D smooth quadrilateral domains. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

7.
建立了边值问题｛div(|Du|^p-2Du) a(|x|f(|x|,u)=0,在B属于R^Nδu／δ λu=-α 在δB｝正对称解的存在性和唯一性,这里B是R^N中的以原点为圆心的单位图。  相似文献   

8.
    
The extended finite element method (XFEM) enables the representation of cracks in arbitrary locations of a mesh. We introduce here a variant of the XFEM rendering an optimally convergent scheme. Its distinguishing features are as follows: (a) the introduction of singular asymptotic crack tip fields with support on only a small region around the crack tip (the enrichment region), (b) only one and two enrichment functions are added for anti‐plane shear and planar problems, respectively and (c) the relaxation of the continuity between the enrichment region and the rest of the domain, and the adoption of a discontinuous Galerkin (DG) method therein. The method is provably stable for any positive value of a stabilization parameter, and by weakly enforcing the continuity between the two regions it eliminates ‘blending elements’ partly responsible for the suboptimal convergence of some early XFEMs. Moreover, the particular choice of enrichment functions results in a surprisingly sparse stiffness matrix that remains reasonably conditioned as the mesh is refined. More importantly, the stress intensity factors can be extracted with a satisfactory accuracy as primary unknowns. Quadrature strategies required for the optimal convergence are also discussed. Finally, the DG method was modified to retain stability based on an inf‐sup condition. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

9.
针对由不同特性物体所组成的多体系统,探讨了能够涵盖各种变形状态和运动形式的广义有限单元模式及其插值函数形式.对于多体接触问题,发展了能够合理描述界面特性的接触力元模型,即采用某种应力插值函数将界面上的相互作用力由接触对上的接触应力来表达,并将接触对上的接触应力当作需满足界面上屈服准则与流动法则等状态控制条件的参变量,将其作为约束条件加入系统控制方程.根据非连续变形系统的分区参变量最小势能变分原理,联立变分驻值条件与参变量的状态控制条件建立了多体系统非连续变形计算力学分析的基本控制方程,将问题最终归结为一个含有自由变量和等式约束条件的线性互补问题,对此发展了数值解法,并进行了多个算例的数值分析.计算结果表明该模型不仅能够对多体系统进行静、动力耦合分析,而且还能够模拟多体系统的变形与应力及接触界面上的接触应力和相对运动等复杂的非线性过程.  相似文献   

10.
    
We construct a method for the parameterization of a class of planar piecewise C2‐curves over a collection of edges in an ambient triangulation. The map from the collection of edges to the curve is the closest‐point projection. A distinguishing feature of the method is that edges in the ambient triangulation need not interpolate the curve. We formulate conditions on the ambient triangulations so that the resulting parameterization over its selected edges is (i) bijective, (ii) maps simple, connected collection of edges to simple, connected components of the curve, and (iii) is C1 within each edge of the collection. These properties of the parameterization make it particularly useful in the construction of high‐order finite element approximation spaces on planar curves immersed in triangulations. We discuss this application and illustrate it with numerical examples. The parameterization method applies to a large class of planar curves, including most ones of interest in engineering and computer graphics applications, and to a large family of triangulations, including acute‐angled triangulations. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

11.
    
We describe an algorithm to recover a boundary-fitting triangulation for a bounded C2-regular domain immersed in a nonconforming background mesh of tetrahedra. The algorithm consists in identifying a polyhedral domain ωh bounded by facets in the background mesh and morphing ωh into a boundary-fitting polyhedral approximation Ωh of Ω. We discuss assumptions on the regularity of the domain, on element sizes and on specific angles in the background mesh that appear to render the algorithm robust. With the distinctive feature of involving just small perturbations of a few elements of the background mesh that are in the vicinity of the immersed boundary, the algorithm is designed to benefit numerical schemes for simulating free and moving boundary problems. In such problems, it is now possible to immerse an evolving geometry in the same background mesh, called a universal mesh, and recover conforming discretizations for it. In particular, the algorithm entirely avoids remeshing-type operations and its complexity scales approximately linearly with the number of elements in the vicinity of the immersed boundary. We include detailed examples examining its performance.  相似文献   

12.
A methodology to connect two or more non-conforming meshes generated from a geometric model is presented. The problem leads to a variational problem subject to constraints. The emphasis of the work is on constraints construction. Two methods, namely, discrete connection and continuous connection, are developed. The relation with hybrid formulations is considered and two application examples are presented.  相似文献   

13.
    
The paper describes the main features of an automatic and three‐dimensional Cartesian mesher specifically designed for compressible inviscid/viscous flow solvers based on an immersed boundary technique. The development of a meshing tool able of conducting non‐isotropic cell refinements is a very tiresome task. The major difficulty is to imagine, at the pre‐design phase, a light but flexible data management, which minimizes the memory and CPUs' resources. In particular, the embedded geometry has to be detected by means of a fast and robust tagging procedure. Cells in proximity of the wall have to be refined in a proper way to adequately solve large flow gradients. Smooth variation of mesh density between differently refined zones has to be guaranteed to increase the flow solver robustness. A procedure to obtain accurate data on the geometry surfaces should be foresee. Here, a robust algorithm is developed to reconstruct a surface triangulation starting from the intersection points among volume cells and the geometry surfaces. The paper attempts to address all the above issues to help the readers in designing their own tools and suggesting them a way forward.  相似文献   

14.
    
A new class of finite elements is described for dealing with non-matching meshes, for which the existing finite elements are hardly efficient. The approach is to employ the moving least-square (MLS) scheme to devise a class of elements with an arbitrary number of nodal points on the parental domain. This approach generally leads to elements with the rational shape functions, which significantly extends the function space of the conventional finite element method. With a special choice of the nodal points and the base functions, the method results in useful elements with the polynomial shape functions for which the C1 continuity breaks down across the boundaries between the subdomains comprising one element. The present scheme possesses an extremely high potential for applications which deal with various problems with discontinuities, such as material inhomogeneity, crack propagation, phase transition and contact mechanics. The effectiveness of the new elements for handling the discontinuities due to non-matching interfaces is demonstrated using appropriate examples. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

15.
自适应笛卡尔网格湍流模拟时,网格数目呈雷诺数的指数方量级增长,带来网格规模“灾难”问题。壁面函数方法可以将壁面第一层网格尺度放宽1~2个量级,是有效的解决办法之一,然而,由于笛卡尔网格的非贴体特性,壁面函数的实现及其与不同湍流模型的匹配是关键。针对这一问题,从经典的Spalding壁面函数出发,采用显式赋初值的方式,解决了壁面摩擦速度迭代求解初值敏感及分层壁面函数速度不光滑问题;基于切应力恒定的壁面条件,构造了多层虚拟网格的切向速度和湍流粘性系数,实现了笛卡尔网格下的壁面函数方法;采用结构系综理论,通过DNS数据标定,构造了湍动能随湍流粘性比率变化的关系式,实现了壁面函数与KDO湍流模型的匹配,并以平板T3b、湍流凸包和跨音速NACA0012翼型绕流为算例,与常用的SST、SA湍流模型壁面函数方法进行了对比分析。结果表明：该文提出的KDO湍流模型壁面函数法准确、可靠,且因其形式简单、计算量相对较小,有较大的推广应用价值。  相似文献   

16.
We apply a combination of the transient scaled boundary finite element method (SBFEM) and quadtree‐based discretization to model dynamic problems at high frequencies. We demonstrate that the current formulation of the SBFEM for dynamics tends to require more degrees of freedom than a corresponding spectral element discretization when dealing with smooth problems on regular domains. Thus, we improve the efficiency of the SBFEM by proposing a novel approach to reduce the number of auxiliary variables for transient analyses. Based on this improved SBFEM, we present a modified meshing procedure, which creates a quadtree mesh purely based on the geometry and allows arbitrary sizes and orders of elements, as well as an arbitrary number of different materials. The discretization of each subdomain is created automatically based on material parameters and the highest frequency of interest. The transition between regions of different properties is straightforward when using the SBFEM. The proposed approach is applied to image‐based analysis with a particular focus on geological models. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

17.
    
A paradigm is developed for generating structured finite element models from solid models by means of implicit surface definitions. The implicit surfaces are defined by radial basis functions. Internal features, such as material interfaces, sliding interfaces and cracks are treated by enrichment techniques developed in the extended finite element method. Methods for integrating the weak form for such models are proposed. These methods simplify the generation of finite element models. Results presented for several examples show that the accuracy of this method is comparable to standard unstructured finite element methods. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

18.
    
Two methods are presented for connecting dissimilar three‐dimensional finite element meshes. The first method combines the concept of master and slave surfaces with the uniform strain approach for finite elements. By modifying the boundaries of elements on a slave surface, corrections are made to element formulations such that first‐order patch tests are passed. The second method is based entirely on constraint equations, but only passes a weaker form of the patch test for non‐planar surfaces. Both methods can be used to connect meshes with different element types. In addition, master and slave surfaces can be designated independently of relative mesh resolutions. Example problems in three‐dimensional linear elasticity are presented. Published in 2000 by John Wiley & Sons, Ltd.  相似文献   

19.
    
This paper gives an explicit geometric interpretation of finite element formulations of the eddy‐current problem. The paper shows, step‐by‐step how the eddy‐current problem can be implemented in a finite element kind of software system exploiting familiar geometric ideas such as lengths of edges, areas of faces, volumes of tetrahedra, and the mutual interconnections between the simplices of a mesh. The approach is a specific case of the so‐called geometric techniques. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

20.
    
  相似文献   

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

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