外加磁场对熔融半导体CZ法拉晶过程中流动场影响的数值计算

数值分析结果表明,外加磁场可以改变熔融半导体中的流型,几千高斯的磁场可以显著地减小熔体的流动,但对温度场影响不大。     

条形药包爆炸作用在土中显著影响范围的数值模拟

爆炸固结法是软土地基处理的新方法。为了分析爆炸作用在土中的变化规律,通过有限元方法对单孔条形药包在土中爆炸的数值模拟,在验证计算结果与实测数据的基础上,对地基中单孔条形药包在上覆压力和药量两组爆炸参数下进行模拟对比分析,以解决在不同上覆压力和药量条件下条形药包爆炸作用在土中的显著影响范围。针对条形药包土中爆炸问题,运用ALE算法解决了炸药网格的畸变问题。研究结果表明,上覆压力发挥作用的主要时间不是在爆炸瞬间,而是在爆炸以后的时间段,主要的贡献是对爆炸产生的孔隙水压的消散作用。而炸药的药量对于布置爆炸孔和竖向排水通道的间距影响具有较大影响。     

流向磁场作用下圆柱绕流的直接数值模拟

绕流是托卡马克装置中液态包层内常见的流动形态,对流场与热量分布有着重要的影响.本文通过直接数值模拟(DNS),研究了不同磁场强度下Re=3900的圆柱绕流,分析了磁场强度对于湍流尾迹的影响.无磁场情况下,直接数值模拟的结果与前人的实验及模拟结果吻合很好.圆柱下游的尾迹中,随着流向距离的增大,流向速度剖面逐渐从U型进化呈V型,并慢慢趋于平缓,这表明尾迹中的流动结构受圆柱影响逐渐减小.圆柱后方两侧的剪切层中,由于Kelvin-Helmholtz不稳定性的影响,可以清晰地看到小尺度剪切层涡的脱落.通过对无磁场的计算结果施加流向磁场,本文计算了哈特曼数(Ha)分别为20, 40和80的工况,以研究磁场效应对于湍流的影响.结果表明磁场较弱时,流动依然呈三维湍流状态.随着磁场增强,近圆柱尾流区受磁场抑制明显,回流区被拉长,剪切层失稳位置向下游转移.圆柱后方的涡结构由于受到竖直方向洛伦兹力的挤压作用,随着哈特曼数的增加尾迹区域逐渐变窄.相比于无磁场情况的涡结构,由于磁场的耗散作用,相应的涡结构尺度变小.该研究不仅扩展了现有磁场下湍流运动的参数范围,对于液态包层的设计及安全运行同样具有重要的理论指导意义和工程应用价值.     

铸造充型流动过程的三维数值模拟

本文以ＳＯＬＡ－ＶＯＦ方法为基础,结合作者提出的三维自由表面边界速度确定方法,实现了铸造充型流动过程的三维数值模拟,模拟结果准确地提供了液体金属在充型过程中的流动形态、充型时间、流量等实用的参数．这表明,通过数值计算的手段,对用实验方法难以观测的铸造充型流动过程进行分析研究是有意义可行的,所得到的参数有助于铸造工艺的设计及优化     

磁场对双扩散液层热毛细对流的影响

通过数值模拟的方法对磁场作用下的双扩散液层热毛细对流进行了研究, 模型中同时考虑了热毛细效应和溶质毛细效应的存在. 研究结果显示, 外部磁场能够有效削弱液层内热毛细对流的强度, 改变热毛细对流的对流结构; 随着磁场强度的增大, 液层内热毛细对流的对流强度逐渐减小, 热质传递过程中扩散效应逐渐得到增强; 最终, 溶质浓度沿水平方向呈梯度分布. 因此, 当磁场强度足够大时能够实现晶体生长中所需的纯扩散条件.     

微重力环境和外加磁场对晶体生长中杂质分凝的影响

