What's the name of this “shooting”-like method for numerical solving of PDE eigenproblems?
Consider a PDE of the following form:
$$left(partial_r^2+frac5rpartial_r+frac4{r^2} hat Lright)Psi(r,p)+(E- V(p)U(r))Psi(r,p)=0,tag1$$
where $hat L$ is a differential operator differentiating only with respect to $p$ (i.e. ignoring $r$) and $V$, $U$ are some functions.
Suppose that boundary conditions are given such that $|Psi(0,p)|<infty$ and $|Psi(infty,p)|<infty$, and some appropriate homogeneous conditions along the $p$ variable, independent from $r$.
I've read about a method for solving it numerically, in particular finding the eigenvalues, which I'll outline here.
First, a the solution $Psi(r,p)$ is expanded into series over a set of orthogonal functions $b_k(y)$:
$$Psi(r,p)=sumlimits_{k=1}^infty R_k(r) b_k(p),tag2$$
where $R_k$ are $r$-dependent coefficients of expansion. These coefficients are collected into a column vector $R(r)$.
Expanding the $p$-dependent operators $V(p)$ and $hat L$ in the same basis functions $b_k$, thus obtaining respectively the matrices $tilde V$ and $tilde L$, we get the following system of ODEs:
$$left(partial_r^2+frac5rpartial_r+frac4{r^2} tilde Lright)R(r)+(I E-tilde V U(r))R(r)=0,tag3$$
where $I$ is the unit matrix
Then we introduce a matrix $F(r)$, which is defined as follows:
$$r R'(r)=F(r)R(r).tag4$$
After inserting it into $(3)$, we get the following nonlinear equation for $F(r)$:
$$rF'(r)+F(r)^2+4F(r)+4tilde L+(IE-tilde V)r^2=0.tag5$$
Now the equation $(5)$ automatically gives us an initial condition:
$$F(0)^2+4F(0)=-4tilde L,tag6$$
from which we can take (with some good definition of square root of a matrix):
$$F(0)=-2I+2sqrt{I-L}.tag7$$
Now (after truncating the series so that we work with finite-dimensional matrices) we can choose some value of $E$, and see if the solution of $(5)$ reaches a pole before we reach the second boundary. If it does, then we've found a possible zero of $R(r)$, and we have to lower $E$ to move the zero farther to the second boundary. Otherwise, $E$ is too small, and we should increase it to move the zero of $R(r)$ closer to the second boundary. Bisecting, we can find the eigenvalue as we do in the shooting method for ODE.
Other eigenvalues could be found by switching to an equation for $G(r)=F(r)^{-1}$ near the pole, and then resuming propagation for $F(r)$ itself.
My question: what is the name of this method? Where can I read more about it, its properties, pitfalls etc.? The source where I first encountered it simply described it similarly to my above description, without any justification of e.g. the choice of initial condition $(7)$ among all possible solutions of $(6)$, or any other explanations.
pde numerical-methods eigenfunctions
add a comment |
Consider a PDE of the following form:
$$left(partial_r^2+frac5rpartial_r+frac4{r^2} hat Lright)Psi(r,p)+(E- V(p)U(r))Psi(r,p)=0,tag1$$
where $hat L$ is a differential operator differentiating only with respect to $p$ (i.e. ignoring $r$) and $V$, $U$ are some functions.
Suppose that boundary conditions are given such that $|Psi(0,p)|<infty$ and $|Psi(infty,p)|<infty$, and some appropriate homogeneous conditions along the $p$ variable, independent from $r$.
I've read about a method for solving it numerically, in particular finding the eigenvalues, which I'll outline here.
First, a the solution $Psi(r,p)$ is expanded into series over a set of orthogonal functions $b_k(y)$:
$$Psi(r,p)=sumlimits_{k=1}^infty R_k(r) b_k(p),tag2$$
where $R_k$ are $r$-dependent coefficients of expansion. These coefficients are collected into a column vector $R(r)$.
Expanding the $p$-dependent operators $V(p)$ and $hat L$ in the same basis functions $b_k$, thus obtaining respectively the matrices $tilde V$ and $tilde L$, we get the following system of ODEs:
$$left(partial_r^2+frac5rpartial_r+frac4{r^2} tilde Lright)R(r)+(I E-tilde V U(r))R(r)=0,tag3$$
where $I$ is the unit matrix
Then we introduce a matrix $F(r)$, which is defined as follows:
$$r R'(r)=F(r)R(r).tag4$$
After inserting it into $(3)$, we get the following nonlinear equation for $F(r)$:
$$rF'(r)+F(r)^2+4F(r)+4tilde L+(IE-tilde V)r^2=0.tag5$$
Now the equation $(5)$ automatically gives us an initial condition:
$$F(0)^2+4F(0)=-4tilde L,tag6$$
from which we can take (with some good definition of square root of a matrix):
$$F(0)=-2I+2sqrt{I-L}.tag7$$
Now (after truncating the series so that we work with finite-dimensional matrices) we can choose some value of $E$, and see if the solution of $(5)$ reaches a pole before we reach the second boundary. If it does, then we've found a possible zero of $R(r)$, and we have to lower $E$ to move the zero farther to the second boundary. Otherwise, $E$ is too small, and we should increase it to move the zero of $R(r)$ closer to the second boundary. Bisecting, we can find the eigenvalue as we do in the shooting method for ODE.
Other eigenvalues could be found by switching to an equation for $G(r)=F(r)^{-1}$ near the pole, and then resuming propagation for $F(r)$ itself.
My question: what is the name of this method? Where can I read more about it, its properties, pitfalls etc.? The source where I first encountered it simply described it similarly to my above description, without any justification of e.g. the choice of initial condition $(7)$ among all possible solutions of $(6)$, or any other explanations.
pde numerical-methods eigenfunctions
add a comment |
Consider a PDE of the following form:
$$left(partial_r^2+frac5rpartial_r+frac4{r^2} hat Lright)Psi(r,p)+(E- V(p)U(r))Psi(r,p)=0,tag1$$
where $hat L$ is a differential operator differentiating only with respect to $p$ (i.e. ignoring $r$) and $V$, $U$ are some functions.
Suppose that boundary conditions are given such that $|Psi(0,p)|<infty$ and $|Psi(infty,p)|<infty$, and some appropriate homogeneous conditions along the $p$ variable, independent from $r$.
I've read about a method for solving it numerically, in particular finding the eigenvalues, which I'll outline here.
First, a the solution $Psi(r,p)$ is expanded into series over a set of orthogonal functions $b_k(y)$:
$$Psi(r,p)=sumlimits_{k=1}^infty R_k(r) b_k(p),tag2$$
where $R_k$ are $r$-dependent coefficients of expansion. These coefficients are collected into a column vector $R(r)$.
Expanding the $p$-dependent operators $V(p)$ and $hat L$ in the same basis functions $b_k$, thus obtaining respectively the matrices $tilde V$ and $tilde L$, we get the following system of ODEs:
$$left(partial_r^2+frac5rpartial_r+frac4{r^2} tilde Lright)R(r)+(I E-tilde V U(r))R(r)=0,tag3$$
where $I$ is the unit matrix
Then we introduce a matrix $F(r)$, which is defined as follows:
$$r R'(r)=F(r)R(r).tag4$$
After inserting it into $(3)$, we get the following nonlinear equation for $F(r)$:
$$rF'(r)+F(r)^2+4F(r)+4tilde L+(IE-tilde V)r^2=0.tag5$$
Now the equation $(5)$ automatically gives us an initial condition:
$$F(0)^2+4F(0)=-4tilde L,tag6$$
from which we can take (with some good definition of square root of a matrix):
$$F(0)=-2I+2sqrt{I-L}.tag7$$
Now (after truncating the series so that we work with finite-dimensional matrices) we can choose some value of $E$, and see if the solution of $(5)$ reaches a pole before we reach the second boundary. If it does, then we've found a possible zero of $R(r)$, and we have to lower $E$ to move the zero farther to the second boundary. Otherwise, $E$ is too small, and we should increase it to move the zero of $R(r)$ closer to the second boundary. Bisecting, we can find the eigenvalue as we do in the shooting method for ODE.
Other eigenvalues could be found by switching to an equation for $G(r)=F(r)^{-1}$ near the pole, and then resuming propagation for $F(r)$ itself.
My question: what is the name of this method? Where can I read more about it, its properties, pitfalls etc.? The source where I first encountered it simply described it similarly to my above description, without any justification of e.g. the choice of initial condition $(7)$ among all possible solutions of $(6)$, or any other explanations.
pde numerical-methods eigenfunctions
Consider a PDE of the following form:
$$left(partial_r^2+frac5rpartial_r+frac4{r^2} hat Lright)Psi(r,p)+(E- V(p)U(r))Psi(r,p)=0,tag1$$
where $hat L$ is a differential operator differentiating only with respect to $p$ (i.e. ignoring $r$) and $V$, $U$ are some functions.
Suppose that boundary conditions are given such that $|Psi(0,p)|<infty$ and $|Psi(infty,p)|<infty$, and some appropriate homogeneous conditions along the $p$ variable, independent from $r$.
I've read about a method for solving it numerically, in particular finding the eigenvalues, which I'll outline here.
First, a the solution $Psi(r,p)$ is expanded into series over a set of orthogonal functions $b_k(y)$:
$$Psi(r,p)=sumlimits_{k=1}^infty R_k(r) b_k(p),tag2$$
where $R_k$ are $r$-dependent coefficients of expansion. These coefficients are collected into a column vector $R(r)$.
Expanding the $p$-dependent operators $V(p)$ and $hat L$ in the same basis functions $b_k$, thus obtaining respectively the matrices $tilde V$ and $tilde L$, we get the following system of ODEs:
$$left(partial_r^2+frac5rpartial_r+frac4{r^2} tilde Lright)R(r)+(I E-tilde V U(r))R(r)=0,tag3$$
where $I$ is the unit matrix
Then we introduce a matrix $F(r)$, which is defined as follows:
$$r R'(r)=F(r)R(r).tag4$$
After inserting it into $(3)$, we get the following nonlinear equation for $F(r)$:
$$rF'(r)+F(r)^2+4F(r)+4tilde L+(IE-tilde V)r^2=0.tag5$$
Now the equation $(5)$ automatically gives us an initial condition:
$$F(0)^2+4F(0)=-4tilde L,tag6$$
from which we can take (with some good definition of square root of a matrix):
$$F(0)=-2I+2sqrt{I-L}.tag7$$
Now (after truncating the series so that we work with finite-dimensional matrices) we can choose some value of $E$, and see if the solution of $(5)$ reaches a pole before we reach the second boundary. If it does, then we've found a possible zero of $R(r)$, and we have to lower $E$ to move the zero farther to the second boundary. Otherwise, $E$ is too small, and we should increase it to move the zero of $R(r)$ closer to the second boundary. Bisecting, we can find the eigenvalue as we do in the shooting method for ODE.
Other eigenvalues could be found by switching to an equation for $G(r)=F(r)^{-1}$ near the pole, and then resuming propagation for $F(r)$ itself.
My question: what is the name of this method? Where can I read more about it, its properties, pitfalls etc.? The source where I first encountered it simply described it similarly to my above description, without any justification of e.g. the choice of initial condition $(7)$ among all possible solutions of $(6)$, or any other explanations.
pde numerical-methods eigenfunctions
pde numerical-methods eigenfunctions
edited 2 days ago
asked Jan 3 at 23:02
Ruslan
3,70421533
3,70421533
add a comment |
add a comment |
0
active
oldest
votes
Your Answer
StackExchange.ifUsing("editor", function () {
return StackExchange.using("mathjaxEditing", function () {
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
});
});
}, "mathjax-editing");
StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "69"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);
StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});
function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
noCode: true, onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});
}
});
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3061123%2fwhats-the-name-of-this-shooting-like-method-for-numerical-solving-of-pde-eige%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
0
active
oldest
votes
0
active
oldest
votes
active
oldest
votes
active
oldest
votes
Thanks for contributing an answer to Mathematics Stack Exchange!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
Use MathJax to format equations. MathJax reference.
To learn more, see our tips on writing great answers.
Some of your past answers have not been well-received, and you're in danger of being blocked from answering.
Please pay close attention to the following guidance:
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3061123%2fwhats-the-name-of-this-shooting-like-method-for-numerical-solving-of-pde-eige%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown