Information Preservation Modeling of Rayleigh-Benard Transition from Thermal Conduction to Convection

Onset and evolution of the Rayleigh-Benard (R-B) convection are investigated using the Information Preservation (IP) method. The information velocity and temperature are updated using the Octant Flux Splitting (OFS) model developed by Masters & Ye based on the Maxwell transport equation suggested by Sun & Boyd. Statistical noise inherent in particle approaches such as the direct simulation Monte Carlo (DSMC) method is effectively reduced by the IP method, and therefore the evolutions from an initial quiescent fluid to a final steady state are shown clearly. An interesting phenomenon is observed: when the Rayleigh number (Ra) exceeds its critical value, there exists an obvious incubation stage. During the incubation stage, the vortex structure clearly appears and evolves, whereas the Nusselt number (Nu) of the lower plate is close to unity. After the incubation stage, the vortex velocity and Nu rapidly increase, and the flow field quickly reaches a steady, convective state. A relation of Nu to Ra given by IP agrees with those given by DSMC, the classical theory and experimental data.

Two bifurcation processes for onset of oscillatory thermocapillary convection in a floating half zone

An overview on the onset of thermocapillary oscillatory convection in a floating half zone is provided, and it is a typical subject in the microgravity sciences related to the space materials science, especially the floating zone processing, and also to the microgravity fluid physics. The main interests are focused around the process for onset of oscillatory thermocapillary convection, which is known also as the bifurcation transition from quasi-steady convection to oscillatory convection. The onset of oscillation depends on a set of critical parameters, such as the Marangoni number, Prandtl number, geometrical parameters, and heat transfer parameters. Recent studies show that, there exists the bifurcation transition from steady and axial symmetric convection to the steady and axial non-symmetric convection before the onset of oscillation in cases of small Prandtl number fluids and in cases of larger Prandtl number fluids of fat liquid bridge with small aspect ratio. The transition process is a strong non-linear process because the velocity deviation has the same order of magnitude as that of an average flow after the onset of oscillation, and unsteady 3-D numerical simulation is suitable to do in depth analysis on strong non-linear process, and leads generally to a better comparison with the experimental results.

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.

Scale-similarity model for Lagrangian velocity correlations in isotropic and stationary turbulence

A scale-similarity model for Lagrangian two-point, two-time velocity correlations LVCs in isotropic turbulence is developed from the Kolmogorov similarity hypothesis. It is a second approximation to the isocontours of LVCs, while the Smith-Hay model is only a first approximation. This model expresses the LVC by its space correlation and a dispersion velocity. We derive the analytical expression for the dispersion velocity from the Navier-Stokes equations using the quasinormality assumption. The dispersion velocity is dependent on enstrophy spectra and shown to be smaller than the sweeping velocity for the Eulerian velocity correlation. Therefore, the Lagrangian decorrelation process is slower than the Eulerian decorrelation process. The data from direct numerical simulation of isotropic turbulence support the scale-similarity model: the LVCs for different space separations collapse into a universal form when plotted against the separation axis defined by the model.

Measurement of velocity profiles in a rectangular microchannel with aspect ratio α = 0.35

In this work, we measured 14 horizontal velocity profiles along the vertical direction of a rectangular microchannel with aspect ratio alpha = h/w = 0.35 (h is the height of the channel and w is the width of the channel) using microPIV at Re = 1.8 and 3.6. The experimental velocity profiles are compared with the full 3D theoretical solution, and also with a Poiseuille parabolic profile. It is shown that the experimental velocity profiles in the horizontal and vertical planes are in agreement with the theoretical profiles, except for the planes close to the wall. The discrepancies between the experimental data and 3D theoretical results in the center vertical plane are less than 3.6%. But the deviations between experimental data and Poiseuille's results approaches 5%. It indicates that 2D Poiseuille profile is no longer a perfect theoretical approximation since a = 0.35. The experiments also reveal that, very near the hydrophilic wall (z = 0.5-1 mu m), the measured velocities are significantly larger than the theoretical velocity based on the no-slip assumption. A proper discussion on some physical effects influencing the near wall velocity measurement is given.

Subgrid scale fluid velocity timescales seen by inertial particles in large-eddy simulation of particle-laden turbulence

