This work could not have been accomplished without the support of many persons. First of all, I would like to express my deep gratitude to Professor… [610112]
UNIVERSITY OF BUCHAREST
HABILITATION THESIS
HOMOGENIZATION RESULTS FOR
HETEROGENEOUS MEDIA
CLAUDIA TIMOFTE
Specialization: Mathematics
Bucharest, 2016
.
Acknowledgments
This work could not have been accomplished without the support of many persons. First
of all, I would like to express my deep gratitude to Professor Horia I. Ene from the Institute
of Mathematics "Simion Stoilow" of the Romanian Academy for his continuous support and
professional guidance. His mentorship was very important for the evolution of my academic
and scientic career.
After completing my Ph.D., I beneted, between 2000 and 2003, from four post-doctoral
fellowships at Complutense University of Madrid (Spain), University of Pisa (Italy), and
Center of Mathematical Modelling, University of Chile, Santiago de Chile (Chile). I had the
chance to work, in stimulating environments, under the supervision of top specialists in ap-
plied mathematics: Professor Enrique Zuazua, Professor Giuseppe Buttazzo, and Professor
Carlos Conca. I was really impressed by their remarkable ability to connect different elds of
research and I want to express my deep gratitude to all of them, for their support and generos-
ity, and for the willingness to share their extraordinary knowledge with me. Also, the research
visits performed at Complutense University of Madrid (Spain), Friedrich-Alexander Univer-
sity Erlangen-N urnberg (Germany), Eindhoven University of Technology (Holland), Taras
Shevchenko National University of Kiev (Ukraine), University of Pavia (Italy), Universit e de
Lorraine, Metz (France), or University of Cantabria, Santander (Spain) gave me the chance of
fruitful interactions with well-known specialists in the eld of homogenization theory: Profes-
sor Jes us Ildefonso D az, Professor Iuliu Sorin Pop, Professor Maria Radu-Neuss, Dr. Renata
Bunoiu, Professor Gennady Sandrakov, Professor Giuseppe Savar e, David G omez-Castro,
Professor Mar a Eugenia P erez, Dr. Delna G omez. I am grateful to all of them, for their
warm hospitality, kindness, and for sharing with me their love for mathematics.
The work presented here represents a collective effort, the fruit of many encounters I
had over the years with many persons and I am fully conscious about their importance at
many steps in my career. It is impossible to me to thank now all the people that I met in
this scientic journey. Therefore, I shall mention here only the co-authors of my papers on
which this thesis is based on: Professor Carlos Conca, Professor Fran cois Murat, Professor
Jes us Ildefonso D az, Professor Amable Li~ n an, Professor Horia I. Ene, Dr. Anca C ap at ^ n a,
Dr. Iulian T ent ea, Dr. Renata Bunoiu. Working together was very important for my
development as a mathematician.
During the last years, I have beneted a lot from inspiring discussions with my colleagues
from the Institute of Mathematics "Simion Stoilow" of the Romanian Academy. I am very
3
4 Acknowledgments
indebted to them for their valuable and constructive suggestions. I would also like to thank to
my colleagues from the Faculty of Physics of the University of Bucharest, who accompanied
me in this transdisciplinary journey, for their support and for the emulating atmosphere they
have always created in our faculty.
There are many other people who helped me at one stage or another in my work and are
not mentioned here. I express my deep gratitude to all of them.
Last, but not least, I am grateful to my family and to my friends for their unwavering
support and understanding.
Bucharest, October 2016 Claudia Timofte
Preface
The aim of this manuscript, prepared to defend my Habilitation thesis , is to give an overview
of my research activity in the eld of homogenization theory, which represents the core of my
scientic work done during the last fteen years.
The thesis, written in English, starts by a short summary in Romanian and a brief
overview of the eld of homogenization and then summarizes, with less details in some proofs
and with some additional hindsights, some of my research works in this eld, performed
after completing my Ph.D. studies. The thesis relies on some of my original contributions
to the applications of the homogenization theory, contained in twenty-ve articles already
published or submitted for publication in international journals. Many of the results in
the publications I selected to support my application are closely related to or motivated by
practical applications to real-life problems.
The results included in this thesis have been obtained alone or in collaboration with
several academic and research institutions from Romania or from abroad. I am grateful to all
my co-authors for their important contribution, for useful advices and friendly discussions.
5
Contents
Preface 5
Rezumat 9
Abstract 11
I Main Scientic Achievements 13
1 Introduction 15
2 Homogenization of reactive
ows in porous media 19
2.1 Upscaling in stationary reactive
ows in porous media . . . . . . . . . . . . . 20
2.1.1 The model problem . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.1.2 The case of a smooth function g. The macroscopic model . . . . . . . 22
2.1.3 The case of a non-smooth function g. The macroscopic model . . . . . 29
2.1.4 Chemical reactions inside the grains of a porous medium . . . . . . . . 34
2.2 Nonlinear adsorption of chemicals in porous media . . . . . . . . . . . . . . . 37
2.2.1 The microscopic model and its weak solvability . . . . . . . . . . . . . 37
2.2.2 The main result . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
2.2.3 The case of a non-smooth boundary condition . . . . . . . . . . . . . . 42
2.2.4 Laplace-Beltrami model with oscillating coefficients . . . . . . . . . . . 43
3 Homogenization results for unilateral problems 49
3.1 Homogenization results for Signorini's type problems . . . . . . . . . . . . . . 50
3.1.1 The microscopic problem . . . . . . . . . . . . . . . . . . . . . . . . . 50
3.1.2 The macroscopic models . . . . . . . . . . . . . . . . . . . . . . . . . . 52
3.2 Homogenization results for elliptic problems in perforated domains with mixed-
type boundary conditions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
3.2.1 Setting of the microscopic problem . . . . . . . . . . . . . . . . . . . . 55
3.2.2 The limit problem . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
7
8 Contents
4 Mathematical models in biology 65
4.1 Homogenization results for ionic transport phenomena in periodic charged media 65
4.1.1 Setting of the problem . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
4.1.2 The homogenized problem . . . . . . . . . . . . . . . . . . . . . . . . . 70
4.2 Multiscale analysis of a carcinogenesis model . . . . . . . . . . . . . . . . . . 73
4.2.1 The microscopic problem . . . . . . . . . . . . . . . . . . . . . . . . . 73
4.2.2 The macroscopic model . . . . . . . . . . . . . . . . . . . . . . . . . . 78
4.2.3 A nonlinear carcinogenesis model involving free receptors . . . . . . . 84
4.3 Homogenization results for the calcium dynamics in living cells . . . . . . . . 86
4.3.1 Setting of the problem . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
4.3.2 The main convergence results . . . . . . . . . . . . . . . . . . . . . . . 92
5 Multiscale modeling of composite media with imperfect interfaces 95
5.1 Multiscale analysis in thermal diffusion problems in composite structures . . 96
5.1.1 Problem setting . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
5.1.2 The main results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
5.2 Diffusion problems with dynamical boundary conditions . . . . . . . . . . . . 100
5.3 Homogenization of a thermal problem with
ux jump . . . . . . . . . . . . . 102
5.3.1 Setting of the problem . . . . . . . . . . . . . . . . . . . . . . . . . . . 102
5.3.2 The macroscopic models . . . . . . . . . . . . . . . . . . . . . . . . . . 105
5.4 Other homogenization problems in composite media with imperfect interfaces 112
II Career Evolution and Development Plans 115
6 Scientic and academic background and research perspectives 117
6.1 Scientic and academic background . . . . . . . . . . . . . . . . . . . . . . . . 117
6.2 Further research directions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 119
6.3 Future plans . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 122
List of publications . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 121
Bibliography . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 123
Rezumat
Lucrarea de fat a, preg atit a pentru obt inerea atestatului de abilitare, cuprinde o select ie a
rezultatelor stiint ice pe care le-am obt inut ^ n domeniul teoriei omogeniz arii dup a dob^ andirea
titlului de doctor ^ n Matematic a. Principala motivat ie din spatele acestui demers o constituie
dorint a de a descrie succint stadiul actual al cunoa sterii ^ n domeniul teoriei omogeniz arii, de
a oferi o imagine de ansamblu asupra contribut iilor mele ^ n acest vast domeniu si de a discuta
c^ ateva probleme deschise si c^ ateva posibile perspective de dezvoltare a viitoarei mele cariere
stiint ice si academice.
Teza se bazeaz a pe o parte dintre contribut iile mele originale la aplicarea teoriei omoge-
niz arii, contribut ii cont inute^ n dou azeci si cinci de articole publicate sau trimse spre publicare
^ n reviste internat ionale cu o larg a audient a, incluz^ and nu doar matematicieni, ci si zicieni,
ingineri si cercet atori din diferite domenii aplicative. Multe dintre rezultatele cont inute ^ n
publicat iile pe care le-am selectat pentru a sprijini aceast a aplicat ie sunt str^ ans legate sau
motivate de probleme de interes practic. Voi ^ ncerca s a ofer o imagine de ansamblu, de sine-
st at atoare, asupra contribut iilor mele si, acolo unde va necesar, s a dau mai multe detalii
care nu sunt prezente ^ n lucr arile corespunz atoare publicate, f ac^ and astfel ca principalele mele
rezultate s a e accesibile unui public mai larg, cu solide cuno stint e generale de matematic a,
dar nu neap arat expert i ^ n domeniul specic al teoriei omogeniz arii. Teza prezint a c^ ateva
dintre cele mai relevante rezultate pe le-am obt inut pe parcursul ultimilor cincisprezece ani
de cercetare efectuat a, singur a sau ^ n colaborare, ^ n patru arii majore: analiza de multiscar a
a proceselor de react ie-difuzie ^ n medii poroase, omogenizarea problemelor unilaterale, mo-
delarea, cu ajutorul sc arilor multiple, a mediilor compozite cu interfet e imperfecte si obt inerea
de modele matematice ^ n biologie si ^ n inginerie.
^In afar a de dou a scurte rezumate, unul ^ n limba rom^ an a si altul ^ n limba englez a, teza
cont ine dou a p art i si o bibliograe cuprinz atoare. Prima parte, structurat a ^ n cinci capitole,
este dedicat a prezent arii principalelor mele realiz ari stiint ice obt inute dup a nalizarea studi-
ilor doctorale. Dup a un capitol introductiv ^ n care este prezentat stadiul actual al cercet arii
^ n domeniul teoriei omogeniz arii si care ofer a cadrul general si o motivat ie pentru activitatea
mea de cercetare post-doctoral a ^ n acest domeniu, al doilea capitol cuprinde c^ ateva dintre
contribut iile mele referitoare la omogenizarea mi sc arilor reactive ^ n medii poroase. Mai pre-
cis, sunt prezentate rezultate originale de omogenizare pentru curgerile reactive neliniare
stat ionare ^ n medii poroase si pentru studiul fenomenelor de adsorbt ie neliniar a ^ n medii
poroase. Capitolul se bazeaz a pe lucr arile [72], [73], [209] si [211]. Al treilea capitol, bazat
9
10 Rezumat
pe articolele [51], [54], [75], [77] si [208], este dedicat omogeniz arii unor probleme unilaterale
^ n domenii perforate. Sunt cuprinse rezultate originale privitoare la omogenizarea unor pro-
bleme de tip Signorini si a unor probleme eliptice cu condit ii mixte pe frontiere ^ n domenii
perforate. Al patrulea capitol cont ine o serie de rezultate recente despre omogenizarea unor
modelele din biologie. Sunt discutate noi modele matematice pentru transportul ionic^ n medii
periodice ^ nc arcate electric, pentru studiul carcinogenezei ^ n celulele vii sau pentru analiza
dinamicii calciului ^ n celulele biologice. Rezultatele prezentate ^ n acest capitol sunt cuprinse
^ n lucr arile [210], [211], [212], [217], [218], [219] si [221]. Ultimul capitol al acestei prime p art i
sintetizeaz a cele mai importante rezultate pe care le-am obt inut ^ n domeniul transferului de
c aldur a ^ n materiale compozite cu interfet e imperfecte si se bazeaz a pe lucr arile [47], [48],
[49], [213], [214] si [215]. Denit iile not iunilor de baz a din teoria omogeniz arii si rezultatele
generale din analiza funct ional a care vor folosite pe parcursul acestei lucr ari pot g asite
^ n [30], [37], [45], [56], [62], [137], [147], [155], [160], [161] si [205].
A doua parte a acestei teze prezint a c^ ateva planuri de dezvoltare si de evolut ie ^ n carier a.
Dup a o scurt a trecere ^ n revist a a parcursului meu stiint ic si academic de p^ an a acum, sunt
prezentate direct iile viitoare de cercetare si c^ ateva planuri de dezvoltare pe termen scurt,
mediu si lung a carierei mele stiint ice si academice.
Lucrarea se ^ ncheie cu o bibliograe cuprinz atoare, menit a s a ilustreze stadiul actual al
cunoa sterii ^ n acest domeniu vast al teoriei omogeniz arii si al aplicat iilor sale.
Principalele mele contribut ii originale cont inute ^ n aceast a tez a pot sintetizate astfel:
efectuarea unui studiu riguros al proceselor neliniare de react ie-difuzie^ n medii poroase,
care includ difuzie, react ii chimice si diferite tipuri de rate de adsorbt ie;
obt inerea unor rezultate de omogenizare pentru probleme unilaterale^ n medii perforate;
elaborarea de noi modele matematice pentru fenomenele de transport ionic ^ n medii
periodice ^ nc arcate electric;
derivarea de noi modele matematice pentru studiul proceselor neliniare de carcinogene-
z a ^ n celulele umane si al dinamicii calciului ^ n celulele vii;
efectuarea unei analize asimptotice riguroase pentru procese de difuzie termic a ^ n
structuri compozite;
ranarea studiului problemelor de difuzie cu condit ii dinamice pe frontier a;
obt inerea de noi modele matematice pentru probleme de difuzie cu salt ^ n
ux.
Rezultatele incluse ^ n aceast a tez a au fost obt inute singur a sau ^ n str^ ans a colaborare cu
mai multe institut ii academice si de cercetare din Rom^ ania sau din str ain atate. Sunt profund
recunosc atoare tuturor co-autorilor mei, Profesor dr. C. Conca, Profesor dr. F. Murat,
Profesor dr. J. I. D az, Profesor dr. A. Li~ n an, Profesor dr. H. I. Ene, Dr. A. C ap at ^ n a, Dr.
I. T ent ea, Dr. R. Bunoiu, pentru o frumoas a colaborare, pentru contribut ia lor important a,
pentru sfaturile utile si pentru discut iile interesante pe care le-am avut de-a lungul anilor.
Toate aceste rezultate ar putea deschide perspective noi si promit atoare pentru dezvolt ari
ulterioare si pentru viitoare colabor ari cu cercet atori reputat i din t ar a si din str ain atate.
Abstract
This manuscript, prepared to defend my Habilitation thesis, summarizes a selection of my
research results obtained in the eld of homogenization theory after defending my Ph.D.
thesis. The main motivation behind this endeavour is to brie
y describe the state of the art
in the eld of homogenization theory, to give an overview of my contributions in this broad
research area and to discuss some open problems and several perspectives I see for my future
scientic and academic career.
The thesis relies on some of my original contributions to the applications of the homoge-
nization theory, contained in twenty-ve articles already published or submitted for publica-
tion in international journals with a broad audience, including not only mathematicians, but
also physicists, engineers, and scientists from various applied elds. Many of the results in
the publications I selected to support my application are closely related to or motivated by
practical applications to real-life problems. I shall try to make a self-contained overview and,
where necessary, to give more details that are not present in the corresponding published
papers, making my main results accessible to an audience with strong, general mathematical
background, but not necessarily experts in the specic eld of homogenization theory.
The thesis is based on some of the most relevant results I obtained during the last fteen
years of research conducted, alone or in collaboration, in four major areas: multiscale analysis
of reaction-diffusion processes in porous media, upscaling in unilateral problems, multiscale
modeling of composite media with imperfect interfaces, and mathematical models in biology
and in engineering. Thus, the homogenization theory and its applications represent the core
of my scientic work done during these last fteen years.
Apart from two short abstracts, one in Romanian and another one in English, the thesis
comprises two parts and a comprehensive bibliography. The rst part, structured into ve
chapters, is devoted to the presentation of my main scientic achievements since the comple-
tion of my Ph.D. thesis. After a brief introductory chapter presenting the state of the art in
the eld of homogenization theory and offering the general framework and a motivation for
my post-doctoral research work in this area, the second chapter is divided in two distinct sec-
tions, summarizing my main contributions to the homogenization of reactive
ows in porous
media. More precisely, some original results for upscaling in stationary nonlinear reactive
ows in porous media and, also, results on nonlinear adsorption phenomena in porous media
are presented. The chapter relies on the papers [72], [73], [209], and [211]. The third chapter
is devoted to the homogenization of some relevant unilateral problems in perforated domains.
11
12 Abstract
More precisely, homogenization results for Signorini's type problems and for elliptic problems
with mixed boundary conditions in perforated media are presented. The chapter is based on
the articles [51], [54], [75], [77], and [208]. The fourth chapter contains some recent results
about homogenized models in biology. New mathematical models for ionic transport phe-
nomena in periodic charged media, for carcinogenesis in living cells or for analyzing calcium
dynamics in biological cells are discussed. The results presented in this chapter are contained
in [210], [211], [212], [217], [218], [219], and [221]. The last chapter of this rst part sum-
marizes the most important results I achieved, alone or in collaboration, in the eld of heat
transfer in composite materials with imperfect interfaces and is based on the papers [47], [48],
[49], [213], [214], and [215]. For the denitions of the basic notions in homogenization theory
and for well-known general results of functional analysis we shall use throughout this thesis,
we refer to [30], [37], [45], [56], [62], [137], [147], [155], [160], [161], and [205].
The second part of this thesis presents some career evolution and development plans.
After a brief review of my scientic and academic background, further research directions
and some future plans on my scientic and academic career are presented. I shall discuss
some short, medium and long term development plans. A brief description of some open
questions I would like to study in the future will be made, as well.
The thesis ends by a comprehensive bibliography, illustrating the state of the art in this
vast eld of homogenization theory and its applications.
My major original contributions contained in this habilitation thesis can be summarized
as follows:
performing a rigorous study of nonlinear reaction-diffusion processes in porous media,
including diffusion, chemical reactions and different types of adsorption rates;
obtaining new homogenization results for unilateral problems in perforated media;
elaborating new mathematical models for ionic transport phenomena in periodic
charged media;
getting original homogenization results for calcium dynamics in living cells;
deriving new nonlinear mathematical models for carcinogenesis in human cells;
performing a rigorous multiscale analysis of some relevant thermal diffusion processes
in composite structures;
rening the study of diffusion problems with dynamical boundary conditions;
obtaining new mathematical models for diffusion problems with
ux jump.
The results included in this thesis have been obtained alone or in close collaboration with
several academic and research institutions from Romania or from abroad. I am grateful to
all my co-authors, Professor C. Conca, Professor F. Murat, Professor J. I. D az, Professor A.
Li~ n an, Professor H. I. Ene, Dr. A. C ap at ^ n a, Dr. I. T ent ea, and Dr. R. Bunoiu, for a nice
collaboration, for their important contribution, for useful advices and interesting discussions.
I hope that all these results might open new and promising perspectives for further
developments and future collaborations with well-known scientists from Romania and from
abroad.
Part I
Main Scientic Achievements
13
Chapter 1
Introduction
In the last decades, there has been an explosive growth of interest in studying the macroscopic
properties of systems having a very complicated microscopic structure. In mechanics, physics,
chemistry, engineering, in material science or in biology, we are often led to consider boundary-
value problems in periodic media exhibiting multiple scales. It is widely recognized that mul-
tiscale techniques represent an essential tool for understanding the macroscopic properties
of such systems having a very complicated microscopic structure. A periodic distribution is
a realistic hypothesis in many situations with practical applications. Typically, in periodic
heterogeneous structures, the physical parameters, such as the electrical or thermal conduc-
tivity or the elastic coefficients, are discontinuous and, moreover, highly oscillating. For
example, in a composite material, constituted by the ne mixing of two or more components,
the physical parameters are obviously discontinuous and they are highly oscillating between
different values characterizing each distinct component. Therefore, the microscopic structure
becomes extremely complicated. If the period of the structure is very small compared to
the domain where we study the given system or, in other words, if the nonhomogeneities
are small compared to the global dimension of the structure, then an asymptotic analysis
becomes necessary. Two scales are important for a suitable description of the considered
structure: one which is comparable with the dimension of the period, called the microscopic
scale, and another one which is of the same order of magnitude as the global dimension of
our system, called the macroscopic scale . The main goal of the homogenization methods is to
pass from the microscopic scale to the macroscopic one; more precisely, using homogenization
methods, we try to describe the macroscopic properties of the nonhomogeneous system in
terms of the properties of its microscopic structure. The nonhomogeneous system is replaced
by a ctitious homogeneous one, with global characteristics which represent a good approx-
imation of the initial system. In this way, we are led to a general framework for obtaining
these macroscale properties, eliminating the huge difficulties related to the explicit determi-
nation of a solution of the microscopic problem. Also, from the point of view of numerical
computation, the homogenized equations will be easier to solve. This is due to the fact that
they are dened on a xed domain and they have, in general, simpler or even constant coeffi-
cients (the so-called effective orhomogenized coefficients ), while the microscopic equations are
15
16 Introduction
dened on a complicated domain, have rapidly oscillating coefficients, and satisfy nonlinear
boundary conditions. Let us remark that the dependence on the real microstructure is given
through the homogenized coefficients.
The study of the macroscopic properties of composite media was initiated by the physicists
Rayleigh, Maxwell, and Einstein. Around 1970, such problems were formulated in such a
way that they became interesting for mathematicians, as well, and this gave rise to a new
mathematical discipline, the homogenization theory . The rst rigorous developments of this
theory appeared with the seminal works of I. Babuka [28], E. De Giorgi and S. Spagnolo [86],
A. Bensoussan, J. L. Lions and G. Papanicolaou [37], and L. Tartar [205]. De Giorgi's notion
of Gamma-convergence marked also an important step in the development of this theory. F.
Murat and L. Tartar (see [178], [179], [180], and [206]) introduced the notion of compensated
compactness, which is an important tool to prove convergence results. A rigorous method,
the two-scale convergence method, was introduced by G. Nguetseng in 1989 [181] and was
further developed by G. Allaire in [1]. An extension to multiscale problems was obtained
by A.I. Ene and J. Saint Jean Paulin [104] and by G. Allaire and M. Briane [2]. In 1990,
T. Arbogast, J. Douglas, and U. Hornung [21] dened a dilation operator in order to study
homogenization problems in a periodic medium with double porosity. An alternative approach
was offered by the Bloch-wave homogenization method [76], which is a high frequency method
that can provide dispersion relations for wave propagation in periodic structures. Recently,
D. Cioranescu, A. Damlamian, P. Donato, and G. Griso combined the dilation technique with
ideas from nite element approximations to give rise to a very general method for studying
classical or multiscale periodic homogenization problems: the periodic unfolding method
(see, e.g., [56]). Let us nally mention that probabilistic and numerical methods, such as
the heterogeneous multiscale method, have been recently developed and successfully applied
to a broad category of problems of both practical and theoretical interest (see [100]). It is
important to emphasize that homogenization theory can be applied to non-periodic media,
as well. To this end, one can use G- orH-convergence techniques. Also, it is possible to deal
with general geometrical settings, without assuming periodicity or randomness.
Homogenization methods have been successfully applied to various problems, such as
the convective-diffusive transport in porous media, nonlinear elasticity problems, the study
of composite polymers, the study of nanocomposite materials, the modeling of interface
phenomena in biology and chemistry, or the problem of obtaining new composite materials
with applications in modern technology. The literature on this subject is vast (see, e.g.,
[61], [74], [63], [73], [75], and the references therein). We also mention here some remarkable
monographs dedicated to the mathematical problems of homogenization: [147], [29], [37],
[160], [164], [184], [198], [62], [69], [102].
Multiscale methods offer multiple possibilities for further developments and for useful
applications in many domains of contemporary science and technology. Their study is one
of the most active and fastest growing areas of modern applied mathematics, and denitely
one of the most interdisciplinary eld of mathematics.
My interest in this broad eld of homogenization theory started after defending my Ph.D.
Introduction 17
thesis at "Simion Stoilow" Institute of Mathematics of the Romanian Academy under the
supervision of Professor Horia I. Ene. I focused on the applications of the homogenization
theory to a wide category of problems arising in physics, chemistry, biology or engineering.
To summarize, my main research interests have been related to the following areas: mul-
tiscale analysis of reaction-diffusion processes in porous media, homogenization results for
unilateral problems, multiscale modeling of composite media with imperfect interfaces, and
mathematical models in biology and in engineering. My research activity in the eld of ho-
mogenization is interdisciplinary in its nature and in the last years I tried to publish my
results in more application-oriented high quality journals, with a broad audience, including
not only mathematicians, but also physicists, engineers, and scientists from various applied
elds, such as biology or geology.
The aim of this manuscript, prepared to defend my Habilitation thesis, is to give an
overview of my research work in the eld of homogenization theory. As a matter of fact,
the homogenization theory and its applications represent the core of my scientic work done
during the last fteen years. Many of the results presented herein are closely related to or
motivated by practical applications to real-life problems.
Apart from two short abstracts, one in Romanian and another one in English, the thesis
comprises two parts and a comprehensive bibliography. The rst part, structured into ve
chapters, is devoted to the presentation of my main scientic achievements since the comple-
tion of my Ph.D. thesis. After a brief introductory chapter presenting the state of the art
in the eld of homogenization theory and offering the general framework and a motivation
for my post-doctoral research work in this area, the second chapter is divided in two distinct
sections, summarizing my main contributions related to the homogenization of reactive
ows
in porous media. More precisely, some original results for upscaling in stationary nonlinear
reactive
ows in porous media and, also, results on nonlinear adsorption phenomena in porous
media are presented. The chapter relies on the papers [72], [73], [209], and [211]. The third
chapter is devoted to the homogenization of some relevant unilateral problems in perforated
domains. More precisely, some homogenization results for Signorini's type problems and for
elliptic problems with mixed boundary conditions in perforated media are presented. The
chapter is based on the articles [51], [54], [75], [77], and [208]. The fourth chapter contains
some recent results about homogenized models in biology. New mathematical models for
ionic transport phenomena in periodic charged media, for carcinogenesis in living cells or for
analyzing calcium dynamics in biological cells are discussed. The results presented in this
chapter are contained in the papers [210], [211], [212], [217], [218], [219], and [221]. The
last chapter of this rst part summarizes the most important results I achieved, alone or in
collaboration, in the eld of heat transfer in composite materials with imperfect interfaces
and is mainly based on the articles [47], [48], [49], [213], [214], and [215].
For the denitions of the basic notions in homogenization theory and for well-known
general results of functional analysis we shall use throughout this thesis, we refer to [30], [37],
[45], [56], [62], [137], [147], [155], [160], [161], and [205].
The second part of this thesis presents some career evolution and development plans.
18 Introduction
After a brief review of my scientic and academic background, further research directions
and some future plans on my scientic and academic career are presented.
The manuscript end by a comprehensive bibliography, illustrating the state of the art in
this vast eld of homogenization theory and its applications.
The thesis relies on some of my original contributions to the applications of homogeniza-
tion theory, contained in twenty-ve articles. The results included in this thesis have been
obtained during the last fteen years of research studies conducted alone or in collaboration
with various research institutions from Romania and from abroad. Let me emphasize that
most of this work is already published or submitted for publication in international journals.
Thus, in this thesis, I shall explicitly use some parts from my own articles, mentioning each
time the precise references to the corresponding original work.
I am the unique author in twelve of the papers on which this thesis is based. I am co-
author, with equal contribution, for the rest thirteen papers on which this thesis is based
(in Mathematics, the academic norm is to list equally contributed authors in alphabetical
order). I gratefully acknowledge the equal contribution of all my co-authors: Professor C.
Conca, Professor F. Murat, Professor J.I. D az, Professor A. Li~ n an, Professor H.I. Ene, Dr.
A. Capatina, Dr. I. T ent ea, and Dr. R. Bunoiu.
My main original contributions in the eld of homogenization theory, contained in this
habilitation thesis, can be summarized as follows (a more detailed description of a selection of
my results in the eld of homogenization theory will be presented in the following chapters):
performing a rigorous study of nonlinear reaction-diffusion processes in porous media,
including diffusion, chemical reactions, and different types of adsorption rates;
obtaining new homogenization results for unilateral problems in perforated media;
elaborating new mathematical models for ionic transport phenomena in periodic
charged media;
getting original homogenization results for calcium dynamics in living cells;
deriving new nonlinear mathematical models for carcinogenesis in human cells;
performing a rigorous multiscale analysis of some relevant thermal diffusion processes
in composite structures;
rening the study of diffusion problems with dynamical boundary conditions;
obtaining new mathematical models for thermal problems with
ux jump.
All these results might open new and promising perspectives for further developments and
future collaborations with well-known academic and research institutions from Romania and
from abroad.
Chapter 2
Homogenization of reactive
ows in
porous media
The problem of obtaining suitable global descriptions for some complex reactive
ows in
porous media was addressed in the literature by using various upscaling methods: heuristic
and empirically based methods, variational methods, stochastic methods, methods based on
homogenization, mixture theories, or volume averaging techniques. Also, the use of numerical
models for studying single-phase or multi-phase
ows in heterogeneous porous media has
received considerable attention in the last decades. However, even with the increases in the
power of computers, the complex multiscale structure of the analyzed media constitutes a
critical problem in the numerical treatment of such models and there is a considerable interest
in the development of upscaled or homogenized models in which the effective properties of
the medium vary on a coarse scale which is suitable for efficient computation and accurately
captures the in
uence of the ne-scale structure on the coarse-scale properties of the solution.
Porous media play an important role in many areas, such as hydrology (groundwater
ow, salt water intrusion into coastal aquifers), geology (petroleum reservoir engineering,
geothermal energy), chemical engineering (packed bed rectors, drying of granular materials),
mechanical engineering (heat exchangers, porous gas burners), the study of industrial ma-
terials (glass ber materials, brick manufacturing). There is an extensive literature on the
determination of the effective properties of heterogeneous porous media (see, e.g., [32] and
[137], and the references therein).
Transport processes in porous media have been extensively studied in last decades by
engineers, geologists, hydrologists, mathematicians, physicists. In particular, mathematical
modeling of chemical reactive
ows through porous media is a topic of huge practical im-
portance in many engineering, physical, chemical, and biological applications. Obtaining
suitable macroscopic laws for the processes in geometrically complex porous media (such as
soil, concrete, rock, or pellets) involving
ow, diffusion, convection, and chemical reactions
is a difficult task. The homogenization theory proves to be a very efficient tool by provid-
ing suitable techniques allowing us to pass from the microscopic scale to the macroscopic
one and to obtain suitable macroscale models. Since the seminal work of G.I. Taylor [207],
19
20 2.1. Upscaling in stationary reactive
ows in porous media
dispersion phenomena in porous media have attracted a lot of attention. There are many
formal or rigorous methods in the literature. We refer to [138] and [139] as one of the rst
works containing rigorous homogenization results for reactive
ows in porous media. By us-
ing the two-scale convergence method, coupled with monotonicity methods and compensated
compactness, the convergence of the homogenization procedure was proven for problems with
nonlinear reactive terms and nonlinear transmission conditions. Since then, many works have
been devoted to the homogenization of reactive transport in porous media (see [3], [25], [32],
[167], [172], [156], [139], [166], [157], [98] and the references therein). For instance, rigorous
homogenization results for reactive
ows with adsorption and desorption at the boundaries
of the perforations, for dominant P eclet numbers and Damkohler numbers, are obtained in
[9], [8], and [171]. For reactive
ows combined with the mechanics of cells, we refer to [146].
Rigorous homogenization techniques for obtaining the effective model for dissolution and pre-
cipitation in a complex porous medium were successfully applied in [157]. Solute transport in
porous media is also a topic of interest for chemists, geologists and environmental scientists
(see, e.g., [6] and [97]). Related problems, such that single or two-phase
ow or miscible
displacement problems were addressed in various papers (see, for instance, [16], [21], [22],
[170]). For an interesting survey on homogenization techniques applied to problems involving
ow, diffusion, convection, and reactions in porous media, we refer to [137].
In this chapter, some applications of the homogenization method to the study of reactive
ows in periodic porous media will be presented. The chapter represents a summary of the
results I obtained in this area, alone or in collaboration, and is based on the papers [72], [73],
[209], and [211].
2.1 Upscaling in stationary reactive
ows in porous media
We shall discuss now, following [72] and [211], some homogenization results for chemical
reactive
ows through porous media. For more details about the chemical aspects involved in
this kind of problems and, also, for some mathematical and historical backgrounds, we refer
to S. N. Antontsev et al. [20], J. Bear [32], J. I. D az [88], [91], [90], and U. Hornung [137]
and the references therein. We shall be concerned with a problem modeling the stationary
reactive
ow of a
uid conned in the exterior of some periodically distributed obstacles,
reacting on the boundaries of the obstacles. More precisely, the challenge in our rst paper
dedicated to this subject, namely [72], consists in dealing with Lipschitz or even non-Lipschitz
continuous reaction rates such as Langmuir or Freundlich kinetics, which, at that time, were
open cases in the literature. Our results represent a generalization of some of the results
in [137]. Using rigorous multiscale techniques, we derive a macroscopic model system for
such elliptic problems modeling chemical reactions on the grains of a porous medium. The
effective model preserves all the relevant information from the microscopic level. The case in
which chemical reactions arise inside the grains of a porous medium will be also discussed.
Also, we shall present some results obtained in [211], where we have analyzed the effective
behavior of the solution of a nonlinear problem arising in the modeling of enzyme catalyzed
Chapter 2. Homogenization of reactive
ows in porous media 21
reactions through the exterior of a domain containing periodically distributed reactive solid
obstacles.
2.1.1 The model problem
We consider an open smooth connected bounded set Ω in Rn(n3) and we insert in it a set
of"-periodically distributed reactive obstacles. In this way, we obtain an open set Ω", called
theexterior domain ;"represents a small parameter related to the characteristic size of the
reactive obstacles. More precisely, let Y= (0;1)nbe the unit cell in Rn. Denote by Fan
open subset of Ywith smooth boundary @Fsuch thatFY. We shall refer to Fas being
the elementary obstacle . We set
Y=YnF:
If"is a real parameter taking values in a sequence of positive numbers converging to zero,
for each"and for any integer vector k2Zn, setF"
kthe translated image of "Fby the vector
k,
F"
k="(k+F):
The setF"
krepresents the obstacles in Rn. Also, let us denote by F"the set of all the obstacles
contained in Ω, i.e.
F"=∪{
F"
kjF"
kΩ; k2Zn}
:
Set Ω"= ΩnF". Hence, Ω"is a periodic domain with periodically distributed obstacles
of size of the same order as the period. We remark that the obstacles do not intersect the
boundary@Ω. Let
S"=[f@F"
kjF"
kΩ; k2Zng:
So,
@Ω"=@Ω[S":
We denote by j!jthe Lebesgue measure of any measurable subset !Rnand, for an
arbitrary function 2L2(Ω"), we denote by e its extension by zero to the whole of Ω. Also,
throughout this thesis, by Cwe denote a generic xed strictly positive constant, whose value
can change from line to line.
The rst problem we present in this section concerns the stationary reactive
ow of a
uid conned in Ω", with concentration u", reacting on the boundary of the obstacles. A
simplied version of this kind of problems can be written as follows:
8
>>><
>>>: Df∆u"=fin Ω";
Df@u"
@="g(u") onS";
u"= 0 on@Ω:(2.1)
Here,is the unit exterior normal to Ω",f2L2(Ω) andS"is the boundary of the exterior
medium Ω nΩ". For simplicity, we assume that the
uid is homogeneous and isotropic,
with a constant diffusion coefficient Df>0. We can treat in a similar manner the more
22 2.1. Upscaling in stationary reactive
ows in porous media
general case in which, instead of considering constant diffusion coefficients, we work with
an heterogeneous medium represented by periodic symmetric bounded matrices which are
assumed to be uniformly coercive.
The semilinear boundary condition imposed on S"in problem (2.1) describes the chemical
reactions which take place locally at the interface between the reactive
uid and the obstacles.
In fact, from a chemical perspective, such a situation represents, equivalently, the effective
reaction on the walls of the chemical reactor between the
uid lling Ω"and a chemical
reactant located inside the rigid solid grains.
For the function g, which is assumed to be given, two representative situations will be
considered: the case in which gis a monotone smooth function satisfying the condition
g(0) = 0 and the case of a maximal monotone graph with g(0) = 0, i.e. the case in which gis
the subdifferential of a convex lower semicontinuous function G. These two general situations
are well illustrated by the following important practical examples:
a)g(v) =v
1 +v; ; > 0 (Langmuir kinetics) (2.2)
and
b)g(v) =jvjp 1v;0<p< 1 (Freundlich kinetics) : (2.3)
The exponent pis called the order of the reaction . We point out that if we assume f0, one
can prove (see, e.g. [90]) that u"0 in Ω nΩ"andu">0 in Ω", althoughu"is not uniformly
positive except in the case in which gis a monotone smooth function satisfying the condition
g(0) = 0, as, for instance, in example a). In fact, since u"is, in practical applications, a
concentration, we can impose suitable conditions on the data to ensure that 0 u"1 (see,
e.g., [89]).
As usual in homogenization, our goal is to obtain a suitable description of the asymptotic
behavior, as the small parameter "tends to zero, of the solution u"of problem (2.1) in such
domains.
2.1.2 The case of a smooth function g. The macroscopic model
Let us deal rst with the case of a smooth function g. We consider that gis a continuously
differentiable function, monotonously non-decreasing and such that g(v) = 0 if and only if
v= 0. Moreover, we suppose that there exist a positive constant Cand an exponent q, with
0q<n= (n 2), such thatdg
dvC(1 +jvjq): (2.4)
We introduce the functional space
V"={
v2H1(Ω")jv= 0 on@Ω}
;
endowed with the norm
∥v∥V"=∥∇v∥L2(Ω"):
Chapter 2. Homogenization of reactive
ows in porous media 23
The weak formulation of problem (2.1) is:
8
><
>:Findu"2V"such that
Df∫
Ω"∇u" ∇φdx+"∫
S"g(u")φd=∫
Ω"fφdx8φ2V":(2.5)
By classical existence results (see [45]), there exists a unique weak solution u"2V"\H2(Ω")
of problem (2.1). This solution being dened only on Ω", we need to extend it to the whole of
Ω to be able to state the convergence result. To this end, let us recall the following well-known
extension result (see [69]):
Lemma 2.1 There exists a linear continuous extension operator
P"2 L(L2(Ω");L2(Ω))\ L(V";H1
0(Ω))
and a positive constant C, independent of ", such that, for any v2V",
∥P"v∥L2(Ω)C∥v∥L2(Ω")
and
∥∇P"v∥L2(Ω)C∥∇v∥L2(Ω"):
Therefore, we have the following Poincar e's inequality in V":
Lemma 2.2 There exists a positive constant C, independent of ", such that
∥v∥L2(Ω")C∥∇v∥L2(Ω")for anyv2V":
The main convergence result for this case is stated in the following theorem.
Theorem 2.3 (Theorem 2.3 in [72]) There exists an extension P"u"of the solution u"of
the variational problem (2.5) such that
P"u"⇀u weakly inH1
0(Ω);
whereuis the unique solution of
8
><
>: n∑
i;j=1qij@2u
@xi@xj+j@Fj
jYjg(u) =finΩ;
u= 0 on@Ω:(2.6)
Here,Q= ((qij))is the homogenized matrix, whose entries are dened as follows:
qij=Df(
ij+1
jYj∫
Y@j
@yidy)
; (2.7)
in terms of the functions i; i= 1;:::;n; solutions of the cell problems
8
>>>><
>>>>: ∆i= 0 inY;
@(i+yi)
@= 0 on@F;
iY periodic.(2.8)
The homogenized matrix Qis constant, symmetric, and positive-denite.
24 2.1. Upscaling in stationary reactive
ows in porous media
Proof. The proof of this theorem is divided into four steps.
First step. Letu"2V"be the solution of the variational problem (2.5) and let P"u"be
the extension of u"inside the obstacles given by Lemma 2.1. Taking φ=u"as a test function
in (2.5), using Schwartz and Poincar e's inequalities, we get
∥P"u"∥H1
0(Ω)C
and, by passing to a subsequence, still denoted by P"u", we can suppose that there exists
u2H1
0(Ω) such that
P"u"⇀u weakly inH1
0(Ω): (2.9)
It remains to determine the limit equation satised by u.
Second step . In order to obtain the limit equation satised by u, we need to pass to the
limit in (2.5). The most delicate part, and, in fact, the main novelty brought by our paper, is
the passage to the limit, in the variational formulation (2.5) of problem (2.1), in the nonlinear
term on the boundary of the grains, i.e. in the second term in the left-hand side of (2.5).
To this end, we introduce, for any h2Ls′(@F), 1s′ 1 , the linear form "
honW1;s
0(Ω)
dened by
⟨"
h;φ⟩="∫
S"h(x
")
φd8φ2W1;s
0(Ω);
with 1=s+ 1=s′= 1. Then (see [61]),
"
h!hstrongly in ( W1;s
0(Ω))′; (2.10)
where
⟨h;φ⟩=h∫
Ωφdx;
with
h=1
jYj∫
@Fh(y) d:
Ifh2L1(@F) or ifhis constant, we have "
h!hstrongly in W 1;1(Ω) and we denote
by"the above introduced measure in the particular case in which h= 1. Notice that in
this casehbecomes1=j@Fj=jYj. We shall prove now that for any φ2 D(Ω) and for any
v"⇀v weakly inH1
0(Ω), one has
φg(v")⇀φg (v) weakly in W1;q
0(Ω); (2.11)
where
q=2n
q(n 2) +n:
To this end, let us remark that
sup∥∇g(v")∥Lq(Ω)<1: (2.12)
Indeed, using the growth condition (2.4) imposed to g, we have
∫
Ω@g
@xi(v")q
dxC∫
Ω(
1 +jv"jqq)@v"
@xiq
dx
Chapter 2. Homogenization of reactive
ows in porous media 25
C(
1 +(∫
Ωjv"jqq
dx)1=
)(∫
Ωj∇v"jqdx)1=
;
where we took
andsuch thatq= 2, 1=
+ 1== 1 andqq
= 2n=(n 2). Notice that
it is from here that we get
q=2n
q(n 2) +n:
Also, due to the fact that 0 q<n= (n 2), it follows that q>1. Since
sup∥v"∥
L2n
n 2(Ω)<1;
we easily get (2.12). Therefore, to obtain (2.11), it remains only to show that
g(v")!g(v) strongly in Lq(Ω): (2.13)
But this convergence is a direct consequence of the following well-known result (see [80], [155]
and [160]):
Theorem 2.4 LetG: ΩR!Rbe a Carath eodory function, i.e.
a) for every vthe function G(;v)is measurable with respect to x2Ω:
b) for every (a.e.) x2Ω, the function G(x;)is continuous with respect to v.
Moreover, if we assume that there exists a positive constant Csuch that
jG(x;v)j C(
1 +jvjr=t)
;
withr1andt <1, then the map v2Lr(Ω)7!G(x;v(x))2Lt(Ω)is continuous in the
strong topologies.
Indeed, since
jg(v)j C(1 +jvjq+1);
applying the above theorem for G(x;v) =g(v),t=qandr= (2n=(n 2)) r′, withr′>0
such thatq+ 1< r=t and using the compact injection H1(Ω),!Lr(Ω) we obtain (2.13).
Now, from (2.10), written for h= 1, and (2.11) written for v"=P"u", we get
⟨";φg(P"u")⟩ !j@Fj
jYj∫
Ωφg(u) dx8φ2 D(Ω) (2.14)
and this completes the proof of this step.
Third step . Let"be the gradient of u"in Ω"and let us denote by e"its extension by
zero to the whole of Ω. Then, e"is bounded in ( L2(Ω))nand, as a consequence, there exists
2(L2(Ω))nsuch that
e"⇀ weakly in (L2(Ω))n: (2.15)
Let us identify now the equation satised by . If we take φ2 D(Ω), from (2.5) we get
∫
Ωe" ∇φdx+"∫
S"g(u")φd=∫
ΩΩ"fφdx (2.16)
26 2.1. Upscaling in stationary reactive
ows in porous media
and we can pass to the limit, with "!0, in all the terms of (2.16). For the rst one, we have
lim
"!0∫
Ωe" ∇φdx=∫
Ω ∇φdx: (2.17)
For the second term, using (2.14), we obtain
lim
"!0"∫
S"g(u")φd=j@Fj
jYj∫
Ωg(u)φdx: (2.18)
It is not difficult to pass to the limit in the right-hand side of (2.16). Since
Ω"f ⇀jYj
jYjfweakly inL2(Ω);
we get
lim
"!0∫
ΩΩ"fφdx=jYj
jYj∫
Ωfφdx: (2.19)
Putting together (2.17)-(2.19), we have
∫
Ω ∇φdx+j@Fj
jYj∫
Ωg(u)φdx=jYj
jYj∫
Ωfφdx8φ2 D(Ω):
Thus,satises
div+j@Fj
jYjg(u) =jYj
jYjfin Ω (2.20)
and it remains now only to identify .
Fourth step. For identifying , we shall use the solutions of the local problems (2.8). For
any xedi= 1;:::;n , we dene
i"(x) ="(
i(x
")
+yi)
8×2Ω"; (2.21)
wherey=x=". By periodicity, it follows that
P"i"⇀x iweakly inH1(Ω): (2.22)
Let"
ibe the gradient of i"in Ω"ande"
ibe the extension by zero of "
iinside the holes.
We have
(
e"
i)
j=(g@i"
@xj)
=(
^@i
@yj(y))
+ijY
and, therefore,
(
e"
i)
j⇀1
jYj(∫
Y@i
@yjdy+jYjij)
=jYj
jYjqijweakly inL2(Ω): (2.23)
We notice that "
isatises{
div"
i= 0 in Ω";
"
i= 0 onS":(2.24)
Chapter 2. Homogenization of reactive
ows in porous media 27
Letφ2 D(Ω). Multiplying the rst equation in (2.24) by φu"and integrating by parts on
Ω", we obtain∫
Ω""
i ∇φu"dx+∫
Ω""
i ∇u"φdx= 0:
Thus,∫
Ωe"
i ∇φP"u"dx+∫
Ω""
i ∇u"φdx= 0: (2.25)
Takingφi"as a test function in (2.5), we get
∫
Ω"(∇u" ∇φ)i"dx+∫
Ω"(∇u" ∇i")φdx+"∫
S"g(u")φi"d=∫
Ω"fφi"dx;
which, using the denitions of e"ande"
i, leads to
∫
Ωe" ∇φP"i"dx+∫
Ω"∇u""
iφdx+"∫
S"g(u")φi"d=∫
ΩfΩ"φP"i"dx:
From (2.25), we obtain
∫
Ωe" ∇φP"i"dx ∫
Ωe"
i ∇φP"i"dx+"∫
S"g(u")φi"d=∫
ΩfΩ"φP"i"dx:(2.26)
We pass now to the limit in (2.26). By using (2.15) and (2.22), we get
lim
"!0∫
Ωe" ∇φP"i"dx=∫
Ω ∇φxidx: (2.27)
From (2.9) and (2.23) we obtain
lim
"!0∫
Ωe"
i ∇φP"u"dx=jYj
jYj∫
Ωqi ∇φudx; (2.28)
whereqiis the vector having the j-component equal to qij.
Since the boundary of Fis of classC2,P"i"2W1;1(Ω) andP"i"!xistrongly in
L1(Ω):Then, since g(P"u")P"i"!g(u)xistrongly in Lq(Ω) andg(P"u")P"i"is bounded
inW1;q(Ω), we have g(P"u")P"i"⇀g(u)xiweakly inW1;q(Ω). Thus,
lim
"!0"∫
S"g(u")φi"d=j@Fj
jYj∫
Ωg(u)φxidx: (2.29)
For the limit of the right-hand side of (2.26), since
Ω"f ⇀jYj
jYjfweakly inL2(Ω);
using (2.22), we have
lim
"!0∫
ΩfΩ"φP"i"dx=jYj
jYj∫
Ωfφx idx: (2.30)
28 2.1. Upscaling in stationary reactive
ows in porous media
Thus, we get
∫
Ω ∇φxidx jYj
jYj∫
Ωqi ∇φudx+j@Fj
jYj∫
Ωg(u)φxidx=jYj
jYj∫
Ωfφx idx: (2.31)
Using Green's formula and equation (2.20), we get
∫
Ω ∇xiφdx+jYj
jYj∫
Ωqi ∇uφdx= 0 in Ω:
The above equality is valid for any φ2 D(Ω) and this implies that
∇xi+jYj
jYjqi ∇u= 0 in Ω: (2.32)
Writing (2.32) by components, differentiating with respect to xi;summing after iand using
(2.19), we are led to
jYj
jYjn∑
i;j=1qij@2u
@xi@xj= div= jYj
jYjf+j@Tj
jYjg(u);
which means that uveries
n∑
i;j=1qij@2u
@xi@xj+j@Fj
jYjg(u) =fin Ω:
Sinceu2H1
0(Ω) (i.e.u= 0 on@Ω) anduis uniquely determined, the whole sequence P"u"
converges to uand Theorem 2.3 is proven.
Remark 2.5 The right scaling "in front of the function gmodeling the contribution of the
nonlinear reactions on the boundary of the grains leads in the limit to the presence of a new
term distributed all over the domain Ω. Also, let us emphasize again that if we assume f0,
the function gin example a)is indeed a particular example of our rst model situation.
Remark 2.6 The results in [72] are obtained for the case n3. All of them are still valid,
under our assumptions, in the case n= 2. Of course, for this case, n=(n 2)has to be
replaced by +1and, hence, (2.4) holds true for 0q<1. The results of this section could
be obtained, under our assumptions, without imposing any growth condition for the nonlinear
functiong(see [209]).
Remark 2.7 In [72], the proof of Theorem 2.3 was done by using the so-called energy method
of L. Tartar (see [205]). We point out that one can use also the recently developed periodic
unfolding method, introduced by Cioranescu, Damlamian, Donato, Griso and Zaki (see, e.g.,
[66], [56], [57], and [64]), which, apart from a signicant simplication in the proof, allows
us to deal with more general media, since we are not forced to use extension operators.
Chapter 2. Homogenization of reactive
ows in porous media 29
2.1.3 The case of a non-smooth function g. The macroscopic model
The case in which the function gappearing in (2.1) is a single-valued maximal monotone
graph in RR, satisfying the condition g(0) = 0, is also treated in [72]. If we denote by
D(g) the domain of g, i.e.D(g) =f2Rjg()̸=∅g, then we suppose that D(g) =R.
Moreover, we assume that gis continuous and there exist C0 and an exponent q, with
0q<n= (n 2), such that
jg(v)j C(1 +jvjq): (2.33)
Notice that the second important practical example b) mentioned above is a particular ex-
ample of such a single-valued maximal monotone graph.
In this case, there exists a lower semicontinuous convex function GfromRto ] 1;+1],
Gproper, i.e. G̸+1, such that gis the subdifferential of G,g=@G. Let
G(v) =∫v
0g(s) ds:
We dene the convex set
K"={
v2V"jG(v)jS"2L1(S")}
:
For a given function f2L2(Ω), the weak solution of the problem (2.1) is also the unique
solution of the following variational inequality:
8
>>><
>>>:ndu"2K"such that
Df∫
Ω"Du"D(v"
