Numerical simulation of pool boiling of a Lennard-Jones liquid

We performed a numerical simulation of pool boiling by a molecular dynamics model. In the simulation, a liquid composed of Lennard-Jones particles in a uniform gravitational field is heated by a heat source at the bottom of the system. The model successfully reproduces the change in regimes of boiling from nucleate boiling to film boiling with the increase of the heat source temperature. We present the pool boiling curve by the model, whose general behavior is consistent with those observed in experiments of pool boiling.     

Fire safety distances for open pool fires

Fire accidents that carry huge loss with them have increased in the previous two decades than at any time in the history. Hence, there is a need for understanding the safety distances from different fires with different fuels. Fire safety distances are computed for different open pool fires. Diesel, gasoline and hexane are used as fuels for circular pool diameters of 0.5 m, 0.7 m and 1.0 m. A large square pool fire of 4 m × 4 m is also conducted with diesel as a fuel. All the prescribed distances in this study are purely based on the thermal analysis. IR camera is used to get the thermal images of pool fires and there by the irradiance at different locations is computed. The computed irradiance is presented with the threshold heat flux limits for human beings.     

Flame characteristics of small-scale pool fires under low gravity environments

The flame characteristics of pool fires such as oscillatory frequency, flame height, and so on differ depending on gravity. To improve our understanding of the effects of gravity on flame characteristics, we conducted an experimental investigation on small-scale pool fires under different gravity levels to normal. The drop tower facility at Hirosaki University in Japan was used to obtain an arbitrary low gravity environment. Our measurements include analyses of oscillatory frequency and flame height under low gravity environments and flow visualization by a particle-track laser-sheet technique (PTLS) combined with a high-speed camera. These measurements revealed that the puffing frequency and flame height declined with decreasing gravity. The puffing frequency was summarized by the relationship between the Strouhal and Froude numbers. The flame height showed a tendency to peak at a certain level of gravity. The behavior of the flame height under low gravity did not correlate to the scaling prediction previously presented, while the normal velocity component to the flame sheet declined with decreasing gravity. The chemical reaction in the flame may be restricted due to decreasing oxygen supply. Consequently, the convective heat feedback to the liquid surface falls, causing the flame height to decline with decreasing gravity.     

Numerical simulation of experiments on turbulent natural convection of heat generating liquid in cylindrical pool

Available experimental data on heat transfer of a melt with volumetric heat generation are analyzed in order to use them for validating the computer codes that describe a core catcher. The problem for CFD simulation of the experiments on heat transfer by laminar and turbulent natural convection is described. Information that can be obtained from experiments for verifying the models of convective heat transfer in a melt is analyzed. The effect of variable viscosity on the integral heat flux is discussed. Calculation results are represented and compared with experimental data on temperature distribution and integral heat transfer. The calculations are in good agreement with the experiment. The results are numerically extrapolated to the range of Rayleigh numbers up to 7 · 10¹⁶. It is concluded that the CFD calculations with the κ-ɛ turbulence model can be used in problems concerned with analysis of melt convection in a core catcher.     

Cross relaxation in liquid methanol

The cross relaxation rate for intermolecular dipole-dipole interactions between methyl protons in liquid methanol at 0 °C was measured and compared to the rate predicted from MD simulations of the experimental system. The experimental and calculated values agree well, even though the translational diffusion coefficient and bulk viscosity of the sample are not well-predicted by the simulations.     

Conditional moment closure modelling of turbulent methanol jet flames

Conditional moment closure (CMC) predictions for a turbulent piloted jet diffusion flame of methanol in air at velocities of 66.2 and 90.3 m s^?1 are presented. Predictions are compared with the experimental joint Raman-Rayleigh-LIF results of Masri et al and laminar flamelet calculations. Three comprehensive chemical mechanisms (SKELETAL, GRI-Mech and SUBGRI) are used to represent the chemistry of the methanol flame. The SKELETAL mechanism shows the best agreement among the various mechanisms employed. It is found that the SUBGRI mechanism reduces computational cost in terms of memory and CPU time without compromising results where the focus is on the main reactive chemistry.

The k-ε-g turbulence model underpredicts the rate of mixing and the predicted flames are somewhat longer than that reported by experiment. In general, the CMC predictions for conditional mean temperature and species mass fractions are very good and show qualitative agreement with experiment. At downstream locations, the overall trends of predicted temperature and species concentration levels are similar to the upstream ones with the latter showing better agreement with the conditional measured levels. CMC predictions show the same order of agreement at higher velocities.

It is believed that the discrepancies on the fuel-rich side may be due to lack of consideration of the conditional fluctuations. The absence of a radially dependent CMC formulation, excluding differential diffusion effects and the inadequacy of the chemical mechanism may also account, partly, for the degree of discrepancy in the predictions.     

Critical heat load upon pool liquid boiling

The theoretical derivation of the formula for calculating heat load, given in the treatise titled “The Theory of Heat and Mass Transfer” (edited by A.I. Leont’ev, Bauman State Technical University, Moscow, 1997) is noted to be incorrect. A corrected deduction of the formula, performed using the dimensional method, is considered. A brief review of the whole problem is given. It is noted that the “second crisis of boiling” is due to the effect of the “spheroidal” state in liquids.     

Numerical modelling of foam Couette flows

A numerical computation based on a tensorial visco-elasto-plastic model based on continuous mechanics is compared to experimental measurements on liquid foams for a bidimensional Couette flow between two glass plates, both in stationary and transient cases. The main features of the model are elasticity up to a plastic yield stress, and viscoelasticity above it. The effect of the friction of the plates is taken into account. The numerical modelling is based on a small set of standard material parameters that are fully characterised. Shear localisation as well as acute transient observations are reproduced and agree with experimental measurements. The plasticity appears to be the fundamental mechanism of the localisation of the flow. Finally, the present approach could be extended from liquid foams to similar materials such as emulsions, colloids or wet granular materials, that exhibit localisation.     