浮区结晶法是生长体单晶的重要方法之一．本文从理论和实验两方面综述了浮区结晶法中熔体内宏观场对体单晶中杂质分布均匀性影响的定性和定量研究,并从改善体单晶杂质分布均匀性的角度出发,总结了目前所采取的两种重要方法：微重力环境和外加磁场的研究现状     

线圈绕Y轴倾斜时方腔内空气磁力与重力耦合对流的数值模拟

数值分析了微重力与重力环境下圆形载流线圈绕Y轴倾斜时方腔内空气热磁对流。磁场计算采用毕奥萨伐定律;控制方程基本变量采用控制容积法离散,求解采用SIMPLE算法。计算过程中Ra数的变化范围为104~105,线圈倾斜角yeuler的变化范围为-90°~90°,磁场力数γ数的变化范围为0~200。获得了空气热磁对流的流场和温...     

微重力环境下Marangoni对流的有限元数值模拟

本文以悬浮区为背景研究液桥中气液交界面上由表面张力所驱动的流体对流。我们采用有限元方法对轴对称、定常运动方程,能量方程,自由面切向、法向应力平衡条件迭代求解,首次考虑了边界形状的影响,得到了自洽的流区位形和流场、温度场、表面压强分布。一般说来,流区自由面呈弯月形。结果表明,只要表面张力数 C_a<1,静态平衡的流区位形就是动态情况的良好近似。本文还分析了C_a数及G_r 数对流区位形的影响,得到了不同 M_a 数及散热条件下的温度场和流函数分布。     

高超声速非定常流动的数值模拟与气动热计算

高超声速飞行器研究中的一个重点问题是飞行器表面的气动加热,它对飞行器的气动、热特性及安全性有重要的影响.受到当前实验技术的限制,地面实验无法准确模拟真实飞行条件,所以采用数值模拟研究气动加热问题成为目前重要的研究手段.本文采用数值方法求解三维N-S方程,得到钝头体再入模型绕流的瞬态流场,驻点温度及表面热流沿轨道变化规律.计算中采用变边界条件模拟沿轨道飞行的非定常性.     

高超声速流动中液滴蒸发对流动参数的影响研究

采用带有非平衡相变模型的数值方法对高超声速流动中液滴蒸发的影响进行了研究,重点探讨了液滴蒸发对斜激波下游区域流场和劈面参数的影响.研究结果表明:液滴在激波后区域的蒸发存在弛豫现象,蒸发弛豫过程与激波强度相关,并影响模型表面气动参数的均匀性;当偏转角从3°增大至16°时,研究发现特别是偏转角为10°时,激波强度刚好能使液滴在整个激波下游区域完全蒸发,此时蒸发引起的流场参数不均性影响最大,与不考虑蒸发时相比,液滴质量含量、压强、静温减小幅度分别达到27.2g/kg、401Pa、54K,该三项参数变化量的比分别达到96.8%、4.8%、14.2%.     

Influence of external magnetic field on the flow field in molten semiconductor of Czochralski crystal growth —A numerical simulation

Numerical results show that an external magnetic field may influence significantly the flow pattern in the molten semiconductor of Czochralski crystal growth. The melt flow could be pronouncedly damped by a magnetic field with the intensity of several thousands Gauss, while the temperature field is affected only in a less extent by the magnetic field. The project is supported by the National Natural Foundation of China.     

Numerical study of buoyancy- and thermocapillary-driven flows in a cavity

Thermocapillary- and buoyancy-driven convection in open cavities with differentially heated endwalls is investigated by numerical solutions of the two-dimensional Navier-Stokes equations coupled with the energy equation. We studied the thermocapillary and buoyancy convection in the cavities, filled with low-Prandtl-number fluids, with two aspect-ratiosA=1 and 4, Grashof number up to 10⁵ and Reynolds number ⋎Re⋎≤10⁴. Our results show that thermocapillary can have a quite significant effect on the stability of a primarily buoyancy-driven flow, as well as on the flow structures and dynamic behavior for both additive effect (i.e., positiveRe) and opposing effect (i.e., negativeRe).     

Influences of Marangoni effect on solidification of curved interfaces with constant curvature

