Falling liquid films: wave patterns and thermocapillary effects
|
|
- Kelly Eustace Cobb
- 6 years ago
- Views:
Transcription
1 Falling liquid films: wave patterns and termocapillary effects Benoit Sceid Cimie-Pysique E.P., Université Libre de Bruxelles C.P. 65/6 Avenue F.D. Roosevelt, Bruxelles - BELGIUM bsceid@ulb.ac.be Abstract Tree-dimensional waves in film flows down a uniformly eated plane are investigated by simulation. For small Reynolds number, regularly spaced rivulets are observed, aligned wit te flow and sustaining two-dimensional waves of larger amplitude and pase speed tan in isotermal conditions. For larger Reynolds number, te picture is similar to te isotermal case and no rivulets appear. Te transition between tese two regimes sows complex cooperative beavior between bot ydrodynamic and termocapillary modes. Key words: Falling film, Hydrodynamic instability, Termocapillary effect Introduction Tree-dimensional (3D) ydrodynamic waves in falling films ave been investigated experimentally by various autors (see for instance [). Experiments by Liu & Gollub [4 and more recently by Park & Nosoko [5 give te clearest picture of te wide penomenology of te interacting 3D waves on isotermal film flows suc as syncronously deformed fronts, subarmonic patterns and orsesoe-like waves. Teoretically, we ave applied a regularization procedure to te weigted integral boundary layer metod providing a two-dimensional modeling tat describe te complex wave patterns observed in experiments up to moderate Reynolds numbers [8. For eated walls, Joo et al. [3 ave modeled rivulet patterns induced by termocapillary effects for small Reynolds number, i.e. using a single evolution equation for te film tickness. However teir simulations experienced finite-time blow-up tat was sown to be intrinsic to te Benney-type equation. Indeed, Ramaswamy et al. [6 ave simulated te full tree-dimensional boundary layer equations and could resolve te wole dynamics of a rivulet formation, from onset to rupture. Neverteless, tey were restricted to a domain size of about one wavelengt in bot streamwise and spanwise directions. We aim in tis paper to propose a two-dimensional model describing te tree-dimensional dynamics of a eated falling film in large-scale domains and to investigate te influence of te Marangoni effect on wave patterns observed in isotermal conditions up to moderate Reynolds number. We will in fact extend to 3D case te model developed by Ruyer-Quil et al. [7 for two-dimensional waves.
2 Matematical formulation of te problem Te flow of a Newtonian liquid down a eated vertical plane is considered. We look for D equations in te streamwise (x) and spanwise (z) coordinates tat mimic te full 3D motion of te fluid. Including te eating of te plane into te regularized tree-field model as detailed in [8 leads to a four-field model independent on te cross-stream coordinate (y), for te film tickness, te streamwise and spanwise flow rates q and p and te temperature field at te interface θ: t = x q z p, (a) [ 9 q δ t q = δ 7 x 7 { [ q 7 xq + δ 8 q z p 7 9 p z q qp z q [ +η 4 q ( x) 9 x q x 3 z p x 43 x p z 73 p xz 6 q xx [ 9 δ t p = δ 7 [ zq z + 3 q ( z ) + 3 p x z 3 q zz xxq + zz q + 7 xzp ( 5 M 4 xθ δ ) 4 q xxθ + 5 } 6 ( xxx + xzz ) ( δ 70 q x + M 5 ) x θ, (b) 56 +η p z 7 p 7 zp 8 7 xp x 4 p ( z) z p z p x q x q z 6 q x p p ( x ) + 3 q x z qp x 5 p 43 z q x 73 q xz 6 p zz 6 6 p xx 5 4 M zθ ( xxz + zzz ), [ ( θ Bθ) 7 δ t θ = 3 + δ Pr 40 ( θ)( xq + z p) 7 + η [( θ 3 )( Pr Bθ ( x ) + ( ) z) +( θ) ( xx + zz ) were te reduced dimensionless numbers are δ = (3Re)/9 Γ /3, B = Bi(3Re) /3, M = + 9 zzp + xx p + 7 xzq (q x θ + p z θ) + x x θ + z z θ (c) + xx θ + zz θ, (d) Ma Γ /3 (3Re) 4/9 and η = (3Re)4/9 Γ /3,
3 based in turn on te usual Reynolds, Biot, Marangoni and Kapitza numbers: Re = g3 N 3ν, αν/3 γ T Bi =, Ma = and Γ = λg/3 ρν 4/3 g /3 σ ρν 4/3 g /3, wit g te gravity acceleration, N te tickness of te uniform film, ν te kinematic viscosity, ρ te density, σ te surface tension, γ = dσ/dt, α te eat transfer coefficient and λ te termal conductivity; Pr = ν/χ is te Prandtl number were χ is te termal diffusivity. Te lengtscales are N and N / η in te wall-normal y and in te in-plane x and z directions, respectively, wile te timescale is ν/g N η. Equation (a) is te mass conservation equation, equations (b,c) express te averaged momentum balances in bot directions x and z, and equation (d) is te energy balance. 3 Computations 3. Small-size domain We first validate our model wit a small-scale simulation investigated by Ramaswamy et al. [6 for small Reynolds numbers, i.e. wen a one-field model can still be valid. Tey demonstrated a mecanism of rivulet formation based solely on instability penomena. Figure sows simulation results, done wit periodic boundary conditions and a simple armonic disturbance of te form (x, z, 0) = + 0. cos(k x x)+ 0. cos(k z z), as initial condition. Te wavenumbers k x and k z ave been cosen to be below tose for te maximal linear growt rate in eac direction, wic are k xmax =0.56 and k z max =0.53. Actually, tese appropriate values allow for interesting secondary flow development. Te initial perturbation creates a troug in te center of te domain (a). Ten, termocapillarity sets in, displacing te fluid from te otter trougs towards te surrounding colder crests. However, te growt rate of te ydrodynamic mode is dominant at initial stage and surface waves develop (b). Te local pase speed being proportional to te square of te local film tickness, te crest travel faster tan te troug (c). In te absence of te mean flow in te spanwise direction, te liquid is displaced laterally due to termocapillarity. Terefore, as time progresses, te tinning of te liquid layer persists and forms a valley surrounded by rivulets aligned wit te flow (d). Tis process is similar to te evolution of a eated tin film on a orizontal substrate. Likewise, it exibits te formation of a secondary rivulet between te main one (e). As found by Boos & Tess [ for orizontal layers, a true cascade of structures takes place in tinner zones (f), precursor to te film rupture. Te last stage before rupture obtained by Ramaswamy et al. wit DNS is at t = 53, wic is in excellent agreement wit Fig.e. However, our simulation runs beyond tis time and exibits finer structures in te film before its rupture (see Fig.g). Te reason wy te full-scale computations wit Navier-Stokes/Fourier equations did not provide suc furter evolution 3
4 (a) t =0 (b) t =5 (c) t = (d) t =90 (e) t = 53 (f) t = 75 Figure : Inception and development of rivulet aligned wit te flow computed wit () arising for Re =/3, Γ = 0, Ma =, Bi = and Pr = 7. Te domain is square of side π/k x wit k x = k z = Te flow direction is from te rigt top towards te left bottom. of te film is likely related to te coice of te number of mes points used in [6. Te autors would ave most probably been able to compute larger time wit refined grid resolution. However, tis would ave been at te expense of computing time, wic demonstrate te great advantage of working wit a model of reduced dimensionality. Indeed, our reduced model involves variables tat are defined only at te film surface and not in te bulk. Actually, te mes points are 3 3 and te time necessary for computing te case of Fig., wit an accuracy of 4, is about one our on a standard desktop computer. Now, let increase te Reynolds number up to Re =, i.e. out of te range of validity of a single evolution equation for te film tickness. Te wavenumbers for te maximal linear growt rate are k xmax 0.6 and k z max 0.34 (approximated values calculated in te long-wave limit), wic yields k x k xmax / wile k z k zmax. Figure sows tat te ydrodynamic mode quickly generates ig amplitude waves (a,b) tat evolves to a solitary-like wave wit preceding capillary ripples (c). However, te termocapillary mode leads te film to rupture before te surface wave saturates. An interesting interaction between te two modes can be pointed out ere: as te termocapillary flow feeds te core of te rivulet, te mean film tickness at te crest increases, and so te local flow rate. Hence, te wave solution 4
5 (a) t =40 (b) t = 400 (c) t = 690 Figure : Same as for Fig. wit Re =. does not saturate and rater follows te cange of te local Reynolds number, by increasing its amplitude and its pase speed. Tis process vanises at t>6 wen te film approaces te wall enoug for te viscous stress to slow down te lateral termocapillary flow. Te ydrodynamic wave and te longitudinal rivulet were found to coexist along a distance tat is about 50 wavelengts ere and tat increases wit Re. 3. Large-size domain Let us turn now to large-scale simulations of a water film at C (Γ = 3375, Pr =7) wit a temperature difference between te vertical wall and te ambient air of 5 C (Ma = 50) and a ig eat transfer coefficient of 00W/cm K(Bi =0.). Te simulations tat follow were started wit a wite noise of amplitude /00 of te flat film tickness. Figure 3 sows again te formation of rivulets due to Marangoni effect: after te development of a parallel wave train (a), drop-like accumulation breaks te D waves into 3D patterns (b,c), precursor of rivulet patterns aligned wit te flow (d,e). As depicted previously for Fig., te liquid ten accumulates into rivulets, wic increases te local Reynolds number and fosters D solitarylike waves of larger amplitude and pase speed tan in isotermal conditions (f). Tis process continues until te rupture of te film, te snapsot of wic is sown in Fig.4a. Similar rivulet pattern, but wit larger wavelengt, is also sown for Re = 5 (b) wile no rivulet develops for still larger Re = (c) were te wave pattern beaves like in isotermal conditions [5. Figure 5 sows finally te pase speed of solitary waves versus Re. For small Reynolds number and small pase speed, i.e. in te drag-gravity regime, inertia effects are small as compared to termocapillary effects wile it is te contrary for large Reynolds number and large pase speed, i.e. in te drag-inertia regime, were inertia is dominant. Te resulting patterns are drastically different (a,c). Te transition between tese two regimes for 4 <Re<6 sows complex cooperative beavior between bot ydrodynamic and termocapillary modes, mostly wen teir 5
6 (a) 0 (b) 500 (c) 700 (d) 900 (e) 0 (f) 500 Figure 3: Water film free surface at different time for Re = for Ma = 50, Bi = 0., Pr = 7 and Γ = Te domain is a square of side π/k x were k x = k z =0.05. Brigt (resp. dark) zones correspond to elevations (resp. depressions). amplitudes are comparable, as illustrated for Re = 5 in Fig.4b. 4 Discussion Te dynamical coupled equations for te film tickness, te surface temperature θ and te streamwise and spanwise flow rates q and p given in () as been sown to be robust and accurate in describing te competition between ydrodynamic waves and termocapillary effect responsible for te rivulet formation in eated falling liquid films. For constant temperature difference across te film, it as been found tat in te drag-gravity regime, quasi-regularly spaced rivulets arise and grow up until rupture (a). Meanwile, te rivulet confines te flow in suc a way tat waves beaves like two-dimensional solitary waves, but of iger flow rate because of te local increase of te Reynolds number. On te contrary, no qualitative influence of te Marangoni effect as been observed in te drag-inertia regime, at least during te 6
7 (a) Re = (b) Re =5 (c) Re = Figure 4: Wave patterns for Ma = 50, Bi =0., Pr = 7 and Γ = (b) (c) c.5 (a) Re Figure 5: Pase speed c of solitary waves versus Re for Ma = 50 (solid line) and Ma = 0 (dased line). Te letters refer to Fig.4. time of te simulations, sowing tat inertia fully dominates te dynamics of te film (c). Finally te transition between tese two regimes is also instructive in terms of rupture wic occurs at only one place and not between all rivulets at te same time. Tis is due to te dispersion of te rivulet size and amplitude, allowing also to observe different kind of waves: D wave trains, D solitary waves and 3D wave structure (b). Having at one s disposal reliable low-dimension model, one is now able to undergo systematic numerical analyzes at low cost (in terms of CPU time) and for large domains. Tis will allow to explore in te wole parameter space te complex 3D interactions between ydrodynamic and termocapillary effects. Notice owever tat, to our knowledge, no experimental data of film flowing down uniformly eated substrates are yet available for comparison purpose. 7
8 Acknowledgment Te autor would like to tank Serafim Kalliadasis and Cristian Ruyer-Quil for fruitful discussions. References [ S.V. Alekseenko, V.E. Nakoryakov, and B.G. Pokusaev. Wave flow in liquid films. Begell House (New York), tird edition, 994. [ W. Boos and A. Tess. Cascade of structures in long-wavelengt marangoni instability. Pys. Fluids, :633, 999. [3 S.W. Joo and S.H. Davis. Instabilities of tree-dimensional teory viscous falling films. J. Fluid Mec., 4:59, 99. [4 J. Liu, J.B. Scneider, and J.P. Gollub. Tree-dimensional instabilities of film flows. Pys. Fluids, 7:55 67, 995. [5 C. D. Park and T. Nosoko. Tree-dimensional wave dynamics on a falling film and associated mass transfer. AICE J., 49():75 77, 03. [6 B. Ramaswamy, S. Krisnamoorty, and S.W. Joo. Tree-dimensional simulation of instabilities and rivulet formation in eated falling films. J. Comp. Pys., 3:70 88, 997. [7 C. Ruyer-Quil, B. Sceid, S. Kalliadasis, M.G. Velarde, and R.K. Zeytounian. Termocapillary long waves in a liquid film flow. part. low dimensional formulation. J. Fluid Mec., 538:99, 05. [8 B. Sceid, C. Ruyer-Quil, and P. Manneville. Wave patterns in film flows: Modelling and tree-dimensional waves. J. Fluid Mec., to appear, 06. 8
Thermocapillary long waves in a liquid film flow. Part 2. Linear stability and nonlinear waves
J. Fluid Mec. (5), vol. 58, pp. 44. c 5 Cambridge University Press doi:.7/s55446 Printed in te United Kingdom Termocapillary long waves in a liquid film flow. Part. Linear stability and nonlinear waves
More informationThin films flowing down inverted substrates: Two dimensional flow
PHYSICS OF FLUIDS, 55 Tin films flowing down inverted substrates: Two dimensional flow Te-Seng Lin and Lou Kondic Department of Matematical Sciences and Center for Applied Matematics and Statistics, New
More informationPrediction of Coating Thickness
Prediction of Coating Tickness Jon D. Wind Surface Penomena CE 385M 4 May 1 Introduction Tis project involves te modeling of te coating of metal plates wit a viscous liquid by pulling te plate vertically
More informationEPJ E. Soft Matter and Biological Physics. EPJ.org. Thin viscoelastic dewetting films of Jeffreys type subjected to gravity and substrate interactions
EPJ E Soft Matter and Biological Pysics EPJ.org your pysics journal Eur. Pys. J. E (9) : DOI./epje/i9-77- Tin viscoelastic dewetting films of Jeffreys type subjected to gravity and substrate interactions
More informationREVIEW LAB ANSWER KEY
REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g
More informationChapters 19 & 20 Heat and the First Law of Thermodynamics
Capters 19 & 20 Heat and te First Law of Termodynamics Te Zerot Law of Termodynamics Te First Law of Termodynamics Termal Processes Te Second Law of Termodynamics Heat Engines and te Carnot Cycle Refrigerators,
More informationHT TURBULENT NATURAL CONVECTION IN A DIFFERENTIALLY HEATED VERTICAL CHANNEL. Proceedings of 2008 ASME Summer Heat Transfer Conference HT2008
Proceedings of 2008 ASME Summer Heat Transfer Conference HT2008 August 10-14, 2008, Jacksonville, Florida USA Proceedings of HT2008 2008 ASME Summer Heat Transfer Conference August 10-14, 2008, Jacksonville,
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationSimulation and verification of a plate heat exchanger with a built-in tap water accumulator
Simulation and verification of a plate eat excanger wit a built-in tap water accumulator Anders Eriksson Abstract In order to test and verify a compact brazed eat excanger (CBE wit a built-in accumulation
More informationGrade: 11 International Physics Olympiad Qualifier Set: 2
Grade: 11 International Pysics Olympiad Qualifier Set: 2 --------------------------------------------------------------------------------------------------------------- Max Marks: 60 Test ID: 12111 Time
More informationDesalination by vacuum membrane distillation: sensitivity analysis
Separation and Purification Tecnology 33 (2003) 75/87 www.elsevier.com/locate/seppur Desalination by vacuum membrane distillation: sensitivity analysis Fawzi Banat *, Fami Abu Al-Rub, Kalid Bani-Melem
More informationA Numerical Scheme for Particle-Laden Thin Film Flow in Two Dimensions
A Numerical Sceme for Particle-Laden Tin Film Flow in Two Dimensions Mattew R. Mata a,, Andrea L. Bertozzi a a Department of Matematics, University of California Los Angeles, 520 Portola Plaza, Los Angeles,
More informationOn the absence of marginal pinching in thin free films
On te absence of marginal pincing in tin free films P. D. Howell and H. A. Stone 6 August 00 Abstract Tis paper concerns te drainage of a tin liquid lamella into a Plateau border. Many models for draining
More information3D Large scale Marangoni convection in liquid films
Eur. Pys. J. B 33, 457 467 (2003) DOI: 0.40/epjb/e2003-0086-3 THE EUROPEAN PHYSICAL JOURNAL B 3D Large scale Marangoni convection in liquid films M. Besteorn,a,A.Pototsky, and U. Tiele 2 Lerstul für Teoretisce
More informationFlow of a Rarefied Gas between Parallel and Almost Parallel Plates
Flow of a Rarefied Gas between Parallel and Almost Parallel Plates Carlo Cercignani, Maria Lampis and Silvia Lorenzani Dipartimento di Matematica, Politecnico di Milano, Milano, Italy 033 Abstract. Rarefied
More informationEffects of Radiation on Unsteady Couette Flow between Two Vertical Parallel Plates with Ramped Wall Temperature
Volume 39 No. February 01 Effects of Radiation on Unsteady Couette Flow between Two Vertical Parallel Plates wit Ramped Wall Temperature S. Das Department of Matematics University of Gour Banga Malda 73
More informationLarge eddy simulation of turbulent flow downstream of a backward-facing step
Available online at www.sciencedirect.com Procedia Engineering 31 (01) 16 International Conference on Advances in Computational Modeling and Simulation Large eddy simulation of turbulent flow downstream
More informationWhy gravity is not an entropic force
Wy gravity is not an entropic force San Gao Unit for History and Pilosopy of Science & Centre for Time, SOPHI, University of Sydney Email: sgao7319@uni.sydney.edu.au Te remarkable connections between gravity
More informationExcluded Volume Effects in Gene Stretching. Pui-Man Lam Physics Department, Southern University Baton Rouge, Louisiana
Excluded Volume Effects in Gene Stretcing Pui-Man Lam Pysics Department, Soutern University Baton Rouge, Louisiana 7083 Abstract We investigate te effects excluded volume on te stretcing of a single DNA
More information6. Non-uniform bending
. Non-uniform bending Introduction Definition A non-uniform bending is te case were te cross-section is not only bent but also seared. It is known from te statics tat in suc a case, te bending moment in
More informationUNIVERSITY OF SASKATCHEWAN Department of Physics and Engineering Physics
UNIVERSITY O SASKATCHEWAN Department of Pysics and Engineering Pysics Pysics 117.3 MIDTERM EXAM Regular Sitting NAME: (Last) Please Print (Given) Time: 90 minutes STUDENT NO.: LECTURE SECTION (please ceck):
More informationNotes: Most of the material in this chapter is taken from Young and Freedman, Chap. 12.
Capter 6. Fluid Mecanics Notes: Most of te material in tis capter is taken from Young and Freedman, Cap. 12. 6.1 Fluid Statics Fluids, i.e., substances tat can flow, are te subjects of tis capter. But
More informationShear Thinning and Hydrodynamic Friction of Viscosity Modifier- Containing Oils. Part II: Impact of Shear Thinning on Journal Bearing Friction
Tribology Letters 208 66:9 ttps://doi.org/0.007/s249-08-040-z ORIGINAL PAPER Sear Tinning and Hydrodynamic Friction of Viscosity Modifier- Containing Oils. Part II: Impact of Sear Tinning on Journal Bearing
More informationFlavius Guiaş. X(t + h) = X(t) + F (X(s)) ds.
Numerical solvers for large systems of ordinary differential equations based on te stocastic direct simulation metod improved by te and Runge Kutta principles Flavius Guiaş Abstract We present a numerical
More informationNumerical analysis of a free piston problem
MATHEMATICAL COMMUNICATIONS 573 Mat. Commun., Vol. 15, No. 2, pp. 573-585 (2010) Numerical analysis of a free piston problem Boris Mua 1 and Zvonimir Tutek 1, 1 Department of Matematics, University of
More informationReynolds number effects in stratified turbulent wakes
Abstract Reynolds number effects in stratified turbulent wakes Qi Zou 1 and Peter Diamessis 2 1 Department of Applied Matematics and Teoretical Pysics, University of Cambridge 2 Scool of Civil and Environmental
More informationThe Basics of Vacuum Technology
Te Basics of Vacuum Tecnology Grolik Benno, Kopp Joacim January 2, 2003 Basics Many scientific and industrial processes are so sensitive tat is is necessary to omit te disturbing influence of air. For
More informationSIMG Solution Set #5
SIMG-303-0033 Solution Set #5. Describe completely te state of polarization of eac of te following waves: (a) E [z,t] =ˆxE 0 cos [k 0 z ω 0 t] ŷe 0 cos [k 0 z ω 0 t] Bot components are traveling down te
More informationFabric Evolution and Its Effect on Strain Localization in Sand
Fabric Evolution and Its Effect on Strain Localization in Sand Ziwei Gao and Jidong Zao Abstract Fabric anisotropy affects importantly te overall beaviour of sand including its strengt and deformation
More informationHall Effcts Eon Unsteady MHD Free Convection Flow Over A Stretching Sheet With Variable Viscosity And Viscous Dissipation
IOSR Journal of Matematics (IOSR-JM) e-issn: 78-578, p-issn: 39-765X. Volume, Issue 4 Ver. I (Jul - Aug. 5), PP 59-67 www.iosrjournals.org Hall Effcts Eon Unsteady MHD Free Convection Flow Over A Stretcing
More informationHydrodynamics of Bounded Vertical Film with Nonlinear Surface Properties
Under consideration for publication in Journal of Colloids and Interface Science Hydrodynamics of Bounded Vertical Film wit Nonlinear Surface Properties A. Homayoun Heidari 1,, Ricard J. Braun,, Amir H.
More informationCFD calculation of convective heat transfer coefficients and validation Part I: Laminar flow Neale, A.; Derome, D.; Blocken, B.; Carmeliet, J.E.
CFD calculation of convective eat transfer coefficients and validation Part I: Laminar flow Neale, A.; Derome, D.; Blocken, B.; Carmeliet, J.E. Publised in: IEA Annex 41 working meeting, Kyoto, Japan Publised:
More informationarxiv: v1 [physics.flu-dyn] 3 Jun 2015
A Convective-like Energy-Stable Open Boundary Condition for Simulations of Incompressible Flows arxiv:156.132v1 [pysics.flu-dyn] 3 Jun 215 S. Dong Center for Computational & Applied Matematics Department
More informationTHE EUROPEAN PHYSICAL JOURNAL E. Modelling thin-film dewetting on structured substrates and templates: Bifurcation analysis and numerical simulations
ur. Pys. J., 7 (3) DOI./epje/i3-9- TH UROPAN PHYSICAL JOURNAL Modelling tin-film dewetting on structured substrates and templates: Bifurcation analysis and numerical simulations U. Tiele,a,L.Brusc,M.Besteorn,andM.Bär
More informationEffect of surface tension on swell-induced surface instability of substrate-confined hydrogel layers
PAPER CREATED USING THE RSC ARTICLE TEMPLATE (VER. 3.) - SEE WWW.RSC.ORG/ELECTRONICFILES FOR DETAILS www.rsc.org/xxxxxx XXXXXXXX 5 5 Effect of surface tension on swell-induced surface instability of substrate-confined
More informationInclined flow of a heated fluid film with temperature dependent fluid properties
Inclined flow of a eated fluid film wit temperature dependent fluid properties N. Gonput, J.P. Pascal and S.J.D. D Alessio Department of Matematics, Ryerson University, Toronto, Ontario, M5B K, Canada
More informationA thin-film equation for viscoelastic liquids of Jeffreys type
Eur. Pys. J. E 7, 373 379 5) DOI.4/epje/i5-6-8 THE EUROPEAN PHYSICAL JOURNAL E A tin-film equation for viscoelastic liquids of Jeffreys type M. Rauscer,, A. Münc 3, B. Wagner 4, and R. Blossey 5,a Max-Planck-Institut
More informationAN EFFICIENT AND ROBUST METHOD FOR SIMULATING TWO-PHASE GEL DYNAMICS
AN EFFICIENT AND ROBUST METHOD FOR SIMULATING TWO-PHASE GEL DYNAMICS GRADY B. WRIGHT, ROBERT D. GUY, AND AARON L. FOGELSON Abstract. We develop a computational metod for simulating models of gel dynamics
More informationStability of stratified two-phase channel flows of Newtonian/non-Newtonian shear-thinning fluids
Stability of stratified two-pase cannel flows of Newtonian/non-Newtonian sear-tinning fluids D. Picci a), I. Barmak b), A. Ullmann, and N. Brauner Scool of Mecanical Engineering, Tel Aviv University, Tel
More informationAPPROXIMATION OF CRYSTALLINE DENDRITE GROWTH IN TWO SPACE DIMENSIONS. Introduction
Acta Mat. Univ. Comenianae Vol. LXVII, 1(1998), pp. 57 68 57 APPROXIMATION OF CRYSTALLINE DENDRITE GROWTH IN TWO SPACE DIMENSIONS A. SCHMIDT Abstract. Te pase transition between solid and liquid in an
More informationComment on Experimental observations of saltwater up-coning
1 Comment on Experimental observations of saltwater up-coning H. Zang 1,, D.A. Barry 2 and G.C. Hocking 3 1 Griffit Scool of Engineering, Griffit University, Gold Coast Campus, QLD 4222, Australia. Tel.:
More informationContinuity and Differentiability Worksheet
Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationGeneviève BELANGÉR, Kristjan KANNIKE, Alexander PUKHOV and Martti RAIDAL
Geneviève BELANGÉR,, Alexander PUKHOV and Martti RAIDAL cuola Normale uperiore di Pisa E-mail: belanger@lapt.cnrs.fr, kristjan.kannike@cern.c, pukov@lapt.cnrs.fr, martti.raidal@cern.c calar dark matter
More informationMVT and Rolle s Theorem
AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationFilm thickness Hydrodynamic pressure Liquid saturation pressure or dissolved gases saturation pressure. dy. Mass flow rates per unit length
NOTES DERITION OF THE CLSSICL REYNOLDS EQTION FOR THIN FIL FLOWS Te lecture presents te derivation of te Renolds equation of classical lubrication teor. Consider a liquid flowing troug a tin film region
More informationThe entransy dissipation minimization principle under given heat duty and heat transfer area conditions
Article Engineering Termopysics July 2011 Vol.56 No.19: 2071 2076 doi: 10.1007/s11434-010-4189-x SPECIAL TOPICS: Te entransy dissipation minimization principle under given eat duty and eat transfer area
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationFirooz Arash (a,b) , which have caused confusion in the normalization of ZEUS data is discussed and resolved. PCAC I. INTRODUCTION
PION STRUCTURE FUNCTION F π 2 In te Valon Model Firooz Aras (a,b) (a) Center for Teoretical Pysics and Matematics, AEOI P.O.Box 365-8486, Teran, Iran (b) Pysics Department, AmirKabir University (Tafres
More informationUndulation instability in two-dimensional foams of magnetic fluid
Eur. Pys. J. B 3, 3 9 (998) THE EUROPEAN PHYSICAL JOURNAL B c EDP Sciences Springer-Verlag 998 Undulation instability in two-dimensional foams of magnetic fluid F. Elias,a,I.Drikis,, A. Cebers,C.Flament,andJ.-C.Bacri
More informationOptimal Shape Design of a Two-dimensional Asymmetric Diffsuer in Turbulent Flow
THE 5 TH ASIAN COMPUTAITIONAL FLUID DYNAMICS BUSAN, KOREA, OCTOBER 7 ~ OCTOBER 30, 003 Optimal Sape Design of a Two-dimensional Asymmetric Diffsuer in Turbulent Flow Seokyun Lim and Haeceon Coi. Center
More informationFinding and Using Derivative The shortcuts
Calculus 1 Lia Vas Finding and Using Derivative Te sortcuts We ave seen tat te formula f f(x+) f(x) (x) = lim 0 is manageable for relatively simple functions like a linear or quadratic. For more complex
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationhal , version 1-29 Jan 2008
ENERGY OF INTERACTION BETWEEN SOLID SURFACES AND LIQUIDS Henri GOUIN L. M. M. T. Box 3, University of Aix-Marseille Avenue Escadrille Normandie-Niemen, 3397 Marseille Cedex 0 France E-mail: enri.gouin@univ-cezanne.fr
More informationCFD calculation of convective heat transfer coefficients and validation Part I: Laminar flow. Annex 41 Kyoto, April 3 rd to 5 th, 2006
CFD calculation of convective eat transfer coefficients and validation Part I: Laminar flow Annex 41 Kyoto, April 3 rd to 5 t, 2006 Adam Neale 1, Dominique Derome 1, Bert Blocken 2 and Jan Carmeliet 2,3
More informationChapter 3 Thermoelectric Coolers
3- Capter 3 ermoelectric Coolers Contents Capter 3 ermoelectric Coolers... 3- Contents... 3-3. deal Equations... 3-3. Maximum Parameters... 3-7 3.3 Normalized Parameters... 3-8 Example 3. ermoelectric
More informationStudy on Dynamic Wetting Behavior in Microgravity Condition Targeted for Propellant Tank
Int. J. Microgravity Sci. Appl., 34 (3) (2017) 340306 DOI: 10.15011//jasma.34.340306 IIIII AMS2016 Proceedings IIIII (Original Paper) Study on Dynamic Wetting Beavior in Microgravity Condition Targeted
More informationJian-Guo Liu 1 and Chi-Wang Shu 2
Journal of Computational Pysics 60, 577 596 (000) doi:0.006/jcp.000.6475, available online at ttp://www.idealibrary.com on Jian-Guo Liu and Ci-Wang Su Institute for Pysical Science and Tecnology and Department
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationINTRODUCTION DEFINITION OF FLUID. U p F FLUID IS A SUBSTANCE THAT CAN NOT SUPPORT SHEAR FORCES OF ANY MAGNITUDE WITHOUT CONTINUOUS DEFORMATION
INTRODUCTION DEFINITION OF FLUID plate solid F at t = 0 t > 0 = F/A plate U p F fluid t 0 t 1 t 2 t 3 FLUID IS A SUBSTANCE THAT CAN NOT SUPPORT SHEAR FORCES OF ANY MAGNITUDE WITHOUT CONTINUOUS DEFORMATION
More informationQuantum Theory of the Atomic Nucleus
G. Gamow, ZP, 51, 204 1928 Quantum Teory of te tomic Nucleus G. Gamow (Received 1928) It as often been suggested tat non Coulomb attractive forces play a very important role inside atomic nuclei. We can
More informationA Reconsideration of Matter Waves
A Reconsideration of Matter Waves by Roger Ellman Abstract Matter waves were discovered in te early 20t century from teir wavelengt, predicted by DeBroglie, Planck's constant divided by te particle's momentum,
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More information3.1 Extreme Values of a Function
.1 Etreme Values of a Function Section.1 Notes Page 1 One application of te derivative is finding minimum and maimum values off a grap. In precalculus we were only able to do tis wit quadratics by find
More informationShrinkage anisotropy characteristics from soil structure and initial sample/layer size. V.Y. Chertkov*
Srinkage anisotropy caracteristics from soil structure and initial sample/layer size V.Y. Certkov* Division of Environmental, Water, and Agricultural Engineering, Faculty of Civil and Environmental Engineering,
More informationTHE MODEL OF DRYING SESSILE DROP OF COLLOIDAL SOLUTION
THE MODEL OF DRYING SESSILE DROP OF COLLOIDAL SOLUTION I.V.Vodolazskaya, Yu.Yu.Tarasevic Astrakan State University, a Tatiscev St., Astrakan, 41456, Russia e-mail: tarasevic@aspu.ru e ave proposed and
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationModels of non-boussinesq lock-exchange flow
Under consideration for publication in J. Fluid Mec. Models of non-boussinesq lock-excange flow R. ROTUNNO, J. B. KLEMP, G. H. BRYAN AND D. J. MURAKI 2 National Center for Atmosperic Researc, Boulder,
More informationSome Review Problems for First Midterm Mathematics 1300, Calculus 1
Some Review Problems for First Midterm Matematics 00, Calculus. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd,
More informationChapter 2 Ising Model for Ferromagnetism
Capter Ising Model for Ferromagnetism Abstract Tis capter presents te Ising model for ferromagnetism, wic is a standard simple model of a pase transition. Using te approximation of mean-field teory, te
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationDistribution of reynolds shear stress in steady and unsteady flows
University of Wollongong Researc Online Faculty of Engineering and Information Sciences - Papers: Part A Faculty of Engineering and Information Sciences 13 Distribution of reynolds sear stress in steady
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationLecture 10: Carnot theorem
ecture 0: Carnot teorem Feb 7, 005 Equivalence of Kelvin and Clausius formulations ast time we learned tat te Second aw can be formulated in two ways. e Kelvin formulation: No process is possible wose
More informationINVERSE TEMPERATURE RECONSTRUCTION IN THERMOCAPILLARY-DRIVEN THIN LIQUID FILMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING, SERIES B Volume 3, Number 3, Pages 85 96 c Institute for Scientific Computing and Information INVERSE TEMPERATURE RECONSTRUCTION IN THERMOCAPILLARY-DRIVEN
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More informationClick here to see an animation of the derivative
Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,
More informationRelativistic nuclear matter in generalized thermo-statistics
Relativistic nuclear matter in generalized termo-statistics K. Miyazaki E-mail: miyazakiro@rio.odn.ne.jp Abstract Te generalized Fermi-Dirac termo-statistics is developed for relativistic nuclear matter.
More informationAnalysis of Static and Dynamic Load on Hydrostatic Bearing with Variable Viscosity and Pressure
Indian Journal of Science and Tecnology Supplementary Article Analysis of Static and Dynamic Load on Hydrostatic Bearing wit Variable Viscosity and Pressure V. Srinivasan* Professor, Scool of Mecanical
More informationNotes on wavefunctions II: momentum wavefunctions
Notes on wavefunctions II: momentum wavefunctions and uncertainty Te state of a particle at any time is described by a wavefunction ψ(x). Tese wavefunction must cange wit time, since we know tat particles
More informationTravelling waves for a thin liquid film with surfactant on an inclined plane
IOP PUBLISHING Nonlinearity (009) 85 NONLINEARITY doi:0.088/095-775///006 Travelling waves for a tin liquid film wit surfactant on an inclined plane Vaagn Manukian and Stepen Scecter Matematics Department,
More informationHydraulic validation of the LHC cold mass heat exchanger tube.
Hydraulic validation o te LHC cold mass eat excanger tube. LHC Project Note 155 1998-07-22 (pilippe.provenaz@cern.c) Pilippe PROVENAZ / LHC-ACR Division Summary Te knowledge o te elium mass low vs. te
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationAN ANALYSIS OF AMPLITUDE AND PERIOD OF ALTERNATING ICE LOADS ON CONICAL STRUCTURES
Ice in te Environment: Proceedings of te 1t IAHR International Symposium on Ice Dunedin, New Zealand, nd t December International Association of Hydraulic Engineering and Researc AN ANALYSIS OF AMPLITUDE
More informationTheoretical Analysis of Flow Characteristics and Bearing Load for Mass-produced External Gear Pump
TECHNICAL PAPE Teoretical Analysis of Flow Caracteristics and Bearing Load for Mass-produced External Gear Pump N. YOSHIDA Tis paper presents teoretical equations for calculating pump flow rate and bearing
More information5.1 We will begin this section with the definition of a rational expression. We
Basic Properties and Reducing to Lowest Terms 5.1 We will begin tis section wit te definition of a rational epression. We will ten state te two basic properties associated wit rational epressions and go
More informationA Modified Distributed Lagrange Multiplier/Fictitious Domain Method for Particulate Flows with Collisions
A Modified Distributed Lagrange Multiplier/Fictitious Domain Metod for Particulate Flows wit Collisions P. Sing Department of Mecanical Engineering New Jersey Institute of Tecnology University Heigts Newark,
More informationComparison of Heat Transfer Conditions in. Tube Bundle Cross-Flow for Different Tube Shapes
Comparison of Heat Transfer Conditions in Tube Bundle Cross-Flow for Different Tube Sapes Dr. Andrej Horvat, researc associate* Jožef Stefan Institute, actor Engineering Division Jamova 39, SI 1001, Ljubljana,
More informationDevelopment of new and validation of existing convection correlations for rooms with displacement ventilation systems
Energy and Buildings 38 (2006) 163 173 www.elsevier.com/locate/enbuild Development of new and validation of existing convection correlations for rooms wit displacement ventilation systems Atila Novoselac
More informationEmpirical models for estimating liquefaction-induced lateral spread displacement
Empirical models for estimating liquefaction-induced lateral spread displacement J.J. Zang and J.X. Zao Institute of Geological & Nuclear Sciences Ltd, Lower Hutt, New Zealand. 2004 NZSEE Conference ABSTRACT:
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More information1. Consider the trigonometric function f(t) whose graph is shown below. Write down a possible formula for f(t).
. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd, periodic function tat as been sifted upwards, so we will use
More informationQuantum Mechanics Chapter 1.5: An illustration using measurements of particle spin.
I Introduction. Quantum Mecanics Capter.5: An illustration using measurements of particle spin. Quantum mecanics is a teory of pysics tat as been very successful in explaining and predicting many pysical
More informationApplication of numerical integration methods to continuously variable transmission dynamics models a
ttps://doi.org/10.1051/ssconf/2018440005 Application of numerical integration metods to continuously variable transmission dynamics models a Orlov Stepan 1 1 Peter te Great St. Petersburg Polytecnic University
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
More informationImplicit Multigrid Computations of Buoyant Drops Through Sinusoidal Constrictions
Metin Muradoglu 1 Department of Mecanical Engineering, Koc University, Rumelifeneri Yolu, Sariyer, Istanbul 34450, Turkey e-mail: mmuradoglu@ku.edu.tr Seckin Gokaltun 2 Computational Science and Engineering
More informationWork and Energy. Introduction. Work. PHY energy - J. Hedberg
Work and Energy PHY 207 - energy - J. Hedberg - 2017 1. Introduction 2. Work 3. Kinetic Energy 4. Potential Energy 5. Conservation of Mecanical Energy 6. Ex: Te Loop te Loop 7. Conservative and Non-conservative
More informationICES REPORT Isogeometric Phase-field Simulation of Boiling
ICES REPORT 15-16 June 015 Isogeometric Pase-field Simulation of Boiling by Ju Liu, Tomas J.R. Huges Te Institute for Computational Engineering and Sciences Te University of Texas at Austin Austin, Texas
More informationJournal of Applied Science and Agriculture. The Effects Of Corrugated Geometry On Flow And Heat Transfer In Corrugated Channel Using Nanofluid
Journal o Applied Science and Agriculture, 9() February 04, Pages: 408-47 AENSI Journals Journal o Applied Science and Agriculture ISSN 86-9 Journal ome page: www.aensiweb.com/jasa/index.tml Te Eects O
More information