MPS-FEM数值分析带自由面的流固耦合问题

随着计算科学的发展,研究人员为探索流固耦合问题的物理机理而提出了众多的数值方法。其中,耦合的移动粒子半隐式方法 MPS(Moving Particle Semi-Implicit method)和有限单元法FEM(Finite Element method)为流固耦合问题的数值仿真工作提供了新的途径。本文所有流场的数值模拟工作均采用课题组自主开发的无网格法求解器MLParticle-SJTU来完成。该求解器在原始的MPS法基础上,对核函数、压力梯度模型、压力泊松方程的求解和自由面判断方式等方面进行了改进。此外,在该求解器框架内,基于FEM法拓展了针对结构场进行求解的功能。首先,对MPS和FEM方法的理论模型及其耦合策略进行了介绍。然后,采用该自研MPS-FEM耦合求解器,数值模拟了溃坝流动对弹性结构的冲击及其相互作用的标准问题。通过将结构变形及自由面波型变化等结果与已发表结果进行对比,验证了该求解器在处理带自由面剧烈变化的粘性流体和柔性变形结构的耦合作用问题上的可行性。     

Effect of coarse particles on the heat transfer in a particle-laden turbulent boundary layer

The temperature distribution in particle-laden turbulent flow, in a flume, was investigated both by DNS and experimentally. Simulations were performed at Re=171 and Pr=5.4 in order to study the interaction between the particle motion and flow turbulence. Two-way coupling was used to obtain various turbulence statistics, the grid resolution was sufficiently fine to resolve all essential turbulent scales. The effect of particle diameter on momentum, heat transfer and particle deposition was considered. The details of particle-turbulence interaction depend on the particle Stokes number and the particle Reynolds number.

The spatial structures of instantaneous flow and temperature fields were visualized. Low frequency small oscillations of deposited particles were observed. It was found that these small deviations from the initial position, caused strong changes in the instantaneous temperature field near the particle.

The experiments provided details of the temperature field on the heated wall close to the particle. In the front of the particle, a sharp increase in heat transfer coefficient was observed. The experimental results agree well with the computational predictions.     

Mixed convective heat transfer characteristics and mechanisms in structured packed beds

The structured packed bed is considered a promising reactor owing to its low pressure drop and good heat transfer performance. In the heat transfer process of thermal storage in packed beds, natural convection plays an important role. To obtain the mixed convective heat transfer characteristics and mechanisms in packed beds, numerical simulations and coupling analyses were carried out in this study on the unsteady process of fluid flow and heat transfer. A three-dimensional model of the flow channel in the packed bed was established, and the Navier–Stokes equations and Laminar model were adopted for the computations. The effects of the driving force on fluid flow around a particle were studied in detail. The differences in velocity and density distributions under different flow directions due to effect of the aiding flow or opposing flow were intuitively demonstrated and quantitatively analyzed. It was found that the driving force strengthens the fluid flow near the particle surface when aiding flow occurs and inhibits the fluid flow when opposing flow occurs. The boundary layer structure was changed by the natural convection, which in turn influences the field synergy angle. For the aiding flow, the coordination between the velocity and density fields is higher than that for the opposing flow. By analysis the effects of physical parameters on mixed convective heat transfer, it is indicated that with an increase in the fluid-solid temperature difference or the particle diameter, or a decrease in the fluid temperature, the strengthening or inhibiting effect of natural convection on the heat transfer became more significant.     

Micro/meso simulation of a fluidized bed in a homogeneous bubbling regime

