Kinetic Study Of The Rayleigh-Benard Flows

The present paper employs the direct simulation Monte Carlo (DSMC) method to study the Rayleigh-Benard flows, where the temperature ratio of the upper to lower plate is fixed to 0.1. For a Knudsen number (Kn) of 0.01, as the Rayleigh number (Ra) increases, the flow changes from the thermal conductive state to the convective state at about Ra=1700, and the calculated relation of heat flux through the lower plate versus Ra is in good agreement with classical experimental and theoretical results. For Kn=0.05, the thermal conductive state remains stable, and the increase of Ra cannot trigger thermal instability.

Complex Transition of Double-Diffusive Convection in a Rectangular Enclosure with Height-to-Length Ratio Equal to 4: Part I

This is the first part of direct numerical simulation (DNS) of double-diffusive convection in a slim rectangular enclosure with horizontal temperature and concentration gradients. We consider the case with the thermal Rayleigh number of 10^5, the Pradtle number of 1, the Lewis number of 2, the buoyancy ratio of composition to temperature being in the range of [0,1], and height-to-width aspect ration of 4. A new 7th order upwind compact scheme was developed for approximation of convective terms, and a three-stage third-order Runge-Kutta method was employed for time advancement. Our DNS suggests that with the buoyancy ratio increasing form 0 to 1, the flow of transition is a complex series changing fromthe steady to periodic, chaotic, periodic, quasi-periodic, and finally back to periodic. There are two types of periodic flow, one is simple periodic flow with single fundamental frequency (FF), and another is complex periodic flow with multiple FFs. This process is illustrated by using time-velocity histories, Fourier frequency spectrum analysis and the phase-space rajectories.

Heat Transfer of Aviation Kerosene at Supercritical Conditions

The heat transfer characteristics of China no. 3 kerosene were investigated experimentally and analytically under conditions relevant to a regenerative cooling system for scramjet applications. A test facility developed for the present study can handle kerosene in a temperature range of 300-1000 K, a pressure range of 2.6-5 MPa, and a mass How rate range of 10-100 g/s. In addition, the test section was uniquely designed such that both the wall temperature and the bulk fuel temperature were measured at the same location along the flowpath. The measured temperature distributions were then used to analytically deduce the local heat transfer characteristics. A 10-component kerosene surrogate was proposed and employed to calculate the fuel thermodynamic and transport properties that were required in the heat transfer analysis. Results revealed drastic changes in the fuel flow properties and heat transfer characteristics when kerosene approached its critical state. Convective heat transfer enhancement was also found as kerosene became supercritical. The heat transfer correlation in the relatively low-fuel-temperature region yielded a similar result to other commonly used jet fuels, such as JP-7 and JP-8, at compressed liquid states. In the high-fuel-temperature region, near and beyond the critical temperature, heat transfer enhancement was observed; hence, the associated correlation showed a more significant Reynolds number dependency.

千瓦级双包层光纤激光器冷却方案设计理论和实验研究

尽管双包层光纤激光器的散热性能好于传统的固体激光器的散热性能,光纤激光器中的热沉积仍然是限制提高其输出功率的重要因素.以双端抽运的400W双包层光纤激光器为实例,定量分析了光纤内的热沉积分布.根据所建立的散热模型,为了确保千瓦级双包层光纤激光器安全稳定的运行,抽运端附近的对流换热系数应大于2.8×10-2W·cm-2K-1.据此设计出高功率双包层光纤激光器抽运端冷却装置并成功应用在激光系统中,获得了千瓦级的激光输出.

Numerical analysis of LEC growth of GaAs with an axial magnetic field

A set of numerical analyses for momentum and heat transfer For a 3 in. (0.075 m) diameter Liquid Encapsulant Czochralski (LEC) growth of single-crystal GaAs with or without all axial magnetic field was carried Out using the finite-element method. The analyses assume a pseudosteady axisymmetric state with laminar floats. Convective and conductive heat transfers. radiative heat transfer between diffuse surfaces and the Navier-Stokes equations for both melt and encapsulant and electric current stream function equations Cor melt and crystal Lire considered together and solved simultaneously. The effect of the thickness of encapsulant. the imposed magnetic field strength as well as the rotation rate of crystal and crucible on the flow and heat transfer were investigated. (C) 2002 Published by Elsevier Science Ltd.

