固-液相变糊状区的格子Boltzmann研究

为深入探讨固-液相变材料在相界面内(糊状区)的流动和传热对固-液相变过程的影响,本文基于"焓-多孔介质"模型,提出研究糊状区的"多相流-多孔介质"两区域复合模型,将糊状区的高含液率视为多相流区域,而低含液率区域视为多孔介质,采用格子Boltzmann方法对固-液相变糊状区的流动和传热过程进行研究.重点讨论了相变温度半径,不同高、低含液率区的分界点γ_(tr)对糊状区的发展、迁移以及对相变过程的影响.研究结果表明糊状区厚度随相变半径的增大而增加,进而影响相变过程中流动和换热;对比复合模型的不同高、低含液率分界点γ_(tr)时和单纯多孔介质模型时的糊状区的边界、流动和换热情况,表明将糊状区高液相率区采用多相流模型进行分析,可以正确反映出高含液率区域内固相微粒随液相流动的特性及其对相变过程的影响。     

半透明介质相变中的辐射变物性问题

半透明介质相变传热在许多工程领域引起关注,其实质是变物性条件下伴随相变的辐射-导热瞬态耦合传热问题,其难点在于相变过程中的辐射变物性问题,尤其是折射率的时空演变问题.本文首次研究相变过程中的变折射率问题.假设相变介质存在固相-糊状-液相区,固/糊界面、液/糊界面假定为半透明漫反射.采用多层辐射传递模型模拟相变介质内的辐...     

基于伪焓法的多孔介质固液相变格子Boltzmann方法

格子Boltzmann方法是求解相变材料(PCM,phase change material)固液相变界面上复杂质量、能量和动量传递的有效方法。本文将伪焓LB模型推广至多孔介质固液相变中,构建相应的LB模型,并通过焓和温度的关系,获得下一时间步的焓、温度和对应液相率。为验证多孔介质固液相变LB模型的正确性,本文分别对一维Stefan问题和二维方腔多孔介质内固液相变问题进行求解。结果表明,在无自然对流作用下(一维Stefan问题),本文所构建模型能精确还原温度分布,与解析解相对误差控制在2%以内。除此之外,本文模型能反映自然对流对固液相变过程的影响,与已有工作差别较小,能准确追踪固液相变界面位置。     

微胶囊化相变悬浮液层流传热强化的参数分析

通过对微胶囊相变悬浮液管内层流恒热流对流换热的数学建模和模拟计算,获得在恒热流加热边界条件下,不同参数匹配时相变糊状区的确切范围以及固液相变两条相界面的数值结果。利用该模型对影响悬浮型固液两相流传热特性的诸多参数进行了比较细致的定量分析,这些参数包括固相体积浓度、斯蒂芬(Ste)数、粒径比、无量纲相变温度区间宽度以及无量纲过冷度等。分析和模拟计算的结果显示微胶囊浓度和斯蒂芬数对管内层流传热具有最重要的影响。     

固液相变过程最优控制数学定解分析及其算法

本文应用分布参数系统控制理论,采用变分法系统地分析,导出了二维固液相变传热过程最优控制问题的数学 定解条件,提出模型离散和数值求解方法,对共轭梯度最优控制算法进行适应性分析。     

相变储热体释热强度研究

对相变储热体以导热为主的释热过程,可以应用一维相变导热模型描述。当进一步简化,认为散热表面的温度低于相变温度且保持不变,不考虑液相区的自然对流和固液两相有明确的分界面,利用Neuman解法求解出释热过程的热强度,并通过实验进行研究。理论与实验结果比较表明,其相对误差不大于5％,从而证明这种方法可作为这类储热体的设计基础。     

有限厚度冻冰相变过程的传热研究

1前言确定相界面的瞬态位置是求解冻冰相变问题的重要部分。目前对于有限区域相变问题还没有可靠的精确解[1]。冻冰过程最简单分析方法是忽略固液两相的显热变化，建立相界面移动的准稳态传热关系式，近似预测固一液相变位置[2-4]。本文将考虑液、固相中显热的变化，对有限厚平板冻冰过程，采用精确解与积分解相结合的方法，分析冰层冻结过程中的固一液相界面、边界热流密度和贮入固、液相内的冷量随冻冰时间的变化规律。2物理一数学模型分析该研究的问题是一块贮冷板的两侧面直接与一个低温蒸发器直接接触，使贮冷板内冻冰贮冷。通常蒸发器…     

分析半透明单晶生长传热过程的边界元方法

论文发展了一个能求解带相变运动界面非定常传热和非线性热物理特性问题的双倒易边界元方法。数值模拟了半透明单晶生长中热过程的一个例子。由于方法是纯边界积分方法,计算量与计算内存都大大减少。获得了单晶生长过程瞬态温度场分布和固液相界面形状时间推进的一些结果。     

