Adaptive Mesh Refinement FEM for Damage Evolution of Heterogeneous Brittle Media

Damage evolution of heterogeneous brittle media involves a wide range of length scales. The coupling between these length scales underlies the mechanism of damage evolution and rupture. However, few of previous numerical algorithms consider the effects of the trans-scale coupling effectively. In this paper, an adaptive mesh refinement FEM algorithm is developed to simulate this trans-scale coupling. The adaptive serendipity element is implemented in this algorithm, and several special discontinuous base functions are created to avoid the incompatible displacement between the elements. Both the benchmark and a typical numerical example under quasi-static loading are given to justify the effectiveness of this model. The numerical results reproduce a series of characteristics of damage and rupture in heterogeneous brittle media.

A Numerical Study for Turbulent Flow and Thermal Influence Over Inhomogenous Canopy of Roughness Elements

A large-eddy simulation with transitional structure function(TSF) subgrid model we previously proposed was performed to investigate the turbulent flow with thermal influence over an inhomogeneous canopy, which was represented as alternative large and small roughness elements. The aerodynamic and thermodynamic effects of the presence of a layer of large roughness elements were modelled by adding a drag term to the three-dimensional Navier-Stokes equations and a heat source/sink term to the scalar equation, respectively. The layer of small roughness elements was simply treated using the method as described in paper (Moeng 1984, J. Atmos Sci. 41, 2052-2062) for homogeneous rough surface. The horizontally averaged statistics such as mean vertical profiles of wind velocity, air temperature, et al., are in reasonable agreement with Gao et al.(1989, Boundary layer meteorol. 47, 349-377) field observation (homogeneous canopy). Not surprisingly, the calculated instantaneous velocity and temperature fields show that the roughness elements considerably changed the turbulent structure within the canopy. The adjustment of the mean vertical profiles of velocity and temperature was studied, which was found qualitatively comparable with Belcher et al. (2003, J Fluid Mech. 488, 369-398)'s theoretical results. The urban heat island(UHI) was investigated imposing heat source in the region of large roughness elements. An elevated inversion layer, a phenomenon often observed in the urban area (Sang et al., J Wind Eng. Ind. Aesodyn. 87, 243-258)'s was successfully simulated above the canopy. The cool island(CI) was also investigated imposing heat sink to simply model the evaporation of plant canopy. An inversion layer was found very stable and robust within the canopy.

Difference schemes on non-uniform mesh and their application

High order accurate schemes are needed to simulate the multi-scale complex flow fields to get fine structures in simulation of the complex flows with large gradient of fluid parameters near the wall, and schemes on non-uniform mesh are desirable for many CFD (computational fluid dynamics) workers. The construction methods of difference approximations and several difference approximations on non-uniform mesh are presented. The accuracy of the methods and the influence of stretch ratio of the neighbor mesh increment on accuracy are discussed. Some comments on these methods are given, and comparison of the accuracy of the results obtained by schemes based on both non-uniform mesh and coordinate transformation is made, and some numerical examples with non-uniform mesh are presented.

Variable power singular interface elements for a crack normal to the bimaterial interface

Zero thickness crack tip interface elements for a crack normal to the interface between two materials are presented. The elements are shown to have the desired r(lambda-1) (0 < lambda < 1) singularity in the stress field at the crack tip and are compatible with other singular elements. The stiffness matrices of the quadratic and cubic interface element are derived. Numerical examples are given to demonstrate the applicability of the proposed interface elements for a crack perpendicular to the bimaterial interface.

Inhomogeneous isoparametric elements for stress analysis of composites

The concept of inhomogeneous element is proposed and the formulations of the inhomogeneous isoparametric elements for stress analysis of four kinds of problems are derived. As an example of applications of the inhomogeneous elements, the stress distribution in a cone-like composite syntheticrope termination is calculated.

Dynamic-response of a laminated plate with friction damping

快速Delaunay逐点插入网格生成算法

