Numerical simulation of g-jitter effects on half floating zone convection under microgravity environment

The g-jitter effects on the thermocapillary convection in liquid bridge of floating half zone were studied by numerical simulation for unsteady and axi-symmetric model in the cylindrical coordinate system. The g-jitter field was given by a steady microgravity field in addition to an oscillatory low-gravity field, and the effects on the flow field, temperature distribution and free surface deformation were analyzed numerically.

Effects of high-frequency vibration on critical marangoni number

The influence of vibration on thermocapillary convection and critical Marangoni number in liquid bridge of half floating zone was discussed for the low frequency range 0.4-1.5 Hz and the intermediate frequency range 2.5-15 Hz in our previous papers. This paper extends the study to high frequency range 15-100Hz. This ground based experiment was completed on the deck of an electromagnetic vibration machine. The results of our experiment shows when the frequency of the applied acceleration is high enough, the amplitude of the time varying part of the temperature response is disappear and the shape of the free surface of the liquid bridge exhibits no fluctuations due to inertia. The critical Marangoni number which is defined to describe the transitions from a peroidical convection in response to vibration to an oscillatory convection due to internal instability is nearly the same as the critical Marangoni number for oscillatory flow in the absence of vibration.

Experimental study on oscillatory thermocapillary convection

The transition processes from steady flow into oscillatory flow in a liquid bridge of the half floating zone are studied experimentally. Two methods of noncontacted diagnoses are developed to measure the distribution of critical Marangoni numbers described by the onset of the oscillation st the free surface of the liquid bridge.The experimental results obtained for both cases of the upper rod heated and the lower rod heated agree with the prediction by Rayleigh's instability theory.The sensitive relations between the relatively fat or slender liquid bridge and the onset of oscillatory convection are also discussed to reveal the insight of the pressure distribution near the free surface. The experiments have been performed in a small liquid bridge, where the Bond number is much smaller than 1, and the results can be used to simulate the experiment in the microgravity environment.

Effects of Thermocapillary Convection on the Growth Process in Industrial 8-inch Silicon Czochralski Growth

Czochralski (CZ) crystal growth process is a widely used technique in manufacturing of silicon crystals and other semiconductor materials. The ultimate goal of the IC industry is to have the highest quality substrates, which are free of point defect, impurities and micro defect clusters. The scale up of silicon wafer size from 200 mm to 300 mm requires large crucible size and more heat power. Transport phenomena in crystal growth processes are quite complex due to melt and gas flows that may be oscillatory and/or turbulent, coupled convection and radiation, impurities and dopant distributions, unsteady kinetics of the growth process, melt crystal interface dynamics, free surface and meniscus, stoichiometry in the case of compound materials. A global model has been developed to simulate the temperature distribution and melt flow in an 8-inch system. The present program features the fluid convection, magnetohydrodynamics, and radiation models. A multi-zone method is used to divide the Cz system into different zones, e.g., the melt, the crystal and the hot zone. For calculation of temperature distribution, the whole system inside the stainless chamber is considered. For the convective flow, only the melt is considered. The widely used zonal method divides the surface of the radiation enclosure into a number of zones, which has a uniform distribution of temperature, radiative properties and composition. The integro-differential equations for the radiative heat transfer are solved using the matrix inversion technique. The zonal method for radiative heat transfer is used in the growth chamber, which is confined by crystal surface, melt surface, heat shield, and pull chamber. Free surface and crystal/melt interface are tracked using adaptive grid generation. The competition between the thermocapillary convection induced by non-uniform temperature distributions on the free surface and the forced convection by the rotation of the crystal determines the interface shape, dopant distribution, and striation pattern. The temperature gradients on the free surface are influenced by the effects of the thermocapillary force on the free surface and the rotation of the crystal and the crucible.

Short-Term Nonlinear Response of Tension Leg Platform in Random Sea Waves

Most of the existing mathematical models for analyzing the dynamic response of TLP are based on explicit or implicit assumptions that motions (translations and rotations) are small magnitude. However, when TLP works in severe adverse conditions, the a priori assumption on small displacements may be inadequate. In such situation, the motions should be regarded as finite magnitude. This paper will study stochastic nonlinear dynamic responses of TLP with finite displacements in random waves. The nonlinearities considered are: large amplitude motions, coupling the six degrees-of-freedom, instantaneous position, instantaneous wet surface, free surface effects and viscous drag force. The nonlinear dynamic responses are calculated by using numerical integration procedure in the time domain. After the time histories of the dynamic responses are obtained, we carry out cycle counting of the stress histories of the tethers with rain-flow counting method to get the stress range distribution.

Floating zone convection

The microgravity research, as a branch of the advanced sciences and a spe- cialized field of high technology, has been made in China since the late 1980's. The research group investigating microgravity fluid physics consisted of our col- leagues and the authors in the Institute of Mechanics of the Chinese Academy of Sciences (CAS), and we pay special attention to the floating zone convection as our first research priority. Now, the research group has expanded and is a part of the National Microgravity Laboratory of the CAS, and the research fields have been extended to include more subjects related to microgravity science. Howev- er, the floating zone convection is still an important topic that greatly holds our research interests.

