Thermocapillary Instabilities of Low Prandtl Number Fluid in a Laterally Heated Vertical Cylinder

2006 ◽  
Vol 128 (6) ◽  
pp. 1228-1235 ◽  
Author(s):  
B. Xu ◽  
X. Ai ◽  
B. Q. Li
Keyword(s):  
Surface Tension ◽  
Free Surface ◽  
Three Dimensional ◽  
Vertical Cylinder ◽  
Reynolds Numbers ◽  
Base Flow ◽  
Matrix Structure ◽  
Various Boundary Conditions ◽  
Banded Matrix ◽  
High Order Finite Difference

Stabilities of surface-tension-driven convection in an open cylinder are investigated numerically. The cylinder is heated laterally through its sidewall and is cooled at free surface by radiation. A seeding crystal at constant temperature is in contact with the free surface. Axisymmetric base flow is solved using the high-order finite difference method. Three-dimensional perturbation is applied to the obtained base flow to determine the critical Marangoni numbers at which the axisymmetry is broken. The eigenvalue matrix equation is solved using linear fractional transformation with banded matrix structure taken into account. Critical Marangoni-Reynolds numbers are obtained at various boundary conditions.

Be´nard-Marangoni Instability in an Open Vertical Cylinder With Lateral Heating

2005 ◽  
Author(s):  
B. Xu ◽  
X. Ai ◽  
B. Q. Li
Keyword(s):  
Convective Flow ◽  
Marangoni Number ◽  
Vertical Cylinder ◽  
Base Flow ◽  
Matrix Structure ◽  
Governing Equations ◽  
Thermocapillary Effect ◽  
Lateral Heating ◽  
Banded Matrix ◽  
Critical Marangoni Number

A linear stability analysis of Rayleigh-Be´nard-Marangoni flow of low Prandtl number fluid contained in an open vertical cylinder is presented. The cylinder is heated laterally and is cooled at top surface by radiation. Governing equations of the flow are solved for axisymmetric base flow using higher order finite difference scheme. Small perturbation was applied to the obtained base flow to determine the critical Marangoni number and Grashof number at which the axisymmetry is broken. The eigenvalue matrix equation is solved using linear fractional transformation with banded matrix structure taken into account. It is found that the thermocapillary effect stabilizes the convective flow driven by buoyancy.

Entrainment and growth of vortical disturbances in the channel-entrance region

2021 ◽  
Vol 927 ◽  
Author(s):  
Pierre Ricco ◽  
Claudia Alvarenga
Keyword(s):  
Initial Conditions ◽  
Experimental Studies ◽  
Three Dimensional ◽  
Analytical Description ◽  
Reynolds Numbers ◽  
Boundary Region ◽  
Shear Layers ◽  
Base Flow ◽  
Open Boundary ◽  
Classical Stability

The entrainment of free-stream unsteady three-dimensional vortical disturbances in the entry region of a channel is studied via matched asymptotic expansions and by numerical means. The interest is in flows at Reynolds numbers where experimental studies have documented the occurrence of intense transient growth, despite the flow being stable according to classical stability analysis. The analytical description of the vortical perturbations at the channel mouth reveals how the oncoming disturbances penetrate into the wall-attached shear layers and amplify downstream. The effects of the channel confinement, the streamwise pressure gradient and the viscous/inviscid interplay between the oncoming disturbances and the boundary-layer perturbations are discussed. The composite perturbation velocity profiles are employed as initial conditions for the unsteady boundary-region perturbation equations. At a short distance from the channel mouth, the disturbance flow is mostly confined within the shear layers and assumes the form of streamwise-elongated streaks, while farther downstream the viscous disturbances permeate the whole channel although the base flow is still mostly inviscid in the core. Symmetrical disturbances exhibit a more significant growth than anti-symmetrical disturbances, the latter maintaining a nearly constant amplitude for several channel heights downstream before growing transiently, a unique feature not reported in open boundary layers. The disturbances are more intense as the frequency decreases or the bulk Reynolds number increases. We compute the spanwise wavelengths that cause the most intense downstream growth and the threshold wall-normal wavelengths below which the perturbations are damped through viscous dissipation.