Due to the wide range of spatial scales and the complex features associated to fluid/solid and solid/solid interactions in a dense fluidized bed, the system can be studied at different length scales, namely micro, meso and macro. In this work, we select a flow configuration relevant of a homogeneous liquid/solid fluidization and compare computed results from Particle Resolved Simulation (PRS) with those from locally averaged Euler/Lagrange simulation. PRS at the micro-scale is carried out by a parallel Distributed Lagrange Multiplier (DLM) solver in the framework of fictitious domain methods (Wachs, 2011a, 2015). For meso-scale simulations, the set of mass and momentum conservation equations is averaged in control volumes encompassing few particles and momentum transfer between the two phases is modeled using appropriate drag laws. Both methods are coupled to a Discrete Element Method (DEM) combined with a soft-sphere contact model to solve the Newton–Euler equations with collisions for the particles in a Lagrangian framework (Wachs et al., 2012). A test case of intermediate size with 2000 spheres is chosen as a sensible compromise between size limitations of the meso-scale model for an appropriate averaging process and computational resources required to run micro-scale simulations. These two datasets yield new insight on momentum transfer at different spatial scales in the flow, and question the validity of certain approximations adopted in the meso-scale model. Results demonstrate an acceptable agreement between the micro- and meso-scale predictions on integral measures as pressure drop and bed height. Investigating more detailed features of the flow, it has been shown that particles fluctuations are considerably suppressed in meso-scale simulations and in particular the particles transverse motion is underestimated, regardless of the selected drag law. The origin of these dependencies is carefully investigated by reconstructing the closure laws based on PRS results and comparing them to the closure laws proposed in the literature.     

Effect of inertial particles with different specific heat capacities on heat transfer in particle-laden turbulent flow

The effect of inertial particles with different specific heat on heat transfer in particle-laden turbulent channel flows is studied using the direct numerical simulation(DNS) and the Lagrangian particle tracking method. The simulation uses a two-way coupling model to consider the momentum and thermal interactions between the particles and turbulence. The study shows that the temperature fields display differences between the particle-laden flow with different specific heat particles and the particle-free flow,indicating that the particle specific heat is an important factor that affects the heat transfer process in a particle-laden flow. It is found that the heat transfer capacity of the particle-laden flow gradually increases with the increase of the particle specific heat. This is due to the positive contribution of the particle increase to the heat transfer. In addition,the Nusselt number of a particle-laden flow is compared with that of a particle-free flow.It is found that particles with a large specific heat strengthen heat transfer of turbulent flow, while those with small specific heat weaken heat transfer of turbulent flow.     

Dealing with extremely large deformation by Nearest-Nodes FEM with algorithm for updating element connectivity

In the recently developed Nearest-Nodes Finite Element Method (NN-FEM), elements are mainly used for numerical integration; while shape functions are constructed in a similar way as in meshless methods. Based on this strategy, NN-FEM inherits major merits from both the classical Finite Element Method and meshless methods. One of them is that NN-FEM is nearly not affected by element distortion. So NN-FEM is more efficient than the classical FEM on dealing with large deformation problems. Nevertheless, NN-FEM still has a requirement on finite element meshes, that is, elements in a mesh are required not to overlap or penetrate to each other, to avoid difficulty in numerical integration. To eliminate overlapped elements, NN-FEM is supplemented with an algorithm for updating element connectivity. With this supplement, NN-FEM is able to deal with extremely large deformation. In updating element connectivity, element nodes are kept not changed and all information associated with nodes are not touched. Therefore, there is no need to transfer solution data, and error introduced by solution transfer is avoided.     

Design sensitivity analysis for the homogenized elasticity tensor of a polymer filled with rubber particles

The main purpose of this work is the computational simulation of the sensitivity coefficients of the homogenized tensor for a polymer filled with rubber particles with respect to the material parameters of the constituents. The Representative Volume Element (RVE) of this composite contains a single spherical particle, and the composite components are treated as homogeneous isotropic media, resulting in an isotropic effective homogenized material. The sensitivity analysis presented in this paper is performed via the provided semi-analytical technique using the commercial FEM code ABAQUS and the symbolic computation package MAPLE. The analytical method applied for comparison uses the additional algebraic formulas derived for the homogenized tensor for a medium filled with spherical inclusions, while the FEM-based technique employs the polynomial response functions recovered from the Weighted Least-Squares Method. The homogenization technique consists of equating the strain energies for the real composite and the artificial isotropic material characterized by the effective elasticity tensor. The homogenization problem is solved using ABAQUS by the application of uniform deformations on specific outer surfaces of the composite RVE and the use of tetrahedral finite elements C3D4. The energy approach will allow for the future application of more realistic constitutive models of rubber-filled polymers such as that of Mullins and for RVEs of larger size that contain an agglomeration of rubber particles.     

Determination of contact parameters for discrete element method simulations of granular systems