The asymptotic method and phase field simulation are applied to study the influences of variation of the surface tension with temperature on the movements of solid-liquid interfaces in the solidification process of spherical and cylindrical nuclei, respectively. Results indicate that Marangoni effect will increase the critical nucleation radius, slow up the movement of interface. The tip speed of dendrite decreases linearly with Marangoni number for melt without convection. The results of phase field simulation are qualitatively in accord with that of asymptotic method.     

用PSM模型模拟聚合物熔体轴对称收缩流动

采用模拟黏弹流体挤出胀大的方法,计算了IUPAC-LDPE熔体经过4：1轴对称收缩流道的流动．计算的相对涡强度、入口校正和献中的结果基本一致,给出的流场也显示出计算结果是合理的．表明该方法能够适用于用积分型PSM模型表征的黏弹流体在收缩流道内的流动模拟．     

环形浅池内硅熔体中的热流体波

为了了解工业Czochralski炉内硅熔体表面轮型的基本特征,对环形浅池内硅熔体的热毛细-浮力对流进行了三维数值模拟,硅熔池内径为15 mm,外径为50 mm,深度为3 mm,熔池外壁被加热,内壁被冷却,底部固壁和顶部自由表面均绝热或者允许一个垂直方向的传热。模拟结果表明,当径向温差较小时,熔池内会产生稳定的单胞热毛细-浮力流动,随着温差的增大,流动将转变为三维振荡流动,在熔体自由表面会出现沿周向运动的轮型,小的垂直方向的热流密度(3W/cm2)对这种振荡流动没有大的影响。同时,讨论了流动和温度波动的特征,并确定了振荡流动的临界条件。     

Objective function choice for control of a thermocapillary flow using an adjoint-based control strategy

The problem of suppressing flow oscillations in a thermocapillary flow is addressed using a gradient-based control strategy. The physical problem addressed is the “open boat” process of crystal growth, the flow in which is driven by thermocapillary and buoyancy effects. The problem is modeled by the two-dimensional unsteady incompressible Navier–Stokes and energy equations under the Boussinesq approximation. The goal of the control is to suppress flow oscillations which arise when the driving forces are such that the flow becomes unsteady. The control is a spatially and temporally varying temperature gradient boundary condition at the free surface. The control which minimizes the flow oscillations is found using a conjugate gradient method, where the gradient of the objective function with respect to the control variables is obtained from solving a set of adjoint equations. The issue of choosing an objective function that can be both optimized in a computationally efficient manner and optimization of which provides control that damps the flow oscillations is investigated. Almost complete suppression of the flow oscillations is obtained for certain choices of the objective function.     

Numerical simulation of the flow field in a thermal plasma reactor

A numerical simulation is presented for a thermal plasma reactor with particle-trajectory model in this paper.Turbulance is considered by using simple SGS model.Thegoverning equations are solved by means of the algorithm of SIMPLER.The calculatedresults give the velocity and the temperature fields within plasma reactor,and thetrajectories of the injected particles.     

OPTICAL DIAGNOSTIC AND MODELING SOLUTION GROWTH PROCESS OF SODIUM CHLORATE CRYSTALS

Both a real time optical interferometric experiment and a numerical simulation of two-dimension non-steady state model were employed to study the growth process of aqueous sodium chlorate crystals. The parameters such as solution concentration distribution, crystal dimensions, growth rate and velocity field were obtained by both experiment and numerical simulation. The influence of earth gravity during crystal growth process was analyzed. A reasonable theory model corresponding to the present experiment is advanced. The thickness of concentration boundary layer was investigated especially. The results from the experiment and numerical simulation match well.     

Numerical Simulation of Plasma Dynamics in a Nonuniform Magnetic Field

An efficient algorithm is proposed enabling numerical simulations of plasma dynamics in a nonuniform magnetic field. The present numerical data are in good agreement with experimental data obtained in a GOL-3 setup and with previous simulations. The experimentally observed effect of fast transfer of energy to ions is confirmed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 1, pp. 35–45, January–February, 2006.