共查询到20条相似文献,搜索用时 109 毫秒
1.
非结构动网格在多介质流体数值模拟中的应用 总被引:1,自引:1,他引:0
采用非结构动网格方法对含多介质的流场进行数值模拟.采用改进的弹簧方法来处理由于边界运动而产生的网格变形.采用基于格心的有限体积方法求解守恒型的ALE(Arbitrary Lagrangiall-Eulerian)方程,控制面通量的计算采用HLLC(Hartem,Lax,van Leer,Contact)方法,采用几何构造的方法使空间达到二阶精度,时间离散采用四阶Runge-Kutta方法.物质界面的处理采用虚拟流体方法.本文对含动边界的激波管、水下爆炸等流场进行数值模拟,取得较好的结果,不同时刻界面的位置和整个扩张过程被准确模拟. 相似文献
2.
3.
欧拉坐标系下具有锐利相界面的可压缩多介质流动数值方法研究 总被引:1,自引:0,他引:1
可压缩多介质流动问题的数值模拟在国防和工业领域内均具有重要的研究价值,诸如武器设计、爆炸安全防护等,通常具有大变形、高度非线性等特点,是一项极具挑战性的研究课题. 本文提出了一种基于 Euler 坐标系的非结构网格、具有锐利相界面的二维和三维守恒型多介质流动数值方法,可用于模拟可压缩流体和弹塑性固体在极端物理条件下的大变形动力学行为. 利用分片线性的水平集函数重构出单纯形网格内分段线性的相界面,并在混合网格内构建出具有多种介质的相界面几何结构,理论上可以处理全局任意种介质、局部 3 种介质的多介质流动问题. 利用传统的有限体积格式来计算单元边界上同种介质间的数值通量,并通过在相界面法向上求解局部一维多介质 Riemann 问题的精确解来计算不同介质间的数值通量,保证了相界面上的通量守恒. 提出了一种非结构网格上的单元聚合算法,消除了由于网格被相界面分割成较小碎片、违反 CFL 条件,进而可能带来数值不稳定的问题. 针对一维多介质 Riemann 问题、激波与气泡相互作用问题、浅埋爆炸问题、空中强爆炸冲击波和典型坑道内冲击波传播问题开展了数值模拟研究,将计算结果与相关的理论、实验结果进行比对,验证了数值方法的正确性和可靠性. 相似文献
4.
可压缩多介质流动问题的数值模拟在国防和工业领域内均具有重要的研究价值,诸如武器设计、爆炸安全防护等,通常具有大变形、高度非线性等特点,是一项极具挑战性的研究课题. 本文提出了一种基于 Euler 坐标系的非结构网格、具有锐利相界面的二维和三维守恒型多介质流动数值方法,可用于模拟可压缩流体和弹塑性固体在极端物理条件下的大变形动力学行为. 利用分片线性的水平集函数重构出单纯形网格内分段线性的相界面,并在混合网格内构建出具有多种介质的相界面几何结构,理论上可以处理全局任意种介质、局部 3 种介质的多介质流动问题. 利用传统的有限体积格式来计算单元边界上同种介质间的数值通量,并通过在相界面法向上求解局部一维多介质 Riemann 问题的精确解来计算不同介质间的数值通量,保证了相界面上的通量守恒. 提出了一种非结构网格上的单元聚合算法,消除了由于网格被相界面分割成较小碎片、违反 CFL 条件,进而可能带来数值不稳定的问题. 针对一维多介质 Riemann 问题、激波与气泡相互作用问题、浅埋爆炸问题、空中强爆炸冲击波和典型坑道内冲击波传播问题开展了数值模拟研究,将计算结果与相关的理论、实验结果进行比对,验证了数值方法的正确性和可靠性. 相似文献
5.
6.
推导了轴对称几何下的MOF(Moment of Fluid)界面重构,将其与多介质ALE方法相耦合,形成MOFMMALE方法,并应用于多介质大变形流动问题的数值模拟研究。数值算例表明,耦合MOF界面重构的多介质ALE方法是求解多介质大变形流动问题的有效手段,并且具有很好的界面精度和分辨率。 相似文献
7.
8.
9.
10.
在均匀网格上求解对流占优问题时,往往会产生数值震荡现象,因此需要局部加密网格来提高解的精度。针对对流占优问题,设计了一种新的自适应网格细化算法。该方法采用流线迎风SUPG(Petrov-Galerkin)格式求解对流占优问题,定义了网格尺寸并通过后验误差估计子修正来指导自适应网格细化,以泡泡型局部网格生成算法BLMG为网格生成器,通过模拟泡泡在区域中的运动得到了高质量的点集。与其他自适应网格细化方法相比,该方法可在同一框架内实现网格的细化和粗化,同时在所有细化层得到了高质量的网格。数值算例结果表明,该方法在求解对流占优问题时具有更高的数值精度和更好的收敛性。 相似文献
11.
Ravi K. Burla Ashok V. Kumar Bhavani V. Sankar 《International Journal of Solids and Structures》2009,46(11-12):2514-2526
A method that uses a structured grid to perform micromechanical analysis for determining effective properties of a composite microstructure is presented. This approach eliminates the need for constructing a mesh that has nodes along the interfaces between constituent materials of the composite. Implicit boundary method is used to ensure that interface conditions are satisfied at the material boundaries. In this method, solution structures for test and trial functions are constructed using approximate step functions such that the interface conditions are satisfied, even if there are no nodes on the material interface boundary. Since a structured grid does not conform to the geometry of the analysis domain, the geometry of the microstructure is defined independently using equations of the interface boundary curves/surfaces. Structured grids that overlap the geometry are easy to generate, and the elements in the grid are regular shaped and undistorted. A numerical example is presented to demonstrate that the proposed solution structure accurately models the solution across material interface, and convergence analysis is performed to show that the method converges as the grid density is increased. Fiber reinforced microstructures are analyzed to compute the effective elastic properties using both 2D and 3D models to show that the results match closely with the ones available in the literature. 相似文献
12.
A stencil-like volume of fluid (VOF) method is proposed for tracking free interface. A stencil on a grid cell is worked out according to the normal direction of the interface, in which only three interface positions are possible in 2D cases, and the interface can be reconstructed by only requiring the known local volume fraction information. On the other hand, the fluid-occupying-length is defined on each side of the stencil, through which a unified fluid-occupying volume model and a unified algorithm can be obtained to solve the interface advection equation. The method is suitable for the arbitrary geometry of the grid cell, and is extendible to 3D cases. Typical numerical examples show that the current method can give "sharp" results for tracking free interface. 相似文献
13.
We develop an Eulerian fixed grid numerical method for calculating multi‐material fluid flows. This approach relates to the class of interface capturing methods. The fluid is treated as a heterogeneous mixture of constituent materials, and the material interface is implicitly captured by a region of mixed cells that have arisen owing to numerical diffusion. To suppress this numerical diffusion, we propose a composite Riemann problem (CRP), which describes the decay of an initial discontinuity in the presence of a contact point between two different fluids, which is located off the initial discontinuity point. The solution to the CRP serves to calculate multi‐material no mixed numerical flux without introducing any material diffusion. We discuss the CRP solution and its implementation in the multi‐material fluid Godunov method. Numerical results show that a simple framework of the CRP greatly improves capturing material interfaces in the Godunov method and reproduces many of the advantages of more complicated interface tracking multi‐material treatments. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
14.
Y.-G. Zhang J. Ballmann P. Behruzi 《Archive of Applied Mechanics (Ingenieur Archiv)》1998,68(7-8):489-500
Summary Explicit finite-difference procedures are proposed to simulate the elastodynamic response of layered media. To avoid large
numerical dispersion and dissipation, different grid sizes are used for different material layers so that the CFL number in
each layer can be set to one in the limit. The master-slave technique is applied to deal with the mismatch of grids along
the material interface. The interface can be either bonded or unbonded. The dynamic contact between two layers along an unbonded
interface is also considered. Numerical examples are presented.
Received 4 August 1997; accepted for publication 24 April 1998 相似文献
15.
This paper presents an approach to the generation of unstructured surface meshes for Computer‐Aided Design (CAD) surface models represented as lists of polygons with minimum user interventions. Stereolithography (STL) data are adopted as an interface between a CAD system and the surface grid generator. STL files often include problems such as overlapping surfaces, gaps, and intersections. They have to be revised quickly and automatically before the surface models are used for the background grid of the surface grid generation. In this paper, we describe an automatic revision method for use with STL‐defined surface models. The advancing front method using geometric features is adopted directly on the modified STL surfaces. The capability of the method is demonstrated for several 3D surface models. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
16.
A numerical method is developed for solving the 3D, unsteady, incompressible Navier–Stokes equations in curvilinear coordinates containing immersed boundaries (IBs) of arbitrary geometrical complexity moving and deforming under forces acting on the body. Since simulations of flow in complex geometries with deformable surfaces require special treatment, the present approach combines a hybrid immersed boundary method (HIBM) for handling complex moving boundaries and a material point method (MPM) for resolving structural stresses and movement. This combined HIBM & MPM approach is presented as an effective approach for solving fluid–structure interaction (FSI) problems. In the HIBM, a curvilinear grid is defined and the variable values at grid points adjacent to a boundary are forced or interpolated to satisfy the boundary conditions. The MPM is used for solving the equations of solid structure and communicates with the fluid through appropriate interface‐boundary conditions. The governing flow equations are discretized on a non‐staggered grid layout using second‐order accurate finite‐difference formulas. The discrete equations are integrated in time via a second‐order accurate dual time stepping, artificial compressibility scheme. Unstructured, triangular meshes are employed to discretize the complex surface of the IBs. The nodes of the surface mesh constitute a set of Lagrangian control points used for tracking the motion of the flexible body. The equations of the solid body are integrated in time via the MPM. At every instant in time, the influence of the body on the flow is accounted for by applying boundary conditions at stationary curvilinear grid nodes located in the exterior but in the immediate vicinity of the body by reconstructing the solution along the local normal to the body surface. The influence of the fluid on the body is defined through pressure and shear stresses acting on the surface of the body. The HIBM & MPM approach is validated for FSI problems by solving for a falling rigid and flexible sphere in a fluid‐filled channel. The behavior of a capsule in a shear flow was also examined. Agreement with the published results is excellent. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
17.
采用连续切片的方法获取准脆性材料的表面图像,利用数字图像处理技术检测材料的细观结构并进行矢量化。通过一种简单的变换,将每一切片矢量化的细观结构转换成单层的三维结构,然后将这些切片连续的细观结构逐层叠加,形成整个试件的三维真实细观结构,并针对准脆性材料图像的特点,编制了能够批量处理数字图像并进行细观结构矢量化的程序,建立了与有限元三维网格模型之间的数据接口,模型数据可直接导入岩石三维破裂过程分析RFPA3D系统中,研究真实细观结构对准脆性材料破坏力学行为的影响。以颗粒材料为例,分析了在单轴受压情况下的三维空间裂纹的产生及扩展过程,计算结果显示颗粒分布与界面显著影响材料的破裂模式。 相似文献
18.
Hideo Koguchi Joviano Antonio da Costa 《International Journal of Solids and Structures》2010,47(22-23):3131-3140
The stress singularity that occurs at a vertex in a joint with a slanted side surface is investigated. The orders of stress singularity at a vertex and at a point on stress singularity lines for various material properties are determined using eigenanalysis. The stress distribution on an interface and the intensity of stress singularity at the vertex are investigated using BEM. It is shown that the order of stress singularity at the vertex in the joints can be reduced by slanting a side surface so as to decrease the angle between the interface and the side surface. The results of BEM analysis reveal that the distribution of stress on the interface is influenced by the slanted side surface. Finally, the 3D intensities of the singularity for stress components which are continuous at the interface are newly defined and determined for various material combinations. 相似文献
19.
We present a fixed‐grid finite element technique for fluid–structure interaction problems involving incompressible viscous flows and thin structures. The flow equations are discretised with isoparametric b‐spline basis functions defined on a logically Cartesian grid. In addition, the previously proposed subdivision‐stabilisation technique is used to ensure inf–sup stability. The beam equations are discretised with b‐splines and the shell equations with subdivision basis functions, both leading to a rotation‐free formulation. The interface conditions between the fluid and the structure are enforced with the Nitsche technique. The resulting coupled system of equations is solved with a Dirichlet–Robin partitioning scheme, and the fluid equations are solved with a pressure–correction method. Auxiliary techniques employed for improving numerical robustness include the level‐set based implicit representation of the structure interface on the fluid grid, a cut‐cell integration algorithm based on marching tetrahedra and the conservative data transfer between the fluid and structure discretisations. A number of verification and validation examples, primarily motivated by animal locomotion in air or water, demonstrate the robustness and efficiency of our approach. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
20.
The objective of this work is to develop a sliding interface method for simulations involving relative grid motion that is fast and efficient and involves no grid deformation, remeshing, or hole cutting. The method is implemented into a parallel, node‐centred finite volume, unstructured viscous flow solver. The rotational motion is accomplished by rigidly rotating the subdomain representing the moving component. At the subdomain interface boundary, the faces along the interface are extruded into the adjacent subdomain to create new volume elements forming a one‐cell overlap. These new volume elements are used to compute a flux across the subdomain interface. An interface flux is computed independently for each subdomain. The values of the solution variables and other quantities for the nodes created by the extrusion process are determined by linear interpolation. The extrusion is done so that the interpolation will maintain information as localized as possible. The grid on the interface surface is arbitrary. The boundary between the two subdomains is completely independent from one another; meaning that they do not have to connect in a one‐to‐one manner and no symmetry or pattern restrictions are placed on the grid. A variety of numerical simulations were performed on model problems and large‐scale applications to examine conservation of the interface flux. Overall solution errors were found to be comparable to that for fully connected and fully conservative simulations. Excellent agreement is obtained with theoretical results and results from other solution methodologies. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献