带嵌件型腔内熔接过程的数值模拟研究

基于充模过程的两相黏弹性流体模型, 采用同位网格有限体积法, 结合浸入边界法和界面追踪的复合水平集流体体积方法实现了带嵌件型腔内充模过程的动态模拟. 基于上述模型和算法模拟了熔体前沿界面及熔接线的动态演化过程, 而且通过线性应力-光学定律得到了熔接线附近的流动诱导应力分布情况; 讨论了熔体温度及模具温度对熔接线区域凝固层厚度的影响. 数值结果表明: 本文提出的方法可用于模拟复杂型腔内的充模过程以及熔接线的自动追踪; 由于聚合物黏弹性熔体流动的复杂性, 当两股熔体相遇后, 熔接线不同位置的应力分布变化较大; 熔体或模具温度越高, 熔接线区域凝固层厚度越薄, 提高熔体或模具温度能够改善甚至消除充模过程中的熔接线.     

黏弹性熔体充模流动诱导残余应力模拟

流动诱导残余应力是塑料制品产生应力开裂以及翘曲变形等现象的重要原因,对成型过程中流动诱导残余应力研究具有重要意义.推导了基于黏弹性eXtended Pom-Pom本构关系的能量方程,进而建立了描述黏弹性流体非等温充模流动的气-液两相模型.用同位网格有限体积法进行了求解,得到了凝固层和剪切速率分布,给出了充填结束时影响制件力学性能的流动诱导残余应力.结果表明,型腔中凝固层的厚度与注射速率有关,注射速率越大,充模时间越短,凝固层越薄.在制品表层紧邻模壁的地方,剪切速率和残余应力几乎为零;在制品次表层的位置,制件内剪切速率和流动残余应力也较高;而在远离模壁的地方,剪切速率和流动残余应力也较小.     

复杂型腔充模过程中分子构型演化的数值模拟

基于充模过程的两相黏弹性流体模型,采用有限体积、浸入边界和复合水平集流体体积方法,数值模拟了聚合熔体在复杂型腔中的充模过程.首先,借助一类特殊函数(R-functions)将基于基本几何体的水平集函数组合成描述复杂型腔的形状水平集函数.然后,采用浸入边界法处理复杂型腔问题,有限体积方法求解熔体控制方程,利用复合水平集流体体积方法对熔体前沿界面进行隐式追踪.基于有限伸展非线性弹性哑铃本构方程模型,计算熔体分子构型张量,通过取向椭圆描述分子的取向及拉伸行为,实现了充模过程中分子构型的可视化.最后,对带有两个圆形嵌件的环状型腔内的充模过程进行数值模拟研究,得到了充模过程中型腔内的温度、应力及分子构型的变化情况,并重点分析了充模速度、熔体温度和模具温度等对分子构型的影响.数值结果表明:本文提出的耦合模型可以成功模拟复杂型腔内充模过程中的温度、应力和分子取向等物理量的动态变化;适当提高注射速度可以增大熔接痕的强度;提升熔体温度和模具温度,可以有效改善甚至消除熔接痕.     

注塑模充模过程动态分析的有限元/控制体积法

塑料熔体在三维薄壁型腔内的流动属于带有运动边界的粘性不可压流体的流动,基于粘性广义Hele-Shaw流动理论,利用控制体积法建立了求解压力场的有限元方程,耦合利用有限元控制体积法自动跟踪熔体的运动边界,实现了充模过程的动态模拟.     

低温冷冻靶微管燃料充注流动冷凝特性

低温冷冻靶是惯性约束核聚变装置的关键部件之一。冷冻靶靶丸位于黑腔内部,需要利用微管完成燃料充注。针对燃料在微管内发生气液相变后流入靶丸这一关键过程,研究了微通道流体与常规流体流动传热的差异。基于微通道特性对流体流动传热方程进行修正,同时建立了气液相变模型,对微充气管内燃料充注过程进行了数值计算分析。得到重力和表面张力的影响,在微通道中,重力作用可以忽略,表面张力起重要作用。得到了微尺度效应包括速度滑移和温度跳跃对流动传热过程的影响。对多种充气管结构进行比较分析,为选型提供指导。通过选择不同进口条件和出口条件,对充注量控制和充注条件选择提供了指导方案,实际充注时需要同时提高进口温度和压力,保证连续可控充注。     

