Variational Assimilation of Discrete Navier-Stokes Equations
|
|
- Lynne McDowell
- 5 years ago
- Views:
Transcription
1 Variational Assimilation of Discrete Navier-Stokes Equations Souleymane.Kadri-Harouna FLUMINANCE, INRIA Rennes-Bretagne Atlantique Campus universitaire de Beaulieu, Rennes, France
2 Outline Discretization of Navier-Stokes Equations Temporal discretization Spatial discretization
3 Outline Discretization of Navier-Stokes Equations Temporal discretization Spatial discretization Variational Assimilation Principle Discrete adjoint Method
4 Outline Discretization of Navier-Stokes Equations Temporal discretization Spatial discretization Variational Assimilation Principle Discrete adjoint Method Preliminary resultats
5 Outline Discretization of Navier-Stokes Equations Temporal discretization Spatial discretization Variational Assimilation Principle Discrete adjoint Method Preliminary resultats Conclusion and Outlook
6 Incompressible Navier-Stokes Equations Cauchy problem for Navier-Stokes: t v ν v + (v )v + p = f, x Ω, t [0, T ], (NS) v = 0, x Ω, t [0, T ], v(0, x) = v 0 (x), x Ω. Unknowns : velocity v(t, x) and pressure p(t, x) Projecting the system (NS) onto H div (Ω) 1 yields: t v = ν v + P[ (v )v + f] (NSP) with P orthogonal projector from (L 2 (Ω)) d to H div (Ω). The pressure p is recovered through the Helmholtz decomposition: p = (v )v + f P[ (v )v + f] 1 (L 2 (Ω)) d divergence-free function space
7 Helmholtz-Decomposition The projector P is explicite in Fourier domain: u (L 2 (Ω)) d, û(ξ) = ξ (1 ξt ξ 2 û(ξ) + ξ ) ξt ξ 2 û(ξ) Thus P(u)(ξ) = (1 ξ ξt ξ 2 ) û(ξ) For space localization and adaptativity: Periodic Anisotropic divergence-free wavelets [Deriaz,Perrier 08] For physical boundary conditions: Anisotropic divergence-free wavelets [Kadri-Harouna,Perrier 10]
8 Outline Discretization of Navier-Stokes Equations Temporal discretization Spatial discretization Variational Assimilation Principle Discrete adjoint Method Preliminary resultats Conclusion and Outlook
9 Temporal discretization Heat kernel integration problem: t v ν v = f, with f = P[ (v )v + f]. Implicite finite difference approximation v(x, nδt) v n : v n+1 v n δt ν 2 (vn+1 + v n ) = f n, Crank-Nicholson O(δt 2 ) Heat kernel factorization (ADI method): ) ) ) (1 α 2 x 2 α 2 y 2 = (1 α (1 2 x 2 α 2 y 2 + O(α 2 )
10 Outline Discretization of Navier-Stokes Equations Temporal discretization Spatial discretization Variational Assimilation Principle Discrete adjoint Method Preliminary resultats Conclusion and Outlook
11 Spatial discretization Semi implicite treatment for the non-linear term: (v n+1/2 )v n+1/2 = 3 2 (vn )v n 1 2 (vn 1 )v n 1 CFL condition: δt C (δx/v max ) 4/3. Scale separation: v(t, x) = j,k d j,k div (t) Ψdiv j,k (x) ODE system on the coefficients [d div j,k (t)]. Galerkin method in space with V j = (V 1 j V 0 j ) (V 0 j V 1 j ). At each time step we need to compute the projector P.
12 Outline Discretization of Navier-Stokes Equations Temporal discretization Spatial discretization Variational Assimilation Principle Discrete adjoint Method Preliminary resultats Conclusion and Outlook
13 Principle of Variational Assimilation Measurements (observations) denoted v k ob, k = 1,, N Discrete dynamical model equation: with L 1/2 v n+1 L 1/2 v n Bn v n 1 2 Bn 1 v n 1 = 0, L 1/2 := 1 δt 2, L 1/2 := 1 + δt 2, Bn v n := P(v n )v n Objective: find the most probable state defined both by the measurements and dynamical equations. Cost function minimization: J(v 0 ) = 1 2 N H v k v k ob 2 δt + α 2 v 0 2, k=1
14 Differentiation operators Let f : E R be a vector (or scalar) function Directional derivative: if the following limit exists d f (v) := lim h 0 f (v + hd) f (v) h Fréchet derivative: if there exist f (v) E such that f (v + u) = f (v) + f (v), u + o( u ) If the gradient of f exists, then: d f (v) = f (v), d
15 Cost function differentiation J(v 0 + hu) J(v 0 ) = N H v k (v 0 + hu) v k ob 2 δt + α 2 v 0 + hu 2 k=1 N H v k v k ob 2 δt α 2 v 0 2 k=1 Rewritten the terms, we get: H v k (v 0 + hu) v k ob 2 H v k v k ob 2 = + H v k (v 0 + hu) + H v k 2v k ob, H v k (v 0 + hu) H v k and v 0 + hu 2 v 0 2 = 2v 0 + hu, hu Thus: N uj(v 0 ) = H v k v k ob, H v k uv k δt + α v 0, u k=1
16 Outline Discretization of Navier-Stokes Equations Temporal discretization Spatial discretization Variational Assimilation Principle Discrete adjoint Method Preliminary resultats Conclusion and Outlook
17 Problem Let us consider the one-dimensional ODE t y = F (t).y, with F (t) a linear operator The continuous adjoint model is t λ = F (t) T.λ Discretizing with an explicit Euler scheme, we get: y n+1 y n = δtf n.y n y n+1 = (1 + δtf n ).y n For which we get Otherwise: y n = (1 + δtf n ) T.y n+1 λ n λ n+1 = δtf T n+1.λ n+1 λ n = (1 + δtf n+1 ) T.λ n+1 (1 + δtf n+1 ) T (1 + δtf n ) T
18 Linear tangent One disturbs the initial condition: ṽ 0 = v 0 + hu. Then, we get: L 1/2 ṽ n+1 L 1/2 ṽ n B n ṽ n 1 2 B n 1 ṽ n 1 = 0. Taking the difference with the non disturbed equation, we have: L 1/2 (ṽ n+1 v n+1 ) L 1/2 (ṽ n v n ) = 3 2 ( B n ṽ n B n v n ) ( B n 1 ṽ n 1 B n 1 v n 1 ) Multiplying with 1/h and taking the limit as h 0, we get: L 1/2 u v n+1 L 1/2 u v n = 3 2 Bn v n u v n Bn 1 v n 1 u v n 1
19 Adjoint variable Taking the inner product of the linear tangent with λ n+1 yields u v n+1, L 1/2 λ n+1 u v n, L 1/2 λ n+1 = u v n, 3 2 B n λ n+1 + u v n 1, 1 2 B n 1 λ n+1 Thus, making identification, the adjoint model is defined as: L 1/2 λ N = F N, L 1/2 λ N 1 L 1/2 λ N B N 1 λ N = F N 1 with F n = H v n (H v n v n ob), 1 n N For 1 n N 2, L 1/2 λ n L 1/2 λ n B n λ n B n λ n+2 = F n J(v 0 )= αv 0 + L 1/2 λ B 0 λ B 0 λ 2
20 Model error 6.8 x L 2 Norm Time Figure: L 2 -norm error
21 Error on real experience Two types of observation H v k v k ob := v k I v k ob, v k ob Optical-Flow H v k v k ob := I k 1 (x + v k ) I k 0 (x)
22 Pseudo observations error (a) True vorticity (b) Estimated vorticity Figure: Optical-flow observation: RMSE=
23 Pseudo observations error Figure: RSE on the vorticity =
24 DFD observations error (a) True vorticity (b) Estimated vorticity Figure: DFD observation: RMSE=0.0696, j = 7.
25 DFD observations error Figure: RSE on the vorticity
26 DFD observations and diffusion (a) True vorticity (b) Estimated vorticity Figure: DFD observation: RMSE=0.0523,ν 2.73E 5, Re E 6.
27 Conclusion and Outlook Navier-Stokes discretization Discrete adjoint models
28 Conclusion and Outlook Navier-Stokes discretization Discrete adjoint models Models with low complexity
29 Conclusion and Outlook Navier-Stokes discretization Discrete adjoint models Models with low complexity Wavelet adaptativity in the simulation Use methods on a dynamic geophysical flow models
VALÉRIE PERRIER AND. Key words. Divergence-free wavelets, Navier-Stokes simulation, physical boundary conditions
DIVERGENCE-FREE WAVELET PROJECTION METHOD FOR INCOMPRESSIBLE VISCOUS FLOW SOULEYMANE KADRI HAROUNA AND VALÉRIE PERRIER Abstract. We present a new wavelet numerical scheme for the discretization of Navier-Stokes
More informationDirect Numerical Simulation of Turbulence using divergence-free Wavelets
Submitted to: SIAM Multiscale Modeling and Simulation Direct Numerical Simulation of Turbulence using divergence-free Wavelets Erwan Deriaz Valérie Perrier July 10, 2008 Abstract We present a numerical
More informationVALÉRIE PERRIER AND. Key words. Divergence-free wavelets, Navier-Stokes simulation, physical boundary conditions
DIVERGENCE-FREE WAVELET PROJECTION METHOD FOR INCOMPRESSIBLE VISCOUS FLOW SOULEYMANE KADRI HAROUNA AND VALÉRIE PERRIER Abstract. We present a new wavelet numerical scheme for the discretization of Navier-Stokes
More informationDivergence-free wavelet projection method for incompressible viscous flow on the square
Divergence-free wavelet proection method for incompressible viscous flow on the square Souleymane Kadri Harouna, Valérie Perrier To cite this version: Souleymane Kadri Harouna, Valérie Perrier. Divergence-free
More informationREMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID
REMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID DRAGOŞ IFTIMIE AND JAMES P. KELLIHER Abstract. In [Math. Ann. 336 (2006), 449-489] the authors consider the two dimensional
More informationEvolution equations with spectral methods: the case of the wave equation
Evolution equations with spectral methods: the case of the wave equation Jerome.Novak@obspm.fr Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris, France in collaboration with
More informationUne méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé
Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé CMCS/IACS Ecole Polytechnique Federale de Lausanne Erik.Burman@epfl.ch Méthodes Numériques
More informationVariational assimilation of meteorological images for motion estimation
Variational assimilation of meteorological images for motion estimation T. Corpetti work with P. Heas, N. Papadakis and E. Mémin IRISA - COSTEL Brest, Journée Statistique Spatio-temporelle Feb., 14 th
More informationStability of explicit numerical schemes for convection-dominated problems
Stability of explicit numerical schemes for convection-dominated problems Erwan Deriaz and Dmitry Kolomenskiy July 20, 2010 Abstract Operating a refined von Neumann stability analysis of the transport
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 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 informationSECOND-ORDER FULLY DISCRETIZED PROJECTION METHOD FOR INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
Tenth MSU Conference on Differential Equations and Computational Simulations. Electronic Journal of Differential Equations, Conference 3 (06), pp. 9 0. ISSN: 07-669. URL: http://ejde.math.txstate.edu or
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 informationDivergence-free Wavelets for Navier-Stokes
Divergence-free Wavelets for Navier-Stokes Erwan Deriaz, Valérie Perrier To cite this version: Erwan Deriaz, Valérie Perrier. Divergence-free Wavelets for Navier-Stokes. 17 - M. novembre 4. 5.
More informationSome recent results for two extended Navier-Stokes systems
Some recent results for two extended Navier-Stokes systems Jim Kelliher UC Riverside Joint work with Mihaela Ignatova (UCR), Gautam Iyer (Carnegie Mellon), Bob Pego (Carnegie Mellon), and Arghir Dani Zarnescu
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 informationVector penalty-projection methods for outflow boundary conditions with optimal second-order accuracy
Vector penalty-projection methods for outflow boundary conditions with optimal second-order accuracy Philippe Angot, Rima Cheaytou To cite this version: Philippe Angot, Rima Cheaytou. Vector penalty-projection
More informationChapter 2 Finite Element Spaces for Linear Saddle Point Problems
Chapter 2 Finite Element Spaces for Linear Saddle Point Problems Remark 2.1. Motivation. This chapter deals with the first difficulty inherent to the incompressible Navier Stokes equations, see Remark
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 informationMidterm Solution
18303 Midterm Solution Problem 1: SLP with mixed boundary conditions Consider the following regular) Sturm-Liouville eigenvalue problem consisting in finding scalars λ and functions v : [0, b] R b > 0),
More informationSpectral Processing. Misha Kazhdan
Spectral Processing Misha Kazhdan [Taubin, 1995] A Signal Processing Approach to Fair Surface Design [Desbrun, et al., 1999] Implicit Fairing of Arbitrary Meshes [Vallet and Levy, 2008] Spectral Geometry
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 informationInequalities of Babuška-Aziz and Friedrichs-Velte for differential forms
Inequalities of Babuška-Aziz and Friedrichs-Velte for differential forms Martin Costabel Abstract. For sufficiently smooth bounded plane domains, the equivalence between the inequalities of Babuška Aziz
More informationFirst, Second, and Third Order Finite-Volume Schemes for Diffusion
First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
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 informationAn Example of Self Acceleration for Incompressible Flows. Key Words: incompressible flows, Euler equation, Navier-Stokes equation, self-acceleration
Bol. Soc. Paran. Mat. (3s.) v. 31 2 (213): 149 159. c SPM ISSN-2175-1188 on line ISSN-378712 in press SPM: www.spm.uem.br/bspm doi:1.5269/bspm.v31i2.1738 An Example of Self Acceleration for Incompressible
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 informationDiffusion Equation and Mean Free Path
Diffusion Equation and Mean Free Path Speaker: Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University (SUNY) Content General Introduction Analytic
More informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationarxiv: v2 [math.na] 16 Aug 2012
CALCULUS ON SURFACES WITH GENERAL CLOSEST POINT FUNCTIONS THOMAS MÄRZ, COLIN B. MACDONALD arxiv:1202.3001v2 [math.na] 16 Aug 2012 Abstract. The Closest Point Method for solving partial differential equations
More informationOn the p-laplacian and p-fluids
LMU Munich, Germany Lars Diening On the p-laplacian and p-fluids Lars Diening On the p-laplacian and p-fluids 1/50 p-laplacian Part I p-laplace and basic properties Lars Diening On the p-laplacian and
More informationON SINGULAR PERTURBATION OF THE STOKES PROBLEM
NUMERICAL ANALYSIS AND MATHEMATICAL MODELLING BANACH CENTER PUBLICATIONS, VOLUME 9 INSTITUTE OF MATHEMATICS POLISH ACADEMY OF SCIENCES WARSZAWA 994 ON SINGULAR PERTURBATION OF THE STOKES PROBLEM G. M.
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 informationModeling, Simulating and Rendering Fluids. Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan
Modeling, Simulating and Rendering Fluids Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan Applications Mostly Hollywood Shrek Antz Terminator 3 Many others Games Engineering Animating Fluids is
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 informationSuboptimal Open-loop Control Using POD. Stefan Volkwein
Institute for Mathematics and Scientific Computing University of Graz, Austria PhD program in Mathematics for Technology Catania, May 22, 2007 Motivation Optimal control of evolution problems: min J(y,
More informationOpenFOAM selected solver
OpenFOAM selected solver Roberto Pieri - SCS Italy 16-18 June 2014 Introduction to Navier-Stokes equations and RANS Turbulence modelling Numeric discretization Navier-Stokes equations Convective term {}}{
More informationNAVIER-STOKES EQUATIONS IN THIN 3D DOMAINS WITH NAVIER BOUNDARY CONDITIONS
NAVIER-STOKES EQUATIONS IN THIN 3D DOMAINS WITH NAVIER BOUNDARY CONDITIONS DRAGOŞ IFTIMIE, GENEVIÈVE RAUGEL, AND GEORGE R. SELL Abstract. We consider the Navier-Stokes equations on a thin domain of the
More informationMETHODS FOR SOLVING MATHEMATICAL PHYSICS PROBLEMS
METHODS FOR SOLVING MATHEMATICAL PHYSICS PROBLEMS V.I. Agoshkov, P.B. Dubovski, V.P. Shutyaev CAMBRIDGE INTERNATIONAL SCIENCE PUBLISHING Contents PREFACE 1. MAIN PROBLEMS OF MATHEMATICAL PHYSICS 1 Main
More informationProper Orthogonal Decomposition (POD) for Nonlinear Dynamical Systems. Stefan Volkwein
Proper Orthogonal Decomposition (POD) for Nonlinear Dynamical Systems Institute for Mathematics and Scientific Computing, Austria DISC Summerschool 5 Outline of the talk POD and singular value decomposition
More informationAccelerating incompressible fluid flow simulations on hybrid CPU/GPU systems
Accelerating incompressible fluid flow simulations on hybrid CPU/GPU systems Yushan Wang 1, Marc Baboulin 1,2, Karl Rupp 3,4, Yann Fraigneau 1,5, Olivier Le Maître 1,5 1 Université Paris-Sud, France 2
More informationTOPICS IN MATHEMATICAL AND COMPUTATIONAL FLUID DYNAMICS
TOPICS IN MATHEMATICAL AND COMPUTATIONAL FLUID DYNAMICS Lecture Notes by Prof. Dr. Siddhartha Mishra and Dr. Franziska Weber October 4, 216 1 Contents 1 Fundamental Equations of Fluid Dynamics 4 1.1 Eulerian
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 informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationNumerical explorations of a forward-backward diffusion equation
Numerical explorations of a forward-backward diffusion equation Newton Institute KIT Programme Pauline Lafitte 1 C. Mascia 2 1 SIMPAF - INRIA & U. Lille 1, France 2 Univ. La Sapienza, Roma, Italy September
More informationLeast-Squares Finite Element Methods
Pavel В. Bochev Max D. Gunzburger Least-Squares Finite Element Methods Spri ringer Contents Part I Survey of Variational Principles and Associated Finite Element Methods 1 Classical Variational Methods
More informationDecay profiles of a linear artificial viscosity system
Decay profiles of a linear artificial viscosity system Gene Wayne, Ryan Goh and Roland Welter Boston University rwelter@bu.edu July 2, 2018 This research was supported by funding from the NSF. Roland Welter
More informationImplicit-explicit exponential integrators
Implicit-explicit exponential integrators Bawfeh Kingsley Kometa joint work with Elena Celledoni MaGIC 2011 Finse, March 1-4 1 Introduction Motivation 2 semi-lagrangian Runge-Kutta exponential integrators
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 informationSome Semi-Markov Processes
Some Semi-Markov Processes and the Navier-Stokes equations NW Probability Seminar, October 22, 2006 Mina Ossiander Department of Mathematics Oregon State University 1 Abstract Semi-Markov processes were
More informationVANISHING VISCOSITY IN THE PLANE FOR NONDECAYING VELOCITY AND VORTICITY
VANISHING VISCOSITY IN THE PLANE FOR NONDECAYING VELOCITY AND VORTICITY ELAINE COZZI Abstract. Assuming that initial velocity and initial vorticity are bounded in the plane, we show that on a sufficiently
More informationExistence and Continuation for Euler Equations
Existence and Continuation for Euler Equations David Driver Zhuo Min Harold Lim 22 March, 214 Contents 1 Introduction 1 1.1 Physical Background................................... 1 1.2 Notation and Conventions................................
More informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
More informationFFTs in Graphics and Vision. The Laplace Operator
FFTs in Graphics and Vision The Laplace Operator 1 Outline Math Stuff Symmetric/Hermitian Matrices Lagrange Multipliers Diagonalizing Symmetric Matrices The Laplacian Operator 2 Linear Operators Definition:
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 informationThe Navier-Stokes Equations with Time Delay. Werner Varnhorn. Faculty of Mathematics University of Kassel, Germany
The Navier-Stokes Equations with Time Delay Werner Varnhorn Faculty of Mathematics University of Kassel, Germany AMS: 35 (A 35, D 5, K 55, Q 1), 65 M 1, 76 D 5 Abstract In the present paper we use a time
More informationZdzislaw Brzeźniak. Department of Mathematics University of York. joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York)
Navier-Stokes equations with constrained L 2 energy of the solution Zdzislaw Brzeźniak Department of Mathematics University of York joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York) LMS
More informationExperience with DNS of particulate flow using a variant of the immersed boundary method
Experience with DNS of particulate flow using a variant of the immersed boundary method Markus Uhlmann Numerical Simulation and Modeling Unit CIEMAT Madrid, Spain ECCOMAS CFD 2006 Motivation wide range
More informationKonstantinos Chrysafinos 1 and L. Steven Hou Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANALYSIS AND APPROXIMATIONS OF THE EVOLUTIONARY STOKES EQUATIONS WITH INHOMOGENEOUS
More informationPartial Differential Equations II
Partial Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Partial Differential Equations II 1 / 28 Almost Done! Homework
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 informationPOD/DEIM 4DVAR Data Assimilation of the Shallow Water Equation Model
nonlinear 4DVAR 4DVAR Data Assimilation of the Shallow Water Equation Model R. Ştefănescu and Ionel M. Department of Scientific Computing Florida State University Tallahassee, Florida May 23, 2013 (Florida
More informationREDUCED-ORDER MODELING OF HIDDEN DYNAMICS
REDUCED-ORDER MODELING OF HIDDEN DYNAMICS Patrick Héas and Cédric Herzet INRIA Centre Rennes - Bretagne Atlantique, Campus universitaire de Beaulieu, 35000 Rennes, France ABSTRACT The objective of this
More informationDISCRETE HODGE HELMHOLTZ
Monografías Matemáticas García de Galdeano 39, 1 10 (2014) DISCRETE HODGE HELMHOLTZ DECOMPOSITION Etienne Ahusborde, Mejdi Azaïez, Jean-Paul Caltagirone, Mark Gerritsma and Antoine Lemoine Abstract. This
More informationHelmholtz-Hodge Decomposition on [0, 1] d by Divergence-free and Curl-free Wavelets
Helmholtz-Hodge Decomposition on [, 1] d by Divergence-free and Curl-free Wavelets Souleymane Kadri Harouna and Valérie Perrier * Springer-Verlag, Computer Science Editorial, Tiergartenstr. 17, 69121 Heidelberg,
More information8 A pseudo-spectral solution to the Stokes Problem
8 A pseudo-spectral solution to the Stokes Problem 8.1 The Method 8.1.1 Generalities We are interested in setting up a pseudo-spectral method for the following Stokes Problem u σu p = f in Ω u = 0 in Ω,
More informationRemarks on the blow-up criterion of the 3D Euler equations
Remarks on the blow-up criterion of the 3D Euler equations Dongho Chae Department of Mathematics Sungkyunkwan University Suwon 44-746, Korea e-mail : chae@skku.edu Abstract In this note we prove that the
More informationGlobal well-posedness and decay for the viscous surface wave problem without surface tension
Global well-posedness and decay for the viscous surface wave problem without surface tension Ian Tice (joint work with Yan Guo) Université Paris-Est Créteil Laboratoire d Analyse et de Mathématiques Appliquées
More informationMULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
More informationLubrication and roughness
Lubrication and roughness Laurent Chupin 1 & Sébastien Martin 2 1 - Institut Camille Jordan - Lyon 2 - Laboratoire de Mathématiques - Orsay GdR CHANT - August 2010 Laurent Chupin (ICJ Lyon) Lubrication
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 informationZdzislaw Brzeźniak. Department of Mathematics University of York. joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York)
Navier-Stokes equations with constrained L 2 energy of the solution Zdzislaw Brzeźniak Department of Mathematics University of York joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York) Stochastic
More informationImplementation of 3D Incompressible N-S Equations. Mikhail Sekachev
Implementation of 3D Incompressible N-S Equations Mikhail Sekachev Navier-Stokes Equations The motion of a viscous incompressible fluid is governed by the Navier-Stokes equations u + t u = ( u ) 0 Quick
More informationCONSTRUCTION AND ANALYSIS OF WEIGHTED SEQUENTIAL SPLITTING FDTD METHODS FOR THE 3D MAXWELL S EQUATIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 5 Number 6 Pages 747 784 c 08 Institute for Scientific Computing and Information CONSTRUCTION AND ANALYSIS OF WEIGHTED SEQUENTIAL SPLITTING
More informationGlobal Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations
Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations C. David Levermore Department of Mathematics and Institute for Physical Science and Technology
More informationA DG-based incompressible Navier-Stokes solver
A DG-based incompressible Navier-Stokes solver Higer-order DG metods and finite element software for modern arcitectures - Bat 2016 Marian Piatkowski 1 Peter Bastian 1 1 Interdisciplinary Center for Scientific
More informationON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS
ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS Abdelhafid Younsi To cite this version: Abdelhafid Younsi. ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS. 4 pages. 212. HAL Id:
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 informationProper Orthogonal Decomposition. POD for PDE Constrained Optimization. Stefan Volkwein
Proper Orthogonal Decomposition for PDE Constrained Optimization Institute of Mathematics and Statistics, University of Constance Joined work with F. Diwoky, M. Hinze, D. Hömberg, M. Kahlbacher, E. Kammann,
More informationExistence of minimizers for the pure displacement problem in nonlinear elasticity
Existence of minimizers for the pure displacement problem in nonlinear elasticity Cristinel Mardare Université Pierre et Marie Curie - Paris 6, Laboratoire Jacques-Louis Lions, Paris, F-75005 France Abstract
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 informationPost-processing of solutions of incompressible Navier Stokes equations on rotating spheres
ANZIAM J. 50 (CTAC2008) pp.c90 C106, 2008 C90 Post-processing of solutions of incompressible Navier Stokes equations on rotating spheres M. Ganesh 1 Q. T. Le Gia 2 (Received 14 August 2008; revised 03
More informationIntroduction to Signal Spaces
Introduction to Signal Spaces Selin Aviyente Department of Electrical and Computer Engineering Michigan State University January 12, 2010 Motivation Outline 1 Motivation 2 Vector Space 3 Inner Product
More informationGEOMETRICAL TOOLS FOR PDES ON A SURFACE WITH ROTATING SHALLOW WATER EQUATIONS ON A SPHERE
GEOETRICAL TOOLS FOR PDES ON A SURFACE WITH ROTATING SHALLOW WATER EQUATIONS ON A SPHERE BIN CHENG Abstract. This is an excerpt from my paper with A. ahalov [1]. PDE theories with Riemannian geometry are
More informationApproximate tensor-product preconditioners for very high order discontinuous Galerkin methods
Approximate tensor-product preconditioners for very high order discontinuous Galerkin methods Will Pazner 1 and Per-Olof Persson 2 1 Division of Applied Mathematics, Brown University, Providence, RI, 02912
More informationContinuum Mechanics Lecture 5 Ideal fluids
Continuum Mechanics Lecture 5 Ideal fluids Prof. http://www.itp.uzh.ch/~teyssier Outline - Helmholtz decomposition - Divergence and curl theorem - Kelvin s circulation theorem - The vorticity equation
More informationADAPTATION OF STANDARD OPTIC FLOW METHODS TO FLUID MOTION
9TH INTERNATIONAL SYMPOSIUM ON FLOW VISUALISATION, 2000 ADAPTATION OF STANDARD OPTIC FLOW METHODS TO FLUID MOTION T. Corpetti 1, E. Mémin 2, P. Pérez 3 Keywords: Fluid motion, optical flow, continuity
More informationA Riemannian Framework for Denoising Diffusion Tensor Images
A Riemannian Framework for Denoising Diffusion Tensor Images Manasi Datar No Institute Given Abstract. Diffusion Tensor Imaging (DTI) is a relatively new imaging modality that has been extensively used
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 2198-5855 On the divergence constraint in mixed finite element methods for incompressible
More informationLagrangian Dynamics & Mixing
V Lagrangian Dynamics & Mixing See T&L, Ch. 7 The study of the Lagrangian dynamics of turbulence is, at once, very old and very new. Some of the earliest work on fluid turbulence in the 1920 s, 30 s and
More informationSeminar on Vector Field Analysis on Surfaces
Seminar on Vector Field Analysis on Surfaces 236629 1 Last time Intro Cool stuff VFs on 2D Euclidean Domains Arrows on the plane Div, curl and all that Helmholtz decomposition 2 today Arrows on surfaces
More informationarxiv: v1 [physics.comp-ph] 30 Sep 2015
On the quasi-unconditional stability of BDF-ADI solvers for the compressible Navier-Stokes equations arxiv:1509.09213v1 [physics.comp-ph] 30 Sep 2015 Oscar P. Bruno and Max Cubillos Abstract The companion
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More informationRegularity and Decay Estimates of the Navier-Stokes Equations
Regularity and Decay Estimates of the Navier-Stokes Equations Hantaek Bae Ulsan National Institute of Science and Technology (UNIST), Korea Recent Advances in Hydrodynamics, 216.6.9 Joint work with Eitan
More informationOptimal control in fluid mechanics by finite elements with symmetric stabilization
Computational Sciences Center Optimal control in fluid mechanics by finite elements with symmetric stabilization Malte Braack Mathematisches Seminar Christian-Albrechts-Universität zu Kiel VMS Worshop
More informationThe Kolmogorov Law of turbulence
What can rigorously be proved? IRMAR, UMR CNRS 6625. Labex CHL. University of RENNES 1, FRANCE Introduction Aim: Mathematical framework for the Kolomogorov laws. Table of contents 1 Incompressible Navier-Stokes
More informationFrom Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray
From Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray C. David Levermore Department of Mathematics and Institute for Physical Science and Technology University of Maryland, College Park
More information