Both linear-spring-dashpot (LSD) and non-linear Hertzian-spring-dnshpot (HSD) contact models are commonly used for the calculation of contact forces in Discrete Element Method (DEM) simulations of granular systems.Despite the popularity of these models, determination of suitable values for the contact parameters of the simulated particles such as stiffness, damping coefficient, coefficient of restitution, and simulation time step,is not altogether obvious.In this work the relationships between these contact parameters for a model system where a particle impacts on a flat base are examined.Recommendations are made concerning the determination of these contact parameters for use in DEM simulations.     

气固两相流介尺度LBM-DEM模型

LBM-DEM耦合方法通常是指一种颗粒流体系统直接数值模拟算法,即是一种不引入经验曳力模型的计算方法,颗粒尺寸通常比计算网格的长度大一个量级,颗粒的受力通过表面的粘性力与压力积分获得,其优点是能描述每个颗粒周围的详细流场,产生详细的颗粒-流体相互作用的动力学信息,可以探索颗粒流体界面的流动、传递和反应的详细信息及两相相互作用的本构关系,但其缺点是计算量巨大,无法应用于真实流化床过程模拟。本文针对气固流化床中的流体以及固体颗粒间的多相流体力学行为,建立了一种稠密气固两相流的介尺度LBMDEM模型,即LBM-DEM耦合的离散颗粒模型,实现在颗粒尺度上流化床的快速离散模拟。该耦合模型采用格子玻尔兹曼方法(LBM)描述气相的流动和传递行为,离散单元法(DEM)用于描述颗粒相的运动,并利用能量最小多尺度(EMMS)曳力解决气固耦合不成熟问题,以提高其模拟精度。通过经典快速流态化的模拟,验证了介尺度LBM-DEM耦合模型的有效性。模拟结果表明介尺度LBM-DEM模型是一种探索实验室规模气固系统的有力手段。     

A thermo-mechanically coupled finite strain model considering inelastic heat generation

The procedure for reuse of finite element method (FEM) programs for heat transfer and structure analysis to solve advanced thermo-mechanical problems is presented as powerful algorithm applicable for coupling of other physical fields (magnetic, fluid flow, etc.). In this case, nonlinear Block-Gauss–Seidel partitioned algorithm strongly couples the heat transfer and structural FEM programs by a component-based software engineering. Component template library provides possibility to exchange the data between the components which solve the corresponding subproblems. The structural component evaluates the dissipative energy induced by inelastic strain. The heat transfer component computes the temperature change due to the dissipation. The convergence is guaranteed by posing the global convergence criterion on the previously locally converged coupled variables. This enables reuse of software and allows the numerical simulation of thermo-sensitive problems.     

Numerical simulation of a dense solid particle flow inside a cyclone separator using the hybrid Euler–Lagrange approach

This paper presents a numerical simulation of the flow inside a cyclone separator at high particle loads. The gas and gas–particle flows were analyzed using a commercial computational fluid dynamics code. The turbulence effects inside the separator were modeled using the Reynolds stress model. The two phase gas–solid particles flow was modeled using a hybrid Euler–Lagrange approach, which accounts for the four-way coupling between phases. The simulations were performed for three inlet velocities of the gaseous phase and several cyclone mass particle loadings. Moreover, the influences of several submodel parameters on the calculated results were investigated. The obtained results were compared against experimental data collected at the in-house experimental rig. The cyclone pressure drop evaluated numerically underpredicts the measured values. The possible reason of this discrepancies was disused.     

NUMERICAL MODELING FOR FLUID-STRUCTURE INTERACTION UNDER BLAST AND IMPACT LOADING RESPONSE

建立了基于位移相等条件的流固耦合数值模拟程序框架,通过串联弱耦合方式对爆炸冲击波作用下流固耦合效应进行数值模拟. 其中非线性固体位移场采用基于Lagrange方法描述的时域间断伽辽金有限元方法进行处理. 基于修正弹簧近似的非结构动网格新技术,非定常流场采用格心格式的有限体积方法进行求解. 数值模拟结果表明: 该文所发展的弱耦合分析程序在流固耦合数值模拟过程中具有可靠的计算精度. 同时,程序对于爆炸强间断气动冲击载荷冲击作用下固体变形引起流场的反射和叠加效应, 以及流场变化引起复杂的固体响应具有良好的耦合求解能力.     

