非牛顿流体在Hele-Shaw模型中的流动

本文研究了非牛顿流体在Hele-Shaw模型中的流动特性,并用摄动法求出了Bingham流体、Power-law流体的速度分布、压力降与流量的关系,分析了屈服应力及幂律指数对物理参数的影响.     

非等温广义牛顿流体的无网格模拟

基于原始变量法,将RPIM推广到非等温广义牛顿流动问题的求解.为了减少未知量的个数,压力和速度采用罚函数方法耦合;同时采用积分降阶技术以保证获得满意的数值解.数值实验结果表明,采用罚函数法处理速度和压力在无网格方法中同样适用,且RPIM用于非等温广义牛顿流动问题的求解时具有易施加本质边界条件、计算精度高和收敛性较快的优点.     

膨胀性非牛顿流体多孔介质流动的模拟分析

在表征体元尺度采用格子Boltzmann方法分析膨胀性非牛顿流体在多孔介质中的流动,基于二阶矩模型在演化方程中引入表征介质阻力的作用力项,求解描述渗流模型的广义Navier-Stokes方程．采用局部法计算形变速率张量,通过循环迭代得到非牛顿粘度和松弛时间．对多孔介质的Poiseuille流动进行分析,通过比较发现结果与孔隙尺度的解析解十分吻合,并且收敛较快,表明方法合理有效．分析了渗透率和幂律指数对速度和压力降的影响,研究结果表明,膨胀性流体的多孔介质流动不符合达西规律,压力降的增加幅度小于渗透率的减小幅度．当无量纲渗透率Da小于10-5时,流道中的速度呈现均匀分布,并且速度分布随着幂律指数的减小趋于平滑．压力降随着幂律指数的增加而增加,Da越大幂律指数对压力降的影响越明显．     

有限应变下蠕变—弹塑性—损伤耦合本构模型的数值模拟

本文提出了一个模拟有限元应变下各向同性蠕变－弹塑性－损伤耦合本构行为的数值模型。模型基于不可逆过程的热动力学和内状态变量理论、以及应力就变本构关系的超弹性形式描述和变形梯度的弹性及非弹性部分的乘式分解。     

电流动聚焦中非牛顿流体射流影响因素的实验研究

在一套流动聚焦装置上加载高压直流电场形成电流动聚焦,并开展了非牛顿流体带电射流的不稳定性特性实验研究。实验在自行设计的装置系统上完成,获得了电流动聚焦中非牛顿流体射流的流动状态,考察了不同控制参数下射流形态的变化。结果表明,由于非牛顿流体具有粘弹性,与牛顿流体相比,非牛顿流体带电射流体现了更复杂的流动特点。这些实验结果为我们理解复杂条件下非牛顿流体射流的流动机理提供了参考,也有助于深入的实验分析和理论研究奠定了基础。     

土体的本构模型和超重力物理模拟

数值模拟和物理模拟是分析土体沉降和稳定性的主要手段. 本构模型作为描述土体应力应变关系的数学表达式, 是数值模拟的基础. 土体具有碎散性, 这一基本物理特性导致了其具有压硬性、摩擦性和剪胀性, 这是土的力学特性区别于金属的主要特征, 在土体的本构模型中必须反映这3个基本特性. 传统土力学将土体的变形和强度分离考虑, 分别采用弹性理论和基于刚塑性模型的极限平衡理论分析, 虽然应用广泛, 但由于不能全面地反映土的基本力学特性, 计算结果的精度常常难以满足定量分析的需要. 剑桥模型作为第一个全面反映压硬性、摩擦性和剪胀性的弹塑性本构模型, 实现了变形和强度的统一, 能较好地描述饱和正常固结黏土的应力应变关系, 被视为是现代土力学的开端; 统一硬化模型通过引入一个独特的硬化参数进一步发展了剑桥模型, 将适用范围扩大到超固结黏土. 作者认为, 未来岩土体本构模型研究的挑战是: 如何考虑岩土体在受力过程中土骨架相变与多场耦合, 以解决目前本构模型尚无法定量分析的能源、交通、环境和水利相关的重大岩土工程问题. 超重力物理模拟具有缩尺效应和缩时效应, 克服了常重力物理模拟中模型的应力水平低于原型的缺点, 特别适用于大尺度、长历时问题的模拟. 相较数值模拟, 超重力物理模拟的优势在于能够检验本构模型的合理性, 揭示本构模型无法描述的未知特性. 最后, 介绍了采用数值模拟和物理模拟联合分析大直径钢管桩水平受荷特性的工程案例.      

