Thermal effects on Rayleigh-Marangoni-Benard instability in a system of superposed fluid and porous layers

Thermal effects of the heat transfer at free surface (represented by Biot number) on the Rayleigh-Marangoni-Benard instability in a system of liquid-porous layers with top free surface are investigated numerically. The results indicate that this thermal effect can evidently lead to the mode transition of convection, which is overlooked in previous works. (C) 2010 Elsevier Ltd. All rights reserved.

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

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

微重力气液两相流动与池沸腾传热

综述了近年来中国科学院微重力重点实验室(国家微重力实验室)完成的一系列微重力气液两相流动与池沸腾传热方面的地基实验、飞行实验和理论研究等方面获得的主要成果.在微重力气液两相流动方面, 提出了半理论Weber数模型用于预测微重力条件下气液两相弹一环状流转换, 并采用Monte Carlo方法, 针对气泡初始尺寸对泡一弹状流转换的影响进行数值研究.通过俄罗斯"和平号"空间站与IL-76失重飞机实验, 获得了微重力下的气液两相流型图, 与此同时在地面利用小尺度毛细管模型模拟了微重力气液两相流动特征.实验测量了微重力气液两相流压降, 并基于微重力流动特性建立了一个泡状流压降关联模型.在微重力池沸腾传热方面, 利用我国返回式卫星完成了两次空间实验, 其中, 第22颗返回式卫星搭载铂丝表面R113池沸腾实验采用控制温度的稳态加热方式, 而实践8号育种卫星搭载平面FC-72池沸腾实验则采用控制加热电压的准稳态加热方式.同时, 还进行了地面常重力和落塔短时微重力条件下的对比实验研究.观察到丝状加热表面微重力时轻微的传热强化现象, 而平板加热表面微重力核态池沸腾低热流时传热强化、高热流时传热恶化.微重力实验中观察到气泡脱落前存在横向运动现象, 据此分析了气泡行为与传热之间关系, 并提出了一个预测丝状加热表面气泡脱落直径的半理论模型.旨在对相关领域的进一步发展和空间两相流系统的应用提供数据及理论支持

Pool Boiling on Wire in Microgravity - What Have We Learned

In the past years, steady pool boiling of degassed R113 on thin platinum wires has been studied systematically in our lab, including experiments in long-term microgravity aboard RS-22, in short-term microgravity in the Drop Tower Beijing / NMLC, and in normal gravity on the ground. Slight enhancement of nucleate boiling heat transfer is observed in microgravity, while dramatic changes of bubble behaviors are much evident. The value of CHF in microgravity is lower than that in normal gravity, but it can be predicted well by the Lienhard-Dhir correlation, although the dimensionless radius in the present case is far beyond its initial application range. The scaling of CHF with gravity is thus much different from the traditional viewpoint. Considering the influence of the Marangoni effects, the different characteristics of bubble behaviors in microgravity have been explained. A new bubble departure model has also been proposed, which can predict the whole observation both in microgravity and in normal gravity.

Thermally induced evolution of phase transformations in gas hydrate sediment

Thermally induced evolution of phase transformations is a basic physical-chemical process in the dissociation of gas hydrate in sediment (GHS). Heat transfer leads to the weakening of the bed soil and the simultaneous establishment of a time varying stress field accompanied by seepage of fluids and deformation of the soil. As a consequence, ground failure could occur causing engineering damage or/and environmental disaster. This paper presents a simplified analysis of the thermal process by assuming that thermal conduction can be decoupled from the flow and deformation process. It is further assumed that phase transformations take place instantaneously. Analytical and numerical results are given for several examples of simplified geometry. Experiments using Tetra-hydro-furan hydrate sediments were carried out in our laboratory to check the theory. By comparison, the theoretical, numerical and experimental results on the evolution of dissociation fronts and temperature in the sediment are found to be in good agreement.

等离子体电弧炉内流动与传热数值模拟