近场动力学与有限元方法耦合求解热传导问题

求解含裂纹等不连续问题一直是计算力学的重点研究课题之一,以偏微分方程为基础的连续介质力学方法处理不连续问题时面临很大的困难. 近场动力学方法是一种基于积分方程的非局部理论,在处理不连续问题时有很大的优越性. 本文提出了求解含裂纹热传导问题的一种新的近场动力学与有限元法的耦合方法. 结合近场动力学方法处理不连续问题的优势以及有限元方法计算效率高的优势,将求解区域划分为两个区域,近场动力学区域和有限元区域. 包含裂纹的区域采用近场动力学方法建模,其他区域采用有限元方法建模. 本文提出的耦合方案实施简单方便,近场动力学区域与有限元区域之间不需要设置重叠区域. 耦合方法通过近场动力学粒子与其域内所有粒子(包括近场动力学粒子和有限元节点)以非局部方式连接,有限元节点与其周围的所有粒子以有限元方式相互作用. 将有限元热传导矩阵和近场动力学粒子相互作用矩阵写入同一整体热传导矩阵中,并采用Guyan缩聚法进一步减小计算量. 分别采用连续介质力学方法和近场动力学方法对一维以及二维温度场算例进行模拟,结果表明,本文的耦合方法具有较高的计算精度和计算效率. 该耦合方案可以进一步拓展到热力耦合条件下含裂纹材料和结构的裂纹扩展问题.      

Flow and temperature field measurements of thermal convection in a small vertical gap using liquid crystals

Thermal convection in a small vertical gap is studied experimentally applying digital particle image velocimetry/thermometry. This optical method enables the simultaneous measurement of two-dimensional flow and temperature fields in a liquid. The principle is based on seeding the liquid flow medium with thermochromic liquid crystal particles. The temperature is measured by the crystal particles which change their reflected colour as function of temperature. The flow velocity is measured by using the same particles as flow tracers. The investigation shall contribute to the understanding of the fluid mechanical behaviour of biological liquids within micro reactor systems. However, the problem is also of fundamental interest as far as heat and mass transfer is concerned. Measured temperature and flow velocity fields are presented and discussed. Presented in part at the 4th Chemnitz/Hamburger Colloquium (CHC) on Microflows, Hamburg, Germany, November 2004.     

Accurate numerical estimation of interphase momentum transfer in Lagrangian–Eulerian simulations of dispersed two-phase flows

The Lagrangian–Eulerian (LE) approach is used in many computational methods to simulate two-way coupled dispersed two-phase flows. These include averaged equation solvers, as well as direct numerical simulations (DNS) and large-eddy simulations (LES) that approximate the dispersed-phase particles (or droplets or bubbles) as point sources. Accurate calculation of the interphase momentum transfer term in LE simulations is crucial for predicting qualitatively correct physical behavior, as well as for quantitative comparison with experiments. Numerical error in the interphase momentum transfer calculation arises from both forward interpolation/approximation of fluid velocity at grid nodes to particle locations, and from backward estimation of the interphase momentum transfer term at particle locations to grid nodes. A novel test that admits an analytical form for the interphase momentum transfer term is devised to test the accuracy of the following numerical schemes: (1) fourth-order Lagrange Polynomial Interpolation (LPI-4), (3) Piecewise Cubic Approximation (PCA), (3) second-order Lagrange Polynomial Interpolation (LPI-2) which is basically linear interpolation, and (4) a Two-Stage Estimation algorithm (TSE). A number of tests are performed to systematically characterize the effects of varying the particle velocity variance, the distribution of particle positions, and fluid velocity field spectrum on estimation of the mean interphase momentum transfer term. Numerical error resulting from backward estimation is decomposed into statistical and deterministic (bias and discretization) components, and their convergence with number of particles and grid resolution is characterized. It is found that when the interphase momentum transfer is computed using values for these numerical parameters typically encountered in the literature, it can incur errors as high as 80% for the LPI-4 scheme, whereas TSE incurs a maximum error of 20%. The tests reveal that using multiple independent simulations and higher number of particles per cell are required for accurate estimation using current algorithms. The study motivates further testing of LE numerical methods, and the development of better algorithms for computing interphase transfer terms.     

