带环状翅片管式相变储热器的数值模拟

相变材料的低导热性是相变储热器传热过程的主要障碍。在相变材料侧进行强化,是目前提升相变储热器蓄热速率的主要方法。本文采用有限容积法对带有环状翅片的管式相变储热器的蓄热过程进行了数值模拟,得到了温度场、相界面随时间的变化。在此基础上,本文对翅片导热系数、翅片厚度和翅片间距等影响储热速率的因素进行了计算和分析。为相变储热器的优化设计提供了一定的参考。     

基于焓法的固-液相变格子Boltzmann模型

潜热储能系统中的相变材料固液相变问题是一个重要的多相流动问题.本文研究了在具有不同类型源项下的焓法格子Boltzmann相变模型在还原温度场分布上的特性.为了消除外力项对温度场的额外影响,本文对源项进行了修正,并通过求解单区域相变问题和自然对流相变问题,验证了由Chapman-Enskog展开得到的焓守恒方程中源项的影响.结果 表明,在单区域相变问题中,改进模型的温度分布更加准确.另外,在有相变的自然对流问题中,没有潜热附加项的改进模型的熔化过程比其他模型略慢,同一时刻的最大液相率与最小液相率之差小于0.01.     

相变存储单元RESET多值存储过程的数值仿真研究

使用数值仿真方法对相变随机存储器存储单元的RESET操作的多值存储过程进行了研究,建立了三维存储单元模型,用有限元法解Laplace方程及热传导方程以模拟电脉冲作用下的存储单元物性变化过程.计算出单元内相变层的相态分布及单元整体电阻,分析了单元内部尺寸变化对多值存储过程及状态的影响.结果表明,通过精确控制输入电脉冲,相变存储单元能够实现4值存储;多值存储状态受单元内相变层厚度及下电极接触尺寸变化的影响较大;存储状态在80℃的环境温度下均可保持10年以上不失效.     

球墨铸铁QT600-3激光相变硬化数值模拟研究

采用商业有限元软件ProCAST,建立了激光相变硬化温度场的3维数值模拟模型,对球墨铸铁QT600-3进行激光相变硬化的数值模拟研究,计算了表面淬火的温度场,根据温度场预测了硬化层的深度和宽度。激光功率在800~1 000 W之间,激光扫描速度在2.000~2.667 mm/s之间,光斑直径在4~5 mm之间时,计算得到的硬化层深度在0.20~0.64 mm之间,硬化层宽度在2.0~3.7 mm之间。球墨铸铁数值模拟结果与试验基本吻合,两者变化趋势一致。     

基于有限单元法和计算流体动力学的光机热集成仿真分析研究

使用计算流体动力学、传热学、固体弹性力学和光学的多物理场耦合集成仿真分析技术完成了红外热像仪中光机的热特性与热设计研究.采用有限单元法(Finite Element Method,FEM)和计算流体动力学(Computational Fluid Dynamics,CFD)得到了光机内部的温度场、应力变形分布及光学件面形...     

利用迈克耳孙干涉仪重建蜡烛火焰温度场

利用教学实验中的迈克耳孙干涉仪,进行蜡烛火焰温度场的测量.对干涉条纹进行了处理和细化,并分别用阿贝尔变换法和环带法对蜡烛火焰的温度场进行了计算,得到蜡烛火焰横截面上温度场的定量结果.     

一维伊辛模型的相变-Ⅱ

找到了一个能够用矩阵法计算的新序参量,此参量既能给出有限数量格点时一维伊辛模型中存在相变,又能给出无限数量格点时相变消失的结果.利用此序参量求出了一个计算相变点的简洁近似公式.     

变折射率介质相变过程辐射导热耦合换热RKPM模拟

基于无网格再生核粒子法(RKMP),采用等效热容法处理相变潜热的吸收和释放,通过求解梯度折射率辐射传输方程处理变折射率问题,建立半透明材料相变过程的辐射导热耦合统一的计算方法。计算了一维晶体和非晶体半透明材料凝固过程温度响应特性,考虑了折射率变化对相变过程的影响,通过和文献值对比,验证了RKPM能够模拟含有相变过程的辐射传输问题,且具有良好的精度。     

有限元与无网格伽辽金耦合法分析二相连续多孔介质

将有限元法和无网格伽辽金法相结合,创造出一种新的数值计算方法,该法既能利用有限元法的优点,又能发挥无网格伽辽金法局部化技术的无单元特性,使其具有计算优势.推导了二相多孔介质的有限元法与无网格伽辽金法耦合的离散方法.给出了两个算例,在计算二维固结方程时,精度较高且处理液体流量准确,是有限元法的有力补充.     

高能激光破岩的传热学特性研究

高能激光破岩中激光与岩石相互作用的过程,是一个复杂的非稳定的传热学问题。利用能量守恒原理,得到了激光破岩的传热方程,并对其相变过程中固液相变界面、液气界面上的焓、导热系数和比热等热物性参数进行了近拟简化,从而建立了激光破岩温度场分析的物理模型。利用Galerkin方法,对激光破岩的温度场模型进行了数值分析求解和实例计算。     

一种改进的无单元伽辽金方法

使用单位分解积分,对传统的无单元伽辽金方法进行改进.有限覆盖和单位分解是单位分解积分的数学基础,对单位分解积分进行了严格证明,并指出使用Shepard函数作为单位分解函数是一个很好的选择.数值实例表明,使用单位分解积分进行数值求积的无单元伽辽金方法是一种真正的无网格方法,与经典的背景网格积分相比具有更高的精度.     

