Part IV: Numerical schemes for the phase-filed model
|
|
- Brent Nelson
- 5 years ago
- Views:
Transcription
1 Part IV: Numerical schemes for the phase-filed model Jie Shen Department of Mathematics Purdue University IMS, Singapore July 29-3, 29
2 The complete set of governing equations Find u, p, (φ, ξ) such that ρ(u t + (u )u) + p = µ( u + t u) λ ( φ φ); u = ; φ t + (u )φ = γ( φ f(φ) + ξ(t)), d φdx = ; dt Ω (ρ, µ) = + φ 2 (ρ, µ ) + φ 2 (ρ 2, µ 2 ). Remark: No energy law for this system. 2
3 Time discretization of the full system ρ nun+ u n µ n D(u n+ ) + p n = expl. nonl. terms; δt ( ρ n ψn+, q) = ( un+ u n, q), δt p n+ = ψ n+ + p n + ν u n+ ; φ n+ φ n γ( C s δt η 2 I)φn+ = γξ k+ + expl. nonl. terms, φ k+ dx = φ k dx; Ω Ω (ρ n+, µ n+ ) = + φn+ 2 (ρ, µ ) + φn+ (ρ 2, µ 2 ). 2 3
4 Summary: Only a sequence of elliptic equations needs to be solved at each time step: The variable-coefficient term (µ n u n+ ) can be approximated by µ u n+ (µ n µ) u n, where µ is some average of µ(φ); 4
5 the same trick can not be applied to ( ρ n ψn+, q) = ( un+ u n, q). δt For problems with relatively small density ratio, one can use the so called Boussinesq approximation: replacing ρ by a constant and accounting the density difference through a gravity force. Volume of φ is conserved for all time. Only limited stability and error analysis were available for the constant density case. 5
6 Second-order accurate schemes can be constructed easily. The above scheme is very easy to implement as the main building blocks are elliptic/fast Poisson solvers. However, two challenges remains: The elliptic solve for the pressure becomes very expensive as the density ratio becomes large. The time step is mainly constraintly by the explict treatment of the surface tension term. 6
7 The new formulation with energy law Let σ = ρ. Find u, p, (φ, ξ) such that σ(σu) t + (ρu )u + (ρu)u + p 2 = µ( u + t u) λ ( φ φ); u = ; φ t + (u )φ = γ( φ f(φ) + ξ(t)), d φdx = ; dt Ω (ρ, µ) = + φ 2 (ρ, µ ) + φ 2 (ρ 2, µ 2 ). 7
8 The gauge-uzawa scheme for the new model σ n+σn+ ũ n+ σ n u n + ρ n (u n )ũ n+ + δt 2 ( (ρn u n ))ũ n+ µ ũ n+ + µ s n + λ ρ n+ ρ n+ =, ũ n+ Ω = ; ( ρ n+ ψn+ ) = ũ n+, n ψ n+ = ; u n+ = ũ n+ + ρ n+ ψn+, s n+ = s n ũ n+ ; 8
9 ρ n+ ρ n + (ũ n+ )ρ n+ γ( ρ n+ g(ρ n+ )) =, δt n ρ n+ Ω =. Thanks to the energy law for the density based phase-field model, we are able to prove that the above scheme is unconditionally stable. In practice, we treat the nonlinear terms explicitly. Second-order version can be constructed. Still need to solve a pressure equation with variable coefficients. 9
10 To avoid solving a pressure equation with variable coefficients, we consider a pressure-stabilized formulation of NSE (Rannacher 92, S. 93): ρ ɛ (u ɛ t + (u ɛ )u ɛ ) = µ u ɛ p ɛ + ρ ɛ f, ρ ɛ t + (u ɛ )ρ ɛ =, with or u ɛ ɛ p ɛ =, u ɛ ɛ p ɛ t =, p ɛ n Ω =, p ɛ t n Ω =. Notice that only a pressure Poisson equation is involved.
11 A new scheme for the density based phase-field model: σ n+σn+ u n+ σ n u n + ρ n (u n )u n+ δt µ n u n+ + p n + λ ρ n+ ρ n+ =, u n+ Ω = ; p n+ = ρ min δt un+, n p n+ Ω = ; ρ n+ ρ n + (u n+ )ρ n+ γ( ρ n+ g(ρ n+ )) =, δt n ρ n+ Ω =.
12 An improved scheme based on pressure-correction: ρ n+ u n+ ρ n u n + ρ n (u n )u n+ + δt 2 ( (ρn u n ))u n+ µ n u n+ + (p n + ψ n ) + λ ρ n+ ρ n+ =, u n+ Ω = ; ψ n+ = χ δt un+, n ψ n+ Ω =, p n+ = p n + ψ n+ ; ρ n+ ρ n δt + (u n+ )ρ n+ γ( ρ n+ g(ρ n+ )) =, n ρ n+ Ω =. 2
13 Several remarks: Second-order schemes based on pressure-correction can be constructed. For problems with large density ratios, the schemes based on pressure-correction are as accurate and much more efficient, when compared with the gauge-uzawa method. 3
14 Moving mesh method r-refinement Goal: redistribute points, through a smooth mapping, to resolve thin interfaces/regions with large gradients. Approach: Given a function u(x; t), it amounts to find a suitable smooth mapping ξ = ξ(x; t) such that v(ξ; t) := u(x(ξ; t), t) is well-behaved u(x) in Ω P u(x) in Ω C 4
15 An effective way to find such a mapping is to min x(ξ;t) Ω c + β 2 v 2 ξ dξ which can be approximated by solving a nonlinear moving mesh PDE (Huang, Ren & Wang, Ceniceros & Hou, Du et al. 6): x d t = τ (a i a j ) ξ i( + β 2 ξ v 2 x ξ j) with a i = x ξ i. i,j= 5
16 With the coordinate transform, the elliptic operator with constant coefficients will become one with variable coefficients, but can still be solved using a suitable problem with constant coefficients as preconditioner. Typically a moving mesh strategy can reduce the number of points needed in each direction by a factor of 3-4, leading to significantly savings, particularly in 3-D. We are currently developing an adaptive (Legendre polynomial based) multi-wavelet method which is capable to further reduce the total number of points. 6
17 7
18 8
19 Two-phase Newtonian flow with large density ratios: air bubble in water 9
20 Figure : ρ =.6 and ρ max = 995 and homogenous viscosity µ =.86 at t =.5,.75,,.5,.75, 2. 2
21 2.5 IP PM NIP GUM 2.5 IP PM NIP GUM 2.5 IP PM NIP GUM IP PM NIP GUM 2.5 IP PM NIP GUM 2.5 IP PM NIP GUM
22 Figure 3: ρ =.6, µ =.86 and ρ 2 = 995, µ 2 =.7977 at t =.5,.75,,.5,.75, 2. 22
23 2.5 IP PM NIP GUM 2.5 IP PM NIP GUM 2.5 IP PM NIP GUM IP PM NIP GUM 2.5 IP PM NIP GUM 2.5 IP PM NIP GUM
24 Figure 5: inverted heart shape in experiment by Akers & Belmonte 6: air bubble rising in a polymeric fluid 24
25 Mixture of nematic liquid crystals and a Newtonian fluid Let φ = representing the liquid crystal drop (represented by director n) and φ = representing the Newtonian fluid. Mixing energy density: f mix (φ, φ) = λ 2 φ 2 + λ 4η 2(φ2 ) 2, Bulk energy density: f bulk = K [ 2 n : ( n)t + ( n 2 ) 2 ] 4δ 2, 25
26 Anchoring energy density: f anch = A 2 (n φ)2, Total energy: F (φ, n, φ, n) = Ω f mix + + φ 2 f bulk + f anch where ( + φ)/2 is the volume fraction of the nematic component. 26
27 Governing equations φ t + v φ = γ 2δF (Cahn-Hilliard) or, δφ φ t + v φ = γ δf (Allen-Cahn) + a Lagrange multiplier, δφ n t + v n = γ δf 2 δn, v =, ( ) v ρ t + v v = p + [µ( v + v T ) + τ e] + g(ρ), where γ is the interfacial mobility and γ 2 determines the relax- 27
28 ation time of n, with [ ] δf δφ = λ 2 φ + φ(φ2 ) η 2 δf δn = K [ ( ) + φ 2 n + 2 f bulk A [(n φ)n], + + φ 2 (n 2 ] )n δ 2 A(n φ) φ, τ e = λ( φ φ) K + φ 2 ( n) ( n)t A(n φ)n φ. 28
29 29 Figure 6: A Newtonian bubble rising in a nematic fluid (Zhou, Feng, Yue, Liu & S. 7)
30 Nematic liquid crystal drops in a Newtonian matrix under simple shear Figure 7: Newtonian drop in a Newtonian matrix: Cahn-Hilliard phase equation 3
31 3 Figure 8: Liquid crystal drop in a Newtonian matrix: Top, Cahn-Hilliard; bottom: Allen-Cahn
32 Figure 9: Liquid crystal drop in a Newtonian matrix (Cahn-Hilliard): left figure with larger viscosity 32
33 Concluding remarks The mixture of two incompressible fluids can be described by a flexible and robust energetic phase field model. An alternative phase-field model based on the density is proposed: it is particularly suitable for problems with large density ratios and it admits an energy law. Some efficient numerical approaches for solving the coupled system are proposed: a stablized semi-implicit scheme for the phase equation; 33
34 new projection type schemes which only require solving regular Poisson equation for the pressure; a moving mesh spectral discretization in space. Advantages and challenges: The proposed numerical schemes are easy to implement as they are all based only on elliptic solves with regular Poisson equation for the pressure. It can be easily extended to handle some non-newtonian flows such as liquid crystal flows and visco-elastic flows,... 34
35 It is possible to treat multiphase (> 2) flows by introducing multiple phase functions. Treating all nonlinear terms explicitly lead to a somewhat restrictive time step constraint. We are currently working on ways to efficiently treat some nonlinear terms implicitly to improve the stability. Thank you! 35
c 2010 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 3, No. 3, pp. 59 79 c 00 Society for Industrial and Applied Mathematics A PHASE-FIELD MODEL AND ITS NUMERICAL APPROXIMATION FOR TWO-PHASE INCOMPRESSIBLE FLOWS WITH DIFFERENT DENSITIES
More informationERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION. Xiaofeng Yang. (Communicated by Jie Shen)
DISCRETE AND CONTINUOUS doi:1.3934/dcdsb.29.11.157 DYNAMICAL SYSTEMS SERIES B Volume 11, Number 4, June 29 pp. 157 17 ERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION Xiaofeng Yang
More informationGiordano Tierra Chica
NUMERICAL METHODS FOR SOLVING THE CAHN-HILLIARD EQUATION AND ITS APPLICABILITY TO MIXTURES OF NEMATIC-ISOTROPIC FLOWS WITH ANCHORING EFFECTS Giordano Tierra Chica Department of Mathematics Temple University
More informationA PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM
A PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM Hyun Geun LEE 1, Jeong-Whan CHOI 1 and Junseok KIM 1 1) Department of Mathematics, Korea University, Seoul
More informationCOMPUTING FRACTIONAL LAPLACIANS ON COMPLEX-GEOMETRY DOMAINS: ALGORITHMS AND SIMULATIONS
SIAM J. SCI. COMPUT. Vol. 39, No. 4, pp. A3 A344 c 7 Society for Industrial and Applied Mathematics COMPUTING FRACTIONAL LAPLACIANS ON COMPLEX-GEOMETRY DOMAINS: ALGORITHMS AND SIMULATIONS FANGYING SONG,
More informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
More informationNUMERICAL APPROXIMATIONS OF ALLEN-CAHN AND CAHN-HILLIARD EQUATIONS
UMERICAL APPROXIMATIOS OF ALLE-CAH AD CAH-HILLIARD EQUATIOS JIE SHE AD XIAOFEG YAG Abstract. Stability analyses and error estimates are carried out for a number of commonly used numerical schemes for the
More informationJournal of Computational Physics
Journal of Computational Physics 305 016 539 556 Contents lists available at ScienceDirect Journal of Computational Physics www.elsevier.com/locate/jcp A decoupled energy stable scheme for a hydrodynamic
More informationA New Approach for the Numerical Solution of Diffusion Equations with Variable and Degenerate Mobility
A New Approach for the Numerical Solution of Diffusion Equations with Variable and Degenerate Mobility Hector D. Ceniceros Department of Mathematics, University of California Santa Barbara, CA 93106 Carlos
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 Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
More informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
More informationFully Discrete Second-Order Linear Schemes for Hydrodynamic Phase Field Models of Binary Viscous Fluid Flows with Variable Densities
Utah State University DigitalCommons@USU Mathematics and Statistics Faculty Publications Mathematics and Statistics -30-08 Fully Discrete Second-Order Linear Schemes for Hydrodynamic Phase Field Models
More informationConvergence of the MAC scheme for incompressible flows
Convergence of the MAC scheme for incompressible flows T. Galloue t?, R. Herbin?, J.-C. Latche??, K. Mallem????????? Aix-Marseille Universite I.R.S.N. Cadarache Universite d Annaba Calanque de Cortiou
More information3D simulations of phase separation with a liquid crystal component
Trabalho apresentado no CNMAC, Gramado - RS, 2016. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics D simulations of phase separation with a liquid crystal component
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationThe Shallow Water Equations
If you have not already done so, you are strongly encouraged to read the companion file on the non-divergent barotropic vorticity equation, before proceeding to this shallow water case. We do not repeat
More informationA Fourier Spectral Moving Mesh Method for the Cahn-Hilliard Equation with Elasticity
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 5, No. 2-4, pp. 582-599 Commun. Comput. Phys. February 2009 A Fourier Spectral Moving Mesh Method for the Cahn-Hilliard Equation with Elasticity W. M. Feng
More informationAn asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations
An asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations Institute of Applied Physics and Computational Mathematics, Beijing NUS, Singapore, March 2-6, 2015 (joint
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
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 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 informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
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 informationA Variational Analysis of a Gauged Nonlinear Schrödinger Equation
A Variational Analysis of a Gauged Nonlinear Schrödinger Equation Alessio Pomponio, joint work with David Ruiz Dipartimento di Meccanica, Matematica e Management, Politecnico di Bari Variational and Topological
More informationHIGHLY EFFICIENT AND ACCURATE NUMERICAL SCHEMES FOR THE EPITAXIAL THIN FILM GROWTH MODELS BY USING THE SAV APPROACH
HIGHLY EFFICIENT AND ACCURATE NUMERICAL SCHEMES FOR THE EPITAXIAL THIN FILM GROWTH MODELS BY USING THE SAV APPROACH QING CHENG, JIE SHEN, XIAOFENG YANG Abstract. We develop in this paper highly efficient,
More informationNavier-Stokes equations
1 Navier-Stokes equations Introduction to spectral methods for the CSC Lunchbytes Seminar Series. Incompressible, hydrodynamic turbulence is described completely by the Navier-Stokes equations where t
More informationNonsmooth Schur Newton Methods and Applications
Nonsmooth Schur Newton Methods and Applications C. Gräser, R. Kornhuber, and U. Sack IMA Annual Program Year Workshop Numerical Solutions of Partial Differential Equations: Fast Solution Techniques November
More informationON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS. Chun Liu and Jie Shen
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 7, Number2, April2001 pp. 307 318 ON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS Chun Liu and Jie Shen Department
More informationNonlinear elasticity and gels
Nonlinear elasticity and gels M. Carme Calderer School of Mathematics University of Minnesota New Mexico Analysis Seminar New Mexico State University April 4-6, 2008 1 / 23 Outline Balance laws for gels
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 informationA phase field model for the coupling between Navier-Stokes and e
A phase field model for the coupling between Navier-Stokes and electrokinetic equations Instituto de Matemáticas, CSIC Collaborators: C. Eck, G. Grün, F. Klingbeil (Erlangen Univertsität), O. Vantzos (Bonn)
More informationChapter 7: Natural Convection
7-1 Introduction 7- The Grashof Number 7-3 Natural Convection over Surfaces 7-4 Natural Convection Inside Enclosures 7-5 Similarity Solution 7-6 Integral Method 7-7 Combined Natural and Forced Convection
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 informationThree-dimensional coarsening dynamics of a conserved, nematic liquid crystal-isotropic fluid mixture
Three-dimensional coarsening dynamics of a conserved, nematic liquid crystal-isotropic fluid mixture Rudimar L. Nós a,, Alexandre M. Roma b, Carlos J. García-Cervera c, Hector D. Ceniceros d a Departamento
More informationNumerics for Liquid Crystals with Variable Degree of Orientation
Numerics for Liquid Crystals with Variable Degree of Orientation Ricardo H. Nochetto, Shawn W. Walker 2 and Wujun hang Department of Mathematics, University of Maryland 2 Department of Mathematics, Louisiana
More informationIntroduction to immersed boundary method
Introduction to immersed boundary method Ming-Chih Lai mclai@math.nctu.edu.tw Department of Applied Mathematics Center of Mathematical Modeling and Scientific Computing National Chiao Tung University 1001,
More informationIncompressible MHD simulations
Incompressible MHD simulations Felix Spanier 1 Lehrstuhl für Astronomie Universität Würzburg Simulation methods in astrophysics Felix Spanier (Uni Würzburg) Simulation methods in astrophysics 1 / 20 Outline
More informationPHASE TRANSITIONS IN SOFT MATTER SYSTEMS
OUTLINE: Topic D. PHASE TRANSITIONS IN SOFT MATTER SYSTEMS Definition of a phase Classification of phase transitions Thermodynamics of mixing (gases, polymers, etc.) Mean-field approaches in the spirit
More informationSome results on the nematic liquid crystals theory
Some results on the nematic liquid crystals theory Marius Paicu University of Bordeaux joint work with Arghir Zarnescu Mathflows 2015, Porquerolles September 17, 2015 Complex fluids: Basic laws Incompressibility:
More informationA Momentum Exchange-based Immersed Boundary-Lattice. Boltzmann Method for Fluid Structure Interaction
APCOM & ISCM -4 th December, 03, Singapore A Momentum Exchange-based Immersed Boundary-Lattice Boltzmann Method for Fluid Structure Interaction Jianfei Yang,,3, Zhengdao Wang,,3, and *Yuehong Qian,,3,4
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
More 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 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 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 informationReport Documentation Page
Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,
More informationZ. Guo a,, P. Lin b,, J. Lowengrub a,, S.M. Wise c,
Mass Conservative and Energy Stable Finite Difference Methods for the Quasi-incompressible Navier-Stokes-Cahn-Hilliard system: Primitive Variable and Projection-Type Schemes Z. Guo a,, P. Lin b,, J. Lowengrub
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
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 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 informationIsogeometric Analysis of the Navier Stokes Cahn Hilliard equations with application to incompressible two-phase flows
Isogeometric Analysis of the Navier Stokes Cahn Hilliard equations with application to incompressible two-phase flows Babak S. Hosseini 1, Stefan Turek 1, Matthias Möller 2 1 TU Dortmund University, Institute
More informationIntroduction to Heat and Mass Transfer. Week 9
Introduction to Heat and Mass Transfer Week 9 補充! Multidimensional Effects Transient problems with heat transfer in two or three dimensions can be considered using the solutions obtained for one dimensional
More informationSpline Element Method for Partial Differential Equations
for Partial Differential Equations Department of Mathematical Sciences Northern Illinois University 2009 Multivariate Splines Summer School, Summer 2009 Outline 1 Why multivariate splines for PDEs? Motivation
More informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationNumerical Approximations for Allen-Cahn Type Phase Field Model of Two-Phase Incompressible Fluids with Moving Contact Lines
Commun. Comput. Phys. doi: 1.48/cicp.OA-16-8 Vol. 1, No. 3, pp. 867-889 March 17 Numerical Approximations for Allen-Cahn Type Phase Field Model of Two-Phase Incompressible Fluids with Moving Contact Lines
More informationFinite difference MAC scheme for compressible Navier-Stokes equations
Finite difference MAC scheme for compressible Navier-Stokes equations Bangwei She with R. Hošek, H. Mizerova Workshop on Mathematical Fluid Dynamics, Bad Boll May. 07 May. 11, 2018 Compressible barotropic
More informationNon-Newtonian Fluids and Finite Elements
Non-Newtonian Fluids and Finite Elements Janice Giudice Oxford University Computing Laboratory Keble College Talk Outline Motivating Industrial Process Multiple Extrusion of Pastes Governing Equations
More informationPriority Programme The Combination of POD Model Reduction with Adaptive Finite Element Methods in the Context of Phase Field Models
Priority Programme 1962 The Combination of POD Model Reduction with Adaptive Finite Element Methods in the Context of Phase Field Models Carmen Gräßle, Michael Hinze Non-smooth and Complementarity-based
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More informationConstrained Minimization and Multigrid
Constrained Minimization and Multigrid C. Gräser (FU Berlin), R. Kornhuber (FU Berlin), and O. Sander (FU Berlin) Workshop on PDE Constrained Optimization Hamburg, March 27-29, 2008 Matheon Outline Successive
More informationNonlinear Solitary Wave Basis Functions for Long Time Solution of Wave Equation
Nonlinear Solitary Wave Basis Functions for Long Time Solution of Wave Equation John Steinhoff Subhashini Chitta Nonlinear Waves --Theory and Applications Beijing - June, 2008 Objective Simulate short
More informationTHE CAHN-HILLIARD EQUATION WITH A LOGARITHMIC POTENTIAL AND DYNAMIC BOUNDARY CONDITIONS
THE CAHN-HILLIARD EQUATION WITH A LOGARITHMIC POTENTIAL AND DYNAMIC BOUNDARY CONDITIONS Alain Miranville Université de Poitiers, France Collaborators : L. Cherfils, G. Gilardi, G.R. Goldstein, G. Schimperna,
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 informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
More informationElmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties. SIF file : phasechange solid-solid
Elmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties 3 6 1. Tb=1750 [K] 2 & 5. q=-10000 [W/m²] 0,1 1 Ω1 4 Ω2 7 3 & 6. α=15 [W/(m²K)] Text=300 [K] 4.
More informationHYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS
1 / 36 HYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS Jesús Garicano Mena, E. Valero Sánchez, G. Rubio Calzado, E. Ferrer Vaccarezza Universidad
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
More informationFinite Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
More informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
More informationNumerical Optimization of Partial Differential Equations
Numerical Optimization of Partial Differential Equations Bartosz Protas Department of Mathematics & Statistics McMaster University, Hamilton, Ontario, Canada URL: http://www.math.mcmaster.ca/bprotas Rencontres
More informationNonlinear Analysis: Modelling and Control, 2008, Vol. 13, No. 4,
Nonlinear Analysis: Modelling and Control, 2008, Vol. 13, No. 4, 513 524 Effects of Temperature Dependent Thermal Conductivity on Magnetohydrodynamic (MHD) Free Convection Flow along a Vertical Flat Plate
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 informationLECTURE 4: GEOMETRIC PROBLEMS
LECTURE 4: GEOMETRIC PROBLEMS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and Related
More informationModelling of interfaces and free boundaries
University of Regensburg Regensburg, March 2009 Outline 1 Introduction 2 Obstacle problems 3 Stefan problem 4 Shape optimization Introduction What is a free boundary problem? Solve a partial differential
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 informationRicci Dark Energy Chao-Jun Feng SUCA, SHNU
Ricci Dark Energy 21 2011.3.18-19 Chao-Jun Feng SUCA, SHNU Outline n Introduction to Ricci Dark Energy (RDE) n Stability and the constraint on the parameter n Age problem alleviated in viscous version
More informationAllen Cahn Equation in Two Spatial Dimension
Allen Cahn Equation in Two Spatial Dimension Yoichiro Mori April 25, 216 Consider the Allen Cahn equation in two spatial dimension: ɛ u = ɛ2 u + fu) 1) where ɛ > is a small parameter and fu) is of cubic
More informationA High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations
A High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations Khosro Shahbazi 1, Paul F. Fischer 2 and C. Ross Ethier 1 1 University of Toronto and 2 Argonne National
More informationFinite Difference Solution of the Heat Equation
Finite Difference Solution of the Heat Equation Adam Powell 22.091 March 13 15, 2002 In example 4.3 (p. 10) of his lecture notes for March 11, Rodolfo Rosales gives the constant-density heat equation as:
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationChapter 4. COSMOLOGICAL PERTURBATION THEORY
Chapter 4. COSMOLOGICAL PERTURBATION THEORY 4.1. NEWTONIAN PERTURBATION THEORY Newtonian gravity is an adequate description on small scales (< H 1 ) and for non-relativistic matter (CDM + baryons after
More informationThe Spectral-Element Method: Introduction
The Spectral-Element Method: Introduction Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Computational Seismology 1 / 59 Outline 1 Introduction 2 Lagrange
More informationEnergy Preserving Numerical Integration Methods
Energy Preserving Numerical Integration Methods Department of Mathematics and Statistics, La Trobe University Supported by Australian Research Council Professorial Fellowship Collaborators: Celledoni (Trondheim)
More informationFEniCS Course. Lecture 6: Incompressible Navier Stokes. Contributors Anders Logg André Massing
FEniCS Course Lecture 6: Incompressible Navier Stokes Contributors Anders Logg André Massing 1 / 11 The incompressible Navier Stokes equations u + u u ν u + p = f in Ω (0, T ] u = 0 in Ω (0, T ] u = g
More informationTitle of your thesis here
Title of your thesis here Your name here A Thesis presented for the degree of Doctor of Philosophy Your Research Group Here Department of Mathematical Sciences University of Durham England Month and Year
More information(1:1) 1. The gauge formulation of the Navier-Stokes equation We start with the incompressible Navier-Stokes equation 8 >< >: u t +(u r)u + rp = 1 Re 4
Gauge Finite Element Method for Incompressible Flows Weinan E 1 Courant Institute of Mathematical Sciences New York, NY 10012 Jian-Guo Liu 2 Temple University Philadelphia, PA 19122 Abstract: We present
More informationLinear Theory of Evolution to an Unstable State
Chapter 2 Linear Theory of Evolution to an Unstable State c 2012 by William Klein, Harvey Gould, and Jan Tobochnik 1 October 2012 2.1 Introduction The simple theory of nucleation that we introduced in
More informationComputers and Mathematics with Applications
Computers and Mathematics with Applications 6 () 77 745 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa Multiphase
More informationNUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART II: ERROR ANALYSIS AND CONVERGENCE OF THE INTERFACE
NUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART II: ERROR ANALYSIS AND CONVERGENCE OF THE INTERFACE XIAOBING FENG AND ANDREAS PROHL Abstract. In this
More informationAn Efficient, Energy Stable Scheme for the Cahn-Hilliard-Brinkman System
An Efficient, Energy Stable Scheme for the Cahn-Hilliard-Brinkman System C. Collins Mathematics Department The University of Tennessee; Knoxville, TN 37996-064 craig.collins@math.utk.edu Jie Shen Mathematics
More informationLecture 1: Water waves and Hamiltonian partial differential equations
Lecture : Water waves and Hamiltonian partial differential equations Walter Craig Department of Mathematics & Statistics Waves in Flows Prague Summer School 208 Czech Academy of Science August 30 208 Outline...
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 informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Partial differential equations arise in a number of physical problems, such as fluid flow, heat transfer, solid mechanics and biological processes. These
More informationA formulation for fast computations of rigid particulate flows
Center for Turbulence Research Annual Research Briefs 2001 185 A formulation for fast computations of rigid particulate flows By N. A. Patankar 1. Introduction A formulation is presented for the direct
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationTurbulent drag reduction by streamwise traveling waves
51st IEEE Conference on Decision and Control December 10-13, 2012. Maui, Hawaii, USA Turbulent drag reduction by streamwise traveling waves Armin Zare, Binh K. Lieu, and Mihailo R. Jovanović Abstract For
More informationEMAFF: MAgnetic Equations with FreeFem++
EMAFF: MAgnetic Equations with FreeFem++ The Grad-Shafranov equation & The Current Hole Erwan DERIAZ Bruno DESPRÉS Gloria FACCANONI š Kirill Pichon GOSTAF Lise-Marie IMBERT-GÉRARD Georges SADAKA Remy
More information