Ice ripple formation at large Reynolds numbers

2012 ◽  
Vol 694 ◽  
pp. 225-251 ◽  
Author(s):  
Carlo Camporeale ◽  
Luca Ridolfi
Keyword(s):  
Free Surface ◽  
Complete Solution ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Thermal Boundary Conditions ◽  
Ripple Formation ◽  
High Reynolds ◽  
Lower Temperature ◽  
The Stability ◽  
Stefan Condition

AbstractA free-surface-induced morphological instability is studied in the laminar regime at large Reynolds numbers ($\mathit{Re}= 1\text{{\ndash}} 1{0}^{3} $) and on sub-horizontal walls ($\vartheta \lt 3{0}^{\ensuremath{\circ} } $). We analytically and numerically develop the stability analysis of an inclined melting–freezing interface bounding a free-surface laminar flow. The complete solution of both the linearized flow field and the heat conservation equations allows the exact derivation of the upper and lower temperature gradients at the interface, as required by the Stefan condition, from which the dispersion relationship is obtained. The eigenstructure is obtained and discussed. Free-surface dynamics appears to be crucial for the triggering of upstream propagating ice ripples, which grow at the liquid–solid interface. The kinematic and the dynamic conditions play a key role in controlling the formation of the free-surface fluctuations; these latter induce a streamline distortion with an increment of the wall-normal velocities and a destabilizing phase shift in the net heat transfer to the interface. Three-dimensional effects appear to be crucial at high Reynolds numbers. The role of inertia forces, vorticity, and thermal boundary conditions are also discussed.

On cusped interfaces

1998 ◽  
Vol 359 ◽  
pp. 313-328 ◽  
Author(s):  
YULII D. SHIKHMURZAEV
Keyword(s):  
Surface Tension ◽  
Free Surface ◽  
Reynolds Numbers ◽  
Present Theory ◽  
Similar Process ◽  
Surface Tension Gradient ◽  
Low Reynolds Numbers ◽  
Moving Contact Line ◽  
Moving Contact ◽  
Rate Of Dissipation

An asymptotic analysis of two-dimensional free-surface cusps associated with flows at low Reynolds numbers is presented on the basis of a model which, in agreement with direct experimental observations, considers this phenomenon as a particular case of an interface formation–disappearance process. The model was derived from first principles and earlier applied to another similar process: the moving contact-line problem. As is shown, the capillary force acting on a cusp from the free surface, which in the classical approach can be balanced by viscous stresses only if the associated rate of dissipation of energy is infinite, in the present theory is always balanced by the force from the surface-tension-relaxation ‘tail’, which stretches from the cusp towards the interior of the fluid. The flow field near the cusp is shown to be regular, and the surface-tension gradient in the vicinity of the cusp, caused and maintained by the external flow, induces and is balanced by the shear stress. Existing approaches to the free-surface cusp description and some relevant experimental aspects of the problem are discussed.

Interactions Between High-Viscosity Droplets Deposited on a Surface: Experiments and Simulations

2012 ◽  
Author(s):  
J. Esmaeelpanah ◽  
A. Dalili ◽  
S. Chandra ◽  
J. Mostaghimi ◽  
H. C. Fan ◽  
...  
Keyword(s):  
Surface Tension ◽  
Free Surface ◽  
Piecewise Linear ◽  
Three Dimensional ◽  
Surface Force ◽  
Line Length ◽  
High Viscosity ◽  
Droplet Shape ◽  
Viscous Liquids ◽  
Straight Lines

