Spinodal decomposition

Last updated
Microstructural evolution under the Cahn-Hilliard equation, demonstrating distinctive coarsening and phase separation. CahnHilliard Animation.gif
Microstructural evolution under the Cahn–Hilliard equation, demonstrating distinctive coarsening and phase separation.

Spinodal decomposition is a mechanism by which a single thermodynamic phase spontaneously separates into two phases (without nucleation). [1] Decomposition occurs when there is no thermodynamic barrier to phase separation. As a result, phase separation via decomposition does not require the nucleation events resulting from thermodynamic fluctuations, which normally trigger phase separation.

Contents

Spinodal decomposition is observed when mixtures of metals or polymers separate into two co-existing phases, each rich in one species and poor in the other. [2] When the two phases emerge in approximately equal proportion (each occupying about the same volume or area), characteristic intertwined structures are formed that gradually coarsen (see animation). The dynamics of spinodal decomposition is commonly modeled using the Cahn–Hilliard equation.

Spinodal decomposition is fundamentally different from nucleation and growth. When there is a nucleation barrier to the formation of a second phase, time is taken by the system to overcome that barrier. As there is no barrier (by definition) to spinodal decomposition, some fluctuations (in the order parameter that characterizes the phase) start growing instantly. Furthermore, in spinodal decomposition, the two distinct phases start growing in any location uniformly throughout the volume, whereas a nucleated phase change begins at a discrete number of points.

Spinodal decomposition occurs when a homogenous phase becomes thermodynamically unstable. An unstable phase lies at a maximum in free energy. In contrast, nucleation and growth occur when a homogenous phase becomes metastable. That is, another biphasic system becomes lower in free energy, but the homogenous phase remains at a local minimum in free energy, and so is resistant to small fluctuations. J. Willard Gibbs described two criteria for a metastable phase: that it must remain stable against a small change over a large area. [3]

History

In the early 1940s, Bradley reported the observation of sidebands around the Bragg peaks in the X-ray diffraction pattern of a Cu-Ni-Fe alloy that had been quenched and then annealed inside the miscibility gap. Further observations on the same alloy were made by Daniel and Lipson, who demonstrated that the sidebands could be explained by a periodic modulation of composition in the <100> directions. From the spacing of the sidebands, they were able to determine the wavelength of the modulation, which was of the order of 100 angstroms (10 nm).

The growth of a composition modulation in an initially homogeneous alloy implies uphill diffusion or a negative diffusion coefficient. Becker and Dehlinger had already predicted a negative diffusivity inside the spinodal region of a binary system, but their treatments could not account for the growth of a modulation of a particular wavelength, such as was observed in the Cu-Ni-Fe alloy. In fact, any model based on Fick's law yields a physically unacceptable solution when the diffusion coefficient is negative.

The first explanation of the periodicity was given by Mats Hillert in his 1955 Doctoral Dissertation at MIT. Starting with a regular solution model, he derived a flux equation for one-dimensional diffusion on a discrete lattice. This equation differed from the usual one by the inclusion of a term, which allowed for the effect of the interfacial energy on the driving force of adjacent interatomic planes that differed in composition. Hillert solved the flux equation numerically and found that inside the spinodal it yielded a periodic variation of composition with distance. Furthermore, the wavelength of the modulation was of the same order as that observed in the Cu-Ni-Fe alloys. [4] [5]

Building on Hillert's work, a more flexible continuum model was subsequently developed by John W. Cahn and John Hilliard, who included the effects of coherency strains as well as the gradient energy term. The strains are significant in that they dictate the ultimate morphology of the decomposition in anisotropic materials. [6] [7] [8]

Cahn–Hilliard model for spinodal decomposition

Free energies in the presence of small amplitude fluctuations, e.g. in concentration, can be evaluated using an approximation introduced by Ginzburg and Landau to describe magnetic field gradients in superconductors. This approach allows one to approximate the free energy as an expansion in terms of the concentration gradient , a vector. Since free energy is a scalar and we are probing near its minima, the term proportional to is negligible. The lowest order term is the quadratic expression , a scalar. Here is a parameter that controls the free energy cost of variations in concentration .

The Cahn–Hilliard free energy is then

where is the bulk free energy per unit volume of the homogeneous solution, and the integral is over the volume of the system.

We now want to study the stability of the system with respect to small fluctuations in the concentration , for example a sine wave of amplitude and wavevector , for the wavelength of the concentration wave. To be thermodynamically stable, the free energy change due to any small amplitude concentration fluctuation , must be positive.