冷热坩埚中电子枪加热金属熔池的数值分析

利用NavierStokes方程在一定的边界条件下求解了电子枪加热长槽形冷、热坩埚中的二维熔池流场,获得了熔池流场的速度和温度分布,并且详细地对比研究了冷、热坩埚中金属熔池的蒸发量与电子枪功率、束宽以及坩埚尺寸的关系.电子枪功率越高、束宽越小(除非小于05mm)、坩埚尺寸越大,则蒸发量越大,电子枪能量的有效利用率也越高.当电子枪束宽较大时,热坩埚的蒸发量较大,而当束宽较小时,冷坩埚的蒸发量较大.计算中对铁、铜、钆和铝等金属熔池进行了数值分析,获得了相似的结果. 关键词： NavierStokes方程 热毛细流 浮力流 金属熔池     

Energy transfer between laser filaments in liquid methanol

We demonstrate energy exchange between two filament-forming femtosecond laser beams in liquid methanol. Our results are consistent with those of previous works documenting coupling between filaments in air; in addition, we identify an unreported phenomenon in which the direction of energy exchange oscillates at increments in the relative pulse delay equal to an optical period (2.6 fs). Energy transfer from one filament to another may be used in remote sensing and spectroscopic applications utilizing femtosecond laser filaments in water and air.     

Numerical modelling of ion transport in flames

This paper presents a modelling framework to compute the diffusivity and mobility of ions in flames. The (n, 6, 4) interaction potential is adopted to model collisions between neutral and charged species. All required parameters in the potential are related to the polarizability of the species pair via semi-empirical formulas, which are derived using the most recently published data or best estimates. The resulting framework permits computation of the transport coefficients of any ion found in a hydrocarbon flame. The accuracy of the proposed method is evaluated by comparing its predictions with experimental data on the mobility of selected ions in single-component neutral gases. Based on this analysis, the value of a model constant available in the literature is modified in order to improve the model's predictions. The newly determined ion transport coefficients are used as part of a previously developed numerical approach to compute the distribution of charged species in a freely propagating premixed lean CH₄/O₂ flame. Since a significant scatter of polarizability data exists in the literature, the effects of changes in polarizability on ion transport properties and the spatial distribution of ions in flames are explored. Our analysis shows that changes in polarizability propagate with decreasing effect from binary transport coefficients to species number densities. We conclude that the chosen polarizability value has a limited effect on the ion distribution in freely propagating flames. We expect that the modelling framework proposed here will benefit future efforts in modelling the effect of external voltages on flames. Supplemental data for this article can be accessed at http://dx.doi.org/10.1080/13647830.2015.1090018.     

Effect of pressure on hydrogen bonding in liquid methanol

We have investigated the effect of pressure on the hydrogen bonding in liquid methanol using Raman spectroscopy. Specifically, we have measured the OH and CO stretching modes and assigned the bands, in agreement with recent IR and crossed molecular beam experiments on methanol clusters. At about 7 to 8 kbar, we note indications that the intrinsic nature of the methanol clusters in our samples has changed. Our results provide support for and extend conclusions derived from Monte Carlo simulations, explain anomalies observed by previous researchers, and provide new insights into general hydrogen-bonding phenomena.     

Geometric modelling of viscosity of copper-containing liquid alloys

In this work, viscosities of ternary Au–Ag–Cu and Al–Cu–Si liquid alloys have been calculated as a function of gold, aluminium and copper compositions for the sections Au–Ag–Cu (x_Ag/x_Cu = 0.543 at 1373 K), Al_x(Cu₅₀–Si₅₀)_(1–x) and Cu_x(Al₅₀–Si₅₀)_(1–x) at 1375 K using Chou’s general solution model, Muggianu, Kohler, Toop, Hillert, Budai et al., Kozlov et al., Schick et al. and Kaptay et al. models. The present study finds that a comparison of the predicted values of viscosities associated with the geometric and physical models indicate good mutual agreement. The Muggianu model indicates the best agreement with the results obtained for Au–Ag–Cu and Al_x–Cu₅₀–Si₅₀ alloy systems and the Kaptay et al. model, which is a physical model, indicates the best agreement with the results obtained for Al₅₀–Cu_x–Si₅₀.     

电子枪加热合金熔池的数值分析

利用数值方法求解了电子枪加热二维合金熔池的流体力学方程组，获得了熔池流场的速度、温度和浓度分布.并且详细地研究了各种特征参量（如雷诺数、普朗特数、格拉晓夫数和李斯特数）对合金熔池流动特性的影响，一般地，雷诺数和普朗特数越小、而格拉晓夫数和李斯特数越大，则熔池中溶剂金属的蒸发速率就越大.在此基础上提出参杂适当浓度的溶质金属可以提高某些难熔难挥发金属的蒸发速率，计算中对参杂了钨、钕、锡和铬等金属的钍、铂、铱和钌等合金熔池进行了数值分析，都相应提高了溶剂金属的蒸发速率. 关键词： 液态合金 毛细流 浮力流 熔池     

Numerical modelling of anisotropy in 4H－SiC MESFET‘s

A numerical model for 4H-SiC MESFET anisotropy is presented in this paper and the device performances, such as breakdown, temperature and transient characteristics, are demonstrated. The simulation results show obvious effects of the anisotropy for 4H-SiC and are in better accordance with the experimental results. The anisotropy for 4H-SiC should be involved in the device design to acquire better performances.