Experimental Study on Benard-Marangoni Convection by PIV and TCL

An experimental investigation of Benard-Marangoni convection has been performed in double immiscible liquid layers of rectangular configuration. The two kinds of liquid are 10cst silicon oil and FC-70 respectively. The velocity fields in the vertical cross-section are obtained by PIV. Flow patterns and/or temperature distributions on the horizontal interface are displayed by using thermal color liquid crystal (TLC), and the velocity distributions on the interface were also obtained with the help of the serial particle image of TCL. The evolution processes of convection are observed in the differential thickness ratio of two liquid layers, and the convection styles are discussed.

Experimental Study of Surface Deformation and Flow Pattern on Buoyant-Thermocapillary Convection

Free surface deformation is one of the most important physical phenomena in fluids with free surface. In the present paper, convection and surface deformation caused by thermocapillary effect in a rectangular cavity were investigated. In ground experiments, the convection was also affected by gravity. The cavity has a horizontal cross section of 52mm×42mm and the thikkness of the liquid layer is 4mm. Temperature difference between two sides of the liquid layer was increased gradually, and the flow in liquid layer will develop from steady to unstable convection. An optical diagnostic system consisting of a revised Michelson interferometer with image processor was developed to study fluid surface deformation in convection, and the displacements of free surface oscillation were determined. PIV technique was adopted to observe the evolution of flow pattern, and the velocity fields were obtained quantitatively. The present experiments demonstrate that surface deformation is quite distinct in buoyant-thermocapillary convection. in order to understand the mechanism of buoyant-thermocapillary convection, not only the hydrothermal wave instability but also the surface wave instability should be discussed.

Experimental Investigation Of Influence Of Interfacial Tension On Convection Of Two-Layer Immiscible Liquid

Benard-Marangoni convections of two-layer fluids heated from the bottom are investigated experimentally with a particle imagine velocimetry. The flows are visualized from the side, and various velocity fields near the onset of convection, such as three-layer vortex convective patterns, are observed when the depth ratio varies in a wide range. A new classification of the convective patterns is proposed with more detail than in previous studies. The analysis of the results indicates that the interface tension greatly influences the motion intensities of the bottom and top layers. The dimensionless wave number increases with the Bond number when the motion in the top layer is not more intense than that in the bottom layer, which agrees with the theoretical prediction.

Prediction of particle distribution in isotropic turbulence by large-eddy simulation

The recent application of large-eddy simulation (LES) to particle-laden turbulence requires that the LES with a subgrid scale (SGS) model could accurately predict particle distributions. Usually, a SGS particle model is used to recover the small-scale structures of velocity fields. In this study, we propose a rescaling technique to recover the effects of small-scale motions on the preferential concentration of inertial particles. The technique is used to simulate particle distribution in isotropic turbulence by LES and produce consistent results with direct numerical simulation (DNS). Key words: particle distribution, particle-laden turbulence, large-eddy simulation, subgrid scale model.

A kinematic subgrid scale model for large-eddy simulation of turbulence generated sound

In the hybrid approach of large-eddy simulation (LES) and Lighthill’s acoustic analogy for turbulence-generated sound, the turbulence source fields are obtained using an LES and the turbulence-generated sound at far fields is calculated from Lighthill’s acoustic analogy. As only the velocity fields at resolved scales are available from the LES, the Lighthill stress tensor, serving as a source term in Lighthill’s acoustic equation, has to be evaluated from the resolved velocity fields. As a result, the contribution from the unresolved velocity fields is missing in the conventional LES. The sound of missing scales is shown to be important and hence needs to be modeled. The present study proposes a kinematic subgrid-scale (SGS) model which recasts the unresolved velocity fields into Lighthill’s stress tensors. A kinematic simulation is used to construct the unresolved velocity fields with the imposed temporal statistics, which is consistent with the random sweeping hypothesis. The kinematic SGS model is used to calculate sound power spectra from isotropic turbulence and yields an improved result: the missing portion of the sound power spectra is approximately recovered in the LES.

Plasma-Arc Technology For The Thermal Treatment Of Chemical Wastes

