|
|
| Line 1: |
Line 1: |
| In [[mathematics]], the '''method of matched asymptotic expansions''' is a common approach to finding an accurate approximation to the solution to an [[equation]], or [[simultaneous equations|system of equations]]. It is particularly used when solving [[singular perturbation|singularly perturbed]] [[differential equation]]s. It involves finding several different approximate solutions, each of which is valid (i.e. accurate) for part of the range of the independent variable, and then combining these different solutions together to give a single approximate solution that is valid for the whole range of values of the independent variable.
| | Golda is what's created on my beginning certification although it is not the name on my beginning certificate. free psychic readings - [http://www.indosfriends.com/profile-253/info/ indosfriends.com] - My wife and I reside in Mississippi and I [http://165.132.39.93/xe/visitors/372912 best psychics] love each day residing right here. He is an order clerk and it's something he truly enjoy. What me and my family members love is bungee leaping but I've been using on new issues recently.<br><br>Look into my page ... email psychic readings ([http://Checkmates.Co.za/index.php?do=/profile-56347/info/ Suggested Web page]) |
| | |
| ==Method overview==
| |
| In a large class of singularly perturbed problems, the [[Domain of a function|domain]] may be divided into two or more subdomains. In one of these, often the largest, the solution is accurately approximated by an [[asymptotic series]] found by treating the problem as a regular [[perturbation theory|perturbation]] (i.e. by setting a relatively small parameter to zero). The other subdomains consist of one or more small areas in which that approximation is inaccurate, generally because the perturbation terms in the problem are not [[negligible]] there. These areas are referred to as transition layers, and as boundary or interior layers depending on whether they occur at the domain boundary (as is the usual case in applications) or inside the domain.
| |
| | |
| An approximation in the form of an asymptotic series is obtained in the transition layer(s) by treating that part of the domain as a separate perturbation problem. This approximation is called the "inner solution," and the other is the "outer solution," named for their relationship to the transition layer(s). The outer and inner solutions are then combined through a process called "matching" in such a way that an approximate solution for the whole domain is obtained.<ref name=verhulst>{{cite book | author=Verhulst, F. | title=Methods and Applications of Singular Perturbations: Boundary Layers and Multiple Timescale Dynamics | publisher=Springer | year=2005 | isbn=0-387-22966-3 }}</ref><ref>{{Cite book | author=Nayfeh, A. H. | title=Perturbation Methods | series=Wiley Classics Library | year=2000 | publisher=Wiley-Interscience | isbn=978-0-471-39917-9 }}</ref><ref>{{Cite book | author=Kevorkian, J.; Cole, J. D. | title=Multiple scale and singular perturbation methods | year=1996 | publisher=Springer | isbn=0-387-94202-5 }}</ref>
| |
| | |
| ==A simple example==
| |
| Consider the [[boundary value problem]]
| |
| | |
| :<math> \epsilon y'' + (1+\epsilon) y' + y = 0,</math> | |
| | |
| where <math>y</math> is a function of independent time variable <math>t</math>, which ranges from 0 to 1, the boundary conditions are <math>y(0)=0</math> and <math>y(1)=1</math>, and <math>\epsilon</math> is a small parameter, such that <math>0<\epsilon\ll 1</math>.
| |
| | |
| ===Outer solution, valid for ''t'' = ''O''(1)===
| |
| | |
| Since <math>\epsilon</math> is very small, our first approach is to treat the equation as a regular perturbation problem, i.e. make the approximation <math>\epsilon=0</math>, and hence find the solution to the problem
| |
| | |
| :<math>y'+y=0.\,</math>
| |
| | |
| Alternatively, consider that when <math>y</math> and <math>t</math> are both of size ''O''(1), the four [[Term (mathematics)|terms]] on the left hand side of the original equation are respectively of sizes ''O''(<math>\epsilon</math>), ''O''(1), ''O''(<math>\epsilon</math>) and ''O''(1). The [[leading-order]] balance on this timescale, valid in the distinguished limit <math>\epsilon \to 0</math>, is therefore given by the second and fourth terms, i.e. <math>y'+y=0.\,</math>
| |
| | |
| This has solution
| |
| | |
| :<math>y=Ae^{-t}\,</math> | |
| | |
| for some constant <math>A</math>. Applying the boundary condition <math>y(0) = 0</math>, we would have <math>A=0</math>; applying the boundary condition <math>y(1) = 1</math>, we would have <math>A=e</math>. It is therefore impossible to satisfy both boundary conditions, so <math>\epsilon=0</math> is not a valid approximation to make across the whole of the domain (i.e. this is a [[singular perturbation]] problem). From this we infer that there must be a boundary layer at one of the endpoints of the domain where <math>\epsilon</math> needs to be included. This region will be where <math>\epsilon</math> is no longer [[negligible]] compared to the independent variable <math>t</math>, i.e. <math>t</math> and <math>\epsilon</math> are of comparable size, i.e. the boundary layer is adjacent to <math>t=0</math>. Therefore the other boundary condition <math>y(1) = 1</math> applies in this outer region, so <math>A=e</math>, i.e. <math>y_O=e^{1-t}\,</math> is an accurate approximate solution to the original boundary value problem in this outer region. It is the leading-order solution.
| |
| | |
| ===Inner solution, valid for ''t'' = ''O''(''ε'')===
| |
| | |
| In the inner region, <math>t</math> and <math>\epsilon</math> are both tiny, but of comparable size, so define the new ''O''(1) time variable <math>\tau = t/\epsilon</math>. Rescale the original boundary value problem by replacing <math>t</math> with <math>\tau\epsilon</math>, and the problem becomes
| |
| | |
| :<math> \frac{1}{\epsilon} y''(\tau ) + \left( {1 + \epsilon } \right)\frac{1}{\epsilon }y'(\tau ) + y(\tau ) = 0,\,</math>
| |
| | |
| which, after multiplying by <math>\epsilon</math> and taking <math>\epsilon = 0</math>, is
| |
| | |
| :<math>y'' + y' = 0. \, </math>
| |
| | |
| Alternatively, consider that when <math>t</math> has reduced to size ''O''(<math>\epsilon</math>), then <math>y</math> is still of size ''O''(1) (using the expression for <math>y_O</math>), and so the four terms on the left hand side of the original equation are respectively of sizes ''O''(<math>\epsilon</math><sup>−1</sup>), ''O''(<math>\epsilon</math><sup>−1</sup>), ''O''(1) and ''O''(1). The [[leading-order]] balance on this timescale, valid in the distinguished limit <math>\epsilon \to 0</math>, is therefore given by the first and second terms, i.e. <math>y'' + y'=0.\,</math>
| |
| | |
| This has solution
| |
| | |
| :<math>y=B-Ce^{-\tau}\,</math>
| |
| | |
| for some constants <math>B</math> and <math>C</math>. Since <math>y(0)=0</math> applies in this inner region, this gives <math>B=C</math>, so an accurate approximate solution to the original boundary value problem in this inner region (it is the leading-order solution) is
| |
| | |
| :<math>y_I = B\left( {1 - e^{ - \tau } } \right)= B\left( {1 - e^{ - t/\epsilon } } \right).\,</math> | |
| | |
| ===Matching===
| |
| | |
| We use matching to find the value of the constant <math>B</math>. The idea of matching is that the inner and outer solutions should agree for values of <math>t</math> in an intermediate (or overlap) region, i.e. where <math>\epsilon \ll t \ll 1</math>. We need the outer limit of the inner solution to match the inner limit of the outer solution, i.e.
| |
| <math>\lim_{\tau \rightarrow \infty} y_I = \lim_{t \to 0} y_O ,\,</math>
| |
| which gives <math>B=e</math>.
| |
| | |
| ===Composite solution===
| |
| | |
| To obtain our final, matched, composite solution, valid on the whole domain, one popular method is the uniform method. In this method, we add the inner and outer approximations and subtract their overlapping value, <math>\,y_\mathrm{overlap}</math>, which would otherwise be counted twice. The overlapping value is the outer limit of the inner boundary layer solution, and the inner limit of the outer solution; these limits were above found to equal <math>e</math>. Therefore, the final approximate solution to this boundary value problem is,
| |
| | |
| :<math>y(t) = y_I + y_O - y_\mathrm{overlap} = e\left( {1 - e^{ - t/\epsilon } } \right) + e^{1 - t} - e = e\left( {e^{ - t} - e^{ - t/\epsilon } } \right).\,</math>
| |
| | |
| Note that this expression correctly reduces to the expressions for <math>y_I</math> and <math>y_O</math> when <math>t</math> is ''O''(<math>\epsilon</math>) and ''O''(1), respectively.
| |
| | |
| ===Accuracy===
| |
| [[Image:Singular_perturbation_convergence.jpg|thumb|right|400px|Convergence of approximations. Approximations and exact solutions, which are indistinguishable at this scale, are shown for various <math>\epsilon</math>. The outer solution is also shown. Note that since the boundary layer becomes narrower with decreasing <math>\epsilon</math>, the approximations [[Limit of a sequence|converge]] to the outer solution [[Pointwise convergence|pointwise]], but not [[Uniform convergence|uniformly]], almost everywhere.]]
| |
| | |
| This final solution satisfies the problem's original differential equation (shown by substituting it and its derivatives into the original equation). Also, the boundary conditions produced by this final solution match the values given in the problem. This implies, due to the uniqueness of the solution, that the matched asymptotic solution is identical to the exact solution up to a constant multiple. This is not necessarily always the case, any remaining terms should go to zero uniformly as <math> \epsilon \rightarrow 0 </math>.
| |
| | |
| Not only does our solution successfully approximately solve the problem at hand, it closely approximates the problem's exact solution. It happens that this particular problem is easily found to have exact solution
| |
| | |
| :<math>y(t) = \frac{{e^{ - t} - e^{ - t/\varepsilon } }}{{e^{ - 1} - e^{ - 1/\varepsilon } }},\,</math>
| |
| | |
| which has the same form as the approximate solution, bar the multiplying constant. Note also that the approximate solution is the first term in a binomial expansion of the exact solution in powers of <math>y(1) = e^{1 - 1/\epsilon }</math>.
| |
| | |
| ===Location of boundary layer===
| |
| Conveniently, we can see that the boundary layer, where <math>y'</math> and <math>y''</math> are large, is near <math>t=0</math>, as we supposed earlier. If we had supposed it to be at the other endpoint and proceeded by making the rescaling <math>\tau = (1 - t)/\epsilon</math>, we would have found it impossible to satisfy the resulting matching condition. For many problems, this kind of trial and error is the only way to determine the true location of the boundary layer.<ref name=verhulst/>
| |
| | |
| ==Harder problems==
| |
| | |
| The problem above is a simple example because it is a single equation with only one dependent variable, and there is one boundary layer in the solution. Harder problems may contain several co-dependent variables in a system of several equations, and/or with several boundary and/or interior layers in the solution.
| |
| | |
| ==See also==
| |
| * [[Multiple-scale analysis]]
| |
| | |
| ==References==
| |
| {{reflist}}
| |
| | |
| [[Category:Differential equations]]
| |
| [[Category:Asymptotic analysis]]
| |
Golda is what's created on my beginning certification although it is not the name on my beginning certificate. free psychic readings - indosfriends.com - My wife and I reside in Mississippi and I best psychics love each day residing right here. He is an order clerk and it's something he truly enjoy. What me and my family members love is bungee leaping but I've been using on new issues recently.
Look into my page ... email psychic readings (Suggested Web page)