首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 218 毫秒
1.
许晓阳  赵雨婷  李家宇  余鹏 《力学学报》2023,55(5):1099-1112
非等温黏弹性流体广泛存在于自然界和工业生产中,准确预测黏弹性流体的非等温流动机理和复杂流变特性有着重要的应用价值.文章提出一种改进的光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法对非等温黏弹性复杂流动进行了数值模拟,其中流体的黏弹特性通过eXtended Pom-Pom本构模型来表征.为了提高模拟结果的精度,采用了一种核函数梯度的修正算法;为了灵活地施加边界条件,发展了边界粒子和虚拟粒子相联合的边界处理方法;为了消除流动过程中的拉伸不稳定性,施加了粒子迁移技术.运用改进SPH方法数值模拟了液滴撞击固壁和F型腔注塑成型问题,通过与Basilisk软件得到的结果进行比较验证了改进SPH方法求解非等温黏弹性流体的有效性.通过利用不同粒子初始间距进行计算,评价了改进SPH方法的数值收敛性.研究了非等温流动相较于等温流动的不同流动特征,深入分析了不同热流变参数对流动过程的影响.数值结果表明,文章提出的改进SPH方法可稳定、准确地描述非等温黏弹性复杂流动的传热机理、复杂流变特性和自由面变化特性.  相似文献   

2.
采用光滑粒子动力学SPH方法建立液滴冲击弹性基底的流固耦合数值模型,给出描述粘性流体和弹性固体运动的SPH离散方程和数值处理格式,引入人工耗散项来抑制标准SPH方法的数值震荡。为模拟液滴的表面张力效应,通过精确检测边界粒子,采用拉格朗日插值方法计算表面法向量和曲率,结合界面理论中的连续表面力CSF方法,建立了适用于自由表面液滴的表面力模型,方形液滴变形的模拟结果与拉普拉斯理论解吻合较好。随后,采用SPH流固耦合模型模拟1.0 mm直径水滴以不同速度(0.2 m/s~3.0 m/s)冲击两种薄板型基底,分析了基底弹性变形对液滴铺展、收缩以及回弹行为的影响。  相似文献   

3.
NUMERICAL SIMULATION OF WIND-BLOWN SAND MOVEMENT BASED ON SPH   总被引:1,自引:0,他引:1  
运用光滑粒子流体动力学(smoothed particles hydrodynamics, SPH)方法对沙粒和气流的相互耦合运动特性进行了分析,研究提出了风沙流的SPH数值方法并进行了数值模拟. 首先提出了风沙流的SPH建模方法和基本理论,建立了风沙流动的SPH数值模拟平台. 其次通过建立风沙流的SPH 模型并施加边界条件,对自然风作用下沙粒的运动情况进行了数值模拟,详细分析了沙粒运动轨迹及特性,最后通过与相关研究成果对比分析,验证了完善后的SPH方法有效性. 通过考虑气流场的可变性,在风沙流SPH计算模型中引入了加载(起风)和卸载(停风)方式,观察并对比分析了沙粒的运动轨迹和特性. 为进一步研究风沙流的实时动态非线性行为提供了SPH理论基础和数值分析方法.  相似文献   

4.
采用光滑粒子动力学SPH(Smoothed Particle Hydrodynamics)方法对三维溃坝流问题进行了数值模拟。为了逼真地模拟出坝内水体与壁面间相互作用而产生的水花飞溅、融合以及近壁面流动等现象,加入了混合长度形式的湍流模型。为了有效地防止粒子穿透固壁,提出了一种新型的适合三维数值模拟的固壁边界处理方法。应用SPH方法对三维溃坝流进行了数值模拟,并分别考虑了未添加障碍物和添加圆柱障碍物两种情形。计算结果表明,改进SPH方法能够精细地捕捉溃坝流在不同时刻的自由液面,并获得稳定而精确的数值结果。  相似文献   

