Mapping Molecular Dynamics to Mesoscopic Models: Challenges at Interfaces
|
|
- Judith Lilian Sutton
- 5 years ago
- Views:
Transcription
1 Mapping Molecular Dynamics to Mesoscopic Models: Challenges at Interfaces Mark Robbins, Johns Hopkins University Shiyi Chen, Shengfeng Cheng, Colin Denniston, Weinan E, Jin Liu, Xiaobo Nie, Peter Thompson Development and Analysis of Multiscale Methods University of Minnesota, Novermber 6, 008 Funding: National Science Foundation Nanomotor Challenge When composition, strain rate, stress or other property changes rapidly compared to correlation length or mean free path, hard to describe with coarse-grained equations. Show examples for fluid/solid and fluid/fluid interfaces, and flows where continuum equations predict singularities Funding: National Science Foundation Nanomotor 1
2 Atomic Interactions Fix Continuum Behavior Continuum models need: Boundary conditions (BC) + Constitutive Relations velocity or stress stress vs. strain (rate) slip, friction, adhesion viscous, elastic, plastic Traditionally assume simple form, fit to experiment: Assume no-slip for fluids, Simple analytic functions simple friction and Linear response gives adhesion laws for solids. viscosity or elastic modulii Molecular simulations determine BC s and constitutive relations. Can answer: Do usual assumptions for BC and constitutive laws work? Down to what scale do continuum equations apply? Is there new mesoscopic behavior between atomic and bulk? Issue: When is it better to find relations on the fly vs. fitting simple relations? Simulations: Generic behavior simple potential Spherical molecules, Lennard-Jones (LJ) interactions: V LJ =4ε ij [(a ij /r) 1 -(a ij /r) 6 ] for r<r c (usually r c = 1/6 a) Add unbreakable FENE bonds to make n-mers Solid walls Usually flat surface of fcc crystal (111) Usually rigid or atoms held to sites by springs Characteristic LJ energy and length depend on species w solid wall, 1 and fluid atoms Vary ε 1, a 1, make fluids miscible or immiscible Vary ε Iw, a Iw, change wetting, flow BC Units: average a (0.3nm), ε, m, and τ=a(m/ε) 1/ (3ps)
3 Continuum Eqs: Newtonian, steady, low Re, Ma 0 Conservation of each species i: ( ρ ) i uiα ρ i -mass density, u iα velocity Conservation of momentum: α σ βα = 0 σ αβ =stress tensor, σ αα /3=-p, p-pressure Constitutive relation for stress: σ x z =ηb η b - viscosity, u=mass-averaged velocity α = Need BC s: u iz =0, σ zz =const u x (0) = slip velocity or length L s? z u ix (0)-u x (0) = surf. diffusion? σ (0)-σ (wall) = surf. stress? L s wall at z=0 Want to apply to nearly molecular scale flows Have fit both sharp interface and mesoscopic BC s u x z u x x Single Fluid Structure Near walls layering & in-plane order Average over layers to generate mesoscopic model and compare with continuum Choice of wall position z=0 is not unique at molecular level. Choose center of first layer of solid atoms Midpoint between solid and fluid layer gives same physics, different numbers - local ρ - layer ave. Wall σ z 3
4 Single Fluid Slip Bulk flow extrapolates (lines) to u(0) u wall (shaded bars) u s u(0)-u wall Define local η from local v η z v x (z)=σ (or finite difference) Then find Navier slip condition w u s = dz z ( u v ) = L sσ / η b b b L s dz ( η b / η 1) w L s > 0 slip, L s < 0 stick L s /a typically ~- 0a can be much bigger for polymers (Thompson & Robbins PRA41, 6830 (1990)) u s u s monomers 16-mers 16-mers exponential decay Wall Induces layering Density vs. height z between walls Deviations from no-slip BC Viscosity bulk within ~ σ Thompson & Robbins, Phys. Rev. A41, 6830 (1990) 4
5 In-plane Order Controls L S ρ(r) Layer Layer 1 S(q) Fluid atoms cluster between wall atoms Atoms in next fluid layer between atoms of first layer Results for all interaction parameters, densities, collapse when plotted against S max =S(G) Quantify with S(q)= ρ(q) Thompson & Robbins PRA41, 6830 (1990)) Boundary Conditions For Mixtures? What average velocity should we apply the boundary condition (& const. relation) to? Mass average u= (ρ 1 u 1 + ρ u )/ ρ Number average? Volume average? Koplik and Banavar, PRL 80, 515 ( 98). v 1 =v at wall, both vanish. Ganesan & Brenner, PRL 8, 1333 ( 99); PRE 61, 6879 ( 00). If true, not a useful B.C. for macroscopic equations: D=0. Denniston and Robbins, PRL 87, ( 01), JCP 15, 1410 ( 06) Diffusion leads to v 1 v, wetting drives flows, no unique choice of average velocity for BC s Unexpected contributions to slip velocity 5
6 Diffusion in Concentration Gradient Maxwell Demon s Wall z y x Lennard-Jones interactions Diffusion bulk on same length as viscosity bulk Slip D higher at wall slip J A =ρ A (u A -u)=-ρd c A stick Simulation times ~ μs. 1 rich Wetting + Conc. Gradient Couette or Poisseuille-like flow between static walls rich z x u s u s u s ~ m/s z/a z/a 6
7 Flux of Individual Species Large diffusive flux J i =ρ i (v i -v) v 1 Const in center since x ρ i, D const. Drops near walls since x ρ i, D drop asymmetric wetting symmetric wetting v v 1 v Generalized Slip Boundary Condition Sharp interface model for slip can only depend on continuum values at wall: σb, x γ and x c 1 γ=surface stress, c 1 =concentration of species 1 L b L f1 f s m us = σ + xγ + D xc 1 ηb ηb c1 c L L s m dz( η / η 1) b dz( η / η)( σ b σ (0)) / γ Usual Navier condition is limiting case x 0 f i depend on relative momentum transfer of two components last term vanishes for some average u BUT which depends on atomic interactions, mass, For any other average can have net flow with no drag x 7
8 Marangoni Surface Stress Surface stress x γ near wall since σ xx =σ yy σ zz =-p b γ dz( σ xx + p) w Force balance gives shear stress BC b wall Δ σ σ ( w) σ = xγ -p Δσ a 3 /ε asymmetric wetting Stress BC obeyed in simulations Gradient in σ xx in fluid drives interfacial flow x γ a 3 /ε L M dz( η / η)( σ σ (0)) / γ b x Test of Slip Boundary Condition Measured slip agrees with calculated u L = L + f γ + s b m 1 s σ x D xc 1 ηb ηb c1 c f Diffusion term =0 for identical atoms Constant L s ~ -0.6a, stick Constant L m ~ 0.7a measured slip calculated L s term x γ term Denniston & Robbins, JCP 15, 1410 ( 06) 8
9 u Flow without drag neutral wetting L L f s b m 1 s = σ + xγ + D xc 1 ηb ηb c1 c P>0 =0 P=0 f 0 if m 1 > m, couple to wall differently, etc. P<0 Mass flux due to diffusion with no pressure gradient. Adding pressure gradient Extra Poiseuille flow with constant L s Denniston & Robbins, PRL 87, ( 01); JCP 15, 1410 ( 06) Optimize Parameters for Nanomotor Efficiency of interfacial stress transfer b dz( η / η 1)( σ σ (0)) / γ b x v x a/τ Velocities ~m/s Δσ a 3 /ε Surface stress ~MPa z/a x γ a 3 /ε 9
10 Summary for Solid-Fluid Interfaces Atomistic flows fit by generalized sharp interface boundary conditions down to scales ~ a few nm Stress balance b wall Δ σ σ ( w) σ = xγ Slip boundary condition L b L f1 f s m us = σ + xγ + D c x 1 ηb ηb c1 c Navier Surface stress Diffusive Can choose average velocity to make last term vanish, but not always mass, number or volume average Considered design and optimization of nanomotors and pumps: u s ~ m/s, stresses ~ MPa Difficult to find on fly because diffusive flux small What about Fluid Interfaces? Mesoscopic Model for Mixtures Conservation of each species i: ( ) α ρ i uiα = 0 ρ i -mass density, u iα velocity Conservation of momentum: Navier-Stokes Equations Constitutive relation for stress: 1 σ αβ = p0δ αβ + η ( ) αuβ + βuα +??? η - viscosity, u=mass-averaged velocity??? stress from variations in free energy F with mass-averaged density ρ and conc. φ= ρ A - ρ B Diffusive flux from gradient in chemical potential μ tφ + α ( φuα ) = Γ μ Lattice-Boltzmann method efficient dynamics that guarantees conservation, but need F and μ 10
11 Phase Field Model for Mixtures If correlation length ξ» σ coarse-grained F as functional of ρ, φ Usually write simplest form with proper symmetry 1 { ψ ( ρ, φ, T ) + ( ) + ( ) + ( } 1 1 φ K ρ ρ Kφ φ Kρφ ρ F = ) Typically ignore K ρφ (odd in φ), assume K ρρ >0 Often fix ρ, expand ψ as quartic polynomial in φ Goal: Obtain ψ and K s from MD None of usual assumptions is good Denniston & Robbins, Physical Review E69, (004) Model Lennard-Jones truncated to give pure repulsion Temperature mostly irrelevant Interaction energy between unlike species greater by factor ε* drives segregation coexistence region 11
12 Finding K From Linear Response Apply perturbation: δμ ρ = μ ρθ sinqx, δμ φ =-μ φθ sinqx Measure resulting ρ q, φ q Lρρ Lρφ ρq μρq = Lρφ Lφφ φq uφq L ρρ = ψ/ ρ + K ρρ q, Find L s linear in q to π/q σ, but K ρρ,k ρφ <0 Usual to add q 4 terms if K ρρ <0 but expect K ρρ 0 as q since cost of forcing atoms onto arbitrarily fine lattice 0 Fit F using only: Fitting F 1) Phase boundaries of coexistence region ) Linear response near coexistence 1
13 Test of Model No information about the interface is used in the fit, but the model reproduces interfacial width and tension γ lines fit symbols -MD Liquid-liquid Interface Small dip in ρ at interface big change in free energy and surface tension 13
14 Fit near coexistence line reproduces Laplace pressure Δp=γ/R Molecular Dynamics Lattice Boltzmann R 1.5 Δp Denniston & Robbins, Phys. Rev. E 69, (003). a/r Dp ÿ1 0 0 s t êm sêr Laplace pressure (sharp interface) prediction Using models like those of Chen, Jasnow and Vinals: an unphysical mode driven by numerical truncation error. Removal of truncation error removes driving, but mode is still present and is easily driven in contact line motion. Δ = γ μ φ R = Δp Recall, isotropic part of pressure was: ( ρμ Φ ) P0 = kbt ψ ρ μφ In a system taking all parameters from molecular dynamics simulation the unphysical mode is gone. Spurious velocities from Numerical truncation errors ~
15 Lattice Boltzmann simulation of contact line motion Concentration and velocity field (same velocity field at top and bottom) Linking Atomistic and Continuum Regions Three overlap regions where solve both continuum and MD Outermost Continuum solution gives MD boundary condition Innermost MD gives continuum boundary condition Middle Two solutions equilibrate independently Fluids: Apply boundary conditions to velocities Solids: Apply boundary conditions to displacements Streamlines in L~0.3mm channel with moving top wall. Atomistic solution in small area (green) removes continuum singularity Model contact region atomistically, elastic deformations with finite-elements, constrain deformations in overlap region 15
16 Hybrid Algorithm Applied to Fluids Continuum: Incompressible Navier-Stokes (Projection method) Atomistic: Molecular dynamics of Lennard-Jones atoms, no-slip Potential: U(r) =4ε((σ/r) 1 -(σ/r) 6 ] ; Units ε, σ N 1 J MD Continuum u J = vi N Continuum MD 1 N NJ J i= 1 J i= 1 v = u ( t) i J NJ 1 DuJ mx && i = Fi Fk + m, NJ k= 1 Dt LJ Fi = Vik. x Potential confines particles at y 3 Insert/remove number of particles equal to net flux k F y Particle Confinement and Mass Flux External force F for y y y3 ( y y) = αp0σ 1 ( y y ) /( y 3 y ) Continuum Maintain a mass flux by introducing particles n(x,t) near y = y mn ( x, t) = Aρu y ( x, y, t) Δt MD Overlap y 3 y y 1 Langevin thermostat for y < y < y3 : LJ V ( rij ) m&& yi = mγy& i + ζ i x j i i ζ ( t) ζ ( t ) = δ δ ( t t ) k TΓ i j ij B MD Werder et al. J. Comp. Phys. 05, 373 (005) claim artifacts, but use different implementation y 0 16
17 Schematic of simulation Moving Wall Continuum Dynamic Couette Flow U Hybrid solution (symbols) tracks full continuum (lines) as a function of time after motion starts Overlap MD Still Wall Δt FE =40Δt MD X. B. Nie, S. Y. Chen, W. N. E and M. O. Robbins, J. Fluid Mech Singular Cavity Flow Continuum approach: Navier-Stokes + no-slip boundary condition (bc) Usually phenomenological no-slip bc has little effect at large scales r Corner flow Molecular scale influences macroscopic forces No-slip boundary condition is discontinuous at corners a, b Stress diverges as 1/r Log divergence in total force on wall y x Only need atomic information near corners Use hybrid method that treats bulk with continuum Navier-Stokes equations, corners with MD 17
18 Coupling in Overlap Region MD Navier Stokes Mean atomic velocity gives boundary condition to NS eqs. Continuum MD 1) Average tangential MD velocity in shadowed bins forced to NS value: Fi 1 && xi = m mn N J Du Fj + Dt J,Fj = x J j= 1 k ) Normal MD velocity constrained by matching mass flux at boundary V Have tested: Agrees with pure MD calculations. Independent of continuum grid 1, 3 and 6σ and specific set of constrained velocities (within MD noise) X.B. Nie, S.Y. Chen and M. R. Robbins, Physics of Fluids 004. LJ jk NS and Hybrid Velocities Near Corners hybrid Effect like slip BC on scale S S is larger of ~σ and U/0.1σ/τ discreteness shear-thinning Hard to use effective Navier BC: spatially varying, nonlinear 18
19 Treating Large Range of Length Scales Problem: Size of atomistic region independent of system size L BUT time to equilibrate NS flow field grows with L. Initial approach limited to L~0.1μm. Solution: Multigrid and time approach Integrate to steady state at each scale with optimum time step. Iterate between scales till self-consistent (~10 times). Result: Size limited only by onset of non-steady, turbulent flow Show results for 0.3mm cavities. > 10 orders of magnitude faster than fully atomistic ~ 0 minutes per iteration Use average over 16 MD representations to accelerate Schematic of Local Refinement M Coarse ->Fine: Prolongation. Fine->Coarse: Restriction. M Flow at each scale reaches steady state at its own characteristic time 19
20 Multiscale Solution for Re=6400 (U=0.068σ/τ) 0.1mm Ten grid levels, largest 56x56, others 64x64, smallest mesh 0.95σ Dashed lines: the regions expanded in successive plots. Final plot MD region Stokes equations bottom corners self-similar under mag. by ~16 (red arrows) This scaling is cut off by atomic structure. Computational time saving more than over fully atomistic. Stress along the moving wall Three regions contribute to force F: Atomistic, Stokes, high Re Re=ρUL/μ Shear stress on wall Re=6400 U=0.7 σ/τ Re= U=0.7 σ/τ Breakdown of Stokes for r<s atomistic or r>r I μ/ρu inertial Little change for r < R I as increase Re by increasing L Large r contribution gives change in F for fixed U, atomic props. 0
21 Total Force on the Moving Wall Re - only parameter in continuum theory Find strong variation with U at fixed Re, atomic model Re=400 U (σ/τ) F μu f Stokes f Re = fs + fstokes + f Re r<s S<r<R I R I <r f S = 4.3 8π = ln( R / ) I S π = Re f S given by assumption that stress saturates at S S= Ut LJ ; f Re is phenomenological fit Summary for Multiscale Fluid Flow We have developed a multiscale hybrid method that can simulate a macro-length scale flow while still resolving the atomistic structure in a small region. Treats mass and heat flux The ability to resolve the stress on all scales enables the first calculation of the drag force on the moving wall in cavity flow. The force depends on three dimensionless numbers: Re=ρUL/μ, I=Uτ/σ and R m =S/r I =ρus/μ Recently extended to time extrapolation X. Nie, S. Chen and M. O. Robbins, Physics of Fluids 16, (004). X. B. Nie, S. Y. Chen, W. N. E and M. O. Robbins, J. Fluid Mech. 500, (004). X. Nie, M. O. Robbins and S. Chen, Phys. Rev. Lett. 96, (006). B. Luan, S. Hyun, J. F. Molinari, N. Bernstein, and Mark O. Robbins, Phys. Rev. E74, (006). J. Liu, S. Chen, X. Nie and M. O. Robbins, J. Comp. Phys. 7, (007). 1
22 Time Extrapolation Scheme J. Liu, S. Chen, X. Nie and M. O. Robbins, Commun. Comp. Phys. 4,179 (008). Couette flow driven by oscillating wall Δx=Δy=15.6σ, Δt FD =10τ, Δt 1 =Δt =0.5 τ Δx=Δy=31.3σ, Δt FD =50τ, Δt 1 =Δt =.5 τ
23 See Delay If No Time Extrapolation 3
Modeling the combined effect of surface roughness and shear rate on slip flow of simple fluids
Modeling the combined effect of surface roughness and shear rate on slip flow of simple fluids Anoosheh Niavarani and Nikolai Priezjev www.egr.msu.edu/~niavaran November 2009 A. Niavarani and N.V. Priezjev,
More informationthe moving contact line
Molecular hydrodynamics of the moving contact line Tiezheng Qian Mathematics Department Hong Kong University of Science and Technology in collaboration with Ping Sheng (Physics Dept, HKUST) Xiao-Ping Wang
More informationBoundary Conditions for the Moving Contact Line Problem. Abstract
Boundary Conditions for the Moving Contact Line Problem Weiqing Ren Courant Institute of Mathematical Sciences, New York University, New York, NY 10012, USA Weinan E Department of Mathematics and PACM,
More informationMD simulation of methane in nanochannels
MD simulation of methane in nanochannels COCIM, Arica, Chile M. Horsch, M. Heitzig, and J. Vrabec University of Stuttgart November 6, 2008 Scope and structure Molecular model for graphite and the fluid-wall
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 informationDensity Functional Modeling of Nanocrystalline Materials
Density Functional Modeling of Nanocrystalline Materials A new approach for modeling atomic scale properties in materials Peter Stefanovic Supervisor: Nikolas Provatas 70 / Part 1-7 February 007 Density
More informationAnalytical and Numerical Study of Coupled Atomistic-Continuum Methods for Fluids
Analytical and Numerical Study of Coupled Atomistic-Continuum Methods for Fluids Weiqing Ren Courant Institute, NYU Joint work with Weinan E Multiscale modeling for two types of problems: Complex fluids
More informationAdvantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations
The Hilltop Review Volume 7 Issue 1 Winter 2014 Article 10 December 2014 Advantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations Tai-Hsien Wu Western Michigan University
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 informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
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 informationFluid Mechanics II Viscosity and shear stresses
Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small
More informationMicrofluidics 1 Basics, Laminar flow, shear and flow profiles
MT-0.6081 Microfluidics and BioMEMS Microfluidics 1 Basics, Laminar flow, shear and flow profiles 11.1.2017 Ville Jokinen Outline of the next 3 weeks: Today: Microfluidics 1: Laminar flow, flow profiles,
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 informationc 2003 International Press
COMM. MATH. SCI. Vol. 1, No. 2, pp. 333 341 c 23 International Press GENERALIZED NAVIER BOUNDARY CONDITION FOR THE MOVING CONTACT LINE TIEZHENG QIAN, XIAO-PING WANG AND PING SHENG Abstract. From molecular
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 informationNumerical modeling of sliding contact
Numerical modeling of sliding contact J.F. Molinari 1) Atomistic modeling of sliding contact; P. Spijker, G. Anciaux 2) Continuum modeling; D. Kammer, V. Yastrebov, P. Spijker pj ICTP/FANAS Conference
More informationThis article was originally published in a journal published by Elsevier, and the attached copy is provided by Elsevier for the author s benefit and for the benefit of the author s institution, for non-commercial
More informationDissipative Particle Dynamics: Foundation, Evolution and Applications
Dissipative Particle Dynamics: Foundation, Evolution and Applications Lecture 4: DPD in soft matter and polymeric applications George Em Karniadakis Division of Applied Mathematics, Brown University &
More informationContinuum Model of Avalanches in Granular Media
Continuum Model of Avalanches in Granular Media David Chen May 13, 2010 Abstract A continuum description of avalanches in granular systems is presented. The model is based on hydrodynamic equations coupled
More informationFluid Mechanics Theory I
Fluid Mechanics Theory I Last Class: 1. Introduction 2. MicroTAS or Lab on a Chip 3. Microfluidics Length Scale 4. Fundamentals 5. Different Aspects of Microfluidcs Today s Contents: 1. Introduction to
More informationHybrid Atomistic-Continuum Methods for Dense Liquids
Hybrid Atomistic-Continuum Methods for Dense Liquids Matej Praprotnik Laboratory for Molecular Modeling National Institute of Chemistry Ljubljana Slovenia TUD Autumn School 2009, Darmstadt, September 24-25,
More informationCoupling Atomistic and Continuum Hydrodynamics
Coupling Atomistic and Continuum Hydrodynamics Matej Praprotnik praprot@cmm.ki.si Laboratory for Molecular Modeling National Institute of Chemistry Ljubljana, Slovenia & Department of Physics Faculty of
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 information1.3 Molecular Level Presentation
1.3.1 Introduction A molecule is the smallest chemical unit of a substance that is capable of stable, independent existence. Not all substances are composed of molecules. Some substances are composed of
More informationMolecular Hydrodynamics of the Moving Contact Line in Two-phase Immiscible Flows
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 1, No. 1, pp. 1-52 Commun. Comput. Phys. February 26 Molecular Hydrodynamics of the Moving Contact Line in Two-phase Immiscible Flows Tiezheng Qian 1,, Xiao-Ping
More informationThe effective slip length and vortex formation in laminar flow over a rough surface
The effective slip length and vortex formation in laminar flow over a rough surface Anoosheh Niavarani and Nikolai V. Priezjev Movies and preprints @ http://www.egr.msu.edu/~niavaran A. Niavarani and N.V.
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 informationExplaining and modelling the rheology of polymeric fluids with the kinetic theory
Explaining and modelling the rheology of polymeric fluids with the kinetic theory Dmitry Shogin University of Stavanger The National IOR Centre of Norway IOR Norway 2016 Workshop April 25, 2016 Overview
More informationInterfacial dynamics
Interfacial dynamics Interfacial dynamics = dynamic processes at fluid interfaces upon their deformation Interfacial rheological properties: elasticity, viscosity, yield stress, Relation between macroscopic
More informationA Control Algorithm for Multiscale Simulations of Liquid Water
A Control Algorithm for Multiscale Simulations of Liquid Water Evangelos M. Kotsalis and Petros Koumoutsakos Computational Science, ETH Zürich, CH-8092, Switzerland petros@ethz.ch Abstract. We present
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 informationCESSATION OF VISCOPLASTIC POISEUILLE FLOW IN A RECTANGULAR DUCT WITH WALL SLIP
8 th GRACM International Congress on Computational Mechanics Volos, 2 July 5 July 205 CESSATION OF VISCOPLASTIC POISEUILLE FLOW IN A RECTANGULAR DUCT WITH WALL SLIP Yiolanda Damianou, George Kaoullas,
More informationCoupled continuum hydrodynamics and molecular dynamics method for multiscale simulation
Coupled continuum hydrodynamics and molecular dynamics method for multiscale simulation Matthew K. BORG 1,, Duncan A. LOCKERBY 2, Jason M. REESE 1 * Corresponding author: Tel.: +44() 141 548 4386; Email:
More informationThe Polymers Tug Back
Tugging at Polymers in Turbulent Flow The Polymers Tug Back Jean-Luc Thiffeault http://plasma.ap.columbia.edu/ jeanluc Department of Applied Physics and Applied Mathematics Columbia University Tugging
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 informationAn Introduction to Theories of Turbulence. James Glimm Stony Brook University
An Introduction to Theories of Turbulence James Glimm Stony Brook University Topics not included (recent papers/theses, open for discussion during this visit) 1. Turbulent combustion 2. Turbulent mixing
More informationDynamics of Glaciers
Dynamics of Glaciers McCarthy Summer School 01 Andy Aschwanden Arctic Region Supercomputing Center University of Alaska Fairbanks, USA June 01 Note: This script is largely based on the Physics of Glaciers
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More informationExercise: concepts from chapter 10
Reading:, Ch 10 1) The flow of magma with a viscosity as great as 10 10 Pa s, let alone that of rock with a viscosity of 10 20 Pa s, is difficult to comprehend because our common eperience is with s like
More informationHydrodynamics, Thermodynamics, and Mathematics
Hydrodynamics, Thermodynamics, and Mathematics Hans Christian Öttinger Department of Mat..., ETH Zürich, Switzerland Thermodynamic admissibility and mathematical well-posedness 1. structure of equations
More informationSuperhydrophobic surfaces: stability of the Cassie-Baxter state and its effect on liquid water slippage
Superhydrophobic surfaces: stability of the Cassie-Baxter state and its effect on liquid water slippage Mauro Chinappi Center for Life Nano Science IIT@Sapienza Talk outlines Talk 1: Solid Molecular Dynamics
More informationMolecular Dynamics. Timothy O Sullivan Supervised by Naida Lacevic & John Sader The University of Melbourne
Molecular Dynamics Timothy O Sullivan Supervised by Naida Lacevic & John Sader The University of Melbourne 1 1 Introduction Molecular Dynamics (MD) is a method of simulating systems of nanoscale volumes,
More information4. The Green Kubo Relations
4. The Green Kubo Relations 4.1 The Langevin Equation In 1828 the botanist Robert Brown observed the motion of pollen grains suspended in a fluid. Although the system was allowed to come to equilibrium,
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 informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
More informationChapter 1: Basic Concepts
What is a fluid? A fluid is a substance in the gaseous or liquid form Distinction between solid and fluid? Solid: can resist an applied shear by deforming. Stress is proportional to strain Fluid: deforms
More informationExcerpt from the Proceedings of the COMSOL Users Conference 2006 Boston
Using Comsol Multiphysics to Model Viscoelastic Fluid Flow Bruce A. Finlayson, Professor Emeritus Department of Chemical Engineering University of Washington, Seattle, WA 98195-1750 finlayson@cheme.washington.edu
More informationTurbulentlike Quantitative Analysis on Energy Dissipation in Vibrated Granular Media
Copyright 011 Tech Science Press CMES, vol.71, no., pp.149-155, 011 Turbulentlike Quantitative Analysis on Energy Dissipation in Vibrated Granular Media Zhi Yuan Cui 1, Jiu Hui Wu 1 and Di Chen Li 1 Abstract:
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 informationNumber of pages in the question paper : 05 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
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 informationDerivation of continuum models for the moving contact line problem based on thermodynamic principles. Abstract
Derivation of continuum models for the moving contact line problem based on thermodynamic principles Weiqing Ren Courant Institute of Mathematical Sciences, New York University, New York, NY 002, USA Weinan
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 informationHeterogeneous multiscale method for the modeling of complex fluids and micro-fluidics
Journal of Computational Physics xxx (4) xxx xxx www.elsevier.com/locate/jcp Heterogeneous multiscale method for the modeling of complex fluids and micro-fluidics Weiqing Ren *, Weinan E Department of
More informationNumber of pages in the question paper : 06 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More informationDiffusive Transport Enhanced by Thermal Velocity Fluctuations
Diffusive Transport Enhanced by Thermal Velocity Fluctuations Aleksandar Donev 1 Courant Institute, New York University & Alejandro L. Garcia, San Jose State University John B. Bell, Lawrence Berkeley
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 comparative study between dissipative particle dynamics and molecular dynamics for simple- and complex-geometry flows
THE JOURNAL OF CHEMICAL PHYSICS 123, 104107 2005 A comparative study between dissipative particle dynamics and molecular dynamics for simple- and complex-geometry flows Eric E. Keaveny, Igor V. Pivkin,
More informationProjektbericht Applikationsserver Phoenix
Guido Fuchs Institut für Strömungsmechanik und Wärmeübertragung, Tu-Wien guido.fuchs@.tuwien.ac.at Projektbericht Applikationsserver Phoenix MOLECULAR DYNAMICS SIMULATION OF MARANGONI FLOW IN A CAVITY
More informationIn the name of Allah the most beneficent the most merciful
In the name of Allah the most beneficent the most merciful Transient flows of Maxwell fluid with slip conditions Authors Dr. T. Hayat & Sahrish Zaib Introduction Non-Newtonian Newtonian fluid Newtonian
More informationarxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 28 Jun 2001
arxiv:cond-mat/665v [cond-mat.mtrl-sci] 28 Jun 2 Matching Conditions in -Continuum Modeling of Materials Weinan E and Zhongyi Huang 2 Department of Mathematics and PACM, Princeton University and School
More informationCHAPTER 9. Microscopic Approach: from Boltzmann to Navier-Stokes. In the previous chapter we derived the closed Boltzmann equation:
CHAPTER 9 Microscopic Approach: from Boltzmann to Navier-Stokes In the previous chapter we derived the closed Boltzmann equation: df dt = f +{f,h} = I[f] where I[f] is the collision integral, and we have
More informationNDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.
CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo
More informationSIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED
1 SIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED WITH LATERAL FRICTIONAL EFFECTS High slope results and comparison with experimental data Nathan Martin1, Ioan Ionescu2, Anne Mangeney1,3 François
More informationModelling of dispersed, multicomponent, multiphase flows in resource industries. Section 3: Examples of analyses conducted for Newtonian fluids
Modelling of dispersed, multicomponent, multiphase flows in resource industries Section 3: Examples of analyses conducted for Newtonian fluids Globex Julmester 017 Lecture # 04 July 017 Agenda Lecture
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 informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationHydrodynamics. l. mahadevan harvard university
Thermodynamics / statistical mechanics : Hydrodynamics N>>1; large number of interacting particles! Hydrodynamics : L >> l; T >> t; long wavelength, slow time - average over (some) microscopic length and
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 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 informationTime dependent slip flow of a micropolar fluid between two parallel plates through state space approach
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 1, Number (016), pp. 1709 17 Research India Publications http://www.ripublication.com/gjpam.htm Time dependent slip flow of a micropolar
More informationQuick Recapitulation of Fluid Mechanics
Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change
More informationMicro-Scale Gas Transport Modeling
Micro-Scale Gas Transport Modeling Continuum & Slip Flow Regimes: Navier-Stokes Equations Slip Boundary Conditions U g U w = ( σ ) σ U Kn n Slip, Transitional & Free Molecular: Direct Simulation Monte
More informationTable of Contents. Preface... xiii
Preface... xiii PART I. ELEMENTS IN FLUID MECHANICS... 1 Chapter 1. Local Equations of Fluid Mechanics... 3 1.1. Forces, stress tensor, and pressure... 4 1.2. Navier Stokes equations in Cartesian coordinates...
More informationStatistical Mechanics of Active Matter
Statistical Mechanics of Active Matter Umberto Marini Bettolo Marconi University of Camerino, Italy Naples, 24 May,2017 Umberto Marini Bettolo Marconi (2017) Statistical Mechanics of Active Matter 2017
More informationCoarse-graining limits in open and wall-bounded dissipative particle dynamics systems
THE JOURNAL OF CHEMICAL PHYSICS 124, 184101 2006 Coarse-graining limits in open and wall-bounded dissipative particle dynamics systems Igor V. Pivkin and George E. Karniadakis a Division of Applied Mathematics,
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
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 informationViscoelasticity. Basic Notions & Examples. Formalism for Linear Viscoelasticity. Simple Models & Mechanical Analogies. Non-linear behavior
Viscoelasticity Basic Notions & Examples Formalism for Linear Viscoelasticity Simple Models & Mechanical Analogies Non-linear behavior Viscoelastic Behavior Generic Viscoelasticity: exhibition of both
More informationarxiv:comp-gas/ v1 28 Apr 1993
Lattice Boltzmann Thermohydrodynamics arxiv:comp-gas/9304006v1 28 Apr 1993 F. J. Alexander, S. Chen and J. D. Sterling Center for Nonlinear Studies and Theoretical Division Los Alamos National Laboratory
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 informationUnsteady Flow of a Newtonian Fluid in a Contracting and Expanding Pipe
Unsteady Flow of a Newtonian Fluid in a Contracting and Expanding Pipe T S L Radhika**, M B Srinivas, T Raja Rani*, A. Karthik BITS Pilani- Hyderabad campus, Hyderabad, Telangana, India. *MTC, Muscat,
More informationAn electrokinetic LB based model for ion transport and macromolecular electrophoresis
An electrokinetic LB based model for ion transport and macromolecular electrophoresis Raffael Pecoroni Supervisor: Michael Kuron July 8, 2016 1 Introduction & Motivation So far an mesoscopic coarse-grained
More informationContents. I Introduction 1. Preface. xiii
Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................
More informationCHAPTER 1 Fluids and their Properties
FLUID MECHANICS Gaza CHAPTER 1 Fluids and their Properties Dr. Khalil Mahmoud ALASTAL Objectives of this Chapter: Define the nature of a fluid. Show where fluid mechanics concepts are common with those
More information2 Law of conservation of energy
1 Newtonian viscous Fluid 1 Newtonian fluid For a Newtonian we already have shown that σ ij = pδ ij + λd k,k δ ij + 2µD ij where λ and µ are called viscosity coefficient. For a fluid under rigid body motion
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 11 Dec 2002
arxiv:cond-mat/0212257v1 [cond-mat.stat-mech] 11 Dec 2002 International Journal of Modern Physics B c World Scientific Publishing Company VISCOUS FINGERING IN MISCIBLE, IMMISCIBLE AND REACTIVE FLUIDS PATRICK
More informationPhysical Modeling of Multiphase flow. Boltzmann method
with lattice Boltzmann method Exa Corp., Burlington, MA, USA Feburary, 2011 Scope Scope Re-examine the non-ideal gas model in [Shan & Chen, Phys. Rev. E, (1993)] from the perspective of kinetic theory
More informationcontact line dynamics
contact line dynamics part 2: hydrodynamics dynamic contact angle? lubrication: Cox-Voinov theory maximum speed for instability corner shape? dimensional analysis: speed U position r viscosity η pressure
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 informationMesoscale fluid simulation of colloidal systems
Mesoscale fluid simulation of colloidal systems Mingcheng Yang Institute of Physics, CAS Outline (I) Background (II) Simulation method (III) Applications and examples (IV) Summary Background Soft matter
More informationChristel Hohenegger A simple model for ketchup-like liquid, its numerical challenges and limitations April 7, 2011
Notes by: Andy Thaler Christel Hohenegger A simple model for ketchup-like liquid, its numerical challenges and limitations April 7, 2011 Many complex fluids are shear-thinning. Such a fluid has a shear
More informationOn the displacement of two immiscible Stokes fluids in a 3D Hele-Shaw cell
On the displacement of two immiscible Stokes fluids in a 3D Hele-Shaw cell Gelu Paşa Abstract. In this paper we study the linear stability of the displacement of two incompressible Stokes fluids in a 3D
More informationSimulation of mixing of heterogeneous HE components
Chapter Simulation of mixing of heterogeneous HE components The majority on high explosives (HEs) used are blend ones. Properties of components differ that produces interaction on the grain scale (mesoprocesses).
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 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 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 information