本文采用磁流体动力学(MHD)模型对直流等离子体自由燃烧电弧和电弧炉内部的流动与传热进行了数值模拟研究. 通过对基于磁矢量势描述的电磁场方程组和流体力学方程组的耦合迭代计算, 求解得到了流体的温度场和速度场等, 计算结果清晰地反映出等离子体电弧的高温阴极射流现象, 并与同行的实验和数值结果进行了对比. 本模拟方法和结果对于电弧炉的工业应用和优化设计有重要的指导意义.

Numerical Simulation of a Low-power Hydrazine Arcjet Thruster

A modeling study is conducted to investigate the plasma flow and heat transfer characteristics of low-power (kW class) arc-heated thrusters (arcjets) with 2:1 hydrogen/nitrogen to simulate decomposed hydrazine as the propellant. The all-speed SIMPLE algorithm is employed to solve the governing equations, which take into account the effects of compressibility, the Lorentz force and Joule heating, as well as the temperature- and pressure-dependence of the gas properties. Typical computed results about the temperature, velocity and Mach number distributions within arcjet thruster are presented for the case with arc current of 9 A and inlet stagnant pressure of 3.3×105 Pa to show the flow and heat transfer characteristics. It is found that the propellant is heated mainly in the near-cathode and constrictor region, with the highest plasma temperature appearing near the cathode tip, and the flow transition from the subsonic to supersonic regime occurs within the constrictor region. The effect of gas viscosity on the plasma flow within arcjet thruster is examined by an additional numerical test using artificially reduced values of gas viscosity. The test results show that the gas viscosity appreciably affects the plasma flow and the performance of the arcjet thruster for the cases with the hydrazine or hydrogen as the propellant. The integrated axial Lorentz force in the thruster nozzle is also calculated and compared with the thrust force of the arcjet thruster. It is found that the integrated axial Lorentz force is much smaller than the thrust force for the low-power arcjet thruster. Modeling results for the NASA 1-kW class arcjet thruster with simulated hydrazine as the propellant are found to be reasonably consistent with available experimental data.

Shallow gas in the muddy sediments of Eckernförde Bay, Germany

The seafloor of central Eckernförde Bay is characterised by soft muddy sediments that contain free methane gas. Bubbles of free gas cause acoustic turbidity which is observed with acoustic remote sensing systems. Repeated surveys with subbottom profiler and side scan sonar revealed an annual period both of depth of the acoustic turbidity and backscatter strength. The effects are delayed by 3–4 months relative to the atmospheric temperature cycle. In addition, prominent pockmarks, partly related to gas seepage, were detected with the acoustic systems. In a direct approach gas concentrations were measured from cores using the gas chromatography technique. From different tests it is concluded that subsampling of a core should start at its base and should be completed as soon as possible, at least within 35 min after core recovery. Comparison of methane concentrations of summer and winter cores revealed no significant seasonal variation. Thus, it is concluded that the temperature and pressure influences upon solubility control the depth variability of acoustic turbidity which is observed with acoustic remote sensing systems. The delay relative to the atmospheric temperature cycle is caused by slow heat transfer through the water column. The atmospheric temperature cycle as ‘exiting function’ for variable gas solubility offers an opportunity for modelling and predicting the depth of the acoustic turbidity. In practice, however, small-scale variations of, e.g., salinity, or gas concentration profile in the sediment impose limits to predictions. In addition, oceanographic influences as mixing in the water column, variable water inflow, etc. are further complications that reduce the reliability of predictions.

In vitro study on the binding of neutral red to bovine serum albumin by molecular spectroscopy

In this paper, the binding of neutral red (NR) to bovine serum albumin (BSA) under physiological conditions has been studied by spectroscopy method including fluorescence, circular dichroism (CD) and Fourier transform infrared (FT-IR) spectroscopy. The Stern-Volmer fluorescence quenching constant (K-SV), binding constant (K-b) and the number of binding sites (It) were measured by fluorescence quenching method. Fluorescence experiments were also performed at different ionic strengths. It was found K-SV was ionic strength dependent, which indicated the electrostatic interactions were part of the binding forces. The distance r between donor (BSA) and acceptor (NR) was obtained according to Foster's non-radiative energy transfer theory. CD spectroscopy and FT-IR spectroscopy were used to investigate the structural information of BSA molecules on the binding of NR, and the results showed no change of BSA conformation in our experimental conditions.

