Two-way coupling model for shock-induced laminar boundary-layer flows of a dusty gas

The present paper describes a numerical two-way coupling model for shock-induced laminar boundary-layer flows of a dust-laden gas and studies the transverse migration of fine particles under the action of Saffman lift force. The governing equations are formulated in the dilute two-phase continuum framework with consideration of the finiteness of the particle Reynolds and Knudsen numbers. The full Lagrangian method is explored for calculating the dispersed-phase flow fields (including the number density of particles) in the regions of intersecting particle trajectories. The computation results show a significant reaction of the particles on the two-phase boundary-layer structure when the mass loading ratio of particles takes finite values.

Mathematical Modeling of Near-Wall Flows of Two-Phase Mixture with Evaporating Droplets

In the framework of the two-continuum approach, using the matched asymptotic expansion method, the equations of a laminar boundary layer in mist flows with evaporating droplets were derived and solved. The similarity criteria controlling the mist flows were determined. For the flow along a curvilinear surface, the forms of the boundary layer equations differ from the regimes of presence and absence of the droplet inertia deposition. The numerical results were presented for the vapor-droplet boundary layer in the neighborhood of a stagnation point of a hot blunt body. It is demonstrated that, due to evaporation, a droplet-free region develops near the wall inside the boundary layer. On the upper edge of this region, the droplet radius tends to zero and the droplet number density becomes much higher than that in the free stream. The combined effect of the droplet evaporation and accumulation results in a significant enhancement of the heat transfer on the surface even for small mass concentration of the droplets in the free stream. 在双连续介质理论框架下,采用匹配渐进展开方法导出并求解了具有蒸发液滴的汽雾流中层流边界层方程,给出了控制汽雾流的相似判据。对于沿曲面的流动,边界层方程的形式取决于是否存在液滴的惯性沉积。给出了热钝体验点附近蒸汽。液滴边界层的数值计算结果。它们表明：由于蒸发,在边界层内近壁处形成了一个无液滴区域；在该区上边界处,液滴半径趋于零而液滴数密度急剧增高。液滴蒸发及聚集的联合效应造成了表面热流的显著增加,甚至在自由来流中液滴质量浓度很低时此效应依然存在。

An improvement on the theoretical model for the large-scale coherent structures in the outer region of a turbulent boundary layer

According to the experimental results, there exist large-scale coherent structures in the outer region of a turbulent boundary layer, which have been studied by many authors.As experimental results, Antonia (1990) showed the phase- aver aged streamlines and isovorticity lines of the large-scale coherent structures in a turbulent boundary layer for different Reynolds numbers. Based on the hydrodynamic stability theory, the 2-D theoretical model for the large-scale structures was proposed by Luo and Zhou, in which the eddy viscosity was defined as a complex function of the position in the normal direction. The theoretical results showed in ref. were in agreement with those in ref. However, there were two problems in the results. One is that in the experimental results, there were divergent focuses between two saddle points in the streamlines, but in the theoretical results, there were centers. The other is that the stretched parts of the isovorticity lines appear at the location of centers in the theoretical results, while in the experimental results they located somewhere between the focuses and saddle points. The reason is that the computations were based on a 2-D model.

Compressible laminar boundary-layer flows of a dusty gas over a semi-infinite flat-plate

The compressible laminar boundary-layer flows of a dilute gas-particle mixture over a semi-infinite flat plate are investigated analytically. The governing equations are presented in a general form where more reasonable relations for the two-phase interaction and the gas viscosity are included. The detailed flow structures of the gas and particle phases are given in three distinct regions : the large-slip region near the leading edge, the moderate-slip region and the small-slip region far downstream. The asymptotic solutions for the two limiting regions are obtained by using a seriesexpansion method. The finite-difference solutions along the whole length of the plate are obtained by using implicit four-point and six-point schemes. The results from these two methods are compared and very good agreement is achieved. The characteristic quantities of the boundary layer are calculated and the effects on the flow produced by the particles are discussed. It is found that in the case of laminar boundary-layer flows, the skin friction and wall heat-transfer are higher and the displacement thickness is lower than in the pure-gas case alone. The results indicate that the Stokes-interaction relation is reasonable qualitatively but not correct quantitatively and a relevant non-Stokes relation of the interaction between the two phases should be specified when the particle Reynolds number is higher than unity.

