Compressibility and collisional effects on thermal instability of a partially ionized medium

Similar documents
Pressure And Entropy Variations Across The Weak Shock Wave Due To Viscosity Effects

Compressibility Effects

Kinetics of Particles. Chapter 3

Chapter 3 Kinematics in Two Dimensions; Vectors

Phys101 Final Code: 1 Term: 132 Wednesday, May 21, 2014 Page: 1

Module 4: General Formulation of Electric Circuit Theory

Chapter 4. Unsteady State Conduction

The Destabilization of Rossby Normal Modes by Meridional Baroclinic Shear

Surface and Contact Stress

The influence of a semi-infinite atmosphere on solar oscillations

February 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA

Suggested reading: Lackmann (2011), Sections

Yeu-Sheng Paul Shiue, Ph.D 薛宇盛 Professor and Chair Mechanical Engineering Department Christian Brothers University 650 East Parkway South Memphis, TN

Flipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System

Short notes for Heat transfer

Flipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System

ChE 471: LECTURE 4 Fall 2003

Electric Current and Resistance

CHAPTER 8b Static Equilibrium Units

Problem 1 Known: Dimensions and materials of the composition wall, 10 studs each with 2.5m high

Free Vibrations of Catenary Risers with Internal Fluid

Dispersion Ref Feynman Vol-I, Ch-31

OF SIMPLY SUPPORTED PLYWOOD PLATES UNDER COMBINED EDGEWISE BENDING AND COMPRESSION

Numerical Simulation of the Thermal Resposne Test Within the Comsol Multiphysics Environment

CLASS XI SET A PHYSICS

Kinematic transformation of mechanical behavior Neville Hogan

Heat and Mass Transfer Analysis of an Unsteady MHD Flow Past an Impulsively Started Vertical Plate in Presence of Thermal Radiation

Advanced Heat and Mass Transfer by Amir Faghri, Yuwen Zhang, and John R. Howell

ON THE EFFECTIVENESS OF POROSITY ON UNSTEADY COUETTE FLOW AND HEAT TRANSFER BETWEEN PARALLEL POROUS PLATES WITH EXPONENTIAL DECAYING PRESSURE GRADIENT

7.0 Heat Transfer in an External Laminar Boundary Layer

A PLETHORA OF MULTI-PULSED SOLUTIONS FOR A BOUSSINESQ SYSTEM. Department of Mathematics, Penn State University University Park, PA16802, USA.

FIELD QUALITY IN ACCELERATOR MAGNETS

Lecture 17: Free Energy of Multi-phase Solutions at Equilibrium

1/2 and e0 e s ' 1+ imm w 4 M s 3 πρ0 r 3 m. n 0 ktr. .Also,since n 0 ktr 1,wehave. 4 3 M sπρ 0 r 3. ktr. 3 M sπρ 0

(1.1) V which contains charges. If a charge density ρ, is defined as the limit of the ratio of the charge contained. 0, and if a force density f

Examiner: Dr. Mohamed Elsharnoby Time: 180 min. Attempt all the following questions Solve the following five questions, and assume any missing data

PHYS 314 HOMEWORK #3

Introduction to Smith Charts

DYNAMIC MODELLING OF N-CARDAN TRANSMISSIONS WITH SHAFTS IN SPATIAL CONFIGURATION. Part II. THE ALGORITHM OF DYNAMIC MODELLING

Lecture 5: Equilibrium and Oscillations

ENGINEERING COUNCIL CERTIFICATE LEVEL THERMODYNAMIC, FLUID AND PROCESS ENGINEERING C106 TUTORIAL 5 THE VISCOUS NATURE OF FLUIDS

CANKAYA UNIVERSITY FACULTY OF ENGINEERING MECHANICAL ENGINEERING DEPARTMENT ME 313 HEAT TRANSFER

Aircraft Performance - Drag

Assume that the water in the nozzle is accelerated at a rate such that the frictional effect can be neglected.

