Lecture 17: Natural (Free) convection

1. Lecture 17: Natural (Free) convection http://www.youtube.com/watch?v=qIlcXp-cKHg 1. Assume that gravitational field affects the flow

2. 2. Statics: 3. We separate the static parts and the non-static contributions: 4. Assume that 5. Linearization of the pressure term 6. Substituting this into Navier-Stokes equation or

3. 7. In incompressible flow, density variations are due to temperature non-isothermalities Finally, the governing equations for thermal convection in Boussinesq approximation Let us non-dimensionalise these equations.

4. Dimensions of the viscosity and thermal conductivity coefficients Dimension of the heat expansion coefficient Gravity acceleration -- length scale (typical size) -- time-scale (also called convective time scale) -- scale of velocity -- scale of pressure -- scale of temperature (Θ is the typical temperature difference) Non-dimensionalisation

5. Non-dimensionalization of the Navier-Stokes equation: or Heat transfer equation, or Non-dimensionalization of the continuity equation is trivial:

6. Non-dimensionalised equations of thermal convection Non-dimensional parameters: -- Grashof number (defines the strength of the forcing (buoyancy) term; characterises intensity of the convective flow) -- Prandtl number (characterises the fluid properties; defines the relative importance of thermal conduction and convection as two mechanisms of heat transfer, Pr>>1: convection dominates, Pr<<1: convection can be disregarded)

7. Czochralski process The Czochralski process is a one of the methods of crystal growth from melt, used to obtain single crystals of semiconductors (e.g. silicon, germanium and gallium arsenide), metals (e.g. palladium, platinum, silver, gold), salts, and synthetic gemstones. When the silicon is fully melted (~15000C), a small seed crystal mounted on the end of a rotating shaft is slowly lowered until it just dips below the surface of the molten silicon. The seed crystal's rod is slowly pulled upwards. Silicon ingot Main difficulty: strong convective flows in the melt. From one hand, convection makes the melt more uniform, which is positive. From another hand, strong flows near the crystallization plane introduce crystal defects. Convection is controlled by rotation/vibration of ingot and crucible, by magnetic field, …

8. Solutal convection Consider incompressible isothermal fluid flow with an admixture. Equation of mass conservation for an impurity: (this equation is valid for not very large C) Diffusion coefficient Variations in density can be caused by admixture: Equations of solutal convection: P.S. If flow is non-isothermal, the thermosolutal (or double-diffusive) convective flow is induced.

9. Lecture 1Flows driven by surface tension gradients Spatial variations in surface tension result in added tangential stresses at an interface, giving rise to fluid motions in the underlying bulk liquids. The motion induced by tangential gradients of surface tension is termed the Marangoni effect. Examples: the camphor ball ‘dances’ on a water surface (http://www.youtube.com/watch?v=Pe88T45VdR8), the calming effect of ‘oil troubled waters’ (http://www.youtube.com/watch?v=00PPPt7EJqo).

10. This motion is induced at an interface, where, besides the normal force, there is another force tangential to the surface, . Adding this force, we obtain the following boundary condition is the unit normal vector directed into medium 1. The surface tension coefficient is function of temperature; the thermocapillary flows are induced. function of concentration; the difusocapillary flows. function of electric charge; the electrocapillary flows.

11. Thermocapillary motion in a thin liquid layer The thermocapillary motion generated in an open rectangular shallow pan with a very thin layer at the bottom. The difference in side wall temperatures results in a temperature gradient along the surface. For liquids, z T2 (cold) T1 (hot) h2 h(x) h1 x The flow in nearly lateral. Any flow nonuniformities at the side walls are small and the flow is essentially 2D. l α α2 α1 x

12. Except at the side walls, which are far removed from the bulk flow, the vertical velocity component is very much smaller than the horizontal component and any effects of the free surface curvature may be neglected. We assume that the liquid layer is thin enough to neglect inertial effects, implying Re<<1 (‘shallow water’ approximation). Governing equations: x-projection of the momentum equation, (1) z-projection, (2) The continuity equation (in integral form) for the fully developed flow is (3)

13. (4) 1. Bottom (z=0): 2. free surface (z=h): Boundary conditions: Or, for our flow, (5) (6) Integration of (1) with the use of (4) and (5) gives (7) Integration of (2) with the use of (6) gives (8)

14. (9) (8) specifies the relation between the pressure gradient and variation of free surface height in the x-direction as By using (7) and (9), condition (3) results in Here, we assumed that α=α1 and h=h1 at x=0. Or, The velocity profile: umax

15. or The requirement Re<<1 may be expressed as The variation of αwith T for liquids is close to linear. For water, Thus with This would give a value Hence, and We did not consider the gravity-driven convection. The gravity between the gravity and capillary forces is defined by the Bond number: Here, -- capillary length

16. For water/air interface, . And, for our configuration, This means that, for the configuration examined, the flows driven by the surface tension effect (Marangoni convection) dominate (in comparison with the gravity-driven convective flows). The surface tension gradients can be important in such very thin layers of mm size or less, or in a reduced gravity environment. For example, a crystal grown from its melt under reduced gravity is governed by convection driven by thermally induced surface tension gradients rather than buoyancy forces.

17. Ludwig Prandtl (4 February 1875, Freising, Upper Bavaria – 15 August 1953, Gottingen) was a German scientist. He was a pioneer in the development of rigorous systematic mathematical analyses which he used to underlay the science of aerodynamics, which have come to form the basis of the applied science of aeronautical engineering. In the 1920s he developed the mathematical basis for the fundamental principles of subsonic aerodynamics in particular. 17

18. Carlo Giuseppe Matteo Marangoni (29 April 1840, Pavia Italy – 14 April 1925, Florence, Italy) was an Italian physicist. He held a position of High School Physics Teacher at the Liceo Dante (Florence) for 45 years, until retirement in 1916. He mainly studied surface phenomena in liquids and also contributed to meteorology. Franz Grashof (July 11, 1826 in Dusseldorf - October 26, 1893 in Karsruhe) was a German engineer. He was a professor of Applied Mechanics at the Technische Hochschule Karlsruhe. He developed some early steam-flow formulas but made no significant contribution to free convection. Joseph Valentin Boussinesq (13 March 1842, Saint-Andre-de-Sangonis, France –Died 19 February 1929, Paris, France) was a French mathematician and physicist who made significant contributions to the theory of hydrodynamics, vibration, light, and heat. 18

19. Jan Czochralski (23 October 1885 – 22 April 1953) was a Polish chemist who invented the Czochralski process, which is used for growing single crystals and in the production of semiconductor wafers. He was educated at Charlottenburg Polytechnic in Berlin, where he specialized in metal chemistry. Czochralski discovered the Czochralski method in 1916, when he accidentally dipped his pen into a crucible of molten tin. He immediately pulled his pen out to discover that a thin thread of solidified metal was hanging from the nib. The nib was replaced by a capillary, and Czochralski verified that the crystallized metal was a single crystal. Czochralski's experiments produced single crystals a millimeter in diameter and up to 150 centimeters long. In 1950, the method was used to grow single germanium crystals, leading to its use in semiconductor production.