5.
光滑粒子流体动力学(SPH)法是一种无网格的拉格朗日效值方法,广泛应用于计算流体领域模拟复杂自由表面流问题.SPH方法的主要缺点就是计算量过大,而基于GPU的并行计算方法可使SPH计算得到有效加速.本文应用基于GPU的SPH并行计算方法研究了二维楔形体的入水砰击问题.数值计算结果与文献中对应的解析解比较一致,验证了基于GPU的SPH方法的精度和可靠性.仿真结果同时显示基于GPU的并行计算方法可使SPH计算速度得到显著提高.  相似文献   

6.
SPH方法Delaunay三角刨分与自由液面重构   总被引:1,自引:1,他引:0  
光滑粒子法(SPH)作为一种拉格朗日型无网格方法,兼具欧拉网格方法和拉格朗日网格方法的优势,已经成功应用于科学和工程的众多领域。SPH方法后处理一般基于无规则分布的粒子,不如网格类方法后处理简便、直接。另外,SPH方法模拟自由液面流动等问题时,通过粒子位置难以重构自由液面的准确位置。发展一种基于Delaunay三角刨分的SPH后处理方法,即先基于SPH粒子位置利用Delaunay三角刨分建立三角网格,然后将粒子信息转化成网格单元/节点信息,从而可以在三角网格上进行后处理,实现基于网格方法的后处理功能,并可以在三角网格上直接提取或重构自由液面。将本文的方法应用到液滴碰撞和溃坝流SPH模拟结果的后处理中,得到了非常好的结果,表明本文的方法有效可靠。  相似文献   

7.
陈飞国  葛蔚 《力学学报》2021,53(9):2357-2373
光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)具有粒子方法的无网格和全拉格朗日特征, 适用于具有界面大变形、不连续性和多物理场的多相流的高精度模拟. SPH方法模拟多相流已有大量报道, 具体的实现方式也大不相同. 本文首先阐述了采用SPH方法模拟流体的基本控制方程, 以及求解过程中需要考虑的流体压力求解、表面张力、固体边界等问题. 整理和总结了基于SPH方法进行多相流模拟的主要实现方式: (1)双流体模型的拉格朗日求解器: 两相离散为两组独立SPH粒子, 并用显式相间作用耦合两相; (2)多相SPH方法: SPH方法对多相流模拟的自然延伸, 相间作用由SPH参数隐式描述; (3) SPH与其他离散方法的耦合: 差异较大的两相各自采用不同离散方法, 发挥不同拉格朗日方法的优点; (4) SPH和基于网格方法的耦合: 网格方法处理简单的单相流动主体, 获得精度和效率间的平衡. 另外, 还在模拟参数物理化等方面论述了与SPH方法模拟多相流相关的一些改进和修正方法, 并在最后讨论和建议了提高多相流SPH模拟效率和精度的措施.   相似文献   

8.
修正SPH方法在自由表面模拟中的应用   总被引:1,自引:0,他引:1  
对传统光滑粒子动力学(SPH)方法进行修正,提出了一种密度初始化方法,同时采用一种新的固壁边界处理方法,模拟溃坝问题。通过液滴旋转和无透空块体溃坝问题的模拟验证了修正SPH方法的有效性和在自由表面模拟中的准确性,分析了密度初始化对流动的影响;数值结果表明,修正SPH方法提高了数值计算的精度和稳定性。最后应用修正SPH方法模拟了有透空块体和挡板紧挨水柱的溃坝现象,比较了有无透空块体两种情况下右端直墙上压力变化情况,结果表明,透空块体可使右端直墙上的压力减小,有无挡板、挡板位置和水柱长高比对溃坝现象有重要影响。  相似文献   

