Asymptotic homogenization

Last updated

In mathematics and physics, homogenization is a method of studying partial differential equations with rapidly oscillating coefficients, [1] [2] [3] such as

Contents

where is a very small parameter and is a 1-periodic coefficient: , .

It turns out that the study of these equations is also of great importance in physics and engineering, since equations of this type govern the physics of inhomogeneous or heterogeneous materials. Of course, all matter is inhomogeneous at some scale, but frequently it is convenient to treat it as homogeneous. A good example is the continuum concept which is used in continuum mechanics. Under this assumption, materials such as fluids, solids, etc. can be treated as homogeneous materials and associated with these materials are material properties such as shear modulus, elastic moduli, etc.

Frequently, inhomogeneous materials (such as composite materials) possess microstructure and therefore they are subjected to loads or forcings which vary on a lengthscale which is far bigger than the characteristic lengthscale of the microstructure. In this situation, one can often replace the equation above with an equation of the form

where is a constant tensor coefficient and is known as the effective property associated with the material in question. It can be explicitly computed as

from 1-periodic functions satisfying:

This process of replacing an equation with a highly oscillatory coefficient with one with a homogeneous (uniform) coefficient is known as homogenization. This subject is inextricably linked with the subject of micromechanics for this very reason.

In homogenization one equation is replaced by another if for small enough , provided in some appropriate norm as .

As a result of the above, homogenization can therefore be viewed as an extension of the continuum concept to materials which possess microstructure. The analogue of the differential element in the continuum concept (which contains enough atom, or molecular structure to be representative of that material), is known as the "Representative Volume Element" [4] in homogenization and micromechanics. This element contains enough statistical information about the inhomogeneous medium in order to be representative of the material. Therefore averaging over this element gives an effective property such as above.

Classical results of homogenization theory [1] [2] [3] were obtained for media with periodic microstructure modeled by partial differential equations with periodic coefficients. These results were later generalized to spatially homogeneous random media modeled by differential equations with random coefficients which statistical properties are the same at every point in space. [5] [6] In practice, many applications require a more general way of modeling that is neither periodic nor statistically homogeneous. For this end the methods of the homogenization theory have been extended to partial differential equations, which coefficients are neither periodic nor statistically homogeneous (so-called arbitrarily rough coefficients). [7] [8]

The method of asymptotic homogenization

Mathematical homogenization theory dates back to the French, Russian and Italian schools. [1] [2] [3] [9] The method of asymptotic homogenization proceeds by introducing the fast variable and posing a formal expansion in :

which generates a hierarchy of problems. The homogenized equation is obtained and the effective coefficients are determined by solving the so-called "cell problems" for the function .

See also

Notes

  1. 1 2 3 Sanchez-Palencia, E. (1980). Non-homogeneous media and vibration theory. Lecture Notes in Physics. 127. Springer Verlag. doi:10.1007/3-540-10000-8. ISBN   978-3-540-10000-3.
  2. 1 2 3 Bakhvalov, N.; Panasenko, G. (1989). Homogenisation: Averaging Processes in Periodic Media. Mathematics and its Applications. Dordrecht: Kluwer. doi:10.1007/978-94-009-2247-1. ISBN   978-94-010-7506-0.
  3. 1 2 3 Bensoussan, A.; Lions, J.L.; Papanicolaou, G. (1978). Asymptotic Analysis for Periodic Structures. Studies in Mathematics and its Applications. Amsterdam: North-Holland. ISBN   0-444-85172-0.
  4. Ostoja-Starzewski, M. (2007). Microstructural randomness and scaling in materials. Modern Mechanics and Mathematics. Chapman and Hall/CRC Press. ISBN   9781584884170.
  5. Kozlov, S.M. (1979). "Homogenization of Random Operators". Mat. Sbornik. 109 (151): 188–202. (English transl.: Math. USSR, Sb. 37:2, 1980, pp. 167-180)
  6. Papanicolaou, G. C.; Varadhan, S.R. (1981). "Boundary Value Problems with Rapidly Oscillating Coefficients" (PDF). Seria Colloq. Math. Society Janos Bolyai. Amsterdam. 27: 835–873.
  7. Berlyand, L.; Owhadi, H. (November 2010). "Flux Norm Approach to Finite Dimensional Homogenization Approximations with Non-Separated Scales and High Contrast". Archive for Rational Mechanics and Analysis. 198 (2): 677–721. arXiv: 0901.1463 . Bibcode:2010ArRMA.198..677B. doi:10.1007/s00205-010-0302-1.
  8. Målqvist, A.; Peterseim, D. (2014). "Localization of elliptic multiscale problems". Mathematics of Computation. 83 (290): 2583–2603. doi: 10.1090/S0025-5718-2014-02868-8 .
  9. Dal Maso, G. (1993). An Introduction to Γ-Convergence. Progress in Nonlinear Differential Equations and Their Applications. Birkhauser. doi:10.1007/978-1-4612-0327-8. ISBN   9780817636791.