Experimental study on two-phase gas-liquid flow patterns at normal and reduced gravity conditions

Experimental studies have been performed for horizontal two-phase air-water flows at normal and reduced gravity conditions in a square cross-section channel. The experiments at reduced gravity are conducted on board the Russian IL-76 reduced gravity airplane. Four flow patterns, namely bubble, slug, slug-annular transition and annular flows, are observed depending on the liquid and gas superficial velocities at both conditions. Semi-theoretical Weber number model is developed to include the shape influence on the slug-annular transition. It is shown that its prediction is in reasonable agreement with the experimental slug-annular transition under both conditions. For the case of two-phase gas-liquid flow with large value of the Froude number, the drift-flux model can predict well the observed boundary between bubble and slug flows.

Two-Phase Flow In Fuel Cells In Short-Term Microgravity Condition

A visual observation of liquid-gas two-phase flow in anode channels of a direct methanol proton exchange membrane fuel cells in microgravity has been carried out in a drop tower. The anode flow bed consisted of 2 manifolds and 11 parallel straight channels. The length, width and depth of single channel with rectangular cross section was 48.0 mm, 2.5 mm and 2.0 mm, respectively. The experimental results indicated that the size of bubbles in microgravity condition is bigger than that in normal gravity. The longer the time, the bigger the bubbles. The velocity of bubbles rising is slower than that in normal gravity because buoyancy lift is very weak in microgravity. The flow pattern in anode channels could change from bubbly flow in normal gravity to slug flow in microgravity. The gas slugs blocked supply of reactants from channels to anode catalyst layer through gas diffusion layer. When the weakened mass transfer causes concentration polarization, the output performance of fuel cells declines.

Modeling, numerical methods, and simulation for particle-fluid two-phase flow problems

In this paper, we study the issues of modeling, numerical methods, and simulation with comparison to experimental data for the particle-fluid two-phase flow problem involving a solid-liquid mixed medium. The physical situation being considered is a pulsed liquid fluidized bed. The mathematical model is based on the assumption of one-dimensional flows, incompressible in both particle and fluid phases, equal particle diameters, and the wall friction force on both phases being ignored. The model consists of a set of coupled differential equations describing the conservation of mass and momentum in both phases with coupling and interaction between the two phases. We demonstrate conditions under which the system is either mathematically well posed or ill posed. We consider the general model with additional physical viscosities and/or additional virtual mass forces, both of which stabilize the system. Two numerical methods, one of them is first-order accurate and the other fifth-order accurate, are used to solve the models. A change of variable technique effectively handles the changing domain and boundary conditions. The numerical methods are demonstrated to be stable and convergent through careful numerical experiments. Simulation results for realistic pulsed liquid fluidized bed are provided and compared with experimental data. (C) 2004 Elsevier Ltd. All rights reserved.

Study of the shock motion in a hypersonic shock system turbulent boundary-layer interaction

Wall pressure fluctuations and surface heat transfer signals have been measured in the hypersonic turbulent boundary layer over a number of compression-corner models. The distributions of the separation shock oscillation frequencies and periods have been calculated using a conditional sampling algorithm. In all cases the oscillation frequency distributions are of broad band, but the most probable frequencies are low. The VITA method is used for deducing large scale disturbances at the wall in the incoming boundary layer and the separated flow region. The results at present showed the existence of coherent structures in the two regions. The zero-cross frequencies of the large scale structures in the two regions are of the same order as that of the separation shock oscillation. The average amplitude of the large scale structures in the separated region is much higher than that in the incoming boundary layer. The length scale of the separation shock motion region is found to increase with the disturbance strength. The results show that the shock oscillation is of inherent nature in the shock wave/turbulent boundary layer interaction with separation. The shock oscillation is considered to be the consequence of the coherent structures in the separated region.

Air flow over a mountain and the convective boundary-layer

