Two Bifurcation Transition Processes in Floating Half Zone Convection of Larger Prandtl Number Fluid

Processes of the onset oscillation in the thermocapillaxy convection under the Earth's gravity are investigated by the numerical simulation and experiments in a floating half zone of large Prandtl number with different volume ratio. Both computational and experimental results show that the steady and axisymmetric convection turns to the oscillatory convection of m=1 for the slender liquid bridge, and to the oscillatory convection before a steady and 3D asymmetric state for the case of a fat liquid bridge. It implies that, there are two critical Marangoni numbers related, respectively, to these two bifurcation transitions for the fat liquid bridge. The computational results agree with the results of ground-based experiments.

Two-Dimensional Kinematic Wave Model of Overland-Flow

A two-dimensional kinematic wave model was developed for simulating runoff generation and flow concentration on an experimental infiltrating hillslope receiving artificial rainfall. Experimental observations on runoff generation and flow concentration on irregular hillslopes showed that the topography of the slope surface controlled the direction and flow lines of overland flow. The model-simulated results satisfactorily compared with experimental observations. The erosive ability of the concentrated flow was found to mainly depend on the ratio of the width and depth of confluent grooves.

High Aspect Ratio (L/D) Riser Viv Prediction Using Wake Oscillator Model

A two-dimensional (2-D) vortex-induced vibration (VIV) prediction model for high aspect ratio (LID) riser subjected to uniform and sheared flow is studied in this paper. The nonlinear structure equations are considered. The near wake dynamics describing the fluctuating nature of vortex shedding is modeled using classical van der Pol equation. A new approach was applied to calibrate the empirical parameters in the wake oscillator model. Compared the predicted results with the experimental data and computational fluid dynamic (CFD) results. Good agreements are observed. It can be concluded that the present model can be used as simple computational tool in predicting some aspects of VIV of long flexible structures. (C) 2008 Elsevier Ltd. All rights reserved.

Experimental Investigation On The Holdup Distribution Of Oil-Water Two-Phase Flow In Horizontal Parallel Tubes

An experimental investigation was conducted to study the holdup distribution of oil and water two-phase flow in two parallel tubes with unequal tube diameter. Tests were performed using white oil (of viscosity 52 mPa s and density 860 kg/m(3)) and tap water as liquid phases at room temperature and atmospheric outlet pressure. Measurements were taken of water flow rates from 0.5 to 12.5 m(3)/h and input oil volume fractions from 3 to 94 %. Results showed that there were different flow pattern maps between the run and bypass tubes when oil-water two-phase flow is found in the parallel tubes. At low input fluid flow rates, a large deviation could be found on the average oil holdup between the bypass and the run tubes. However, with increased input oil fraction at constant water flow rate, the holdup at the bypass tube became close to that at the run tube. Furthermore, experimental data showed that there was no significant variation in flow pattern and holdup between the run and main tubes. In order to calculate the holdup in the form of segregated flow, the drift flux model has been used here.

Large Scale Numerical Simulations For Multi-Phase Fluid Dynamics With Moving Interfaces

This short communication presents our recent studies to implement numerical simulations for multi-phase flows on top-ranked supercomputer systems with distributed memory architecture. The numerical model is designed so as to make full use of the capacity of the hardware. Satisfactory scalability in terms of both the parallel speed-up rate and the size of the problem has been obtained on two high rank systems with massively parallel processors, the Earth Simulator (Earth simulator research center, Yokohama Kanagawa, Japan) and the TSUBAME (Tokyo Institute of Technology, Tokyo, Japan) supercomputers.

Thermovibrational Instability Of Rayleigh-Marangoni-Benard Convection In Two-Layer Fluid Systems

The thermovibrational instability of Rayleigh-Marangoni-Benard convection in a two-layer system under the high-frequency vibration has been investigated by linear instability analysis in the present paper. General equations for the description of the convective flow and within this framework, the generalized Boussinesq approximation are formulated. These equations are dealt with using the averaging method. The theoretical analysis results show that the high-frequency thermovibrations can change the Marangoni-Benard convection instabilities as well as the oscillatory gaps of the Rayleigh-Marangoni-Benard convection in two-layer liquid systems. It is found that vertical high-frequency vibrations can delay convective instability of this system, and damp the convective flow down. (C) 2007 COSPAR. Published by Elsevier Ltd. All rights reserved.

Derivation Of A Nonlinear Reynolds Stress Model Using Renormalization Group Analysis And Two-Scale Expansion Technique

Adopting Yoshizawa's two-scale expansion technique, the fluctuating field is expanded around the isotropic field. The renormalization group method is applied for calculating the covariance of the fluctuating field at the lower order expansion. A nonlinear Reynolds stress model is derived and the turbulent constants inside are evaluated analytically. Compared with the two-scale direct interaction approximation analysis for turbulent shear flows proposed by Yoshizawa, the calculation is much more simple. The analytical model presented here is close to the Speziale model, which is widely applied in the numerical simulations for the complex turbulent flows.

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.

Validity ranges of interfacial wave theories in a two-layer fluid system

