Method of dominant balance

Short description: Solution of a simplified form of an equation

In mathematics, the method of dominant balance approximates the solution to an equation by solving a simplified form of the equation containing 2 or more of the equation's terms that most influence (dominate) the solution and excluding terms contributing only small modifications to this approximate solution. Following an initial solution, iteration of the procedure may generate additional terms of an asymptotic expansion providing a more accurate solution.

An early example of the dominant balance method is the Newton polygon method. Newton developed this method to find an explicit approximation for an algebraic function. Newton expressed the function as proportional to the independent variable raised to a power, retained only the lowest-degree polynomial terms (dominant terms), and solved this simplified reduced equation to obtain an approximate solution. Dominant balance has a broad range of applications, solving differential equations arising in fluid mechanics, plasma physics, turbulence, combustion, nonlinear optics, geophysical fluid dynamics, and neuroscience.

Asymptotic relations

The functions $f(z)$ and $g(z)$ of parameter or independent variable $z$ and the quotient $f(z)/g(z)$ have limits as $z$ approaches the limit $L$.

The function $f(z)$ is much less than $g(z)$ as $z$ approaches $L$, written as $f(z)\ll g(z)(z\to L)$, if the limit of the quotient $f(z)/g(z)$ is zero as $z$ approaches $L$.

The relation $f(z)$ is lower order than $g(z)$ as $z$ approaches $L$, written using little-o notation $f(z)=o(g(z))(z\to L)$, is identical to the $f(z)$ is much less than $g(z)$ as $z$ approaches $L$ relation.

The function $f(z)$ is equivalent to $g(z)$ as $z$ approaches $L$, written as $f(z)\sim g(z)(z\to L)$, if the limit of the quotient $f(z)/g(z)$ is 1 as $z$ approaches $L$.

This result indicates that the zero function, $f(z)=0$ for all values of $z$, can never be equivalent to any other function.

Asymptotically equivalent functions remain asymptotically equivalent under integration if requirements related to convergence are met. There are more specific requirements for asymptotically equivalent functions to remain asymptotically equivalent under differentiation.

Equation properties

An equation's approximate solution is $s(z)$ as $z$ approaches limit $L$. The equation's terms that may be constants or contain this solution are $T_0(s),T_1(s),\dots ,T_n(s)$. If the approximate solution is fully correct, the equation's terms sum to zero in this equation: $$T_0(s)+T_1(s)+\dots +T_n(s)=0.$$ For distinct integer indices $i,j$, this equation is a sum of 2 terms and a remainder $R_{ij}(s)$ expressed as $$\begin{array}{rl}& T_i(s)+T_j(s)+R_{ij}(s)=0\\ & R_{ij}(s)=\sum _{\genfrac{}{}{0ex}{}{k=0}{k\ne i,k\ne j}}^n{T}_k(s).\end{array}$$

Balance equation terms $T_i(s)$ and $T_j(s)$ means find the function $s(z)$ that solves the reduced equation $T_i(s)+T_j(s)=0$ with $T_i(s)\ne 0$ and $T_j(s)\ne 0$. This solution $s(z)$ is consistent if terms $T_i(s)$ and $T_j(s)$ are dominant; dominant means the remaining equation terms $R_{ij}(s)$ are much less than terms $T_i(s)$ and $T_j(s)$ as $z$ approaches $L$. A consistent solution that balances two equation terms may generate an accurate approximation to the full equation's solution for $z$ values approaching $L$. Approximate solutions arising from balancing different terms of an equation may generate distinct approximate solutions e.g. inner and outer layer solutions.

In some cases, the reduced equation has coefficients containing a parameter $ϵ$, but the parameter's exponent is different for the 2 terms. In this case, the reduced equation requires a scale transformation $z={ϵ}^{\lambda }\tilde{z}$ of the variable $z$ using scaling exponent $\lambda$. The dominant balance method selects a scaling exponent $\lambda$ to generate transformed terms, each transformed term now containing the parameter with the same exponent. Distinct pairs of terms and scaling transformations lead to distinct simplified equations called distinguished limits, and each equation may lead to a distinct solution. The Kruskal-Newton diagram facilitates identifying the required scale transformation needed for dominant balance of algebraic and differential equations.