One century record of contamination by polycyclic aromatic hydrocarbons and polychlorinated biphenyls in core sediments from the southern Yellow Sea

Sixteen polycyclic aromatic hydrocarbons (PAHs) and 28 polychlorinated biphenyls (PCBs) were measured at a 2-cm interval in a core sample from the middle of the southern Yellow Sea for elucidating their historical variations in inflow and sources. The chronology was obtained using the Pb-210 method. PAHs concentrations decreased generally with depth and two climax values occurred in 14-16 cm and 20-22 cm layers, demonstrating that the production and usage of PAHs might reach peaks in the periods of 1956-1962 and 1938-1944. The booming economy and the navy battles of the Second World War might explain why the higher levels were detected in the two layers. The result of principal component analysis (PCA) revealed that PAHs were primarily owing to the combustion product. Down-cored variation of PCB concentrations was complex. Higher concentrations besides the two peaks being the same as PAHs were detected from 4 to 8 cm, depositing from 1980 to 1992, which probably resulted from the disposal of the out-dated PCB-containing equipment. The average Cl percentage of PCBs detected was similar to that of the mixture of Aroclor 1254 and 1242, suggesting they might origin from the dielectrical and heat-transfer fluid. The total organic carbon (TOC) content played a prevalent role in the adsorption of high molecular weight PAHs (>= 4-ring), while no obvious relationship among total PCBs, the concentration of congeners, and TOC was found.

Observation and modeling for terrestrial processes in alpine meadow

The water-heat transfer process between land and atmosphere in Haibei alpine meadow area has been systematically observed. A multi-layer coupling model for land-atmosphere interaction was presented with special attention paid to the moisture transfer in leaf stomata under unsaturated condition. A profound investigation on the physical process of turbulent transfer inside the vegetation has been performed with a revised formula of water absorption for root system. The present model facilitates the study of vertically distributed physical variables in detail. Numerical simulation was conducted according to the transfer process of Kinesia humility meadow in the area of Haibei Alpine Meadow Ecosystem Station, CAS. The calculated results agree well with observation.

温度、压力耦合背景下多相流体运移的数学模型研究

The Mathematical modeling of multiphase fluid flow is an important aspect of basin simulation, and also is a topic of geological frontier. Based on coupling relation of temperature, pressure and fluid flow, this dissertation discusses the modeling which conform to geological regularities of fluid migration. The modeling that is multi-field and multiphase includes heat transport equation, pressure evolvement equation, solution transport equation and fluid transport equation. The finite element method is effective numerical calculation methods. Author applies it to solve modeling and implements the finite element program, and the modeling is applied to Ying-Qiong Basin. The channels of fluid vertical migration are fault, fracture and other high penetrability area. In this thesis, parallel fracture model and columnar channel model have been discussed, and a characteristic time content and a characteristic space content been obtained to illustrate the influences of stratigraphic and hydrodynamic factors on the process. The elliptoid fracture model is established and its approximately solution in theory is gotten. Three kinds of modeling are applied to analyze the transient variation process of fluid pressure in the connected permeable formations. The elliptoid fracture model is the most similar geology model comparing with the other fracture models so the research on this fracture model can enhance the understanding to fluid pressure. In the non-hydrodynamic condition, because of the difference between water density and nature gas density, nature gas can migrate upon by float force. A one-dimension mathematical model of nature gas migration by float force is established and also applied to analyze the change in the saturation of gas. In the process of gas migration its saturation is non-continuous. Fluid flow is an important factor which influences the distribution of the temperature-field, the change of temperature can influence fluid property (including density, viscidity, and solubility),a nd the temperature field has coupling relations to the fluid pressure field. In this dissertation one-dimension and two-dimension thermal convection modeling is developed and also applied to analyze convective and conductive heat transfer. Author has established one-dimension and two-dimension mathematical modeling in which fluid is a mixture of water and nature gas based on the coupling relation between temperature and pressure, discussed mixture fluid convection heat transfer in different gas saturation, and analyzed overpressure form mechanism. Based on geothermal abnormity and pore pressure distribution in Dongfong 1-1, Yinggehai Basin, South China Sea, one-dimension mathematical modeling of coupling temperature and pressure is established. The modeling simulates the process that fluid migrates from deep to shallow and overpressure forms in shallow. When overpressure is so large that fractures appear and overpressure is released. As deep fluid flow to shallow, the high geothermal then forms in shallow. Based on the geological characteristics in Ya13-1, two-dimension mathematical modeling of coupling temperature and pressure is established. Fluid vertically flows in fault and then laterally migrates in reservoir. The modeling simulates the geothermal abnormity and pore pressure distribution in reservoir.