使用新型物态方程的超高速碰撞物质点法模拟

为更准确地对超高速碰撞进行数值模拟、获得与实验结果相似度更高的碎片云形态,利用分子动力学方法求解材料的冷能、冷压,并结合Grover定标律方程,建立了一种表达形式简洁、可处理相变影响的新型物态方程,并代入自编柱坐标物质点法计算程序,使用新型物态方程计算所得的碎片云与使用Mie-Grüneisen、Tillotson等传统物态方程的计算结果相比,在尺寸、形态方面均能够与实验结果更好地吻合,证明了新型物态方程的有效性。

     

广义双剪粘弹塑性本构模型在渤海海冰动力学中的应用

在连续介质力学基础上建立了一个广义双剪粘弹塑性海冰动力学本构模型。该模型在海冰屈服前采用Kelvin-Vogit粘弹性模型,考虑中间主应力和静水压力对海冰屈服的影响选用广义双剪应力屈服准则作为海冰屈服判据,屈服后采用相关联的正则流动法则。采用该本构模型对渤海海冰动力过程进行了48小时数值模拟,讨论了辽东湾海冰的厚度、密集度、冰速和主应力的分布规律,其中海冰厚度分布与卫星遥感资料符合良好,从而有效地验证了该广义双剪粘弹塑性本构模型在海冰动力学中的可靠性。     

颗粒材料柱体崩塌物质点法数值模拟

为了研究颗粒材料崩塌的运动规律和堆积特性,采用物质点法对颗粒材料柱体崩塌试验进行数值模拟,并将模拟结果与试验结果进行对比验证。对颗粒材料柱体崩塌过程中颗粒的流动特性(滑动距离、堆积高度、速度、能量和动能通量的演化)进行了分析。进一步探究了颗粒材料柱体高宽比对颗粒流动能通量的影响,从而反应颗粒材料柱体崩塌过程中颗粒流的破坏能力。颗粒材料柱体高宽比越大,颗粒材料柱体外侧边缘颗粒速度越大,其溃散的程度更加强烈,并且滑动距离和动能均在增大。对于动能通量的分布,水平方向越靠近初始颗粒材料柱体,动能通量越大。     

基于内时大变形弹塑性本构模型的回弹和二次屈服分析

采用了一个由简单机构建立的材料热力学一致性弹塑性大变形本构模型，该模型满足热力学基本的不可逆约束条件，并考虑到了储存在残余微应力场中的能量，用相应的算法，编制了一个用户子程序UMAT，并将其嵌入到商用软件ABAQUS中，对圆柱形工件在弹性工具中镦粗和闭合模镦粗时的回弹和二次屈服进行了分析。比较模拟和实验所测得的载荷一位移关系，可以看出在镦粗的第一阶段发生塑性变形，工件在轴向加载过程中产生的横向变形受到了由于模具的弹性变形引起的径向约束，载荷迅速增加；第二阶段轴向压力减少，模具回弹使得工件的径向变形得以恢复，但随着塑性应变的增加，发生工作强化；第三阶段模具约束解除，工作的径向尺寸增加。用该本构模型对闭合模镦粗过程的计算，表明计算结果与实验结果是一致的。     

A study on non-steady flow of upper-convected maxwell fluid in tube

The transient response of an upper-convected Maxwell fluid flow in a circular tube is analysed by variational approach of Kantorovich and the method of finite difference. The solution of the variational method is in agreement with the numerical results by the difference schemes. The results show that the method of Kantorovich is suitable for the study of non-steady flow of non-Newtonian fluids and the effect of elasticity of the fluid has an influence on the non-steady flow. project supported by National Natural Science Foundation of China     

Coupling between finite element method and material point method for problems with extreme deformation

As a Lagrangian meshless method, the material point method (MPM) is suitable for dynamic problems with extreme deformation, but its efficiency and accuracy are not as good as that of the finite element method (FEM) for small deformation problems. Therefore, an algorithm for the coupling of FEM and MPM is proposed to take advantages of both methods. Furthermore, a conversion scheme of elements to particles is developed. Hence, the material domain is firstly discretized by finite elements, and then the distorted elements are automatically converted into MPM particles to avoid element entanglement. The interaction between finite elements and MPM particles is implemented based on the background grid in MPM framework. Numerical results are in good agreement with experimental data and the efficiency of this method is higher than that of both FEM and MPM.     

极端变形问题的物质点法研究进展

