Conservation Laws of Surfactant Transport Equations
|
|
- Angelica Shonda Chambers
- 5 years ago
- Views:
Transcription
1 Conservation Laws of Surfactant Transport Equations Alexei Cheviakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada Winter 2011 CMS Meeting Dec. 10, 2011 A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
2 Outline 1 Surfactants Brief Overview Applications Surfactant Transport Equations 2 Conservation Laws Definition and Examples Applications 3 Direct Construction Method for Conservation Laws The Idea Completeness 4 Conservation Laws of Examples of Surfactant Dynamics Equations The Convection Case The Convection-Diffusion Case 5 Conclusions A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
3 Outline 1 Surfactants Brief Overview Applications Surfactant Transport Equations 2 Conservation Laws Definition and Examples Applications 3 Direct Construction Method for Conservation Laws The Idea Completeness 4 Conservation Laws of Examples of Surfactant Dynamics Equations The Convection Case The Convection-Diffusion Case 5 Conclusions A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
4 Surfactants - Brief Overview Surfactants Surfactant = Surface active agent. Act as detergents, wetting agents, emulsifiers, foaming agents, and dispersants. Consist of a hydrophobic group (tail) and a hydrophilic group (head). A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
5 Surfactants - Brief Overview Surfactants Surfactant = Surface active agent. Act as detergents, wetting agents, emulsifiers, foaming agents, and dispersants. Consist of a hydrophobic group (tail) and a hydrophilic group (head). Sodium lauryl sulfate: Saponification: e.g., triglyceride (fat) + base sodium stearate (soap) A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
6 Surfactants - Brief Overview Surfactants Surfactant = Surface active agent. Act as detergents, wetting agents, emulsifiers, foaming agents, and dispersants. Consist of a hydrophobic group (tail) and a hydrophilic group (head). Hydrophilic groups (heads) can have various properties: A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
7 Surfactants - Applications Surfactant molecules adsorb at phase separation interfaces. Stabilization of growth of bubbles / droplets. Creation of emulsions of insoluble substances. Multiple industrial and medical applications. A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
8 Surfactants - Applications Can form micelles, double layers, etc. A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
9 Surfactants - Applications Soap bubbles... A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
10 Surfactant Transport Equations Derivation: Can be derived as a special case of multiphase flows with moving interfaces and contact lines: [Y.Wang, M. Oberlack, 2011] Illustration: Y. Wang, M. Oberlack Fig. 1 The investigated material domain B with its outer boundary B consisting of three-phase subdomains B (i) (i = 1, 2, 3) surrounded by their outer boundaries B (i) and phase interfaces S (i) A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
11 Surfactant Transport Equations (ctd.) n = 0 u Parameters Surfactant concentration c = c(x, t). Flow velocity u(x, t). Two-phase interface: phase separation surface Φ(x, t) = 0. Unit normal: n = Φ Φ. A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
12 Surfactant Transport Equations (ctd.) n = 0 u Surface gradient Surface projection tensor: p ij = δ ij n i n j. Surface gradient operator: s = p = (δ ij n i n j ) x. j Surface Laplacian: s F = (δ ij n i n j ) ( (δ x j ik n i n k ) F ). x k A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
13 Surfactant Transport Equations (ctd.) n = 0 u Governing equations Incompressibility condition: u = 0. Fluid dynamics equations: Euler or Navier-Stokes. Interface transport by the flow: Φ t + u Φ = 0. Surfactant transport equation: c t + u i c x i cn i n j u i x j α(δ ij n i n j ) x j ( (δ ik n i n k ) c ) = 0. x k A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
14 Surfactant Transport Equations (ctd.) n = 0 u Fully conserved form Specific numerical methods (e.g., discontinuous Galerkin) require the system to be written in a fully conserved form. Straightforward for continuity, momentum, and interface transport equations. Can the surfactant transport equation be written in the conserved form? c t + u i c u i cn x i i n j α(δ x j ij n i n j ) ( (δ x j ik n i n k ) c ) = 0. x k A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
15 Outline 1 Surfactants Brief Overview Applications Surfactant Transport Equations 2 Conservation Laws Definition and Examples Applications 3 Direct Construction Method for Conservation Laws The Idea Completeness 4 Conservation Laws of Examples of Surfactant Dynamics Equations The Convection Case The Convection-Diffusion Case 5 Conclusions A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
16 Conservation Laws: Some Definitions Conservation law: Given: some model. Independent variables: x = (t, x, y,...); dependent variables: u = (u, v,...). A conservation law: a divergence expression equal to zero D tθ(x, u,...) + D xψ 1 (x, u,...) + D y Ψ 2 (x, u,...) + = 0. Time-independent: D i Ψ i div x,y,...ψ = 0. Time-dependent: A conserved quantity: D tθ + div x,y,...ψ = 0. D t Θ dv = 0. V A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
17 A Simple PDE Example Example: small string oscillations, 1D wave equation Independent variables: x, t; dependent: u(x, t). u tt = c 2 u xx, c 2 = T /ρ. Conservation of momentum: Divergence expression: D t(ρu t) D x(tu x) = 0. Arises from a multiplier: D t(ρu t) D x(tu x) = Λ(u tt c 2 u xx) = ρ(u tt c 2 u xx) = 0; Conserved quantity: total momentum M = ρu tdx = const. A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
18 A Simple PDE Example Example: small string oscillations, 1D wave equation Independent variables: x, t; dependent: u(x, t). u tt = c 2 u xx, c 2 = T /ρ. Conservation of Energy: ( ρu 2 Divergence expression: D t t 2 + Tu2 x 2 Arises from a multiplier: ( ρu 2 D t t 2 + Tu2 x 2 ) D x(tu tu x) = 0. ) D x(tu tu x) = ρu t(u tt c 2 u xx) = 0; Conserved quantity: total ebergy E = ( ρut Tu2 x 2 ) dx = const. A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
19 Applications of Conservation Laws ODEs Constants of motion. Integration. PDEs Direct physical meaning. Constants of motion. Differential constraints (div B = 0, etc.). Analysis: existence, uniqueness, stability. Nonlocally related PDE systems, exact solutions. Potentials, stream functions, etc.: V = (u, v), div V = u x + v y = 0, { u = Φy, v = Φ x. An infinite number of conservation laws can indicate integrability / linearization. Numerical simulation. A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
20 Outline 1 Surfactants Brief Overview Applications Surfactant Transport Equations 2 Conservation Laws Definition and Examples Applications 3 Direct Construction Method for Conservation Laws The Idea Completeness 4 Conservation Laws of Examples of Surfactant Dynamics Equations The Convection Case The Convection-Diffusion Case 5 Conclusions A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
21 The Idea of the Direct Construction Method Definition The Euler operator with respect to U j : E U j = D U j i U j i + + ( 1) s D i1... D is U j i 1...i s +, j = 1,..., m. Consider a general DE system R σ [u] = R σ (x, u, u,..., k u) = 0, variables x = (x 1,..., x n ), u = u(x) = (u 1,..., u m ). σ = 1,..., N with Theorem The equations E U j F (x, U, U,..., s U) 0, hold for arbitrary U(x) if and only if j = 1,..., m for some functions Ψ i (x, U,...). F (x, U, U,..., s U) D i Ψ i A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
22 The Idea of the Direct Construction Method Consider a general DE system R σ [u] = R σ (x, u, u,..., k u) = 0, σ = 1,..., N with variables x = (x 1,..., x n ), u = u(x) = (u 1,..., u m ). Direct Construction Method [Anco, Bluman (1997,2002)] Specify dependence of multipliers: Λ σ = Λ σ(x, U,...), σ = 1,..., N. Solve the set of determining equations E U j (Λ σr σ ) 0, j = 1,..., m, for arbitrary U(x) (off of solution set!) to find all such sets of multipliers. Find the corresponding fluxes Φ i (x, U,...) satisfying the identity Λ σr σ D i Φ i. Each set of fluxes and multipliers yields a local conservation law D i Φ i (x, u,...) = 0, holding for all solutions u(x) of the given PDE system. A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
23 Completeness Completeness of the Direct Construction Method For the majority of physical DE systems (in particular, all systems in solved form), all conservation laws follow from linear combinations of equations! Λ σr σ D i Φ i. A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
24 Outline 1 Surfactants Brief Overview Applications Surfactant Transport Equations 2 Conservation Laws Definition and Examples Applications 3 Direct Construction Method for Conservation Laws The Idea Completeness 4 Conservation Laws of Examples of Surfactant Dynamics Equations The Convection Case The Convection-Diffusion Case 5 Conclusions A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
25 CLs of the Surfactant Dynamics Equations: The Convection Case Governing equations (α = 0) R 1 = ui x i = 0, R 2 = Φ t + (ui Φ) x i = 0, R 3 = c t + u i c x i cn i n j u i x j = 0. A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
26 CLs of the Surfactant Dynamics Equations: The Convection Case Governing equations (α = 0) Multiplier ansatz R 1 = ui x i = 0, R 2 = Φ t + (ui Φ) x i = 0, R 3 = c t + u i c x i cn i n j u i x j = 0. Λ i = Λ i (t, x, Φ, c, u, Φ, c, u, 2 Φ, 2 c, 2 u). A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
27 CLs of the Surfactant Dynamics Equations: The Convection Case Governing equations (α = 0) Multiplier ansatz R 1 = ui x i = 0, R 2 = Φ t + (ui Φ) x i = 0, R 3 = c t + u i c x i cn i n j u i x j = 0. Λ i = Λ i (t, x, Φ, c, u, Φ, c, u, 2 Φ, 2 c, 2 u). Conservation Law Determining Equations E u j (Λ σ R σ ) = 0, j = 1,..., 3; E Φ (Λ σ R σ ) = 0; E c(λ σ R σ ) = 0. A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
28 CLs of the Surfactant Dynamics Equations: The Convection Case (ctd.) Governing equations (α = 0) R 1 = ui x i = 0, R 2 = Φ t + (ui Φ) x i = 0, R 3 = c t + u i c x i cn i n j u i x j = 0. Principal Result 1 (multipliers) There exist an infinite family of multiplier sets with Λ 3 0, i.e., essentially involving c. Family of conservation laws with Λ 3 = Φ K(Φ, c Φ ). A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
29 CLs of the Surfactant Dynamics Equations: The Convection Case (ctd.) Governing equations (α = 0) R 1 = ui x i = 0, R 2 = Φ t + (ui Φ) x i = 0, R 3 = c t + u i c x i cn i n j u i x j = 0. Principal Result 1 (divergence expressions) Usual form: Material form: G(Φ, c Φ ) + ( ) u i G(Φ, c Φ ) = 0. t x i d G(Φ, c Φ ) = 0. dt A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
30 CLs of the Surfactant Dynamics Equations: The Convection Case (ctd.) Governing equations (α = 0) R 1 = ui x i = 0, R 2 = Φ t + (ui Φ) x i = 0, R 3 = c t + u i c x i cn i n j u i x j = 0. Simplest conservation law with c-dependence Can take G(Φ, c Φ ) = c Φ. t (c Φ ) + ( ) u i c Φ = 0. x i A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
31 The Convection-Diffusion Case Governing equations (α 0) R 1 = ui x i = 0, R 2 = Φ t + (ui Φ) = 0, x i R 3 = c t + u i c x i cn i n j u i x j α(δ ij n i n j ) x j ( (δ ik n i n k ) c ) = 0. x k Conservation Law Determining Equations E u j (Λ σ R σ ) = 0, j = 1,..., 3; E Φ (Λ σ R σ ) = 0; E c(λ σ R σ ) = 0. A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
32 The Convection-Diffusion Case (ctd.) Governing equations (α = 0) R 1 = ui x i = 0, R 2 = Φ t + (ui Φ) = 0, x i R 3 = c t + u i c x i cn i n j u i x j α(δ ij n i n j ) x j ( (δ ik n i n k ) c ) = 0. x k Principal Result 2 (multipliers) Λ 1 = ΦF(Φ) Φ 1 ( x j Λ 2 = F(Φ) Φ 1 ( x j Λ 3 = F(Φ) Φ, ( c Φ ) ) 2 Φ cn x j i n j, x i x ( j c Φ ) ) 2 Φ cn x j i n j, x i x j A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
33 The Convection-Diffusion Case (ctd.) Governing equations (α = 0) R 1 = ui x i = 0, R 2 = Φ t + (ui Φ) = 0, x i R 3 = c t + u i c x i cn i n j u i x j α(δ ij n i n j ) x j ( (δ ik n i n k ) c ) = 0. x k Principal Result 2 (divergence expressions) An infinite family of conservation laws: (c F(Φ) Φ ) + ( ) A i F(Φ) Φ = 0, t x i where A i = cu i α ( (δ ik n i n k ) c ), i = 1, 2, 3, x k and F is an arbitrary sufficiently smooth function. A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
34 The Convection-Diffusion Case (ctd.) Governing equations (α = 0) R 1 = ui x i = 0, R 2 = Φ t + (ui Φ) = 0, x i R 3 = c t + u i c x i cn i n j u i x j α(δ ij n i n j ) x j ( (δ ik n i n k ) c ) = 0. x k Simplest conservation law with c-dependence Can take F(Φ) = 1: (c Φ ) + ( ) A i Φ = 0. t x i A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
35 Conclusions Some conclusions Conservation laws can be obtained systematically through the Direct Construction Method, which employs multipliers and Euler operators. The method is implemented in a symbolic package GeM for Maple. Surfactant dynamics equations can be written in a fully conserved form. Infinite families of c-dependent conservation laws for cases with and without diffusion. A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
36 Conclusions Some conclusions Conservation laws can be obtained systematically through the Direct Construction Method, which employs multipliers and Euler operators. The method is implemented in a symbolic package GeM for Maple. Surfactant dynamics equations can be written in a fully conserved form. Infinite families of c-dependent conservation laws for cases with and without diffusion. Open problems Higher order conservation laws? Can we get more by including fluid dynamics equations? A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
37 Some references Anco, S.C. and Bluman, G.W. (1997). Direct construction of conservation laws from field equations. Phys. Rev. Lett. 78, Anco, S.C. and Bluman, G.W. (2002). Direct construction method for conservation laws of partial differential equations. Part I: Examples of conservation law classifications. Eur. J. Appl. Math. 13, Bluman, G.W., Cheviakov, A.F., and Anco, S.C. (2010). Applications of Symmetry Methods to Partial Differential Equations. Springer: Applied Mathematical Sciences, Vol Cheviakov, A.F. (2007). GeM software package for computation of symmetries and conservation laws of differential equations. Comput. Phys. Commun. 176, Kallendorf, C., Cheviakov, A.F., Oberlack, M., and Wang, Y. Conservation Laws of Surfactant Transport Equations. Submitted, A. Cheviakov (Math&Stat) Cons. laws for Surfactant Transport Equations 10 Dec / 22
Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations.
Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations. Alexey Shevyakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon,
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
More informationConservation Laws of Fluid Dynamics Models
Conservation Laws of Fluid Dynamics Models Prof. Alexei Cheviakov (Alt. English spelling: Alexey Shevyakov) Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada June
More informationA Recursion Formula for the Construction of Local Conservation Laws of Differential Equations
A Recursion Formula for the Construction of Local Conservation Laws of Differential Equations Alexei Cheviakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada December
More informationConservation Laws for Nonlinear Equations: Theory, Computation, and Examples
Conservation Laws for Nonlinear Equations: Theory, Computation, and Examples Alexei Cheviakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada Seminar February 2015
More informationNatural States and Symmetry Properties of. Two-Dimensional Ciarlet-Mooney-Rivlin. Nonlinear Constitutive Models
Natural States and Symmetry Properties of Two-Dimensional Ciarlet-Mooney-Rivlin Nonlinear Constitutive Models Alexei Cheviakov, Department of Mathematics and Statistics, Univ. Saskatchewan, Canada Jean-François
More informationSymbolic Computation of Nonlocal Symmetries and Nonlocal Conservation Laws of Partial Differential Equations Using the GeM Package for Maple
Symbolic Computation of Nonlocal Symmetries and Nonlocal Conservation Laws of Partial Differential Equations Using the GeM Package for Maple Alexei F. Cheviakov Abstract The use of the symbolic software
More informationFree energy concept Free energy approach LBM implementation Parameters
BINARY LIQUID MODEL A. Kuzmin J. Derksen Department of Chemical and Materials Engineering University of Alberta Canada August 22,2011 / LBM Workshop OUTLINE 1 FREE ENERGY CONCEPT 2 FREE ENERGY APPROACH
More informationA comparison of conservation law construction approaches for the two-dimensional incompressible Mooney-Rivlin hyperelasticity model
JOURNAL OF MATHEMATICAL PHYSICS 56, 505 (05) A comparison of conservation law construction approaches for the two-dimensional incompressible Mooney-Rivlin hyperelasticity model A. F. Cheviakov and S. St.
More informationConstruction of Conservation Laws: How the Direct Method Generalizes Noether s Theorem
Proceedings of 4th Workshop Group Analysis of Differential Equations & Integrability 2009, 1 23 Construction of Conservation Laws: How the Direct Method Generalizes Noether s Theorem George W. BLUMAN,
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationThe Orchestra of Partial Differential Equations. Adam Larios
The Orchestra of Partial Differential Equations Adam Larios 19 January 2017 Landscape Seminar Outline 1 Fourier Series 2 Some Easy Differential Equations 3 Some Not-So-Easy Differential Equations Outline
More informationNonlinear Equations for Finite-Amplitude Wave Propagation in Fiber-Reinforced Hyperelastic Media
Nonlinear Equations for Finite-Amplitude Wave Propagation in Fiber-Reinforced Hyperelastic Media Alexei F. Cheviakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada
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 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 informationREGULARIZATION AND BOUNDARY CONDITIONS FOR THE 13 MOMENT EQUATIONS
1 REGULARIZATION AND BOUNDARY CONDITIONS FOR THE 13 MOMENT EQUATIONS HENNING STRUCHTRUP ETH Zürich, Department of Materials, Polymer Physics, CH-8093 Zürich, Switzerland (on leave from University of Victoria,
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 informationThe Euler Equations and Non-Local Conservative Riccati Equations
The Euler Equations and Non-Local Conservative Riccati Equations Peter Constantin Department of Mathematics The University of Chicago November 8, 999 The purpose of this brief note is to present an infinite
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 informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
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 informationModeling using conservation laws. Let u(x, t) = density (heat, momentum, probability,...) so that. u dx = amount in region R Ω. R
Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...) so that u dx = amount in region R Ω. R Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...)
More informationA symmetry-based method for constructing nonlocally related partial differential equation systems
A symmetry-based method for constructing nonlocally related partial differential equation systems George W. Bluman and Zhengzheng Yang Citation: Journal of Mathematical Physics 54, 093504 (2013); doi:
More informationScientific Computing I
Scientific Computing I Session 11: Basics of Partial Differential Equations Tobias Weinzierl Winter 2010/2011 Session 11: Basics of Partial Differential Equations, Winter 2010/2011 1 Motivation Instead
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Partial Differential Equations (Lecture 1, Week 1) Markus Schmuck Department of Mathematics and Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh
More informationTRANSPORT IN POROUS MEDIA
1 TRANSPORT IN POROUS MEDIA G. ALLAIRE CMAP, Ecole Polytechnique 1. Introduction 2. Main result in an unbounded domain 3. Asymptotic expansions with drift 4. Two-scale convergence with drift 5. The case
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 informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) 1D Heat Equation: Derivation Current Semester 1 / 19 Introduction The derivation of the heat
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 informationThe CG1-DG2 method for conservation laws
for conservation laws Melanie Bittl 1, Dmitri Kuzmin 1, Roland Becker 2 MoST 2014, Germany 1 Dortmund University of Technology, Germany, 2 University of Pau, France CG1-DG2 Method - Motivation hp-adaptivity
More informationDiffusion - The Heat Equation
Chapter 6 Diffusion - The Heat Equation 6.1 Goal Understand how to model a simple diffusion process and apply it to derive the heat equation in one dimension. We begin with the fundamental conservation
More informationSymmetry Reductions of (2+1) dimensional Equal Width. Wave Equation
Authors: Symmetry Reductions of (2+1) dimensional Equal Width 1. Dr. S. Padmasekaran Wave Equation Asst. Professor, Department of Mathematics Periyar University, Salem 2. M.G. RANI Periyar University,
More informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
More informationModule 4: "Surface Thermodynamics" Lecture 21: "" The Lecture Contains: Effect of surfactant on interfacial tension. Objectives_template
The Lecture Contains: Effect of surfactant on interfacial tension file:///e /courses/colloid_interface_science/lecture21/21_1.htm[6/16/2012 1:10:36 PM] Surface Thermodynamics: Roles of Surfactants and
More information0.2. CONSERVATION LAW FOR FLUID 9
0.2. CONSERVATION LAW FOR FLUID 9 Consider x-component of Eq. (26), we have D(ρu) + ρu( v) dv t = ρg x dv t S pi ds, (27) where ρg x is the x-component of the bodily force, and the surface integral is
More informationThis note presents an infinite-dimensional family of exact solutions of the incompressible three-dimensional Euler equations
IMRN International Mathematics Research Notices 2000, No. 9 The Euler Equations and Nonlocal Conservative Riccati Equations Peter Constantin This note presents an infinite-dimensional family of exact solutions
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationWeek 6 Notes, Math 865, Tanveer
Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,
More informationLOCAL WELL-POSEDNESS FOR AN ERICKSEN-LESLIE S PARABOLIC-HYPERBOLIC COMPRESSIBLE NON-ISOTHERMAL MODEL FOR LIQUID CRYSTALS
Electronic Journal of Differential Equations, Vol. 017 (017), No. 3, pp. 1 8. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu LOCAL WELL-POSEDNESS FOR AN ERICKSEN-LESLIE S
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 informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 21 May 2001 with E. A. Spiegel
More informationControl of Interface Evolution in Multi-Phase Fluid Flows
Control of Interface Evolution in Multi-Phase Fluid Flows Markus Klein Department of Mathematics University of Tübingen Workshop on Numerical Methods for Optimal Control and Inverse Problems Garching,
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 informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel
More informationSymmetry Properties of Two-Dimensional Ciarlet-Mooney-Rivlin Constitutive Models in Nonlinear Elastodynamics
Symmetry Properties of Two-Dimensional Ciarlet-Mooney-Rivlin Constitutive Models in Nonlinear Elastodynamics A. F. Cheviakov a, J.-F. Ganghoffer b a Department of Mathematics and Statistics, University
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 informationMagnetized Target Fusion: Insights from Mathematical Modelling
: Insights from Mathematical Modelling, UBC, PhD Candidate - Applied Math Brian Wetton, UBC, Supervisor Michael Ward, UBC, Committee Member Sandra Barsky, General Fusion, Committee Member November 28,
More informationA NONLINEAR GENERALIZATION OF THE CAMASSA-HOLM EQUATION WITH PEAKON SOLUTIONS
A NONLINEAR GENERALIZATION OF THE CAMASSA-HOLM EQUATION WITH PEAKON SOLUTIONS STEPHEN C. ANCO 1, ELENA RECIO 1,2, MARÍA L. GANDARIAS2, MARÍA S. BRUZÓN2 1 department of mathematics and statistics brock
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 information2. Conservation Equations for Turbulent Flows
2. Conservation Equations for Turbulent Flows Coverage of this section: Review of Tensor Notation Review of Navier-Stokes Equations for Incompressible and Compressible Flows Reynolds & Favre Averaging
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationUNIVERSITY OF OSLO. Faculty of Mathematics and Natural Sciences. MEK4300/9300 Viscous flow og turbulence
UNIVERSITY OF OSLO Faculty of Mathematics and Natural Sciences Examination in: Day of examination: Friday 15. June 212 Examination hours: 9. 13. This problem set consists of 5 pages. Appendices: Permitted
More informationexamples of equations: what and why intrinsic view, physical origin, probability, geometry
Lecture 1 Introduction examples of equations: what and why intrinsic view, physical origin, probability, geometry Intrinsic/abstract F ( x, Du, D u, D 3 u, = 0 Recall algebraic equations such as linear
More informationConnecting Euler and Lagrange systems as nonlocally related systems of dynamical nonlinear elasticity
Arch. Mech., 63, 4, pp. 1 20, Warszawa 2011 Connecting Euler and Lagrange systems as nonlocally related systems of dynamical nonlinear elasticity G. BLUMAN 1), J. F. GANGHOFFER 2) 1) Mathematics Department
More informationNSW Higher School Certificate Senior Science 9.2 Lifestyle Chemistry
NSW Higher School Certificate Senior Science 9.2 Lifestyle Chemistry Section 2 Cleaning Products 9.2 Lifestyle Chemistry Section 2 ::: Cleaning Products 9.2.2 A wide range of cleaning products are made
More informationMultidimensional partial differential equation systems: Nonlocal symmetries, nonlocal conservation laws, exact solutions
JOURNAL OF MATHEMATICAL PHYSICS 51, 103522 2010 Multidimensional partial differential equation systems: Nonlocal symmetries, nonlocal conservation laws, exact solutions Alexei F. Cheviakov 1,a and George
More informationDiscretization of PDEs and Tools for the Parallel Solution of the Resulting Systems
Discretization of PDEs and Tools for the Parallel Solution of the Resulting Systems Stan Tomov Innovative Computing Laboratory Computer Science Department The University of Tennessee Wednesday April 4,
More informationFrom a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction
From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy
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 informationNumerical modelling of phase change processes in clouds. Challenges and Approaches. Martin Reitzle Bernard Weigand
Institute of Aerospace Thermodynamics Numerical modelling of phase change processes in clouds Challenges and Approaches Martin Reitzle Bernard Weigand Introduction Institute of Aerospace Thermodynamics
More informationChapter 3 Burgers Equation
Chapter 3 Burgers Equation One of the major challenges in the field of comple systems is a thorough understanding of the phenomenon of turbulence. Direct numerical simulations (DNS) have substantially
More informationThe Euler equations in fluid mechanics
The Euler equations in fluid mechanics Jordan Bell jordan.bell@gmail.com Department of Mathematics, niversity of Toronto April 14, 014 1 Continuity equation Let Ω be a domain in R n and let C (Ω R); perhaps
More informationPotential symmetry and invariant solutions of Fokker-Planck equation in. cylindrical coordinates related to magnetic field diffusion
Potential symmetry and invariant solutions of Fokker-Planck equation in cylindrical coordinates related to magnetic field diffusion in magnetohydrodynamics including the Hall current A. H. Khater a,1,
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 informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
More informationVariable Exponents Spaces and Their Applications to Fluid Dynamics
Variable Exponents Spaces and Their Applications to Fluid Dynamics Martin Rapp TU Darmstadt November 7, 213 Martin Rapp (TU Darmstadt) Variable Exponent Spaces November 7, 213 1 / 14 Overview 1 Variable
More informationIMPACTS OF SIGMA COORDINATES ON THE EULER AND NAVIER-STOKES EQUATIONS USING CONTINUOUS/DISCONTINUOUS GALERKIN METHODS
Approved for public release; distribution is unlimited IMPACTS OF SIGMA COORDINATES ON THE EULER AND NAVIER-STOKES EQUATIONS USING CONTINUOUS/DISCONTINUOUS GALERKIN METHODS Sean L. Gibbons Captain, United
More informationWeak-Strong Uniqueness of the Navier-Stokes-Smoluchowski System
Weak-Strong Uniqueness of the Navier-Stokes-Smoluchowski System Joshua Ballew University of Maryland College Park Applied PDE RIT March 4, 2013 Outline Description of the Model Relative Entropy Weakly
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationPast, present and space-time
Past, present and space-time Arnold Reusken Chair for Numerical Mathematics RWTH Aachen Utrecht, 12.11.2015 Reusken (RWTH Aachen) Past, present and space-time Utrecht, 12.11.2015 1 / 20 Outline Past. Past
More informationUNIVERSITY of LIMERICK
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4607 SEMESTER: Autumn 2012-13 MODULE TITLE: Introduction to Fluids DURATION OF
More informationNotes on theory and numerical methods for hyperbolic conservation laws
Notes on theory and numerical methods for hyperbolic conservation laws Mario Putti Department of Mathematics University of Padua, Italy e-mail: mario.putti@unipd.it January 19, 2017 Contents 1 Partial
More informationMcGill University April 20, Advanced Calculus for Engineers
McGill University April 0, 016 Faculty of Science Final examination Advanced Calculus for Engineers Math 64 April 0, 016 Time: PM-5PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer Student
More informationMathematical Theory of Non-Newtonian Fluid
Mathematical Theory of Non-Newtonian Fluid 1. Derivation of the Incompressible Fluid Dynamics 2. Existence of Non-Newtonian Flow and its Dynamics 3. Existence in the Domain with Boundary Hyeong Ohk Bae
More informationThe Distribution Function
The Distribution Function As we have seen before the distribution function (or phase-space density) f( x, v, t) d 3 x d 3 v gives a full description of the state of any collisionless system. Here f( x,
More informationInfo. No lecture on Thursday in a week (March 17) PSet back tonight
Lecture 0 8.086 Info No lecture on Thursday in a week (March 7) PSet back tonight Nonlinear transport & conservation laws What if transport becomes nonlinear? Remember: Nonlinear transport A first attempt
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) Key Concepts Current Semester 1 / 25 Introduction The purpose of this section is to define
More informationThe Origins of Surface and Interfacial Tension
The Origins of Surface and Interfacial Tension Imbalance of intermolecular forces exists at the liquid-air interface γ la= the surface tension that exists at the liquid-air interface Suppose we have a
More informationGroup analysis, nonlinear self-adjointness, conservation laws, and soliton solutions for the mkdv systems
ISSN 139-5113 Nonlinear Analysis: Modelling Control, 017, Vol., No. 3, 334 346 https://doi.org/10.15388/na.017.3.4 Group analysis, nonlinear self-adjointness, conservation laws, soliton solutions for the
More informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2011 Randall J. LeVeque Applied Mathematics University of Washington January 3, 2011 R.J. LeVeque, University of Washington AMath 574,
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 informationLast time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler:
Lecture 7 18.086 Last time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler: U j,n+1 t U j,n = U j+1,n 2U j,n + U j 1,n x 2 Expected accuracy: O(Δt) in time,
More informationOn Reduction and Q-conditional (Nonclassical) Symmetry
Symmetry in Nonlinear Mathematical Physics 1997, V.2, 437 443. On Reduction and Q-conditional (Nonclassical) Symmetry Roman POPOVYCH Institute of Mathematics of the National Academy of Sciences of Ukraine,
More information2. Conservation of Mass
2 Conservation of Mass The equation of mass conservation expresses a budget for the addition and removal of mass from a defined region of fluid Consider a fixed, non-deforming volume of fluid, V, called
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 informationÉquation de Burgers avec particule ponctuelle
Équation de Burgers avec particule ponctuelle Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France 7 juin 2010 En collaboration avec B. Andreianov, F. Lagoutière et T. Takahashi Nicolas Seguin
More informationCONNECTIONS BETWEEN A CONJECTURE OF SCHIFFER S AND INCOMPRESSIBLE FLUID MECHANICS
CONNECTIONS BETWEEN A CONJECTURE OF SCHIFFER S AND INCOMPRESSIBLE FLUID MECHANICS JAMES P. KELLIHER Abstract. We demonstrate connections that exists between a conjecture of Schiffer s (which is equivalent
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationEntropy and Relative Entropy
Entropy and Relative Entropy Joshua Ballew University of Maryland October 24, 2012 Outline Hyperbolic PDEs Entropy/Entropy Flux Pairs Relative Entropy Weak-Strong Uniqueness Weak-Strong Uniqueness for
More informationEULERIAN DERIVATIONS OF NON-INERTIAL NAVIER-STOKES EQUATIONS
EULERIAN DERIVATIONS OF NON-INERTIAL NAVIER-STOKES EQUATIONS ML Combrinck, LN Dala Flamengro, a div of Armscor SOC Ltd & University of Pretoria, Council of Scientific and Industrial Research & University
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationAn analogue of Rionero s functional for reaction-diffusion equations and an application thereof
Note di Matematica 7, n., 007, 95 105. An analogue of Rionero s functional for reaction-diffusion equations and an application thereof James N. Flavin Department of Mathematical Physics, National University
More informationThe incompressible Navier-Stokes equations in vacuum
The incompressible, Université Paris-Est Créteil Piotr Bogus law Mucha, Warsaw University Journées Jeunes EDPistes 218, Institut Elie Cartan, Université de Lorraine March 23th, 218 Incompressible Navier-Stokes
More informationConservation and dissipation principles for PDEs
Conservation and dissipation principles for PDEs Modeling through conservation laws The notion of conservation - of number, energy, mass, momentum - is a fundamental principle that can be used to derive
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 informationTravelling Wave Solutions and Conservation Laws of Fisher-Kolmogorov Equation
Gen. Math. Notes, Vol. 18, No. 2, October, 2013, pp. 16-25 ISSN 2219-7184; Copyright c ICSRS Publication, 2013 www.i-csrs.org Available free online at http://www.geman.in Travelling Wave Solutions and
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 informationIntroduction to PDEs and Numerical Methods Lecture 1: Introduction
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Lecture 1: Introduction Dr. Noemi Friedman, 28.10.2015. Basic information on the course Course
More information