THE TOPOLOGY OF SURFACE SKIN FRICTION AND VORTICITY FIELDS IN WALL-BOUNDED FLOWS

Computational modeling techniques

Lab 11 LRC Circuits, Damped Forced Harmonic Motion

37 Maxwell s Equations

SGP - TR - 30 PROCEEDINGS FOURTH WORKSHOP GEOTHERMAL RESERVOIR ENGINEERING. Editors. December13-15, , 1978 SGP - TR - 30 CONF

CHEM-443, Fall 2013, Section 010 Midterm 2 November 4, 2013

Curvature Effects on Thermal Buckling Load of DWCNT Under Axial Compression Force

Analysis of Hydrodynamics and Heat Transfer in a Thin Liquid Film Flowing Over a Rotating Disk by the Integral Method

Unit code: H/ QCF level: 5 Credit value: 15 OUTCOME 3 - STATIC AND DYNAMIC FLUID SYSTEMS TUTORIAL 3 - VISCOSITY

Entropy. Chapter The Clausius Inequality and Entropy

Revision: August 19, E Main Suite D Pullman, WA (509) Voice and Fax

A Simple Set of Test Matrices for Eigenvalue Programs*

On Boussinesq's problem

Study Guide Physics Pre-Comp 2013

Part a: Writing the nodal equations and solving for v o gives the magnitude and phase response: tan ( 0.25 )

Materials Engineering 272-C Fall 2001, Lecture 7 & 8 Fundamentals of Diffusion

Chapter 2 GAUSS LAW Recommended Problems:

Chapter 30. Inductance

General Chemistry II, Unit I: Study Guide (part I)

39th International Physics Olympiad - Hanoi - Vietnam Theoretical Problem No. 1 /Solution. Solution

CHAPTER 3 INEQUALITIES. Copyright -The Institute of Chartered Accountants of India

A Few Basic Facts About Isothermal Mass Transfer in a Binary Mixture

Quantum Harmonic Oscillator, a computational approach

ZVS Boost Converter. (a) (b) Fig 6.29 (a) Quasi-resonant boost converter with M-type switch. (b) Equivalent circuit.

PHYS College Physics II Final Examination Review

Effects of piezo-viscous dependency on squeeze film between circular plates: Couple Stress fluid model

ENGI 4430 Parametric Vector Functions Page 2-01

14. Which shows the direction of the centripetal force acting on a mass spun in a vertical circle?

Review Problems 3. Four FIR Filter Types

Spring Pendulum with Dry and Viscous Damping

More Tutorial at

14. Which shows the direction of the centripetal force acting on a mass spun in a vertical circle?

SAFE HANDS & IIT-ian's PACE EDT-04 (JEE) Solutions

Thermodynamics Partial Outline of Topics

GAUSS' LAW E. A. surface

1 Course Notes in Introductory Physics Jeffrey Seguritan

SPH3U1 Lesson 06 Kinematics

ELECTRON CYCLOTRON HEATING OF AN ANISOTROPIC PLASMA. December 4, PLP No. 322

22.54 Neutron Interactions and Applications (Spring 2004) Chapter 11 (3/11/04) Neutron Diffusion

Chapter 6. Dielectrics and Capacitance

3. Mass Transfer with Chemical Reaction

CHAPTER 4 Dynamics: Newton s Laws of Motion /newtlaws/newtltoc.html

Dead-beat controller design

1 The limitations of Hartree Fock approximation

Calculating the optimum pressure and temperature for vacancy minimization from theory; Niobium is an example. Jozsef Garai

12. Shallow water equations with rotation Poincaré waves. Considering now motions with L<<R, we can write the equations of motion in

Optimization of frequency quantization. VN Tibabishev. Keywords: optimization, sampling frequency, the substitution frequencies.

AP CHEMISTRY CHAPTER 6 NOTES THERMOCHEMISTRY

Engineering Approach to Modelling Metal THz Structures

Chapter 5: Force and Motion I-a

1. Transformer A transformer is used to obtain the approximate output voltage of the power supply. The output of the transformer is still AC.

