多相流局部混合型质点网格法

提出模拟多相流的局部混合型质点网格法，该方法能稳定地模拟高浓度流体颗粒两相 流. 在每一个颗粒团尺度的欧拉网格下(本文称之为欧拉微元)，基于Lagrangian 追踪原理，可直接估计体积内颗粒的总量，从而准确求出欧拉微元和控制容积内颗粒的 浓度(即颗粒在容积内的体积含量). 同时，假设在新的时间步下，颗粒在欧拉微元里充分 混合，形成新颗粒团. 作者对竖直和倾斜容器中单粒径颗粒沉降和竖直容器中双粒径颗粒双 峰悬浮液沉降过程进行了计算，结果与实测数据相符.     

高速冲击过程数值分析的再生核质点法

利用新型的无网格方法-再生核质点方法对高速冲击过程进行数值模拟，给出了控制 方程，分析中引入Bordner-Partom本构模型来实现材料高速冲击条件下的大应变和高应变 率的特性，提出了新的接触面处理方法，介绍了速度配点法来实现接触面及其它本质边界的 速度条件. 数值模拟表明，再生核质点法能够准确模拟高速冲击过程，具有较高的计算精度.     

节理岩体的爆破松动机理

利用应力波理论对节理岩体的爆破松动机理进行了探讨,同时利用动力有限元数值模拟方法验证了理论分析结果。研究表明,爆破荷载的加载阶段对岩体的挤压蓄能和快速卸载阶段的弹性回复是岩体产生后续松动位移的根本原因;节理岩体的爆破松动在开挖面的法向表现为岩体的剪切错动,在平行于开挖面方向表现为结构面的张开;拉格朗日坐标下的动力有限元数值方法,如ANSYS/LS-DYNA程序,能够较好地用来模拟节理岩体的爆破松动过程。     

基于隐式扩散的直接力格式浸没边界格子Boltzmann方法

采用浸没边界格子Boltzmann(immersed boundary-lattice Boltzmann,IB-LB)模型执行动边界绕流数值模拟时,信息交互界面和边界力计算格式直接影响流动求解器的数值精度和计算效率.基于隐式扩散界面,一种改进的直接力格式IB-LB模型被提出.边界力表达式基于欧拉/拉格朗日变量同一性准...     

碎冰区海冰与船舶结构相互作用的离散元分析

采用离散元模型对碎冰区海冰与船舶结构的相互作用进行了数值模拟。碎冰由三维圆盘单元构成,并考虑了其在海流作用下的浮力、拖曳力、附加质量;采用单元组合的方式构造了船体的简易模型。通过海冰与船体单元间的接触判断和接触力计算确定海冰与船体结构之间的相互作用。采用以上离散单元模型对不同冰速、冰厚、冰块尺寸、密集度条件下海冰与船体的动力作用过程进行了数值计算,并对比分析了不同海冰参数下冰对船体作用力的幅值。结果表明:冰对船体的作用力随冰厚、流速、密集度、冰块尺寸的增加而增大。本文工作可为冰区船舶的安全运行和结构设计提供参考。     

模拟结构倒塌破坏的有限质点法

有限质点法是以向量式力学为基础的新兴结构分析方法,本文将其应用于冲击荷载作用下的网壳结构的倒塌破坏模拟中。以空间杆单元为例建立了有限质点法的基本方程,推导了求解几何和材料非线性问题的基本公式。为计算断裂问题,建立了空间杆单元的断裂准则和断裂模型,发展了有限质点法进行断裂分析的基本算法。通过对某双层网壳冲击荷载下破坏过程的模拟和分析,验证了该方法在结构倒塌破坏过程模拟中的有效性和适用性。     

海冰动力学中本构模型研究的若干进展