Rarefied gas dynamics: fundamentals, simulations and micro flows

谈镐生文集

本文集收入了著名力学家、应用数学家、中国科学院院士谈镐生先生在流体力学、稀薄气体动力学和应用数学研究领域的论文和研究报告26篇，谈镐生先生倡导和支持力学基础研究的有关文章和论述21篇，谈先生的学术活动和生活图片多幅，以及谈镐生先生生平，最后附有谈先生生平年表。

 

 

Boiling heat transfer enhancement by using micro-pin-finned surface for electronics cooling

For efficiently cooling electronic components with high heat flux, experiments were conducted to study the flow boiling heat transfer performance of FC-72 over square silicon chips with the dimensions of 10 × 10 × 0.5 mm3. Four kinds of micro-pin-fins with the dimensions of 30 × 60, 30 × 120, 50 × 60, 50 × 120 μm2 (thickness, t × height, h) were fabricated on the chip surfaces by the dry etching technique for enhancing boiling heat transfer. A smooth surface was also tested for comparison. The experiments were made at three different fluid velocities (0.5, 1 and 2 m/s) and three different liquid subcoolings (15, 25 and 35 K). The results were compared with the previous published data of pool boiling. All micro-pin-fined surfaces show a considerable heat transfer enhancement compared with a smooth surface. Flow boiling can remarkably decrease wall superheat compared with pool boiling. At the velocities lower than 1 m/s, the micro-pin-finned surfaces show a sharp increase in heat flux with increasing wall superheat. For all surfaces, the maximum allowable heat flux, qmax, for the normal operation of LSI chips increases with fluid velocity and subcooling. For all micro-pin-finned surfaces, the wall temperature at the critical heat flux (CHF) is less than the upper limit for the reliable operation of LSI chips, 85◦C. The largest value of qmax can reach nearly 148 W/cm2 for micro-pin-finned chips with the fin height of 120 μm at the fluid velocity of 2 m/s and the liquid subcooling of 35 K. The perspectives for the boiling heat transfer experiment of the prospective micro-pin-finned sur- faces, which has been planned to be made in the Drop Tower Beijing/NMLC in the future, are also presented.

Effect of interfacial heat transfer on the onset of oscillatory convection in liquid bridge

In present study, effect of interfacial heat transfer with ambient gas on the onset of oscillatory convection in a liquid bridge of large Prandtl number on the ground is systematically investigated by the method of linear stability analyses. With both the constant and linear ambient air temperature distributions, the numerical results show that the interfacial heat transfer modifies the free-surface temperature distribution directly and then induces a steeper temperature gradient on the middle part of the free surface, which may destabilize the convection. On the other hand, the interfacial heat transfer restrains the temperature disturbances on the free surface, which may stabilize the convection. The two coupling effects result in a complex dependence of the stability property on the Biot number. Effects of melt free-surface deformation on the critical conditions of the oscillatory convection were also investigated. Moreover, to better understand the mechanism of the instabilities, rates of kinetic energy change and "thermal" energy change of the critical disturbances were investigated (C) 2009 Elsevier Ltd. All rights reserved.

Turbulent buoyant jets into stratified or flowing ambient fluids

Theoretical and experimental studies were made on two classes of buoyant jet problems, namely:

1) an inclined, round buoyant yet in a stagnant environment with linear density-stratification;

2) a round buoyant jet in a uniform cross stream of homogenous density.

Using the integral technique of analysis, assuming similarity, predictions can be made for jet trajectory, widths, and dilution ratios, in a density-stratified or flowing environment. Such information is of great importance in the design of disposal systems for sewage effluent into the ocean or waste gases into the atmosphere.

The present study of a buoyant jet in a stagnant environment has extended the Morton type of analysis to cover the effect of the initial angle of discharge. Numerical solutions have been presented for a range of initial conditions. Laboratory experiments were conducted for photographic observations of the trajectories of dyed jets. In general the observed jet forms agreed well with the calculated trajectories and nominal half widths when the value of the entrainment coefficient was taken to be α = 0.082, as previously suggested by Morton.

The problem of a buoyant jet in a uniform cross stream was analyzed by assuming an entrainment mechanism based upon the vector difference between the characteristic jet velocity and the ambient velocity. The effect of the unbalanced pressure field on the sides of the jet flow was approximated by a gross drag term. Laboratory flume experiments with sinking jets which are directly analogous to buoyant jets were performed. Salt solutions were injected into fresh water at the free surface in a flume. The jet trajectories, dilution ratios and jet half widths were determined by conductivity measurements. The entrainment coefficient, α, and drag coefficient, C_d, were found from the observed jet trajectories and dilution ratios. In the ten cases studied where jet Froude number ranged from 10 to 80 and velocity ratio (jet: current) K from 4 to 16, α varied from 0.4 to 0.5 and C_d from 1.7 to 0.1. The jet mixing motion for distance within 250D was found to be dominated by the self-generated turbulence, rather than the free-stream turbulence. Similarity of concentration profiles has also been discussed.

The unsteady forces on slender delta wing hydrofoils oscillating in heave