We may expand about the average composition co as follows:

and for the perturbation the free energy change is

When this is integrated over the volume , the gives zero, while and integrate to give . So, then [9]

As , thermodynamic stability requires that the term in brackets be positive. The is always positive but tends to zero at small wavevectors, large wavelengths. Since we are interested in macroscopic fluctuations, , stability requires that the second derivative of the free energy be positive. When it is, there is no spinodal decomposition, but when it is negative, spinodal decomposition will occur. Then fluctuations with wavevectors become spontaneously unstable, where the critical wave number is given by:

which corresponds to a fluctuations above a critical wavelength

Dynamics of spinodal decomposition when molecules move via diffusion

Spinodal decomposition can be modeled using a generalized diffusion equation: [10] [11] [12]

for the chemical potential and the mobility. As pointed out by Cahn, this equation can be considered as a phenomenological definition of the mobility M, which must by definition be positive. [13] It consists of the ratio of the flux to the local gradient in chemical potential. The chemical potential is a variation of the free energy and when this is the Cahn–Hilliard free energy this is [10]

and so

and now we want to see what happens to a small concentration fluctuation - note that now it has time dependence as a wavevector dependence. Here is a growth rate. If then the perturbation shrinks to nothing, the system is stable with respect to small perturbations or fluctuations, and there is no spinodal decomposition. However, if then the perturbation grows and the system is unstable with respect to small perturbations or fluctuations: There is spinodal decomposition.

Substituting in this concentration fluctuation, we get

This gives the same expressions for the stability as above, but it also gives an expression for the growth rate of concentration perturbations

which has a maximum at a wavevector

So, at least at the beginning of spinodal decomposition, we expect the growing concentrations to mostly have this wavevector.

Phase diagram

This type of phase transformation is known as spinodal decomposition, and can be illustrated on a phase diagram exhibiting a miscibility gap. Thus, phase separation occurs whenever a material transition into the unstable region of the phase diagram. The boundary of the unstable region sometimes referred to as the binodal or coexistence curve, is found by performing a common tangent construction of the free-energy diagram. Inside the binodal is a region called the spinodal, which is found by determining where the curvature of the free-energy curve is negative. The binodal and spinodal meet at the critical point. It is when a material is moved into the spinodal region of the phase diagram that spinodal decomposition can occur. [14]

The free energy curve is plotted as a function of composition for a temperature below the convolute temperature, T. Equilibrium phase compositions are those corresponding to the free energy minima. Regions of negative curvature (∂2f/∂c2 < 0 ) lie within the inflection points of the curve (∂2f/∂c2 = 0 ) which are called the spinodes. Their locus as a function of temperature defines the spinodal curve. For compositions within the spinodal, a homogeneous solution is unstable against infinitesimal fluctuations in density or composition, and there is no thermodynamic barrier to the growth of a new phase. Thus, the spinodal represents the limit of physical and chemical stability.

To reach the spinodal region of the phase diagram, a transition must take the material through the binodal region or the critical point. Often phase separation will occur via nucleation during this transition, and spinodal decomposition will not be observed. To observe spinodal decomposition, a very fast transition, often called a quench, is required to move from the stable to the spinodal unstable region of the phase diagram.

In some systems, ordering of the material leads to a compositional instability and this is known as a conditional spinodal, e.g. in the feldspars. [15] [16] [17] [18] [19]

Coherency strains

For most crystalline solid solutions, there is a variation of lattice parameters with the composition. If the lattice of such a solution is to remain coherent in the presence of a composition modulation, mechanical work has to be done to strain the rigid lattice structure. The maintenance of coherency thus affects the driving force for diffusion. [13] [20] [21] [22]

Consider a crystalline solid containing a one-dimensional composition modulation along the x-direction. We calculate the elastic strain energy for a cubic crystal by estimating the work required to deform a slice of material so that it can be added coherently to an existing slab of cross-sectional area. We will assume that the composition modulation is along the x' direction and, as indicated, a prime will be used to distinguish the reference axes from the standard axes of a cubic system (that is, along the <100>). [11]

Let the lattice spacing in the plane of the slab be ao and that of the undeformed slice a. If the slice is to be coherent after the addition of the slab, it must be subjected to a strain ε in the z' and y' directions which is given by:

In the first step, the slice is deformed hydrostatically in order to produce the required strains to the z' and y' directions. We use the linear compressibility of a cubic system 1 / ( c11 + 2 c12 ) where the c's are the elastic constants. The stresses required to produce a hydrostatic strain of δ are therefore given by:

The elastic work per unit volume is given by:

where the ε's are the strains. The work performed per unit volume of the slice during the first step is therefore given by:

In the second step, the sides of the slice parallel to the x' direction are clamped and the stress in this direction is relaxed reversibly. Thus, εz' = εy' = 0. The result is that:

The net work performed on the slice in order to achieve coherency is given by:

or

The final step is to express c1'1' in terms of the constants referred to the standard axes. From the rotation of axes, we obtain the following:

where l, m, n are the direction cosines of the x' axis and, therefore the direction cosines of the composition modulation. Combining these, we obtain the following:

The existence of any shear strain has not been accounted for. Cahn considered this problem, and concluded that shear would be absent for modulations along <100>, <110>, <111> and that for other directions the effect of shear strains would be small. It then follows that the total elastic strain energy of a slab of cross-sectional area A is given by:

We next have to relate the strain δ to the composition variation. Let ao be the lattice parameter of the unstrained solid of the average composition co. Using a Taylor series expansion about co yields the following:

in which

where the derivatives are evaluated at co. Thus, neglecting higher-order terms, we have:

Substituting, we obtain:

This simple result indicates that the strain energy of a composition modulation depends only on the amplitude and is independent of the wavelength. For a given amplitude, the strain energy WE is proportional to Y. Consider a few special cases.

For an isotropic material:

so that:

This equation can also be written in terms of Young's modulus E and Poisson's ratio υ using the standard relationships:

Substituting, we obtain the following:

For most metals, the left-hand side of this equation

is positive, so that the elastic energy will be a minimum for those directions that minimize the term: l2m2 + m2n2 + l2n2. By inspection, those are seen to be <100>. For this case:

the same as for an isotropic material. At least one metal (molybdenum) has an anisotropy of the opposite sign. In this case, the directions for minimum WE will be those that maximize the directional cosine function. These directions are <111>, and

As we will see, the growth rate of the modulations will be a maximum in the directions that minimize Y. These directions, therefore, determine the morphology and structural characteristics of the decomposition in cubic solid solutions.

Rewriting the diffusion equation and including the term derived for the elastic energy yields the following:

or

which can alternatively be written in terms of the diffusion coefficient D as:

The simplest way of solving this equation is by using the method of Fourier transforms.

Fourier transform

The motivation for the Fourier transformation comes from the study of a Fourier series. In the study of a Fourier series, complicated periodic functions are written as the sum of simple waves mathematically represented by sines and cosines. Due to the properties of sine and cosine, it is possible to recover the amount of each wave in the sum by an integral. In many cases it is desirable to use Euler's formula, which states that e2πiθ = cos 2πθ + i sin 2πθ, to write Fourier series in terms of the basic waves e2πiθ, with the distinct advantage of simplifying many unwieldy formulas.

The passage from sines and cosines to complex exponentials makes it necessary for the Fourier coefficients to be complex-valued. The usual interpretation of this complex number is that it gives both the amplitude (or size) of the wave present in the function and the phase (or the initial angle) of the wave. This passage also introduces the need for negative "frequencies". (E.G. If θ were measured in seconds then the waves e2πiθ and e−2πiθ would both complete one cycle per second—but they represent different frequencies in the Fourier transform. Hence, frequency no longer measures the number of cycles per unit time, but is closely related.)

If A(β) is the amplitude of a Fourier component of wavelength λ and wavenumber β = 2π/λ the spatial variation in composition can be expressed by the Fourier integral: [13]

in which the coefficients are defined by the inverse relationship:

Substituting, we obtain on equating coefficients:

This is an ordinary differential equation that has the solution:

in which A(β) is the initial amplitude of the Fourier component of wave wavenumber β and R(β) defined by:

or, expressed in terms of the diffusion coefficient D:

In a similar manner, the new diffusion equation:

has a simple sine wave solution given by:

where is obtained by substituting this solution back into the diffusion equation as follows:

For solids, the elastic strains resulting from coherency add terms to the amplification factor as follows:

where, for isotropic solids:

,