充模过程中熔接痕的改进光滑粒子动力学方法模拟与预测

提出了一种黏弹性流体的改进光滑粒子动力学(SPH)方法以试探性地模拟和预测黏弹性FENE-P熔体充模过程中熔接痕的形态演化. 首先基于SPH方法建立了聚合物流动的宏微观耦合模型, 同时提出了黏弹性流体的改进SPH离散格式. 随后, 通过模拟一些基准算例验证了改进的SPH方法模拟聚合物宏微观耦合问题的有效性及收敛性, 以及所提出的黏弹性温度模型的有效性. 最后, 模拟了环型腔内的充模过程, 试探性地展示了充模过程中微观分子的变形过程. 同时采用顺序热流道技术模拟了多浇口C形腔内的充模过程, 并与其他数值结果做比较. 数值结果表明: 对于大制件多浇口充模过程, 顺序热流道技术能够改善甚至消除充模过程中的熔接痕.     

急冷条件下Cu-Pb偏晶合金的相分离研究

研究了Cu-Pb偏晶合金的急冷快速凝固和组织形成规律，并通过将金属熔体的热传导方程和Navier-Stokes方程相耦合, 理论分析了合金熔体的冷却速率、液固相变时间等物理参量与液相分离之间的相关性. 研究结果表明，在急冷快速凝固条件下，熔体的快速冷却对偏晶合金组织形成的影响要比熔体内部液相流动的影响更为显著. 快速凝固使液相分离受到抑制，Cu-Pb偏晶合金均可获得均匀的微观组织结构. 随着冷速的增大，晶粒尺寸明显减小，凝固组织显著细化，晶体形态由粗大枝晶向均匀细小的等轴晶过渡. 提高冷却速率，缩短液固相变时间是重力场中抑制液相分离、获得均匀偏晶组织结构的重要条件. 关键词： 偏晶合金 快速凝固 液相分离 微观结构     

充液率对小型重力热管传热特性影响实验研究

开展了矩形槽道铝-乙醇小型重力热管的传热特性的实验研究,分析讨论了充液率对壁面温度分布、气液两相分布、热阻等热管传热性能的影响。研究表明,充液率对高热负荷工况的两相流动状态和传热特性有显著影响。两相脉动是高热负荷工况小型重力热管特有的两相流动现象.低充液率时,液塞易被气流冲破形成环状流,壁面温度几乎无波动。中等充液率时,在蒸汽和液塞的交替冲刷作用下,热管各段壁面温度均表现出脉动特性。高充液率时,液塞脉动速度的减小削弱了液塞对壁面的冲刷作用,壁面温度未出现明显波动.并且,中等充液率工况下气液两相的快速脉动增强了热管的传热性能,使得均温性和传热极限均优于低充液率和高充液率的情况.     

微通道内光热相变驱动流体运动特性研究

本文通过可视化实验,对微通道内光热效应致相变驱动流体运动特性进行了研究,通过红外聚焦激光跟随微通道内液柱气液界面进行加热持续产生的蒸发冷凝-聚合过程对液柱进行连续驱动。实验研究了激光功率、加热点位置对相变过程中的界面行为、冷凝液滴分布、驱动速率的影响规律。结果表明,激光功率越高,光斑距离界面越近,液柱蒸发速率越大,蒸汽浓度高,冷凝液滴分布越密集,驱动流体流动速度越快。     

黏弹介质包裹的充液管道中低频轴对称波传播特性

研究埋地充液管道中低频轴对称波传播特性。将土壤考虑为黏弹介质,结合Kennard薄壳方程和Kelvin-Voigt线性黏弹性模型,引入土壤载荷矩阵,推导出土-管滑移情形下流体主导波和管壁压缩波的相速度表达式。通过数值模拟计算得到流体主导波和管壁压缩波的频散和衰减曲线并进行可靠性验证,分析两种波引起的管壁径向位移之比,讨论厚径比和品质因子对流体主导波传播的影响。结果表明,黏弹介质对流体主导波和管壁压缩波的相速度影响较小,但对衰减影响较大;流体主导波对管壁径向位移有较大的影响,是泄露噪声传播的主要载体;厚径比越大,流体主导波的相速度越大,衰减越小;而品质因子越大,流体主导波的频散和衰减都越小。研究结果可为埋地充液管道的泄漏检测提供一定的理论参考。      