A combined numerical and experimental investigation of coalescence of droplets of highly viscous liquids dropped on a surface has been carried out. Droplets of 87 wt% glycerin-in-water solutions with viscosity 110 centistokes were deposited sequentially in straight lines onto a flat, solid steel plate and droplet impact photographed. Impacting droplets spread on the surface until liquid surface tension and viscosity overcame inertial forces and the droplets recoiled, eventually reaching equilibrium. Droplet center-to-center distance was varied and droplet line length was measured from photographs. As droplet spacing was increased there was less interaction between the droplets. A three dimensional parallel code has been developed to simulate fluid flow and free surface interaction by solving the continuity, momentum and volume-of-fluid (VOF) equations. The two-step projection method was employed to solve the governing equations for the whole domain including both liquid and air phases. The continuum-surface-force (CSF) scheme was applied to model surface tension and the piecewise-linear-interface-construction (PLIC) technique used to reconstruct the free surface. Computer generated images of impacting droplets modeled droplet shape evolution correctly and compared well with photographs taken during experiments. Accurate predictions were obtained for droplet line length during spreading and at equilibrium.

Numerical simulations of nonlinear axisymmetric flows with a free surface

1987 ◽  
Vol 178 ◽  
pp. 195-219 ◽  
Author(s):  
Douglas G. Dommermuth ◽  
Dick K. P. Yue
Keyword(s):  
Free Surface ◽  
Numerical Simulations ◽  
Three Dimensional ◽  
Vertical Cylinder ◽  
Free Surface Flows ◽  
Boundary Integral ◽  
The Body ◽  
Special Treatment ◽  
Computational Domain ◽  
Vortex Sheets

A numerical method is developed for nonlinear three-dimensional but axisymmetric free-surface problems using a mixed Eulerian-Lagrangian scheme under the assumption of potential flow. Taking advantage of axisymmetry, Rankine ring sources are used in a Green's theorem boundary-integral formulation to solve the field equation; and the free surface is then updated in time following Lagrangian points. A special treatment of the free surface and body intersection points is generalized to this case which avoids the difficulties associated with the singularity there. To allow for long-time simulations, the nonlinear computational domain is matched to a transient linear wavefield outside. When the matching boundary is placed at a suitable distance (depending on wave amplitude), numerical simulations can, in principle, be continued indefinitely in time. Based on a simple stability argument, a regriding algorithm similar to that of Fink & Soh (1974) for vortex sheets is generalized to free-surface flows, which removes the instabilities experienced by earlier investigators and eliminates the need for artificial smoothing. The resulting scheme is very robust and stable.For illustration, three computational examples are presented: (i) the growth and collapse of a vapour cavity near the free surface; (ii) the heaving of a floating vertical cylinder starting from rest; and (iii) the heaving of an inverted vertical cone. For the cavity problem, there is excellent agreement with available experiments. For the wave-body interaction calculations, we are able to obtain and analyse steady-state (limit-cycle) results for the force and flow field in the vicinity of the body.

scholarly journals Experiments and Simulations of Free-Surface Flow behind a Finite Height Rigid Vertical Cylinder

Fluids ◽  
2021 ◽  
Vol 6 (10) ◽  
pp. 367
Author(s):  
Valentin Ageorges ◽  
Jorge Peixinho ◽  
Gaële Perret ◽  
Ghislain Lartigue ◽  
Vincent Moureau
Keyword(s):  
Free Surface ◽  
Numerical Study ◽  
Three Dimensional ◽  
Vertical Cylinder ◽  
Free Surface Flow ◽  
Air Entrainment ◽  
Surface Flow ◽  
Finite Height ◽  
Wave Patterns ◽  
Different Diameters