Lecture 7: Damped and Driven Oscillations

Fundamental Concepts in Structural Plasticity

Phys102 Final-061 Zero Version Coordinator: Nasser Wednesday, January 24, 2007 Page: 1

making triangle (ie same reference angle) ). This is a standard form that will allow us all to have the X= y=

ECE 546 Lecture 02 Review of Electromagnetics

Transcription:

Pram~na, Vl. I0, N. 3, March 978, pp. 267-272, printed in India. Cmpressibility and cllisinal effects n thermal instability f a partially inized medium R C SHARMA and K C SHARMA Department f Mathematics, Himachal Pradesh University, Simla 7 005 MS received 23 June 977 Abstract. The thermal instability f a finitely cnducting hydrmagnetic cmpsite and cmpressible medium is studied t include the frictinal effects with neutrals. The effect f cmpressibility is fund t be stabilizing. In cntrast t the nnscillatry mdes fr (Cp/g)~> in the absence f a magnetic field; C~,/3 and g being specific heat at cnstant pressure, unifrm adverse temperature gradient and acceleratin due t gravity respectively, the presence f magnetic field intrduces scillatry mdes in the system. The verstable case is als discussed. The magnetic field is fund t have a stabilizing effect n the system fr (Cp/g)~>. Keywrds. Cmpressibility; cllisins; thermal instability.. Intrductin The prblem f thermal instability f fluids, bth in hydrdynamics and hydrmagnetics, has been treated in detail by Chandrasekhar (96). The Bussinesq apprximatin has been used thrughut. The equatins gverning the system becme quite cmplicated when the fluids are cmpressible. Spiegel and Vernis (960) made the fllwing assumptins: (i) the fluctuatins in density, pressure and temperature intrduced due t mtin, d nt exceed their ttal static variatins and (ii) the depth f fluid layer is much less than the scale height as defined by them. Using the abve assumptins, Spiegel and Vernis (960) have fund the flw equatins t be the same as fr incmpressible fluids except that the static temperature gradient is replaced by its excess ver the adiabatic. Fllwing Hans (968) the medium has been idealized as a cmpsite mixture f a hydrmagnetic (inized) cmpnent and a neutral cmpnent, the tw interacting thrugh mutual cuisinal (frictinal) effects. As it is quite frequent that the medium is nt fully inized and may be permeated with neutral atms, the abve idealizatin f the medium as a cmpsite mixture was called fr. Hans (968) als fund that these cllisins have stabilizing effect n the Rayleigh-Taylr instability. Sharma (976) studied the thermal hydrmagnetic instability f a partially inized medium. The bject f the present paper is t study the cmpressibility and cllisinal effects n thermal instability f a cmpsite medium. 267