9.
应用有限元(FE)-光滑粒子流体动力学(SPH)耦合法模拟了具有自由表面的不可压流体与结构的相互作用问题.流体和结构分别采用SPH法和有限元法同时求解,两者在交界面处的相互作用通过接触算法进行处理.为了避免隐式计算压力,通过引入人工压缩率,将不可压流体近似为人工可压缩流体.采用FE-SPH耦合法对弹性板在随时间变化的水压作用下的变形以及倒塌水柱冲击弹性结构两个问题进行了模拟.模拟结果与实验结果以及其他已有数值结果符合良好,说明FE-SPH耦合法用于流体与结构相互作用问题的模拟是可行和有效的.  相似文献   

10.
杨秋足  徐绯  王璐  杨扬 《力学学报》2019,51(3):730-742
多相流界面存在密度、黏性等物理场间断,直接采用传统光滑粒子水动力学(smoothedparticle hydrodynamics,SPH)方法进行数值模拟,界面附近的压力和速度存在震荡.一套基于黎曼解能够处理大密度比的多相流SPH计算模型被提出,该模型利用黎曼解在处理接触间断问题方面的优势,将黎曼解引入到SPH多相流计算模型中,为了能够准确求解多相流体物理黏性、减小黎曼耗散,对黎曼形式的SPH动量方程进行了改进,又将Adami固壁边界与黎曼单侧问题相结合来施加多相流SPH固壁边界,同时模型中考虑了表面张力对小尺度异相界面的影响,该模型没有添加任何人工黏性、人工耗散和非物理人工处理技术,能够反应多相流真实物理黏性和物理演变状态.采用该模型首先对三种不同粒子间距离散下方形液滴震荡问题进行了数值模拟,验证了该模型在处理异相界面的正确性和模型本身的收敛性;后又通过对Rayleigh--Taylor不稳定、单气泡上浮、双气泡上浮问题进行了模拟计算,结果与文献对比吻合度高,异相界面捕捉清晰,结果表明,本文改进的多相流SPH模型能够稳定、有效的模拟大密度比和黏性比的多相流问题.   相似文献   

11.
In this work, a corrected symmetric and periodic density reinitialized SPH (CSPDR‐SPH) method is proposed and extended to simulate the viscoelastic free surface flows based on the Phan–Thien–Tanner model. The improvements mainly lie in deriving a corrected symmetric kernel gradient, and combining it with a periodic density reinitialization procedure. In addition, a simple artificial viscosity and a simple artificial stress form are adopted. Thus, the CSPDR‐SPH method has higher accuracy and better stability than the SPH method, and conserves both linear and angular momentums. The consistency and convergence of the CSPDR‐SPH method are justified by approximating a function in one and two dimensions. The merits of CSPDR‐SPH method are demonstrated by several benchmarks. The simple flow in a two‐dimensional channel is investigated to show the capability of the CSPDR‐SPH method to simulate the viscoelastic free surface flow. Then the CSPDR‐SPH method is extended to simulate the impacting drop problem. Numerical results show that the CSPDR‐SPH method can precisely capture the viscoelastic free surface. The Reynolds number, Weissenberg number and elongation parameter have remarkable effect on the flows. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

12.
This work is concerned with the development of a numerical method capable of simulating two-dimensional viscoelastic free surface flows governed by the non-linear constitutive equation PTT (Phan-Thien–Tanner). In particular, we are interested in flows possessing moving free surfaces. The fluid is modelled by a marker-and-cell type method and employs an accurate representation of the fluid surface. Boundary conditions are described in detail and the full free surface stress conditions are considered. The PTT equation is solved by a high order method which requires the calculation of the extra-stress tensor on the mesh contour. The equations describing the numerical technique are solved by the finite difference method on a staggered grid. In order to validate the numerical method fully developed flow in a two-dimensional channel was simulated and the numerical solutions were compared with known analytic solutions. Convergence results were obtained throughout by using mesh refinement. To demonstrate that complex free surface flows using the PTT model can be computed, extrudate swell and a jet flowing onto a rigid plate were simulated.  相似文献   