Related Research Articles

Laplaces equation Second order partial differential equation

In mathematics and physics, Laplace's equation is a second-order partial differential equation named after Pierre-Simon Laplace who first studied its properties. This is often written as

In mathematics, a recurrence relation is an equation that recursively defines a sequence or multidimensional array of values, once one or more initial terms are given; each further term of the sequence or array is defined as a function of the preceding terms.

Del, or nabla, is an operator used in mathematics, in particular in vector calculus, as a vector differential operator, usually represented by the nabla symbol . When applied to a function defined on a one-dimensional domain, it denotes its standard derivative as defined in calculus. When applied to a field, it may denote the gradient of a scalar field, the divergence of a vector field, or the curl (rotation) of a vector field, depending on the way it is applied.

Heat equation

In mathematics and physics, the heat equation is a certain partial differential equation. Solutions of the heat equation are sometimes known as caloric functions. The theory of the heat equation was first developed by Joseph Fourier in 1822 for the purpose of modeling how a quantity such as heat diffuses through a given region.

Fokker–Planck equation Partial differential equation

In statistical mechanics, the Fokker–Planck equation is a partial differential equation that describes the time evolution of the probability density function of the velocity of a particle under the influence of drag forces and random forces, as in Brownian motion. The equation can be generalized to other observables as well. It is named after Adriaan Fokker and Max Planck, and is also known as the Kolmogorov forward equation, after Andrey Kolmogorov, who independently discovered the concept in 1931. When applied to particle position distributions, it is better known as the Smoluchowski equation, and in this context it is equivalent to the convection–diffusion equation. The case with zero diffusion is known in statistical mechanics as the Liouville equation. The Fokker–Planck equation is obtained from the master equation through Kramers–Moyal expansion.

In mathematics, the Laplace operator or Laplacian is a differential operator given by the divergence of the gradient of a function on Euclidean space. It is usually denoted by the symbols ∇·∇, 2 or Δ. In a Cartesian coordinate system, the Laplacian is given by the sum of second partial derivatives of the function with respect to each independent variable. In other coordinate systems, such as cylindrical and spherical coordinates, the Laplacian also has a useful form. Informally, the Laplacian Δf(p) of a function f at a point p measures by how much the average value of f over small spheres or balls centered at p deviates from f(p).

In the calculus of variations, a field of mathematical analysis, the functional derivative relates a change in a functional to a change in a function on which the functional depends.

In mathematics, a Green's function is the impulse response of an inhomogeneous linear differential operator defined on a domain with specified initial conditions or boundary conditions.

Blochs theorem

In condensed matter physics, Bloch's theorem states that solutions to the Schrödinger equation in a periodic potential take the form of a plane wave modulated by a periodic function. Mathematically, they are written:

In mathematics, a linear differential equation is a differential equation that is defined by a linear polynomial in the unknown function and its derivatives, that is an equation of the form

Geometrical optics, or ray optics, is a model of optics that describes light propagation in terms of rays. The ray in geometric optics is an abstraction useful for approximating the paths along which light propagates under certain circumstances.