where E is Young's modulus of elasticity, ν is Poisson's ratio, and η is the linear strain per unit composition difference. For anisotropic solids, the elastic term depends on the direction in a manner that can be predicted by elastic constants and how the lattice parameters vary with composition. For the cubic case, Y is a minimum for either (100) or (111) directions, depending only on the sign of the elastic anisotropy.

Thus, by describing any composition fluctuation in terms of its Fourier components, Cahn showed that a solution would be unstable concerning to the sinusoidal fluctuations of a critical wavelength. By relating the elastic strain energy to the amplitudes of such fluctuations, he formalized the wavelength or frequency dependence of the growth of such fluctuations, and thus introduced the principle of selective amplification of Fourier components of certain wavelengths. The treatment yields the expected mean particle size or wavelength of the most rapidly growing fluctuation.

Thus, the amplitude of composition fluctuations should grow continuously until a metastable equilibrium is reached with preferential amplification of components of particular wavelengths. The kinetic amplification factor R is negative when the solution is stable to the fluctuation, zero at the critical wavelength, and positive for longer wavelengths—exhibiting a maximum at exactly times the critical wavelength.

Consider a homogeneous solution within the spinodal. It will initially have a certain amount of fluctuation from the average composition which may be written as a Fourier integral. Each Fourier component of that fluctuation will grow or diminish according to its wavelength.