从20世纪70年代, 人们在不同尺度下建立了一系列的海冰动力学本构模型用于海冰数值模拟和预测. 将目前应用的海冰动力学本构模型分为弹塑性、黏塑性、各向异性和颗粒流体动力学中的黏弹塑性模型4类,并分别讨论了各类模型的特点和适用范围. 尤其对在极区及副极区大、中尺度下广泛应用的黏塑性及其改进的本构模型进行了深入地分析. 最后得出: 在大、中尺度下建立海冰动力学的黏弹塑性本构模型是提高海冰动力学计算精度的有效途径; 将小尺度下采用黏弹塑性本构模型的海冰颗粒流体动力学的计算时效进一步改进后, 可在大、中尺度下对海冰的重叠、堆积特性进行模拟; 进一步开展海冰动力学的尺度效应研究, 进行海冰本构模型的实验验证, 并建立不同尺度模型间的相互联系也应是海冰动力学本构模型研究的重要内容.      

无网格MPM法三维爆炸焊接数值模拟

物质点法MPM(Material Point Method)是无网格方法之一.它是在质点网格法(PIC)基础上发展而来的一种新数值方法,它利用了欧拉法和拉格朗日法两者的优点,计算物质点在冲击载荷下的应力和应交;通过物质点来跟踪材料体的变形和破损,而在整个计算过程中背景网格始终固定不变,避免了重新划分网格.本文应用MPM法计算三维爆炸焊接问题,在爆轰载荷作用下的飞板和基板的金属动态变形过程进行了三维数值模拟,并且对飞板的碰撞点速度和爆轰压力变化进行了计算分析.     

薄膜断裂和穿透的向量式有限元分析及应用

向量式有限元基于牛顿运动定律,通过质点描述和向量分析来求解整体结构的动力响应。首先,给出了向量式有限元三角形薄膜单元的基本理论,进而针对薄膜结构的断裂和穿透破坏过程,提出相应解决方案。对于薄膜断裂问题,采用Mises应力状态变量达到失效应力限值作为断裂判据,通过质点分裂方式将相连单元的对应节点断开,并对分裂后的新质点进行状态更新来模拟其断裂过程;对于薄膜穿透问题,则同时结合碰撞和断裂过程模拟来实现。在此基础上,编制了薄膜结构的断裂和穿透求解计算程序,算例分析表明,程序可很好地完成薄膜结构的大变形大转动、断裂和穿透等不连续行为的模拟,验证了理论及程序的可靠性和有效性,体现了本文方法进行薄膜结构复杂不连续行为分析的优势。     

碎冰区海冰与船舶结构相互作用的离散元分析

采用离散元模型对碎冰区海冰与船舶结构的相互作用进行了数值模拟。碎冰由三维圆盘单元构成,并考虑了其在海流作用下的浮力、拖曳力、附加质量;采用单元组合的方式构造了船体的简易模型。通过海冰与船体单元间的接触判断和接触力计算确定海冰与船体结构之间的相互作用。采用以上离散单元模型对不同冰速、冰厚、冰块尺寸、密集度条件下海冰与船体的动力作用过程进行了数值计算,并对比分析了不同海冰参数下冰对船体作用力的幅值。结果表明：冰对船体的作用力随冰厚、流速、密集度、冰块尺寸的增加而增大。本文工作可为冰区船舶的安全运行和结构设计提供参考。     

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

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

Cahn–Hilliard modeling of particles suspended in two‐phase flows

In this paper, we present a model for the dynamics of particles suspended in two‐phase flows by coupling the Cahn–Hilliard theory with the extended finite element method (XFEM). In the Cahn–Hilliard model the interface is considered to have a small but finite thickness, which circumvents explicit tracking of the interface. For the direct numerical simulation of particle‐suspended flows, we incorporate an XFEM, in which the particle domain is decoupled from the fluid domain. To cope with the movement of the particles, a temporary ALE scheme is used for the mapping of field variables at the previous time levels onto the computational mesh at the current time level. By combining the Cahn–Hilliard model with the XFEM, the particle motion at an interface can be simulated on a fixed Eulerian mesh without any need of re‐meshing. The model is general, but to demonstrate and validate the technique, here the dynamics of a single particle at a fluid–fluid interface is studied. First, we apply a small disturbance on a particle resting at an interface between two fluids, and investigate the particle movement towards its equilibrium position. In particular, we are interested in the effect of interfacial thickness, surface tension, particle size and viscosity ratio of two fluids on the particle movement towards its equilibrium position. Finally, we show the movement of a particle passing through multiple layers of fluids. Copyright © 2011 John Wiley & Sons, Ltd.     