268 R C Sharma and K C Sharma 2. Frmulatin f the prblem and dispersin relatin Cnsider an infinite hrizntal cmpressible and cmpsite layer f thickness d cnsisting f a finitely cnducting hydrmagnetic fluid f density p, permeated with neutrals f density Pa, acted n by a unifrm vertical magnetic field H(0, 0, H) and gravity frce g(0, 0, --g). This layer is heated frm belw such that a steady adverse temperature gradient/3(-----[dt/dz[) is maintained. Bth the inized fluid and neutral atms are assumed t behave like cntinuum fluids. Fllwing Spiegel and Vernis (960), Sharma (976), the basic linearized equatins gverning the mtin f cmpressible and cmpsite medium are: _ pcgq=_ V 8P + g Sp + pv V~q q-~rr (ST h) H -t- pav (qa--q), () 8t ~gq d = _ vc (qa -- q), (2) cgt ~gh --V (q H) q- ~ ~7~h, (3) ~gt V.q:O, V.h=O, (4) ~O_[[36._~_]wt ~ +K A~O, (5) /gt ~ C t, I where q(u, v, w), h(h~, h, hz), 0, 8p and 8p dente respectively the perturbatins in velcity, magnetic field H, temperature T, density p and pressure p; g[ C, t~, v (-=--/~/pm), K',, (:, '/pmcp), vc, q~ and ~ stand fr the adiabatic gradient, the viscsity, the kinematic viscsity, the thermal cnductivity, the thermal diffusivity, the cllisinal frequency between the tw cmpnents f the cmpsite medium, the velcity f the neutral cmpnent and the resistivity respectively, am(= a, say) is the cefficient f thermal expansin and pm (= p, say) is the density f the inized medium. In writing eq. (2), it is assumed that the effects n the neutral cmpnent resulting frm the fields f gravity and pressure are neglected and that the neutral particles are nncnducting. Analyzing in terms f nrmal mdes, we seek slutins f the abve equatins, whse dependence n space crdinates x, y, z and time t is f the frm [w, 0, hz] = [W(z), O(z), K(z)] exp (ik~xf ikyyj-nt), (6) where k,, k, are hrizntal wave numbers f the harmnic disturbance, k~=kx2+k, ~ and n is the frequency. Eliminating qd between eqs (I) and (2) and using expressin (6), eqs (I)-(5) give,c dg/, i W --a)dk- (DS--a 2) ( D ~ -- a ~ -- + a"~v~-~vl W + 4rrpv gadz a 20 = 0, (7) V

Thermal instability f cmpressible and cmpsite medium 269 (D2--a~--V~a)K=--(H-~-) DW, (8) (D~--a~--pl~) @ =- -- fl-- g IV, (9) K where a==kd, ~,~-nd~/v, px=v/k, pz~v]~l, %=pa]p, D=d/dz and x, y, z stand fr the crdinates in the new unit f length d. Use has als been made f the Bussinesq equatin f state 8p= -- apo. Eliminating O and K between eqs (7)-(9), we get ( D~--a~)( D2--a~--pz~) [(~-a'-,, ~,)(D'-a'-,, l + +vd,l~! a vcda/v ~ -- QD'] W = -- R(~.~)(D'--a'--p,g)a'W, (0) where R--~gafld4/vtc is the Rayleigh number, Q:Had~/4,rpv77 is the Chandrasekhar number and G:Cdt]g. Cnsider the case in which bth the bundaries are free and the medium adjining the fluid is nncnducting. The bundary cnditins apprpriate fr the prblem are (Chandrasekhar 96): W=D~W=O, 0=0 X=0 and h is cntinuus at z---o and. () In the absence f any surface current, the tangential cmpnents f magnetic field are cntinuus. Hence the bundary cnditins in additin t () are DK=O, n the bundaries. Using the abve bundary cnditins, it can be shwn that all the even rder derivatives f W must vanish fr z=o and and hence the prper slutin f W characterizing the lwest mde is (2) W= Wshl*r z, (3) where W 0 is a cnstant. Substituting (3) in eq. (0), we btain the dispersin relatin Rx =(~)I( + b)( + b-l-p~-~) b( +b+p.~)

270 R C Sharma and K C Sharma where a 2 _ Q = R and b =--. 3. The statinary cnvectin When instability sets in as rdinary cnvectin, the marginal state will be characterized by =0 and eq. (4) reduces t (5) Fr the fixed value f Q, let the nndimensinal number G accunting fr the cmpressibility effects be als kept as fixed, then we find that where Rc and Rc dente respectively the critical Rayleigh number in the absence and presence f cmpressibility. The effect f cmpressibility is, thus, t pstpne the nset f thermal instability. Hence we btain a stabilizing effect f cmpressibility. The cases G< and G= crrespnd t negative and infinite values f critical Rayleigh numbers in the presence f cmpressibility which are nt relevant fr the present prblem. 4. Stability f the system and nn-scillatry mdes Multiplying eq. (7) by W*, the cmplex cnjugate f W, and using eqs (8) and (9), we get s ~ a~/~ ~ '7 (8 + v, ~*,) i~ + + g-~..~, ~/ i~ + ff-~p~ + Cp a Ka z (i 6 + Px * 6) = 0, (7),,( -- G) where I~= f (ID~WI~+2a~IDWI~+a4WI~)dz, ~=f (I D~KI ~ + 2a~l DKI ~ + a~ IKI ~) dz,

Thermal instabibty f cmpressible and cmpsite medium 27 ~=f (I D Kl~ a'l Kl')dz, z~= f (ID~l' a'l~l')dz, ~= f (I[~)az, (is) which are all psitive definite and ~* is the cmplex cnjugate f ~. Putting ~ -----,, + i - t and then equating the real and imaginary parts f eq. (7), we btain 6 + (~'+" d~/~) {~"+~"+~" (" d~/") (I +%)} Z~+ ~ (I~+p~, Z~) (~, + v~d2/v) 2+ et ~ 4~pv Cp ~,ca ~ t r + p~ c,,/6) = 0, +,,(I_G) vs (9) and [~ ' + a, (red'/v) (I + %) + (v d'/v)" (I + %) + ~,' z~ (~, + v~ d~/v) ~ + ~fl P2/4 q- Cp a Ka 9 ] 4~pv 7(6:~ p~ z, :. (20) It fllws frm eq. (20) that if G>I and if the magnetic field is absent, cry=0, which means that the scillatry mdes are nt allwed fr G>I and in the absence f magnetic field. The presence f magnetic field, in cntrast t nnscillatry mdes fr G > in the absence f magnetic field, intrduces scillatry mdes in the system. 5. The verstable ease In this sectin we cnsider the pssibility f whether instability may ccur as an verstability. Put G,/~-=igl, it being remembered that, may be cmplex. Since fr verstability, we wish t determine the critical Rayleigh number fr the nset f instability via a state f pure scillatins, it suffices t find cnditins fr which (4) will admit f slutin with, real. Equatin (4) becmes I G \ [- I ( II b( + b + ipagl) (2)

272 R C Sharma and K C Sharma Equating real and imaginary parts f eq. (2) and eliminating R~ between them, we btain (t,~2,:~/a ~) [p:0 ~, + 0 +b)( +pb (,:~/d'~)] ~: + (l+b) [p~ ( +p+%) +Pl (l+b)(:rid ~) %vc + ( +b)" (l+pb (,::/a') + Q~ (p~-p,) (,::/a,)]~' + (l+b) [( +b)uvc 9 (I+p+%) + Q: (p:-pz)vc ~] = 0. (22) Equatin (22) is quadratic in ~,2 and des nt allw any f its rts t be such that Re(c,a) is psitive s that vl is imaginary when Pl >~P2. Hence Pl >~P2 i.e. ~ ~<~, (23) is a sufficient cnditin fr the nnexistence f verstability. The cnditin (23) is the same as in the absence f cmpressibility as well as frictinal (cllisinal) effects with neutrals n thermal instability (Chandrasekhar 96). Thus K ~<~7 is a sufficient cnditin fr the nnexistence f verstability which hlds bth in the presence r absence (Chandrasekhar 96) f cmpressibility and cllisinal effects with neutrals n thermal instability f a cmpsite and cmpressible medium. T study the effect f magnetic field n thermal instability f a cmpressible and cmpsite medium, we examine the nature f dr/dq. It fllws frm eq. (4) that (24) The imaginary part f eq. (24) equated t zer gives pl =,02. (25) Equating real parts f eq. (24) and substituting (25) in it, we btain dr _ (26) It fllws frm eq. (26) that dr:/dq is psitive if G>I. Hence fr G>I, the Rayleigh number increases as the magnetic field increases, shwing the stabilizing effect f magnetic field. The cases G< and G:-I are nt allwed here as thse wuld mean negative and infinite Rayleigh numbers respectively which are nt relevant in the present discussin. References Chandrasekhar S 96 Hydrdynamic and ttydrmagnetic Stability (Oxfrd: Clarendn Press) Chap. 4 Hans H K 968 Nucl. Fusin 8 89 Sharma R C 976 Physica C8 99 Spiegel E A and Vernis G 960 Astrphys. J. 3 442