聚合物充填过程中裹气现象的数值模拟

针对聚合物充填过程中的裹气现象,采用一种有限元（FEM）-间断有限元（DG）耦合算法对其进行数值模拟。对于自由运动界面,采用水平集（Level Set）方法进行捕捉;用XPP（eXtended Pom-Pom）本构模型来描述黏弹性流体的流变行为。采用有限元-间断有限元耦合算法求解统一的流场方程,并采用隐式间断有限元求解XPP本构方程、Level Set及其重新初始化方程。数值结果与文献中的实验结果及模拟结果吻合较好,验证了数值方法的稳定性及准确性。分析带有非规则嵌件型腔内,注射速度及浇口尺寸对裹气现象的影响,裹气容易出现在较高注射速度及较小浇口的情形。     

An enthalpy-temperature hybrid method for solving phase-change problems and its application to polymer pyrolysis and ignition

A numerical approach based on the enthalpy method is proposed for solving generalized phase-change problems. The method is applied to predict pyrolysis and ignition of polymeric combustible materials. In contrast to the traditional approach, here both enthalpy and temperature are treated as independent variables, and the conservation equations are solved simultaneously in conjunction with the constitutive equations. Also, the formulation of the constitutive equations for the phase change is not necessarily the same for all of the possible phases, but can be chosen independently according to the characteristics of the physical problem and the requirements of the numerical analysis of each respective phase. Thus with this new approach, which we refer to as the enthalpy-temperature hybrid method, the enthalpy method is applicable to the generalized phase-change problems regardless of the form of the constitutive equations. The proposed method is first applied to a one-dimensional classical freezing problem for verification. It is found that the numerical results for the temperature history and the position of the phase-change interface agree well with the analytic solution existing in the literature. The method is then applied to the numerical simulation of the pyrolysis and ignition of a composite material with a polymer as the matrix and fibreglass as the filling material. Three models of oxygen distribution in the molten layer are considered to explore the melting and oxygen effects on the polymer pyrolysis. Numerical calculation shows that high oxygen concentrations in the molten layer enhance the pyrolysis reaction, resulting in a larger amount of pyrolysate, but in lower surface temperatures of the sample. It also shows that the distribution of oxygen in the molten layer has a strong effect on the pyrolysate rate, and therefore on ignition and combustion of the polymers. Comparison with available experimental data indicates that a model of oxygen distribution in the molten layer that is limited to a thin layer near the surface best describes the ignition process for a homogeneously blended polypropylene/fibreglass composite.     

复杂型腔充模中纤维取向的动态模拟

基于描述短纤维增强复合材料充模过程的气-固-液三相模型及同位网格有限体积法, 实现了纤维增强复合材料沿复杂型腔水平中面充模过程的动态模拟. 不仅得到了界面位置、各物理量的信息, 而且得到了纤维在型腔中的运动情况(包括纤维的平动和取向). 结果表明, 与沿型腔厚度方向纤维取向的表层-芯层结构不同, 纤维沿型腔水平中面的取向与型腔结构有关, 入口处纤维取向环绕型腔入口, 沿水平或竖直方向纤维取向与来流方向垂直, 型腔拐角处纤维取向指向拐点.     

Heat transfer in MHD flow of dusty viscoelastic (Walters’ liquid model-B) stratified fluid in porous medium under variable viscosity

The paper investigates the effects of heat transfer in MHD flow of viscoelastic stratified fluid in porous medium on a parallel plate channel inclined at an angle θ. A laminar convection flow for incompressible conducting fluid is considered. It is assumed that the plates are kept at different temperatures which decay with time. The partial differential equations governing the flow are solved by perturbation technique. Expressions for the velocity of fluid and particle phases, temperature field, Nusselt number, skin friction and flow flux are obtained within the channel. The effects of various parameters like stratification factor, magnetic field parameter, Prandtl number on temperature field, heat transfer, skin friction, flow flux, velocity for both the fluid and particle phases are displayed through graphs and discussed numerically.     