Plasma-arc technology was developed to dispose of chemical wastes from a chemical plant by the Institute of Mechanics, Chinese Academy of Sciences (CAS-IMECH). A pilot plant system with this technology was constructed to destroy two types of chemical wastes. The system included shredding, mixing, and feeding subsystems, a plasma-arc reactor of 150 kW, an off-gas burning subsystem, and a scrubbing subsystem. The additives (CaO, SiO2, and Fe) were added into the reactor to form vitrified slag and capture the hazardous elements. The molten slag was quickly quenched to form an amorphous glassy structure. A direct current (DC) experimental facility of 30kW with plasma-arc technology was also set up to study the pyrolysis process in the laboratory, and the experimental results showed the cooling speed is the most important factor for good vitrified structure of the slag. According to previous tests, the destruction and removal efficiency (DRE) for these chemical wastes was more than 99.999%, and the polychlorinated biphenyls (PCBs) concentration in the solid residues was in the range of 1.28 to 12.9mg/kg, which is far below the Chinese national emission limit for the hazardous wastes. A simplified electromagneto model for numerical simulation was developed to predict the temperature and velocity fields. This model can make satisfactory maximum temperature and velocity distributions in the arc region, as well as the results by the magneto hydrodynamic approach.

Detailed properties of high redshift galaxies

Galaxies evolve throughout the history of the universe from the first star-forming sources, through gas-rich asymmetric structures with rapid star formation rates, to the massive symmetrical stellar systems observed at the present day. Determining the physical processes which drive galaxy formation and evolution is one of the most important questions in observational astrophysics. This thesis presents four projects aimed at improving our understanding of galaxy evolution from detailed measurements of star forming galaxies at high redshift.

We use resolved spectroscopy of gravitationally lensed z ≃ 2 - 3 star forming galaxies to measure their kinematic and star formation properties. The combination of lensing with adaptive optics yields physical resolution of ≃ 100 pc, sufficient to resolve giant Hii regions. We find that ~ 70 % of galaxies in our sample display ordered rotation with high local velocity dispersion indicating turbulent thick disks. The rotating galaxies are gravitationally unstable and are expected to fragment into giant clumps. The size and dynamical mass of giant Hii regions are in agreement with predictions for such clumps indicating that gravitational instability drives the rapid star formation. The remainder of our sample is comprised of ongoing major mergers. Merging galaxies display similar star formation rate, morphology, and local velocity dispersion as isolated sources, but their velocity fields are more chaotic with no coherent rotation.

We measure resolved metallicity in four lensed galaxies at z = 2.0 − 2.4 from optical emission line diagnostics. Three rotating galaxies display radial gradients with higher metallicity at smaller radii, while the fourth is undergoing a merger and has an inverted gradient with lower metallicity at the center. Strong gradients in the rotating galaxies indicate that they are growing inside-out with star formation fueled by accretion of metal-poor gas at large radii. By comparing measured gradients with an appropriate comparison sample at z = 0, we demonstrate that metallicity gradients in isolated galaxies must flatten at later times. The amount of size growth inferred by the gradients is in rough agreement with direct measurements of massive galaxies. We develop a chemical evolution model to interpret these data and conclude that metallicity gradients are established by a gradient in the outflow mass loading factor, combined with radial inflow of metal-enriched gas.

We present the first rest-frame optical spectroscopic survey of a large sample of low-luminosity galaxies at high redshift (L < L*, 1.5 < z < 3.5). This population dominates the star formation density of the universe at high redshifts, yet such galaxies are normally too faint to be studied spectroscopically. We take advantage of strong gravitational lensing magnification to compile observations for a sample of 29 galaxies using modest integration times with the Keck and Palomar telescopes. Balmer emission lines confirm that the sample has a median SFR ∼ 10 M_sun yr^−1 and extends to lower SFR than has been probed by other surveys at similar redshift. We derive the metallicity, dust extinction, SFR, ionization parameter, and dynamical mass from the spectroscopic data, providing the first accurate characterization of the star-forming environment in low-luminosity galaxies at high redshift. For the first time, we directly test the proposal that the relation between galaxy stellar mass, star formation rate, and gas phase metallicity does not evolve. We find lower gas phase metallicity in the high redshift galaxies than in local sources with equivalent stellar mass and star formation rate, arguing against a time-invariant relation. While our result is preliminary and may be biased by measurement errors, this represents an important first measurement that will be further constrained by ongoing analysis of the full data set and by future observations.

