Formal Asymptotic Homogenization

Similar documents
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media

Periodic homogenization and effective mass theorems for the Schrödinger equation

Solving x 2 Dy 2 = N in integers, where D > 0 is not a perfect square. Keith Matthews

Digital Control Systems State Feedback Control

TRANSPORT IN POROUS MEDIA

Qualifying Examination

A BRIEF INTRODUCTION TO HOMOGENIZATION AND MISCELLANEOUS APPLICATIONS

HOMOGENIZATION OF A CONVECTIVE, CONDUCTIVE AND RADIATIVE HEAT TRANSFER PROBLEM

1 Overview. 2 A Characterization of Convex Functions. 2.1 First-order Taylor approximation. AM 221: Advanced Optimization Spring 2016

CONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence

[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,

Homogenization Method and Multiscale Modeling. Adrian Muntean Vladimír Chalupecký

Ordinary Differential Equations II

/ / MET Day 000 NC1^ INRTL MNVR I E E PRE SLEEP K PRE SLEEP R E

Module Two: Differential Calculus(continued) synopsis of results and problems (student copy)

Introduction to gradient descent

ALGEBRAIC GEOMETRY HOMEWORK 3

Computational Optimization. Convexity and Unconstrained Optimization 1/29/08 and 2/1(revised)

Method of Lagrange Multipliers

JUST THE MATHS SLIDES NUMBER 9.6. MATRICES 6 (Eigenvalues and eigenvectors) A.J.Hobson

Additional questions for Chapter 1 - solutions: 1. See Varian, p.20. Notice that Varian set a1=a2. 2. See Varian, p.20.

Stochastic homogenization 1

Upscaling Wave Computations

Extra Problems and Examples

FFT-based Galerkin method for homogenization of periodic media

Second Order Elliptic PDE

Math Exam 2, October 14, 2008

Goals: Second-order Linear Equations Linear Independence of Solutions and the Wronskian Homogeneous DEs with Constant Coefficients

MATH443 PARTIAL DIFFERENTIAL EQUATIONS Second Midterm Exam-Solutions. December 6, 2017, Wednesday 10:40-12:30, SA-Z02

Lecture Notes 6: Dynamic Equations Part B: Second and Higher-Order Linear Difference Equations in One Variable

Homogenization in Elasticity

Introduction to the Calculus of Variations

Second order corrector in the homogenization of a conductive-radiative heat transfer problem

STEP Support Programme. Pure STEP 1 Questions

THE STATISTICAL SECOND-ORDER TWO-SCALE METHOD FOR HEAT TRANSFER PERFORMANCES OF RANDOM POROUS MATERIALS WITH INTERIOR SURFACE RADIATION


Learning Target: I can sketch the graphs of rational functions without a calculator. a. Determine the equation(s) of the asymptotes.

Crew of25 Men Start Monday On Showboat. Many Permanent Improvements To Be Made;Project Under WPA

10. Rank-nullity Definition Let A M m,n (F ). The row space of A is the span of the rows. The column space of A is the span of the columns.

Homogenization Theory

Linear Normed Spaces (cont.) Inner Product Spaces

1.7. Stability and attractors. Consider the autonomous differential equation. (7.1) ẋ = f(x),

MEMORY EFFECT IN HOMOGENIZATION OF VISCOELASTIC KELVIN-VOIGT MODEL WITH TIME DEPENDENT COEFFICIENTS

G H. Extended Unit Tests A L L. Higher Still Advanced Higher Mathematics. (more demanding tests covering all levels) Contents. 3 Extended Unit Tests

Homogenization with stochastic differential equations

Nonlinear equations and optimization

Math Linear Algebra II. 1. Inner Products and Norms

1. Nonlinear Equations. This lecture note excerpted parts from Michael Heath and Max Gunzburger. f(x) = 0

Spectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem

First and Second Order ODEs

98 Cours d analyse September 1, Chapter V. Determination of continuous functions of a single variable that satisfy certain conditions.

Examination paper for TMA4180 Optimization I

Continuous dependence estimates for the ergodic problem with an application to homogenization

ODE classification. February 7, Nasser M. Abbasi. compiled on Wednesday February 07, 2018 at 11:18 PM

Functions & Graphs. Section 1.2

Homogenization of micro-resonances and localization of waves.

An Overly Simplified and Brief Review of Differential Equation Solution Methods. 1. Some Common Exact Solution Methods for Differential Equations

Lecture 7: Advanced Coupling & Mixing Time via Eigenvalues

17.2 Nonhomogeneous Linear Equations. 27 September 2007

Midterm Solution

Lecture 4: Convex Functions, Part I February 1

Abstract Vector Spaces and Concrete Examples

Local properties of plane algebraic curves

Ph.D. Katarína Bellová Page 1 Mathematics 2 (10-PHY-BIPMA2) EXAM - Solutions, 20 July 2017, 10:00 12:00 All answers to be justified.

Cheat Sheet on Linearization

Lecture Notes on PDEs

Further Mathematical Methods (Linear Algebra)

Wave Equation. Finite difference Calculation in Tesseral. Content. 1 Introduction. Finite difference solution of Wave Equation in Tesseral

Functional Analysis, Math 7320 Lecture Notes from August taken by Yaofeng Su

Introduction to Dynamical Systems

Finite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations

Partial Differential Equations for Engineering Math 312, Fall 2012

Simple Examples on Rectangular Domains

STABILITY. Phase portraits and local stability

Orthonormal Systems. Fourier Series

Basic concepts in DT systems. Alexandra Branzan Albu ELEC 310-Spring 2009-Lecture 4 1

Preliminary draft only: please check for final version

The spectral action for Dirac operators with torsion

Core Mathematics C34

A very short introduction to the Finite Element Method

Basic Theory of Linear Differential Equations

DO NOT OPEN THIS QUESTION BOOKLET UNTIL YOU ARE TOLD TO DO SO

Chapter 3 Second Order Linear Equations

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

Math 2a Prac Lectures on Differential Equations

Schrödinger operators exhibiting parameter-dependent spectral transitions

Integral Representation Formula, Boundary Integral Operators and Calderón projection

ASYMPTOTIC HOMOGENISATION IN STRENGTH AND FATIGUE DURABILITY ANALYSIS OF COMPOSITES

Homogenization of a Conductive, Convective and Radiative Heat Transfer Problem in a Heterogeneous Domain

CENTER MANIFOLD AND NORMAL FORM THEORIES

Jose Luis Blázquez Salcedo

Short-time Fourier transform for quaternionic signals

Generalized Gradient Descent Algorithms

0.1 Problems to solve

CS60021: Scalable Data Mining. Large Scale Machine Learning

Homomorphism Testing

Analysis IV : Assignment 3 Solutions John Toth, Winter ,...). In particular for every fixed m N the sequence (u (n)

Homogenization and Multiscale Modeling

Chapter 2 Subspaces of R n and Their Dimensions

A Note on the Central Limit Theorem for a Class of Linear Systems 1

Transcription:

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 allow to move away the oscillations and successfully lead to upscaled model equations and effective coefficients, under the basic assumption of periodicity The introduced expansions are only formal; analysis work is required to justify the validity of the two-scale expansion

Model Problem: Setting of the problem Consider a model problem of diffusion or conductivity in a periodic medium (e.g., an heterogeneous domain obtained by mixing periodically two different phases, one being the matrix and the other the inclusions) Let a domain Ω (a bounded open set in R d, d 1) be periodic with period diam(ω) Denote Y = (0,1) d. For a given source term f(x) L 2 (Ω) scalar function defined in Ω, consider the following BVP: Stationary diffusion in perforated media { (a(x) u (x)) = f(x), x Ω u (x) = 0, x Ω As before, assume a (x) = a ( x ), where a( ) is a Y-periodic function, bounded and s.t. d α, β > 0 : α ζ 2 a ij(y)ζ iζ j β ζ 2, ζ R d, y Y i=1 RK. Here we don t make any assumptions on symmetry of a( ) Function u (x) H 1 (Ω) is the unknown potential (temperature) (1)

Asymptotic Homogenization Method Remark: Well-possedness of (1) follows from assumption f(x) L 2 (Ω) and the Lax-Milgram Lemma:! u (x) H 1 0(Ω) : u H 1 (Ω) C, where C > 0 is independent of Definition (homogenization) Averaging the solution of (1) and finding the effective properties of Ω is called homogenization Ansatz: We seek for the solution of (1) in the form of the series: u = k u k (2) k=0

Asymptotic Homogenization Method (cont.) We search for the behavior of the system as 0 Since coefficient of (1) are Y-periodic in x, it is reasonable to assume that all terms of (2) are periodic functions of x For x Ω, introduce y = x Y, then variables x,y are called x slow micro global y fast macro local Hence, assume u (x) = u 0(x,y)+u 1(x,y)+ 2 u 2(x,y)+..., (3) where u k (x,y) is Y-periodic in the second variable y Now scale separation is exploited: we treat x and y as independent variables Using = x + 1 y (4) i.e. for G (x) := g ( x, ) x one has ( G (x) = xg x, x ) + 1 ( y= x yg x, x ) (5) y= x

Asymptotic Homogenization Method (cont.) Introduce the following notations for operators: A 0 := y [a(y) y] A 1 := y [a(y) x] x [a(y) y] A 2 := x [a(y) x] Then substituting (3) into (1), using (4) and (5), and collecting terms of the same powers of, we obtain (0) A 0u 0 := 0 (1) A 0u 1 := A 1u 0 (6) (2) A 0u 2 := A 2u 0 A 1u 1 +f(x) supplied with periodicity condition of u i, i = 0,1,2, w.r.t. y-variable Problem (6-0) with periodicity condition for u 0 yields u 0(x,y) = u 0(x) (7) Note that y is the variable here, while x plays role of a parameter

Asymptotic Homogenization Method (cont.) Co back to (6-1) with periodicity condition for u 1. With (7) it becomes A 0u 1 = y [a(y) xu 0(x)] (8) The RHS of (8) is a multiple of xu 0(x), hence, (by linearity) one may seek for the solution of (6-1) in the form u 1(x,y) = λ(y) xu 0(x)+û 1(x), (9) where the function λ(y) = (λ 1(y),...,λ d (y)) R d is vector-valued Substituting (9) into (6-1) we obtain equations for λ i(y), i = 1,...,d with periodicity of λ i(y), called cell problems (will be given below). But now let us turn to (6-2) for which we apply the following Lemma (Compatibility condition) Let f(y) L 2 #(Y) then there exists a unique u H 1 #(Y) solving { (A(y) w(y)) = f, y Y iff f(y)dy = 0 Y u is Y periodic

Asymptotic Homogenization Method (cont.) Application of this compatibility condition to the RHS of (6-2) yields { } x a(y)[ yu 1(y)+ xu 0(x)]dy = f(x), x Ω Y u 0(x) = 0, x Ω (10)