Chemical etching of a GaSb crystal incorporated with Mn grown by the Bridgman method under microgravity conditions

A GaSb crystal incorporated with Mn has been grown by the Bridgman method on the Polizon facility onboard the FOTON-M3 spacecraft. Structural defects and growth striations have been successfully revealed by the chemical etching method. By calculating various parameters of the convection, the striation patterns can be explained, and the critical value of the Taylor number, which characterizes the convective condition of the rotating magnetic field induced azimuthal flow, was shown. The stresses generated during crystal growth can be reflected by the observations of etch pit distribution and other structural defects. Suggestions for improving the space experiment to improve the quality of the crystal are given.

超临界压力下航空煤油水平管内对流换热特性数值研究

采用RNG(renormalization group)k-ε湍流模型和近壁区的Wolfstein一方程模型对超临界压力下大庆RP-3航空煤油在水平圆管内的流动和换热特性进行了数值研究.超临界压力下,由于航空煤油在拟临界点附近热物性的剧烈变化,浮升力将引起显著的二次流动.二次流动使得水平圆管的下表面湍流强度和对流换热增强,而上表面的湍流强度和对流换热减弱.最后分析了两种水平管内对流换热受浮升力影响判别标准在超临界流体中的适用性

超临界压力下航空煤油流动与传热特性试验

本文针对超燃冲压发动机再生冷却系统运行条件,实验研究了大庆RP-3煤油在超临界压力下的流动和传热特性,目的在于研究煤油在各种工况下的流动参数变化以及对流传热规律. 煤油通过二级煤油加热/输运系统加热,试验的煤油压力约2.6M~5.0MPa,油温约300~800 K. 相应的壁面热流密度为10~300 kW/m2. 通过油温与壁温的同步测量,结合非定常传热分析,获得了超临界压力下、亚/超临界温度范围内煤油的流动和传热物理参量的变化曲线

对流扩散方程的三阶迎风格式的数值摄动高精度重构

利用高智提出的数值摄动算法,把求解对流扩散方程常用三阶迎风格式(3-UDS)(粘性项和对流项分别用二阶中心格式和3-UDS离散)进行了高精度重构,包括使用离散单元内所有节点的全域重构和分别使用上下游节点的上下游重构,得到两类新的更高阶精度迎风差分格式,称为高的迎风差分格式(记作GUDS)。讨论了GUDS的数学性质,GUDS比原来的3-UDS精度显著提高;全域重构的GUDS和3-UDS均为条件稳定,一些上下游重构GUDS为绝对稳定。本文通过稳定性分析和四个算例(一维常系数、变系数、非线性及二维变系数对流扩散方程)的计算证实了GUDS的优良性质。上下游重构GUDS为避免在3-UDS中使用人工粘性提供了一条有效途径,适合于求解高Reynolds数线性和非线性问题。

对流扩散方程的绝对稳定高阶中心差分格式

将作者提出的数值摄动算法改进为区分离散单元内上游和下游并分别对通量进行高精度重构的双重数值摄动算法,与原(单重)摄动算法相比,双重摄动算法既提高了格式精度又明显扩大了格式的稳定域范围,利用双重摄动算法,即分别利用上游和下游基点变量的摄动重构将高阶流体力学关系及迎风机制耦合进二阶中心格式之中,由此构建了对流扩散方程的对网格Reynolds数的任意值均稳定(绝对稳定)高精度(四阶和八阶精度)三基点中心TVD差分格式,通过解析分析以及3个算例计算证实了构建格式的优良性能;3个算例包括一维线性、非线性(Burgers方程)和二维变系数对流扩散方程,数值计算表明:构建的格式在粗网格下不振荡,构建格式在粗网格时的最大误差L∞和均方误差L2与二阶中心格式在细网格时的相应误差一致,对线性方程,构建格式在细网格下可达到L2精度阶

Thermocapillary Migration of Deformable Bubbles at Moderate to Large Marangoni Number in Microgravity

