TRANSPORT IN POROUS MEDIA
|
|
- Pierce Tucker
- 5 years ago
- Views:
Transcription
1 1 TRANSPORT IN POROUS MEDIA G. ALLAIRE CMAP, Ecole Polytechnique 1. Introduction 2. Main result in an unbounded domain 3. Asymptotic expansions with drift 4. Two-scale convergence with drift 5. The case of bounded domains Ecole CEA-EDF-INRIA, Décembre 20
2 2 -I- INTRODUCTION Based on a joint work with A. Mikelic and A. Piatnitski. Infinite porous medium: (connected) fluid part Ω f, solid part Ω s = R n \ Ω f. Saturated incompressible single phase flow in Ω f and a single solute. The unknown is the concentration u in the fluid. convection diffusion in the fluid: u τ + b yu div y (D y u) = 0 in Ω f (0, T ), Given incompressible and steady-state velocity: div b = 0 in Ω f no-flux boundary condition on the pore boundaries: b n = 0 and D y u n = 0 on Ω f (0, T ),
3 3 Scaling We want to upscale this model, so we define a large macroscopic scale ǫ 1 and we choose a long time scale of order ǫ 2 (parabolic or diffusion scaling) x = ǫy and t = ǫ 2 τ. We define u ǫ (t, x) = u(τ, y) which is a solution of u ǫ t + 1 ǫ b ǫ x u ǫ div x (D ǫ x u ǫ ) = 0 in Ω ǫ (0, T) with T = ǫ 2 T. u ǫ (x, 0) = u 0 (x), x Ω ǫ, D ǫ x u ǫ n = 0 on Ω ǫ (0, T)
4 4 Periodicity assumption ε Ω Periodic unit cell Y = (0, 1) n = Y O with fluid part Y Periodic (infinite) porous media x Ω ǫ y Y ( x ) Stationary incompressible periodic flow b ǫ (x) = b ǫ Y and b n = 0 on O ( x ) Periodic symmetric coercive diffusion D ǫ (x) = D ǫ with div y b = 0 in
5 5 Another approach to scaling: dimensional analysis We write the same equations with dimensional constants denoted by * : c t + b x c div x (D x c ) = 0 in Ω f (0, T ), D x c n = 0 on Ω f (0, T ),
6 6 Dimensional analysis (ctd.) We adimensionalize the equations as follows: Characteristic lengthscale L R and timescale T R. Period l << L R : we introduce a small parameter ǫ = l L R. Characteristic velocity b R. Characteristic concentration c R. Characteristic diffusivity D R. New adimensionalized variables and functions: x = x, t = t, b ǫ (x, t) = b (x, t ), D = D, u ǫ = c L R T R b R D R c R
7 7 Dimensional analysis (ctd.) Dimensionless equation u ǫ t + b RT R b ǫ x u ǫ D RT R L R L 2 div x (D x u ǫ ) = 0 in Ω ǫ (0, T) R and D x u ǫ n = 0 on Ω ǫ (0, T). We choose a diffusion timescale, i.e., we assume T R = L2 R D R. Péclet number: Pe = L Rb R D R. We assume Pe = ǫ 1. u ǫ t + Pe b ǫ x u ǫ div x (D x u ǫ ) = 0 in Ω ǫ (0, T)
8 8 Microscopic model u ǫ t + 1 ǫ b ǫ x u ǫ div x (D ǫ x u ǫ ) = 0 in Ω ǫ (0, T) u ǫ (x, 0) = u 0 (x), x Ω ǫ, D ǫ x u ǫ n = 0 on Ω ǫ (0, T) Assumptions: Stationary incompressible periodic flow div y b = 0 in Y, b n = 0 on O Periodic symmetric coercive diffusion D Goal of homogenization: find the effective diffusion tensor.
9 9 -II- MAIN RESULT Theorem. The solution u ǫ satisfies u ǫ (t, x) u (t, x b ǫ ) t with the effective drift b = 1 Y b(y) dy Y and u the solution of the homogenized problem u t div (A u) = 0 in R n (0, T) u(t = 0, x) = u 0 (x) in R n
10 10 Precise convergence u ǫ (t, x) = u (t, x b ǫ ) t + r ǫ (t, x) with T lim ǫ 0 0 R n r ǫ (t, x) 2 dt dx = 0,
11 11 Homogenized diffusion tensor A ij = D(e i + y χ i ) (e j + y χ j )dy Y where χ i (y), 1 i n, are solutions of the cell problems b(y) y χ i div y (D(y) ( y χ i + e i )) = (b b(y)) e i D(y) ( y χ i + e i ) n = 0 on O y χ i (y) Y -periodic in Y Remark that the value of b is exactly the compatibility condition for the existence of χ i.
12 12 Define ũ ǫ (t, x) = u Equivalent homogenized equation (t, x b ǫ t ). Then, it is solution of ũ ǫ t + 1 ǫ b ũ ǫ div (A ũ ǫ ) = 0 in R n (0, T) ũ ǫ (t = 0, x) = u 0 (x) in R n
13 13 -III- TWO-SCALE ANSATZ WITH DRIFT To motivate our result, let us start with a formal process. Standard two-scale asymptotic expansions should be modified to introduce an unknown large drift b R n u ǫ (t, x) = + i=0 ( ǫ i u i t, x b t ǫ, x ), ǫ with u i (t, x, y) a function of the macroscopic variable x and of the periodic microscopic variable y Y = (0, 1) n.
14 14 We plug these ansatz in the system of equations and use the usual chain rule derivation ( (u i t, x b t ǫ, x )) = ( ) ( ǫ 1 y u i + x u i t, x b t ǫ ǫ, x ), ǫ plus a new contribution t ( (u i t, x b t ǫ, x )) ǫ = u i t ǫ 1 b x u }{{} i ( t, x, x ǫ new term )
15 15 Fredholm alternative in the unit cell Lemma. The boundary value problem b(y) y χ div y (D(y) y χ) = f D(y) y χ n + g = 0 on O y χ(y) Y -periodic in Y admits a unique solution χ H 1 (Y )/R (up to an additive constant), if and only if f(y) dy + g(y) ds = 0. Y O
16 16 Variational formulation of the cell problem b(y) y χ(y)φ(y) dy + D y χ(y) y φ(y) dy = Y Y f(y)φ(y) dy + g(y)φ(y) ds Y O for any test function φ H 1 (Y ). Coercive bilinear form on the orthogonal subspace to its kernel R.
17 17 Cascade of equations Equation of order ǫ 2 : b(y) y u 0 div y (D(y) y u 0 ) = 0 in Y D(y) y u 0 n = 0 on O y u 0 (t, x, y) Y -periodic We deduce u 0 (t, x, y) u(t, x)
18 18 Equation of order ǫ 1 : b x u 0 + b(y) ( x u 0 + y u 1 ) div y (D(y) ( x u 0 + y u 1 )) = 0 in Y D ( x u 0 + y u 1 ) n = 0 on O y u 1 (t, x, y) Y -periodic We deduce u 1 (t, x, y) = n i=1 u 0 x i (t, x)χ i (y)
19 19 Cell problem b(y) y χ i div y (D(y) ( y χ i + e i )) = (b b(y)) e i in Y D(y) ( y χ i + e i ) n = 0 on O y χ i (y) Y -periodic The compatibility condition (Fredholm alternative) for the existence of χ i is b = 1 Y b(y) dy Y
20 20 Equation of order ǫ 0 : u t b x u 1 + b ( x u 1 + y u 2 ) div y (D ( x u 1 + y u 2 )) div x (D( y u 1 + x u)) = 0 in Y D ( x u 1 + y u 2 ) n = 0 on O y u 2 (t, x, y) Y -periodic Compatibility condition for the existence of u 2 homogenized problem. (Recall that u 1 is linear in x u.) u t div (A u) = 0 in R n (0, T) u(t = 0, x) = u 0 (x) in R n,
21 21 -IV- RIGOROUS PROOF The proof is made of 3 steps 1. A priori estimates. 2. Passing to the limit by two-scale convergence with drift. 3. Strong convergence.
22 22 A priori estimates Assume u 0 L 2 (R n ). For any final time T > 0, there exists a constant C > 0 that does not depend on ǫ such that u ǫ L (0,T;L 2 (Ω ǫ )) + u ǫ L 2 ((0,T) Ω ǫ ) C Proof. Multiply the fluid equation by u ǫ and integrate by parts to get the usual parabolic estimates.
23 23 Usual two-scale convergence Proposition. Let w ǫ be a bounded sequence in L 2 (R n ). Up to a subsequence, there exist a limit w(x, y) L 2 (R n T n ) such that w ǫ two-scale converges to w in the sense that ( lim w ǫ (x)φ x, x ) dx = ǫ 0 R ǫ n for all functions φ(x, y) L 2 (R n ; C(T n )). R n T n w(x, y)φ(x, y) dxdy
24 24 Two-scale convergence with drift Proposition (Marusic-Paloka, Piatnitski). Let V R N be a given drift velocity. Let w ǫ be a bounded sequence in L 2 ((0, T) R n ). Up to a subsequence, there exist a limit w 0 (t, x, y) L 2 ((0, T) R n T n ) such that w ǫ two-scale converges with drift weakly to w 0 in the sense that T ( lim w ǫ (t, x)φ t, x + V ǫ 0 0 R ǫ t, x ) dt dx = ǫ n T w 0 (t, x, y)φ(t, x, y) dt dxdy 0 R n T n for all functions φ(t, x, y) L 2 ((0, T) R n ; C(T n )).
25 25 Lemma. Let φ(t, x, y) L 2 ( (0, T) T N ; C c (R N ) ). Then T lim ǫ 0 0 R N (t, φ x + V ( x ) ) 2 ǫ t, dt dx = ǫ T 0 R N T N φ(t, x, y) 2 dt dxdy. Proof. Change of variables x = x + V ǫ t T 0 R N (t, φ x + V ( x ǫ t, ǫ ) ) 2 dt dx = T 0 R N ( φ t, x, x ǫ V ) 2 ǫ 2 t dt dx We mesh R N with cubes of size ǫ, R N = i Z Yi ǫ with Yi ǫ = x ǫ i R N ( φ t, x, x = i Z ǫ V ) 2 ǫ 2 t dx = i Y ǫ i ǫ N T N φ (x ǫ i, y) 2 dy + o(1) = + (0, ǫ)n ( φ t, x ǫ i, x ǫ V ) 2 ǫ 2 t dx + o(1) Ω Y φ (x, y) 2 dxdy + o(1)
26 26 Passing to the limit We multiply the equation by an oscillating test function with drift V = b Ψ ǫ = φ (t, x + Vǫ ) ( t + ǫ φ 1 t, x + V ǫ t, x ) ǫ and we use the two-scale convergence with drift to get the homogenized equation.
27 27 STRONG CONVERGENCE We use the notion of strong two-scale convergence with drift. Proposition. If w ǫ (t, x) two-scale converges with drift weakly to w 0 (t, x, y) (assumed to be smooth enough) and lim w ǫ L2 ((0,T) R ǫ 0 n ) = w 0 L2 ((0,T) R n T n ), then it converges strongly in the sense that w ǫ(t, x) w 0 (t, x b ǫ t, x ) ǫ lim ǫ 0 L 2 ((0,T) R n ) = 0
28 28 -V- NUMERICAL RESULTS Numerical computations done with FreeFem++ in 2-d for circular obstacles. The velocity b(y) is generated by solving the following filtration problem in the fluid part Y of the unit cell Y y p y b = e i in Y ; div y b = 0 in Y ; b = 0 p, b are Y periodic. on O;
29 29
30 Log-log plot of the longitudinal dispersion A 11 as a function of the local Péclet s number (asymptotic slope 1.7).
31 31 -VI- THE CASE OF BOUNDED DOMAINS Consider now a bounded domain Ω with a Dirichlet boundary condition on Ω. To study the long time behavior we consider the eigenvalue problem 1 ( x ) ( x ) ǫ b u ǫ div(d u ǫ ) = λ ǫ u ǫ in Ω ǫ ǫ ǫ u ǫ = 0 on Ω Ω ǫ, ( x ) D u ǫ n = 0 on Ω ǫ \ Ω ǫ where λ ǫ and u ǫ are the first eigenvalue and eigenfunction (which exists by the Krein-Rutman theorem). Assumptions: Stationary incompressible periodic flow div y b = 0 in Y, b n = 0 on O Periodic symmetric coercive diffusion D
32 Typical expected behavior of the first eigenfunction: because of the large drift it behaves like a boundary layer.
33 33 Following Bensoussan-Lions-Papanicolaou and Capdeboscq, for θ R N, we introduce the spectral cell problem div y (D(y) y ψ θ ) + b(y) y ψ θ = λ(θ)ψ θ in Y D(y) y ψ θ n = 0 on O y ψ θ (y)e θ y Y -periodic where λ(θ) and ψ θ are the first eigenvalue and eigenfunction. We also introduce the adjoint spectral cell problem div y (D(y) y ψθ) div y (b(y)ψθ) = λ(θ)ψθ D(y) y ψθ n = 0 on O y ψθ (y)e+θ y Y -periodic in Y The first eigenfunctions ψ θ and ψθ can be chosen positive and normalized by ψ θ (y)e θ y 2 dy = 1 and ψ θ (y)ψθ(y) dy = 1 T N T N
34 34 Lemma. The function θ λ(θ) is strictly concave from R N into R and admits a maximum λ which is obtained for a unique θ = θ. Denoting ψ = ψ θ and ψ = ψ θ, the vector field b(y) = ψ ψ b(y) + ψ D y ψ (y) ψ D y ψ (y) satisfies div y b = 0 in T N, T N b(y) dy = 0.
35 35 Change of unknown Define a new unknown function ũ ǫ (x) = u ǫ(x) ( ψ x ) ǫ ( and multiply the equation by ψ x ) ǫ. We obtain 1 x ) ũ ǫ div( ǫ b( ǫ D ( x ) ( x ) ũ ǫ ) = µ ǫ (ψ ψ ǫ ) ǫ ũ ǫ = 0 on Ω Ω ǫ, ( x ) D ũ ǫ n = 0 on Ω ǫ \ Ω ǫ ũ ǫ in Ω ǫ with D(y) = ψ (y)ψ (y)d(y) and µ ǫ = λ ǫ λ(θ ) ǫ 2
36 36 Homogenization Theorem 1. Since the new velocity b is divergence-free and has zero average, classical periodic homogenization can be applied lim µ ǫ = µ, ǫ 0 ũ ǫ ũ weakly in H 1 0(Ω) where (µ, ũ) is the first eigencouple of div( D ũ) = µũ ũ = 0 on Ω in Ω with D the usual homogenized matrix for D.
37 37 Conclusion Theorem 2. The asymptotic behavior of the first eigencouple of the original problem is λ ǫ = λ(θ ( ) ǫ 2 + µ + o(1), u ǫ (x) e + θ x x ) ǫ φ ũ(x) ǫ where φ (y) = ψ θ (y)e θ y is a Y -periodic function. Remark. The direction of concentration θ is different from the drift b = Y b(y) dy!
38 38 (a) t=0 (b) t=0. (c) t=0.02 (d) t=0.03 (e) t=0.04 (f) t=0.05 Isovalues of u ǫ for various t in the parabolic case (computations by I. Pankratova)
39 39 (g) t=0 (h) t=0.1 (i) t=0.2 (j) t=0.3 (k) t=0.4 (l) t=0.5 Rescaled isovalues of u ǫ for larger times t in the parabolic case (computations by I. Pankratova)
CONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More informationPeriodic homogenization and effective mass theorems for the Schrödinger equation
Periodic homogenization and effective mass theorems for the Schrödinger equation Grégoire Allaire September 5, 2006 Abstract The goal of this course is to give an introduction to periodic homogenization
More informationHOMOGENIZATION OF A CONVECTIVE, CONDUCTIVE AND RADIATIVE HEAT TRANSFER PROBLEM
Homogenization of a heat transfer problem 1 G. Allaire HOMOGENIZATION OF A CONVECTIVE, CONDUCTIVE AND RADIATIVE HEAT TRANSFER PROBLEM Work partially supported by CEA Grégoire ALLAIRE, Ecole Polytechnique
More informationHomogenization of reactive flows in porous media
Homogenization of reactive flows in porous media HUTRIDURGA RAMAIAH Harsha Abstract In this article we study reactive flows through porous media. We use the method of two-scale asymptotic expansions with
More informationA BRIEF INTRODUCTION TO HOMOGENIZATION AND MISCELLANEOUS APPLICATIONS
ESAIM: PROCEEDINGS, September 22, Vol. 37, p. 1-49 E. Cancès and S. Labbé, Editors A BRIEF INTRODUCTION TO HOMOGENIZATION AND MISCELLANEOUS APPLICATIONS Grégoire Allaire 1 Abstract. This paper is a set
More informationHOMOGENIZATION OF ELECTROKINETIC FLOWS IN POROUS MEDIA: THE ROLE OF NON-IDEALITY.
Ion transport in porous media 1 G. Allaire HOMOGENIZATION OF ELECTROKINETIC FLOWS IN POROUS MEDIA: THE ROLE OF NON-IDEALITY. Grégoire ALLAIRE, Robert BRIZZI, CMAP, Ecole Polytechnique Jean-François DUFRECHE,
More informationTWO-SCALE CONVERGENCE ON PERIODIC SURFACES AND APPLICATIONS
TWO-SCALE CONVERGENCE ON PERIODIC SURFACES AND APPLICATIONS Grégoire ALLAIRE Commissariat à l Energie Atomique DRN/DMT/SERMA, C.E. Saclay 91191 Gif sur Yvette, France Laboratoire d Analyse Numérique, Université
More informationHomogenization of micro-resonances and localization of waves.
Homogenization of micro-resonances and localization of waves. Valery Smyshlyaev University College London, UK July 13, 2012 (joint work with Ilia Kamotski UCL, and Shane Cooper Bath/ Cardiff) Valery Smyshlyaev
More informationLINEAR FLOW IN POROUS MEDIA WITH DOUBLE PERIODICITY
PORTUGALIAE MATHEMATICA Vol. 56 Fasc. 2 1999 LINEAR FLOW IN POROUS MEDIA WITH DOUBLE PERIODICITY R. Bunoiu and J. Saint Jean Paulin Abstract: We study the classical steady Stokes equations with homogeneous
More informationTHE STOKES SYSTEM R.E. SHOWALTER
THE STOKES SYSTEM R.E. SHOWALTER Contents 1. Stokes System 1 Stokes System 2 2. The Weak Solution of the Stokes System 3 3. The Strong Solution 4 4. The Normal Trace 6 5. The Mixed Problem 7 6. The Navier-Stokes
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part III Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationHomogenization of reactive flows in porous media and competition between bulk and surface diffusion
Homogenization of reactive flows in porous media and competition between bulk and surface diffusion G. Allaire, H. Hutridurga Abstract In this work, we study the convection and diffusion of a solute in
More informationTWO-SCALE ASYMPTOTIC EXPANSION : BLOCH-FLOQUET THEORY IN PERIODIC STRUCTURES
TWO-SCALE ASYMPTOTIC EXPANSION : BLOCH-FLOQUET THEORY IN PERIODIC STRUCTURES JACK ARBUNICH 1. Setting of Periodic Structures Our aim is to study an application of Bloch-Floquet Theory in the multiscale
More informationColloids transport in porous media: analysis and applications.
Colloids transport in porous media: analysis and applications. Oleh Krehel joint work with Adrian Muntean and Peter Knabner CASA, Department of Mathematics and Computer Science. Eindhoven University of
More informationContinuous dependence estimates for the ergodic problem with an application to homogenization
Continuous dependence estimates for the ergodic problem with an application to homogenization Claudio Marchi Bayreuth, September 12 th, 2013 C. Marchi (Università di Padova) Continuous dependence Bayreuth,
More informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
More informationHT HT HOMOGENIZATION OF A CONDUCTIVE AND RADIATIVE HEAT TRANSFER PROBLEM, SIMULATION WITH CAST3M
Proceedings of HT2005 2005 ASME Summer Heat Transfer Conference July 17-22, 2005, San Francisco, California, USA HT2005-72193 HT2005-72193 HOMOGENIZATION OF A CONDUCTIVE AND RADIATIVE HEAT TRANSFER PROBLEM,
More informationECOLE POLYTECHNIQUE. Two Asymptotic Models for Arrays of Underground Waste Containers CENTRE DE MATHÉMATIQUES APPLIQUÉES UMR CNRS 7641
ECOLE POLYTECHNIQUE CENTRE DE MATHÉMATIQUES APPLIQUÉES UMR CNRS 764 928 PALAISEAU CEDEX FRANCE. Tél: 69 33 46. Fax: 69 33 46 46 http://www.cmap.polytechnique.fr/ Two Asymptotic Models for Arrays of Underground
More informationFree energy estimates for the two-dimensional Keller-Segel model
Free energy estimates for the two-dimensional Keller-Segel model dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine in collaboration with A. Blanchet (CERMICS, ENPC & Ceremade) & B.
More informationDurham Research Online
Durham Research Online Deposited in DRO: 23 April 2018 Version of attached le: Accepted Version Peer-review status of attached le: Peer-reviewed Citation for published item: Cooper, S. (2018 'Quasi-periodic
More informationHomogenization of a Conductive, Convective and Radiative Heat Transfer Problem in a Heterogeneous Domain
Homogenization of a Conductive, Convective and Radiative Heat Transfer Problem in a Heterogeneous Domain Grégoire Allaire, Zakaria Habibi To cite this version: Grégoire Allaire, Zakaria Habibi. Homogenization
More informationHomogenization of the Schrodinger equation with a time oscillating potential
Homogenization of the Schrodinger equation with a time oscillating potential Grégoire Allaire, M. Vanninathan To cite this version: Grégoire Allaire, M. Vanninathan. Homogenization of the Schrodinger equation
More informationLecture 12: Detailed balance and Eigenfunction methods
Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility), 4.2-4.4 (explicit examples of eigenfunction methods) Gardiner
More informationHomogenization and Multiscale Modeling
Ralph E. Showalter http://www.math.oregonstate.edu/people/view/show Department of Mathematics Oregon State University Multiscale Summer School, August, 2008 DOE 98089 Modeling, Analysis, and Simulation
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More informationTransport equations with disparate advection fields. Application to the gyrokinetic models in plasma physics
Transport equations with disparate advection fields. Application to the gyrokinetic models in plasma physics Mihai Bostan (July 22, 2 Abstract The subject matter of this paper concerns the asymptotic regimes
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationUniform estimates for Stokes equations in domains with small holes and applications in homogenization problems
Uniform estimates for Stokes equations in domains with small holes and applications in homogenization problems Yong Lu Abstract We consider the Dirichlet problem for the Stokes equations in a domain with
More informationSecond order corrector in the homogenization of a conductive-radiative heat transfer problem
Second order corrector in the homogenization of a conductive-radiative heat transfer problem Grégoire Allaire, Zakaria Habibi To cite this version: Grégoire Allaire, Zakaria Habibi. Second order corrector
More informationSmilansky-Solomyak model with a δ -interaction
Smilansky-Solomyak model with a δ -interaction Jiří Lipovský University of Hradec Králové, Faculty of Science jiri.lipovsky@uhk.cz joint work with P. Exner Ostrava, April 6, 2018 Jiří Lipovský Smilansky-Solomyak
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 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 information1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω.
1 The Stokes System The motion of a (possibly compressible) homogeneous fluid is described by its density ρ(x, t), pressure p(x, t) and velocity v(x, t). Assume that the fluid is barotropic, i.e., the
More informationApplications of the periodic unfolding method to multi-scale problems
Applications of the periodic unfolding method to multi-scale problems Doina Cioranescu Université Paris VI Santiago de Compostela December 14, 2009 Periodic unfolding method and multi-scale problems 1/56
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
More informationSOME PROBLEMS IN SCALING UP THE SOURCE TERMS IN AN UNDERGROUND WASTE REPOSITORY. A.Bourgeat
SOME PROBLEMS IN SCALING UP THE SOURCE TERMS IN AN UNDERGROUND WASTE REPOSITORY A.Bourgeat Oberwolfach Conference - October 31, November 4; 2005 - Reactive Flow and Transport in Complex Systems Alain Bourgeat
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 informationand finally, any second order divergence form elliptic operator
Supporting Information: Mathematical proofs Preliminaries Let be an arbitrary bounded open set in R n and let L be any elliptic differential operator associated to a symmetric positive bilinear form B
More informationExit times of diffusions with incompressible drifts
Exit times of diffusions with incompressible drifts Andrej Zlatoš University of Chicago Joint work with: Gautam Iyer (Carnegie Mellon University) Alexei Novikov (Pennylvania State University) Lenya Ryzhik
More informationQuantitative Homogenization of Elliptic Operators with Periodic Coefficients
Quantitative Homogenization of Elliptic Operators with Periodic Coefficients Zhongwei Shen Abstract. These lecture notes introduce the quantitative homogenization theory for elliptic partial differential
More informationIV. Transport Phenomena Lecture 18: Forced Convection in Fuel Cells II
IV. Transport Phenomena Lecture 18: Forced Convection in Fuel Cells II MIT Student (and MZB) As discussed in the previous lecture, we are interested in forcing fluid to flow in a fuel cell in order to
More informationHomogenization limit for electrical conduction in biological tissues in the radio-frequency range
Homogenization limit for electrical conduction in biological tissues in the radio-frequency range Micol Amar a,1 Daniele Andreucci a,2 Paolo Bisegna b,2 Roberto Gianni a,2 a Dipartimento di Metodi e Modelli
More informationKPP Pulsating Traveling Fronts within Large Drift
KPP Pulsating Traveling Fronts within Large Drift Mohammad El Smaily Joint work with Stéphane Kirsch University of British olumbia & Pacific Institute for the Mathematical Sciences September 17, 2009 PIMS
More informationTHE L 2 -HODGE THEORY AND REPRESENTATION ON R n
THE L 2 -HODGE THEORY AND REPRESENTATION ON R n BAISHENG YAN Abstract. We present an elementary L 2 -Hodge theory on whole R n based on the minimization principle of the calculus of variations and some
More informationMath The Laplacian. 1 Green s Identities, Fundamental Solution
Math. 209 The Laplacian Green s Identities, Fundamental Solution Let be a bounded open set in R n, n 2, with smooth boundary. The fact that the boundary is smooth means that at each point x the external
More informationThe Factorization Method for Inverse Scattering Problems Part I
The Factorization Method for Inverse Scattering Problems Part I Andreas Kirsch Madrid 2011 Department of Mathematics KIT University of the State of Baden-Württemberg and National Large-scale Research Center
More informationChapter 13. Eddy Diffusivity
Chapter 13 Eddy Diffusivity Glenn introduced the mean field approximation of turbulence in two-layer quasigesotrophic turbulence. In that approximation one must solve the zonally averaged equations for
More informationDRIFT OF SPECTRALLY STABLE SHIFTED STATES ON STAR GRAPHS
DRIFT OF SPECTRALLY STABLE SHIFTED STATES ON STAR GRAPHS ADILBEK KAIRZHAN, DMITRY E. PELINOVSKY, AND ROY H. GOODMAN Abstract. When the coefficients of the cubic terms match the coefficients in the boundary
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationTOPICS IN HOMOGENIZATION OF DIFFERENTIAL EQUATIONS 7
TOPICS IN HOMOGENIZATION OF DIFFERENTIAL EQUATIONS 7. Two-scale asymptotic epansions method We will illustrate this method by applying it to an elliptic boundary value problem for the unknown u : æ R given
More informationMATH 425, HOMEWORK 3 SOLUTIONS
MATH 425, HOMEWORK 3 SOLUTIONS Exercise. (The differentiation property of the heat equation In this exercise, we will use the fact that the derivative of a solution to the heat equation again solves the
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationDiffusion on the half-line. The Dirichlet problem
Diffusion on the half-line The Dirichlet problem Consider the initial boundary value problem (IBVP) on the half line (, ): v t kv xx = v(x, ) = φ(x) v(, t) =. The solution will be obtained by the reflection
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 informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationStability of an abstract wave equation with delay and a Kelvin Voigt damping
Stability of an abstract wave equation with delay and a Kelvin Voigt damping University of Monastir/UPSAY/LMV-UVSQ Joint work with Serge Nicaise and Cristina Pignotti Outline 1 Problem The idea Stability
More informationGLOBAL ATTRACTIVITY IN A CLASS OF NONMONOTONE REACTION-DIFFUSION EQUATIONS WITH TIME DELAY
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 17, Number 1, Spring 2009 GLOBAL ATTRACTIVITY IN A CLASS OF NONMONOTONE REACTION-DIFFUSION EQUATIONS WITH TIME DELAY XIAO-QIANG ZHAO ABSTRACT. The global attractivity
More information1 Introduction We will consider traveling waves for reaction-diusion equations (R-D) u t = nx i;j=1 (a ij (x)u xi ) xj + f(u) uj t=0 = u 0 (x) (1.1) w
Reaction-Diusion Fronts in Periodically Layered Media George Papanicolaou and Xue Xin Courant Institute of Mathematical Sciences 251 Mercer Street, New York, N.Y. 10012 Abstract We compute the eective
More informationProblem set 3: Solutions Math 207B, Winter Suppose that u(x) is a non-zero solution of the eigenvalue problem. (u ) 2 dx, u 2 dx.
Problem set 3: Solutions Math 27B, Winter 216 1. Suppose that u(x) is a non-zero solution of the eigenvalue problem u = λu < x < 1, u() =, u(1) =. Show that λ = (u ) 2 dx u2 dx. Deduce that every eigenvalue
More informationHomogenization of the compressible Navier Stokes equations in domains with very tiny holes
Homogenization of the compressible Navier Stokes equations in domains with very tiny holes Yong Lu Sebastian Schwarzacher Abstract We consider the homogenization problem of the compressible Navier Stokes
More informationON A PRIORI ERROR ANALYSIS OF FULLY DISCRETE HETEROGENEOUS MULTISCALE FEM
MULTISCALE MODEL. SIMUL. Vol. 4, No. 2, pp. 447 459 c 2005 Society for Industrial and Applied Mathematics ON A PRIORI ERROR ANALYSIS OF FULLY DISCRETE HETEROGENEOUS MULTISCALE FEM ASSYR ABDULLE Abstract.
More informationThe KPP minimal speed within large drift in two dimensions
The KPP minimal speed within large drift in two dimensions Mohammad El Smaily Joint work with Stéphane Kirsch University of British Columbia & Pacific Institute for the Mathematical Sciences Banff, March-2010
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
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 informationUniform estimates for Stokes equations in domains with small holes and applications in homogenization problems
Uniform estimates for Stokes equations in domains with small holes and applications in homogenization problems Yong Lu Abstract We consider the Dirichlet problem for the Stokes equations in a domain with
More informationFormal Asymptotic Homogenization
September 21 25, 2015 Formal asymptotic homogenization Here we present a formal asymptotic technique (sometimes called asymptotic homogenization) based on two-spatial scale expansions These expansions
More information4.7 Dispersion in an oscillatory shear flow
Lecture notes in Fluid ynamics.63j/.0j) by Chiang C. Mei, MIT, Spring, 007 4-6dispersion.tex March, 007 [Refs]:. Aris:. Fung, Y. C. Biomechanics 4.7 ispersion in an oscillatory shear flow Relevant to the
More informationADJOINT METHODS FOR OBSTACLE PROBLEMS AND WEAKLY COUPLED SYSTEMS OF PDE. 1. Introduction
ADJOINT METHODS FOR OBSTACLE PROBLEMS AND WEAKLY COPLED SYSTEMS OF PDE F. CAGNETTI, D. GOMES, AND H.V. TRAN Abstract. The adjoint method, recently introduced by Evans, is used to study obstacle problems,
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationKrein-Rutman Theorem and the Principal Eigenvalue
Chapter 1 Krein-Rutman Theorem and the Principal Eigenvalue The Krein-Rutman theorem plays a very important role in nonlinear partial differential equations, as it provides the abstract basis for the proof
More informationReal Variables # 10 : Hilbert Spaces II
randon ehring Real Variables # 0 : Hilbert Spaces II Exercise 20 For any sequence {f n } in H with f n = for all n, there exists f H and a subsequence {f nk } such that for all g H, one has lim (f n k,
More informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More informationBoundary conditions. Diffusion 2: Boundary conditions, long time behavior
Boundary conditions In a domain Ω one has to add boundary conditions to the heat (or diffusion) equation: 1. u(x, t) = φ for x Ω. Temperature given at the boundary. Also density given at the boundary.
More informationControl, Stabilization and Numerics for Partial Differential Equations
Paris-Sud, Orsay, December 06 Control, Stabilization and Numerics for Partial Differential Equations Enrique Zuazua Universidad Autónoma 28049 Madrid, Spain enrique.zuazua@uam.es http://www.uam.es/enrique.zuazua
More informationNonlinear and Nonlocal Degenerate Diffusions on Bounded Domains
Nonlinear and Nonlocal Degenerate Diffusions on Bounded Domains Matteo Bonforte Departamento de Matemáticas, Universidad Autónoma de Madrid, Campus de Cantoblanco 28049 Madrid, Spain matteo.bonforte@uam.es
More informationSTRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD
Shape optimization 1 G. Allaire STRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD G. ALLAIRE CMAP, Ecole Polytechnique with F. JOUVE and A.-M. TOADER 1. Introduction 2. Setting of the problem 3. Shape differentiation
More informationLocal smoothing and Strichartz estimates for manifolds with degenerate hyperbolic trapping
Local smoothing and Strichartz estimates for manifolds with degenerate hyperbolic trapping H. Christianson partly joint work with J. Wunsch (Northwestern) Department of Mathematics University of North
More informationREACTION-DIFFUSION EQUATIONS FOR POPULATION DYNAMICS WITH FORCED SPEED II - CYLINDRICAL-TYPE DOMAINS. Henri Berestycki and Luca Rossi
Manuscript submitted to Website: http://aimsciences.org AIMS Journals Volume 00, Number 0, Xxxx XXXX pp. 000 000 REACTION-DIFFUSION EQUATIONS FOR POPULATION DYNAMICS WITH FORCED SPEED II - CYLINDRICAL-TYPE
More informationExit times of diffusions with incompressible drift
Exit times of diffusions with incompressible drift Gautam Iyer, Carnegie Mellon University gautam@math.cmu.edu Collaborators: Alexei Novikov, Penn. State Lenya Ryzhik, Stanford University Andrej Zlatoš,
More informationHomogenization Theory
Homogenization Theory Sabine Attinger Lecture: Homogenization Tuesday Wednesday Thursday August 15 August 16 August 17 Lecture Block 1 Motivation Basic Ideas Elliptic Equations Calculation of Effective
More informationKramers formula for chemical reactions in the context of Wasserstein gradient flows. Michael Herrmann. Mathematical Institute, University of Oxford
eport no. OxPDE-/8 Kramers formula for chemical reactions in the context of Wasserstein gradient flows by Michael Herrmann Mathematical Institute, University of Oxford & Barbara Niethammer Mathematical
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 informationLECTURE 3 NUMERICAL METHODS OF HOMOGENIZATION
CEA-EDF-INRIA school on homogenization, 13-16 December 20 LECTURE 3 NUMERICAL METHODS OF HOMOGENIZATION Grégoire ALLAIRE Ecole Polytechnique CMAP 91128 Palaiseau, France The goal of this lecture is to
More informationControllability of linear PDEs (I): The wave equation
Controllability of linear PDEs (I): The wave equation M. González-Burgos IMUS, Universidad de Sevilla Doc Course, Course 2, Sevilla, 2018 Contents 1 Introduction. Statement of the problem 2 Distributed
More informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
More informationAnalysis of PDE models of molecular motors. Benoît Perthame based on works with P. E. Souganidis
Analysis of PDE models of molecular motors Benoît Perthame based on works with P. E. Souganidis Asymmetry and motion Except human constructions movement always uses asymmetry Asymmetry and motion See M.
More informationSection 12.6: Non-homogeneous Problems
Section 12.6: Non-homogeneous Problems 1 Introduction Up to this point all the problems we have considered are we what we call homogeneous problems. This means that for an interval < x < l the problems
More informationSome Notes on Elliptic Regularity
Some Notes on Elliptic Regularity Here we consider first the existence of weak solutions to elliptic problems of the form: { Lu = f, u = 0, (1) and then we consider the regularity of such solutions. The
More informationSome recent results on controllability of coupled parabolic systems: Towards a Kalman condition
Some recent results on controllability of coupled parabolic systems: Towards a Kalman condition F. Ammar Khodja Clermont-Ferrand, June 2011 GOAL: 1 Show the important differences between scalar and non
More informationCRYSTAL DISSOLUTION AND PRECIPITATION IN POROUS MEDIA: L 1 -CONTRACTION AND UNIQUENESS
DSCRETE AND CONTNUOUS Website: www.amsciences.org DYNAMCAL SYSTEMS SUPPLEMENT 2007 pp. 1013 1020 CRYSTAL DSSOLUTON AND PRECPTATON N POROUS MEDA: L 1 -CONTRACTON AND UNQUENESS T. L. van Noorden,. S. Pop
More informationMathematical modelling of collective behavior
Mathematical modelling of collective behavior Young-Pil Choi Fakultät für Mathematik Technische Universität München This talk is based on joint works with José A. Carrillo, Maxime Hauray, and Samir Salem
More informationConservation Laws of Surfactant Transport Equations
Conservation Laws of Surfactant Transport Equations Alexei Cheviakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada Winter 2011 CMS Meeting Dec. 10, 2011 A. Cheviakov
More informationON THE SPECTRUM OF THE DIRICHLET LAPLACIAN IN A NARROW STRIP
ON THE SPECTRUM OF THE DIRICHLET LAPLACIAN IN A NARROW STRIP LEONID FRIEDLANDER AND MICHAEL SOLOMYAK Abstract. We consider the Dirichlet Laplacian in a family of bounded domains { a < x < b, 0 < y < h(x)}.
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 informationanalysis for transport equations and applications
Multi-scale analysis for transport equations and applications Mihaï BOSTAN, Aurélie FINOT University of Aix-Marseille, FRANCE mihai.bostan@univ-amu.fr Numerical methods for kinetic equations Strasbourg
More informationOn Schrödinger equations with inverse-square singular potentials
On Schrödinger equations with inverse-square singular potentials Veronica Felli Dipartimento di Statistica University of Milano Bicocca veronica.felli@unimib.it joint work with Elsa M. Marchini and Susanna
More informationOn quasi-richards equation and finite volume approximation of two-phase flow with unlimited air mobility
On quasi-richards equation and finite volume approximation of two-phase flow with unlimited air mobility B. Andreianov 1, R. Eymard 2, M. Ghilani 3,4 and N. Marhraoui 4 1 Université de Franche-Comte Besançon,
More informationEnhanced resolution in structured media
Enhanced resolution in structured media Eric Bonnetier w/ H. Ammari (Ecole Polytechnique), and Yves Capdeboscq (Oxford) Outline : 1. An experiment of super resolution 2. Small volume asymptotics 3. Periodicity
More informationComposite Membranes. 8th August 2006
Composite Membranes Aldo Lopes UFMG, aldoelopes@hotmail.com Russell E. Howes, Brigham Young University rhowes@byu.edu Cynthia Shepherd, California State University - Northridge Cynthia.Shepherd.87@csun.edu
More information