飞机结冰过程的传热研究

基于液/固相变的基本理论,对飞机结冰过程中冰层和液膜的传热特性进行了研究,建立了考虑固/液相区传热特性及固/液、气/液界面同时移动条件下的传热模型,并采用该模型对飞机结冰过程中的冰层和水膜生长特征进行了分析.研究表明,所建模型考虑了固、液相区的传热过程及液/固相变的非线性特征,对冰层生长速率及冰型特征有着较好的预测作用.     

合金材料激光表面重熔快速凝固过程数值模拟

根据合金固液盯变多区域性特点，用数值方法求解了适用于多区域相变的单相统一控制方程组，并据此对１Ｃｒ１８ＮｉｏＴｉ不锈钢激光表面重熔熔池的快速凝固过程进行了理论分析。     

Phase field simulation of the columnar dendritic growth and microsegregation in a binary alloy

This paper applies a phase field model for polycrystalline solidification in binary alloys to simulate the formation and growth of the columnar dendritic array under the isothermal and constant cooling conditions. The solidification process and microsegregation in the mushy zone are analysed in detail. It is shown that under the isothermal condition solidification will stop after the formation of the mushy zone, but dendritic coarsening will progress continuously, which results in the decrease of the total interface area. Under the constant cooling condition the mushy zone will solidify and coarsen simultaneously. For the constant cooling solidification, microsegregation predicted by a modified Brody- Flemings model is compared with the simulation results. It is found that the Fourier number which characterizes microsegregation is different for regions with different microstructures. Dendritic coarsening and the larger area of interface should account for the enhanced Fourier number in the region with well developed second dendritic arms.     

磁场对二元合金凝固过程中糊状层稳定性的影响

利用线性稳定性方法研究了外加磁场对二元合金凝固过程中糊状层稳定性的影响,且模型同时考虑了温度场、浓度场和流动的耦合作用.利用计算得出的色散关系式分析了磁场对糊状层稳定性的影响,其中包括直接模式和振荡模式.给出了不同情况下外加磁场对糊状层稳定性的影响,发现磁洛伦兹力可以减小由浮力引起的失稳效应.振荡模式下外加磁场对糊状层产生稳定作用,但直接模式下外加磁场对糊状层的稳定作用具有不确定性.本文所给出结果为工业中利用外加磁场改善产品的质量提供了重要的理论参考.     

定向凝固宏微观热质传递及双扩散流作用

针对氯化铵水溶液定向凝固实验,从枝晶微尺度到宏观尺度耦合模拟了相变热质传递过程,并考察了双扩散流的影响.研究结果表明,NH4CI-H2O溶液的枝晶微界面浓度变化与非平衡态的谢尔方程趋于一致.两相区内液相率较高,易于流动的产生;热浮升力形成逆时针回流,溶质浮升力形成顺时针回流;热浮升力作用大于溶质浮升力.但双扩散流几乎不影响凝固沿着结晶室高度方向进行的一维定向特性.数值模拟中低估了两相区的渗透率,是造成两相区厚度模拟结果小于实验值的原因.     

The effect of ultrasonic processing on solidification microstructure and heat transfer in stainless steel melt

The heat transfer in the ultrasonic processing of stainless steel melt is studied in this thesis. The temperature field is simulated when the metal melt is treated with and without ultrasound. In order to avoid the erosion of high temperature melt, ultrasound was introduced from the bottom of melt. It is found that the temperature of melt apparently increases when processed with ultrasound, and the greater the ultrasonic power is, the higher the melt temperature will be; ultrasonic processing can reduce the temperature gradient, leading to more uniform temperature distribution in the melt. The solidification speed is obviously brought down due to the introduction of ultrasound during solidification, with the increasing of ultrasonic power, the melt temperature rises and the solidification speed decreases; as without ultrasound, the interface of solid and mushy zone is arc-shaped, so is the interface of liquid and mushy zone, with ultrasound, the interface of solid and mushy zone is still arc-shaped, but the interface of liquid and mushy zone is almost flat. The simulation results of temperature field are verified in experiment, which also indicates that the dendrite growth direction is in accord with thermal flux direction. The effect of ultrasonic treatment, which improves with the increase of treating power, is in a limited area due to the attenuation of ultrasound.     

Solute distribution in KNbO3 melt-solution and its effect on dendrite growth during rapid solidification

