Modelling of a free-surface ferrofluid flow
|
|
- Elinor Malone
- 6 years ago
- Views:
Transcription
1 Modelling of a free-surface ferrofluid flow M. Habera a,, J. Hron a a Mathematical Institute, Faculty of Mathematics and Physics, Charles University in Prague, Soolovsa 83, Prague Abstract The Cauchy s stress tensor of a ferrofluid exposed to an external magnetic field is subject to additional magnetic terms. For a linearly magnetizable medium, the terms results in interfacial magnetic force acting on the ferrofluid boundaries. This force changes the characteristics of many free-surface ferrofluid phenomena. The aim of this wor is to implement this force into the incompressible Navier-Stoes equations and propose a numerical method to solve them. The interface of ferrofluid is traced with the use of the characteristic level-set method and additional reinitialization step assures conservation of its volume. Incompressible Navier-Stoes equations are formulated for a divergence-free velocity fields while discrete interfacial forces are treated with continuous surface force model. Velocity-pressure coupling is implemented via the projection method. To predict the magnetic force effect quantitatively, Maxwell s equations for magnetostatics are solved in each time step. Finite element method is utilized for the spatial discretization. At the end of the wor, equilibrium droplet shape are compared to nown experimental results. 1. Introduction Free surface fluid flows and processes involved in a fluid behavior fascinated scientists since the very beginning of the scientific history. Problems as a breaup of a liquid jet, droplet formation and merging, rising bubbles etc. still lacs deeper understanding because of a complex and nonlinear equations governing such phenomena. In addition, they play a role in many industrial processes: fuel injection, fibre spinning, in-jet printing, etc. All these phenomena become even more attractive in terms of ferrohydrodynamics. Ferrofluid reacts to a magnetic field and changes its shape and rheology. This wor attempts to present a physical model and numerical method succesful in capturing multiphasefreesurface) ferrofluid flow with special focus on the interface development. 2. Theory and methods 2.1. Balance of momentum for ferrofluid Balance of momentum for continuous medium in Eulerian setting reads ρu) t + [u ] ρu) = T + f, 1) with ut, x) : I Ω R d being Eulerian velocity, a dimension d, T a Cauchy s stress tensor and ft, x) : I Ω R d volume density forces. Corresponding author addresses: habera@arlin.mff.cuni.cz M. Habera), hron@arlin.mff.cuni.cz J. Hron) The continuous medium Ω := Ω air t) Ω ferr t), t I = [0, T ], Ω C 0,1, in this wor is composed of two immiscible phases - surrounding fluid, e.g. air, Ω air t) and ferrofluid, Ω ferr t). So called constitutive relation must be specified. Here, we adopt several assumptionssimplifications) on the behaviour and properties of the medium, namely 1. the Ω fluid is incompressible, newtonian, the stress tensor without magnetic field T n = ˆT n D) = mi + 2ηD, with D := ) 1 2 u + u T the rate of deformation tensor, m the Lagrange multiplier pressure ) enforcing the incompressibility constraint, η dynamic viscosity, 2. the ferrofluid Ω ferr is linearly magnetizable, isotropic and homogenous, 3. there are no macroscopic electric currents and ferrofluid is non-conductive, we wor in field of magnetostatics, 4. the only effect of magnetic field on ferrofluid is additional magnetic stress tensor T m = µ 2 H 2 I + µh H. 2) With the assumptions above the total stress tensor T = ˆTD, H) = T n + T m. The surface tension effect must be included in multiphase flows. We chose so called continuous surface force, where surface forces are approximated with volume body forces concentrated on the interface. This allows us to write f = f s + f g, f s = σκδn, f g = ρge z, Preprint submitted to Elsevier July 14, 2016
2 where σ is surface tension coefficient, κ curvature, δ an infinitely smooth representative of the Dirac delta distribution concentrated on the interface, n outward normal to the interface, density ρ and gravitational acceleration g. Moreover, according to [8] there exists a tensor T s, such that 2.2. Magnetostatics f s = T s, T s = σδ I n n). As assumed above, the total stress tensor is function of magnetic field intensity H, it must be found a priori, before the momentum balance is solved. Magnetostatic equations and constitutive relation for linearly magnetizable medium are B = 0, H = 0, B = µh. 3) 2.3. Diffused interface, level-set function So called characteristic level-set function is constructed to trac the domains Ω air t), Ω ferr t) and their boundaries It is a regularized characteristic function of Ω ferr t), such that { C 0 x Ω air, distx, Ω air ) > /2, Ω), = 1 x Ω ferr, distx, Ω ferr ) > /2, and χ Ωferr L2 Ω) 0 as 0+. The epsilon parameter is the thicness of the diffused interface. With the aid of this function all discontinuous piecewise constant physical quantities are approximated, therefore ρ = ρ ferr ρ1 air, η = η ferr η1 air, µ = µ ferr µ1 air, 4) for density, dynamic viscosity and absolute permeability respectively. Outer normal to the implicitly represented interface is approximated as n Advection Since the evolution of boundaries Ω ferr t) and Ω air t) is the essence of this method, one must advect/transport the level-set function in the velocity field u. It is achieved solving t + u = 0. 5) This distorts the profile of level-set function and deteriorates it s property as regularized characteristic function of the ferrofluid phase. Level-set function is therefore reinitialized similarly to [5], [7], [6] τ + [ 1 )nt 0, x)] =, 6) note the change in the term on right-hand-side. tae full laplacian instead of normal derivative. We Numerical methodology Generally, evolution equations are here discretized in time with finite differences FDM). Finite element FEM) method is used to discretize in spatial variable. In the sense of FDM we denote a discrete solution at -th time level with upper index so u x) ut, x), etc. Because Cauchy s stress tensor is function of magnetic field, we must see H before the momentum balance is solved. For the Ω domain is simply connected the potential character of magnetic field intensity allows us to define H =: ξ. Wea formulation for the magnetostatics problem 3) is: find ξ 1 V := H0 Ω) such that for all s V µ 1 s ξ 1 n ds + µ 1 ξ 1, s ) = 0 7) Ω with u, v) := Ω uvdx being the standard L2 Ω) scalar product. The term ξ 1 n = H 1 n is substitued with boundary conditions. The balance of momentum 1) contains in each time step two unnown quantities, velocity u and pressure m enforcing the incompressibility constraint u = 0. To overcome this difficulty Chorin s projection method is used [9], [5], [7]. First find tentative u V := H 1 0 Ω; R d ), such that for all v V ρ u u 1), v ) ρ 1 u u, v ) m 1, v ) + η u + u Re ) T ), v ) + ) I 1 We 1 2 ), v ) 1 µ 1 Mg 2 H 2 ) 1, v + µ 1 H H, v ) = 0 8) Mg where Re, We are Reynolds and Weber numbers. The dimensionless number Mg := ρ ref u 2 ref µ ref is often called Magnetic Reynolds number. Href 2 Then, in so called pressure correction step find m V that for all q V u, q ) + m ref x ref u ref ρ ref m m 1 ) ρ 1 ), q = 0. Finally, correct the velocity to have zero divergence by finding u V that for all v V u u, v ) m ref m m 1 ) ) + x ref u ref ρ ref ρ 1, v = 0. 10) We use implicit Euler schema for the advection and reinitialization equation. Dimensionless wea formulation 9)
3 for the advection problem 5) is: find V := H 1 0 Ω) such that for all ϕ V 1 1, ϕ ) u ref x, u ϕ ) = 0. 11) ref Dimensionless wea formulation for the reinitialization 6): find n V such that for all ϕ V n τ n 1, ϕ ) n 1 n ) ), ϕ + n, ϕ) = 0. 12) In the Figure 1 we give a simple overview of the method used. Wea formulations and mesh generation are implemented into FEniCS [3]. Magnetostatics, 7) Balance of momentum, 8), 9), 10) 3. Results Advection, 11) Reinitialization, 12) next time step, 1 Figure 1: Schema of developed numerical method. p m h Figure 2: Convergence of pressure inside a droplet in equilibrium. y/cm x/cm Simple model validations Here we give a simple validation of the developed model and methodology in trivial static cases. In all examples the dimension d = 2. First, consider an ideally circular ferrofluid droplet with radius r = 2 mm placed in the center of square computational domain Ω. From well nown Young-Laplace equation there is a jump in pressures, p = σ R, for droplet with radius R and surface tension σ. Magnetic field is not assumed. We depicted in Figure 2 a ratio of the theoretical pressure difference and the incompressibility constraint p jump, m, computed from our model. This ratio is evaluated for several values of mesh cell diameter, h. As a mesh is more refined and h 0 + we have p m 1. Another simple example is a magnetic field intensity H in a domain with ferrofluid droplet. Let have a droplet with radius r = 1 mm, permeability of the surrounding phase be the permeability of free space, µ air = 4π 10 7 N A 2 and susceptibility of the droplet be χ = 1. Boundary conditions are set in 7) so that H n = 1 A m 1 at top and bottom boundaries. Note, that due to applied magnetic field, the magnetic intensity increases at the top and bottom tips of the droplet. The magnetic stress tensor 2) yields to the force which elongate the droplet. 3 Figure 3: Magnetic field intensity H in A m 1 with ferrofluid droplet present Equillibrium ferrofluid droplet shape Let have ferrofluid droplet of initial radius r = 1 mm and set the boundary conditions H n = 12 A m 1. The droplet elongates [10], [13] due to the effect of magnetic field. Developed model could capture this behaviour. In equillibrium state the level-set function of ferrofluid phase is depicted in Figure 4 and magnetic field intensity in Figure 5. Parameters are inspired by [10] where experimental results are presented. Their experiment is compared to our model in Figure Ferrofluid dripping Ferrofluid dripping is experimentally studied in [1] and [2]. It is noted, that magnetic field substantially changes equilibrium of droplet formation and formed droplet and evolution dynamics are affected by magnetic field. Here we give a simple numerical example of such effect. Pendant droplet profile is depicted in Figure 6.
4 Figure 6: Pendant droplet profile. Magnetic field H n = 4 A m 1 and χ = 1 in the left picture, while no magnetic field is present in the right picture. Profiles are taen at time 33 ms after simulation started. Figure 4: Equillibrium droplet shape. Boundary conditions H n = 12 A m 1 and χ = Surface tension is set to σ = 13.5 mn m 1. Viscosity and density do not play role in equillibrium state, since u 0. Experimental result from [10] is included in actual dimensions in dotted rectangle. y/cm x/cm time evolution is done solving advection equation. To preserve volume of ferrofluid level-set function is reinitialized. In the last part we give wea formulations of derived equations and show some simple numerical experiments. The possible extensions of this wor could be the implementation of non-newtonian constitutive laws, with viscosity depending on magnetic field, etc. [4] Another numerical experiments are needed to verify the validity of the model. Dynamics of droplet pinch-off process and influence of magnetic field on it are now experimentally studied [1], [2], so simulations of this phenomena are be welcomed. Forces due to surface tension and magnetic field are written in terms of tensors, wea formulation maes it possible to per-partes derivatives and less regularity is needed for physical quantities. Several wors with finite-volume approach were validated in recent years [10], [11], [12], [13] on the other hand, this is the first time finite-elements are used, according to authors best nowledge. Figure 5: Magnetic field for equillibrium droplet shape. Boundary conditions H n = 12 A m 1 and χ = Surface tension is set to σ = 13.5 mn m Conclusions The main goal of this wor was to formulate the equations of ferrohydrodynamics and to put together advanced mathematical techniques, so simple free-surface ferrofluid problems could be attaced. In the beginning a ferrofluid is treated in terms of continuum mechanics and several assumptions lead us to the balance law form. Continuous surface force method is utilized for surface tension force. The implicit description of the interface is adopted and 4
5 References [1] HABERA, Michal, Michal FABIAN, Maria SVIKOVA and Milan TIMKO. The influence of magnetic field on free surface ferrofluid flow. Magnetohydrodynamics. 2013, 3-4). ISSN X. [2] FABIAN, Michal, Peter BURDA, Maria SVIKOVA and Robert HUNADY. The Influence of Magnetic Field on the Separation of Droplets from Ferrofluid Jet. Journal of Magnetism and Magnetic materials - special ICMF14 issuesubmitted), [3] ALNAES, M.S., A. LOGG, K A. MARDAL, O. SKAVHAUG and H.P. LANGTANGEN. Unified framewor for finite element assembly. International Journal of Computational Science and Engineering [online]. 2009, 44). DOI: /ijcse [4] ODENBACH, S., L. M. POP a A. Yu. ZUBAREV. Rheological properties of magnetic fluids and their microstructural bacground. GAMM-Mitteilungen [online]. 2007, 301), DOI: /gamm ISSN [5] OLSSON, Elin and Gunilla KREISS. A conservative level set method for two phase flow. Journal of Computational Physics [online]. 2005, 2101): DOI: /j.jcp [6] SUSSMAN, Mar, Emad FATEMI, Peter SMEREKA and Stanley OSHER. An improved level set method for incompressible two-phase flows. Computers & Fluids [online]. 1998, 275-6): DOI: /s ) [7] OLSSON, Elin, Gunilla KREISS and Sara ZAHEDI. A conservative level set method for two phase flow II. Journal of Computational Physics [online]. 2007, 2251): DOI: /j.jcp [8] LAFAURIE, Bruno, Carlo NARDONE, Ruben SCAR- DOVELLI, Stphane ZALESKI and Gianluigi ZANETTI. Modelling Merging and Fragmentation in Multiphase Flows with SURFER. Journal of Computational Physics [online]. 1994, 1131): DOI: /jcph [9] CHORIN, Alexandre Joel. Numerical solution of the Navier- Stoes equations. Mathematics of Computation [online]. 1968, 22104): DOI: /s [10] AFKHAMI, S., A. J. TYLER, Y. RENARDY, M. RENARDY, T. G. St. PIERRE, R. C. WOODWARD and J. S. RIF- FLE. Deformation of a hydrophobic ferrofluid droplet suspended in a viscous medium under uniform magnetic fields. Journal of Fluid Mechanics [online]. 2010, 663: DOI: /s [11] KI, Hyungson. Level set method for two-phase incompressible flows under magnetic fields. Computer Physics Communications [online]. 2010, 1816): DOI: /j.cpc [12] SHI, Dongxiao, Qincheng BI and Rongqi ZHOU. Numerical Simulation of a Falling Ferrofluid Droplet in a Uniform Magnetic Field by the VOSET Method. Numerical Heat Transfer, Part A: Applications [online]. 2014, 662): DOI: / [13] LAVROVA, O., V. POLEVIKOV and L. TOBISKA. Equilibrium shapes of a ferrofluid drop. PAMM [online]. 2005, 51): DOI: /pamm
Deformation of a hydrophobic ferrofluid droplet suspended in a viscous medium under uniform magnetic fields
Under consideration for publication in J. Fluid Mech. 1 Deformation of a hydrophobic ferrofluid droplet suspended in a viscous medium under uniform magnetic fields S. AFKHAMI 1,, A. J. TYLER 3, Y. RENARDY,
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 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 informationDeformation of a hydrophobic ferrofluid droplet suspended in a viscous medium under uniform magnetic fields
J. Fluid Mech. (010), vol. 663, pp. 358 384. c Cambridge University Press 010 doi:10.1017/s0011010003551 Deformation of a hydrophobic ferrofluid droplet suspended in a viscous medium under uniform magnetic
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 informationINDEX 363. Cartesian coordinates 19,20,42, 67, 83 Cartesian tensors 84, 87, 226
INDEX 363 A Absolute differentiation 120 Absolute scalar field 43 Absolute tensor 45,46,47,48 Acceleration 121, 190, 192 Action integral 198 Addition of systems 6, 51 Addition of tensors 6, 51 Adherence
More informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
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 informationFinite Element Techniques for the Numerical Simulation of Two-Phase Flows with Mass Transport
Finite Element Techniques for the Numerical Simulation of Two-Phase Flows with Mass Transport Christoph Lehrenfeld and Arnold Reusken Preprint No. 413 December 2014 Key words: Two-phase flow, mass transport,
More informationFEM-Level Set Techniques for Multiphase Flow --- Some recent results
FEM-Level Set Techniques for Multiphase Flow --- Some recent results ENUMATH09, Uppsala Stefan Turek, Otto Mierka, Dmitri Kuzmin, Shuren Hysing Institut für Angewandte Mathematik, TU Dortmund http://www.mathematik.tu-dortmund.de/ls3
More informationDevelopment and analysis of a Lagrange-Remap sharp interface solver for stable and accurate atomization computations
ICLASS 2012, 12 th Triennial International Conference on Liquid Atomization and Spray Systems, Heidelberg, Germany, September 2-6, 2012 Development and analysis of a Lagrange-Remap sharp interface solver
More informationNonlinear oscillations and rotations of a liquid droplet
Nonlinear oscillations and rotations of a liquid droplet Tadashi Watanabe watanabe.tadashi66@jaea.go.jp Abstract Nonlinear oscillations and rotations of a liquid droplet are simulated numerically by solving
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 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 informationFERROHYDRODYNAMICS R. E. ROSENSWEIG. DOVER PUBLICATIONS, INC. Mineola, New York
FERROHYDRODYNAMICS R. E. ROSENSWEIG DOVER PUBLICATIONS, INC. Mineola, New York CONTENTS Preface 1 Introduction 1.1 Scope of ferrohydrodynamics 1.2 Ferromagnetic solids 1.3 Magnetic fluids 1.4 Ferromagnetic
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 informationA comparison of viscoelastic stress wakes for 2D and 3D Newtonian drop deformation in a viscoelastic matrix under shear
A comparison of viscoelastic stress wakes for 2D and 3D Newtonian drop deformation in a viscoelastic matrix under shear S. Afkhami, P. Yue, and Y. Renardy Department of Mathematics, 460 McBryde Hall, Virginia
More informationSimulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework
Simulating Interfacial Tension of a Falling Drop in a Moving Mesh Framework Anja R. Paschedag a,, Blair Perot b a TU Berlin, Institute of Chemical Engineering, 10623 Berlin, Germany b University of Massachusetts,
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 informationShear-induced rupturing of a viscous drop in a Bingham liquid
J. Non-Newtonian Fluid Mech. 95 (2000) 235 251 Shear-induced rupturing of a viscous drop in a Bingham liquid Jie Li, Yuriko Y. Renardy Department of Mathematics and ICAM, Virginia Polytechnic Institute
More informationLes Houches School of Foam: Rheology of Complex Fluids
Les Houches School of Foam: Rheology of Complex Fluids Andrew Belmonte The W. G. Pritchard Laboratories Department of Mathematics, Penn State University 1 Fluid Dynamics (tossing a coin) Les Houches Winter
More information- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum)
2.20 - Marine Hydrodynamics, Spring 2005 Lecture 4 2.20 - Marine Hydrodynamics Lecture 4 Introduction Governing Equations so far: Knowns Equations # Unknowns # density ρ( x, t) Continuity 1 velocities
More informationAn immersed boundary method for interfacial flows with insoluble surfactant
An immersed boundary method for interfacial flows with insoluble surfactant Ming-Chih Lai Yu-Hau Tseng Huaxiong Huang Abstract In this paper, an immersed boundary method is proposed for the simulation
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 informationA MULTIPHASE FLUID-SOLID MODEL BASED ON THE LEVEL SET METHOD
Ninth International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 10-12 December 2012 A MULTIPHASE FLUID-SOLID MODEL BASED ON THE LEVEL SET METHOD James E. HILTON
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 informationNumerical Modeling of Fluid Interface Phenomena
Numerical Modeling of Fluid Interface Phenomena SARA ZAHEDI Avhandling som med tillstånd av Kungliga Tekniska högskolan framlägges till offentlig granskning för avläggande av teknologie licentiatexamen
More informationNavier-Stokes Equation: Principle of Conservation of Momentum
Navier-tokes Equation: Principle of Conservation of Momentum R. hankar ubramanian Department of Chemical and Biomolecular Engineering Clarkson University Newton formulated the principle of conservation
More informationIntroduction to Fluid Mechanics
Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the
More informationTurbulence Models for Flows with Free Surfaces and Interfaces
AIAA JOURNAL Vol. 44, No. 7, July 2006 Turbulence Models for Flows with Free Surfaces and Interfaces Ebrahim Shirani, Ali Jafari, and Nasser Ashgriz University of Toronto, Toronto, Ontario M5S 3G8, Canada
More informationSimulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang,
Simulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang, yfwang09@stanford.edu 1. Problem setting In this project, we present a benchmark simulation for segmented flows, which contain
More informationx j r i V i,j+1/2 r Ci,j Ui+1/2,j U i-1/2,j Vi,j-1/2
Merging of drops to form bamboo waves Yuriko Y. Renardy and Jie Li Department of Mathematics and ICAM Virginia Polytechnic Institute and State University Blacksburg, VA -, U.S.A. May, Abstract Topological
More informationLevel Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows
Level Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows Abstract Maged Ismail Claremont Graduate University Level Set and Phase Field methods are well-known interface-capturing
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 informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationDetailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and Boundary Condition
ILASS-Americas 25th Annual Conference on Liquid Atomization and Spray Systems, Pittsburgh, PA, May 23 Detailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and
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 informationHeat Transfer Benchmark Problems Verification of Finite Volume Particle (FVP) Method-based Code
PROCEEDING OF 3 RD INTERNATIONAL CONFERENCE ON RESEARCH, IMPLEMENTATION AND EDUCATION OF MATHEMATICS AND SCIENCE YOGYAKARTA, 16 17 MAY 2016 Heat Transfer Benchmark Problems Verification of Finite Volume
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More informationDirect Numerical Simulation of Single Bubble Rising in Viscous Stagnant Liquid
Direct Numerical Simulation of Single Bubble Rising in Viscous Stagnant Liquid Nima. Samkhaniani, Azar. Ajami, Mohammad Hassan. Kayhani, Ali. Sarreshteh Dari Abstract In this paper, direct numerical simulation
More informationMagnetic Drug Targeting in Cancer Therapy
Magnetic Drug Targeting in Cancer Therapy SOLVED WITH COMSOL MULTIPHYSICS 3.5a COPYRIGHT 2008. All right reserved. No part of this documentation may be photocopied or reproduced in any form without prior
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 informationA COMPUTATIONAL STUDY OF THE COALESCENCE PROCESS BETWEEN A DROP AND AN INTERFACE IN AN ELECTRIC FIELD
6th International Conference on CFD in the Oil & Gas, Metallurgical and Process Industries SINTEF / NTNU, Trondheim, Norway 1 12 June 28 CFD28-78 A COMPUTATIONAL STUDY OF THE COALESCENCE PROCESS BETWEEN
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationA PARALLEL LEVEL-SET APPROACH FOR TWO-PHASE FLOW PROBLEMS WITH SURFACE TENSION IN THREE SPACE DIMENSIONS
A PARALLEL LEVEL-SET APPROACH FOR TWO-PHASE FLOW PROBLEMS WITH SURFACE TENSION IN THREE SPACE DIMENSIONS ROBERTO CROCE, MICHAEL GRIEBEL, AND MARC ALEXANDER SCHWEITZER Abstract. In this paper we present
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationFinite volume method for two-phase flows using level set formulation
Finite volume method for two-phase flows using level set formulation Peter Frolkovič,a,1, Dmitry Logashenko b, Christian Wehner c, Gabriel Wittum c a Department of Mathematics and Descriptive Geometry,
More informationFlow Field and Oscillation Frequency of a Rotating Liquid Droplet
Flow Field and Oscillation Frequency of a Rotating Liquid Droplet TADASHI WATANABE Center for Computational Science and e-systems Japan Atomic Energy Agency (JAEA) Tokai-mura, Naka-gun, Ibaraki-ken, 319-1195
More informationAn Extended Finite Element Method for a Two-Phase Stokes problem
XFEM project An Extended Finite Element Method for a Two-Phase Stokes problem P. Lederer, C. Pfeiler, C. Wintersteiger Advisor: Dr. C. Lehrenfeld August 5, 2015 Contents 1 Problem description 2 1.1 Physics.........................................
More information( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #
Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue
More informationOblique Drop Impact on Deep and Shallow Liquid
1 2 3 4 5 6 7 8 9 11 10 12 13 14 15 16 17 18 19 20 21 22 23 24 25 Commun. Comput. Phys. doi: 10.4208/cicp.XXX.XXX Oblique Drop Impact on Deep and Shallow Liquid B. Ray 1, G. Biswas 1,2, and A. Sharma 3
More informationPart IV: Numerical schemes for the phase-filed model
Part IV: Numerical schemes for the phase-filed model Jie Shen Department of Mathematics Purdue University IMS, Singapore July 29-3, 29 The complete set of governing equations Find u, p, (φ, ξ) such that
More informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
More information3. FORMS OF GOVERNING EQUATIONS IN CFD
3. FORMS OF GOVERNING EQUATIONS IN CFD 3.1. Governing and model equations in CFD Fluid flows are governed by the Navier-Stokes equations (N-S), which simpler, inviscid, form is the Euler equations. For
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationAA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow
AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit
More informationFluid flows through unsaturated porous media: An alternative simulation procedure
Engineering Conferences International ECI Digital Archives 5th International Conference on Porous Media and Their Applications in Science, Engineering and Industry Refereed Proceedings Summer 6-24-2014
More informationABSTRACT NUMERICAL SIMULATION OF LIQUID - VAPOR INTERFACES IN THE SHARP INTERFACE LIMIT
ABSTRACT Title of Thesis: NUMERICAL SIMULATION OF LIQUID - VAPOR INTERFACES IN THE SHARP INTERFACE LIMIT Vamsee K. Yerramilli, Master of Science, 2005 Thesis directed by: Dr. Elias Balaras Department of
More informationEfficient simulation techniques for incompressible two-phase flow
3D-Surface Engineering für Werkzeugsysteme der Blechformteilefertigung - Erzeugung, Modellierung, Bearbeitung - Efficient simulation techniques for incompressible two-phase flow O. Mierka, O. Ouazzi, T.
More informationMath 575-Lecture Viscous Newtonian fluid and the Navier-Stokes equations
Math 575-Lecture 13 In 1845, tokes extended Newton s original idea to find a constitutive law which relates the Cauchy stress tensor to the velocity gradient, and then derived a system of equations. The
More informationSimulations with Dispersive Mixtures Model
Simulations with Dispersive Mixtures Model P.J.A. Janssen, P.D. Anderson & G.W.M. Peters MT03.05 May 12, 2003 1 Introduction The combination of two different materials in a blend can yield interesting
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 informationComputational Fluid Dynamics 2
Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationBoundary Conditions in Fluid Mechanics
Boundary Conditions in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University The governing equations for the velocity and pressure fields are partial
More informationA Permeability Function for Brinkman s Equation
Proceedings of the th WSEAS Int. Conf. on MATHEMATICAL METHOS, COMPUTATIONAL TECHNIQUES AN INTELLIGENT SYSTEMS A Permeability Function for Brinman s Equation M.H. HAMAN, M.T. KAMEL, and H.I. SIYYAM 3,
More informationHeight functions for applying contact angles to 2D VOF simulations
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids (2007) Published online in Wiley InterScience (www.interscience.wiley.com). DOI: 0.002/fld.65 Height functions for applying
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More informationDynamics of magnetic chains in a shear flow under the influence of a uniform magnetic field
Dynamics of magnetic chains in a shear flow under the influence of a uniform magnetic field Tae Gon Kang, 1, a) Martien A. Hulsen, 2 and Jaap M.J. den Toonder 2 1) School of Aerospace and Mechanical Engineering,
More informationInternational Journal of Multiphase Flow
International Journal of Multiphase Flow 52 (2013) 22 34 Contents lists available at SciVerse ScienceDirect International Journal of Multiphase Flow journal homepage: www.elsevier.com/locate/ijmulflow
More informationFVM for Fluid-Structure Interaction with Large Structural Displacements
FVM for Fluid-Structure Interaction with Large Structural Displacements Željko Tuković and Hrvoje Jasak Zeljko.Tukovic@fsb.hr, h.jasak@wikki.co.uk Faculty of Mechanical Engineering and Naval Architecture
More informationNumerical simulation of breakup of a viscous drop in simple shear flow through a volume-of-fluid method
Numerical simulation of breakup of a viscous drop in simple shear flow through a volume-of-fluid method Jie Li, Yuriko Y. Renardy, and Michael Renardy Citation: Physics of Fluids (1994-present) 12, 269
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 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 informationNUMERICAL INVESTIGATION OF THERMOCAPILLARY INDUCED MOTION OF A LIQUID SLUG IN A CAPILLARY TUBE
Proceedings of the Asian Conference on Thermal Sciences 2017, 1st ACTS March 26-30, 2017, Jeju Island, Korea ACTS-P00786 NUMERICAL INVESTIGATION OF THERMOCAPILLARY INDUCED MOTION OF A LIQUID SLUG IN A
More informationContinuum Mechanics Fundamentals
Continuum Mechanics Fundamentals James R. Rice, notes for ES 220, 12 November 2009; corrections 9 December 2009 Conserved Quantities Let a conseved quantity have amount F per unit volume. Examples are
More informationarxiv:chao-dyn/ v2 24 May 1996
Two-dimensional Navier Stokes simulation of deformation and break up of liquid patches Stéphane Zaleski arxiv:chao-dyn/9502016v2 24 May 1996 Jie Li Sauro Succi Laboratoire de Modélisation en Mécanique,
More informationSimulation of a Pressure Driven Droplet Generator
Simulation of a Pressure Driven Droplet Generator V. Mamet* 1, P. Namy 2, N. Berri 1, L. Tatoulian 1, P. Ehouarn 1, V. Briday 1, P. Clémenceau 1 and B. Dupont 1 1 DBV Technologies, 2 SIMTEC *84 rue des
More informationInternational Journal of Multiphase Flow
International Journal of Multiphase Flow 36 (010) 333 34 Contents lists available at ScienceDirect International Journal of Multiphase Flow journal homepage: www.elsevier.com/locate/ijmulflow A coupled
More informationChapter 1 Fluid Characteristics
Chapter 1 Fluid Characteristics 1.1 Introduction 1.1.1 Phases Solid increasing increasing spacing and intermolecular liquid latitude of cohesive Fluid gas (vapor) molecular force plasma motion 1.1.2 Fluidity
More informationA SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations
A SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations Hervé Lemonnier DM2S/STMF/LIEFT, CEA/Grenoble, 38054 Grenoble Cedex 9 Ph. +33(0)4 38 78 45 40 herve.lemonnier@cea.fr, herve.lemonnier.sci.free.fr/tpf/tpf.htm
More informationdynamics of f luids in porous media
dynamics of f luids in porous media Jacob Bear Department of Civil Engineering Technion Israel Institute of Technology, Haifa DOVER PUBLICATIONS, INC. New York Contents Preface xvii CHAPTER 1 Introduction
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 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 informationThe effect of insoluble surfactant at dilute concentration on drop breakup under shear with inertia
PHYSICS OF FLUIDS VOLUME 16, NUMBER 1 JANUARY 2004 The effect of insoluble surfactant at dilute concentration on drop breakup under shear with inertia M. A. Drumright-Clarke Southwest Research Institute,
More informationNonlinear Deformation of a Ferrofluid Droplet in a Uniform Magnetic Field
pubs.acs.org/ Nonlinear Deformation of a Ferrofluid Droplet in a Uniform Magnetic Field Gui-Ping Zhu, Nam-Trung Nguyen,*, R. V. Ramanujan, and Xiao-Yang Huang School of Mechanical and Aerospace Engineering,
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 informationA Numerical Study of an Injection-Compression Molding Process by using a Moving Grid Bambang Arip Dwiyantoro
Applied Mechanics and Materials Online: 2012-12-13 IN: 1662-7482, ols. 249-250, pp 472-476 doi:10.4028/www.scientific.net/amm.249-250.472 2013 Trans Tech Publications, witzerland A Numerical tudy of an
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 1 June 2006 1.30 to 4.30 PAPER 76 NONLINEAR CONTINUUM MECHANICS Attempt FOUR questions. There are SIX questions in total. The questions carry equal weight. STATIONERY
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 informationCFD in COMSOL Multiphysics
CFD in COMSOL Multiphysics Mats Nigam Copyright 2016 COMSOL. Any of the images, text, and equations here may be copied and modified for your own internal use. All trademarks are the property of their respective
More information1 Modeling Immiscible Fluid Flow in Porous Media
Excerpts from the Habilitation Thesis of Peter Bastian. For references as well as the full text, see http://cox.iwr.uni-heidelberg.de/people/peter/pdf/bastian_habilitationthesis.pdf. Used with permission.
More informationA level set projection model of lipid vesicles in general flows
, A level set projection model of lipid vesicles in general flows D. Salac a, M. Miksis a a Northwestern University, Engineering Sciences and Applied Mathematics, Evanston, IL, 60208 Abstract A new numerical
More informationCourse Syllabus: Continuum Mechanics - ME 212A
Course Syllabus: Continuum Mechanics - ME 212A Division Course Number Course Title Academic Semester Physical Science and Engineering Division ME 212A Continuum Mechanics Fall Academic Year 2017/2018 Semester
More informationA numerical analysis of drop impact on liquid film by using a level set method
Journal of Mechanical Science and Technology 25 (10) (2011) 2567~2572 wwwspringerlinkcom/content/1738-494x DOI 101007/s12206-011-0613-7 A numerical analysis of drop impact on liquid film by using a level
More informationReview of fluid dynamics
Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence
More informationLevel Set Tumor Growth Model
Level Set Tumor Growth Model Andrew Nordquist and Rakesh Ranjan, PhD University of Texas, San Antonio July 29, 2013 Andrew Nordquist and Rakesh Ranjan, PhD (University Level Set of Texas, TumorSan Growth
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 informationSurface phase separation and flow in a simple model of drops and vesicles
Surface phase separation and flow in a simple model of drops and vesicles Tutorial Lecture 4 John Lowengrub Department of Mathematics University of California at Irvine Joint with J.-J. Xu (UCI), S. Li
More informationDEVELOPING A MICRO-SCALE MODEL OF SOIL FREEZING
Proceedings of ALGORITMY 2016 pp. 234 243 DEVELOPING A MICRO-SCALE MODEL OF SOIL FREEZING A. ŽÁK, M. BENEŠ, AND T.H. ILLANGASEKARE Abstract. In this contribution, we analyze thermal and mechanical effects
More information