淬火过程中具有非线性表面系数考虑相变时温度场的有限元分析

温度场的计算对淬火过程中热应力和热应变的分析有较大影响，对淬火后试件的残余应力和微观结构分析也有较大影响·本文以４２ＣｒＭｏ钢圆柱体试件为研究对象，从相变等温动力学ＴＴＴ图，给出了连续冷却相变动力学ＣＣＴ图的数学模拟计算式，计算了钢淬火过程中相变组织成分百分比，将热物性系数处理为温度和相变体积百分比的函数，用有限单元法计算了淬火过程中具有非线性表面系数考虑相变时的温度场·并建立了相应的泛函，为以后计算淬火的热应力和热应变做好准备工作     

膛口流场动力学过程数值研究

采用基于ALE方法的动网格及嵌入网格技术,运用有限体积方法,结合二阶精度Roe格式,对弹丸由高压气体驱动从静止状态加速至超音速,射出膛口到完全飞离初始流场的整个过程进行了数值模拟．根据数值结果,详细讨论了初始流场、火药燃气流场的形成与发展以及与弹丸的耦合和相互作用过程,揭示了在这一变化过程中激波与激波、激波与漩涡、激波与弹丸等的相互作用以及激波衍射、聚焦等对弹丸加速的影响．     

环形蓄热片相变材料固化时热传导问题的Green函数瞬时解

当环形蓄热片为相变材料所固化时,应用Green函数法分析其瞬时的温度分布.在一个柱形壳体中存储该相变材料(PCM),蓄热片上固体相变材料的厚度随时间而变化,并由Megerlin法求得.找到模型用Bessel方程构成的分析解,讨论了3种不同类型的边界条件,给出了解析解和数值解的比较.结果表明,所有情况下得到的蓄热片温度分布精度是令人满意的.     

金属材料热处理过程中的瞬态温度场与解态相变的数值计算方法

本文首先将Kirchhoff变换推广到导热系数为温度的多项式的非定常非线性热传导问题.并用分析方法确定热传导问题的边界条件.其次提出以孕育期叠加法并引用线性混合法则来模拟金属热处理过程的多相瞬态相变,较为简便地确定相变的开始时间、相变的种类及相变组织的数量.最后利用三维双重边界元法分析工件多种形式的热处理全过程.算例的数值计算结果表明本文方法行之有效.     

含湿相变粗糙多孔材质的热质耦合分形研究

多孔材质复杂的内部结构和含湿状态对传热和传质特性有着重要意义,其热质耦合传递过程广泛存在于能源开发和工程隔热等领域。不同于在多孔材质理想状态下对传热和传质特性的单方面分析,该文将孔道的分布参数、粗糙表面、含湿状态和相变等因素考虑进去,运用分形理论推导出了含湿相变粗糙表面多孔材质的渗流系数和耦合等效导热系数的表达式。结果表明,渗流系数与面积分形维数、含湿饱和度呈正相关,与相对粗糙度、迂曲分形维数呈负相关;耦合等效导热系数与渗流系数、相变量呈正相关,与相对粗糙度呈负相关。此外,结果还表明,相变量以及相变引起的气体膨胀压强差对热质耦合传递也有着重要影响。     

激光熔凝加工中瞬时温度场数值模拟

以激光熔凝表面强韧化处理为背景,应用空间弹塑性有限单元和高精度数值算法同时考虑材料组织性能的变化模拟工件的温度场,主要研究激光熔凝加工中瞬时温度场数值模拟,同时考虑相变潜热的影响,为第二步热应力场及残余应力的数值模拟做准备．用算例验证了模型的正确性,最后给出了激光熔凝加工不同时刻温度场分布．     

单向拉伸镍钛合金带从奥氏体到马氏体的相变分析

单向拉伸镍钛合金带中从奥氏体到马氏体的相变已在实验中观测到,并被看作为局部变形进行了数值模拟．该文采用相变理论对其进行分析,考虑了两相界面处变形梯度的跳跃以及Maxwell关系,导出了相变的控制方程．相变分析归结为寻求载荷的最小值,使在该值下控制方程具有唯一的、物理上可以接受的实数解．控制方程被数值求解,证明了该唯一解确实存在．相变的Maxwell 应力,马氏体相与奥氏体相内的应力与应变,以及相边界的倾角都可求出,并与实验所观测到的结果相吻合．     

垃圾填埋气体渗流过程中压力分布的滑脱解

考虑气体滑脱效应条件下,建立了垃圾填埋场气体渗流的数学模型．采用摄动法及积分变换法对模型进行解析求解,定量研究填埋气体的压力分布特征,并结合监测数据对模型参数与模型可靠性进行灵敏度分析和验证．结果表明:滑脱效应对气体渗流有较大影响,考虑滑脱效应条件的气体压力小于未考虑滑脱效应条件的气体压力,且得到的滑脱解与实测数值吻合较好．因此,研究填埋气体渗流时不能忽略滑脱效应．这不仅对于填埋场释放气体控制系统的设计和管理提供理论依据,而且可为低渗透油气藏工程开发过程中试井数据的确定提供科学的技术支持．     

