Homogenization and Multiscale Modeling
|
|
- Myles Moore
- 6 years ago
- Views:
Transcription
1 Ralph E. Showalter Department of Mathematics Oregon State University Multiscale Summer School, August, 2008 DOE Modeling, Analysis, and Simulation of Multiscale Preferential Flow with Małgorzata Peszyńska
2 Introduction Two-scale approximation The Homogenized problem Origins Homogenization is a collection of methods to approximate a heterogeneous problem by a homogeneous one. Example (Elliptic) A ɛ (u ɛ ) = f : (a ɛ ij (x)uɛ x i ) xj = f, where a ɛ ij (x) = a ij(x/ɛ), and a ij ( ) are 1-periodic inr n. Exact model at microscale replaced by homogenized model a ij (x/ɛ) with constant ã ij
3 Introduction Two-scale approximation The Homogenized problem First Issues Question Can we find a limiting problem, Au = f, whose solution u characterizes the limit, lim ɛ 0 u ɛ = u? Example (Elliptic) A(u) = f : (ã ij u xi ) xj = f, where ã ij are constant. This is the homogenized equation with effective coefficients. Question Is A of the same type as A ɛ? The answer is frequently No!.
4 Introduction Two-scale approximation The Homogenized problem An Elliptic Problem Let Y be the unit cube in R N. Assume a( ) is Y -periodic. We want to approximate the solution to the singular problem u ɛ H0 1 (Ω) : a(x/ɛ) u ɛ (x) ϕ(x) dx = Ω Ω F(x)ϕ(x) dx for all ϕ H 1 0 (Ω). (1)
5 Introduction Two-scale approximation The Homogenized problem The approximation We seek an approximation u ɛ (x) = u(x, x/ɛ) + ɛu(x, x/ɛ) + O(ɛ 2 ), x Ω, in which each u(x, y) and U(x, y) is Y -periodic. The gradient is given (formally) by u ɛ (x) = x u(x, x/ɛ) + 1 ɛ yu(x, x/ɛ) + y U(x, x/ɛ) + O(ɛ).
6 Introduction Two-scale approximation The Homogenized problem The solution u ɛ has a bounded gradient, so y u(x, y) = 0 and u ɛ (x) = u(x) + ɛu(x, x/ɛ) + O(ɛ 2 ) (2a) u ɛ (x) = u(x) + y U(x, x/ɛ) + O(ɛ) (2b) Substitute (2) into (1) with a test function of the same form ϕ(x) + ɛφ(x, x/ɛ) where Φ(x, y) is Y -periodic for each x Ω.
7 Introduction Two-scale approximation The Homogenized problem The variational form Two-scale limits ɛ 0 give the system u H0 1 (Ω), U L 2 (Ω, H#(Y 1 )) : a(y)( u(x) + y U(x, y)) ( ϕ(x) + y Φ(x, y)) dy dx Ω Y = F (x)ϕ(x) dx for all ϕ H0 1 (Ω), Φ L 2 (Ω, H#(Y 1 )). (3) Ω
8 Introduction Two-scale approximation The Homogenized problem The homogenized system Decouple the system: U L 2 (Ω, H#(Y 1 )) : a(y)( y U(x, y) + u(x)) y Φ(x, y) dy dx = 0 Ω Y for all Φ L 2 (Ω, H 1 #(Y )). (4a) u H0 1 (Ω) : a(y)( u(x) + y U(x, y)) dy ϕ(x) dx = F (x)ϕ(x) dx Ω Y Ω for all ϕ H0 1 (Ω). (4b)
9 Introduction Two-scale approximation The Homogenized problem The Cell problem The local problem (4a) is equivalent to requiring U(x, ) H#(Y 1 ) : a(y)( y U(x, y) + u(x)) y Φ(y) dy = 0 for all Φ H#(Y 1 ). Y Define ω i (y) for each 1 i N to be the solution of the cell problem ω i H#(Y 1 ) : a(y)( y ω i (y)+e i ) y Φ(y) dy = 0 for all Φ H#(Y 1 ). Y where e i is the indicated coordinate vector in R N. (5)
10 Introduction Two-scale approximation The Homogenized problem The Homogenized problem By linearity the solution is given by i=n U(x, y) = i u(x)ω i (y) i=1 This is substituted into (4b) to obtain u H0 1 (Ω) : ã ij i u(x) j ϕ(x) dx = Ω Ω F(x)ϕ(x) dx for all ϕ H 1 0 (Ω). (6) where the constant coefficients are given by ã ij = a(y) ( δ ij + j ω i (y) ) dy. (7) Y
11 Introduction Two-scale approximation The Homogenized problem The approximate solution u ɛ (x) = u(x) + ɛ u(x) (ω 1 (x/ɛ), ω 2 (x/ɛ),..., ω N (x/ɛ)) + O(ɛ 2 ), u ɛ (x) = u(x) + u(x) ( ω 1 (x/ɛ), ω 2 (x/ɛ),..., ω N (x/ɛ)) + O(ɛ) of the singular elliptic problem.
12 Introduction Two-scale approximation The Homogenized problem Homogenization Alain Bensoussan, Jacques-Louis Lions, and George Papanicolaou. Asymptotic analysis for periodic structures, volume 5 of Studies in Mathematics and its Applications. North-Holland Publishing Co., Amsterdam, Enrique Sánchez-Palencia. Nonhomogeneous media and vibration theory, volume 127 of Lecture Notes in Physics. Springer-Verlag, Berlin, V. V. Jikov, S. M. Kozlov, and O. A. Oleinik. Homogenization of differential operators and integral functionals. Springer-Verlag, Berlin, 1994.
13 ... the model u(x) = K(x) p(x), u(x) = 0 c(x, t) φ + (u(x)c(x, t) D(u(x)) c(x, t)) t = 0 K fast, K slow u fast, u slow D fast, D slow
14 The Parabolic Equation The Exact Model... with fine-scale coefficients φ(x) c t D(x) c = 0, x Ω, with D(x) = D slow on Ω slow and = D fast on Ω fast.... or in transmission form φ α c α t D α c α = 0, x Ω α, α = fast, slow with the interface conditions on Ω slow Ω fast : c fast = c slow, D fast c fast ν = D slow c slow ν
15 Homogenization The coefficients are D = [D fast, D slow ] : φ c t D c = 0 The fine-scale geometry is eliminated replaced by the constant effective coefficient D. The upscaled limit is of the same type... a single equation. The fast and slow regions are coupled by gradients of the solution (flux). Accurate only in this low contrast case.
16 ... the algorithm Exact model at microscale Homogenization replaced by homogenized model D(x) = D slow, D fast with constant D Compute homogenized coefficient D = D(D fast, D slow ) D jk = 1 D jk (y)(δ jk + k ω j (y))da Ω 0 Ω 0 { D(y) ωj (y) = (D(y)e j ), ω j (y) is Ω 0 periodic.
17 ... the algorithm Exact model at microscale Homogenization replaced by homogenized model D(x) = D slow, D fast with constant D Compute homogenized coefficient D = D(D fast, D slow ) D jk = 1 D jk (y)(δ jk + k ω j (y))da Ω 0 Ω 0 { D(y) ωj (y) = (D(y)e j ), ω j (y) is Ω 0 periodic. It works well for problems with low contrast D fast /D slow
18 The Obstacle Problem... a special case The coefficients are D = [D fast, 0] : φ 0 c0 t D 0 c 0 = 0 The slow region Ω slow is impervious.
19 Classical homogenization is not adequate for time-dependent problems with large contrast D fast /D slow Low contrast: classical homogenization High contrast: double porosity models u ɛ 0 + u ɛ 0 C local averages u ɛ 0 +ɛ u ɛ C local averages versus special averages
20 The Double-porosity model The coefficients are D = [D fast, ɛ 2 0 D slow] : φ 0 c t + i χ i q i (t) D 0 c = 0, q i (t) = 1 D slow c i νds, ˆΩ i Γ i c i φ slow D slow c i = 0, c i Γi = 1 c(x)da. t ˆΩ i ˆΩ i The upscaled model is a system, highly parallel. The coupling from the cell to global equation is via gradients. The coupling from the global equation to the cell is via the values. The contribution of the cells is a secondary storage. Accurate only in this high contrast case. It will not couple any advective effects to the local cell.
21 Double-porosity Micro-structure Model Exact model replaced by two layers + Global equation, x Ω φ 0 c t + χ i q i (t) D 0 c = 0 i q i (t) = Π 0,i (D slow c i (, t) ν) Cell problem at each x Ω i c i φ slow t D slow c i = 0 c i Γi = Π 0,i ( c)(t)
22 Notes Homogenization The coefficients in the global equation are precisely those of the obstacle problem. The cell input to the global equation is q i (t) = 1 φ slow c i (x, t) dx ˆΩ i t Ω i... the rate of secondary storage. The spatially-constant input to the cell problem will not transmit any advective transport. We shall more tightly couple the cell by replacing the constant coupling with an affine coupling to the surrounding fast medium. This will provide a gradient coupling.
23 Affine approximations Π 1 AVERAGE: Π 0 f = 1 Ω 0 Ω 0 f (x)da; assume here Ω 0 = 1. Denote x C - center of mass of Ω 0. General affine approximation f (x) Π 1 f := m + n x, x Ω 0 Choice of m, n Taylor (f C 1 (Ω 0 )) about midpoint f (x) f (x C ) + f (x C )(x x C ) L 2 (Ω 0 )-projection onto affines that is: (f, v) Ω0 = (m + n x, v) Ω0, affine v H 1 (Ω 0 ) projection: f (x) Π 1 f := Π 0 f + Π 0 f (x x C )
24 Summary: the affine coupling Π i and its dual Π i Affine coupling Π i : H 1 (Ω) H 1 (ˆΩ i ) Π i (w)(x) Π 0 w + Π 0 ( w) (x x C ) Its dual Π i : H 1 (ˆΩ i ) H 1 (Ω) pointwise Π i (q)(x) = χ i(x)m 0 i (q) χ i (x)m 1 i (q)
25 Double-porosity Model with secondary flux φ 0 c t + i χ i q i (t) D 0 c = 0, φ slow c i t D slow c i = 0, c i Γi = 1 ˆΩ i cda + 1 c da (x x C ). ˆΩ i ˆΩ i ˆΩ i q i (t) = 1 D slow c i νds χ i(x) D slow c i ν (x x C )ds ˆΩ i Γ i ˆΩ i Γ i This is t cells. (secondary storage) (secondary flux) from the local
26 Computational experiments at microscale GOAL: reproduce qualitatively experimental results
27 Porous Media Homogenization Ulrich Hornung, editor. Homogenization and Porous Media, volume 6 of Interdisciplinary Applied Mathematics. Springer-Verlag, New York, Peszynska, M.; Showalter, R. E. Multiscale Elliptic-Parabolic Systems for Flow and Transport, Electron. J. Differential Equations 2007, No. 147, 30 pp. (electronic). Brendan Zinn, Lucy C. Meigs, Charles F. Harvey, Roy Haggerty, Williams J. Peplinski, and Claudius Freiherr von Schwerin. Experimental visualization of solute transport and mass transfer processes in two-dimensional conductivity fields with connected regions of high conductivity. Environ Sci. Technol., 38: , 2004.
Nonlocal models of transport in multiscale porous media:something old and som. something borrowed...
Nonlocal models of transport in multiscale porous media: something old and something new, something borrowed... Department of Mathematics Oregon State University Outline 1 Flow and transport in subsurface,
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 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 informationTRANSPORT IN POROUS MEDIA
1 TRANSPORT IN POROUS MEDIA G. ALLAIRE CMAP, Ecole Polytechnique 1. Introduction 2. Main result in an unbounded domain 3. Asymptotic expansions with drift 4. Two-scale convergence with drift 5. The case
More 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 informationUpscaling Wave Computations
Upscaling Wave Computations Xin Wang 2010 TRIP Annual Meeting 2 Outline 1 Motivation of Numerical Upscaling 2 Overview of Upscaling Methods 3 Future Plan 3 Wave Equations scalar variable density acoustic
More informationDOUBLE-DIFFUSION MODELS FROM A HIGHLY-HETEROGENEOUS MEDIUM
DOUBLE-DIFFUSION MODELS FROM A HIGHL-HETEROGENEOUS MEDIUM R.E. SHOWALTER AND D.B. VISARRAGA Abstract. A distributed microstructure model is obtained by homogenization from an exact micro-model with continuous
More informationFinite Volume Discretization of a Class of Parabolic Equations with Discontinuous and Oscillating Coefficients
Applied Mathematical Sciences, Vol, 8, no 3, - 35 HIKARI Ltd, wwwm-hikaricom https://doiorg/988/ams876 Finite Volume Discretization of a Class of Parabolic Equations with Discontinuous and Oscillating
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part IV Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationSchur Complement Technique for Advection-Diffusion Equation using Matching Structured Finite Volumes
Advances in Dynamical Systems and Applications ISSN 973-5321, Volume 8, Number 1, pp. 51 62 (213) http://campus.mst.edu/adsa Schur Complement Technique for Advection-Diffusion Equation using Matching Structured
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 informationThe Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition
The Dirichlet boundary problems for second order parabolic operators satisfying a Martin Dindos Sukjung Hwang University of Edinburgh Satellite Conference in Harmonic Analysis Chosun University, Gwangju,
More informationHomogenization of the Transmission Eigenvalue Problem for a Periodic Media
Homogenization of the Transmission Eigenvalue Problem for a Periodic Media Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work
More informationHeat Transfer in a Medium in Which Many Small Particles Are Embedded
Math. Model. Nat. Phenom. Vol. 8, No., 23, pp. 93 99 DOI:.5/mmnp/2384 Heat Transfer in a Medium in Which Many Small Particles Are Embedded A. G. Ramm Department of Mathematics Kansas State University,
More informationMicropolar fluid flow with rapidly variable initial conditions
Mathematical Communications 3998), 97 24 97 Micropolar fluid flow with rapidly variable initial conditions Nermina Mujaković Abstract. In this paper we consider nonstationary D-flow of a micropolar viscous
More informationNUMERICAL JUSTIFICATION OF ASYMPTOTIC EDDY CURRENTS MODEL FOR HETEROGENEOUS MATERIALS
Proceedings of ALGORITMY 2009 pp. 212 218 NUMERICAL JUSTIFICATION OF ASYMPTOTIC EDDY CURRENTS MODEL FOR HETEROGENEOUS MATERIALS V. MELICHER AND J.BUŠA JR. Abstract. We study electromagnetic properties
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 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 informationLECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI
LECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI Georgian Technical University Tbilisi, GEORGIA 0-0 1. Formulation of the corresponding
More informationHeterogeneous media: Case studies in non-periodic averaging and perspectives in designing packaging materials
Heterogeneous media: Case studies in non-periodic averaging and perspectives in designing packaging materials Department of Mathematics and Computer Science, Karlstad University, Sweden Karlstad, January
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 informationHETEROGENEOUS MULTISCALE METHOD IN EDDY CURRENTS MODELING
Proceedings of ALGORITMY 2009 pp. 219 225 HETEROGENEOUS MULTISCALE METHOD IN EDDY CURRENTS MODELING JÁN BUŠA, JR. AND VALDEMAR MELICHER Abstract. The induction of eddy currents in a conductive piece is
More informationOperator Upscaling for the Wave Equation
Operator Upscaling for the Wave Equation Tetyana Vdovina Susan E. Minkoff UMBC), Oksana Korostyshevskaya Department of Computational and Applied Mathematics Rice University, Houston TX vdovina@caam.rice.edu
More informationPeaceman and Thiem well models or how to remove a logarithmic singularity fr. your numerical solution
Peaceman and Thiem well models or how to remove a logarithmic singularity from your numerical solution Department of Mathematics Oregon State University AMC, 3/2/2007 Current work supported by NSF-0511190
More informationApplications of Mathematics
Applications of Mathematics Jeanette Silfver On general two-scale convergence and its application to the characterization of G-limits Applications of Mathematics, Vol. 52 (2007, No. 4, 285--302 Persistent
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 informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationMath 240 Calculus III
Calculus III Summer 2015, Session II Monday, August 3, 2015 Agenda 1. 2. Introduction The reduction of technique, which applies to second- linear differential equations, allows us to go beyond equations
More informationThe Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition
The Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition Sukjung Hwang CMAC, Yonsei University Collaboration with M. Dindos and M. Mitrea The 1st Meeting of
More informationSOLUTION OF POISSON S EQUATION. Contents
SOLUTION OF POISSON S EQUATION CRISTIAN E. GUTIÉRREZ OCTOBER 5, 2013 Contents 1. Differentiation under the integral sign 1 2. The Newtonian potential is C 1 2 3. The Newtonian potential from the 3rd Green
More informationStochastic Volatility and Correction to the Heat Equation
Stochastic Volatility and Correction to the Heat Equation Jean-Pierre Fouque, George Papanicolaou and Ronnie Sircar Abstract. From a probabilist s point of view the Twentieth Century has been a century
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 informationMultiscale Diffusion Modeling in Charged and Crowded Biological Environments
Multiscale Diffusion Modeling in Charged and Crowded Biological Environments Andrew Gillette Department of Mathematics University of Arizona joint work with Pete Kekenes-Huskey (U. Kentucky) and J. Andrew
More informationROLE OF PORE-SCALE HETEROGENEITY ON REACTIVE FLOWS IN POROUS MATERIALS: VALIDITY OF THE CONTINUUM REPRESENTATION OF REACTIVE TRANSPORT
ROLE OF PORE-SCALE HETEROGENEITY ON REACTIVE FLOWS IN POROUS MATERIALS: VALIDITY OF THE CONTINUUM REPRESENTATION OF REACTIVE TRANSPORT PETER C. LICHTNER 1, QINJUN KANG 1 1 Los Alamos National Laboratory,
More informationModeling of two-phase flow in fractured porous media on unstructured non-uniform coarse grids
Modeling of two-phase flow in fractured porous media on unstructured non-uniform coarse grids Jørg Espen Aarnes and Vera Louise Hauge SINTEF ICT, Deptartment of Applied Mathematics Applied Mathematics
More informationHomogenization theory and the assessment of extreme field values in composites with random microstructure Robert Lipton 1
SIAM J. Applied Mathematics, 65, 2004, pp. 475-493 Homogenization theory and the assessment of extreme field values in composites with random microstructure Robert Lipton Abstract. Suitable macroscopic
More informationHomogenization with stochastic differential equations
Homogenization with stochastic differential equations Scott Hottovy shottovy@math.arizona.edu University of Arizona Program in Applied Mathematics October 12, 2011 Modeling with SDE Use SDE to model system
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 informationNOTES ON SCHAUDER ESTIMATES. r 2 x y 2
NOTES ON SCHAUDER ESTIMATES CRISTIAN E GUTIÉRREZ JULY 26, 2005 Lemma 1 If u f in B r y), then ux) u + r2 x y 2 B r y) B r y) f, x B r y) Proof Let gx) = ux) Br y) u r2 x y 2 Br y) f We have g = u + Br
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 informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY SPRING SEMESTER 207 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationMultiscale method and pseudospectral simulations for linear viscoelastic incompressible flows
Interaction and Multiscale Mechanics, Vol. 5, No. 1 (2012) 27-40 27 Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows Ling Zhang and Jie Ouyang* Department of
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part II Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationCONVERGENCE 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 informationWEYL S LEMMA, ONE OF MANY. Daniel W. Stroock
WEYL S LEMMA, ONE OF MANY Daniel W Stroock Abstract This note is a brief, and somewhat biased, account of the evolution of what people working in PDE s call Weyl s Lemma about the regularity of solutions
More informationUniversity of Minnesota, Minneapolis
Coarse Graining of Electric Field Interactions with Materials Prashant Kumar Jha Ph.D., CMU prashant.j16o@gmail.com Adviser Dr. Kaushik Dayal Funded by Army Research Office Research Talk University of
More informationEECS 598: Statistical Learning Theory, Winter 2014 Topic 11. Kernels
EECS 598: Statistical Learning Theory, Winter 2014 Topic 11 Kernels Lecturer: Clayton Scott Scribe: Jun Guo, Soumik Chatterjee Disclaimer: These notes have not been subjected to the usual scrutiny reserved
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 informationAsymptotic Behavior of Waves in a Nonuniform Medium
Available at http://pvamuedu/aam Appl Appl Math ISSN: 1932-9466 Vol 12, Issue 1 June 217, pp 217 229 Applications Applied Mathematics: An International Journal AAM Asymptotic Behavior of Waves in a Nonuniform
More informationIntroduction to PDEs and Numerical Methods Lecture 1: Introduction
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Lecture 1: Introduction Dr. Noemi Friedman, 28.10.2015. Basic information on the course Course
More 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 information17.8 Nonhomogeneous Linear Equations We now consider the problem of solving the nonhomogeneous second-order differential
ADAMS: Calculus: a Complete Course, 4th Edition. Chapter 17 page 1016 colour black August 15, 2002 1016 CHAPTER 17 Ordinary Differential Equations 17.8 Nonhomogeneous Linear Equations We now consider the
More informationApplications of Mathematics
Applications of Mathematics Liselott Flodén; Marianne Olsson Homogenization of some parabolic operators with several time scales Applications of Mathematics, Vol. 52 2007, No. 5, 431--446 Persistent URL:
More informationInverse Scattering Theory: Transmission Eigenvalues and Non-destructive Testing
Inverse Scattering Theory: Transmission Eigenvalues and Non-destructive Testing Isaac Harris Texas A & M University College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with: F. Cakoni, H.
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 informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
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 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 informationOn some weighted fractional porous media equations
On some weighted fractional porous media equations Gabriele Grillo Politecnico di Milano September 16 th, 2015 Anacapri Joint works with M. Muratori and F. Punzo Gabriele Grillo Weighted Fractional PME
More informationMultivariable Calculus
2 Multivariable Calculus 2.1 Limits and Continuity Problem 2.1.1 (Fa94) Let the function f : R n R n satisfy the following two conditions: (i) f (K ) is compact whenever K is a compact subset of R n. (ii)
More informationMath 220A - Fall 2002 Homework 5 Solutions
Math 0A - Fall 00 Homework 5 Solutions. Consider the initial-value problem for the hyperbolic equation u tt + u xt 0u xx 0 < x 0 u t (x, 0) ψ(x). Use energy methods to show that the domain of dependence
More informationDetermination of thin elastic inclusions from boundary measurements.
Determination of thin elastic inclusions from boundary measurements. Elena Beretta in collaboration with E. Francini, S. Vessella, E. Kim and J. Lee September 7, 2010 E. Beretta (Università di Roma La
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 informationTHE STATISTICAL SECOND-ORDER TWO-SCALE METHOD FOR HEAT TRANSFER PERFORMANCES OF RANDOM POROUS MATERIALS WITH INTERIOR SURFACE RADIATION
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING, SERIES B Volume 4, Number 2, Pages 151 166 c 2013 Institute for Scientific Computing and Information THE STATISTICAL SECOND-ORDER TWO-SCALE METHOD
More informationRandom Homogenization and Convergence to Integrals with respect to the Rosenblatt Process
andom Homogenization and Convergence to Integrals with respect to the osenblatt Process Yu Gu Guillaume Bal September 7, Abstract This paper concerns the random fluctuation theory of a one dimensional
More informationAN INVERSE PROBLEM FOR THE WAVE EQUATION WITH A TIME DEPENDENT COEFFICIENT
AN INVERSE PROBLEM FOR THE WAVE EQUATION WITH A TIME DEPENDENT COEFFICIENT Rakesh Department of Mathematics University of Delaware Newark, DE 19716 A.G.Ramm Department of Mathematics Kansas State University
More informationTWO-SCALE CONVERGENCE OF A MODEL FOR FLOW IN A PARTIALLY FISSURED MEDIUM
Electronic Journal of Differential Equations, Vol. 19991999), No., pp. 1. ISSN: 17-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu ejde.math.unt.edu login: ftp) TWO-SCALE
More informationABOUT SOME TYPES OF BOUNDARY VALUE PROBLEMS WITH INTERFACES
13 Kragujevac J. Math. 3 27) 13 26. ABOUT SOME TYPES OF BOUNDARY VALUE PROBLEMS WITH INTERFACES Boško S. Jovanović University of Belgrade, Faculty of Mathematics, Studentski trg 16, 11 Belgrade, Serbia
More informationPRACTICE PROBLEMS. Please let me know if you find any mistakes in the text so that i can fix them. 1. Mixed partial derivatives.
PRACTICE PROBLEMS Please let me know if you find any mistakes in the text so that i can fix them. 1.1. Let Show that f is C 1 and yet How is that possible? 1. Mixed partial derivatives f(x, y) = {xy x
More informationON THE ASYMPTOTIC BEHAVIOR OF ELLIPTIC PROBLEMS IN PERIODICALLY PERFORATED DOMAINS WITH MIXED-TYPE BOUNDARY CONDITIONS
Bulletin of the Transilvania University of Braşov Series III: Mathematics, Informatics, Physics, Vol 5(54) 01, Special Issue: Proceedings of the Seventh Congress of Romanian Mathematicians, 73-8, published
More informationIntegral Equations in Electromagnetics
Integral Equations in Electromagnetics Massachusetts Institute of Technology 6.635 lecture notes Most integral equations do not have a closed form solution. However, they can often be discretized and solved
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationPeriodic homogenization for inner boundary conditions with equi-valued surfaces: the unfolding approach
Periodic homogenization for inner boundary conditions with equi-valued surfaces: the unfolding approach Alain Damlamian Université Paris-Est, Laboratoire d Analyse et de Mathématiques Appliquées, CNRS
More informationFrequency functions, monotonicity formulas, and the thin obstacle problem
Frequency functions, monotonicity formulas, and the thin obstacle problem IMA - University of Minnesota March 4, 2013 Thank you for the invitation! In this talk we will present an overview of the parabolic
More informationLecture No 2 Degenerate Diffusion Free boundary problems
Lecture No 2 Degenerate Diffusion Free boundary problems Columbia University IAS summer program June, 2009 Outline We will discuss non-linear parabolic equations of slow diffusion. Our model is the porous
More informationarxiv: v1 [math.ap] 18 Jan 2019
manuscripta mathematica manuscript No. (will be inserted by the editor) Yongpan Huang Dongsheng Li Kai Zhang Pointwise Boundary Differentiability of Solutions of Elliptic Equations Received: date / Revised
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 225 Estimates of the second-order derivatives for solutions to the two-phase parabolic problem
More informationLecture Notes #10. The "paradox" of finite strain and preferred orientation pure vs. simple shear
12.520 Lecture Notes #10 Finite Strain The "paradox" of finite strain and preferred orientation pure vs. simple shear Suppose: 1) Anisotropy develops as mineral grains grow such that they are preferentially
More informationDouble-diffusion models from a highly-heterogeneous medium
J. Math. Anal. Appl. 295 24) 191 21 www.elsevier.com/locate/jmaa Double-diffusion models from a highly-heterogeneous medium R.E. Showalter a, and D.B. Visarraga b a Department of Mathematics, Oregon State
More informationAM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs
AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference
More informationDistributed Microstructure Models of Porous Media
Distributed Microstructure Models of Porous Media Ralph E. Showalter Abstract. Laminar flow through through fissured or otherwise highly inhomogeneous media leads to very singular initial-boundary-value
More information2. Second-order Linear Ordinary Differential Equations
Advanced Engineering Mathematics 2. Second-order Linear ODEs 1 2. Second-order Linear Ordinary Differential Equations 2.1 Homogeneous linear ODEs 2.2 Homogeneous linear ODEs with constant coefficients
More informationEfficient Methods for Homogenization of Random Heterogeneous Materials
Northwestern Polytechnical University Efficient Methods for Homogenization of Random Heterogeneous Materials Yufeng Nie Department of Applied Mathematics Northwestern Polytechnical University yfnie@nwpu.edu.cn
More informationEFFECTIVE CHARACTERISTICS OF POROUS MEDIA AS A FUNCTION OF POROSITY LEVEL
AMS Subject Classification Index: 74Q20 EFFECTIVE CHARACTERISTICS OF POROUS MEDIA AS A FUNCTION OF POROSITY LEVEL Irini DJERAN-MAIGRE and Sergey V.KUZNETSOV INSA de Lyon, URGC 34 Avenue des Arts, 69621
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 informationCONVERGENCE OF A MULTISCALE FINITE ELEMENT METHOD FOR ELLIPTIC PROBLEMS WITH RAPIDLY OSCILLATING COEFFICIENTS
MATHEMATICS OF COMPUTATION Volume 68, Number 227, Pages 913 943 S 0025-5718(99)01077-7 Article electronically published on March 3, 1999 CONVERGENCE OF A MULTISCALE FINITE ELEMENT METHOD FOR ELLIPTIC PROBLEMS
More informationA Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements
W I S S E N T E C H N I K L E I D E N S C H A F T A Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements Matthias Gsell and Olaf Steinbach Institute of Computational Mathematics
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 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 informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationON QUALITATIVE PROPERTIES OF SOLUTIONS OF QUASILINEAR ELLIPTIC EQUATIONS WITH STRONG DEPENDENCE ON THE GRADIENT
GLASNIK MATEMATIČKI Vol. 49(69)(2014), 369 375 ON QUALITATIVE PROPERTIES OF SOLUTIONS OF QUASILINEAR ELLIPTIC EQUATIONS WITH STRONG DEPENDENCE ON THE GRADIENT Jadranka Kraljević University of Zagreb, Croatia
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 informationT. Arbogast, Zhen Tao, and Hailong Xiao, Multiscale mortar mixed methods for heterogeneous elliptic problems, in Recent Advances in Scientific
T. Arbogast, Zhen Tao, and Hailong Xiao, Multiscale mortar mixed methods for heterogeneous elliptic problems, in Recent Advances in Scientific Computing and Applications, Jichun Li, Hongtao Yang, and Eric
More informationDynamical systems with Gaussian and Levy noise: analytical and stochastic approaches
Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches Noise is often considered as some disturbing component of the system. In particular physical situations, noise becomes
More informationHomogeniza*ons in Perforated Domain. Ki Ahm Lee Seoul Na*onal University
Homogeniza*ons in Perforated Domain Ki Ahm Lee Seoul Na*onal University Outline 1. Perforated Domain 2. Neumann Problems (joint work with Minha Yoo; interes*ng discussion with Li Ming Yeh) 3. Dirichlet
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
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 informationAhmed Mohammed. Harnack Inequality for Non-divergence Structure Semi-linear Elliptic Equations
Harnack Inequality for Non-divergence Structure Semi-linear Elliptic Equations International Conference on PDE, Complex Analysis, and Related Topics Miami, Florida January 4-7, 2016 An Outline 1 The Krylov-Safonov
More informationWentzell Boundary Conditions in the Nonsymmetric Case
Math. Model. Nat. Phenom. Vol. 3, No. 1, 2008, pp. 143-147 Wentzell Boundary Conditions in the Nonsymmetric Case A. Favini a1, G. R. Goldstein b, J. A. Goldstein b and S. Romanelli c a Dipartimento di
More informationSpace-Time Nonconforming Optimized Schwarz Waveform Relaxation for Heterogeneous Problems and General Geometries
Space-Time Nonconforming Optimized Schwarz Waveform Relaxation for Heterogeneous Problems and General Geometries Laurence Halpern, Caroline Japhet, and Jérémie Szeftel 3 LAGA, Université Paris XIII, Villetaneuse
More information