Separation of variables
Navier–Stokes differential equations used to simulate airflow around an obstruction.
In mathematics, separation of variables (also known as the Fourier method) is any of several methods for solving ordinary and partial differential equations, in which algebra allows one to rewrite an equation so that each of two variables occurs on a different side of the equation.
Ordinary differential equations (ODE)
Suppose a differential equation can be written in the form
which we can write more simply by letting :
As long as h(y) ≠ 0, we can rearrange terms to obtain:
so that the two variables x and y have been separated. dx (and dy) can be viewed, at a simple level, as just a convenient notation, which provides a handy mnemonic aid for assisting with manipulations. A formal definition of dx as a differential (infinitesimal) is somewhat advanced.
Some who dislike Leibniz's notation may prefer to write this as
but that fails to make it quite as obvious why this is called "separation of variables". Integrating both sides of the equation with respect to , we have
because of the substitution rule for integrals.
If one can evaluate the two integrals, one can find a solution to the differential equation. Observe that this process effectively allows us to treat the derivative as a fraction which can be separated. This allows us to solve separable differential equations more conveniently, as demonstrated in the example below.
(Note that we do not need to use two constants of integration, in equation (1) as in
because a single constant is equivalent.)
Population growth is often modeled by the differential equation
where is the population with respect to time , is the rate of growth, and is the carrying capacity of the environment.
Separation of variables may be used to solve this differential equation.
To evaluate the integral on the left side, we simplify the fraction
and then, we decompose the fraction into partial fractions
Thus we have
- Let .
Therefore, the solution to the logistic equation is
To find , let and . Then we have
Noting that , and solving for A we get
Partial differential equations
The method of separation of variables is also used to solve a wide range of linear partial differential equations with boundary and initial conditions, such as heat equation, wave equation, Laplace equation and Helmholtz equation.
Consider the one-dimensional heat equation.The equation is
The boundary condition is homogeneous, that is
Let us attempt to find a solution which is not identically zero satisfying the boundary conditions but with the following property: u is a product in which the dependence of u on x, t is separated, that is:
Substituting u back into equation and using the product rule,
Since the right hand side depends only on x and the left hand side only on t, both sides are equal to some constant value − λ. Thus:
We will now show that solutions for X(x) for values of λ ≤ 0 cannot occur:
Suppose that λ < 0. Then there exist real numbers B, C such that
From (2) we get
and therefore B = 0 = C which implies u is identically 0.
Suppose that λ = 0. Then there exist real numbers B, C such that
From (7) we conclude in the same manner as in 1 that u is identically 0.
Therefore, it must be the case that λ > 0. Then there exist real numbers A, B, C such that
From (7) we get C = 0 and that for some positive integer n,
This solves the heat equation in the special case that the dependence of u has the special form of (3).
where Dn are coefficients determined by initial condition.
Given the initial condition
we can get
This is the sine series expansion of f(x). Multiplying both sides with and integrating over [0,L] result in
Suppose the equation is nonhomogeneous,
with the boundary condition the same as (2).
Expand h(x,t), u(x,t) and f(x,t) into
where hn(t) and bn can be calculated by integration, while un(t) is to be determined.
If the boundary condition is nonhomogeneous, then the expansion of (9) and (10) is no longer valid. One has to find a function v that satisfies the boundary condition only, and subtract it from u. The function u-v then satisfies homogeneous boundary condition, and can be solved with the above method.
In orthogonal curvilinear coordinates, separation of variables can still be used, but in some details different from that in Cartesian coordinates. For instance, regularity or periodic condition may determine the eigenvalues in place of boundary conditions. See spherical harmonics for example.
The matrix form of the separation of variables is the Kronecker sum.
where and are 1D discrete Laplacians in the x- and y-directions, correspondingly, and are the identities of appropriate sizes. See the main article Kronecker sum of discrete Laplacians for details.
- A. D. Polyanin, Handbook of Linear Partial Differential Equations for Engineers and Scientists, Chapman & Hall/CRC Press, Boca Raton, 2002. ISBN 1-58488-299-9.
- Tyn Myint-U, Lokenath Debnath (2007). Linear Partial Differential Equations for Scientists and Engineers. Boston, MA. ISBN 978-0-8176-4393-5. Retrieved 2011-03-29.
- Teschl, Gerald (2012). Ordinary Differential Equations and Dynamical Systems. Graduate Studies in Mathematics 140. Providence, RI: American Mathematical Society. ISBN 978-0-8218-8328-0.
- Hazewinkel, Michiel, ed. (2001), "Fourier method", Encyclopedia of Mathematics, Springer, ISBN 978-1-55608-010-4
- John Renze, Eric W. Weisstein, Separation of variables (Differential Equation) at MathWorld
- Methods of Generalized and Functional Separation of Variables at EqWorld: The World of Mathematical Equations
- Examples of separating variables to solve PDEs
- "A Short Justification of Separation of Variables"