For differential equation solutions containing an irregular singularity, the leading behavior is the first term of an asymptotic series solution that remains when the independent variable $z$ approaches an irregular singularity $L$. The controlling factor is the fastest changing part of the leading behavior. It is advised to "show that the equation for the function obtained by factoring off the dominant balance solution from the exact solution itself has a solution that varies less rapidly than the dominant balance solution."

Algorithm

The input is the set of equation terms and the limit L. The output is the set of approximate solutions. For each pair of distinct equation terms $T_i(s),T_j(s)$ the algorithm applies a scale transformation if needed, balances the selected terms by finding a function that solves the reduced equation and then determines if this function is consistent. If the function balances the terms and is consistent, the algorithm adds the function to the set of approximate solutions, otherwise the algorithm rejects the function. The process is repeated for each pair of distinct equation terms.

Inputs Set of equation terms $\{T_0(s),T_1(s),\dots ,T_n(s)\}$ and limit $L$

Output Set of approximate solutions $\{s_0(z),s_1(z),\dots \}$

1. For each pair of distinct equation terms $T_i(s),T_j(s)$ do:
   1. Apply a scale transformation if needed.
   2. Solve the reduced equation: $T_i(s)+T_j(s)=0$ with $T_i(s)\ne 0$ and $T_j(s)\ne 0$.
   3. Verify consistency: $R_{ij}(s)\ll T_i(s)(z\to L)$ and $R_{ij}(s)\ll T_j(s)(z\to L).$
   4. If function $s(z)$ is consistent and solves the reduced equation, add this function to the set of approximate solutions, otherwise reject the function.

Improved accuracy

The method may be iterated to generate additional terms of an asymptotic expansion to provide a more accurate solution. Iterative methods such as the Newton-Raphson method may generate a more accurate solution. A perturbation series, using the approximate solution as the first term, may also generate a more accurate solution.

Examples

Algebraic function

The dominant balance method will generate an approximate solution for variable $s$ in algebraic equation $1-16s+ϵs^5=0$ that depends on parameter $ϵ$. This example computes asymptotic values as $ϵ$ approaches zero because this example's goal is finding solutions for small $ϵ$ parameter values. The function $s=s(ϵ)$ denotes how the equation's solution depends on the equation's parameter $ϵ$. The solution set is enumerated as $\{s_0(ϵ),s_1(ϵ),s_2(ϵ),s_3(ϵ),s_4(ϵ)\}$.

Input

The set of equation terms is $\{1,-16s,ϵs^5\}$ and asymptotic computations are as $ϵ$ approaches zero.

First term pair

1. Select the terms $1$ and $-16s$.
2. Rewriting the selected terms as ${ϵ}^{0}1$ and $-{ϵ}^{0}16s$, the parameter's exponents are the same (0,0) for both terms so no scale transformation needed.
3. Solve the reduced equation: $1-16s=0,s(ϵ)=\frac{1}{16}$.
4. Verify consistency: $ϵs^5\ll 1(ϵ\to 0),ϵs^5\ll 16s(ϵ\to 0)$ for $s(ϵ)=\frac{1}{16}.$
5. Add this function to the set of approximate solutions: $s_0(ϵ)=\frac{1}{16}$.

Second term pair

1. Select the terms $-16s$ and $ϵs^5$.
2. Rewriting the selected terms as $-{ϵ}^{0}16s$ and ${ϵ}^1{s}^5$, the parameter's exponents are different (0,1) so a scale transformation is needed. Apply the scale transformation $s={ϵ}^{-1/4}\tilde{s}$. The parameter exponents of the transformed terms $-{ϵ}^{-1/4}16\tilde{s}$ and ${ϵ}^{-1/4}{\tilde{s}}^5=0$ are the same (-1/4,-1/4). The transformed equation becomes ${ϵ}^{1/4}-16\tilde{s}+{\tilde{s}}^5=0$.
3. Solve the reduced equation: $-16\tilde{s}+{\tilde{s}}^5=0,\tilde{s}=2,-2,2i,-2i$.
4. Verify consistency: ${ϵ}^{1/4}\ll 16\tilde{s}(ϵ\to 0),{ϵ}^{1/4}\ll {\tilde{s}}^5(ϵ\to 0)$ for $\tilde{s}=2,-2,2i,-2i.$
5. Add these functions to the set of approximate solutions:

$$s_1(ϵ)=\frac{2}{{ϵ}^{1/4}},s_2(ϵ)=\frac{-2}{{ϵ}^{1/4}},s_3(ϵ)=\frac{2i}{{ϵ}^{1/4}},s_4(ϵ)=\frac{-2i}{{ϵ}^{1/4}}.$$

Third term pair

1. Select the terms $1$ and $ϵs^5$.
2. Rewriting the selected terms as ${ϵ}^{0}1$ and ${ϵ}^1{s}^5$, the parameter's exponents are different (0,1) so a scale transformation is needed. Apply the scale transformation $s={ϵ}^{-1/5}\tilde{s}$. The parameter exponents of the transformed terms ${ϵ}^{0}1$ and ${ϵ}^0{\tilde{s}}^5$ are the same (0,0). The transformed equation is $1-16{ϵ}^{-1/5}\tilde{s}+{\tilde{s}}^5=0.$
3. Solve the reduced equation: $1+{\tilde{s}}^5=0,\tilde{s}=(-1{)}^{1/5}.$
4. The function is not consistent: $-16{ϵ}^{-1/5}\tilde{s}\gg 1(ϵ\to 0),{ϵ}^{-1/5}\tilde{s}\gg {\tilde{s}}^5(ϵ\to 0)$ for $\tilde{s}=(-1{)}^{1/5}.$
5. Reject this function: $s={ϵ}^{-1/5}(-1{)}^{1/5}.$

Output

The set of approximate solutions has 5 functions: $$\{\frac{1}{16},\frac{2}{{ϵ}^{1/4}},\frac{-2}{{ϵ}^{1/4}},\frac{2i}{{ϵ}^{1/4}},\frac{-2i}{{ϵ}^{1/4}}\}.$$

Perturbation series solution

The approximate solutions are the first terms in the perturbation series solutions.

$$\begin{array}{rl}& s_0(ϵ)=\frac{1}{16}+\frac{1}{16777216}{ϵ}^1+\frac{5}{17592186044416}{ϵ}^2+\dots ,\\ & s_1(ϵ)=\frac{2}{{ϵ}^{1/4}}-\frac{1}{64}-\frac{5}{16384}{ϵ}^{\frac{1}{4}}-\frac{5}{524288}{ϵ}^{\frac{1}{2}}-\dots ,\\ & s_2(ϵ)=-\frac{2}{{ϵ}^{1/4}}-\frac{1}{64}+\frac{5}{16384}{ϵ}^{\frac{1}{4}}-\frac{5}{524288}{ϵ}^{\frac{1}{2}}+\dots ,\\ & s_3(ϵ)=\frac{2i}{{ϵ}^{1/4}}-\frac{1}{64}+\frac{5i}{16384}{ϵ}^{\frac{1}{4}}+\frac{5}{524288}{ϵ}^{\frac{1}{2}}-\dots \\ & s_4(ϵ)=-\frac{2i}{{ϵ}^{1/4}}-\frac{1}{64}-\frac{5i}{16384}{ϵ}^{\frac{1}{4}}+\frac{5}{524288}{ϵ}^{\frac{1}{2}}+\dots ,\\ \end{array}$$

Differential equation

The differential equation $z^3{w}^{\prime \prime }-w=0$ is known to have a solution with an exponential leading term. The transformation $w(z)=e^{s(z)}$ leads to the differential equation $1-z^3(s^{\prime }{)}^2-z^3{s}^{\prime \prime }=0$. The dominant balance method will find an approximate solution as $z$ approaches zero because this example's focus is for solutions for small $z$ values. Scaled transformations will not be used because there are no equation parameters. Two consecutive applications of the dominant balance method will generate the 2 leading terms of an asymptotic series.

Input

The set of equation terms is $\{1,-z^3(s^{\prime }{)}^2,-z^3{s}^{\prime \prime }\}$ and the limit is zero.

