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

本文针对超燃冲压发动机再生冷却系统运行条件,实验研究了大庆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.

Chain propagation kinetics on melt grafting reaction

Graft copolymerization in the molten state is of fundamental importance as a probe of chemical modification and reactive compatibilization. However, few grafting kinetic studies on reactive extrusion have been carried out because of the inherent difficulties, as expected. In this work, we have studied chain propagation kinetics on melt grafting using pre-irradiated linear low density polyethylene (LLDPE) and three monomers, acrylic acid (AA), methacrylic acid (MAA), and methyl methacrylate (MMA), as the model system. We measured the apparent chain propagation rate coefficients of grafting (k(p,g)) and homopolymerization (k(p,h)) at an initial stage for the melt grafting by FT-IR spectroscopy and electron spin resonance spectroscopy. It was observed that the convective mixing affected the rate coefficients. The magnitude of k(p,h) and k(p,g) were in the same order, but k(p,h) was slightly larger than k(p,g) The k(p,g) of the three grafting systems increased in the order: LLDPE/MMA < LLDPE/MAA < LLDPE/AA. These results are explained in terms of phase separation, solubility, and inherent reactivity of the monomer.

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

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.

Drainage induced convection rolls in foams

Hutzler, S., Cox, S.J., Janiaud, E. and Weaire, D. (2007) Drainage induced convection rolls in foams. Colloids and Surfaces A: Physicochemical and Engineering Aspects Volume 309, Issues 1-3, 1 November 2007, Pages 33-37 A Collection of Papers Presented at the 6th Eufoam Conference, Potsdam, Germany, 2-6 July, 2006 Sponsorship: European Space Agency (14914/02/NL/SH, 14308/00/NL/SG) (AO-99-031) CCN 002 MAP Project AO-99-075); Science Foundation Ireland (RFP 05/RFP/PHY0016); Royal Society; UWA Learned Societies.

HEATING IN VASCULAR TISSUE AND FLOW-THROUGH TISSUE PHANTOMS INDUCED BY FOCUSED ULTRASOUND

High intensity focused ultrasound (HIFU) can be used to control bleeding, both from individual blood vessels as well as from gross damage to the capillary bed. This process, called acoustic hemostasis, is being studied in the hope that such a method would ultimately provide a lifesaving treatment during the so-called "golden hour", a brief grace period after a severe trauma in which prompt therapy can save the life of an injured person. Thermal effects play a major role in occlusion of small vessels and also appear to contribute to the sealing of punctures in major blood vessels. However, aggressive ultrasound-induced tissue heating can also impact healthy tissue and can lead to deleterious mechanical bioeffects. Moreover, the presence of vascularity can limit one’s ability to elevate the temperature of blood vessel walls owing to convective heat transport. In an effort to better understand the heating process in tissues with vascular structure we have developed a numerical simulation that couples models for ultrasound propagation, acoustic streaming, ultrasound heating and blood cooling in Newtonian viscous media. The 3-D simulation allows for the study of complicated biological structures and insonation geometries. We have also undertaken a series of in vitro experiments, in non-uniform flow-through tissue phantoms, designed to provide a ground truth verification of the model predictions. The calculated and measured results were compared over a range of values for insonation pressure, insonation time, and flow rate; we show good agreement between predictions and measurements. We then conducted a series of simulations that address two limiting problems of interest: hemostasis in small and large vessels. We employed realistic human tissue properties and considered more complex geometries. Results show that the heating pattern in and around a blood vessel is different for different vessel sizes, flow rates and for varying beam orientations relative to the flow axis. Complete occlusion and wall- puncture sealing are both possible depending on the exposure conditions. These results concur with prior clinical observations and may prove useful for planning of a more effective procedure in HIFU treatments.

Scale-wise evolution of rainfall probability density functions fingerprints the rainfall generation mechanism