Large-eddy simulation (LES) has emerged as a promising tool for simulating turbulent flows in general and, in recent years,has also been applied to the particle-laden turbulence with some success (Kassinos et al., 2007). The motion of inertial particles is much more complicated than fluid elements, and therefore, LES of turbulent flow laden with inertial particles encounters new challenges. In the conventional LES, only large-scale eddies are explicitly resolved and the effects of unresolved, small or subgrid scale (SGS) eddies on the large-scale eddies are modeled. The SGS turbulent flow field is not available. The effects of SGS turbulent velocity field on particle motion have been studied by Wang and Squires (1996), Armenio et al. (1999), Yamamoto et al. (2001), Shotorban and Mashayek (2006a,b), Fede and Simonin (2006), Berrouk et al. (2007), Bini and Jones (2008), and Pozorski and Apte (2009), amongst others. One contemporary method to include the effects of SGS eddies on inertial particle motions is to introduce a stochastic differential equation (SDE), that is, a Langevin stochastic equation to model the SGS fluid velocity seen by inertial particles (Fede et al., 2006; Shotorban and Mashayek, 2006a; Shotorban and Mashayek, 2006b; Berrouk et al., 2007; Bini and Jones, 2008; Pozorski and Apte, 2009).However, the accuracy of such a Langevin equation model depends primarily on the prescription of the SGS fluid velocity autocorrelation time seen by an inertial particle or the inertial particle–SGS eddy interaction timescale (denoted by $\delt T_{Lp}$ and a second model constant in the diffusion term which controls the intensity of the random force received by an inertial particle (denoted by C_0, see Eq. (7)). From the theoretical point of view, dTLp differs significantly from the Lagrangian fluid velocity correlation time (Reeks, 1977; Wang and Stock, 1993), and this carries the essential nonlinearity in the statistical modeling of particle motion. dTLp and C0 may depend on the filter width and particle Stokes number even for a given turbulent flow. In previous studies, dTLp is modeled either by the fluid SGS Lagrangian timescale (Fede et al., 2006; Shotorban and Mashayek, 2006b; Pozorski and Apte, 2009; Bini and Jones, 2008) or by a simple extension of the timescale obtained from the full flow field (Berrouk et al., 2007). In this work, we shall study the subtle and on-monotonic dependence of $\delt T_{Lp}$ on the filter width and particle Stokes number using a flow field obtained from Direct Numerical Simulation (DNS). We then propose an empirical closure model for $\delta T_{Lp}$. Finally, the model is validated against LES of particle-laden turbulence in predicting single-particle statistics such as particle kinetic energy. As a first step, we consider the particle motion under the one-way coupling assumption in isotropic turbulent flow and neglect the gravitational settling effect. The one-way coupling assumption is only valid for low particle mass loading.

The near wall velocity measurements in microchannels with different diameter particles by microPIV/PTV. 8th International symposium on PIV

In near wall measurements with microPIV/PTV, whether seeding particles can be effectively used to detect local fluid velocity is a crucial problem. This talk presents our recent measurements in microchannels [1][2]. Based on measured velocity profiles with 200nm and 50nm in pure water, we found that the measured velocity profiles are agreed with the theoretical values in the middle of channel, but large deviations between measured data and theoretical prediction appear close to wall (0.25mm velocity deviation appeared in shear flows.

Numerical Study on Critical Wedge Angle of Cellular Detonation Reflections

The critical wedge angle (CWA) for the transition from regular reflection (RR) to Mach reflection (MR) of a cellular detonation wave is studied numerically by an improved space-time conservation element and solution element method together with a two-step chemical reaction model. The accuracy of that numerical way is verified by simulating cellular detonation reflections at a 19.3∘ wedge. The planar and cellular detonation reflections over 45∘–55∘ wedges are also simulated. When the cellular detonation wave is over a 50∘ wedge, numerical results show a new phenomenon that RR and MR occur alternately. The transition process between RR and MR is investigated with the local pressure contours. Numerical analysis shows that the cellular structure is the essential reason for the new phenomenon and the CWA of detonation reflection is not a certain angle but an angle range.

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.

Experimental Study Of Vortex-Induced Vibrations Of A Cylinder Near A Rigid Plane Boundary In Steady Flow

In this study, the vortex-induced vibrations of a cylinder near a rigid plane boundary in a steady flow are studied experimentally. The phenomenon of vortex-induced vibrations of the cylinder near the rigid plane boundary is reproduced in the flume. The vortex shedding frequency and mode are also measured by the methods of hot film velocimeter and hydrogen bubbles. A parametric study is carried out to investigate the influences of reduced velocity, gap-to-diameter ratio, stability parameter and mass ratio on the amplitude and frequency responses of the cylinder. Experimental results indicate: (1) the Strouhal number (St) is around 0.2 for the stationary cylinder near a plane boundary in the sub-critical flow regime; (2) with increasing gap-to-diameter ratio (e (0)/D), the amplitude ratio (A/D) gets larger but frequency ratio (f/f (n) ) has a slight variation for the case of larger values of e (0)/D (e (0)/D > 0.66 in this study); (3) there is a clear difference of amplitude and frequency responses of the cylinder between the larger gap-to-diameter ratios (e (0)/D > 0.66) and the smaller ones (e (0)/D < 0.3); (4) the vibration of the cylinder is easier to occur and the range of vibration in terms of V (r) number becomes more extensive with decrease of the stability parameter, but the frequency response is affected slightly by the stability parameter; (5) with decreasing mass ratio, the width of the lock-in ranges in terms of V (r) and the frequency ratio (f/f (n) ) become larger.

Group velocity of probe light pulse in an open lambda system

The group velocity of the probe light pulse (GVPLP) propagating through an open Lambda-type atomic system with a spontaneously generated coherence is investigated when the weak probe and strong driving light fields have different frequencies. It is found that adjusting the detuning or Rabi frequency of the probe light field can realize switching of the GVPLP from subluminal to superluminal. Changing the relative phase between the probe and driving light. elds or atomic exit and injection rates can lead to GVPLP varying in a wider range, but cannot induce transformation of the property of the GVPLP. The absolute value of the GVPLP always increases with Rabi frequency of the driving light field increasing. For subluminal and superluminal propagation, the system always exhibits the probe absorption, and GVPLP is mainly determined by the slope of the steep dispersion.

Controllable group velocity of the probe light in a Lambda-type system with two fold levels

The group velocities of the probe laser field are studied in a A-type system where one lower state has two fold levels coupled by a control field. It is found that the interaction of double dark states leads to controllable group velocity of the probe field in this system. It can be easily realized, due to the interacting double dark resonances, that one of the group velocities at transparency positions is much slower than the other by tuning the control field to be off resonance. In particular, when the control field is on resonance. we can obtain two equal slow group velocities with a broader EIT width, which provides potential applications in quantum storage and retrieval of light. (c) 2005 Elsevier B.V. All rights reserved.

Influence of Coulomb potential for photoionization of H atoms in an elliptically polarized laser field: Velocity gauge versus length gauge

We theoretically study the influence of Coulomb potential for photoionization of hydrogen atoms in an intense laser field with elliptical polarization. The total ionization rates, photoelectron energy spectra, and photoelectron angular distributions are calculated with the Coulomb-Volkov wave functions in the velocity gauge and compared with those calculated in the length gauge as well as those calculated with the Volkov wave functions. By comparing the results obtained by the Coulomb-Volkov and Volkov wave functions, we find that for linear polarization the influence of Coulomb potential is obvious for low-energy photoelectrons, and as the photoelectron energy and/or the laser intensity increase, its influence becomes smaller. This trend, however, is not so clear for the case of elliptical polarization. We also find that the twofold symmetry in the photoelectron angular distributions for elliptical polarization is caused by the cooperation of Coulomb potential and interference of multiple transition channels. About the gauge issue, we show that the difference in the photoelectron angular distributions obtained by the velocity and length gauges becomes rather obvious for elliptical polarization, while the difference is generally smaller for linear polarization.

Study of noncollinear quasi-phase matching optical parametric amplification with group-velocity matching based on periodically poled crystals

The properties of noncollinear optical parametric amplification (NOPA) based on quasi-phase matching of periodically poled crystals are investigated, under the condition that the group velocity matching (GVM) of the signal and idler pulses is satisfied. Our study focuses on the dependence of the gain spectrum upon the noncollinear angle, crystal temperature, and crystal angle with periodically poled KTiOPO4 (PPKTP), periodically poled LiNbO3 (PPLN), and periodically poled LiTaO3 (PPLT), and the NOPA gain properties of the three crystals are compared. Broad gain bandwidth exists above 85 nm at a signal wavelength of 800 nm with a 532 nm pump pulse, with proper noncollinear angle and grating period at a fixed temperature for GVM. Deviation from the group-velocity-matched noncollinear angle can be compensated by accurately tuning the crystal angle or temperature with a fixed grating period for phase matching. Moreover, there is a large capability of crystal angle tuning.