物质点法采用了拉格朗日和欧拉双重描述,将物体离散为一组在固定于空间的网格(欧拉描述)中运动的质点(拉格朗日描述),有效地综合了拉格朗日法和欧拉法的优点,是分析冲击爆炸等极端变形问题的一种有效方法. 本文系统总结了本课题组近年来针对冲击、爆炸和流固耦合等极端变形问题在物质点法算法方面的研究进展,并简要介绍了开发的三维显式并行物质点法数值仿真软件MPM3D及其在超高速碰撞、侵彻、爆炸、边坡失效、金属切削和流固耦合等问题中的应用.     

Material point method of modelling and simulation of reacting flow of oxygen

Aerospace vehicles are continually being designed to sustain flight at higher speeds and higher altitudes than previously attainable. At hypersonic speeds, gases within a flow begin to chemically react and the fluid's physical properties are modified. It is desirable to model these effects within the Material Point Method (MPM). The MPM is a combined Eulerian–Lagrangian particle-based solver that calculates the physical properties of individual particles and uses a background grid for information storage and exchange. This study introduces chemically reacting flow modelling within the MPM numerical algorithm and illustrates a simple application using the AeroElastic Material Point Method (AEMPM) code. The governing equations of reacting flows are introduced and their direct application within an MPM code is discussed. A flow of 100% oxygen is illustrated and the results are compared with independently developed computational non-equilibrium algorithms. Observed trends agree well with results from an independently developed source.     

A three-dimensional fictitious domain method for incompressible fluid flow problems

A new Galerkin finite element method for the solution of the Navier–Stokes equations in enclosures containing internal parts which may be moving is presented. Dubbed the virtual finite element method, it is based upon optimization techniques and belongs to the class of fictitious domain methods. Only one volumetric mesh representing the enclosure without its internal parts needs to be generated. These are rather discretized using control points on which kinematic constraints are enforced and introduced into the mathematical formulation by means of Lagrange multipliers. Consequently, the meshing of the computational domain is much easier than with classical finite element approaches. First, the methodology will be presented in detail. It will then be validated in the case of the two-dimensional Couette cylinder problem for which an analytical solution is available. Finally, the three-dimensional fluid flow inside a mechanically agitated vessel will be investigated. The accuracy of the numerical results will be assessed through a comparison with experimental data and results obtained with a standard finite element method. © 1997 John Wiley & Sons, Ltd.     

A coupled immersed boundary‐lattice Boltzmann method with smoothed point interpolation method for fluid‐structure interaction problems

The immersed boundary‐lattice Boltzmann method has been verified to be an effective tool for fluid‐structure interaction simulation associated with thin and flexible bodies. The newly developed smoothed point interpolation method (S‐PIM) can handle the largely deformable solids owing to its softened model stiffness and insensitivity to mesh distortion. In this work, a novel coupled method has been proposed by combining the immersed boundary‐lattice Boltzmann method with the S‐PIM for fluid‐structure interaction problems with large‐displacement solids. The proposed method preserves the simplicity of the lattice Boltzmann method for fluid solvers, utilizes the S‐PIM to establish the realistic constitutive laws for nonlinear solids, and avoids mesh regeneration based on the frame of the immersed boundary method. Both two‐ and three‐dimensional numerical examples have been carried out to validate the accuracy, convergence, and stability of the proposed method in consideration of comparative results with referenced solutions.     

A novel finite point method for flow simulation

A novel finite point method is developed to simulate flow problems. The mashes in the traditional numerical methods are supplanted by the distribution of points in the calculation domain. A local interpolation based on the properties of Taylor series expansion is used to construct an approximation for unknown functions and their derivatives. An upwind‐dominated scheme is proposed to efficiently handle the non‐linear convection. Comparison with the finite difference solutions for the two‐dimensional driven cavity flow and the experimental results for flow around a cylinder shows that the present method is capable of satisfactorily predicting the flow separation characteristic. The present algorithm is simple and flexible for complex geometric boundary. The influence of the point distribution on computation time and accuracy of results is included. Copyright © 2002 John Wiley & Sons, Ltd.     

Solving time‐dependent PDEs using the material point method,a case study from gas dynamics

The material point method (MPM) developed by Sulsky and colleagues is currently being used to solve many challenging problems involving large deformations and/or fragementations with some success. In order to understand the properties of this method, an analysis of the considerable computational properties of MPM is undertaken in the context of model problems from gas dynamics. The MPM method in the form used here is shown both theoretically and computationally to have first‐order accuracy for a standard gas dynamics test problem. Copyright © 2009 John Wiley & Sons, Ltd.