© 2010 by the American Geophysical Union.The cross-scale probabilistic structure of rainfall intensity records collected over time scales ranging from hours to decades at sites dominated by both convective and frontal systems is investigated. Across these sites, intermittency build-up from slow to fast time-scales is analyzed in terms of heavy tailed and asymmetric signatures in the scale-wise evolution of rainfall probability density functions (pdfs). The analysis demonstrates that rainfall records dominated by convective storms develop heavier-Tailed power law pdfs toward finer scales when compared with their frontal systems counterpart. Also, a concomitant marked asymmetry build-up emerges at such finer time scales. A scale-dependent probabilistic description of such fat tails and asymmetry appearance is proposed based on a modified q-Gaussian model, able to describe the cross-scale rainfall pdfs in terms of the nonextensivity parameter q, a lacunarity (intermittency) correction and a tail asymmetry coefficient, linked to the rainfall generation mechanism.

Lightning mapping observation of a terrestrial gamma-ray flash

We report the observation with the North Alabama Lightning Mapping Array (LMA) related to a terrestrial gamma-ray flash (TGF) detected by RHESSI on 26 July 2008. The LMA data explicitly show the TGF was produced during the initial development of a compact intracloud (IC) lightning flash between a negative charge region centered at about 8.5 km above sea level (-22C temperature level) a higher positive region centered at 13 km, both confined to the convective core of an isolated storm in close proximity to the RHESSI footprint. After the occurrence of an LMA source with a high peak power (26 kW), the initial lightning evolution caused an unusually large IC current moment that became detectable 2 ms after the first LMA source and increased for another 2 ms, during which the burst of gamma-rays was produced. This slowly building current moment was most likely associated with the upward leader progression, which produced an uncommonly large IC charge moment change (+90 Ckm) in 3 ms while being punctuated by a sequence of fast discharge. These observations suggest that the leader development may be involved in the TGF production. Copyright © 2010 by the American Geophysical Union.

Using fractal downscaling of satellite precipitation products for hydrometeorological applications

The objective of spatial downscaling strategies is to increase the information content of coarse datasets at smaller scales. In the case of quantitative precipitation estimation (QPE) for hydrological applications, the goal is to close the scale gap between the spatial resolution of coarse datasets (e.g., gridded satellite precipitation products at resolution L × L) and the high resolution (l × l; L»l) necessary to capture the spatial features that determine spatial variability of water flows and water stores in the landscape. In essence, the downscaling process consists of weaving subgrid-scale heterogeneity over a desired range of wavelengths in the original field. The defining question is, which properties, statistical and otherwise, of the target field (the known observable at the desired spatial resolution) should be matched, with the caveat that downscaling methods be as a general as possible and therefore ideally without case-specific constraints and/or calibration requirements? Here, the attention is focused on two simple fractal downscaling methods using iterated functions systems (IFS) and fractal Brownian surfaces (FBS) that meet this requirement. The two methods were applied to disaggregate spatially 27 summertime convective storms in the central United States during 2007 at three consecutive times (1800, 2100, and 0000 UTC, thus 81 fields overall) from the Tropical Rainfall Measuring Mission (TRMM) version 6 (V6) 3B42 precipitation product (~25-km grid spacing) to the same resolution as the NCEP stage IV products (~4-km grid spacing). Results from bilinear interpolation are used as the control. A fundamental distinction between IFS and FBS is that the latter implies a distribution of downscaled fields and thus an ensemble solution, whereas the former provides a single solution. The downscaling effectiveness is assessed using fractal measures (the spectral exponent β, fractal dimension D, Hurst coefficient H, and roughness amplitude R) and traditional operational scores statistics scores [false alarm rate (FR), probability of detection (PD), threat score (TS), and Heidke skill score (HSS)], as well as bias and the root-mean-square error (RMSE). The results show that both IFS and FBS fractal interpolation perform well with regard to operational skill scores, and they meet the additional requirement of generating structurally consistent fields. Furthermore, confidence intervals can be directly generated from the FBS ensemble. The results were used to diagnose errors relevant for hydrometeorological applications, in particular a spatial displacement with characteristic length of at least 50 km (2500 km2) in the location of peak rainfall intensities for the cases studied. © 2010 American Meteorological Society.