A Finite Element Approximation of a Coupled Model for Two Turbulent Fluids
|
|
- Augustine Cobb
- 5 years ago
- Views:
Transcription
1 A Finite Element Approximation of a Coupled Model for Two Turbulent Fluids Driss Yakoubi: GIREF, Université Laval, Québec. Canada joint works with Tomas Chacon Rebollo: Laboratoire Jacques-Louis Lions, UPMC, Paris, France and Departamento EDAN, Universidad de Sevilla, Spain Journée FreeFEM++ December 7-10 Paris
2 Outline Introduction Continuous Scheme Finite Element Analysis Numerical results
3 Model equation It includes: A one-equation simplified turbulent model for each fluid: Only one statistic of turbulent is considered, the turbulent kinetic energy (TKE). Modeling of TKE generation of ocean-atmosphere interface by quadratic law Modeling of friction at interface by manning law
4 Model equation: RANS Reynolds Averaged Navier-Stokes r ( i (k i )ru i )+rp i = f i in, r u i = 0 in, r ( i (k i )rk i ) = i (k i ) ru i 2 in, u i = 0 on i, k i = 0 on i, i (k i )@ ni u i p i n i + apple i (u i u j ) u i u j = 0 on, 1 apple i 6= j apple 2, k i = u 1 u 2 2 on.
5 Model equation Open bounded sets of either convex or with boundary that: i, Data: R = i [ 2 C 1,1 are the viscosity and diffusion coefficients, such i 2 W 1,1 (R) ( apple i (`) apple 1, apple i (`) apple 1, and 0 i(`) apple 2, 0 i(`) apple 2, > 0, apple i > 0 f i 2 L 2 ( ) d are the friction coefficients. are the source terms Specific Purposes: Devise stable numerical methods to compute steady states, as simplified equilibra states of the ocean-atmosphere system.
6 Weak formulation of continuous model Functional spaces: Velocity spaces: X i = {v i 2 H 1 ( ); v i = 0 on i} Pressure spaces: L 2 0( )={q i 2 L 2 ( ); such that q i =0} TKE spaces: Y i = {k i 2 W 1,r0 ( ); k i = 0 on i} 1 r + 1 =1, and r>d r0
7 Weak formulation of continuous model Find (u i,p i,k i ) 2 X i L 2 ( ) W 1,r0 ( ) such that, for all (v i,q i, ' i ) 2 X i L 2 ( ) W 1,r 0 ( ), i (k i )ru i : rv i (r v i ) p i + apple i u i u j (u i u j ) v i d = f i v i (r u i ) q i =0, k i = 0 on i, k i = u i u j 2 on, and i(k i )rk i r' i = i (k i ) ru i 2 ' i.
8 Weak formulation of continuous model Since u i 2 X i then its trace on belongs to H 1 2 ( ),! L 3 ( ) d,then u i u j (u i u j ) v i d is well defined. Analysis of this model by [Bernardi, Chacon, Lewandowski and Murat] Difficulties: 1. ru i 2 2 L 1 ( ) d, 2. Coupling fuilds by: u i u j (u i u j ) v i, and k i = u i u j 2 on, 3. Coupling Eqs by terms: ru i 2, i ( ), and i ( )
9 Analysis of continuous scheme a i ki n ; u n+1 i, rv i + b i v i,p n+1 i + apple i u n+1 i u n+1 j u n+1 i u n+1 j v i = i(k n i )rk n+1 i r' i = i (ki n ) ru n+1 i 2 ' i f i v i, With: k n+1 i = 0 on i, and k n+1 i = u n+1 1 u n on. The iterative scheme is contractive if The turbulent diffusion is large enough with respect to the data, and The iterates enough: u n+1 i n and kn+1 i remains bounded in norms smooth n
10 Analysis of continuous scheme Theorem [T. Chacon Rebollo, S. Del Pino & DY] Assume that the sequences (u n remain bounded i ) n and (kn i ) n in W 1,3+" ( by M. i ) d and W 1,3 ( ) Then, there exists a constant C depending only on data such tha if the the iterative K scheme = C is contractive in the sens that: < 1 2X i=1 2X i=1 u n+1 i k n+1 i u n i k n i 2 1, apple K 2 1, apple K 2X i=1 2X i=1 k n i k n 1 i k n i k n 1 i 2 1,, and 2 1,.
11 Finite Element Approximation Bernardi-Chacon-Lewandowski-Murat: (Numer. Math, 2004) analysis of 2D F.E solution The velocity-pressure is discretized by the Mini-Element on both are supposed to be compatible on The velocity spaces X ih interface in the sens that the trace spaces ih = {v ih, for v ih 2 X ih }, i =1, 2 Theorem [Bernardi-Chacon-Lewandowski-Murat] In the 2D case, there exists a subsequence of the solutions (u 1h,p 1h,k 1h ), (u 2h,p 2h,k 2h )that converge strongly in: H 1 ( 1 ) 2 L 2 0( 1 ) H s ( 1 ) H 1 ( 2 ) 2 L 2 0( 2 ) H s ( 2 ), to a solution of model problem. are equal. for 0 apple s<1/2
12 Finite Element Approximation The problems we face now are: To build a more constructive solution schem To analyse the 3D. Discrete spaces: X i,h = nv i,h 2 C 0 d, 8K 2 Ti,h, 2 P 2 (K) d o \ X i M i,h = q i,h 2 L 2 ( ), 8K 2 T i,h,q i,h K 2 P 1 (K) K i,h = `i,h 2 C 0, 8K 2 T i,h, `i,h K 2 P 2 (K), `i,h i =0. The family (X i,h,m i,h ) h>0, for i =1, 2satisfy the Brezzi-Fortin: 8q i,h 2 M i,h, b i (v i,h,q i,h ) sup v i,h 2X i,h v i,h i,hkq i,h k 0, i
13 Finite Element Approximation Standard interpolation: S i,h : H 1 ( ) C 0! K i,h k i! S i,h (k i ). i,h : X i! X i,h such that v i! i,h (v i ). L i,h : H ( )! i,h such that L i,h (W )=(S i,h W ). This is in the sens that for all k i 2 H 1 ( ) C 0, the trace on of the interpolate S i,h (k i ) coincides with the interpolate of the trace L i,h (k i ).
14 Finite Element Approximation 1. Discrete algorithms: a i ki,h n ; un+1, v i,h i,h and 8q i,h 2 M i,h b i +b i v i,h,p n+1 i,h +apple i u n+1 i,h u n+1 j,h u n+1 i,h,q i,h =0, u n+1 i,h u n+1 j,h v i,h d = f i v i,h, 2. k n+1 i,h = 0 on i, k n+1 i,h = un+1 1,h u n+1 i(k n i,h) rk n+1 ih r' ih = 2,h 2 on, i (k n i,h) ru n+1 i,h 2 ' i,h, 8' i,h 2 K 0 i,h, where K 0 i,h = K i,h \ W 1,r 0 ( )
15 Finite Element Approximation Theorem [T. Chacon Rebollo & DY] Assume that the sequences (u n i ) n, (un ih ) n, (kn i ) n, and (kn ih ) n remain bounded in by M. W 1,3+" ( ) d and W 1,3 ( ) Then there exists a constants depending only on data such that: 2X k n i,h k n i 2 1, + u n i,h u n i 2 1, + kp n i,h p n i k 2 i=1 apple cn 1 Were Furthermore if, discrete scheme converges, and its limit is a solution= " > 1 2, 8" > 0. of the continuous c 1 < model. n + c h2 + h +( 2 + 1) h
16 Finite Element Approximation Keys of the proof: Use convenient choices of test functions: 2X For instance, to obtain the estimate for 1. We introduce the following space: i=1 k n+1 ih W i,h = ' i,h 2 C 0 (@ ); 8e 2 E i,h, ' i,h e 2 P 2 (e). 2. We introduce the lifting operator: R ih : W ih 7! K ih, (R ih (' ih i = ' ih, 8' ih 2 W ih. kr i,h (' i,h )k W 1,p ( ) apple ck' i,h k W 1 1/p,p (@ ). 3. and set the test function TKE: k n i 2 1, ' ih = `n+1 i,h R (`n+1 i,h i,h ), where `n+1 i,h = kn+1 i,h S i,h (k n+1 i ).
17 Finite Element Approximation Keys of the proof: Due to the friction term, it is needed to estimate: h u n+1 1,h u n+1 1 i h u n+1 2,h u n+1 2 u n+1 1,h + un+1 1 This is done using Grisvard s Theorem: i u n+1 2,h + un+1 2 H 1/2 00 ( ) Assume that s a bounded Lipschitz-continuous open subset of R d. Let s, s 1,s 2 0 and p, p 1,p 2 2 [1, +1) such that s 1 s, s 2 s and either 1 s 1 + s 2 s d p p 2 p 0, s i s>d p i p or s 1 + s 2 s>d 0, s i s d 1 p p 2 p 1 1 p i p Then the mapping (u, v) 7! uv is a continuous bilinear map from W s 1,p 1 ( ) W s 2,p 2 ( ) to W s,p ( ).
18 Numerical tests by FreeFEM++ We have tested ou iterative scheme for the data: (Atmosphere:) 1 =[0, 5] [0, 1] [0, 1], (Ocean:) 2 =[0, 5] [0, 1] [ 1, 0]. Turbulent diffusions: (realistic values) 1 (k 1 )= p k 1, Friction coefficients: (realistic values) i ( ) = i ( ). 1(k 2 )= p k 2 apple i = 10 3, = Boundary data: u 1 =(1, 0, 0) on y =1 u i =(0, 0, 0) on the remaining k i = 0 /.
19 Numerical tests by FreeFEM++ Discretization P2-P1 for velocity-pressure Discretization P2 for TKE Computations with FreeFEM++ In order to verify the convergence order, we use different size meshes.
20 Numerical tests by FreeFEM++: Results The algorithm converges to steady state with a rate > 0.25, in agreement with theoritical analysis: 2X i=1 k i,h k i 1, i + u i,h u i 1, i = E h apple c Where = 3 2 We set h = 3 3+" > 1 2, 8" > 0. log Eh E h/2 log(2) Mesh size h + h /2 + h 1/2 h h ---- h/ h/ h/ h/16 not yet! [in progress]
21 Numerical tests by FreeFEM++ Wind-induced flow on swimming pool (Atmosphere:) 1 =[0, 10 4 ] [0, ] [0, 500], The ocean domain is defined by: Horizontal dimensions (m): Bathemetry (m): 8 >< >:! =[0, 10 4 ] [0, ] 50 if 0 apple x apple x x 10 3 if apple x apple if apple x apple 10 4
22 We test the formation of the up-welling effect, due to the interaction between wind-tension and Coriolis forces. Numerical tests by FreeFEM++ Ocean Domain
23 Numerical tests by FreeFEM++ Methodology The system of equation we solve in t u i +(u i r) u i + ( u i,y,u i,x, 0) r ( i (k i )ru i )+rp i = f i in, r u i = 0 t k i + u i rk i r ( i (k i )rk i ) = i (k i ) ru i 2 in, u i = 0 on i, k i = 0 on i, i (k i )@ ni u i p i n i + apple i (u i u j ) u i u j = 0 on, 1 apple i 6= j apple 2, k i = u 1 u 2 2 on.
24 Numerical tests by FreeFEM++ Results:
25 Numerical tests by FreeFEM++ Velocity fields: z=-5m z = - 30 m z = - 10 m z = - 55 m
26 Numerical tests by FreeFEM++ Velocity fields: Movies
Simulation of free surface fluids in incompressible dynamique
Simulation of free surface fluids in incompressible dynamique Dena Kazerani INRIA Paris Supervised by Pascal Frey at Laboratoire Jacques-Louis Lions-UPMC Workshop on Numerical Modeling of Liquid-Vapor
More informationADAPTATION OF THE REYNOLDS STRESS TURBULENCE MODEL FOR ATMOSPHERIC SIMULATIONS
ADAPTATION OF THE REYNOLDS STRESS TURBULENCE MODEL FOR ATMOSPHERIC SIMULATIONS Radi Sadek 1, Lionel Soulhac 1, Fabien Brocheton 2 and Emmanuel Buisson 2 1 Laboratoire de Mécanique des Fluides et d Acoustique,
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 informationEXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 147, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE AND REGULARITY OF SOLUTIONS FOR
More informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationUniform inf-sup condition for the Brinkman problem in highly heterogeneous media
Uniform inf-sup condition for the Brinkman problem in highly heterogeneous media Raytcho Lazarov & Aziz Takhirov Texas A&M May 3-4, 2016 R. Lazarov & A.T. (Texas A&M) Brinkman May 3-4, 2016 1 / 30 Outline
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 Two-Phase Solid/Fluid Model for Dense Granular Flows including Dilatancy Effects
23-26th May 2016, GdR EGRIN, Piriac-sur-Mer A Two-Phase Solid/Fluid Model for Dense Granular Flows including Dilatancy Effects Anne Mangeney1,1 F. Bouchut2, E. Fernandez-Nieto3, G. Narbona-Reina3, and
More informationSome recent developments on ROMs in computational fluid dynamics
Some recent developments on ROMs in computational fluid dynamics F. Ballarin 1, S. Ali 1, E. Delgado 2, D. Torlo 1,3, T. Chacón 2, M. Gómez 2, G. Rozza 1 1 mathlab, Mathematics Area, SISSA, Trieste, Italy
More informationGlobal well-posedness of the primitive equations of oceanic and atmospheric dynamics
Global well-posedness of the primitive equations of oceanic and atmospheric dynamics Jinkai Li Department of Mathematics The Chinese University of Hong Kong Dynamics of Small Scales in Fluids ICERM, Feb
More informationFiltered scheme and error estimate for first order Hamilton-Jacobi equations
and error estimate for first order Hamilton-Jacobi equations Olivier Bokanowski 1 Maurizio Falcone 2 2 1 Laboratoire Jacques-Louis Lions, Université Paris-Diderot (Paris 7) 2 SAPIENZA - Università di Roma
More informationSpectral discretization of Darcy s equations coupled with the heat equation
Spectral discretization of Darcy s equations coupled with the heat equation Christine Bernardi, Sarra Maarouf, Driss Yakoubi To cite this version: Christine Bernardi, Sarra Maarouf, Driss Yakoubi. Spectral
More informationModel Studies on Slag-Metal Entrainment in Gas Stirred Ladles
Model Studies on Slag-Metal Entrainment in Gas Stirred Ladles Anand Senguttuvan Supervisor Gordon A Irons 1 Approach to Simulate Slag Metal Entrainment using Computational Fluid Dynamics Introduction &
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 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 informationChapter 10. Solids and Fluids
Chapter 10 Solids and Fluids Surface Tension Net force on molecule A is zero Pulled equally in all directions Net force on B is not zero No molecules above to act on it Pulled toward the center of the
More informationLeast-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations
Least-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations by Wilhelm Heinrichs Universität Duisburg Essen, Ingenieurmathematik Universitätsstr.
More informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationJuan Vicente Gutiérrez Santacreu Rafael Rodríguez Galván. Departamento de Matemática Aplicada I Universidad de Sevilla
Doc-Course: Partial Differential Equations: Analysis, Numerics and Control Research Unit 3: Numerical Methods for PDEs Part I: Finite Element Method: Elliptic and Parabolic Equations Juan Vicente Gutiérrez
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
More informationNumerical control and inverse problems and Carleman estimates
Numerical control and inverse problems and Carleman estimates Enrique FERNÁNDEZ-CARA Dpto. E.D.A.N. - Univ. of Sevilla from some joint works with A. MÜNCH Lab. Mathématiques, Univ. Blaise Pascal, C-F 2,
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationSecond-order finite-element projection method for 3D flows
Second-order finite-element projection method for 3D flows J.-L. Guermond I and L. Quartapelle 2 1 LIMSI-CNRS, BP 133, 91403 Orsay, France, (guermond@limsi.fr) 2 Dipartimento di Fisica del Politecnico
More informationEffect of near-wall treatments on airflow simulations
Proceedings of 29 International Conference on Computational Methods for Energy Engineering and Environment: ICCM3E. Sousse, Tunisia, 2-22 November, 29, pp. Effect of near-wall treatments on airflow simulations
More informationSplitting Techniques in the Face of Huge Problem Sizes: Block-Coordinate and Block-Iterative Approaches
Splitting Techniques in the Face of Huge Problem Sizes: Block-Coordinate and Block-Iterative Approaches Patrick L. Combettes joint work with J.-C. Pesquet) Laboratoire Jacques-Louis Lions Faculté de Mathématiques
More informationRome - May 12th Université Paris-Diderot - Laboratoire Jacques-Louis Lions. Mean field games equations with quadratic
Université Paris-Diderot - Laboratoire Jacques-Louis Lions Rome - May 12th 2011 Hamiltonian MFG Hamiltonian on the domain [0, T ] Ω, Ω standing for (0, 1) d : (HJB) (K) t u + σ2 2 u + 1 2 u 2 = f (x, m)
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 informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationSobolev Spaces 27 PART II. Review of Sobolev Spaces
Sobolev Spaces 27 PART II Review of Sobolev Spaces Sobolev Spaces 28 SOBOLEV SPACES WEAK DERIVATIVES I Given R d, define a multi index α as an ordered collection of integers α = (α 1,...,α d ), such that
More informationR09. d water surface. Prove that the depth of pressure is equal to p +.
Code No:A109210105 R09 SET-1 B.Tech II Year - I Semester Examinations, December 2011 FLUID MECHANICS (CIVIL ENGINEERING) Time: 3 hours Max. Marks: 75 Answer any five questions All questions carry equal
More informationcentrifugal acceleration, whose magnitude is r cos, is zero at the poles and maximum at the equator. This distribution of the centrifugal acceleration
Lecture 10. Equations of Motion Centripetal Acceleration, Gravitation and Gravity The centripetal acceleration of a body located on the Earth's surface at a distance from the center is the force (per unit
More informationShih-sen Chang, Yeol Je Cho, and Haiyun Zhou
J. Korean Math. Soc. 38 (2001), No. 6, pp. 1245 1260 DEMI-CLOSED PRINCIPLE AND WEAK CONVERGENCE PROBLEMS FOR ASYMPTOTICALLY NONEXPANSIVE MAPPINGS Shih-sen Chang, Yeol Je Cho, and Haiyun Zhou Abstract.
More informationLES of turbulent shear flow and pressure driven flow on shallow continental shelves.
LES of turbulent shear flow and pressure driven flow on shallow continental shelves. Guillaume Martinat,CCPO - Old Dominion University Chester Grosch, CCPO - Old Dominion University Ying Xu, Michigan State
More informationLecture 7 Boundary Layer
SPC 307 Introduction to Aerodynamics Lecture 7 Boundary Layer April 9, 2017 Sep. 18, 2016 1 Character of the steady, viscous flow past a flat plate parallel to the upstream velocity Inertia force = ma
More informationTheoretical and numerical results for a chemo-repulsion model with quadratic production
Theoretical and numerical results for a chemo-repulsion model with quadratic production F. Guillén-Gonzalez, M. A. Rodríguez-Bellido & and D. A. Rueda-Gómez Dpto. Ecuaciones Diferenciales y Análisis Numérico
More informationNumerical simulations of heat transfer in plane channel flow
Numerical simulations of heat transfer in plane channel flow Najla EL GHARBI 1, 3, a, Rafik ABSI 2, b and Ahmed BENZAOUI 3, c 1 Renewable Energy Development Center, BP 62 Bouzareah 163 Algiers, Algeria
More informationd v 2 v = d v d t i n where "in" and "rot" denote the inertial (absolute) and rotating frames. Equation of motion F =
Governing equations of fluid dynamics under the influence of Earth rotation (Navier-Stokes Equations in rotating frame) Recap: From kinematic consideration, d v i n d t i n = d v rot d t r o t 2 v rot
More informationTurbulence. 2. Reynolds number is an indicator for turbulence in a fluid stream
Turbulence injection of a water jet into a water tank Reynolds number EF$ 1. There is no clear definition and range of turbulence (multi-scale phenomena) 2. Reynolds number is an indicator for turbulence
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationGeneralized Newtonian Fluid Flow through a Porous Medium
Generalized Newtonian Fluid Flow through a Porous Medium V.J. Ervin Hyesuk Lee A.J. Salgado April 24, 204 Abstract We present a model for generalized Newtonian fluid flow through a porous medium. In the
More informationNUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD
NUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD V. G. Guedes a, G. C. R. Bodstein b, and M. H. Hirata c a Centro de Pesquisas de Energia Elétrica Departamento de Tecnologias
More informationAN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract
nd Workshop on CFD Uncertainty Analysis - Lisbon, 19th and 0th October 006 AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION Alfredo Iranzo 1, Jesús Valle, Ignacio Trejo 3, Jerónimo
More informationSimulations for Enhancing Aerodynamic Designs
Simulations for Enhancing Aerodynamic Designs 2. Governing Equations and Turbulence Models by Dr. KANNAN B T, M.E (Aero), M.B.A (Airline & Airport), PhD (Aerospace Engg), Grad.Ae.S.I, M.I.E, M.I.A.Eng,
More informationPressure Losses for Fluid Flow Through Abrupt Area. Contraction in Compact Heat Exchangers
Pressure Losses for Fluid Flow Through Abrupt Area Contraction in Compact Heat Exchangers Undergraduate Research Spring 004 By Bryan J. Johnson Under Direction of Rehnberg Professor of Ch.E. Bruce A. Finlayson
More informationModel adaptation in hierarchies of hyperbolic systems
Model adaptation in hierarchies of hyperbolic systems Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France February 15th, 2012 DFG-CNRS Workshop Nicolas Seguin (LJLL, UPMC) 1 / 29 Outline of the
More informationA simple kinetic equation of swarm formation: blow up and global existence
A simple kinetic equation of swarm formation: blow up and global existence Miroslaw Lachowicz, Henryk Leszczyński, Martin Parisot To cite this version: Miroslaw Lachowicz, Henryk Leszczyński, Martin Parisot.
More informationGlobal Carleman inequalities and theoretical and numerical control results for systems governed by PDEs
Global Carleman inequalities and theoretical and numerical control results for systems governed by PDEs Enrique FERNÁNDEZ-CARA Dpto. E.D.A.N. - Univ. of Sevilla joint work with A. MÜNCH Lab. Mathématiques,
More informationA LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES
A LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES P. HANSBO Department of Applied Mechanics, Chalmers University of Technology, S-4 96 Göteborg, Sweden E-mail: hansbo@solid.chalmers.se
More informationESS Turbulence and Diffusion in the Atmospheric Boundary-Layer : Winter 2017: Notes 1
ESS5203.03 - Turbulence and Diffusion in the Atmospheric Boundary-Layer : Winter 2017: Notes 1 Text: J.R.Garratt, The Atmospheric Boundary Layer, 1994. Cambridge Also some material from J.C. Kaimal and
More informationGENERALISATION OF THE TWO-SCALE MOMENTUM THEORY FOR COUPLED WIND TURBINE/FARM OPTIMISATION
25 th National Symposium on Wind Engineering, Tokyo, Japan, 3-5 December 2018 第 25 回風工学シンポジウム (2018) GENERALISATION OF THE TWO-SCALE MOMENTUM THEORY FOR COUPLED WIND TURBINE/FARM OPTIMISATION Takafumi
More informationNumerical methods for Mean Field Games: additional material
Numerical methods for Mean Field Games: additional material Y. Achdou (LJLL, Université Paris-Diderot) July, 2017 Luminy Convergence results Outline 1 Convergence results 2 Variational MFGs 3 A numerical
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 informationOPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS
PORTUGALIAE MATHEMATICA Vol. 59 Fasc. 2 2002 Nova Série OPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS J. Saint Jean Paulin and H. Zoubairi Abstract: We study a problem of
More informationLecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2)
Lecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2) The ABL, though turbulent, is not homogeneous, and a critical role of turbulence is transport and mixing of air properties, especially in the
More information.u= 0 ρ( u t +(u. )u)= ρ g p+.[µ( u+ t u)]
THETIS is a numerical simulation tool developed by University of Bordeaux. It is a versatile code to solve different problems: fluid flows, heat transfers, scalar transports or porous mediums. The potential
More informationA posteriori error estimates applied to flow in a channel with corners
Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationAcceleration of a Domain Decomposition Method for Advection-Diffusion Problems
Acceleration of a Domain Decomposition Method for Advection-Diffusion Problems Gert Lube 1, Tobias Knopp 2, and Gerd Rapin 2 1 University of Göttingen, Institute of Numerical and Applied Mathematics (http://www.num.math.uni-goettingen.de/lube/)
More informationChapter 7 The Time-Dependent Navier-Stokes Equations Turbulent Flows
Chapter 7 The Time-Dependent Navier-Stokes Equations Turbulent Flows Remark 7.1. Turbulent flows. The usually used model for turbulent incompressible flows are the incompressible Navier Stokes equations
More informationCFD Analysis of Forced Convection Flow and Heat Transfer in Semi-Circular Cross-Sectioned Micro-Channel
CFD Analysis of Forced Convection Flow and Heat Transfer in Semi-Circular Cross-Sectioned Micro-Channel *1 Hüseyin Kaya, 2 Kamil Arslan 1 Bartın University, Mechanical Engineering Department, Bartın, Turkey
More informationMean Field Games on networks
Mean Field Games on networks Claudio Marchi Università di Padova joint works with: S. Cacace (Rome) and F. Camilli (Rome) C. Marchi (Univ. of Padova) Mean Field Games on networks Roma, June 14 th, 2017
More informationLocal null controllability of the N-dimensional Navier-Stokes system with N-1 scalar controls in an arbitrary control domain
Local null controllability of the N-dimensional Navier-Stokes system with N-1 scalar controls in an arbitrary control domain Nicolás Carreño Université Pierre et Marie Curie-Paris 6 UMR 7598 Laboratoire
More informationA Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract
A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin
More informationME 305 Fluid Mechanics I. Part 8 Viscous Flow in Pipes and Ducts. Flow in Pipes and Ducts. Flow in Pipes and Ducts (cont d)
ME 305 Fluid Mechanics I Flow in Pipes and Ducts Flow in closed conduits (circular pipes and non-circular ducts) are very common. Part 8 Viscous Flow in Pipes and Ducts These presentations are prepared
More informationCFD Based Optimization of a Vertical Axis Wind Turbine
SCHOOL OF ENGINEERING INSTITUTE FOR ENERGY SYSTEMS CFD Based Optimization of a Vertical Axis Wind Turbine Emmanouil Falagkaris Supervisory Team Prof. David Ingram, Dr Ignazio Maria Viola, Dr Daniel Friedrich
More informationPseudo-divergence-free element free Galerkin method for incompressible fluid flow
International Workshop on MeshFree Methods 23 Pseudo-divergence-free element free Galerkin method for incompressible fluid flow Y. Vidal () and A. Huerta (2) Abstract: Incompressible modelling in finite
More informationAn Equal-order DG Method for the Incompressible Navier-Stokes Equations
An Equal-order DG Method for the Incompressible Navier-Stokes Equations Bernardo Cockburn Guido anschat Dominik Schötzau Journal of Scientific Computing, vol. 40, pp. 188 10, 009 Abstract We introduce
More informationA Dykstra-like algorithm for two monotone operators
A Dykstra-like algorithm for two monotone operators Heinz H. Bauschke and Patrick L. Combettes Abstract Dykstra s algorithm employs the projectors onto two closed convex sets in a Hilbert space to construct
More informationInverse function theorems
Inverse function theorems Ivar Ekeland CEREMADE, Université Paris-Dauphine Ravello, May 2011 Ivar Ekeland (CEREMADE, Université Paris-Dauphine) Inverse function theorems Ravello, May 2011 1 / 18 Ekeland
More informationRoad to Turbulence in Stratified Flow
Road to Turbulence in Stratified Flow (Exaple: Boundary flow with f =0) Lainar Flow U H Turbulent Flow I II III IV V u u p u i i i u b t x x x x b g 0 i N 0 0 g z Turbulence occurs when: II UL R e R eyonlds
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 informationUNIT 4 FORCES ON IMMERSED BODIES. Lecture-01
1 UNIT 4 FORCES ON IMMERSED BODIES Lecture-01 Forces on immersed bodies When a body is immersed in a real fluid, which is flowing at a uniform velocity U, the fluid will exert a force on the body. The
More informationA posteriori error analysis of an augmented mixed finite element method for Darcy flow
A posteriori error analysis of an augmented mixed finite element method for Darcy flow Tomás P. Barrios, J. Manuel Cascón and María González Abstract We develop an a posteriori error analysis of residual
More informationTutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace
Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and
More informationA two-fluid model of turbulent two-phase flow for simulating turbulent stratified flows
Ocean Engineering 30 (2003) 153 161 www.elsevier.com/locate/oceaneng A two-fluid model of turbulent two-phase flow for simulating turbulent stratified flows Y.M. Shen a,, C.-O. Ng b, A.T. Chwang b a State
More informationOptimized Surface Acoustic Waves Devices With FreeFem++ Using an Original FEM/BEM Numerical Model
Optimized Surface Acoustic Waves Devices With FreeFem++ Using an Original FEM/BEM Numerical Model P. Ventura*, F. Hecht**, Pierre Dufilié*** *PV R&D Consulting, Nice, France Laboratoire LEM3, Université
More informationAn iterative algorithm for nonlinear wavelet thresholding: Applications to signal and image processing
An iterative algorithm for nonlinear wavelet thresholding: Applications to signal and image processing Marie Farge, LMD-CNRS, ENS, Paris Kai Schneider, CMI, Université de Provence, Marseille Alexandre
More informationRegularization in Reproducing Kernel Banach Spaces
.... Regularization in Reproducing Kernel Banach Spaces Guohui Song School of Mathematical and Statistical Sciences Arizona State University Comp Math Seminar, September 16, 2010 Joint work with Dr. Fred
More informationEquivalent electrostatic capacitance Computation using FreeFEM++
Equivalent electrostatic capacitance Computation using FreeFEM++ P. Ventura*, F. Hecht** *PV R&D Consulting, Nice, France **Laboratoire Jacques Louis Lions, Université Pierre et Marie Curie, Paris, France
More informationSHARP BOUNDS FOR PROBABILITIES WITH GIVEN SHAPE INFORMATION
R u t c o r Research R e p o r t SHARP BOUNDS FOR PROBABILITIES WITH GIVEN SHAPE INFORMATION Ersoy Subasi a Mine Subasi b András Prékopa c RRR 4-006, MARCH, 006 RUTCOR Rutgers Center for Operations Research
More informationTomás Chacón Rebollo 1, Roger Lewandowski 2 and Eliseo Chacón Vera 1
Mathematical Modelling and Numerical Analysis M2AN, Vol. 34, N o 3, 2000, pp. 525 537 Modélisation Mathématique et Analyse Numérique ANALYSIS OF THE HYDROSTATIC APPROXIMATION IN OCEANOGRAPHY WITH COMPRESSION
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationME 305 Fluid Mechanics I. Chapter 8 Viscous Flow in Pipes and Ducts
ME 305 Fluid Mechanics I Chapter 8 Viscous Flow in Pipes and Ducts These presentations are prepared by Dr. Cüneyt Sert Department of Mechanical Engineering Middle East Technical University Ankara, Turkey
More information7.4 The Saddle Point Stokes Problem
346 CHAPTER 7. APPLIED FOURIER ANALYSIS 7.4 The Saddle Point Stokes Problem So far the matrix C has been diagonal no trouble to invert. This section jumps to a fluid flow problem that is still linear (simpler
More informationEfficient numerical solution of the Biot poroelasticity system in multilayered domains
Efficient numerical solution of the Biot poroelasticity system Anna Naumovich Oleg Iliev Francisco Gaspar Fraunhofer Institute for Industrial Mathematics Kaiserslautern, Germany, Spain Workshop on Model
More informationA posteriori error estimation in the FEM
A posteriori error estimation in the FEM Plan 1. Introduction 2. Goal-oriented error estimates 3. Residual error estimates 3.1 Explicit 3.2 Subdomain error estimate 3.3 Self-equilibrated residuals 3.4
More informationCompressible Flow Modeling Occurring in a Depressurization Process
Compressible Flow Modeling Occurring in a Depressurization Process Frédéric Viry 1, Vincent Bruyere 1, Thomas Paris 2, Patrick Namy 1 1 SIMTEC, 8 rue Duployé, Grenoble, 38100, France, +33 9 53 51 45 60
More informationGéraldine Pichot 1. Introduction
ESAIM: PROCEEDINGS, October 2007, Vol. 22, 132-139 Gabriel Caloz & Monique Dauge, Editors PROCEEDINGS IN FLOW MODELLING IN THE SURROUNDINGS OF THE END PART OF A TRAWL Géraldine Pichot 1 Abstract. This
More informationA note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations
A note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations Bernardo Cockburn Guido anschat Dominik Schötzau June 1, 2007 Journal of Scientific Computing, Vol. 31, 2007, pp.
More informationElectrokinetic effects in the breakup of electrified jets: a Volume-Of-Fluid numerical study
Electrokinetic effects in the breakup of electrified jets: a Volume-Of-Fluid numerical study J. M. Lopez-Herrera 1, A. M. Gañan-Calvo 1, S. Popinet 2, M. A. Herrada 1 1.-University of Sevilla, Spain. 2.-Université
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More informationAbstract. 1 Introduction
Mass Conservation in shallow water storm surge computations Jeffrey P. Laible Department of Civil and Environmental Engineering University of Vermont Burlington, Vermont 05405 E-mail: laible@emba. uvm.
More informationWater Artificial Circulation For Eutrophication Control
Water Artificial Circulation For Eutrophication Control Francisco Javier Fernández (fjavier.fernandez@cud.uvigo.es) Centro Universitario de la Defensa, Escuela Naval Militar, 36920 Marn, Spain Joint work
More informationTHE SKOROKHOD OBLIQUE REFLECTION PROBLEM IN A CONVEX POLYHEDRON
GEORGIAN MATHEMATICAL JOURNAL: Vol. 3, No. 2, 1996, 153-176 THE SKOROKHOD OBLIQUE REFLECTION PROBLEM IN A CONVEX POLYHEDRON M. SHASHIASHVILI Abstract. The Skorokhod oblique reflection problem is studied
More informationFINITE ELEMENT APPROXIMATION OF ELLIPTIC DIRICHLET OPTIMAL CONTROL PROBLEMS
Numerical Functional Analysis and Optimization, 28(7 8):957 973, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0163-0563 print/1532-2467 online DOI: 10.1080/01630560701493305 FINITE ELEMENT APPROXIMATION
More informationSpectral discretization of a model for organic pollution in waters
Spectral discretization of a model for organic pollution in waters Christine Bernardi, Driss Yakoubi To cite this version: Christine Bernardi, Driss Yakoubi. Spectral discretization of a model for organic
More informationC e n t r u m v o o r W i s k u n d e e n I n f o r m a t i c a
C e n t r u m v o o r W i s k u n d e e n I n f o r m a t i c a Modelling, Analysis and Simulation Modelling, Analysis and Simulation Bilinear forms for the recovery-based discontinuous Galerkin method
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationNONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS
NONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS J.-L. GUERMOND 1, Abstract. This paper analyzes a nonstandard form of the Stokes problem where the mass conservation
More information