This paper reports that the rapid solidification of mixed Li 2 B 4 O 7 and KNbO 3 melted in a Pt loop heater has been performed experimentally by the method of quenching, and various morphologies of KNbO 3 crystals have been observed in different regions of the quenched melt-solution. Dendrites were formed in the central region where mass transfer is performed by diffusion, whereas polygonal crystals with smooth surface grew in the marginal region where convection dominates mass transport. Based on measurement of KNbO 3 concentration along crystal interface by electronic probe analysis, it finds the variety of crystal morphologies, which is the result of different solute distributions: in the central region the inhomogeneity of solute concentration is much sharper and morphological instability is easier to take place; nevertheless in the marginal region the concentration homogeneity has been greatly enhanced by convection which prevents the occurrence of morphological instability. Additional solute distribution in the melt along the primary dendrite trunk axis as well as that in mushy zones has also been determined. Results show that the solute concentration in the liquid increases linearly with distance from the trunk tip and more solutes were found to be concentrated in mushy zones. The closer the mushy zone is to trunk tip, the lower the solute concentration will be there.     

温度梯度区域熔化作用下熔池迁移的元胞自动机模拟

本文采用耦合凝固和熔化效应的二维元胞自动机(cellular automaton, CA)模型,对温度梯度区域熔化(temperature gradient zone melting, TGZM)效应引起的熔池在固液两相区中的迁移现象进行模拟研究.模拟分析了抽拉速度、熔池初始位置、温度梯度和合金成分等因素对TGZM动力学的影响,并将模拟结果与解析模型的预测结果进行比较验证.通过模拟发现,在温度梯度作用下,熔池总是向着高温方向迁移;当抽拉速度低于或高于临界抽拉速度时,熔池朝向移动的液相线或固相线迁移;对于给定的抽拉速度,位于糊状区内临界位置以上的熔池会迁移进入液相,而位于临界位置以下的熔池会逐步靠近固相线.此外,温度梯度越高,合金成分越低,熔池的迁移速度越快.     

Aggregation and chimney formation during the solidification of ammonium chloride

Experiments study large-scale pattern formation during the growth of ammonium chloride (NH4Cl) from solution in a thin (Hele-Shaw) geometry. In particular a solid-liquid mixture ("mushy layer") forms in which growing solid NH4Cl crystals form a solid network interspersed with liquid. There are different ways that the mushy layer can be formed, however. If the cell is heated from below and cooled from above, thermal convection generates large-scale recirculating flows that carry seed crystals from the upper (cold) boundary to the (warmer) side and bottom boundaries. Ballistic deposition of these seed crystals leads to aggregation patterns with significant voids (filled with liquid) with a wide range of length scales. If the cell is cooled from below with a warm environment, the solid NH4Cl grows dendritically without deposition, resulting in a compact mushy layer. Plume convection within this mushy layer produces one or two well-defined "chimneys." If the environment is cool (comparable to the liquidus temperature of the solution), the mushy layer forms by a combination of dendritic growth and ballistic deposition, resulting in a more permeable mushy layer and enhanced chimney formation. The effects of ballistic deposition are enhanced if the cell is tipped, in which case the voids reappear. Plume convection and chimney formation are dramatically enhanced in this case. Additional experiments are done in which fluid flows in the system are enhanced artificially to verify that enhancements in chimney formation are due primarily to the aggregation process, and not to the increases in fluid flows due to thermal and compositional convection.     

Influence of void ratio on phase change of thermal energy storage for heat pipe receiver

In this paper, influence of void ratio on phase change of thermal storage unit for heat pipe receiver under microgravity is numerically simulated. Accordingly, mathematical model is set up. Numerical method is offered. The liquid fraction distribution of thermal storage unit of heat pipe receiver is shown. Numerical results are compared with experimental ones in Japan. Numerical results show that void cavity prevents the process of phase change greatly. PCMmelts slowly during sunlight periods and freezes slowly during eclipse periods as the void ratio increases. The utility ratio of PCM during both sunlight periods and eclipse periods decreases obviously as the void ratio increases. The thermal resistance of void cavity is much higher than that of PCMcanister wall. Void cavity prevents the heat transfer between PCM zone and canister wall.     

Focal shift of apodized truncated hyperbolic-cosine-Gaussian beam

Focal shift of hyperbolic-cosine-Gaussian beams induced by pure phase apodizer was investigated theoretically. The pure phase apodizer consists of three concentric zones: center circle zone, inner annular zone and outer annular zone. And the phase variance of the inner annular zone is adjustable. Results show that intensity peak moves far from optical aperture and then shrinks sharply for certain radii of zones with increasing phase variance of the inner annular zone. Simultaneously, one new intensity peak occurs near optical aperture, moves far from the optical aperture, and then becomes intensity maximum peak, and repeats the evolution process of the former intensity peak. Tunable focal shift occurs with focal switch. Decreasing the phase variance can change the move direction of the intensity peaks. In addition, the maximum distance between the two intensity peaks can be altered by beam parameters of cosh parts, and the distance value increases and then decreases with increasing beam parameters of cosh parts for certain radii of zones. Tunable focal shift is also discussed to construct optical tweezers.