First term pair

1. Select $1$ and $-z^3(s^{\prime }{)}^2$.
2. The scale transformation is not required.
3. Solve the reduced equation: $1-z^3(s^{\prime }{)}^2=0,s(z)=\pm 2z^{-1/2}$
4. Verify consistency: $z^3{s}^{\prime \prime }\ll 1(z\to 0),z^3{s}^{\prime \prime }\ll z^3(s^{\prime }{)}^2(z\to 0)$ for $s(z)=\pm 2z^{-1/2}.$
5. Add these 2 functions to the set of approximate solutions: $s_{+}(z)=+2z^{-1/2},s_{-}(z)=-2z^{-1/2}.$

Second term pair

1. Select $1$ and $-z^3{s}^{\prime \prime }$
2. The scale transformation is not required.
3. Solve the reduced equation: $1-z^3{s}^{\prime \prime }=0,s(z)=\frac{1}{2}{z}^{-1}$
4. The function is not consistent: $z^3(s^{\prime }{)}^2\gg 1(z\to 0),z^3(s^{\prime }{)}^2\gg z^3{s}^{\prime \prime }(z\to 0)$ for $s(z)=\frac{1}{2}{z}^{-1}.$
5. Reject this function: $s(z)=\frac{1}{2}{z}^{-1}.$.

Third term pair

1. Select $-z^3(s^{\prime }{)}^2$ and $-z^3{s}^{\prime \prime }$.
2. The scale transformation is not required.
3. Solve the reduced equation: $z^3(s^{\prime }{)}^2+z^3{s}^{\prime \prime }=0,s(z)=\ln z$.
4. The function is not consistent: $1\gg z^3(s^{\prime }{)}^2(z\to 0)$ and $1\gg z^3{s}^{\prime \prime }(z\to 0)$ for $s(z)=\ln z.$
5. Reject this function: $s(z)=\ln z.$

Output

The set of approximate solutions has 2 functions: $$\{+2z^{-1/2},-2z^{-1/2}\}.$$

Find 2-term solutions

Using the 1-term solution, a 2-term solution is $$s_{2\pm }(z)=\pm 2z^{-1/2}+s(z).$$ Substitution of this 2-term solution into the original differential equation generates a new differential equation: $$\begin{array}{rl}1-z^3(s_{2\pm }^{\prime }{)}^2-z^3{s}_{2\pm }^{\prime \prime }& =0\\ \pm 1\mp \frac{4}{3}z{s}^{\prime }+\frac{2}{3}{z}^{5/2}(s^{\prime }{)}^2+\frac{2}{3}{z}^{5/2}{s}^{\prime \prime }& =0.\end{array}$$

Input

The set of equation terms is $\{\pm 1,\mp \frac{4}{3}z{s}^{\prime },\frac{2}{3}{z}^{5/2}(s^{\prime }{)}^2,\frac{2}{3}{z}^{5/2}{s}^{\prime \prime }\}$ and the limit is zero.

First term pair

1. Select $1$ and $-\frac{4}{3}z{s}^{\prime }$.

2. The scale transformation is not required.

3. Solve the reduced equation: $1-\frac{4}{3}z{s}^{\prime }=0,s(z)=\frac{3}{4}\ln z$.

4. Verify consistency:

$\frac{2}{3}{z}^{5/2}(s^{\prime }{)}^2+\frac{2}{3}{z}^{5/2}{s}^{\prime \prime }\ll 1(z\to 0),\text{for}s(z)=\frac{3}{4}\ln z$

$\frac{2}{3}{z}^{5/2}(s^{\prime }{)}^2+\frac{2}{3}{z}^{5/2}{s}^{\prime \prime }\ll \frac{4}{3}z{s}^{\prime }(z\to 0)\text{for}s(z)=\frac{3}{4}\ln z.$

5. Add these functions to the set of approximate solutions:

$s_{2+}(z)=+2z^{-1/2}+\frac{3}{4}\ln z$

$s_{2-}(z)=-2z^{-1/2}+\frac{3}{4}\ln z$.

Other term pairs

For other term pairs, the functions that solve the reduced equations are not consistent.