We present a study of composite rest-frame ultraviolet spectra of Lyman break galaxies at z = 4 and discuss implications for the distribution of neutral outflowing gas in the circumgalactic medium. In general we find similar spectroscopic trends to those found at z = 3 by earlier surveys. In particular, absorption lines which trace neutral gas are weaker in less evolved galaxies with lower stellar masses, smaller radii, lower luminosity, less dust, and stronger Lyα emission. Typical galaxies are thus expected to have stronger Lyα emission and weaker low-ionization absorption at earlier times, and we indeed find somewhat weaker low-ionization absorption at higher redshifts. In conjunction with earlier results, we argue that the reduced low-ionization absorption is likely caused by lower covering fraction and/or velocity range of outflowing neutral gas at earlier epochs. This result has important implications for the hypothesis that early galaxies were responsible for cosmic reionization. We additionally show that fine structure emission lines are sensitive to the spatial extent of neutral gas, and demonstrate that neutral gas is concentrated at smaller galactocentric radii in higher redshift galaxies.

The results of this thesis present a coherent picture of galaxy evolution at high redshifts 2 ≲ z ≲ 4. Roughly 1/3 of massive star forming galaxies at this period are undergoing major mergers, while the rest are growing inside-out with star formation occurring in gravitationally unstable thick disks. Star formation, stellar mass, and metallicity are limited by outflows which create a circumgalactic medium of metal-enriched material. We conclude by describing some remaining open questions and prospects for improving our understanding of galaxy evolution with future observations of gravitationally lensed galaxies.

Dynamic Stall on Vertical Axis Wind Turbines

In this study the dynamics of flow over the blades of vertical axis wind turbines was investigated using a simplified periodic motion to uncover the fundamental flow physics and provide insight into the design of more efficient turbines. Time-resolved, two-dimensional velocity measurements were made with particle image velocimetry on a wing undergoing pitching and surging motion to mimic the flow on a turbine blade in a non-rotating frame. Dynamic stall prior to maximum angle of attack and a leading edge vortex development were identified in the phase-averaged flow field and captured by a simple model with five modes, including the first two harmonics of the pitch/surge frequency identified using the dynamic mode decomposition. Analysis of these modes identified vortical structures corresponding to both frequencies that led the separation and reattachment processes, while their phase relationship determined the evolution of the flow.

Detailed analysis of the leading edge vortex found multiple regimes of vortex development coupled to the time-varying flow field on the airfoil. The vortex was shown to grow on the airfoil for four convection times, before shedding and causing dynamic stall in agreement with 'optimal' vortex formation theory. Vortex shedding from the trailing edge was identified from instantaneous velocity fields prior to separation. This shedding was found to be in agreement with classical Strouhal frequency scaling and was removed by phase averaging, which indicates that it is not exactly coupled to the phase of the airfoil motion.

The flow field over an airfoil undergoing solely pitch motion was shown to develop similarly to the pitch/surge motion; however, flow separation took place earlier, corresponding to the earlier formation of the leading edge vortex. A similar reduced-order model to the pitch/surge case was developed, with similar vortical structures leading separation and reattachment; however, the relative phase lead of the separation mode, corresponding to earlier separation, necessitated that a third frequency to be incorporated into the reattachment mode to provide a relative lag in reattachment.

Finally, the results are returned to the rotating frame and the effects of each flow phenomena on the turbine are estimated, suggesting kinematic criteria for the design of improved turbines.

Investigation of air-core vortex at hydraulic intakes

Hydrodynamic properties of the surface vortex have been investigated. Based on the Navier-Stokes equations, three sets of the new formulations for the tangential velocity distributions are derived, and verified against the experimental measurements in the literature. It is shown that one modification greatly improves the agreement with the experimental data. Physical model experiments were carried out to study the intake vortex related to the Xiluodu hydropower project. The velocity fields were measured using the Particle Tracking Velocimetry (PTV) technique. The proposed equation for tangential velocity distribution is applied to the Xiluodu project with the solid boundary being considered by the method of images. Good agreement has been observed between the formula prediction and the experimental observation. © 2010 Publishing House for Journal of Hydrodynamics.