Using the level-set method and the continuum interface model, the axisymmetric thermocapillary migration of gas bubbles in an immiscible bulk liquid with a temperature gradient at moderate to large Marangoni number is simulated numerically. Constant material properties of the two phases are assumed. Steady state of the motion can always be reached. The terminal migration velocity decreases monotonously with the increase of the Marangoni number due to the wrapping of isotherms around the front surface of the bubble. Good agreements with space experimental data and previous theoretical and numerical studies in the literature are evident. Slight deformation of bubble is observed, but no distinct influence on the motion occurs. It is also found that the influence of the convective transport of heat inside bubbles cannot be neglected at finite Marangoni number, while the influence of the convective transport of momentum inside bubbles may be actually negligible.

Linear spatio-temporal instability analysis of ice growth under a falling water film

A linear spatio-temporal stability analysis is conducted for the ice growth under a falling water film along an inclined ice plane. The full system of linear stability equations is solved by using the Chebyshev collocation method. By plotting the boundary curve between the linear absolute and convective instabilities (AI/CI) of the ice mode in the parameter plane of the Reynolds number and incline angle, it is found that the linear absolute instability exists and occurs above a minimum Reynolds number and below a maximum inclined angle. Furthermore, by plotting the critical Reynolds number curves with respect to the inclined angle for the downstream and upstream branches, the convectively unstable region is determined and divided into three parts, one of which has both downstream and upstream convectively unstable wavepackets and the other two have only downstream or upstream convectively unstable wavepacket. Finally, the effect of the Stefan number and the thickness of the ice layer on the AI/CI boundary curve is investigated.

超临界压力下航空煤油圆管流动和传热的数值研究

对超临界压力下大庆RP-3航空煤油在小管道内的流动、传热过程进行了数值研究. 湍流模拟采用了 RNG k-ε两方程模型和Wolfstein一方程模型结合的两层模型;同时, 采用煤油的10组分替代模型以及NIST Supertrapp程序库对大庆3号航空煤油的热物理和输运特性进行了确定. 圆管传热的计算条件为:入口压力4 MPa, 入口温度300 K, 质量流量范围:0.06~0.12 kg/s, 壁面热流密度范围:300~700 kW/m~2. 计算结果显示, 煤油的流动和传热特性比水、二氧化碳等简单化合物复杂得多. 在超临界压力下, 煤油的吸热升温导致其热物理特性以及流动特性均发生剧烈变化, 其中, 雷诺数沿管道方向上升了至少一个量级, 而普朗特数下降了一个量级. 在加热开始段, 煤油的对流传热系数迅速上升;当壁面温度超过其拟临界温度后, 对流传热系数略有所回落;随着煤油温度的进一步上升, 传热系数又得到明显增强. 计算表明, 煤油对流换热特性的变化与煤油复杂的高温热物理特性以及湍流流动在近壁区的增强和抑止有关

Short and long waves over a muddy seabed

The available experimental results have shown that in time-periodic motion the rheology of fluid mud displays complex viscoelastic behaviour. Based on the measured rheology of fluid mud from two field sites, we study the interaction of water waves and fluid mud by a two-layered model in which the water above is assumed to be inviscid and the mud below is viscoelastic. As the fluid-mud layer in shallow seas is usually much thinner than the water layer above, the sharp contrast of scales enables an approximate analytical theory for the interaction between fluid mud and small-amplitude waves with a narrow frequency band. It is shown that at the leading order and within a short distance of a few wavelengths, wave pressure from above forces mud motion below. Over a Much longer distance, waves are modified by the accumulative dissipation in mud. At the next order, infragravity waves owing to convective inertia (or radiation stresses) are affected indirectly by mud motion through the slow modulation of the short waves. Quantitative predictions are made for mud samples of several concentrations and from two different field sites.

Combined Effects of Hot Curing Conditions and Reaction Heat on Rubber Vulcanization Efficiency and Vulcanizate Uniformity

A mathematical model of the chemical kinetics of silicone rubber Vulcanization is developed, with the thermal effects being computed using the increment method, and the hot Vulcanization process estimated with the finite element method. The results show that the reaction heat of rubber vulcanization is important for energy saving, and that a proper curing medium temperature is important when considering both vulcanization efficiency and vulcanizate uniformity. The results also indicate that increases in the forced convective heat transfer coefficient have no significant effect above a certain level. The validity of the numerical model is indirectly proven by comparison with existing data.