Output

The set of approximate solutions has 2 functions: $$\{+2z^{-1/2}+\frac{3}{4}\ln z,-2z^{-1/2}+\frac{3}{4}\ln z\}.$$

Asymptotic expansion

The next iteration generates a 3-term solution $s_{3\pm }(z)=\pm 2z^{-1/2}+\frac{3}{4}\ln (z)+h(z)$ with $h(z)\ll 1(z\to 0)$ and this means that a power series expansion can represent the remainder of the solution. The dominant balance method generates the leading term to this asymptotic expansion with constant $A$ and expansion coefficients determined by substitution into the full differential equation:

$w(z)=A{z}^{3/4}{e}^{\pm 2z^{-1/2}}\left(\sum _{n=0}^m{a}_n{z}^{n/2}\right)$

$a_{n+1}=\pm \frac{(n-1/2)(n+3/2)a_n}{4(n+1)}.$

A partial sum of this non-convergent series generates an approximate solution. The leading term corresponds to the Liouville-Green (LG) or Wentzel–Kramers–Brillouin (WKB) approximation.

References

• Bender, C.M.; Orszag, S.A. (1999). Advanced Mathematical Methods for Scientists and Engineers. Springer. ISBN 0-387-98931-5. 
• Callaham, Jared L.; Koch, James V.; Brunton, Bingni W.; Kutz, J. Nathan; Brunton, Steven L. (2021). "Learning dominant physical processes with data-driven balance models" (in en). Nature Communications 12 (1): 1016. doi:10.1038/s41467-021-21331-z. ISSN 2041-1723. PMID 33589607. Bibcode: 2021NatCo..12.1016C. 
• Christensen, Chris (1996). "Newton's Method for Resolving Affected Equations" (in en). The College Mathematics Journal 27 (5): 330–340. doi:10.1080/07468342.1996.11973804. ISSN 0746-8342. https://www.tandfonline.com/doi/abs/10.1080/07468342.1996.11973804. 
• de Bruijn, N. G. (1981), Asymptotic Methods in Analysis, Dover Publications, ISBN 9780486642215, https://books.google.com/books?id=Oqj9AgAAQBAJ 
• Fishaleck, T.; White, R.B. (2008). "Technical Report: The Use of Kruskal-Newton Diagrams for Differential Equations". Princeton Plasma Physics Laboratory (Princeton, NJ: U.S. Department of Energy Office of Scientific and Technical Information) (PPL-4289): 1–29. doi:10.2172/960287. https://www.osti.gov/biblio/960287. 
• Hinch, E. J. (1991) (in en). Perturbation Methods. Cambridge University Press. ISBN 978-0-521-37897-0. https://books.google.com/books?id=mqY4ZM0BWwIC. 
• Kruskal, M.D. (1962). "Technical Report: Asymptotology, Report MATT 160". Princeton Plasma Physics Laboratory (Princeton, NJ: Princeton University): 1–32. https://w3.pppl.gov/~hammett/work/2009/Kruskal62_Asymptotology_MATT160_ocr.pdf. 
• Neu, John C. (2015) (in en). Singular Perturbation in the Physical Sciences. American Mathematical Soc.. ISBN 978-1-4704-2555-5. https://books.google.com/books?id=rv8dCwAAQBAJ. 
• Olver, Frank William John. (1974) (in en). Introduction to Asymptotics and Special Functions. New York: Academic Press. ISBN 0-12-525856-9. https://books.google.com/books?id=rLjiBQAAQBAJ. 
• Paulsen, William (2013) (in en). Asymptotic Analysis and Perturbation Theory. CRC Press. ISBN 978-1-4665-1512-3. https://books.google.com/books?id=Z0DSBQAAQBAJ. 
• Rozman, Michael (2020). "Perturbation methods". University of Connecticut. https://www.phys.uconn.edu/~rozman/Courses/P2400_20F/downloads/syllabus.pdf. 
• White, R. B. (2010) (in en). Asymptotic Analysis of Differential Equations. World Scientific. ISBN 978-1-84816-607-3. https://books.google.com/books?id=F1OwONpC-N8C. 

See also

• Asymptotic analysis