Free Mesh Method: fundamental conception, algorithms and accuracy study

The finite element method (FEM) has been commonly employed in a variety of fields as a computer simulation method to solve such problems as solid, fluid, electro-magnetic phenomena and so on. However, creation of a quality mesh for the problem domain is a prerequisite when using FEM, which becomes a major part of the cost of a simulation. It is natural that the concept of meshless method has evolved. The free mesh method (FMM) is among the typical meshless methods intended for particle-like finite element analysis of problems that are difficult to handle using global mesh generation, especially on parallel processors. FMM is an efficient node-based finite element method that employs a local mesh generation technique and a node-by-node algorithm for the finite element calculations. In this paper, FMM and its variation are reviewed focusing on their fundamental conception, algorithms and accuracy.     

Efficient Analysis of Millimeter Wave Ferrite Circulators by the GMRES Iterative Algorithm

The finite element method (FEM) combined with the perfectly matched layers (PML) is given for simulation of waveguide ferrite circulators. The generalized minimal residual (GMRES) iterative method is applied to solve such sparse large non-symmetric system of linear equations resulting from the use of edge-based finite element method. The formulation of FEM and the algorithm of GMRES method are described in detail. The reflection and insertion losses of millimeter wave waveguide circulator are analyzed and the results are compared with those obtained from literature.     

Finite element method for conserved phase fields: Stress-mediated diffusional phase transformation

Phase-field models with conserved phase-field variables result in a 4th order evolution partial differential equation (PDE). When coupled with the usual 2nd order thermo-mechanics equations, such problems require special treatment. In the past, the finite element method (FEM) has been successfully applied to non-conserved phase fields, governed by a 2nd order PDE. For higher order equations, the convergence of the standard Galerkin FEM requires that the interpolation functions belong to a higher continuity class.     

Analysis of Complex Dielectric Phase Gratings by Hybrid FEM-Gaussian Mode Expansion Technique

Hybrid technique of finite element method (FEM) and Gaussian mode expansion method is implemented for the analysis of transmission property of dielectric phase gratings used in holographic power combining systems. It takes advantages of the versatility of FEM and the simplicity of the mode expansion method. Apparently, this hybrid method is well suited for the analysis of phase gratings with complex profiles, which provide more degrees of structural freedom than those of classical rectangular groove gratings. The transmitted fields of an array of two and three fundamental Gaussian beams through the multi-stepped and curved gratings are calculated by the hybrid method. Numerical results show that much less phase distortion can be fulfilled for the reconstructed pseudo-plane wave by the modified phase gratings, especially the curved ones. This is advantageous for improving the combining efficiency in spatial millimeter wave and sub-millimeter wave power combining systems.     

微波元件与加速结构的有限元分析

研究了有限元方法在二维均匀结构和轴对称加速结构中的应用,采用了带宽优化技术和子空间求解特征值方法,并给出了部分例子,结果表明该方法精度高、速度快.这些方法可直接用于三维程序、并行计算.     

激光辐照受拉铝板破坏行为的时空多尺度数值模拟

推导耦合过渡区内参变量信息交换的元/网格动量传递多尺度算法,建立离散元与有限元耦合时空多尺度计算模型,并应用于激光辐照下受拉铝板破坏行为的数值模拟中.通过对比有限元计算模型、空间多尺度计算模型与时空多尺度计算模型在激光辐照下受拉铝板破坏算例的模拟结果,验证离散元与有限元耦合时空多尺度计算模型的准确性和数值计算高效率优势.使用该多尺度计算模型从宏观和细观尺度对铝板破坏行为进行数值模拟,模拟结果与实验结果基本一致.     

An element-free Galerkin(EFG) method for numerical solution of the coupled Schrdinger-KdV equations

The present paper deals with the numerical solution of the coupled Schrdinger-KdV equations using the elementfree Galerkin(EFG) method which is based on the moving least-square approximation.Instead of traditional mesh oriented methods such as the finite difference method(FDM) and the finite element method(FEM),this method needs only scattered nodes in the domain.For this scheme,a variational method is used to obtain discrete equations and the essential boundary conditions are enforced by the penalty method.In numerical experiments,the results are presented and compared with the findings of the finite element method,the radial basis functions method,and an analytical solution to confirm the good accuracy of the presented scheme.     

Comparison and implementation of the various numerical methods used for calculating transmission loss in silencer systems

Issues concerning the design and use of large-scale silencers are more prevalent today then ever before. With the increased use of large industrial machinery (such as gas turbines) and the increase in public awareness and concern for noise control, the desire to be able to properly design silencers for specific applications is increasing. Even today, most silencer design is performed by simply modifying existing designs without full confidence of the new performance characteristics. Due to the size and expense of these silencers, it would be beneficial to have means to predict the insertion loss (IL) or transmission loss (TL) characteristics at the design stage. To properly accomplish this, many factors such as geometry, absorptive material properties, flow effects, break out noise, and self-generated noise must be considered. The use of the finite element method (FEM) and the boundary element method (BEM) can aid in the prediction and design. This paper examines three of the different methods used in calculation of TL values; namely the “traditional” laboratory method, the 4-pole transfer matrix method and the 3-point method. A comparison of these methods based on such criteria as accuracy, computation time, and ease of use was conducted. In addition, the idiosyncrasies and problems encountered during implementation are presented. The conclusions were that the FEM is better suited for this kind of application and that the 3-point method was the fastest method and was easier to use than the 4-pole method.