地幔对流理论及其在解释地球物理观测的应用

With the improving of mantle convection theory, the developing of computing method and increasing of the measurement data, we can numerically simulate more clearly about the effects on some geophysical observed phenomenons such as the global heat flow and global lithospheric stress field in the Earth's surface caused by mantle convection, which is the primary mechanism for the transport of heat from the Earth's deep interior to its surface and the underlying force mechanism of dynamics in the Earth.Chapter 1 reviews the historical background and present research state of mantle convection theory.In Chapter 2, the basic conception of thermal convection and the basic theory about mantle flow.The effects on generation and distribution of global lithospheric stres s field induced by mantle flow are the subject of Chapter 3. Mantle convection causes normal stress and tangential stresses at the bottom of the lithosphere, and then the sublithospheric stress field induces the lithospheric deformation as sixrface force and results in the stress field within the lithosphere. The simulation shows that the agreement between predictions and observations is good in most regions. Most of subduction zones and continental collisions are under compressive. While ocean ridges, such as the east Pacific ridge, the Atlantic ridge and the east African rift valley, are under tensile. And most of the hotspots preferentially occur in regions where calculated stress is tensile. The calculated directions of the most compressive principal horizontal stress are largely in accord with that of the observation except for some regions such as the NW-Pacifie subduction zone and Qinghai-Tibet Plateau, in which the directions of the most compressive principal horizontal stress are different. It shows that the mantel flow plays an important role in causing or affecting the large-scale stress field within the lithosphere.The global heat flow simulation based on a kinematic model of mantle convection is given in Chapter 4. Mantle convection velocities are calculated based on the internal loading theory at first, the velocity field is used as the input to solve the thermal problem. Results show that calculated depth derivatives of the near surface temperature are closely correlated to the observed surface heat flow pattern. Higher heat flow values around midocean ridge systems can be reproduced very well. The predicted average temperature as a function of function of depth reveals that there are two thermal boundary layers, one is close to the surface and another is close to the core-mantle boundary, the rest of the mantle is nearly isothermal. Although, in most of the mantle, advection dominates the heat transfer, the conductive heat transfer is still locally important in the boundary layers and plays an important role for the surface heat flow pattern. The existence of surface plates is responsible for the long wavelength surface heat flow pattern.In Chapter 5, the effects on present-day crustal movement in the China Mainland resulted from the mantle convection are introduced. Using a dynamic method, we present a quantitative model for the present-day crustal movement in China. We consider not only the effect of the India-Eurasia collision, the gravitational potential energy difference of the Tibet Plateau, but also the contribution of the shear traction on the bottom of the lithosphere induced by the global mantle convection. The comparison between our results and the velocity field obtained from the GPS observation shows that our model satisfactorily reproduces the general picture of crustal deformation in China. Numerical modeling results reveal that the stress field on the base of the lithosphere induced by the mantle flow is probably a considerable factor that causes the movement and deformation of the lithosphere in continental China with its eflfcet focuing on the Eastern China A numerical research on the small-scale convection with variable viscosity in the upper mantle is introduced in Chapter 6. Based on a two-dimensional model, small-scale convection in the mantle-lithosphere system with variable viscosity is researched by using of finite element method. Variation of viscosity in exponential form with temperature is considered in this paper The results show that if viscosity is strongly temperature-dependent, the upper part of the system does not take a share in the convection and a stagnant lid, which is identified as lithosphere, is formed on the top of system because of low temperature and high viscosity. The calculated surface heat flow, topography and gravity anomaly are associated well with the convection pattern, namely, the regions with high heat flow and uplift correspond to the upwelling flow, and vice versa.In Chapter 7, we give a brief of future research subject: The inversion of lateral density heterogeneity in the mantle by minimizing the viscous dissipation.