对插入形心的Delaunay逐点插入算法,提出按单元可插度分组的双向链表组数据结构,避免了对最大可插度单元的搜索。采用了邻接单元搜索、双向链表存储、随机方向搜索、邻接旋转、几何量继承等技术,使算法的计算时间与生成单元数近似呈线性关系,时间复杂度达到O(N1.05),N为生成单元数。算例表明,在一台AMD Athlon3200+(主频2.0 GHz)PC上,该算法的四面体单元生成速度达50 000个/s以上。

Computation of flow over a rotating body on unstructured Chimera mesh

Flow around moving boundary is ubiquitous in engineering applications. To increse the efficienly of the algorithm to handle moving boundaries is still a major challenge in Computational Fluid Dynamics (CFD). The Chimera grid method is one type of method to handle moving boundaries. A concept of domain de-composition has been proposed in this paper. In this method, sub-domains are meshed independently and governing equations are also solved separately on them. The Chimera grid method was originally used only on structured (curvilinear) meshes. However, in a problem which involves both moving boundary and complex geometry, the number of sub-domains required in a traditional (structured) Chimera method becomes fairly large. Thus the time required in the interior boundary locating, link-building and data exchanging also increases. The use of unstructured Chimera grid can reduce the time consumption significantly by the reduction of domain(block) number. Generally speaking, unstructured Chimera grid method has not been developed. In this paper, a well-known pressure correction scheme - SIMPLEC is modified and implemented on unstructured Chimera mesh. A new interpolation scheme regarding the pressure correction is proposed to prevent the possible decoupling of pressure. A moving-mesh finite volume approach is implemented in an inertial reference frame. This approach is then used to compute incompressible flow around a rotating circular and elliptic cylinder. These numerical examples demonstrate the capability of the proposed scheme in handling moving boundaries. The numerical results are in good agreement with other experimental and computational data in literature. The method proposed in this paper can be efficiently applied to more challenge cases such as free-falling objects or heavy particles in fluid.

An Improved CE/SE Scheme for Numerical Simulation of Gaseous and Two-Phase Detonations

A new structure of solution elements and conservation elements based on rectangular mesh was pro- posed and an improved space-time conservation element and solution element (CE/SE) scheme with sec- ond-order accuracy was constructed. Furthermore, the application of improved CE/SE scheme was extended to detonation simulation. Three models were used for chemical reaction in gaseous detonation. And a two-fluid model was used for two-phase (gas–droplet) detonation. Shock reflections were simu- lated by the improved CE/SE scheme and the numerical results were compared with those obtained by other different numerical schemes. Gaseous and gas–droplet planar detonations were simulated and the numerical results were carefully compared with the experimental data and theoretical results based on C–J theory. Mach reflection of a cellular detonation was also simulated, and the numerical cellular pat- terns were compared with experimental ones. Comparisons show that the improved CE/SE scheme is clear in physical concept, easy to be implemented and high accurate for above-mentioned problems.

Numerical simulation of condensation of bubbles under microgravity conditions by moving mesh method in the double-staggered grid

A numerical 2D method for simulation of two-phase flows including phase change under microgravity conditions is presented in this paper, with a level set method being coupled with the moving mesh method in the double-staggered grid systems. When the grid lines bend very much in a curvilinear grid, great errors may be generated by using the collocated grid or the staggered grid. So the double-staggered grid was adopted in this paper. The level set method is used to track the liquid-vapor interface. The numerical analysis is fulfilled by solving the Navier-Stokes equations using the SIMPLER method, and the surface tension force is modeled by a continuum surface force approximation. A comparison of the numerical results obtained with different numerical strategies shows that the double-staggered grid moving-mesh method presented in this paper is more accurate than that used previously in the collocated grid system. Based on the method presented in this paper, the condensation of a single bubble in the cold water under different level of gravity is simulated. The results show that the condensation process under the normal gravity condition is different from the condensation process under microgravity conditions. The whole condensation time is much longer under the normal gravity than under the microgravity conditions.