In classical electromagnetism, reciprocity refers to a variety of related theorems involving the interchange of time-harmonic electric current densities (sources) and the resulting electromagnetic fields in Maxwell's equations for time-invariant linear media under certain constraints. Reciprocity is closely related to the concept of Hermitian operators from linear algebra, applied to electromagnetism.

Maxwell stress tensor

The Maxwell stress tensor is a symmetric second-order tensor used in classical electromagnetism to represent the interaction between electromagnetic forces and mechanical momentum. In simple situations, such as a point charge moving freely in a homogeneous magnetic field, it is easy to calculate the forces on the charge from the Lorentz force law. When the situation becomes more complicated, this ordinary procedure can become impossibly difficult, with equations spanning multiple lines. It is therefore convenient to collect many of these terms in the Maxwell stress tensor, and to use tensor arithmetic to find the answer to the problem at hand.

Covariant formulation of classical electromagnetism

The covariant formulation of classical electromagnetism refers to ways of writing the laws of classical electromagnetism in a form that is manifestly invariant under Lorentz transformations, in the formalism of special relativity using rectilinear inertial coordinate systems. These expressions both make it simple to prove that the laws of classical electromagnetism take the same form in any inertial coordinate system, and also provide a way to translate the fields and forces from one frame to another. However, this is not as general as Maxwell's equations in curved spacetime or non-rectilinear coordinate systems.

Mathematical descriptions of the electromagnetic field Formulations of electromagnetism

There are various mathematical descriptions of the electromagnetic field that are used in the study of electromagnetism, one of the four fundamental interactions of nature. In this article, several approaches are discussed, although the equations are in terms of electric and magnetic fields, potentials, and charges with currents, generally speaking.

Plane wave expansion method (PWE) refers to a computational technique in electromagnetics to solve the Maxwell's equations by formulating an eigenvalue problem out of the equation. This method is popular among the photonic crystal community as a method of solving for the band structure of specific photonic crystal geometries. PWE is traceable to the analytical formulations, and is useful in calculating modal solutions of Maxwell's equations over an inhomogeneous or periodic geometry. It is specifically tuned to solve problems in a time-harmonic forms, with non-dispersive media.

Spinodal decomposition

Spinodal decomposition occurs when one thermodynamic phase spontaneously separates into two phases. Decomposition occurs in the absence of nucleation because certain fluctuations in the system reduce the free energy. As a result, the phase change occurs immediately. There is no waiting, as there typically is when there is a nucleation barrier.

Matrix representation of Maxwells equations

In electromagnetism, a branch of fundamental physics, the matrix representations of the Maxwell's equations are a formulation of Maxwell's equations using matrices, complex numbers, and vector calculus. These representations are for a homogeneous medium, an approximation in an inhomogeneous medium. A matrix representation for an inhomogeneous medium was presented using a pair of matrix equations. A single equation using 4 × 4 matrices is necessary and sufficient for any homogeneous medium. For an inhomogeneous medium it necessarily requires 8 × 8 matrices.

In optics, the Ewald–Oseen extinction theorem, sometimes referred to as just "extinction theorem", is a theorem that underlies the common understanding of scattering. It is named after Paul Peter Ewald and Carl Wilhelm Oseen, who proved the theorem in crystalline and isotropic media, respectively, in 1916 and 1915. Originally, the theorem applied to scattering by an isotropic dielectric objects in free space. The scope of the theorem was greatly extended to encompass a wide variety of bianisotropic media.

Magnetic diffusion refers to the motion of magnetic field lines, typically in the presence of a conducting solid or fluid such as a plasma. The motion of magnetic fields is described by the magnetic diffusion equation and is due primarily to induction and diffusion of magnetic fields through the material. The magnetic diffusion equation is a partial differential equation commonly used in physics. Understanding the phenomenon is essential to magnetohydrodynamics and has important consequences in astrophysics, geophysics, and electrical engineering.

References