单（多）井抽灌对浅部地温场的影响研究

The technique of energy extraction using groundwater source heat pumps, as a sustainable way of low-grade thermal energy utilization, has widely been used since mid-1990's. Based on the basic theories of groundwater flow and heat transfer and by employing two analytic models, the relationship of the thermal breakthrough time for a production well with the effect factors involved is analyzed and the impact of heat transfer by means of conduction and convection, under different groundwater velocity conditions, on geo-temperature field is discussed.A mathematical model, coupling the equations for groundwater flow with those for heat transfer, was developed. The impact of energy mining using a single well system of supplying and returning water on geo-temperature field under different hydrogeological conditions, well structures, withdraw-and-reinjection rates, and natural groundwater flow velocities was quantitatively simulated using the finite difference simulator HST3D. Theoretical analyses of the simulated results were also made. The simulated results of the single well system indicate that neither the permeability nor the porosity of a homogeneous aquifer has significant effect on the temperature of the production segment provided that the production and injection capability of each well in the aquifers involved can meet the designed value. If there exists a lower permeable interlayer, compared with the main aquifer, between the production and injection segments, the temperature changes of the production segment will decrease. The thicker the interlayer and the lower the interlayer permeability, the longer the thermal breakthrough time of the production segment and the smaller the temperature changes of the production segment. According to the above modeling, it can also be found that with the increase of the aquifer thickness, the distance between the production and injection screens, and/or the regional groundwater flow velocity, and/or the decrease of the production-and-reinjection rate, the temperature changes of the production segment decline. For an aquifer of a constant thickness, continuously increase the screen lengths of production and injection segments may lead to the decrease of the distance between the production and injection screens, and the temperature changes of the production segment will increase, consequently.According to the simulation results of the single well system, the parameters, that can cause significant influence on heat transfer as well as geo-temperature field, were chosen for doublet system simulation. It is indicated that the temperature changes of the pumping well will decrease as the aquifer thickness, the distance between the well pair and/or the screen lengths of the doublet increase. In the case of a low permeable interlayer embedding in the main aquifer, if the screens of the pumping and the injection wells are installed respectively below and above the interlayer, the temperature changes of the pumping well will be smaller than that without the interlay. The lower the permeability of the interlayer, the smaller the temperature changes. The simulation results also indicate that the lower the pumping-and-reinjection rate, the greater the temperature changes of the pumping well. It can also be found that if the producer and the injector are chosen reasonably, the temperature changes of the pumping well will decline as the regional groundwater flow velocity increases. Compared with the case that the groundwater flow direction is perpendicular to the well pair, if the regional flow is directed from the pumping well to the injection well, the temperature changes of the pumping well is relatively smaller.Based on the above simulation study, a case history was conducted using the data from an operating system in Beijing. By means of the conceptual model and the mathematical model, a 3-D simulation model was developed and the hydrogeological parameters and the thermal properties were calibrated. The calibrated model was used to predict the evolution of the geo-temperature field for the next five years. The simulation results indicate that the calibrated model can represent the hydrogeological conditions and the nature of the aquifers. It can also be found that the temperature fronts in high permeable aquifers move very fast and the radiuses of temperature influence are large. Comparatively, the temperature changes in clay layers are smaller and there is an obvious lag of the temperature changes. According to the current energy mining load, the temperature of the pumping wells will increase by 0.7°C at the end of the next five years. The above case study may provide reliable base for the scientific management of the operating system studied.