Instabilities in the oscillatory flow of a complex fluid.

The dynamics of a fluid in a vertical tube, subjected to an oscillatory pressure gradient, is studied experimentally for both a Newtonian and a viscoelastic shear-thinning fluid. Particle image velocimetry is used to determine the two-dimensional velocity fields in the vertical plane of the tube axis, in a range of driving amplitudes from 0.8 to 2.5 mm and of driving frequencies from 2.0 to 11.5 Hz. The Newtonian fluid exhibits a laminar flow regime, independent of the axial position, in the whole range of drivings. For the complex fluid, instead, the parallel shear flow regime exhibited at low amplitudes [Torralba, Phys. Rev. E 72, 016308 (2005)] becomes unstable at higher drivings against the formation of symmetric vortices, equally spaced along the tube. At even higher drivings the vortex structure itself becomes unstable, and complex nonsymmetric structures develop. Given that inertial effects remain negligible even at the hardest drivings (Re < 10(-1)), it is the complex rheology of the fluid that is responsible for the instabilities observed. The system studied represents an interesting example of the development of shear-induced instabilities in nonlinear complex fluids in purely parallel shear flow.

Emerging modes of collective cell migration induced by geometrical constraints.

The role of geometrical confinement on collective cell migration has been recognized but has not been elucidated yet. Here, we show that the geometrical properties of the environment regulate the formation of collective cell migration patterns through cell-cell interactions. Using microfabrication techniques to allow epithelial cell sheets to migrate into strips whose width was varied from one up to several cell diameters, we identified the modes of collective migration in response to geometrical constraints. We observed that a decrease in the width of the strips is accompanied by an overall increase in the speed of the migrating cell sheet. Moreover, large-scale vortices over tens of cell lengths appeared in the wide strips whereas a contraction-elongation type of motion is observed in the narrow strips. Velocity fields and traction force signatures within the cellular population revealed migration modes with alternative pulling and/or pushing mechanisms that depend on extrinsic constraints. Force transmission through intercellular contacts plays a key role in this process because the disruption of cell-cell junctions abolishes directed collective migration and passive cell-cell adhesions tend to move the cells uniformly together independent of the geometry. Altogether, these findings not only demonstrate the existence of patterns of collective cell migration depending on external constraints but also provide a mechanical explanation for how large-scale interactions through cell-cell junctions can feed back to regulate the organization of migrating tissues.

Time-resolved particle image velocimetry within the nozzle of a drop-on-demand printhead

Time-resolved particle image velocimetry (PIV) has been performed inside the nozzle of a commercially available inkjet print-head to obtain the time-dependent velocity waveform. A printhead with a single transparent nozzle 80 μm in orifice diameter was used to eject single droplets at a speed of 5 m/s. An optical microscope was used with an ultra-high-speed camera to capture the motion of particles suspended in a transparent liquid at the center of the nozzle and above the fluid meniscus at a rate of half a million frames per second. Time-resolved velocity fields were obtained from a fluid layer approximately 200 μm thick within the nozzle for a complete jetting cycle. A Lagrangian finite-element numerical model with experimental measurements as inputs was used to predict the meniscus movement. The model predictions showed good agreement with the experimental results. This work provides the first experimental verification of physical models and numerical simulations of flows within a drop-on-demand nozzle. © 2012 Society for Imaging Science and Technology.

Numerical simulation of a bubble rising in shear-thinning fluids