The investigations described herein are both experimental and theoretical. An experimental technique is described by which the models tested could be oscillated sinusoidally in heave. The apparatus used to gather the unsteady lift, drag and pitching moment data is also described.

The models tested were two flat delta wings with apex angles of 15° and 30° and they had sharp leading edges to insure flow separation. The models were fabricated from 0.25 inch aluminum plate and were approximately one foot in length.

Three distinct types of flow were investigated: 1) fully wetted, 2) ventilated and 3) planing. The experimental data are compared with existing theories for steady motions in the case of fully wetted delta wings. Ventilation measurements, made only for the 30° model at 20° angle of attack, of lift and drag are presented.

A correction of the theory proposed by M.P. Tulin for high speed planing of slender bodies is presented and it is extended to unsteady motions. This is compared to the experimental measurements made at 6° and 12° angle of attack for the two models previously described.

This is the first extensive measurement of unsteady drag for any shape wing, the first measurement of unsteady planing forces, the first quantitative documentation of unstable oscillations near a free surface, and the first measurements of the unsteady forces on ventilated delta wings. The results of these investigations, both theoretical and experimental, are discussed and further investigations suggested.

Numerical simulations of unsteady wet steam flows with non-equilibrium condensation in the nozzle and the steam turbine

Numerical techniques for non-equilibrium condensing flows are presented. Conservation equations for homogeneous gas-liquid two-phase compressible flows are solved by using a finite volume method based on an approximate Riemann solver. The phase change consists of the homogeneous nucleation and growth of existing droplets. Nucleation is computed with the classical Volmer-Frenkel model, corrected for the influence of the droplet temperature being higher than the steam temperature due to latent heat release. For droplet growth, two types of heat transfer model between droplets and the surrounding steam are used: a free molecular flow model and a semi-empirical two-layer model which is deemed to be valid over a wide range of Knudsen number. The computed pressure distribution and Sauter mean droplet diameters in a convergent-divergent (Laval) nozzle are compared with experimental data. Both droplet growth models capture qualitatively the pressure increases due to sudden heat release by the non-equilibrium condensation. However the agreement between computed and experimental pressure distributions is better for the two-layer model. The droplet diameter calculated by this model also agrees well with the experimental value, whereas that predicted by the free molecular model is too small. Condensing flows in a steam turbine cascade are calculated at different Mach numbers and inlet superheat conditions and are compared with experiments. Static pressure traverses downstream from the blade and pressure distributions on the blade surface agree well with experimental results in all cases. Once again, droplet diameters computed with the two-layer model give best agreement with the experiments. Droplet sizes are found to vary across the blade pitch due to the significant variation in expansion rate. Flow patterns including oblique shock waves and condensation-induced pressure increases are also presented and are similar to those shown in the experimental Schlieren photographs. Finally, calculations are presented for periodically unsteady condensing flows in a low expansion rate, convergent-divergent (Laval) nozzle. Depending on the inlet stagnation subcooling, two types of self-excited oscillations appear: a symmetric mode at lower inlet subcooling and an asymmetric mode at higher subcooling. Plots of oscillation frequency versus inlet sub-cooling exhibit a hysteresis loop, in accord with observations made by other researchers for moist air flow. Copyright © 2006 by ASME.

Shallow water equations: Magic and limitations

D Liang from Cambridge University explains the shallow water equations and their applications to the dam-break and other steep-fronted flow modeling. They assume that the horizontal scale of the flow is much greater than the vertical scale, which means the flow is restricted within a thin layer, thus the vertical momentum is insignificant and the pressure distribution is hydrostatic. The left hand sides of the two momentum equations represent the acceleration of the fluid particle in the horizontal plane. If the fluid acceleration is ignored, then the two momentum equations are simplified into the so-called diffusion wave equations. In contrast to the SWEs approach, it is much less convenient to model floods with the Navier-Stokes equations. In conventional computational fluid dynamics (CFD), cumbersome treatments are needed to accurately capture the shape of the free surface. The SWEs are derived using the assumptions of small vertical velocity component, smooth water surface, gradual variation and hydrostatic pressure distribution.

Hole nucleation and growth induced by crystallization and microphase separation of thin semicrystalline diblock copolymer films

We have investigated the hole nucleation and growth induced by crystallization of thin crystalline-coil diblock copolymer films. Semicrystalline rodlike assemblies from neutral/selective binary solvent are used as seeds to nucleate crystallization at temperatures above the glass transition temperature (T-g) but below melting point (T-m). The crystallization of nanorods drives neighboring copolymer chains to diffuse into the growing nanorods. Depletion of copolymer chains yields hole nucleation and growth at the edge of the nanorods. Simultaneously, the polymer chains unassociated into the nanorods were oriented by induction from the free surface and the substrate, leading to limitation of the hole depth to the lamellar spacing, similar to20 nm. The holes, as well as the nanorods, grow as t(alpha), where t is the annealing time and a crossover in the exponent a. is found. The orientation and stretching of the copolymer chains by the surface and interface are believed to accelerate the crystallization, and in turn, the latter accelerates the growth rate of the holes. At T > T-m, the grains melt and the copolymer chains relax and flow into the first layer of the film.