Coupling Eulerian‐Lagrangian method of air‐particle two‐phase flow with population balance equations to simulate the evolution of vehicle exhaust plume

In this paper, we present a new numerical scheme to describe the dynamic evolution of multiphase polydisperse systems in terms of time, space, and properties by coupling the Eulerian‐Lagrangian method for air‐particle two‐phase flow and population balance equations to describe particle property evolution due to microbehaviors (eg, aggregation, breakage, and growth). This coupling scheme was used to comprehensively simulate the two‐phase flow structure, particle size spectrum, particle number, and volume concentrations. These were characterized by a high‐resolution particle tracking using the Lagrangian approach and the high precision of moments of the particle size spectrum by solving the population balance equation with the quadrature method of moments. The algorithm of the coupling scheme was incorporated into the open source computational fluid dynamics software OpenFOAM to simulate the dynamic evolution of vehicle exhaust plume. The impacts of vehicle velocity, exhaust temperature, and aggregation efficiency on the distribution of auto exhaust particles in space and changes in their properties were analyzed. The results indicate that the particle number concentration, volume concentration, and average diameter of particles in the vehicle exhaust plume could be strongly affected by the plume structure and flow properties.     

Averaging of particle data from phase Doppler anemometry in unsteady two-phase flow: Validation by numerical simulation

A novel post-processing algorithm is proposed to correct statistical bias observed in the treatment of time series obtained by a phase Doppler anemometer (PDA) at flow locations with variable particle velocity and concentration. Extensive properties of each validated particle are weighted with their inverse measuring (validation) volume to account for the procedure of particle sampling and fluctuations in the particle concentration. To compensate for the short characteristic length of the validation volume, the properties of particles are expressed by properties of fields of particle groups, using a local averaging time. A window shift and a decorrelation scheme are applied on the fields to increase their frequency resolution. This algorithm has been tested on numerical time series, provided by an Eulerian/Lagrangian code representing a gas/solids flow past a bluff body. Moments and spectral estimates of concentration and velocity of particle groups were successfully validated by the numerical simulation using the PDA data algorithm and control volume averaging. The control volume was much larger than the PDA validation volume, but the centre positions of the two volumes were identical.     

Two-phase micro- and macro-time scales in particle-laden turbulent channel flows

The micro-and macro-time scales in two-phaseturbulent channel flows are investigated using the direct numerical simulation and the Lagrangian particle trajectorymethods for the fluid-and the particle-phases,respectively.Lagrangian and Eulerian time scales of both phases are calculated using velocity correlation functions.Due to flowanisotropy,micro-time scales are not the same with the theoretical estimations in large Reynolds number(isotropic) turbulence.Lagrangian macro-time scales of particle-phaseand of fluid-phase seen by particles are both dependent onparticle Stokes number.The fluid-phase Lagrangian integral time scales increase with distance from the wall,longerthan those time scales seen by particles.The Eulerian integral macro-time scales increase in near-wall regions but decrease in out-layer regions.The moving Eulerian time scalesare also investigated and compared with Lagrangian integraltime scales,and in good agreement with previous measurements and numerical predictions.For the fluid particles themicro Eulerian time scales are longer than the Lagrangianones in the near wall regions,while away from the walls themicro Lagrangian time scales are longer.The Lagrangianintegral time scales are longer than the Eulerian ones.Theresults are useful for further understanding two-phase flowphysics and especially for constructing accurate predictionmodels of inertial particle dispersion.     

Numerical investigation of ice breaking by a high-pressure bubble based on a coupled BEM-PD model