经验正交函数与遗传算法结合的副热带高压位势场非线性模型反演

针对副热带高压的动力预报模型难以准确构建的困难,基于T106数值预报产品500 hPa位势高度场序列,用经验正交函数（EOF）分解方法对位势场序列进行了时、空分解,引入了动力系统重构思想,以EOF分解的空间模态的时间系数序列作为动力模型变量,用遗传算法全局搜索和并行计算优势,进行了动力模型参数的优化反演,建立了客观合理的非线性动力模型．通过对动力模型积分和EOF的时、空重构,实现了副热带高压的中、长期预报．试验结果表明,本文反演的动力模型的副热带高压预报效果优于常规的数值预报产品,该研究工作为副热带高压等复杂天气系统和要素场预报提供了新的方法思路和技术途径．     

VAPEX技术中溶剂气体传质规律现状与分析

VAPEX(蒸气萃取)技术是一种新型的稠油开采技术,VAPEX技术具有温室气体排放量低、开采原油时所需能量少和开发成本低等优点,对经济高效提高稠油采收率具有重要意义.综述VAPEX技术的发展历程,介绍溶剂气体注入到稠油样品中求解溶剂扩散系数的模型和溶剂分散系数模型,详细分析了溶剂气体进入到模型样品中体积、系统压力、溶解度、界面张力等变化对传质系数和扩散系数影响,溶剂气体在VAPEX过程中的分散系数,概述了VAPEX过程的连接路径和腔室变化情况和溶剂气体在稠油和沥青中的传质机理,并针对目前存在的问题,提出了VAPEX技术在稠油开采领域的前景方向及展望.     

Numerical simulation of immersion quenching process of an engine cylinder head

In this article, we present the numerical simulations of a real cylinder head quench cooling process employing a newly developed boiling phase change model using the commercial CFD code AVL-FIRE v8.5. Separate computational domains constructed for the solid and liquid regions are numerically coupled at the interface of the solid–liquid boundaries using the AVL-Code-Coupling-Interface (ACCI) feature. The boiling phase change process triggered by the dipping hot metal and the ensuing two-phase flow is handled using an Eulerian two-fluid method. Multitude of flow features such as vapor pocket generation, bubble clustering and their disposition, are captured very effectively during the computation, in addition to the variation of the temperature pattern within the solid region. A comparison of the registered temperature readings at different monitoring locations with the numerical results generates an overall very good agreement and indicates the presence of intense non-uniformity in the temperature distribution within the solid. Overall, the predictive capability of the new boiling model is well demonstrated for real-time quenching applications.     

Solution of non-linear thermal transient problems by a new adaptive time-step method in quenching process

Quenching process is a thermo-elastic-plasticity problem with a high material non-linearity. The numerical oscillation is likely caused in the simulation of quenching process. In order to avoid the numerical oscillation and improve the calculation accuracy of temperature and phase-transformation fields in the quenching process, a new self-adaptive time-step size method is presented. The method can adjust the time-step size according to the maximum and minimum differences of temperature fields between the previous simulation step and the current simulation step. FEM software for evaluating the temperature, stress/strain and phase-transformation is also developed. A cooling example with numerical analytical results and a quenching example with experiment results are used to verify the calculation accuracy of this software. Five methods including the method in this paper, two constant time-step sizes and two geometric proportion time-step sizes are applied to simulate the quenching process of a 40Cr steel cylinder, respectively. A comparison of the simulation results shows that, the method presented in this paper can effectively avoid the numerical oscillation, ensure the calculation accuracy and cost less calculation time.     

A direct method for determining the mass of gases in ICF shell

In the present paper, the defects of dew point method for measuring the mass of gas filled in ICF shells are analyzed. An accurate state equation for gas D2 is deduced from Benedict-Webb-Rubin (BWR) equation and experimental data in planar phase. A direct method to determine gas mass in ICF shells via measuring the temperature and pressure outside the shells and solving the equation of state by numerical method is proposed. It overcomes the theoretical defects of dew point method and the complexities of equipment. In the present method, the state equation can be improved by more accurately measuring P-V-T values of gas D2, so the measuring precision of the mass of gas in the shells can also be improved. The present method is effective for treating mix gases filled in the shells as well. The errors between the computational results and experimental data are very small. Some cases in the filling process are predicted, and the proper temperature and pressure for filling gases effectively are also suggested.     

Two-phase laminar axisymmetric jet flow: Explicit,implicit, and split-operator approximations