Elastic properties of a magnetic fluid with an air cavity retained by levitation forces

The paper describes the process of an air cavity rising in a magnetic fluid filling a tube with a bottom, transport, and retention of the cavity by magnetic levitation forces. The elastic and dissipative properties of a vibratory system with an inertial element that is a column of a magnetic fluid over an air cavity are considered. The possibility of using a transported air cavity as a movable reflector for a sound wave is evaluated.     

Aluminum integral foams with tailored density profile by adapted blowing agents

The goal of the present work is the variation of the structure of aluminum integral foams regarding the thickness of the integral solid skin as well as the density profile. A modified die casting process, namely integral foam molding, is used in which an aluminum melt and blowing agent particles (magnesium hydride MgH₂) are injected in a permanent steel mold. The high solidification rates at the cooled walls of the mold lead to the formation of a solid skin. In the inner region, hydrogen is released by thermal decomposition of MgH₂ particles. Thus, the pore formation takes place parallel to the continuing solidification of the melt. The thickness of the solid skin and the density profile of the core strongly depend on the interplay between solidification velocity and kinetics of hydrogen release. By varying the melt and blowing agent properties, the structure of integral foams can be systematically changed to meet the requirements of the desired field of application of the produced component.     

Continuity and momentum equations for moist atmospheres

The moist atmosphere with occurring precipitation is considered to be a multiphase fluid composed of dry air, water vapor and hydrometeors. These compositions move with different velocities: they take a macroscopic motion with a reference velocity and a relative motion with a velocity deviated from the reference velocity. The reference velocity can be chosen as the velocities of dry air, a gas mixture and the total air mixture. The budget equations of continuity and momentum are formulated in the three reference-velocity frames. It is shown that the resulting equations are dependent on the chosen reference velocity. The diffusive flux due to compositions moving with velocities deviated from the reference velocity and the internal sources due to the phase transitions of water substances result in additional source terms in continuity and momentum equations. A continuity equation of the total mass is conserved and free of diffusive flux divergence if the reference velocity is referred to the velocity of the total air mixture. However, continuity equations in the dry-air and gasmixture frames are not conserved due to the mass diffusive flux divergence. The diffusive flux introduces additional source terms in the momentum equation. In the dry-air frame, the diffusive flux of water substances and the phase transitions of water substances contribute to the change of the total momentum. The additional sources of total momentum in the frame of a gas mixture are associated with the diffusive flux of hydrometeors, the phase transitions of hydrometeors and the gasmixture diffusive flux. In the frame of total air mixture, the contribution to the total momentum comes from the diffusive flux of all atmospheric compositions instead of the phase transitions. The continuity and momentum equations derived here are more complicated than the traditional model equations. With increasing computing power, it becomes possible to simulate atmospheric processes with these sophisticated equations. It is helpful to the improvement of precipitation forecast.     

An immersed boundary energy-based method for incompressible viscoelasticity

Incompressible viscoelastic materials are prevalent in biological applications. In this paper we present a method for incompressible viscoelasticity in which the elasticity of the material is described in Lagrangian form (i.e. in material coordinates), and Eulerian (spatial) coordinates are used for the equations of motion and to enforce the incompressibility condition. The elastic forces are computed directly from an energy functional without the use of stress tensors, and the immersed boundary method is used to communicate between Lagrangian and Eulerian variables. The method is first applied to a warm-up problem, in which a viscoelastic incompressible material fills a two-dimensional periodic domain. For this problem, we study convergence of the velocity field, the deformation map, and the Eulerian force density. The numerical results indicate that the velocity field and deformation map converge strongly at second order and the Eulerian force density converges weakly at second order. Incompressibility is well maintained, as indicated by area conservation in this 2D problem. Finally, the method is applied to a three-dimensional fluid–structure interaction problem with two different materials: an isotropic neo-Hookean model and an anisotropic fiber-reinforced model.