自然对流情形下激光辐照液体贮箱的机理

利用激光辐照靶目标时,被辐照部位可能是液体贮箱。通过实验测量与数值模拟的紧密结合,揭示了液体处于自然对流状态时激光辐照下贮箱侧壁温升及液体速度场的演化规律。结果表明：激光辐照初期,铝板中心温升率较高;随着壁面附近液体温度的升高,光斑附近速度边界层内的最大流速增大,传热强度亦增大,导致铝板温升率降低;当铝板吸收的激光能量能够基本被水的对流带走时,铝板中心的温升率趋于零。     

On the Fully Implicit Solution of a Phase-Field Model for Binary Alloy Solidification in Three Dimensions

A fully implicit numerical method, based upon a combination of adaptively refined hierarchical meshes and geometric multigrid, is presented for the simulation of binary alloy solidification in three space dimensions. The computational techniques are presented for a particular mathematical model, based upon the phase-field approach, however, their applicability is of greater generality than for the specific phase-field model used here. In particular, an implicit second order time discretization is combined with the use of second order spatial differences to yield a large nonlinear system of algebraic equations as each time step. It is demonstrated that these equations may be solved reliably and efficiently through the use of a nonlinear multigrid scheme for locally refined grids. In effect, this paper presents an extension of earlier research in two space dimensions (J. Comput. Phys., 225 (2007), pp. 1271-1287) to fully three-dimensional problems. This extension is validated against earlier two-dimensional results and against some of the limited results available in three dimensions, obtained using an explicit scheme. The efficiency of the implicit approach and the multigrid solver are then demonstrated and some sample computational results for the simulation of the growth of dendrite structures are presented.     

复杂多孔介质腔体内自然对流换热的数值研究

采用曲线坐标系下压力与速度耦合的SIMPLEC算法,数值研究复杂多孔介质腔体内的自然对流换热问题.腔体的曲面温度分别保持恒定,上下表面绝热.在曲线坐标系中用有限容积法离散方程,并采用Brinkman扩展达西模型及局部非热平衡模型求解,综合研究Rayleigh数,Darcy数、孔隙率等参数对腔体内自然对流换热的影响.计算结果表明:Rayleigh数和Darcy数的影响最大而孔隙率的影响很小,同时存在使得腔体内换热达到最强的最佳纵横比.     

考虑Soret和Dufour效应的方腔内双扩散自然对流格子Boltzmann模拟

采用格子Boltzmann方法,考虑Soret和Dufour效应,对内置高浓度发热圆的方腔内部双扩散自然对流现象进行数值模拟.高浓度发热圆位于方腔中心,四周壁面均为低温低浓度.在该模型中,用三个独立的LBGK方程分别模拟速度场、温度场和浓度场,并通过Boussinesq近似将它们耦合起来.分析Soret数和Dufour数对方腔内部双扩散自然对流的影响,得到流线图、等温线图、等浓度线图、发热圆表面平均Nusselt数和平均Sherwood数.结果表明：Soret和Dufour效应对方腔内双扩散自然对流影响明显,不能忽略.     

Laminar Natural Convection in a Square Cavity with a Partition on the Heated Vertical Wall

Abstract

The laminar natural convection in an air-filled square cavity with a partition on the heated vertical wall was experimentally investigated. Temperature measurements and flow visualizations were performed for cases with heated and cooled vertical walls (corresponding to a global Grashof number Gr _H of approximately 1.4 × 10⁸) and non dimensional top wall temperatures θ _T of 0.57 (insulated) to 2.3. Experiments were performed with an aluminum partition with non dimensional height H _P /H of 0.0625 and 0.125 attached to the heated vertical wall at y/H = 0.65 and 0.95. The blockage effect and/or the thermal effect of the partition resulted in changes to the temperature and flow fields but were mainly limited to the vicinity of the partition. For the cases with the heated top wall, the change in the height of the partition at y/H = 0.95 resulted in changes to the ambient temperature outside the boundary layer due to the reduction of the size of the recirculating flow in the corner region. The changes in the partition height and the top wall temperature affected the blockage effect of the partition, resulting in the local Nusselt number near the corner region to be affected. The local Nusselt number over most of the heated vertical wall of the partitioned cavity (y/H < 0.7) was correlated to the local Rayleigh number in the form Nu = C · Ra ⁿ .     