The motion of a single bubble rising freely in quiescent non-Newtonian viscous fluids was investigated experimentally and computationally. The non-Newtonian effects in the flow of viscous inelastic fluids are modeled by the Carreau theological model. An improved level set approach for computing the incompressible two-phase flow with deformable free interface is used. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The simulation results demonstrate that the algorithm is robust for shear-thinning liquids with large density (rho(1)/rho(g) up to 10(3)) and high viscosity (eta(1)/eta(g) up to 10(4)). The comparison of the experimental measurements of terminal bubble shape and velocity with the computational results is satisfactory. It is shown that the local change in viscosity around a bubble greatly depends on the bubble shape and the zero-shear viscosity of non-Newtonian shear-thinning liquids. The shear-rate distribution and velocity fields are used to elucidate the formation of a region of large viscosity at the rear of a bubble as a result of the rather stagnant flow behind the bubble. The numerical results provide the basis for further investigations, such as the numerical simulation of viscoelastic fluids. (C) 2010 Elsevier B.V. All rights reserved.

Analysis of current induced by long internal solitary waves in stratified ocean

An approximate theoretical expression for the current induced by long internal solitary waves is presented when the ocean is continuously or two-layer stratified. Particular attention is paid to characterizing velocity fields in terms of magnitude, flow components, and their temporal evolution/spatial distribution. For the two-layer case, the effects of the upper/lower layer depths and the relative layer density difference upon the induced current are further studied. The results show that the horizontal components are basically uniform in each layer with a shear at the interface. In contrast, the vertical counterparts vary monotonically in the direction of the water depth in each layer while they change sign across the interface or when the wave peak passes through. In addition, though the vertical components are generally one order of magnitude smaller than the horizontal ones, they can never be neglected in predicting the heave response of floating platforms in gravitationally neutral balance. Comparisons are made between the partial theoretical results and the observational field data. Future research directions regarding the internal wave induced flow field are also indicated.

边坡工程的极限分析上限法

Theory of limit analysis include upper bound theorem and lower bound theorem. To deal with slope stability analysis by limit analysis is to approximate the real solution from upper limit and lower limit. The most used method of limit analysis is upper bound theorem, therefore it is often applied to slope engineering in many cases. Although upper bound approach of limit analysis can keep away from vague constitutive relation and complex stress analyses, it also can obtain rigorous result. Assuming the critical surface is circular slip surface, two kinematically admissible velocity fields for perpendicular slice method and radial slice method can be established according to the limit analysis of upper bound theorem. By means of virtual work rate equation and strength reduction method, the upper-bound solution of limit analysis for homogeneous soil slope can be obtained. A log-spiral rotational failure mechanism for homogeneous slope is discussed from two different conditions which represent the position of shear crack passing the toe and below the toe. In the dissertition, the author also establishes a rotational failure mechanics with combination of different logarithmic spiral arcs. Furthermore, the calculation formula of upper bound solution for inhomogeneous soil slope stability problem can be deduced based on the upper bound approach of rigid elements. Through calculating the external work rate caused by soil nail, anti-slide pile, geotechnological grid and retaining wall, the upper bound solution of safety factor of soil nail structure slope, slip resistance of anti-slide pile, critical height of reinforced soil slope and active earth pressure of retaining wall can be obtained by upper bound limit analysis method. Taking accumulated body slope as subject investigated, with study on the limit analysis method to calculate slope safety factor, the kinematically admissible velocity fields of perpendicular slice method for slope with broken slip surface is proposed. Through calculating not only the energy dissipation rate produced in the broken slip surfaces and the vertical velocity discontinuity, but also the work rate produced by self-weight and external load, the upper bound solution of slope with broken slip surface is deduced. As a case study, the slope stability of the Sanmashan landslide in the area of the Three Gorges reservoir is analyzed. Based on the theory of limit analysis, the upper bound solution for rock slope with planar failure surface is obtained. By means of virtual work-rate equation, energy dissipation caused by dislocation of thin-layer and terrane can be calculated; furthermore, the formulas of safety factor for upper bound approach of limit analysis can be deduced. In the end, a new computational model of stability analysis for anchored rock slope is presented after taking into consideration the supporting effect of rock-bolts, the action of seismic force and fissure water pressure. By using the model, not only the external woke-rate done by self-weight, seismic force, fissure water pressure and anchorage force but also the internal energy dissipation produced in the slip surface and structural planes can be totally calculated. According to the condition of virtual work rate equation in limit state, the formula of safety factor for upper bound limit analysis can be deduced.