13.
14.
15.
In this paper, an incompressible smoothed particle hydrodynamics (SPH) method is presented to solve unsteady free-surface flows. Both Newtonian and viscoelastic fluids are considered. In the case of viscoelastic fluids, both the Maxwell and Oldroyd-B models are investigated. The proposed SPH method uses a Poisson pressure equation to satisfy the incompressibility constraints. The solution algorithm is an explicit predictor-corrector scheme and employs an adaptive smoothing length based on density variations. To alleviate the numerical difficulties encountered when fluid is highly stretched, an artificial stress term is incorporated into the momentum equation which reduces the risk of unrealistic fractures in the material. Two challenging test cases, the impacting drop and the jet buckling problems, are solved to demonstrate the capability of the proposed scheme in handling viscoelastic flows with complex free surfaces. The jet buckling test case was solved for a wide range of Weissenberg numbers. It was shown that in all cases the method is stable and fairly accurate and agrees well with the available data.  相似文献   

16.
The kernel gradient free (KGF) smoothed particle hydrodynamics (SPH) method is a modified finite particle method (FPM) which has higher order accuracy than the conventional SPH method. In KGF‐SPH, no kernel gradient is required in the whole computation, and this leads to good flexibility in the selection of smoothing functions and it is also associated with a symmetric corrective matrix. When modeling viscous incompressible flows with SPH, FPM or KGF‐SPH, it is usual to approximate the Laplacian term with nested approximation on velocity, and this may introduce numerical errors from the nested approximation, and also cause difficulties in dealing with boundary conditions. In this paper, an improved KGF‐SPH method is presented for modeling viscous, incompressible fluid flows with a novel discrete scheme of Laplacian operator. The improved KGF‐SPH method avoids nested approximation of first order derivatives, and keeps the good feature of ‘kernel gradient free’. The two‐dimensional incompressible fluid flow of shear cavity, both in Euler frame and Lagrangian frame, are simulated by SPH, FPM, the original KGF‐SPH and improved KGF‐SPH. The numerical results show that the improved KGF‐SPH with the novel discrete scheme of Laplacian operator are more accurate than SPH, and more stable than FPM and the original KGF‐SPH. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

17.
In this paper, a corrected particle method based on the smoothed particle hydrodynamics (SPH) method with high-order Taylor expansion (CSPH-HT) for solving the vis-coelastic flow is proposed and investigated. The validity and merits of the CSPH-HT method are first tested by solv-ing the nonlinear high order Kuramoto-Sivishinsky equation and simulating the drop stretching, respectively. Then the flow behaviors behind two stationary tangential cylinders of polymer melt, which have been received little attention, are investigated by the CSPH-HT method. Finally, the CSPH-HT method is extended to the simulation of the filling process of the viscoelastic fluid. The numerical results show that the CSPH-HT method possesses higher accuracy and stability than other corrected SPH methods and is more reliable than other corrected SPH methods.  相似文献   

18.
The gridless smoothed particle hydrodynamics (SPH) method is now commonly used in computational fluid dynamics (CFD) and appears to be promising in predicting complex free‐surface flows. However, increasing flow complexity requires appropriate approaches for taking account of turbulent effects, whereas some authors are still working without any turbulence closure in SPH. A review of recently developed turbulence models adapted to the SPH method is presented herein, from the simplistic point of view of a one‐equation model involving mixing length to more sophisticated (and thus realistic) models like explicit algebraic Reynolds stress models (EARSM) or large eddy simulation (LES). Each proposed model is tested and validated on the basis of schematic cases for which laboratory data, theoretical or numerical solutions are available in the general field of turbulent free‐surface incompressible flows (e.g. open‐channel flow and schematic dam break). They give satisfactory results, even though some progress should be made in the future in terms of free‐surface influence and wall conditions. Recommendations are given to SPH users to apply this method to the modelling of complex free‐surface turbulent flows. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

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

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