多孔介质方腔内混合对流格子Boltzmann模拟

采用格子Boltzmann方法,对多孔介质方腔内的混合对流现象进行研究.方腔内部中心有一发热圆,径宽比D/L=0.4,冷流体从方腔左下角入口流进,从方腔左上角出口流出,四周壁面绝热.在普朗特数Pr=0.71和格拉晓夫数Gr=1.4×10⁴时,分析理查德森数Ri和达西数Da对发热圆表面平均努赛尔数Nu的影响.结果表明：Ri数位于10^-3~10范围内,Nu随Ri的增大而减小.Da越大,Ri对Nu的影响越显著;Da数位于10^-5~10^-2范围内,强制对流占主导的情况下（Ri≤0.1）,Nu随着多孔介质的Da的增大而增大.自然对流占主导的情况下（Ri=10）,Nu对Da的变化不敏感.     

不连续冷源布局多孔介质内热流耦合LBM数值模拟

为研究不连续冷源边界对内置发热体多孔介质方腔内传热及流动的影响,采用格子Boltzmann方法对REV尺度下多孔介质方腔内的自然对流进行计算,并研究瑞利数(Ra)、达西数(Da)、孔隙度对多孔介质方腔内传热流动的影响.发现Da对方腔内的流体流型影响很大,Da为10^-4时,多孔介质方腔内只有一个涡流,而Da为10^-2时,方腔内有两个涡流.增大Ra、Da、孔隙度可以提高冷源壁面的平均努赛尔数(Nu),增强散热效果,孔隙率对平均Nu影响程度和Da的大小有关.当冷源布置在壁面上方,壁面的平均Nu随Ra的增加剧烈变化,方腔处于高Ra条件下时,将冷源布置在边界的上方可以提高散热效果.6种布置方案中Case 6的散热效果最好.     

Growth of Ordered Domains in a Computer Model Alloy with Lattice Misfit

We study via Monte Carlo simulations the influence of elastic interactions on the ordering and decomposition of a two-dimensional model binary alloy with antiferromagnetic nearest and ferromagnetic next nearest neighbor type interactions following a quench into the coexistence region. The elastic interaction leads to the development of a platelet morphology for the segregated ordered and disordered regions. A length scale characterizing the coarsening process follows a law of the type R=a+bt ^1/3 with the growth b decreasing with the amount of ordered phase; this appears to be due to the presence of anti-phase boundaries between neighboring domains ordered on different sublattices which are difficult to eliminate. The application of uniaxial external stress results in rafting of the domains. Many of the simulation results are in agreement with experimentally observed effects in nickel-base superalloys.     

Development and validation of a new LBM-MRT hybrid model with enthalpy formulation for melting with natural convection

This article presents a new model to simulate melting with natural convection of a phase change material. For the phase change problem, the enthalpy formulation is used. Energy equation is solved by a finite difference method, whereas the fluid flow is solved by the multiple relaxation time (MRT) lattice Boltzmann method. The model is first verified and validated using the data from the literature. Then, the model is applied to a tall brick filled with a fatty acid eutectic mixture and the results are presented. The main results are (1) the spatial convergence rate is of second order, (2) the new model is validated against data from the literature and (3) the natural convection plays an important role in the melting process of the fatty acid mixture considered in our work.     

Structure and orientation of an intermetallic phase in a W-Ni-Co alloy