Euler–Lagrange/DEM simulation of wood gasification in a bubbling fluidized bed reactor

We present an Euler–Lagrange method for the simulation of wood gasification in a bubbling fluidized bed. The gas phase is modeled as a continuum using the 2D Navier–Stokes equations and the solid phase is modeled by a Discrete Element Method (DEM) using a soft-sphere approach for the particle collision dynamic. Turbulence is included via a Large-Eddy approach using the Smagorinsky sub-grid model. The model takes into account detailed gas phase chemistry, zero-dimensional modeling of the pyrolysis and gasification of each individual particle, particle shrinkage, and heat and mass transfer between the gas phase and the particulate phase. We investigate the influence of wood feeding rate and compare exhaust gas compositions and temperature results obtained with the model against experimental data of a laboratory scale bubbling fluidized bed reactor.     

On the numerical convergence and performance of different spatial discretization techniques for transient elastodynamic wave propagation problems

We present a systematic investigation of several discretization approaches for transient elastodynamic wave propagation problems. This comparison includes a Finite Difference, a Finite Volume, a Finite Element, a Spectral Element and the Scaled Boundary Finite Element Method. Numerical examples are given for simple geometries with normalized parameters, for heterogeneous materials as well as for structures with arbitrarily shaped material interfaces. General conclusions regarding the accuracy of the methods are presented. Based on the essential numerical examples an expansion of the results to a wide range of problems and thus to numerous fields of application is possible.     

Influence of heat transfer and fluid flow on crack growth in multilayered porous/dense materials using XFEM: Application to Solid Oxide Fuel Cell like material design

An advanced numerical model is developed to investigate the influence of heat transfer and fluid flow on crack propagation in multi-layered porous materials. The fluid flow, governed by the Navier–Stokes and Darcy’s law, is discretized with the nonconforming Crouzeix–Raviart (CR) finite element method. A combination of Discontinuous Galerkin (DG) and Multi-Point Flux Approximation (MPFA) methods is used to solve the advection–diffusion heat transfer equation in the flow channel and in the fluid phase within the porous material. The crack is assumed to affect only the heat diffusion within the porous layer, therefore a time splitting technique is used to solve the heat transfer in the fluid and the solid phases separately. Thus, within the porous material, the crack induces a discontinuity of the temperature at the crack surfaces and a singularity of the flux at the crack tip. Conduction in the solid phase is solved using the eXtended Finite Element Method (XFEM) to better handle the discontinuities and singularities caused by the cracks. The XFEM is also used to solve the thermo-mechanical problem and to track the crack propagation. The multi-physics model is implemented then validated for the transient regime, this necessitated a post processing treatment in which, the stress intensity factors (SIF) are computed for each time step. The SIFs are then used in the crack propagation criterion and the crack orientation angle. The methodology seems to be robust accurate and the computational cost is reduced thanks to the XFEM.     

Effect of temperature gradient within a solid particle on the rotation and oscillation modes in solid-dispersed two-phase flows

Liquid–solid two-phase flow with heat transfer is simulated, and the effect of temperature gradient within a solid particle on the particle behaviour and heat transfer is studied. The interaction between fluid and particles is considered with our original immersed solid approach on a rectangular grid system. The local heat flux at the fluid–solid interface is described with an anisotropic heat conductivity matrix, and the governing equation of temperature is time-updated with an implicit treatment for the diffusion term. The method is applied to a 2-D natural convection flow of a relatively low Rayleigh number including multiple particles. Heat transfer and particle behaviours are studied for different solid heat conductivities (ratio to the fluid conductivity ranging between 10⁻³ and 10³) and solid volume fractions. Under a condition of relatively low heat conductivity ratio, the particles show a simple circulating flow. By increasing the heat conductivity ratio, a transition of the particulate flow is observed to oscillation mode around the domain centre due to the buoyancy force as a restitution force. The oscillation period is found to vary with the heat conductivity ratio, and it is related to the time scales for the heat transfer via fluid and solid.