An hybrid Eulerian-Lagrangian numerical scheme is developed for a two-phase problem and four finite-difference schemes are compared. For this purpose, the problem of hydrodynamic and thermal interactions between a fuel spray and a mixing region of two laminar, unconfined axisymmetric jets is formulated in terms of a set of parabolic differential equations for the gas phase and a set of Lagrangian ordinary differential equations for the condensed phase. Consistent, second-order accurate hybrid numerical schemes, with the exception of the explicit scheme with an accuracy between linear and quadratic, are used to solve these equations. The subset of gas-phase equations has been solved by four different numerical methods: a predictor-corrector explicit method, a sequential implicit method, a block implicit method, and a symmetric operator-splitting method. The subsystem of liquid-phase equations is solved along the droplet trajectories by a second-order Runge-Kutta scheme. The computations have been made to predict the hydro-dynamic and thermal mixing regions of the gas phase as well as the trajectories of each individual group of droplets. In addition, the size, velocity and temperature associated with each group are predicted along these trajectories. The relative merits of the above four difference-schemes are discussed by constructing effectiveness curves. At low error tolerances, the sequential implicit method gives the best results, where for large error tolerances, the explicit and operator splitting give better results. The block implicit scheme is the least effective at all accuracy requirements.     

The use of implicit flux limiting schemes in the simulation of the drying process: A new maximum flow sensor applied to phase mobilities

In this work, a previously proposed two-dimensional wood drying computational model is modified so that an analysis of the use of flux limiters for reducing numerical dispersion in the drying kinetics can be undertaken. In particular, a new sensor based on the ratio of phase fluxes is used to flux limit the liquid and gas phase mobility tensors. An extremely important physical phenomenon that arises during high temperature drying, and one that is exacerbated by the use of the flux limiter, relates to the medium becoming fully saturated. To overcome the numerical convergence problems associated with this phenomenon, the concept of a fixed compressible phase is introduced within the model and the averaged air density is used as one of the primary solution variables. A comparison of the simulation results for high temperature drying will be presented for three different case studies. In case one, a completely isotropic sample is considered while in case two, a slight anisotropy (1:2) is introduced for the modelling of a radial–transverse wood cross-section. In case three, for a wood sample modelled in the radial–longitudinal cross-section, a strongly anisotropic medium (1:1000) is analysed. It is found in all cases that when the new sensor is used during flux limiting, superior results to upstream weighting are obtained. In fact, it is possible with this method to attain the accuracy offered by upstream weighting on a much finer mesh at a considerable reduction of at least a factor of 5 in the overall computation time.     

Numerical simulation of gas-assisted injection molding using CLSVOF method

It is a typical gas-liquid two phase flow phenomenon that gas penetrates the polymer melt in gas-assisted injection molding (GAIM) process. Numerical simulation is now playing an important role in GAIM, in which the accurate simulation of moving interface is of great importance. The level set (LS) method is a popular interface tracking method, but it does not ensure naturally mass-conservation. In order to improve the mass-conservation of LS method, a coupled level-set and volume-of-fluid (CLSVOF) method with mass-correction is presented for the numerical simulations of interfacial flows in GAIM. The performance of this CLSVOF method is demonstrated by two numerical tests including the three-dimensional deformation field test and the dam break problems. Finally the CLSVOF method is employed to simulate the 3D moving interfaces in GAIM, including gas-melt interface and the melt-front interface. The influences of melt temperature and gas delay time are also analyzed detailedly. As a case study, the processes that gas penetrates the polymer melt in complex cavities are also simulated using this method, and the simulation results are in agreement with those obtained by other researchers.     

Numerical analysis of quenching – Heat conduction in metallic materials

Metallic materials present a complex behavior during heat treatment processes. In a certain temperature range, change of temperature induces a phase transformation of metallic structure, which alters physical properties of the material. Indeed, measurements of specific heat and conductivity show strong temperature-dependence during processes such as quenching of steel. Several mathematical models, as solid mixtures and thermal–mechanical coupling, for problems of heat conduction in metallic materials, have been proposed. In this work, we take a simpler approach without thermal–mechanical coupling of deformation, by considering the nonlinear temperature-dependence of thermal parameters as the sole effect due to those complex behaviors. The above discussion of phase transformation of metallic materials serves only as a motivation for the strong temperature-dependence as material properties. In general, thermal properties of materials do depend on the temperature, and the present formulation of heat conduction problem may be served as a mathematical model when the temperature-dependence of material parameters becomes important. For this mathematical model we present the error estimate using the finite element method for the continuous-time case.     

Structure and crystallization of bulk amorphous Pd41Ni10Cu28P21 alloy

Bulk amorphous Pd41Ni10Cu28P21 alloy has been prepared by water quenching method. The system shows excellent glass forming ability (GFA) with a high value of reduced glass transition temperature Trg, 0.714. Structural analyses indicate that the Pd41Ni10Cu28P21 alloy has a dense packing structure closer to "frozen liquid" than that of amorphous Pd40Ni40P20 alloy. Experiments on crystallization reveal that several crystalline phases simultaneously precipitate in the early part of crystallization. Below 710 K, a metastable phase forms, and subsequently disappears at elevated temperatures. In addition, the influence of partial substitute of Cu for Ni on GFA has been discussed with regard to thermodynamics and kinetics.