The aim of this investigation is crystal structure determination of an intermetallic phase formed in a W-Ni-Co alloy during a heat-treatment carried out at a temperature of 800°C. This intermetallic phase is expected to play a critical role on the final microstructure (fine tungsten particles in an FCC matrix that is present in between large tungsten grains) and thereby, on the properties of the alloy. 92W-5.3Ni-2.7Co alloy was prepared through powder metallurgy route (liquid phase sintering) followed by heat-treatment at 800°C for 5?h. The intermetallic phase formed at this temperature was characterised using transmission and scanning electron microscopes. The intermetallic phase was found to have orthorhombic crystal structure with Pnam (62) space group as determined using automated diffraction tomography along with precession electron diffraction. Chemical analysis in TEM suggested that the intermetallic phase is based on stoichiometry (Co,Ni)₂W. Orientation imaging of the phase was also carried out in TEM and EBSD to understand its evolution. Equiaxed or lath morphology of the intermetallic phase was found to depend on the crystallographic orientation relationship of the phase with the tungsten grains and the matrix phase.     

Phase-field simulation of dendritic growth in a binary alloy with thermodynamics data

This paper simulates the dendrite growth process during non-isothermal solidification in the Al-Cu binary alloy by using the phase-field model. The heat transfer equation is solved simultaneously. The thermodynamic and kinetic parameters are directly obtained from existing database by using the Calculation of Phase Diagram （CALPHAD） method. The effects of the latent heat and undercooling on the dendrite growth, solute and temperature profile during the solidification of binary alloy are investigated. The results indicate that the dendrite growing morphologies could be simulated realistically by linking the phase-field method to CALPHAD. The secondary arms of solidification dendritic are better developed with the increase of undercooling. Correspondingly, the tip speed and the solute segregation in solid-liquid interface increase, but the tip radius decreases.     

Monte Carlo simulations of magnetic particle suspensions with a simple assessment method for the particle overlap between magnetic spheroids

We have developed a simple assessment method for the overlap between spheroidal particles, which neither requires the complex manipulation of vectors and matrices that is indispensable in the ordinary methods, nor is based on a model potential. Moreover, we have developed an evaluation method for the interaction energy arising from the overlap of the steric layer coating spheroidal particles. This is based on a sphere-connected particle model, but some modifications are introduced in order to express an appropriate repulsive interaction energy at the deepest overlapping position. We have investigated the phase change in a magnetic spheroidal particle suspension for a two-dimensional system by means of Monte Carlo simulations. In the case of no external magnetic field, if the magnetic particle-particle interaction is sufficiently strong to favour cluster formation, long raft-like clusters tend to be formed in a dilute situation. With decreasing values of area fraction, a chain-like structure in a dense situation transforms into a raft-like structure within a narrow range of the particle area fraction. Similarly, the raft-like clusters are preferred in a weak applied magnetic field, but an increase in the field strength induces a phase change from a raft-like into a chain-like structure.

Highlights of the present paper:

1. A simple assessment method has been proposed for the overlap between two spheroidal particles.

2. The particle overlap assessment is free from a complex mathematical manipulation regarding vectors and matrices.

3. A modified sphere-connected model has been proposed in order to more accurately evaluate a repulsive interaction due to the overlap of the steric layers coating spheroidal particles.

4. 2D Monte Carlo simulations have been performed to elucidate the phenomenon of a phase change by magnetic spheroidal particles on a material plane surface.

5. A phase change between a raft-like and a chain-like aggregate structure is able to be controlled by the area fraction of particles and an external magnetic field.

     

相干点源照明时消球差光学系统的像场结构

利用矢量傅里叶变换和稳相法获得了相干点源位于光轴上任意位置时，消球差光学系统像场结构的积分表达式，详细研究了使用线偏振光照明时像平面上大物方孔径角对像场结构的影响.模拟计算表明，在像空间垂直于光轴的平面上，如果物方孔径角较大，磁场分布绕光轴旋转90°后不再与电场分布相同，电场能量密度、磁能量密度和玻印亭矢量分布的等高线始终近似为椭圆，并且物方孔径角是导致玻印亭矢量分布失去圆对称性的主要因素.同时当使用小像方孔径角时，电场能量密度分布形状的长轴方向垂直于物空间电场的振动方向；随着像方孔径角逐渐增大，电场能量密度分布形状的长轴方向将逐渐变为与物空间电场的振动方向相同.这些结论完全不同于以前理论所预测的结果. 关键词： 成像系统 像场结构 矢量傅里叶变换 稳相法