First-order quasilinear PDE system analysis
$begingroup$
For $(t,r)in[0,infty)times[0,1]$ let be the following PDE's system
$$dot{vartheta}= w' +wvartheta' $$
$$dot{w}= vartheta'+ww' $$
along with initial conditions (i.c.)
$$w(0,r)=w_i(r)qquad vartheta(0,r)=vartheta_i(r) $$
and the boundary condition (b.c)
$$boxed{w(t,1)=0};Rightarrow;w(t,1)=dot{w}(t,1)=0;Rightarrow;boxed{vartheta'(t,1)=0} $$
My question is:
- What is the minimum extra b.c. needed for the problem to be well posed?
- Is there any specific analytical solution, if not a generic one?
- How can one modify the equations so that no shock waves occur? (I mean some modification that does not change the type and number of b.c. needed.)
Remark: If range $[0,1]$ specified for $r$ is interpreted as a sphere of radious $R=1$ then $r=0$ is just an interior point. So the fields $w$ and $vartheta$ at this point should not be restricted by any b.c. just as it happens for any other interior point.
Above all: What would be a good source (lecture notes, book, webpage, etc.) for studying this kind of system without having to revise all of PDE's theory?
I would appreciate any help.
PS: I use "dot" for $partial_t$ and "prime" for $partial_r$.
pde systems-of-equations boundary-value-problem hyperbolic-equations
$endgroup$
add a comment |
$begingroup$
For $(t,r)in[0,infty)times[0,1]$ let be the following PDE's system
$$dot{vartheta}= w' +wvartheta' $$
$$dot{w}= vartheta'+ww' $$
along with initial conditions (i.c.)
$$w(0,r)=w_i(r)qquad vartheta(0,r)=vartheta_i(r) $$
and the boundary condition (b.c)
$$boxed{w(t,1)=0};Rightarrow;w(t,1)=dot{w}(t,1)=0;Rightarrow;boxed{vartheta'(t,1)=0} $$
My question is:
- What is the minimum extra b.c. needed for the problem to be well posed?
- Is there any specific analytical solution, if not a generic one?
- How can one modify the equations so that no shock waves occur? (I mean some modification that does not change the type and number of b.c. needed.)
Remark: If range $[0,1]$ specified for $r$ is interpreted as a sphere of radious $R=1$ then $r=0$ is just an interior point. So the fields $w$ and $vartheta$ at this point should not be restricted by any b.c. just as it happens for any other interior point.
Above all: What would be a good source (lecture notes, book, webpage, etc.) for studying this kind of system without having to revise all of PDE's theory?
I would appreciate any help.
PS: I use "dot" for $partial_t$ and "prime" for $partial_r$.
pde systems-of-equations boundary-value-problem hyperbolic-equations
$endgroup$
$begingroup$
What do you mean by "dot" and "prime"? I imagine it is $tfrac{d}{dt}$ and $tfrac{d}{dr}$ but you really should write it in your question.
$endgroup$
– Giuseppe Negro
Dec 19 '18 at 14:13
add a comment |
$begingroup$
For $(t,r)in[0,infty)times[0,1]$ let be the following PDE's system
$$dot{vartheta}= w' +wvartheta' $$
$$dot{w}= vartheta'+ww' $$
along with initial conditions (i.c.)
$$w(0,r)=w_i(r)qquad vartheta(0,r)=vartheta_i(r) $$
and the boundary condition (b.c)
$$boxed{w(t,1)=0};Rightarrow;w(t,1)=dot{w}(t,1)=0;Rightarrow;boxed{vartheta'(t,1)=0} $$
My question is:
- What is the minimum extra b.c. needed for the problem to be well posed?
- Is there any specific analytical solution, if not a generic one?
- How can one modify the equations so that no shock waves occur? (I mean some modification that does not change the type and number of b.c. needed.)
Remark: If range $[0,1]$ specified for $r$ is interpreted as a sphere of radious $R=1$ then $r=0$ is just an interior point. So the fields $w$ and $vartheta$ at this point should not be restricted by any b.c. just as it happens for any other interior point.
Above all: What would be a good source (lecture notes, book, webpage, etc.) for studying this kind of system without having to revise all of PDE's theory?
I would appreciate any help.
PS: I use "dot" for $partial_t$ and "prime" for $partial_r$.
pde systems-of-equations boundary-value-problem hyperbolic-equations
$endgroup$
For $(t,r)in[0,infty)times[0,1]$ let be the following PDE's system
$$dot{vartheta}= w' +wvartheta' $$
$$dot{w}= vartheta'+ww' $$
along with initial conditions (i.c.)
$$w(0,r)=w_i(r)qquad vartheta(0,r)=vartheta_i(r) $$
and the boundary condition (b.c)
$$boxed{w(t,1)=0};Rightarrow;w(t,1)=dot{w}(t,1)=0;Rightarrow;boxed{vartheta'(t,1)=0} $$
My question is:
- What is the minimum extra b.c. needed for the problem to be well posed?
- Is there any specific analytical solution, if not a generic one?
- How can one modify the equations so that no shock waves occur? (I mean some modification that does not change the type and number of b.c. needed.)
Remark: If range $[0,1]$ specified for $r$ is interpreted as a sphere of radious $R=1$ then $r=0$ is just an interior point. So the fields $w$ and $vartheta$ at this point should not be restricted by any b.c. just as it happens for any other interior point.
Above all: What would be a good source (lecture notes, book, webpage, etc.) for studying this kind of system without having to revise all of PDE's theory?
I would appreciate any help.
PS: I use "dot" for $partial_t$ and "prime" for $partial_r$.
pde systems-of-equations boundary-value-problem hyperbolic-equations
pde systems-of-equations boundary-value-problem hyperbolic-equations
edited Dec 20 '18 at 15:25
Harry49
7,55431342
7,55431342
asked Dec 19 '18 at 14:11
dkstackdkstack
204
204
$begingroup$
What do you mean by "dot" and "prime"? I imagine it is $tfrac{d}{dt}$ and $tfrac{d}{dr}$ but you really should write it in your question.
$endgroup$
– Giuseppe Negro
Dec 19 '18 at 14:13
add a comment |
$begingroup$
What do you mean by "dot" and "prime"? I imagine it is $tfrac{d}{dt}$ and $tfrac{d}{dr}$ but you really should write it in your question.
$endgroup$
– Giuseppe Negro
Dec 19 '18 at 14:13
$begingroup$
What do you mean by "dot" and "prime"? I imagine it is $tfrac{d}{dt}$ and $tfrac{d}{dr}$ but you really should write it in your question.
$endgroup$
– Giuseppe Negro
Dec 19 '18 at 14:13
$begingroup$
What do you mean by "dot" and "prime"? I imagine it is $tfrac{d}{dt}$ and $tfrac{d}{dr}$ but you really should write it in your question.
$endgroup$
– Giuseppe Negro
Dec 19 '18 at 14:13
add a comment |
1 Answer
1
active
oldest
votes
$begingroup$
This system can be written in the quasi-linear form $u_t + A(u), u_r = 0$ where $u = (vartheta,w)^top$ and
$$
A(u) = -begin{pmatrix}w & 1 \ 1 & wend{pmatrix}
$$
is a symmetric matrix. The above PDE system is a first-order symmetric hyperbolic quasi-linear system. Thus, the matrix $A(u)$ can be diagonalized with real eigenvalues $lambda in lbrace-1-w, 1-wrbrace$. Here, the two characteristic fields associated with the eigenvalues are both genuinely nonlinear, i.e. the gradient $lambda_u$ of the eigenvalue $lambda$ is not orthogonal to the corresponding eigenvector (1). Since this system is of first order, the boundary value problem is locally well-posed for small data if we know the function $vartheta_1$ such that $vartheta(t,1) = vartheta_1(t)$ along the boundary $r=1$. Indeed, we can verify that there is no zero eigenvalue on the boundary $r=1$ (which is the case here), so that information will propagate from the boundary.
Concerning further results such as analytical solutions, it would be advantageous to write this system in conservation form $u_t + f(u)_r = 0$, but I doubt that this is possible. The present type of system produces discontinuous solutions in finite time even if the boundary data is smooth. To prevent this, a dissipative right-hand side of the form $epsilon u_{rr}$ or $-epsilon u$ can be added.
(1) P.D. Lax (1973): Hyperbolic Systems of Conservation Laws and the Mathematical Theory of Shock Waves, SIAM. doi:10.1137/1.9781611970562
$endgroup$
add a comment |
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%2f3046437%2ffirst-order-quasilinear-pde-system-analysis%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
1 Answer
1
active
oldest
votes
1 Answer
1
active
oldest
votes
active
oldest
votes
active
oldest
votes
$begingroup$
This system can be written in the quasi-linear form $u_t + A(u), u_r = 0$ where $u = (vartheta,w)^top$ and
$$
A(u) = -begin{pmatrix}w & 1 \ 1 & wend{pmatrix}
$$
is a symmetric matrix. The above PDE system is a first-order symmetric hyperbolic quasi-linear system. Thus, the matrix $A(u)$ can be diagonalized with real eigenvalues $lambda in lbrace-1-w, 1-wrbrace$. Here, the two characteristic fields associated with the eigenvalues are both genuinely nonlinear, i.e. the gradient $lambda_u$ of the eigenvalue $lambda$ is not orthogonal to the corresponding eigenvector (1). Since this system is of first order, the boundary value problem is locally well-posed for small data if we know the function $vartheta_1$ such that $vartheta(t,1) = vartheta_1(t)$ along the boundary $r=1$. Indeed, we can verify that there is no zero eigenvalue on the boundary $r=1$ (which is the case here), so that information will propagate from the boundary.
Concerning further results such as analytical solutions, it would be advantageous to write this system in conservation form $u_t + f(u)_r = 0$, but I doubt that this is possible. The present type of system produces discontinuous solutions in finite time even if the boundary data is smooth. To prevent this, a dissipative right-hand side of the form $epsilon u_{rr}$ or $-epsilon u$ can be added.
(1) P.D. Lax (1973): Hyperbolic Systems of Conservation Laws and the Mathematical Theory of Shock Waves, SIAM. doi:10.1137/1.9781611970562
$endgroup$
add a comment |
$begingroup$
This system can be written in the quasi-linear form $u_t + A(u), u_r = 0$ where $u = (vartheta,w)^top$ and
$$
A(u) = -begin{pmatrix}w & 1 \ 1 & wend{pmatrix}
$$
is a symmetric matrix. The above PDE system is a first-order symmetric hyperbolic quasi-linear system. Thus, the matrix $A(u)$ can be diagonalized with real eigenvalues $lambda in lbrace-1-w, 1-wrbrace$. Here, the two characteristic fields associated with the eigenvalues are both genuinely nonlinear, i.e. the gradient $lambda_u$ of the eigenvalue $lambda$ is not orthogonal to the corresponding eigenvector (1). Since this system is of first order, the boundary value problem is locally well-posed for small data if we know the function $vartheta_1$ such that $vartheta(t,1) = vartheta_1(t)$ along the boundary $r=1$. Indeed, we can verify that there is no zero eigenvalue on the boundary $r=1$ (which is the case here), so that information will propagate from the boundary.
Concerning further results such as analytical solutions, it would be advantageous to write this system in conservation form $u_t + f(u)_r = 0$, but I doubt that this is possible. The present type of system produces discontinuous solutions in finite time even if the boundary data is smooth. To prevent this, a dissipative right-hand side of the form $epsilon u_{rr}$ or $-epsilon u$ can be added.
(1) P.D. Lax (1973): Hyperbolic Systems of Conservation Laws and the Mathematical Theory of Shock Waves, SIAM. doi:10.1137/1.9781611970562
$endgroup$
add a comment |
$begingroup$
This system can be written in the quasi-linear form $u_t + A(u), u_r = 0$ where $u = (vartheta,w)^top$ and
$$
A(u) = -begin{pmatrix}w & 1 \ 1 & wend{pmatrix}
$$
is a symmetric matrix. The above PDE system is a first-order symmetric hyperbolic quasi-linear system. Thus, the matrix $A(u)$ can be diagonalized with real eigenvalues $lambda in lbrace-1-w, 1-wrbrace$. Here, the two characteristic fields associated with the eigenvalues are both genuinely nonlinear, i.e. the gradient $lambda_u$ of the eigenvalue $lambda$ is not orthogonal to the corresponding eigenvector (1). Since this system is of first order, the boundary value problem is locally well-posed for small data if we know the function $vartheta_1$ such that $vartheta(t,1) = vartheta_1(t)$ along the boundary $r=1$. Indeed, we can verify that there is no zero eigenvalue on the boundary $r=1$ (which is the case here), so that information will propagate from the boundary.
Concerning further results such as analytical solutions, it would be advantageous to write this system in conservation form $u_t + f(u)_r = 0$, but I doubt that this is possible. The present type of system produces discontinuous solutions in finite time even if the boundary data is smooth. To prevent this, a dissipative right-hand side of the form $epsilon u_{rr}$ or $-epsilon u$ can be added.
(1) P.D. Lax (1973): Hyperbolic Systems of Conservation Laws and the Mathematical Theory of Shock Waves, SIAM. doi:10.1137/1.9781611970562
$endgroup$
This system can be written in the quasi-linear form $u_t + A(u), u_r = 0$ where $u = (vartheta,w)^top$ and
$$
A(u) = -begin{pmatrix}w & 1 \ 1 & wend{pmatrix}
$$
is a symmetric matrix. The above PDE system is a first-order symmetric hyperbolic quasi-linear system. Thus, the matrix $A(u)$ can be diagonalized with real eigenvalues $lambda in lbrace-1-w, 1-wrbrace$. Here, the two characteristic fields associated with the eigenvalues are both genuinely nonlinear, i.e. the gradient $lambda_u$ of the eigenvalue $lambda$ is not orthogonal to the corresponding eigenvector (1). Since this system is of first order, the boundary value problem is locally well-posed for small data if we know the function $vartheta_1$ such that $vartheta(t,1) = vartheta_1(t)$ along the boundary $r=1$. Indeed, we can verify that there is no zero eigenvalue on the boundary $r=1$ (which is the case here), so that information will propagate from the boundary.
Concerning further results such as analytical solutions, it would be advantageous to write this system in conservation form $u_t + f(u)_r = 0$, but I doubt that this is possible. The present type of system produces discontinuous solutions in finite time even if the boundary data is smooth. To prevent this, a dissipative right-hand side of the form $epsilon u_{rr}$ or $-epsilon u$ can be added.
(1) P.D. Lax (1973): Hyperbolic Systems of Conservation Laws and the Mathematical Theory of Shock Waves, SIAM. doi:10.1137/1.9781611970562
edited Dec 20 '18 at 17:26
answered Dec 20 '18 at 17:10
Harry49Harry49
7,55431342
7,55431342
add a comment |
add a comment |
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.
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%2f3046437%2ffirst-order-quasilinear-pde-system-analysis%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
$begingroup$
What do you mean by "dot" and "prime"? I imagine it is $tfrac{d}{dt}$ and $tfrac{d}{dr}$ but you really should write it in your question.
$endgroup$
– Giuseppe Negro
Dec 19 '18 at 14:13