In the present paper, we endeavor to accomplish a diagram, which demarcates the validity ranges for interfacial wave theories in a two-layer system, to meet the needs of design in ocean engineering. On the basis of the available solutions of periodic and solitary waves, we propose a guideline as principle to identify the validity regions of the interfacial wave theories in terms of wave period T, wave height H, upper layer thickness d(1), and lower layer thickness d(2), instead of only one parameter-water depth d as in the water surface wave circumstance. The diagram proposed here happens to be Le Mehautes plot for free surface waves if water depth ratio r = d(1)/d(2) approaches to infinity and the upper layer water density rho(1) to zero. On the contrary, the diagram for water surface waves can be used for two-layer interfacial waves if gravity acceleration g in it is replaced by the reduced gravity defined in this study under the condition of sigma = (rho(2) - rho(1))/rho(2) -> 1.0 and r > 1.0. In the end, several figures of the validity ranges for various interfacial wave theories in the two-layer fluid are given and compared with the results for surface waves.

A dynamic coupling model for hybrid atomistic-continuum computations

A dynamic coupling model is developed for a hybrid atomistic-continuum computation in micro- and nano-fluidics. In the hybrid atomistic-continuum computation, a molecular dynamics (MD) simulation is utilized in one region where the continuum assumption breaks down and the Navier-Stokes (NS) equations are used in another region where the continuum assumption holds. In the overlapping part of these two regions, a constrained particle dynamics is needed to couple the MD simulation and the NS equations. The currently existing coupling models for the constrained particle dynamics have a coupling parameter, which has to be empirically determined. In the present work, a novel dynamic coupling model is introduced where the coupling parameter can be calculated as the computation progresses rather than inputing a priori. The dynamic coupling model is based on the momentum constraint and exhibits a correct relaxation rate. The results from the hybrid simulation on the Couette flow and the Stokes flow are in good agreement with the data from the full MD simulation and the solutions of the NS equations, respectively. (c) 2007 Elsevier Ltd. All rights reserved.

Attractor crisis and bursting in a fluid flow with two no-slip directions

Characteristic burtsing behavior is observed in a driven, two-dimensional viscous flow, confined to a square domain and subject to no-slip boundaries. Passing a critical parameter value, an existing chaotic attractor undergoes a crisis, after which the flow initially enters a transient bursting regime. Bursting is caused by ejections from and return to a limited subdomain of the phase space, whereas the precrisis chaotic set forms the asymptotic attractor of the flow. For increasing values of the control parameter the length of the bursting regime increases progressively. Passing another critical parameter value, a second crisis leads to the appearance of a secondary type of bursting, of very large dynamical range. Within the bursting regime the flow then switches in irregular intervals from the primary to the secondary type of bursting. Peak enstrophy levels for both types of bursting are associated to the collapse of a primary vortex into a quadrupolar state.

A coupled model for multiphase fluid flow and sedimentation deformation in oil reservoir and its numerical simulation

A mathematical model for coupled multiphase fluid flow and sedimentation deformation is developed based on fluid-solid interaction mechanism. A finite difference-finite element numerical approach is presented. The results of an example show that the fluid-solid coupled effect has great influence on multiphase fluid flow and reservoir recovery performances, and the coupled model has practical significance for oilfield development.

Potential energy model for turbulent entrainment across density interface in a linearly stratified fluid

A potential energy model is developed for turbulent entrainment in the absence of mean shear in a linearly stratified fluid. The relation between the entrainment distance D and the time t and the relation between dimensionless entrainment rate E and the local Richardson number are obtained. An experiment is made for examination. The experimental results are in good agreement with the model, in which the dimensionless entrainment distance D is given by DBAR = A(i)(SBAR)-1/4(fBAR)1/2(tBAR)1/8, where A(i) is the proportional coefficient, S is the dimensionless stroke, fBAR is the dimensionless frequency of the grid oscillation, tBAR the dimensionless time.

Growth of zero-mean-shear turbulent-mixed layer with suspended particles in a 2-layer fluid system

The growth behaviour of zero-mean-shear turbulent-mixed layer containing suspended solid particles has been studied experimentally and analysed theoretically in a two-layer fluid system. The potential model for estimating the turbulent entrainment rate of the mixed layer has also been suggested, including the results of the turbulent entrainment for pure two-layer fluid. The experimental results show that the entrainment behaviour of a mixed layer with the suspended particles is well described by the model. The relationship between the entrainment distance and the time, and the variation of the dimensionless entrainment rate E with the local Richardson number Ri1 for the suspended particles differ from that for the pure two-layer fluid by the factors-eta-1/5 and eta-1, respectively, where eta = 1 + sigma-0-DELTA-rho/DELTA-rho-0.

The wave model of mesothermal plasma near wakes and Korteweg-de Vries equation

The stationary two-dimensional (x, z) near wakes behind a flat-based projectile which moves at a constant mesothermal speed (V∞) along a z-axis in a rarefied, fully ionized, plasma is studied using the wave model previously proposed by one of the authors (VCL). One-fluid theory is used to depict the free expansion of ambient plasma into the vacuum produced behind a fast-moving projectile. This nonstationary, one-dimensional (x, t) flow which is approximated by the K-dV equation can be transformed, through substitution, t=z/V∞, into a stationary two-dimensional (x, z) near wake flow seen by an observer moving with the body velocity (V∞). The initial value problem of the K-dV equation in (x, t) variables is solved by a specially devised numerical method. Comparisons of the present numerical solution for the asymptotically small and large times with available analytical solutions are made and found in satisfactory agreements.