NUMERICAL INVESTIGATION OF SURFACE DISTURBANCE ON LIQUID METAL FALLING FILM FLOWS IN A MAGNETIC FIELD GRADIENT D. Gao, N.B. Morley
|
|
- Camilla Garrett
- 5 years ago
- Views:
Transcription
1 MAGNETOHYDRODYNAMICS Vol. 4 (4), No., pp NUMERICAL INVESTIGATION OF SURFACE DISTURBANCE ON LIQUID METAL FALLING FILM FLOWS IN A MAGNETIC FIELD GRADIENT D. Gao, N.B. Morley Department of Mechanical and Aerospace Engineering University of California, Los Angeles, CA 995 The surface stability characters of falling film flows on an inclined plane in a spanwise magnetic field with constant streamwise gradient are investigated by the direct numerical simulation of the Navier Stokes equations and magnetic induction equation based on the Volume of Fluid (VOF) method for free surface tracking. The spatial disturbance, a small rectangular rise or a sinusoidal bump, is imposed on the otherwise fully developed flat MHD film. The free surface behaviors under the disturbances are calculated based on a D MHD model for the film flow down an inclined plane. It is found that the surface disturbance is rapidly enhanced in streamwise etent and in amplitude, with the traveling speed of disturbance largely reduced. These MHD surface stability features are discussed and eplained in terms of the induced electrical currents and thereby the MHD drag due to the magnetic field gradient. Introduction. Interest in liquid metal film flow arises from the possibility of utilizing this type of flow in future magnetic confined fusion reactors for protection of solid structures from the thermonuclear plasma []. This application places increasing demand on understanding the stability of the magnetohydrodynamic (MHD) falling thin film down an inclined plane. For the surface stability of ordinary films without magnetic field, there are etensive theoretical studies from linear theory [], weak non-linear [3], to the advanced highly non-linear theory [4]. However, with the presence of a magnetic field the comple interaction between liquid metal motion and magnetic field through the Lorentz force J B significantly increases the difficulty for theoretical handling, thus most analyses are limited to linear theory scopes. For eample, Hsieh [5] analyzed the stability of conducting films down an inclined plane in a uniform transverse magnetic field (perpendicular to the plane) and came to the conclusion that the vertical field acts to stabilize film flows. Some works [6, 7] considered thin film flows in a uniform spanwise magnetic field and studied the Hartmann layer effects on the velocity profile. To the authors knowledge, the literature in nonlinear analysis of MHD film flow instability is scarce. The purpose of this study is to eplore the effects of a field gradient on the surface stability behaviors under some finite-sized initial disturbances. There are past works that describe the gradient effects for closed-channel flows [8, 9]. But, as we know, the research for field gradient effects on the surface stability of film flow is still very limited. MHD free surface flows are more difficult to analyze than the closed-channel flows, because the induced currents depend on the positions of unrestrained surfaces. In this work, we consider a spanwise magnetic field oriented in the z-direction, and assume the field has a constant field gradient in streamwise through the complete running length. The thin film flows are studied in the twodimensional geometry (, y) and with the velocity (u, v). The configuration of the two-dimensional film flow is shown in Fig.. The wall is electrically nonconducting. 99
2 h u g B z (Spanwise) η y (Transverse) (Streamwise) Fig.. Geometry of the D MHD thin film in a spanwise magnetic field (here η is the surface variation from the fully developed flat film height h ). β Quite different from the Hartmann type problem [8, ], in which the induced currents are located in the cross sections perpendicular to the streamwise, the currents induced by the field gradient encircle in the planes parallel to the streamwise. The Hartmann layers [8, ] here have been removed to infinity, and are neglected in this study, so that the D film flows under a constant streamwise field gradient can be studied by the D-MHD model presented in our previous work []. In that work, the analytic solution to the equilibrium state is obtained and a preliminary linear stability analysis is performed with the approimation of undisturbed induced currents. The linear analysis, though, is not applicable to finite amplitude disturbances, and the MHD stability is beyond the capability of present nonlinear theory. Hence, a direct numerical simulation of two-dimensional MHD film flows is desired. Since the deformation of the film surface under a disturbance may be dramatic, a simulation method with no topological constraint is preferable. We choose to simulate the time-dependent governing equations and fully non-linear free surface boundary conditions by a finite volume scheme based on the Volume of Fluid (VOF) method to track free surfaces. VOF methods are widely used in the free surface simulations [, 3]. The governing equations and numerical methods are described, but we will put emphasis on the results and discussion for understanding field gradient induction effects on the film instability under the imposed disturbances. The spatial disturbances, a rectangular rise or a sinusoidal bump, are imposed on the otherwise fully developed flat films, and the evolutions of the disturbances are calculated.. Governing equations and numerical model. We consider a thin liquid metal layer of constant density ρ, dynamic viscosity µ, electrical conductivity σ e, magnetic permeability µ m, draining down a plate inclined at the angle β to the horizontal in a magnetic environment, as shown in Fig.. The velocity field V = ( u(, y, t), v(, y, t) ) and the del operator = ( /, / y, ) are considered two-dimensional in the -y plane. The applied magnetic field B = (,,B(, y, t) ) is aligned in the spanwise z-direction and can in general vary both temporally and spatially in the -y plane. The induced currents in turn build their own magnetic field b, which is also aligned spanwise. The field gradient interacts with liquid metal movement, producing currents [8, 9], and thereby the Lorentz force in the -y plane. Thus the field gradient effect on flow can be studied by the D model [] (other variants of this infinite and/or aisymmetric D model can be seen in [4]). We actually use scalar variables B and b to represent the applied and induced spanwise magnetic fields, because only one direction of field is involved
3 here. The incompressible liquid metal flow is described by the time-dependent electromagnetic induction equation, momentum equations and mass conservation [] as: b +(V )b = t Rm b (V )B B t () V t +(V )V = p + Re V + Fr ĝ + m4 ( bẑ) (B + b)ẑ ReRm () V =. (3) Besides the standard hydrodynamic parameters: the Reynolds number Re and the Froude number Fr, we also have the MHD parameters: the magnetic Reynolds Rm and m. The dimensionless group m is actually equal to the square root of the Hartmann number based on a different choice of B. They are defined respectively by: Re = u h ν, Fr = u (gh ), ( B σe h ρν Rm = u h σ e µ m, m = ) /, where ĝ, p are the unit gravity acceleration and pressure respectively; u, h,and B are the characteristic values, respectively, for velocity, length and magnetic field. Time is normalized by h /u, and pressure by ρu. The interest here in the study of stability prompts the selection of h, u as the film height and average velocity of a fully developed film. Induction equation () reveals that the source term of the induced magnetic field, or in other words, the electrical currents, comes from the spatial or temporal variation of the applied field. It does not depend at all on the absolute value of the applied field. Unlike the Aitov analysis [6], which studied the Hartmann layer effect on the mean flow velocity profile, this model is proposed to study the field gradient effect in the absense of the Hartmann layers due to the infinite width of film flow. In the above governing equations, we use the induced field b to represent the induced magnetic field, where, actually, the contours of b work as the streamlines of induced currents, since we have the following relation for the current density: J = bẑ (4) where J is normalized by B /(h µ m ). The induced electrical currents are confined in the liquid metal, and the gas side is considered non-conducting, thus the integral current density with respect to y over the film height is zero. Hence, the MHD boundary conditions at the bottom solid wall and free surface can be written as: b(y =)= (5) b(y = h) b(y =)= (6) When the liquid is draining downstream, the liquid-gas interface h[, t] changes with the streamwise location and time, which is defined by the kinematic equation v(y = h) = h + u(y = h) h t, (7) This equation associates the surface movement to the velocity components at the interface. At a free surface the normal stress is balanced by the capillary force,
4 and the shear stress by the gradient of surface tension along the surface. For the constant surface tension σ, the dynamic conditions at the interface can be written in inde form as T ij ˆn j ˆn i = κwe, T ijˆt j ˆn i = (8) T ij = pδ ij + Re ( ju i + i u j ), We = σ ρu h, (9) where δ ij is the delta function, κ is the curvature of the interface, and We is the Weber number. The unit outward normal vector, tangential vector, and curvature of the interface are given by ˆn = ( h/, )[ +( h/ ) ] /, () ˆt = (, h/ )[ +( h/ ) ] /, () κ = ( h/ )[ +( h/ ) ] 3/. () In the computer code written to model this general D MHD film flow problem, the free surface location is tracked by the VOF method []. The dynamic conditions at free surface are implemented via the continuum surface force (CSF) model [5] combined with the VOF method. In the VOF method, a scalar variable f[, y, t] is defined as the volume fraction of liquid in a control volume, which is advected in a flow field. The VOF method consists of two parts: an interface reconstruction algorithm to determine an interface from the given volume fraction data, and a VOF advection algorithm to determine the volume fraction data after a time step. The first part is required by the second part for maintaining mass conservation. The split operator method from Puckett et al. [3] is used for VOF advection, and the free surface of each cell is constructed as a line-segment. The numerical scheme is designed for the two-dimensional MHD flows with free surfaces, and is built on an Eulerian staggered grid. The velocity field is obtained by solving the Navier Stokes equations and the magnetic induction equation. The Navier Stokes equations are solved by a projection method. The projection method is divided into two steps. First, the intermediate velocity field is determined from the diffusion-convection equations without regarding the solenoidal nature of the velocity; second, the pressure Poisson equation is solved implicitly, and then the pressure field is used to correct the intermediate velocity field to recover the divergence-free nature. The magnetic induction equation is solved implicitly in a fashion similar to the pressure Poisson equation, not causing additional CFL constraint. Complete descriptions of the VOF method can be found in the references [, 6]. The numerical scheme has been validated in regard to the Navier Stokes solver by simulating thin water film flows [6, 7], and magnetic mechanism implementation by computing MHD films and comparing with analytical solution and other numerical results [6].. Solution to a fully developed MHD film with a constant streamwise field gradient. This paper specifically studies the special case of film flow stability in a spanwise magnetic field with a constant streamwise gradient in field strength. Such a field can be epressed dimensionally as B z = C +B,where B z / = C. To highlight the field gradient effect, B is chosen to be equal to the product of field gradient multiplying characteristic length, B = Ch/ (the factor of is used to be consistent with [], where it is adopted to simplify the form of the analytic solutions). With this definition of B and value of B then the current
5 source term in Eq. reduces to: (V )B B =u, (3) t and the MHD force term in Eq. reduces to: m 4 ( bẑ) (B + b)ẑ = m4 ReRm ReRm (b + + B ) b, (4) where B, and for the rest of the work, is the normalized value of the y-intercept of the linearly varying B field at =. The dimensionless applied field is actually B z = + B. This specific special case is chosen because an analytic solution for the equilibrium flat film can be obtained [] and is written as: Re sin β u = Fr m ( cosmsinh my cos(m my)+coshmsin my cosh(m my) sinh m +sinm sinmy cosh my(sin m sinh m cos m cosh m) + sinh m +sinm ) cosmy sinh my(sin m sinh m +cosmcosh m) + sinh m +sinm ReRm sin β b = Fr ( m 4 cosmcosh my sin(m my)+coshmcos my sinh(m my) + sinh m +sinm cosmy sinh my(cos m cosh m sin m sinh m) + sinh m +sinm ) sinmy cosh my(cos m cosh m +sinmsinh m) + sinh m +sinm (5) (6) p = m4 ReRm b( + B )+ cos β ( y) (7) Fr The relationship between Fr and Re is Re sin β 4m 3 (sin m +sinhm) Fr = (8) +3cosm +3coshm 8cosmcosh m The dimensional average velocity is given by ( ) gh u av = sin β νm 3 ( ) (9) cos m +coshm cosmcosh m +(cosm cosh m) sin m +sinhm The solution reduces to ordinary film solution U = 3 (y y )andu av = gh sin β/3ν if there is no field gradient m. The solution shows very definitely the origin of M-shaped velocity profile due to the field gradient induction effect working on the flow via pressurization. The solution also shows a dependence of side layer thickness on Ha / when the Ha definition is based on field gradient not on field magnitude. The flowrate, decreasing with the increase of Ha, is evidence of the MHD drag. 3
6 t = 3 4 t = t = 3 4 t =6 3 4 Fig.. The development of a small rectangular rise on the MHD film at m =. 3. Numerical results and discussion. We intend to roughly eamine the responses of MHD films to various disturbance conditions and magnetic conditions. Here, the parameters for the vertical (β =9 ) MHD film are chosen as: Re =, We =., Rm =., and different m values based on the strength of the field gradient. The field is assumed to be increasing in strength from the inlet at =, and various initial values of B are considered. The small value of We insteadofwe=ischoseninordertohighlight the magnetic effect and also to suppress noise caused probably by some unphysical reasons. The film is initiated by its fully developed velocity profile given by equation (5), the inlet is fied at the same velocity profile. The no-slip boundary condition at the solid-liquid interface and continuous outflow are used for the film flow. The computational area is 4 initial film depths in the -direction and in the y-direction. The mesh size is (, y) =(.,.). 4
7 t = Fig. 3. Development of the bump on the ordinary film without MHD effect. 3.. Small rise on the flat film surface. In this first test, a small square wave (broadband in frequency content) is imposed on the liquid surface at point =4,andB = is assumed. The small rise on the surface is created by two cells.5 y higher than the fully developed flat MHD film. The velocities for the two cells are equal to the surface velocity of the flat film calculated from Eq.(5). The free surface configurations after a few time units are shown in Fig.. Without surface tension, the disturbance is quickly stimulated and amplified around the initial disturbance to emerge as sinusoidal bumps. The film selects a short wave to develop around the small rectangular disturbance, and this short wave grows in amplitude as propagating downstream. In contrast, on the counterpart ordinary film with Re =, We =.andwithout magnetic effect, we found that the small disturbance is swept downstream without perceivable growth within the computation domain. 3.. Sinusoidal bump at m =,, 3. We have seen after the broad-banded small disturbance develops to finite amplitude, it quickly ehibits a finite sinusoidal bump shape. So we will directly impose an upper-half of single period spatial disturbance of the form. sin[π/] on the free surface of the otherwise fully developed vertical film flows, and then observe the development of the disturbance at various magnetic gradient conditions. The velocity for the finite amplitude bump is also calculated from the respective velocity profile of Eq.(5) with the distance of respective cell center to the bottom wall substituted in the relation. Figs. 3 5 show the development of initial disturbance on fully developed vertical films at m =,, 3withB =. Without magnetic effect, Fig. 3 shows the wave propagates downstream at phase speed about.7 without apparent growth in amplitude during the simulation period. The bump is shifted downstream by the dominant inertia at this level of Re =. No other waves are observed ahead or behind the waves, most energy is confined to the fundamental one. In contrast, the disturbance development on the MHD film flow at m =,as shown in Fig. 4, and at m = 3, as shown in Fig. 5, ehibits strong nonlinear dynamics. The original disturbance first spreads over a larger region, then is rapidly amplified. The initial sinusoidal shape of the bump quickly degrades to irregular disturbance. Accompanying the spreading process, the rear part of disturbance is increased in amplitude more rapidly than the front part. The growth rate of amplitude increases appreciably with m. Although the preliminary linear stability analysis given in [] points out that the instability is epected to increase with m for long waves, there must be another factor directly causing the growth and spreading of the disturbance since the velocities of base film for m =andm =, as shown in Fig. 6, are very similar. The small difference of the velocity profile, which was the only factor included in the linear stability analysis of [], should not lead to such a change in stability behaviors. 5
8 t = 3 4 t = t = 3 4 t =6 3 4 Fig. 4. The development of the bump on the MHD film at m =. The traveling velocity of surface disturbance decreases apparently with the increase of m. The disturbance front travels faster than the tail causing the elongated disturbance region. From the simulation result, the tail traveling velocities are estimated.7,.7, respectively for m =,, 3. Referring to Fig. 6 for the solutions to the equilibrium films, the surface liquid velocity increases with m, approaching at m = 3, while being.5 for an ordinary film. So we can see for MHD flows that the surface disturbance travels much slower than even the surface velocity of the base film flow. This seems contradictory to the results for ordinary surface waves [8, 9] without magnetic effect, where the surface wave travels at least at the surface fluid velocity. There must be a mechanism distinctive from the gravity, viscosity for the MHD flow, which opposes the disturbance propagation. It is the MHD drag produced by induction effect between the magnetic field gradient and the fluid motion, which is often the dominant mechanism for MHD flows. 6
9 t = 3 4 t = t = 3 4 t =6 3 4 Fig. 5. The development of the bump on the MHD film at m =3. The MHD effect must be understood from the induced current inside a flow. For eample, induced b contours for m =caseatt = 8 are shown in Fig. 7. Incidently, we can see here that the induced magnetic field is substantially smaller than the applied field B, thus its contribution to the total magnetic field is negligible. But more to the point here, though, the b contours represent the streamlines of electric currents. The directions of currents are downstream near the surface and upstream near the bottom wall. The induced currents follow the curvature of free surface. For the equilibrium film the currents are straight horizontal lines, thus the Lorentz force is perpendicular to the surface, and it is actually balanced by pressure resulting in a distorted M-shape velocity profile. The MHD drag effect for the equilibrium film is epressed by deceasing flowrate with the increase of m values. However, when a disturbance stands out of flat surface, the Lorentz force points to the upstream-downward direction on the front edge, to downstream-downward direction on the tailing edge. The backward force is larger than the forward force 7
10 ..5 U..5 m =3 m = m = Fig. 6. Velocity profiles of the fully developed flat flow y due to the increase applied field, resulting in a net backward resistant force on the outstanding disturbance. The larger the field gradient, the larger the MHD resistance. As a result, the disturbance traveling is slowed down by the drag. In addition, when the MHD forces on the front and tail slopes cannot be balanced by other mechanisms, they then act accumulating liquid there, making the growth in amplitude and change in shape. If the field gradient is large, the net drag force nearly does not allow the disturbance to spread, all the energy is mainly used to increase the amplitude. Once large disturbances arise from the surface in strongly MHD dominant conditions, such as m =, 3 considered here, the Lorentz force tends to destabilize the surface rapidly Effect of absolute values of B. The absolute value of B does not enter the solutions of velocity and b profiles for the fully developed flat film. In other words, whatever the initial value of B,thesamevelocityandb fields for flat film will result. The pressure solution Eq.(7) shows the electromagnetic force incurred from uniformly adding to or subtracting a value from the applied field can be easily balanced by pressure for equilibrium film since the current is strictly flowing in the streamwise direction. Additionally, in [] it is shown after utilizing the shallow water approimation to eliminate inertial and viscous terms in the y-equation of motion, that the velocity of the flow itself should be independent of the absolute value of B as well. But in the finite disturbance cases analyzed here there is an effect on the wave propagation and stability arising from the absolute values of B because the pressure has potential influences on the free surface deformation. 3 4 Fig. 7. b contours at form =film. 8
11 t = 3 4 t = t = 3 4 t =6 3 4 Fig. 8. Development of the bump on the MHD films for B increasing from. We calculated two cases: one with the start value of B =, and another with, for m = and other hydrodynamic parameters and initial disturbances equal to those in the preceding section. The evolution of the surface shape for these two new cases are shown in Figs. 8 and 9. Comparing Fig. 8 to the result in Fig. 4 from the previous section with B =, we do see some effect on the flow behavior. While the main body of the disturbance appears to be about the same length and traveling roughly at the same speed, we do see some small upstream and downstream wavelets that indicate the effect is spreading more rapidly for the B = case. But more dramatic is the increase in the wave amplitude with time (roughly twice higher at t = ) in the tail region of the disturbance with B =. The presence of a stronger field accelerates the disturbance growth rate in the later times, growing out of the confines of the computational space by t = 4. At earlier times the wave amplitudes of the B =andb =cases 9
12 t = 3 4 t = t = 3 4 t =6 3 4 Fig. 9. Development of the bump on the MHD films for B increasing from. are similar in magnitude, with the acceleration taking place as the disturbance becomes more non-linear. This fact leads support to the idea that the dependence on the absolute strength of B is a non-linear effect, neglected in the shallow water approimation used in []. The induced current streamlines for the B =case are shown in Fig. a and can be compared to Fig. 7 at B =. Contrary to the B =case,b = ehibits less amplitude growth of the initial surface disturbance and seems to hold its initial sinusoidal shape for a longer time. The leading edge of the disturbance, however, does seems to be spread faster downstream, elongating the disturbance especially after crossing the = point where the field actually passes through zero and changes the direction. The induced current streamlines for the B = case are shown in Fig. b and show the presence of a second current cell closing in the elongated disturbance.
13 (a) 3 4 (b) 3 4 Fig.. b contours at t = 8 for the case of increasing B from (a), and (b). 4. Conclusions. We have investigated the surface stability behaviors of vertical liquid metal film flows in a constant streamwise field gradient by numerical computation. The destabilizing role of a field gradient is clearly shown by the simulation results of spatial disturbance on the surface of a thin MHD film flow with very small surface tension. The small disturbance on a flat MHD film rapidly develops to a finite amplitude short wavelength disturbance. There is no tendency for the finite amplitude wave to stay in a sinusoidal shape. They tend to choose their own irregular shapes. The disturbance spreads to large region and grows in amplitude. The traveling speed is largely reduced by the MHD drag even to lower than the surface liquid velocity of a corresponding flat film. The effect of the induced currents, via the direct Lorentz force on the disturbed surface, but not the flow velocity profile play the key role for understanding the surface instability behaviors. The induced currents follow the free surface curvature. The preliminary linear stability analysis [] without consideration of disturbed current lines thus is no applicable to finite amplitude short waves although it may be applicable to infinitesimal long waves. The larger anti-streamwise component of the Lorentz force on the front slope of disturbance and the smaller forward force on the tailing slope due to the increase in field strength along the flow result in the net backward MHD drag force. The MHD drag slows down the traveling of surface disturbance and destabilizes the surface. Acknowledgements. The authors would like to gratefully acknowledge the support of U.S. Department of Energy through Grand No. DE-FG3-86ER53, and the support of Professor Mohamed Abdou at UCLA. The first author wishes to thank Professor Dhir at UCLA for his advice and instruction.
14 REFERENCES [] M.A. Abdou, A. Ying, N. Morley, et al. On the eploration of innovative concepts for fusion chamber technology APEX Interim Report Overview. Fusion Engineering and Design, vol. 54 (), pp [] C.-S. Yih. Stability of liquid flow down an inclined plane. Phys. Fluids, vol.6 (963), pp [3] B.J. Benney. Long waves in liquid films. J. Math. Phys., vol. 45 (966), pp [4] H.-C. Chang. Wave evolution on a falling film. Annual Review of Field Mechanics., vol. 6 (994), pp [5] D.Y. Hsieh. Stability of a conducting liquid flowing down an inclined plane in a magnetic field. Phys. Fluids, vol.8 (965), pp [6] T.N. Aitov, E.M. Kirillina, A.V. Tananaev. Stability of the flow of a thin layer of liquid metal in a coplanar magnetic field. Magnetohydrodynamics, vol. 4 (988), pp. 5. [7] N.B. Morley, M.S. Tillack. Eamination of stability calculations for liquid metal film flows in a coplanar magnetic field. Magnetohydrodynamics, vol. 9 (993), pp [8] H. Branover. Magnetohydrodynamic Flows in Ducts (Israel University Press, Jerusalem, Israel, 978). [9] O. Lielausis. Liquid-metal magnetohydrodynamics. In: Atomic Energy Review. (International Atomic Energy, Vienna, 975) vol. 3, pp [] J.C.R. Hunt, J.A. Shercliff. Magnetohydrodynamics at high Hartmann number. In: Annual Review of Fluid Mechanics. (Palo Alto, CA, USA: Annual Review Inc.) vol. 3 (97), pp [] D. Gao, N.B. Morley. Equilibrium and initial linear stability analysis of liquid metal falling film flows in a varying spanwise magnetic field. Magnetohydrodynamics, vol. 38 (), no. 4, pp [] D.B. Kothe, R.C. Mjolsness, M.D. Torrey. RIPPLE: A Computer Program for Incompressible Flows with Free Surfaces (LA-7-MS, Los Alamos National Laboratory, 99). [3] E.G. Puckett, A.S. Almgren, J.B. Bell, et al. A high-order projection method for tracking fluid interfaces in variable density incompressible flows. Journal of Computational Physics, vol. 3 (997), pp [4] N.B. Morley, S. Smolentsev, D. Gao. Modeling infinite/aisymmetric liquid metal magnetohydrodynamic free surface flows. Invited paper at the 6th International Symposium for Fusion Nuclear Technology (San Diego, USA Apr ). Fusion Engineering and Design, vol (), pp [5] J.U. Brackbill, D.B. Kothe, C. Zemach. A continnum method for modeling surface tension. Journal of Computational Physics, vol. (99), pp
15 [6] D. Gao. Numerical Simulation of Surface Wave Dynamics of Liquid Metal MHD Flow on an Inclined Plane in a Magnetic Field with Spatial Variation. Ph.D. Dissertation, Mechanical Engineering Department, University of California Los Angeles, 3. [7] D. Gao, N.B. Morley, V. Dhir. Numerical simulation of wavy falling film flows using the VOF method. Journal of Computational Physics, vol. 9 (3), pp [8] B.E. Anshus. On the asymptotic solution to the falling film stability problem. Ind. Eng. Chem. Fundam., vol. (97), pp [9] J.M. Floryan, S.H. Davis, R.E. Kelly. Instability of a liquid film flowing down a slightly inclined plane. Physics of Fluids, vol.3 (987), pp Received
Understanding Magnetic Field Gradient Effect From a Liquid Metal Droplet Movement
Donghong Gao* Neil B. Morley Vijay Dhir Department of Mechanical and Aerospace Engineering, University of California, Los Angeles, CA 90095 Understanding Magnetic Field Gradient Effect From a Liquid Metal
More informationNumerical study of liquid metal film flows in a varying spanwise. magnetic field
Numerical study of liquid metal film flows in a varying spanwise magnetic field D. Gao, N.B. Morley, V. Dhir Department of Mechanical and Aerospace Engineering, University of California, Los Angeles, CA
More information3D MHD Free Surface Fluid Flow Simulation Based on. Magnetic-Field Induction Equations
3D MHD Free Surface Fluid Flow Simulation Based on Magnetic-Field Induction Equations H.L. HUANG, A. YING, M. A. ABDOU Mechanical and Aerospace Engineering Department, UCLA, Los Angeles, CA, 90095 Huang@fusion.ucla.edu
More informationMagnetohydrodynamic Flow of a Liquid Metal in a Curved Circular Duct subject to the Effect of an External Magnetic Field
Paper 85 Civil-Comp Press, 2012 Proceedings of the Eighth International Conference on Engineering Computational Technology, B.H.V. Topping, (Editor), Civil-Comp Press, Stirlingshire, Scotland Magnetohydrodynamic
More informationSTABILITY ANALYSIS FOR BUOYANCY-OPPOSED FLOWS IN POLOIDAL DUCTS OF THE DCLL BLANKET. N. Vetcha, S. Smolentsev and M. Abdou
STABILITY ANALYSIS FOR BUOYANCY-OPPOSED FLOWS IN POLOIDAL DUCTS OF THE DCLL BLANKET N. Vetcha S. Smolentsev and M. Abdou Fusion Science and Technology Center at University of California Los Angeles CA
More informationVALIDATION CASE RESULTS FOR 2D AND 3D MHD SIMULATIONS
VALIDATION CASE RESULTS FOR 2D AND 3D MHD SIMULATIONS Ming-Jiu Ni 1, 3, Ramakanth Munipalli 2, Neil B Morley 1, Peter Huang 2, Mohamed A. Abdou 1 1 MAE Dept., University of California, Los Angeles, CA,
More informationLiquid Metal Flow Control Simulation at Liquid Metal Experiment
Modestov DOI:10.1088/1741-4326/aa8bf4 FIP/P4-38 Liquid Metal Flow Control Simulation at Liquid Metal Experiment M. Modestov 1,2, E. Kolemen 3, E. P. Gilson 3, J. A. Hinojosa 1, H. Ji 3, R. P. Majeski 3,
More informationNumerical Solution for Coupled MHD Flow Equations in a Square Duct in the Presence of Strong Inclined Magnetic Field
International Journal of Advanced Research in Physical Science (IJARPS) Volume 2, Issue 9, September 2015, PP 20-29 ISSN 2349-7874 (Print) & ISSN 2349-7882 (Online) www.arcjournals.org Numerical Solution
More informationLINEAR STABILITY ANALYSIS FOR THE HARTMANN FLOW WITH INTERFACIAL SLIP
MAGNETOHYDRODYNAMICS Vol. 48 (2012), No. 1, pp. 147 155 LINEAR STABILITY ANALYSIS FOR THE HARTMANN FLOW WITH INTERFACIAL SLIP N. Vetcha, S. Smolentsev, M. Abdou Fusion Science and Technology Center, UCLA,
More informationApplication of the immersed boundary method to simulate flows inside and outside the nozzles
Application of the immersed boundary method to simulate flows inside and outside the nozzles E. Noël, A. Berlemont, J. Cousin 1, T. Ménard UMR 6614 - CORIA, Université et INSA de Rouen, France emeline.noel@coria.fr,
More informationFLUID MECHANICS. Chapter 9 Flow over Immersed Bodies
FLUID MECHANICS Chapter 9 Flow over Immersed Bodies CHAP 9. FLOW OVER IMMERSED BODIES CONTENTS 9.1 General External Flow Characteristics 9.3 Drag 9.4 Lift 9.1 General External Flow Characteristics 9.1.1
More informationVORTICITY FIELD EVOLUTION IN A FORCED WAKE. Richard K. Cohn Air Force Research Laboratory Edwards Air Force Base, CA 92524
Proceedings of the st International Symposium on Turbulence and Shear Flow Phenomena, Santa Barbara, CA, Sep. 5, 999, Eds. Banerjee, S. and Eaton, J. K., pp. 9-96. VORTICITY FIELD EVOLUTION IN A FORCED
More informationarxiv: v4 [physics.comp-ph] 21 Jan 2019
A spectral/hp element MHD solver Alexander V. Proskurin,Anatoly M. Sagalakov 2 Altai State Technical University, 65638, Russian Federation, Barnaul, Lenin prospect,46, k2@list.ru 2 Altai State University,
More informationNonlinear shape evolution of immiscible two-phase interface
Nonlinear shape evolution of immiscible two-phase interface Francesco Capuano 1,2,*, Gennaro Coppola 1, Luigi de Luca 1 1 Dipartimento di Ingegneria Industriale (DII), Università di Napoli Federico II,
More informationVORTICITY GENERATION IN NON-UNIFORM MHD FLOWS
The 15 th Riga and 6 th PAMIR Conference on Fundamental and Applied MHD Invited Lectures VORTICITY GENERATION IN NON-UNIFORM MHD FLOWS S. Cuevas 1,S.Smolentsev 2, M. Abdou 2 1 Center for Energy Research,
More informationHEAT TRANSFER ENHANCEMENT USING VORTEX PROMOTERS IN MAGNETO-HYDRO-DYNAMIC FLOWS
Eleventh International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 7-9 December 2015 HEAT TRANSFER ENHANCEMENT USING VORTEX PROMOTERS IN MAGNETO-HYDRO-DYNAMIC FLOWS
More informationCHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW
CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW 4.1 Introduction Boundary layer concept (Prandtl 1904): Eliminate selected terms in the governing equations Two key questions (1) What are the
More informationMHD Analysis of Dual Coolant Pb-17Li Blanket for ARIES-CS
MHD Analysis of Dual Coolant Pb-17Li Blanket for ARIES-CS C. Mistrangelo 1, A. R. Raffray 2, and the ARIES Team 1 Forschungszentrum Karlsruhe, 7621 Karlsruhe, Germany, mistrangelo@iket.fzk.de 2 University
More informationSolution of Partial Differential Equations
Solution of Partial Differential Equations Introduction and Problem Statement Combination of Variables R. Shankar Subramanian We encounter partial differential equations routinely in transport phenomena.
More informationExam in Fluid Mechanics 5C1214
Eam in Fluid Mechanics 5C1214 Final eam in course 5C1214 13/01 2004 09-13 in Q24 Eaminer: Prof. Dan Henningson The point value of each question is given in parenthesis and you need more than 20 points
More informationMAGNETOHYDRODYNAMIC FLOW OF VISCOUS ELECTRICALLY CONDUCTING INCOMPRESSIBLE FLUID THROUGH VERTICAL PLATES SUBJECTED TO INCLINED MAGNETIC FIELD.
International Journal of Scientific and Technical Research in Engineering (IJSTRE) www.ijstre.com Volume 3 Issue 6 ǁ August 2018. MAGNETOHYROYNAMIC FLOW OF VISCOUS ELECTRICALLY CONUCTING INCOMPRESSIBLE
More informationFluid Mechanics Prof. T.I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay. Lecture - 17 Laminar and Turbulent flows
Fluid Mechanics Prof. T.I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay Lecture - 17 Laminar and Turbulent flows Welcome back to the video course on fluid mechanics. In
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationStatus of DCLL Manifold Modeling
Z Status of DCLL Manifold Modeling Presented by Ming Jiu NI Prepared by Ming Jiu NI and Neil B Morley US ITER TBM Meeting, Feb. 14 15, UCLA Contributors from UCLA: Karim Messadek, Sergey Smolentsev, Alice
More informationCFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE
GANIT J. Bangladesh Math. Soc. (ISSN 1606-3694) 31 (2011) 33-41 CFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE Sharmina Hussain Department of Mathematics and Natural Science BRAC University,
More informationREE Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics
REE 307 - Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics 1. Is the following flows physically possible, that is, satisfy the continuity equation? Substitute the expressions for
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More informationUNIT II CONVECTION HEAT TRANSFER
UNIT II CONVECTION HEAT TRANSFER Convection is the mode of heat transfer between a surface and a fluid moving over it. The energy transfer in convection is predominately due to the bulk motion of the fluid
More informationMicroscopic Momentum Balance Equation (Navier-Stokes)
CM3110 Transport I Part I: Fluid Mechanics Microscopic Momentum Balance Equation (Navier-Stokes) Professor Faith Morrison Department of Chemical Engineering Michigan Technological University 1 Microscopic
More informationTurbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing.
Turbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing. Thus, it is very important to form both a conceptual understanding and a quantitative
More information7. Basics of Turbulent Flow Figure 1.
1 7. Basics of Turbulent Flow Whether a flow is laminar or turbulent depends of the relative importance of fluid friction (viscosity) and flow inertia. The ratio of inertial to viscous forces is the Reynolds
More informationHartmann Flow in a Rotating System in the Presence of Inclined Magnetic Field with Hall Effects
Tamkang Journal of Science and Engineering, Vol. 13, No. 3, pp. 243 252 (2010) 243 Hartmann Flow in a Rotating System in the Presence of Inclined Magnetic Field with Hall Effects G. S. Seth, Raj Nandkeolyar*
More informationValidation analyses of advanced turbulence model approaches for stratified two-phase flows
Computational Methods in Multiphase Flow VIII 361 Validation analyses of advanced turbulence model approaches for stratified two-phase flows M. Benz & T. Schulenberg Institute for Nuclear and Energy Technologies,
More informationTransition to turbulence in plane Poiseuille flow
Proceedings of the 55th Israel Annual Conference on Aerospace Sciences, Tel-Aviv & Haifa, Israel, February 25-26, 2015 ThL2T5.1 Transition to turbulence in plane Poiseuille flow F. Roizner, M. Karp and
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationInternational Journal of Mathematical Archive-5(5), 2014, Available online through ISSN
International Journal of Mathematical Archive-5(5), 014, 135-14 Available online through www.ijma.info ISSN 9 5046 STUDY ON LIQUID METAL MHD FLOW THROUGH A TRIANGULAR DUCT UNDER THE ACTION OF VERY STRONG
More informationUNIT IV BOUNDARY LAYER AND FLOW THROUGH PIPES Definition of boundary layer Thickness and classification Displacement and momentum thickness Development of laminar and turbulent flows in circular pipes
More informationCENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer
CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer You are assigned to design a fallingcylinder viscometer to measure the viscosity of Newtonian liquids. A schematic
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationof Friction in Fluids Dept. of Earth & Clim. Sci., SFSU
Summary. Shear is the gradient of velocity in a direction normal to the velocity. In the presence of shear, collisions among molecules in random motion tend to transfer momentum down-shear (from faster
More informationMHD flow and heat transfer due to a linearly stretching sheet. with induced magnetic field: Exact solution. Tarek M. A.
MHD flow and heat transfer due to a linearly stretching sheet with induced magnetic field: Exact solution Tarek M. A. El-Mistikawy Dept. Eng. Math. & Phys., Faculty of Engineering, Cairo University, Giza
More informationMacroscopic plasma description
Macroscopic plasma description Macroscopic plasma theories are fluid theories at different levels single fluid (magnetohydrodynamics MHD) two-fluid (multifluid, separate equations for electron and ion
More informationA semi-implicit finite volume implementation of the CSF method for treating surface tension in interfacial flows
A semi-implicit finite volume implementation of the CSF method for treating surface tension in interfacial flows M. Raessi, M. Bussmann*, and J. Mostaghimi Department of Mechanical and Industrial Engineering,
More informationDrop Impact on a Wet Surface: Computational Investigation of Gravity and Drop Shape
Drop Impact on a Wet Surface: Computational Investigation of Gravity and Drop Shape MURAT DINC and DONALD D. GRAY Department of Civil and Environmental Engineering West Virginia University P.O. Box 6103,
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationWall turbulence with arbitrary mean velocity profiles
Center for Turbulence Research Annual Research Briefs 7 Wall turbulence with arbitrary mean velocity profiles By J. Jiménez. Motivation The original motivation for this work was an attempt to shorten the
More informationShell Balances in Fluid Mechanics
Shell Balances in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University When fluid flow occurs in a single direction everywhere in a system, shell
More informationAbsorption of gas by a falling liquid film
Absorption of gas by a falling liquid film Christoph Albert Dieter Bothe Mathematical Modeling and Analysis Center of Smart Interfaces/ IRTG 1529 Darmstadt University of Technology 4th Japanese-German
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationA surfactant-conserving volume-of-fluid method for interfacial flows with insoluble surfactant
A surfactant-conserving volume-of-fluid method for interfacial flows with insoluble surfactant Ashley J. James Department of Aerospace Engineering and Mechanics, University of Minnesota John Lowengrub
More informationDepartment of Mechanical Engineering
Department of Mechanical Engineering AMEE401 / AUTO400 Aerodynamics Instructor: Marios M. Fyrillas Email: eng.fm@fit.ac.cy HOMEWORK ASSIGNMENT #2 QUESTION 1 Clearly there are two mechanisms responsible
More informationMA3D1 Fluid Dynamics Support Class 5 - Shear Flows and Blunt Bodies
MA3D1 Fluid Dynamics Support Class 5 - Shear Flows and Blunt Bodies 13th February 2015 Jorge Lindley email: J.V.M.Lindley@warwick.ac.uk 1 2D Flows - Shear flows Example 1. Flow over an inclined plane A
More informationExercise 5: Exact Solutions to the Navier-Stokes Equations I
Fluid Mechanics, SG4, HT009 September 5, 009 Exercise 5: Exact Solutions to the Navier-Stokes Equations I Example : Plane Couette Flow Consider the flow of a viscous Newtonian fluid between two parallel
More informationCorresponding Author: Kandie K.Joseph. DOI: / Page
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 13, Issue 5 Ver. 1 (Sep. - Oct. 2017), PP 37-47 www.iosrjournals.org Solution of the Non-Linear Third Order Partial Differential
More informationOscillatory MHD Mixed Convection Boundary Layer Flow of Finite Dimension with Induced Pressure Gradient
Journal of Applied Fluid Mechanics, Vol. 9, No., pp. 75-75, 6. Available online at www.jafmonline.net, ISSN 75-57, EISSN 75-65. DOI:.8869/acadpub.jafm.68.5.876 Oscillatory MHD Mixed Convection Boundary
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationThe distortion observed in the bottom channel of Figure 1 can be predicted from the full transport equation, C t + u C. y D C. z, (1) x D C.
1 8. Shear Dispersion. The transport models and concentration field solutions developed in previous sections assume that currents are spatially uniform, i.e. u f(,y,). However, spatial gradients of velocity,
More information2.3 The Turbulent Flat Plate Boundary Layer
Canonical Turbulent Flows 19 2.3 The Turbulent Flat Plate Boundary Layer The turbulent flat plate boundary layer (BL) is a particular case of the general class of flows known as boundary layer flows. The
More informationChapter 6 Laminar External Flow
Chapter 6 aminar Eternal Flow Contents 1 Thermal Boundary ayer 1 2 Scale analysis 2 2.1 Case 1: δ t > δ (Thermal B.. is larger than the velocity B..) 3 2.2 Case 2: δ t < δ (Thermal B.. is smaller than
More informationG. C. Hazarika 2 Department of Mathematics Dibrugarh University, Dibrugarh
Effects of Variable Viscosity and Thermal Conductivity on Heat and Mass Transfer Flow of Micropolar Fluid along a Vertical Plate in Presence of Magnetic Field Parash Moni Thakur 1 Department of Mathematics
More informationFigure 11.1: A fluid jet extruded where we define the dimensionless groups
11. Fluid Jets 11.1 The shape of a falling fluid jet Consider a circular orifice of a radius a ejecting a flux Q of fluid density ρ and kinematic viscosity ν (see Fig. 11.1). The resulting jet accelerates
More informationLecture 8: Tissue Mechanics
Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials
More informationOpposed Flow Impact on Flame Spread Above Liquid Fuel Pools
Opposed Flow Impact on Flame Spread Above Liquid s Jinsheng Cai, Feng Liu, and William A. Sirignano Department of Mechanical and Aerospace Engineering University of California, Irvine, CA 92697-3975 Abstract
More informationFluid equations, magnetohydrodynamics
Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics
More informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationHomework #4 Solution. μ 1. μ 2
Homework #4 Solution 4.20 in Middleman We have two viscous liquids that are immiscible (e.g. water and oil), layered between two solid surfaces, where the top boundary is translating: y = B y = kb y =
More informationMHD problems in free liquid surfaces as plasma-facing materials in magnetically confined reactors
Fusion Engineering and Design 6/62 (2002) 223/229 www.elsevier.com/locate/fusengdes MHD problems in free liquid surfaces as plasma-facing materials in magnetically confined reactors I. Konkashbaev, A.
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
More informationViscous non-linear theory of Richtmyer-Meshkov Instability. Abstract
Viscous non-linear theory of Richtmyer-Meshkov Instability Pierre Carles and Stéphane Popinet Laboratoire de Modélisation en Mécanique, Université Pierre et Marie Curie, Case 162, 4 place Jussieu, 75252
More informationUnit operations of chemical engineering
1 Unit operations of chemical engineering Fourth year Chemical Engineering Department College of Engineering AL-Qadesyia University Lecturer: 2 3 Syllabus 1) Boundary layer theory 2) Transfer of heat,
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
More information*Corresponding Author: Surajit Dutta, Department of Mathematics, C N B College, Bokakhat, Golaghat, Assam, India
International Journal of Scientific and Innovative Mathematical Research (IJSIMR) Volume 6, Issue, 8, PP -6 ISSN 347-37X (Print) & ISSN 347-34 (Online) DOI: http://dx.doi.org/.43/347-34.6 www.arcjournals.org
More informationNumerical Study of Steady MHD Plane Poiseuille Flow and Heat Transfer in an Inclined Channel
Numerical Study of Steady MHD Plane Poiseuille Flow and Heat Transfer in an Inclined Channel Muhim Chutia Department of Mathematics, Mariani College, Assam-785634, India ABSTRACT: In this paper, a numerical
More informationFluid Dynamics for Ocean and Environmental Engineering Homework #2 Viscous Flow
OCEN 678-600 Fluid Dynamics for Ocean and Environmental Engineering Homework #2 Viscous Flow Date distributed : 9.18.2005 Date due : 9.29.2005 at 5:00 pm Return your solution either in class or in my mail
More informationLaminar Boundary Layers. Answers to problem sheet 1: Navier-Stokes equations
Laminar Boundary Layers Answers to problem sheet 1: Navier-Stokes equations The Navier Stokes equations for d, incompressible flow are + v ρ t + u + v v ρ t + u v + v v = 1 = p + µ u + u = p ρg + µ v +
More informationTurbulent drag reduction by streamwise traveling waves
51st IEEE Conference on Decision and Control December 10-13, 2012. Maui, Hawaii, USA Turbulent drag reduction by streamwise traveling waves Armin Zare, Binh K. Lieu, and Mihailo R. Jovanović Abstract For
More informationMAGNETOHYDRODYNAMICS
Chapter 6 MAGNETOHYDRODYNAMICS 6.1 Introduction Magnetohydrodynamics is a branch of plasma physics dealing with dc or low frequency effects in fully ionized magnetized plasma. In this chapter we will study
More informationSimple examples of MHD equilibria
Department of Physics Seminar. grade: Nuclear engineering Simple examples of MHD equilibria Author: Ingrid Vavtar Mentor: prof. ddr. Tomaž Gyergyek Ljubljana, 017 Summary: In this seminar paper I will
More informationMagnetohydrodynamics Stability of a Compressible Fluid Layer Below a Vacuum Medium
Mechanics and Mechanical Engineering Vol. 12, No. 3 (2008) 267 274 c Technical University of Lodz Magnetohydrodynamics Stability of a Compressible Fluid Layer Below a Vacuum Medium Emad E. Elmahdy Mathematics
More informationCriteria for locally fully developed viscous flow
1 Criteria for locally fully developed viscous flow Ain A. Sonin, MIT October 00 Contents 1. Locally fully developed flow.. Criteria for locally fully developed flow. 3 3. Criteria for constant pressure
More informationInertial Microfluidic Physics
Electronic Supplementary Material (ESI) for Lab on a Chip. This journal is The Royal Society of Chemistry 2014 Inertial Microfluidic Physics Hamed Amini 1,2, Wonhee Lee 3, Dino Di Carlo 1,2* 1. Department
More informationCFD SIMULATIONS OF FLOW, HEAT AND MASS TRANSFER IN THIN-FILM EVAPORATOR
Distillation Absorption 2010 A.B. de Haan, H. Kooijman and A. Górak (Editors) All rights reserved by authors as per DA2010 copyright notice CFD SIMULATIONS OF FLOW, HEAT AND MASS TRANSFER IN THIN-FILM
More informationENTROPY GENERATION IN HEAT AND MASS TRANSFER IN POROUS CAVITY SUBJECTED TO A MAGNETIC FIELD
HEFAT 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 6 8 July Malta ENTROPY GENERATION IN HEAT AND MASS TRANSFER IN POROUS CAVITY SUBJECTED TO A MAGNETIC FIELD Nawaf
More informationROTATING OSCILLATORY MHD POISEUILLE FLOW: AN EXACT SOLUTION
Kragujevac J. Sci. 35 (23) 5-25. UDC 532.527 ROTATING OSCILLATORY MHD POISEUILLE FLOW: AN EXACT SOLUTION Krishan Dev Singh Wexlow Bldg, Lower Kaithu, Shimla-73, India e-mail: kdsinghshimla@gmail.com (Received
More informationFluid Mechanics. du dy
FLUID MECHANICS Technical English - I 1 th week Fluid Mechanics FLUID STATICS FLUID DYNAMICS Fluid Statics or Hydrostatics is the study of fluids at rest. The main equation required for this is Newton's
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationReduction of parasitic currents in the DNS VOF code FS3D
M. Boger a J. Schlottke b C.-D. Munz a B. Weigand b Reduction of parasitic currents in the DNS VOF code FS3D Stuttgart, March 2010 a Institut für Aerodynamik und Gasdynamik, Universität Stuttgart, Pfaffenwaldring
More informationThe Effect Of MHD On Laminar Mixed Convection Of Newtonian Fluid Between Vertical Parallel Plates Channel
The Effect Of MH On Laminar Mixed Convection Of Newtonian Fluid Between Vertical Parallel Plates Channel Rasul alizadeh,alireza darvish behanbar epartment of Mechanic, Faculty of Engineering Science &
More informationMYcsvtu Notes HEAT TRANSFER BY CONVECTION
www.mycsvtunotes.in HEAT TRANSFER BY CONVECTION CONDUCTION Mechanism of heat transfer through a solid or fluid in the absence any fluid motion. CONVECTION Mechanism of heat transfer through a fluid in
More informationNumerical Simulation of Core- Annular Flow in a Curved Pipe
Numerical Simulation of Core- Annular Flow in a Curved Pipe Simulation of two phase flow with OpenFOAM Report Number:2625(MEAH:277) Master of Science Thesis Process and Energy Numerical Simulation of
More informationDetailed 3D modelling of mass transfer processes in two phase flows with dynamic interfaces
Detailed 3D modelling of mass transfer processes in two phase flows with dynamic interfaces D. Darmana, N.G. Deen, J.A.M. Kuipers Fundamentals of Chemical Reaction Engineering, Faculty of Science and Technology,
More informationReynolds stress analysis of EMHD-controlled wall turbulence. Part I. Streamwise forcing
Reynolds stress analysis of EMHD-controlled wall turbulence. Part I. Streamwise forcing Catherine H. Crawford, and George Em Karniadakis Citation: Physics of Fluids 9, 788 (1997); doi: 10.1063/1.869210
More informationMomentum (Newton s 2nd Law of Motion)
Dr. Nikos J. Mourtos AE 160 / ME 111 Momentum (Newton s nd Law of Motion) Case 3 Airfoil Drag A very important application of Momentum in aerodynamics and hydrodynamics is the calculation of the drag of
More informationGENERALISATION OF THE TWO-SCALE MOMENTUM THEORY FOR COUPLED WIND TURBINE/FARM OPTIMISATION
25 th National Symposium on Wind Engineering, Tokyo, Japan, 3-5 December 2018 第 25 回風工学シンポジウム (2018) GENERALISATION OF THE TWO-SCALE MOMENTUM THEORY FOR COUPLED WIND TURBINE/FARM OPTIMISATION Takafumi
More informationSelf-Organization of Plasmas with Flows
Self-Organization of Plasmas with Flows ICNSP 2003/ 9/10 Graduate School of Frontier Sciences,, National Institute for Fusion Science R. NUMATA, Z. YOSHIDA, T. HAYASHI ICNSP 2003/ 9/10 p.1/14 Abstract
More informationNumerical simulation of wave breaking in turbulent two-phase Couette flow
Center for Turbulence Research Annual Research Briefs 2012 171 Numerical simulation of wave breaking in turbulent two-phase Couette flow By D. Kim, A. Mani AND P. Moin 1. Motivation and objectives When
More informationStorm Water Best Management Practice: Development of Debris Filtering Structure for Supercritical Flow
Storm Water Best Management Practice: Development of Debris Filtering Structure for Supercritical Flow Jungseok Ho 1, Todd Marti 2, and Julie Coonrod 3 1 Department of Civil Engineering, University of
More informationBasic Fluid Mechanics
Basic Fluid Mechanics Chapter 6A: Internal Incompressible Viscous Flow 4/16/2018 C6A: Internal Incompressible Viscous Flow 1 6.1 Introduction For the present chapter we will limit our study to incompressible
More information