Because of the maximum in R as a function of wavelength, those components of the fluctuation with times the critical wavelength will grow fastest and will dominate. This "principle of selective amplification" depends on the initial presence of these wavelengths but does not critically depend on their exact amplitude relative to other wavelengths (if the time is large compared with (1/R). It does not depend on any additional assumptions, since different wavelengths can coexist and do not interfere with one another.

Limitations of this theory would appear to arise from this assumption and the absence of an expression formulated to account for irreversible processes during phase separation which may be associated with internal friction and entropy production. In practice, frictional damping is generally present and some of the energy is transformed into thermal energy. Thus, the amplitude and intensity of a one-dimensional wave decrease with distance from the source, and for a three-dimensional wave, the decrease will be greater.

Dynamics in k-space

In the spinodal region of the phase diagram, the free energy can be lowered by allowing the components to separate, thus increasing the relative concentration of a component material in a particular region of the material. The concentration will continue to increase until the material reaches the stable part of the phase diagram. Very large regions of material will change their concentration slowly due to the amount of material that must be moved. Very small regions will shrink away due to the energy cost of maintaining an interface between two dissimilar component materials. [23] [24] [25]

To initiate a homogeneous quench a control parameter, such as temperature, is abruptly and globally changed. For a binary mixture of -type and -type materials, the Landau free-energy

is a good approximation of the free energy near the critical point and is often used to study homogeneous quenches. The mixture concentration is the density difference of the mixture components, the control parameters which determine the stability of the mixture are and , and the interfacial energy cost is determined by .

Diffusive motion often dominates at the length-scale of spinodal decomposition. The equation of motion for a diffusive system is

where is the diffusive mobility, is some random noise such that , and the chemical potential is derived from the Landau free-energy:

We see that if , small fluctuations around have a negative effective diffusive mobility and will grow rather than shrink. To understand the growth dynamics, we disregard the fluctuating currents due to , linearize the equation of motion around and perform a Fourier transform into -space. This leads to

which has an exponential growth solution:

Since the growth rate is exponential, the fastest growing angular wavenumber

will quickly dominate the morphology. We now see that spinodal decomposition results in domains of the characteristic length scale called the spinodal length:

The growth rate of the fastest-growing angular wave number is

where is known as the spinodal time.

The spinodal length and spinodal time can be used to nondimensionalize the equation of motion, resulting in universal scaling for spinodal decomposition.

Spinodal Architected Materials

Example design of a spinodal architected material. Spinodal-architecture-example.stl
Example design of a spinodal architected material.

Spinodal phase decomposition has been used to generate architected materials by interpreting one phase as solid, and the other phase as void. These spinodal architected materials present interesting mechanical properties, such as high energy absorption, [26] insensitivity to imperfections, [27] superior mechanical resilience, [28] and high stiffness-to-weight ratio. [29] Furthermore, by controlling the phase separation, i.e., controlling the proportion of materials, and/or imposing preferential directions in the decompositions, one can control the density, and preferential directions effectively tuning the strength, weight, and anisotropy of the resulting architected material. [30] Another interesting property of spinodal materials is the capability to seamlessly transition between different classes, orientations, and densities, [30] thereby enabling the manufacturing of effectively multi-material structures. [31]

Related Research Articles

In the general theory of relativity, the Einstein field equations relate the geometry of spacetime to the distribution of matter within it.

In physics, the S-matrix or scattering matrix relates the initial state and the final state of a physical system undergoing a scattering process. It is used in quantum mechanics, scattering theory and quantum field theory (QFT).

In physics, a wave vector is a vector used in describing a wave, with a typical unit being cycle per metre. It has a magnitude and direction. Its magnitude is the wavenumber of the wave, and its direction is perpendicular to the wavefront. In isotropic media, this is also the direction of wave propagation.

In differential geometry, the four-gradient is the four-vector analogue of the gradient from vector calculus.

In differential geometry, a tensor density or relative tensor is a generalization of the tensor field concept. A tensor density transforms as a tensor field when passing from one coordinate system to another, except that it is additionally multiplied or weighted by a power W of the Jacobian determinant of the coordinate transition function or its absolute value. A tensor density with a single index is called a vector density. A distinction is made among (authentic) tensor densities, pseudotensor densities, even tensor densities and odd tensor densities. Sometimes tensor densities with a negative weight W are called tensor capacity. A tensor density can also be regarded as a section of the tensor product of a tensor bundle with a density bundle.

<span class="mw-page-title-main">Voigt effect</span>

The Voigt effect is a magneto-optical phenomenon which rotates and elliptizes linearly polarised light sent into an optically active medium. The effect is named after the German scientist Woldemar Voigt who discovered it in vapors. Unlike many other magneto-optical effects such as the Kerr or Faraday effect which are linearly proportional to the magnetization, the Voigt effect is proportional to the square of the magnetization and can be seen experimentally at normal incidence. There are also other denominations for this effect, used interchangeably in the modern scientific literature: the Cotton–Mouton effect and magnetic-linear birefringence, with the latter reflecting the physical meaning of the effect.

The Newman–Penrose (NP) formalism is a set of notation developed by Ezra T. Newman and Roger Penrose for general relativity (GR). Their notation is an effort to treat general relativity in terms of spinor notation, which introduces complex forms of the usual variables used in GR. The NP formalism is itself a special case of the tetrad formalism, where the tensors of the theory are projected onto a complete vector basis at each point in spacetime. Usually this vector basis is chosen to reflect some symmetry of the spacetime, leading to simplified expressions for physical observables. In the case of the NP formalism, the vector basis chosen is a null tetrad: a set of four null vectors—two real, and a complex-conjugate pair. The two real members often asymptotically point radially inward and radially outward, and the formalism is well adapted to treatment of the propagation of radiation in curved spacetime. The Weyl scalars, derived from the Weyl tensor, are often used. In particular, it can be shown that one of these scalars— in the appropriate frame—encodes the outgoing gravitational radiation of an asymptotically flat system.

<span class="mw-page-title-main">Gravitational lensing formalism</span>

In general relativity, a point mass deflects a light ray with impact parameter by an angle approximately equal to

A synchronous frame is a reference frame in which the time coordinate defines proper time for all co-moving observers. It is built by choosing some constant time hypersurface as an origin, such that has in every point a normal along the time line and a light cone with an apex in that point can be constructed; all interval elements on this hypersurface are space-like. A family of geodesics normal to this hypersurface are drawn and defined as the time coordinates with a beginning at the hypersurface. In terms of metric-tensor components , a synchronous frame is defined such that

Equatorial Rossby waves, often called planetary waves, are very long, low frequency water waves found near the equator and are derived using the equatorial beta plane approximation.

In quantum field theory, a non-topological soliton (NTS) is a soliton field configuration possessing, contrary to a topological one, a conserved Noether charge and stable against transformation into usual particles of this field for the following reason. For fixed charge Q, the mass sum of Q free particles exceeds the energy (mass) of the NTS so that the latter is energetically favorable to exist.

<span class="mw-page-title-main">Plate theory</span> Mathematical model of the stresses within flat plates under loading

In continuum mechanics, plate theories are mathematical descriptions of the mechanics of flat plates that draw on the theory of beams. Plates are defined as plane structural elements with a small thickness compared to the planar dimensions. The typical thickness to width ratio of a plate structure is less than 0.1. A plate theory takes advantage of this disparity in length scale to reduce the full three-dimensional solid mechanics problem to a two-dimensional problem. The aim of plate theory is to calculate the deformation and stresses in a plate subjected to loads.

<span class="mw-page-title-main">Uflyand-Mindlin plate theory</span>

The Reissner–Mindlin theory of plates is an extension of Kirchhoff–Love plate theory that takes into account shear deformations through-the-thickness of a plate. The theory was proposed in 1951 by Raymond Mindlin. A similar, but not identical, theory in static setting, had been proposed earlier by Eric Reissner in 1945. Both theories are intended for thick plates in which the normal to the mid-surface remains straight but not necessarily perpendicular to the mid-surface. The Reissner-Mindlin theory is used to calculate the deformations and stresses in a plate whose thickness is of the order of one tenth the planar dimensions while the Kirchhoff–Love theory is applicable to thinner plates.

In accelerator physics, the term acceleration voltage means the effective voltage surpassed by a charged particle along a defined straight line. If not specified further, the term is likely to refer to the longitudinal effective acceleration voltage.

Q-vectors are used in atmospheric dynamics to understand physical processes such as vertical motion and frontogenesis. Q-vectors are not physical quantities that can be measured in the atmosphere but are derived from the quasi-geostrophic equations and can be used in the previous diagnostic situations. On meteorological charts, Q-vectors point toward upward motion and away from downward motion. Q-vectors are an alternative to the omega equation for diagnosing vertical motion in the quasi-geostrophic equations.

Heat transfer physics describes the kinetics of energy storage, transport, and energy transformation by principal energy carriers: phonons, electrons, fluid particles, and photons. Heat is thermal energy stored in temperature-dependent motion of particles including electrons, atomic nuclei, individual atoms, and molecules. Heat is transferred to and from matter by the principal energy carriers. The state of energy stored within matter, or transported by the carriers, is described by a combination of classical and quantum statistical mechanics. The energy is different made (converted) among various carriers. The heat transfer processes are governed by the rates at which various related physical phenomena occur, such as the rate of particle collisions in classical mechanics. These various states and kinetics determine the heat transfer, i.e., the net rate of energy storage or transport. Governing these process from the atomic level to macroscale are the laws of thermodynamics, including conservation of energy.

Vasiliev equations are formally consistent gauge invariant nonlinear equations whose linearization over a specific vacuum solution describes free massless higher-spin fields on anti-de Sitter space. The Vasiliev equations are classical equations and no Lagrangian is known that starts from canonical two-derivative Frønsdal Lagrangian and is completed by interactions terms. There is a number of variations of Vasiliev equations that work in three, four and arbitrary number of space-time dimensions. Vasiliev's equations admit supersymmetric extensions with any number of super-symmetries and allow for Yang–Mills gaugings. Vasiliev's equations are background independent, the simplest exact solution being anti-de Sitter space. It is important to note that locality is not properly implemented and the equations give a solution of certain formal deformation procedure, which is difficult to map to field theory language. The higher-spin AdS/CFT correspondence is reviewed in Higher-spin theory article.

<span class="mw-page-title-main">Phonovoltaic</span>

A phonovoltaic (pV) cell converts vibrational (phonons) energy into a direct current much like the photovoltaic effect in a photovoltaic (PV) cell converts light (photon) into power. That is, it uses a p-n junction to separate the electrons and holes generated as valence electrons absorb optical phonons more energetic than the band gap, and then collects them in the metallic contacts for use in a circuit. The pV cell is an application of heat transfer physics and competes with other thermal energy harvesting devices like the thermoelectric generator.

In statistical mechanics, Lee–Yang theory, sometimes also known as Yang–Lee theory, is a scientific theory which seeks to describe phase transitions in large physical systems in the thermodynamic limit based on the properties of small, finite-size systems. The theory revolves around the complex zeros of partition functions of finite-size systems and how these may reveal the existence of phase transitions in the thermodynamic limit.

In physical oceanography and fluid mechanics, the Miles-Phillips mechanism describes the generation of wind waves from a flat sea surface by two distinct mechanisms. Wind blowing over the surface generates tiny wavelets. These wavelets develop over time and become ocean surface waves by absorbing the energy transferred from the wind. The Miles-Phillips mechanism is a physical interpretation of these wind-generated surface waves.
Both mechanisms are applied to gravity-capillary waves and have in common that waves are generated by a resonance phenomenon. The Miles mechanism is based on the hypothesis that waves arise as an instability of the sea-atmosphere system. The Phillips mechanism assumes that turbulent eddies in the atmospheric boundary layer induce pressure fluctuations at the sea surface. The Phillips mechanism is generally assumed to be important in the first stages of wave growth, whereas the Miles mechanism is important in later stages where the wave growth becomes exponential in time.

References

  1. Binder, K (1987-07-01). "Theory of first-order phase transitions". Reports on Progress in Physics. 50 (7): 783–859. doi:10.1088/0034-4885/50/7/001. ISSN   0034-4885. S2CID   250879624.
  2. Gennes, Pierre-Gilles de. (1979). Scaling concepts in polymer physics. Ithaca, N.Y.: Cornell University Press. ISBN   0-8014-1203-X. OCLC   4494721.
  3. Gibbs, J.W., Scientific Papers of J Willard Gibbs, 2 vols. Bumstead, H. A., and Van Name, R. G., eds. (Dover, New York, 1961) ISBN   0-918024-77-3
  4. Hillert, M., A Theory of Nucleation for Solid Metallic Solutions, Sc. D. Thesis (MIT, 1955)
  5. Hillert, M (1961). "A solid-solution model for inhomogeneous systems". Acta Metallurgica. 9 (6). Elsevier BV: 525–535. doi:10.1016/0001-6160(61)90155-9. ISSN   0001-6160.
  6. Cahn, John W (1961). "On spinodal decomposition". Acta Metallurgica. 9 (9). Elsevier BV: 795–801. doi:10.1016/0001-6160(61)90182-1. ISSN   0001-6160.
  7. Cahn, John W (1962). "On spinodal decomposition in cubic crystals". Acta Metallurgica. 10 (3). Elsevier BV: 179–183. doi:10.1016/0001-6160(62)90114-1. ISSN   0001-6160.
  8. Cahn, John W (1962). "Coherent fluctuations and nucleation in isotropic solids". Acta Metallurgica. 10 (10). Elsevier BV: 907–913. doi:10.1016/0001-6160(62)90140-2. ISSN   0001-6160.
  9. Cahn, John W.; Hilliard, John E. (1958). "Free Energy of a Nonuniform System. I. Interfacial Free Energy". The Journal of Chemical Physics. 28 (2). AIP Publishing: 258–267. Bibcode:1958JChPh..28..258C. doi:10.1063/1.1744102. ISSN   0021-9606.
  10. 1 2 Bray, A. J. (2002-03-01). "Theory of phase-ordering kinetics". Advances in Physics. 51 (2): 481–587. arXiv: cond-mat/9501089 . Bibcode:2002AdPhy..51..481B. doi:10.1080/00018730110117433. ISSN   0001-8732. S2CID   218646292.
  11. 1 2 Hilliard, J.E., Spinodal Decomposition, in Phase Transformations p. 497 (American Society of Metals, Metals Park, 1970)
  12. Bray, A. J. (1994). "Theory of phase ordering kinetics". Physica A: Statistical Mechanics and Its Applications. 194 (1): 41–52. arXiv: cond-mat/9501089 . doi:10.1016/0378-4371(93)90338-5. ISSN   0378-4371.
  13. 1 2 3 Cahn, J.W., Spinodal Decomposition, 1967 Institute of Metals Lecture, Trans. Met. Soc. AIME, Vol. 242, p. 168 (1968)
  14. Jones, Richard A. L. (2004) [2002]. Soft Condensed Matter. Oxford University Press. p. 33. ISBN   978-0-19-850589-1 . Retrieved 2007-10-22.
  15. Cook, H.E (1973). "A lattice model of structural and dislocation transformations". Acta Metallurgica. 21 (10). Elsevier BV: 1431–1444. doi:10.1016/0001-6160(73)90092-8. ISSN   0001-6160.
  16. Cook, H.E (1973). "On the nature of the omega transformation". Acta Metallurgica. 21 (10). Elsevier BV: 1445–1449. doi:10.1016/0001-6160(73)90093-x. ISSN   0001-6160.
  17. Cook, H.E (1975). "On first-order structural phase transitions—I. General considerations of pre-transition and nucleation phenomena". Acta Metallurgica. 23 (9). Elsevier BV: 1027–1039. doi:10.1016/0001-6160(75)90107-8. ISSN   0001-6160.
  18. Suzuki, T . and Wuttig, M., Analogy between spinodal decomposition and martensitic transformation, Acta Met., Vol. 23, p.1069 (1975)
  19. Carpenter, M. A. (1981). "A "conditional spinodal" within the peristerite miscibility gap of plagioclase feldspars" (PDF). American Mineralogist. 66: 553–560.
  20. de Fontaine, D (1969). "An approximate criterion for the loss of coherency in modulated structures". Acta Metallurgica. 17 (4). Elsevier BV: 477–482. doi:10.1016/0001-6160(69)90029-7. ISSN   0001-6160.
  21. Cook, H.E; De Fontaine, D; Hilliard, J.e (1969). "A model for diffusion on cubic lattices and its application to the early stages of ordering". Acta Metallurgica. 17 (6). Elsevier BV: 765–773. doi:10.1016/0001-6160(69)90083-2. ISSN   0001-6160.
  22. Cook, H.E; de Fontaine, D (1969). "On the elastic free energy of solid solutions—I. Microscopic theory". Acta Metallurgica. 17 (7). Elsevier BV: 915–924. doi:10.1016/0001-6160(69)90112-6. ISSN   0001-6160.
  23. De Fontaine, D. (1970). "Mechanical instabilities in the b.c.c. lattice and the beta to omega phase transformation". Acta Metallurgica. 18 (2). Elsevier BV: 275–279. doi:10.1016/0001-6160(70)90035-0. ISSN   0001-6160.
  24. Cook, H.E.; De Fontaine, D. (1971). "On the elastic free energy of solid solutions—II. Influence of the effective modulus on precipitation from solution and the order-disorder reaction". Acta Metallurgica. 19 (7). Elsevier BV: 607–616. doi:10.1016/0001-6160(71)90013-7. ISSN   0001-6160.
  25. De Fontaine, D; Paton, N.E; Williams, J.C (1971). "The omega phase transformation in titanium alloys as an example of displacement controlled reactions". Acta Metallurgica. 19 (11). Elsevier BV: 1153–1162. doi:10.1016/0001-6160(71)90047-2. ISSN   0001-6160.
  26. Guell Izard, Anna; Bauer, Jens; Crook, Cameron; Turlo, Vladyslav; Valdevit, Lorenzo (November 2019). "Ultrahigh Energy Absorption Multifunctional Spinodal Nanoarchitectures". Small. 15 (45): 1903834. doi:10.1002/smll.201903834. PMID   31531942. S2CID   202672680.
  27. Hsieh, Meng-Ting; Endo, Bianca; Zhang, Yunfei; Bauer, Jens; Valdevit, Lorenzo (1 April 2019). "The mechanical response of cellular materials with spinodal topologies". Journal of the Mechanics and Physics of Solids. 125: 401–419. arXiv: 1904.06733 . Bibcode:2019JMPSo.125..401H. doi:10.1016/j.jmps.2019.01.002. S2CID   118645530.
  28. Portela, Carlos M.; Vidyasagar, A.; Krödel, Sebastian; Weissenbach, Tamara; Yee, Daryl W.; Greer, Julia R.; Kochmann, Dennis M. (2020-03-17). "Extreme mechanical resilience of self-assembled nanolabyrinthine materials". Proceedings of the National Academy of Sciences. 117 (11): 5686–5693. Bibcode:2020PNAS..117.5686P. doi: 10.1073/pnas.1916817117 . ISSN   0027-8424. PMC   7084143 . PMID   32132212.
  29. Zheng, Li; Kumar, Siddhant; Kochmann, Dennis M. (September 2021). "Data-driven topology optimization of spinodoid metamaterials with seamlessly tunable anisotropy". Computer Methods in Applied Mechanics and Engineering. 383: 113894. arXiv: 2012.15744 . Bibcode:2021CMAME.383k3894Z. doi:10.1016/j.cma.2021.113894. S2CID   229923529.
  30. 1 2 Kumar, Siddhant; Tan, Stephanie; Zheng, Li; Kochmann, Dennis M. (December 2020). "Inverse-designed spinodoid metamaterials". npj Computational Materials. 6 (1): 73. Bibcode:2020npjCM...6...73K. doi: 10.1038/s41524-020-0341-6 . hdl: 20.500.11850/418523 . S2CID   219315003.
  31. Senhora, Fernando V.; Sanders, Emily D.; Paulino, Glaucio H. (27 April 2022). "Optimally-Tailored Spinodal Architected Materials for Multiscale Design and Manufacturing". Advanced Materials. 34 (26): 2109304. Bibcode:2022AdM....3409304S. doi: 10.1002/adma.202109304 . PMID   35297113. S2CID   247498613.

Further reading