In this paper, by combining the boundary element method (BEM) and peridynamics (PD), a bubble-ice interaction model is established, which can investigate the dynamic interactions between a high-pressure bubble and an ice plate with particular focus on the mechanical behaviors of ice breaking. The bubble dynamics are solved by BEM based on the potential flow theory. Ice cracks initiation and propagation are simulated by the bond-based peridynamics which is validated by a three-point bending test. The fluid–structure interaction (FSI) is achieved by matching the normal velocity and hydrodynamic loads at the fluid–structure interface. To validate the proposed FSI model, an experiment is carried out in which an oscillating bubble is generated under an ice plate by underwater discharge system. The whole interaction process is captured by a Phantom V711 high-speed camera. Qualitative agreements are achieved between the numerical and experimental results. The underlying mechanism of cracks initiation, propagation, branching, and coalescence of the ice plate is found to highly depend on three parameters, i.e., bubble–ice distance, ice thickness and bubble size. The present study is expected to provide further assists in the understanding of ice breaking problems.     

Micro-rheology of dense particulate flows: Application to immersed avalanches

We rely here on a non-smooth contact dynamics (NSCD) approach to treat particle collisions in a direct numerical simulation of a dense particulate flow. Interactions between particles are considered by a non-smooth formulation of particle dynamics at the microscopic scale, which enables one to straightforwardly implement complex contact laws. The hydrodynamic coupling is achieved by a distributed Lagrange multiplier/fictitious domain (DLM/FD) method. As a preliminary step, the relevance of our NSCD-DLM/FD method is assessed by comparing results of 2D sedimentation simulations with those obtained with a usual molecular dynamics collision model. Then, we use it to investigate how a fully immersed granular packing collapses depending on its initial particle volume fraction, providing clues on the micro-rheology of dense particulate flows.     

Evaluation of the equilibrium Eulerian approach for the evolution of particle concentration in isotropic turbulence

In the current work, the accuracy of the equilibrium Eulerian approach in evolving the particulate concentration field is evaluated by comparing it against the Lagrangian approach, for varying particle response time and terminal velocity. In particular, we compare the statistics of preferential accumulation and gravitational settling of particles in a cubic box of isotropic turbulence. Twelve simulations corresponding to four values of nondimensional particle response time, τ_p=0.05, 0.1, 0.2, 0.4, and three values of nondimensional terminal velocity, |V_s|=0.5,2,4 are considered. The equilibrium Eulerian approach obviates the need to solve additional governing equations for the particle velocity field. It, however, involves evolution of the particle concentration field using the equilibrium Eulerian velocity field. A spectral diffusion term is included in the particle concentration equation to provide an essentially non-oscillatory behavior to the solution. There is good agreement between the equilibrium Eulerian and Lagrangian statistics for small particles. With increasing particle size, the equilibrium Eulerian approach tends to somewhat overestimate particle preferential concentration in regions of excess strain-rate over rotation-rate compared to the Lagrangian approach. Over the entire range of parameters considered, the equilibrium approach provides a good approximation to the actual mean and rms fluctuating settling velocities of the particle.     

A three‐dimensional vortex particle‐in‐cell method for vortex motions in the vicinity of a wall

A new vortex particle‐in‐cell method for the simulation of three‐dimensional unsteady incompressible viscous flow is presented. The projection of the vortex strengths onto the mesh is based on volume interpolation. The convection of vorticity is treated as a Lagrangian move operation but one where the velocity of each particle is interpolated from an Eulerian mesh solution of velocity–Poisson equations. The change in vorticity due to diffusion is also computed on the Eulerian mesh and projected back to the particles. Where diffusive fluxes cause vorticity to enter a cell not already containing any particles new particles are created. The surface vorticity and the cancellation of tangential velocity at the plate are related by the Neumann conditions. The basic framework for implementation of the procedure is also introduced where the solution update comprises a sequence of two fractional steps. The method is applied to a problem where an unsteady boundary layer develops under the impact of a vortex ring and comparison is made with the experimental and numerical literature. Copyright © 2001 John Wiley & Sons, Ltd.