In this paper, the analytical model coupling the convective boundary layer (CBL) with the free atmosphere developed by Qi and Fu (1992) is improved. And by this improved model, the interaction between airflow over a mountain and the CBL is further discussed. The conclusions demonstrate: (1) The perturbation potential temperatures in the free atmosphere can counteract the effect of orographic thermal forcing through entraining and mixing in the CBL. If u(M)BAR > u(F)BAR, the feedback of the perturbation potential temperatures in the free atmosphere is more important than orographic thermal forcing, which promotes the effect of interfacial waves. If u(M)BAR < u(F)BAR, orographic thermal forcing is more important, which makes the interfacial height and the topographic height identical in phase, and the horizontal speeds are a maximum at the top of the mountain. (2) The internal gravity waves propagating vertically in the free atmosphere cause a strong downslope wind to become established above the lee slope in the CBL and result in the hydraulic jump at the top of the CBL. (3) With the CBL deepening, the interfacial gravity waves induced by the potential temperature jump at the top of the CBL cause the airflow in the CBL to be subcritical.

Two-phase Flow in Fuel Cells in Short-term Microgravity Condition

A visual observation of liquid-gas two-phase flow in anode channels of a direct methanol proton exchange membrane fuel cells in microgravity has been carried out in a drop tower. The anode flow bed consisted of 2 manifolds and 11 parallel straight channels. The length, width and depth of single channel with rectangular cross section was 48.0 mm, 2.5 mm and 2.0 mm, respectively. The experimental results indicated that the size of bubbles in microgravity condition is bigger than that in normal gravity. The longer the time, the bigger the bubbles. The velocity of bubbles rising is slower than that in normal gravity because buoyancy lift is very weak in microgravity. The flow pattern in anode channels could change from bubbly flow in normal gravity to slug flow in microgravity. The gas slugs blocked supply of reactants from channels to anode catalyst layer through gas diffusion layer. When the weakened mass transfer causes concentration polarization, the output performance of fuel cells declines.

Direct numerical simulation of hypersonic boundary layer transition over a blunt cone with a small angle of attack

The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [Li et al., AIAA J. 46, 2899(2008)], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 0° ≤ θ ≤ 30° (θ = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ θ ≤ 30°, which leads to the delayed transition location. Very low frequency waves VLFWs� are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [ Li et al., AIAA J. 46, 2899 (2008) ], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 20° ≤ θ ≤ 30° (θ = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ θ ≤ 30°, which leads to the delayed transition location. Very low frequency waves (VLFWs) are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.

Synthesis of Ruthenium Dioxide Nanoparticles by a Two-Phase Route and Their Electrochemical Properties

Dissolvable, size- and shape-controlled ruthenium dioxide nanoparticles are successfully achieved through a two-phase route. The influence of reaction time, temperature, and monomer concentration and the nature of capping agents on the morphologies of nanoparticles are studied through transmission electron microscopy (TEM). A possible mechanism for the formation and growth of nanoparticles is also involved. X-ray powder diffraction (XRD) confirms the amorphous structure for as-prepared ruthenium dioxide nanoparticles. Samples are immobilized by simple dip-coating on a current collector, and the cyclic voltammetry measurement is utilized to investigate their electrochemical properties. The specific capacitance of one sample can teach as high as 840 F g(-1), which reveals the promising application potential to electrochemical capacitors.

Pressure Drop of Bubbly Two-Phase Flow in a Square Channel at Reduced Gravity

A new set of experimental pressure drop data, collected aboard the Russian IL-76MDK, is reported for bubbly airwater two-phase flow in a square channel with a cross-sectional area of 12x 12mm(2). The present data are compared to several frequently used empirical models, e.g. homogeneous model, Lockhart-Martinelli-Chisholm correlation and Friedel's model. It is shown that the predictions of the models mentioned above are generally not satisfied. A new homogeneous model is developed based on the analysis of the characteristics of bubbly two-phase flow at reduced gravity. It is tested with the present data and other data collected by other researchers in circular pipes. Some questions related to the present model are also discussed. (C) 2002 COSPAR. Published by Elsevier Science Ltd. All rights reserved.

Direct Numerical Simulation of a Spatially Evolving Supersonic Turbulent Boundary Layer at Ma=6

Direct numerical simulation is carried out for a spatially evolving supersonic turbulent boundary layer at free-stream Mach number 6. To overcome numerical instability, the seventh-order WENO scheme is used for the convection terms of Navier-Stokes equations, and fine mesh is adopted to minimize numerical dissipation. Compressibilty effects on the near-wall turbulent kinetic energy budget are studied. The cross-stream extended self-similarity and scaling exponents including the near-wall region are studied. In high Mach number flows, the coherence vortex structures are arranged to be smoother and streamwised, and the hair-pin vortices are less likely to occur.