We present the results of a combined experimental and numerical study of the free-surface flow behind a finite height rigid vertical cylinder. The experiments measure the drag and the wake angle on cylinders of different diameters for a range of velocities corresponding to 30,000 <Re< 200,000 and 0.2<Fr<2 where the Reynolds and Froude numbers are based on the diameter. The three-dimensional large eddy simulations use a conservative level-set method for the air-water interface, thus predicting the pressure, the vorticity, the free-surface elevation and the onset of air entrainment. The deep flow looks like single phase turbulent flow past a cylinder, but close to the free-surface, the interaction between the wall, the free-surface and the flow is taking place, leading to a reduced cylinder drag and the appearance of V-shaped surface wave patterns. For large velocities, vortex shedding is suppressed in a layer region behind the cylinder below the free surface. The wave patterns mostly follow the capillary-gravity theory, which predicts the crest lines cusps. Interestingly, it also indicates the regions of strong elevation fluctuations and the location of air entrainment observed in the experiments. Overall, these new simulation results, drag, wake angle and onset of air entrainment, compare quantitatively with experiments.

Instabilities of dynamic thermocapillary liquid layers with magnetic fields

1999 ◽  
Vol 382 ◽  
pp. 87-108 ◽  
Author(s):  
T. E. MORTHLAND ◽  
J. S. WALKER
Keyword(s):  
Magnetic Field ◽  
Free Surface ◽  
Metal Layer ◽  
Three Dimensional ◽  
Previous Treatment ◽  
Magnetic Reynolds Number ◽  
Base Flow ◽  
Linear Temperature ◽  
Linear Temperature Gradient ◽  
Critical Marangoni Number

This paper presents a linear-stability analysis for the transition from a steady, two-dimensional thermocapillary convection in a liquid-metal layer to a periodic, three-dimensional flow involving hydrothermal waves which propagate in the directionzw normal to the plane of the base flow. There is a uniform magnetic field applied parallel to the free surface in the plane of the base flow, and there is a linear temperature gradient along the free surface in the base flow. The ratio of the layer's length to its depth, 2L, is large. The magnetic Reynolds number is small.A key parameter is λ, the ratio of the large Hartmann number based on depth to L. The value of λ increases as either the magnetic field strength is increased or L is decreased. The results for very small values of λ agree with the results of a previous treatment of this instability without a magnetic field. As λ is increased, the critical Marangoni number and the wavenumber for the hydrothermal rolls both increase. For large values of λ, the base flow and the hydrothermal waves are confined to a free-surface layer with O(λ−2) dimensionless thickness.

scholarly journals Transient growth in the flow past a three-dimensional smooth roughness element

2013 ◽  
Vol 724 ◽  
pp. 642-670 ◽  
Author(s):  
S. Cherubini ◽  
M. D. De Tullio ◽  
P. De Palma ◽  
G. Pascazio
Keyword(s):  
Boundary Layer ◽  
Boundary Layer Flow ◽  
Roughness Element ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Base Flow ◽  
Energy Gain ◽  
Roughness Elements ◽  
Layer Flow ◽  
Optimal Disturbances

AbstractThis work provides a global optimization analysis, looking for perturbations inducing the largest energy growth at a finite time in a boundary-layer flow in the presence of smooth three-dimensional roughness elements. Amplification mechanisms are described which can bypass the asymptotical growth of Tollmien–Schlichting waves. Smooth axisymmetric roughness elements of different height have been studied, at different Reynolds numbers. The results show that even very small roughness elements, inducing only a weak deformation of the base flow, can localize the optimal disturbance characterizing the Blasius boundary-layer flow. Moreover, for large enough bump heights and Reynolds numbers, a strong amplification mechanism has been recovered, inducing an increase of several orders of magnitude of the energy gain with respect to the Blasius case. In particular, the highest value of the energy gain is obtained for an initial varicose perturbation, differently to what found for a streaky parallel flow. Optimal varicose perturbations grow very rapidly by transporting the strong wall-normal shear of the base flow, which is localized in the wake of the bump. Such optimal disturbances are found to lead to transition for initial energies and amplitudes considerably smaller than sinuous optimal ones, inducing hairpin vortices downstream of the roughness element.

Sign in / Sign up

Export Citation Format

Share Document