Analytical solution for convection diffusion equation
$begingroup$
Convection-diffusion equation is:
$frac{partial u}{partial t} + frac{partial u}{partial x} = 0.01frac{partial^2 u}{partial x^2}$
Inital conditon is:
$u(x,0) = sin(x)$ over the domain 0 to $2pi$ with periodic boundary conditon that is $u(0,t) =u(2pi,t)$.
Second boundary condition is clamped at $x = 0 $ that is $frac{du}{dx} = 0$
I would like to know
1) Whether this problem has an exact solution? if so please prove the solution.
2) Can any symbolic computing software like Maple, Mathematica, Matlab can solve this problem analytically?
3) Please provide some good tutorial (external links) for finding the analytical solution of the advection-diffusion equation.
Please don't provide a numerical solution because this problem is a toy problem in numerical methods.
pde maple mathematica
asked Jan 3 at 6:22
Arun Govind Neelan AArun Govind Neelan A
17027
$endgroup$
add a comment |
$begingroup$
Convection-diffusion equation is:
$frac{partial u}{partial t} + frac{partial u}{partial x} = 0.01frac{partial^2 u}{partial x^2}$
Inital conditon is:
$u(x,0) = sin(x)$ over the domain 0 to $2pi$ with periodic boundary conditon that is $u(0,t) =u(2pi,t)$.
Second boundary condition is clamped at $x = 0 $ that is $frac{du}{dx} = 0$
I would like to know
1) Whether this problem has an exact solution? if so please prove the solution.
2) Can any symbolic computing software like Maple, Mathematica, Matlab can solve this problem analytically?
3) Please provide some good tutorial (external links) for finding the analytical solution of the advection-diffusion equation.
Please don't provide a numerical solution because this problem is a toy problem in numerical methods.
pde maple mathematica
asked Jan 3 at 6:22
Arun Govind Neelan AArun Govind Neelan A
17027
$endgroup$
$begingroup$
Note, if $u$ solves the above convection-diffusion then $h(x, t)=u(x+t, t)$ solves begin{align} partial_t h =partial_t u + partial_x u= .01 partial_{xx} u = .01partial_{xx}h. end{align} Hence $h$ solves the heat equation.
$endgroup$
– Jacky Chong
Jan 3 at 7:14
$begingroup$
$u(x,t) = exp(-nu t)sin(x-ct)$. Is this the solution to this equation? I found "some" match with the numerical solution but not satisfied. I think I made some mistake. Solution following the physics but I'm not sure whether it is correct or not
$endgroup$
– Arun Govind Neelan A
Jan 3 at 7:29
$begingroup$
In addition to that grid refinement study shows a good match with analytical and numerical result. Could you tell me any way to cross check it? The solution basis is good but I don't know the impact of the boundary condition on the solution and how to cross check it.
$endgroup$
– Arun Govind Neelan A
Jan 3 at 7:43
add a comment |
$begingroup$
Convection-diffusion equation is:
$frac{partial u}{partial t} + frac{partial u}{partial x} = 0.01frac{partial^2 u}{partial x^2}$
Inital conditon is:
$u(x,0) = sin(x)$ over the domain 0 to $2pi$ with periodic boundary conditon that is $u(0,t) =u(2pi,t)$.
Second boundary condition is clamped at $x = 0 $ that is $frac{du}{dx} = 0$
I would like to know
1) Whether this problem has an exact solution? if so please prove the solution.
2) Can any symbolic computing software like Maple, Mathematica, Matlab can solve this problem analytically?
3) Please provide some good tutorial (external links) for finding the analytical solution of the advection-diffusion equation.
Please don't provide a numerical solution because this problem is a toy problem in numerical methods.
pde maple mathematica
asked Jan 3 at 6:22
Arun Govind Neelan AArun Govind Neelan A
17027
$endgroup$
Convection-diffusion equation is:
$frac{partial u}{partial t} + frac{partial u}{partial x} = 0.01frac{partial^2 u}{partial x^2}$
Inital conditon is:
$u(x,0) = sin(x)$ over the domain 0 to $2pi$ with periodic boundary conditon that is $u(0,t) =u(2pi,t)$.
Second boundary condition is clamped at $x = 0 $ that is $frac{du}{dx} = 0$
I would like to know
1) Whether this problem has an exact solution? if so please prove the solution.
2) Can any symbolic computing software like Maple, Mathematica, Matlab can solve this problem analytically?
3) Please provide some good tutorial (external links) for finding the analytical solution of the advection-diffusion equation.
Please don't provide a numerical solution because this problem is a toy problem in numerical methods.
pde maple mathematica
pde maple mathematica
asked Jan 3 at 6:22
Arun Govind Neelan AArun Govind Neelan A
17027
asked Jan 3 at 6:22
Arun Govind Neelan AArun Govind Neelan A
17027
edited Jan 3 at 7:05
Arun Govind Neelan A
asked Jan 3 at 6:22
Arun Govind Neelan AArun Govind Neelan A
17027
asked Jan 3 at 6:22
Arun Govind Neelan AArun Govind Neelan A
17027
asked Jan 3 at 6:22
Arun Govind Neelan AArun Govind Neelan A
17027
17027
$begingroup$
Note, if $u$ solves the above convection-diffusion then $h(x, t)=u(x+t, t)$ solves begin{align} partial_t h =partial_t u + partial_x u= .01 partial_{xx} u = .01partial_{xx}h. end{align} Hence $h$ solves the heat equation.
$endgroup$
– Jacky Chong
Jan 3 at 7:14
$begingroup$
$u(x,t) = exp(-nu t)sin(x-ct)$. Is this the solution to this equation? I found "some" match with the numerical solution but not satisfied. I think I made some mistake. Solution following the physics but I'm not sure whether it is correct or not
$endgroup$
– Arun Govind Neelan A
Jan 3 at 7:29
$begingroup$
In addition to that grid refinement study shows a good match with analytical and numerical result. Could you tell me any way to cross check it? The solution basis is good but I don't know the impact of the boundary condition on the solution and how to cross check it.
$endgroup$
– Arun Govind Neelan A
Jan 3 at 7:43
add a comment |
$begingroup$
Note, if $u$ solves the above convection-diffusion then $h(x, t)=u(x+t, t)$ solves begin{align} partial_t h =partial_t u + partial_x u= .01 partial_{xx} u = .01partial_{xx}h. end{align} Hence $h$ solves the heat equation.
$endgroup$
– Jacky Chong
Jan 3 at 7:14
$begingroup$
$u(x,t) = exp(-nu t)sin(x-ct)$. Is this the solution to this equation? I found "some" match with the numerical solution but not satisfied. I think I made some mistake. Solution following the physics but I'm not sure whether it is correct or not
$endgroup$
– Arun Govind Neelan A
Jan 3 at 7:29
$begingroup$
In addition to that grid refinement study shows a good match with analytical and numerical result. Could you tell me any way to cross check it? The solution basis is good but I don't know the impact of the boundary condition on the solution and how to cross check it.
$endgroup$
– Arun Govind Neelan A
Jan 3 at 7:43
$begingroup$
Note, if $u$ solves the above convection-diffusion then $h(x, t)=u(x+t, t)$ solves begin{align} partial_t h =partial_t u + partial_x u= .01 partial_{xx} u = .01partial_{xx}h. end{align} Hence $h$ solves the heat equation.
$endgroup$
– Jacky Chong
Jan 3 at 7:14
$begingroup$
Note, if $u$ solves the above convection-diffusion then $h(x, t)=u(x+t, t)$ solves begin{align} partial_t h =partial_t u + partial_x u= .01 partial_{xx} u = .01partial_{xx}h. end{align} Hence $h$ solves the heat equation.
$endgroup$
– Jacky Chong
Jan 3 at 7:14
$begingroup$
$u(x,t) = exp(-nu t)sin(x-ct)$. Is this the solution to this equation? I found "some" match with the numerical solution but not satisfied. I think I made some mistake. Solution following the physics but I'm not sure whether it is correct or not
$endgroup$
– Arun Govind Neelan A
Jan 3 at 7:29
$begingroup$
$u(x,t) = exp(-nu t)sin(x-ct)$. Is this the solution to this equation? I found "some" match with the numerical solution but not satisfied. I think I made some mistake. Solution following the physics but I'm not sure whether it is correct or not
$endgroup$
– Arun Govind Neelan A
Jan 3 at 7:29
$begingroup$
In addition to that grid refinement study shows a good match with analytical and numerical result. Could you tell me any way to cross check it? The solution basis is good but I don't know the impact of the boundary condition on the solution and how to cross check it.
$endgroup$
– Arun Govind Neelan A
Jan 3 at 7:43
$begingroup$
In addition to that grid refinement study shows a good match with analytical and numerical result. Could you tell me any way to cross check it? The solution basis is good but I don't know the impact of the boundary condition on the solution and how to cross check it.
$endgroup$
– Arun Govind Neelan A
Jan 3 at 7:43
add a comment |
1 Answer
1
active
oldest
votes
$begingroup$
Let's rewrite
$$ 2ku_t = u_{xx} - 2ku_x $$
where $k=50$. Separation of variables $u(x,t) = X(x)T(t)$ gives
$$ 2kfrac{T'}{T} = frac{X''-2kX}{X} = -lambda $$
The $X$ part has the characteristic polynomial
$$ r^2 - 2kr + lambda = 0 $$
or $r = k pm sqrt{k^2-lambda} $
Due to the periodicity, we need the roots to have an imaginary part, therefore $k^2-lambda = -mu^2$, or $lambda = mu^2+k^2$. This makes
$$ X(x) = e^{kx}big[Asin(mu x) + Bcos(mu x) big] $$
The periodicity condition $X(x) = X(x+2pi)$ requires $mu = n$ where $n=1,2,3$.
The other boundary condition is $X'(0)=0$. Solving this gives $mu A + kB = 0$, we can write
$$ X_n(x) = e^{kx}big[ksin(nx) - ncos(nx)big] $$
up to a multiplicative constant.
Solving for $T(t)$ and using superposition, we have the general solution
$$ u(x,t) = sum_{n=1}^infty c_n e^{-frac{n^2+k^2}{2k}t}e^{kx}big[ksin(nx) - ncos(nx)big] $$
The initial condition is equivalent to solving for $c_n$ such that
$$ sum_{n=1}^infty c_nbig[ksin(nx) - ncos(nx)big] = e^{-kx}sin x $$
Using orthogonality, we have
$$ c_n = frac{int_0^{2pi} e^{-kx}sin x big[ksin(nx) - ncos(nx)big] dx}{int_0^{2pi} big[ksin(nx) - ncos(nx)big]^2 dx} $$
answered Jan 3 at 11:00
DylanDylan
12.6k31026
$endgroup$
$begingroup$
Sorry for the late response. I think the solution may not a good one. When I was plotting the result, I can't see any advection (travelling wave) in the solution. Based on the physics $sin(x)$ should move from one end to other and diffuse over the time.
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:51
$begingroup$
I have learnt solving advection and diffusion equation in engineering maths book, I havn't learn solving advection-diffusion equation. Could you refer some text book for that?
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:54
$begingroup$
Could you tell $u(x,t) =exp(-0.01 t) sin(x-t)$ is the solution for this. I got this solution by combining heat equation and wave equation solution and it is matching with numerical solution. Please note that the problem is little bit ill-posed because of gradient boundary conditon.
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:59
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%2f3060304%2fanalytical-solution-for-convection-diffusion-equation%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$
Let's rewrite
$$ 2ku_t = u_{xx} - 2ku_x $$
where $k=50$. Separation of variables $u(x,t) = X(x)T(t)$ gives
$$ 2kfrac{T'}{T} = frac{X''-2kX}{X} = -lambda $$
The $X$ part has the characteristic polynomial
$$ r^2 - 2kr + lambda = 0 $$
or $r = k pm sqrt{k^2-lambda} $
Due to the periodicity, we need the roots to have an imaginary part, therefore $k^2-lambda = -mu^2$, or $lambda = mu^2+k^2$. This makes
$$ X(x) = e^{kx}big[Asin(mu x) + Bcos(mu x) big] $$
The periodicity condition $X(x) = X(x+2pi)$ requires $mu = n$ where $n=1,2,3$.
The other boundary condition is $X'(0)=0$. Solving this gives $mu A + kB = 0$, we can write
$$ X_n(x) = e^{kx}big[ksin(nx) - ncos(nx)big] $$
up to a multiplicative constant.
Solving for $T(t)$ and using superposition, we have the general solution
$$ u(x,t) = sum_{n=1}^infty c_n e^{-frac{n^2+k^2}{2k}t}e^{kx}big[ksin(nx) - ncos(nx)big] $$
The initial condition is equivalent to solving for $c_n$ such that
$$ sum_{n=1}^infty c_nbig[ksin(nx) - ncos(nx)big] = e^{-kx}sin x $$
Using orthogonality, we have
$$ c_n = frac{int_0^{2pi} e^{-kx}sin x big[ksin(nx) - ncos(nx)big] dx}{int_0^{2pi} big[ksin(nx) - ncos(nx)big]^2 dx} $$
answered Jan 3 at 11:00
DylanDylan
12.6k31026
$endgroup$
$begingroup$
Sorry for the late response. I think the solution may not a good one. When I was plotting the result, I can't see any advection (travelling wave) in the solution. Based on the physics $sin(x)$ should move from one end to other and diffuse over the time.
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:51
$begingroup$
I have learnt solving advection and diffusion equation in engineering maths book, I havn't learn solving advection-diffusion equation. Could you refer some text book for that?
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:54
$begingroup$
Could you tell $u(x,t) =exp(-0.01 t) sin(x-t)$ is the solution for this. I got this solution by combining heat equation and wave equation solution and it is matching with numerical solution. Please note that the problem is little bit ill-posed because of gradient boundary conditon.
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:59
add a comment |
$begingroup$
Let's rewrite
$$ 2ku_t = u_{xx} - 2ku_x $$
where $k=50$. Separation of variables $u(x,t) = X(x)T(t)$ gives
$$ 2kfrac{T'}{T} = frac{X''-2kX}{X} = -lambda $$
The $X$ part has the characteristic polynomial
$$ r^2 - 2kr + lambda = 0 $$
or $r = k pm sqrt{k^2-lambda} $
Due to the periodicity, we need the roots to have an imaginary part, therefore $k^2-lambda = -mu^2$, or $lambda = mu^2+k^2$. This makes
$$ X(x) = e^{kx}big[Asin(mu x) + Bcos(mu x) big] $$
The periodicity condition $X(x) = X(x+2pi)$ requires $mu = n$ where $n=1,2,3$.
The other boundary condition is $X'(0)=0$. Solving this gives $mu A + kB = 0$, we can write
$$ X_n(x) = e^{kx}big[ksin(nx) - ncos(nx)big] $$
up to a multiplicative constant.
Solving for $T(t)$ and using superposition, we have the general solution
$$ u(x,t) = sum_{n=1}^infty c_n e^{-frac{n^2+k^2}{2k}t}e^{kx}big[ksin(nx) - ncos(nx)big] $$
The initial condition is equivalent to solving for $c_n$ such that
$$ sum_{n=1}^infty c_nbig[ksin(nx) - ncos(nx)big] = e^{-kx}sin x $$
Using orthogonality, we have
$$ c_n = frac{int_0^{2pi} e^{-kx}sin x big[ksin(nx) - ncos(nx)big] dx}{int_0^{2pi} big[ksin(nx) - ncos(nx)big]^2 dx} $$
answered Jan 3 at 11:00
DylanDylan
12.6k31026
$endgroup$
$begingroup$
Sorry for the late response. I think the solution may not a good one. When I was plotting the result, I can't see any advection (travelling wave) in the solution. Based on the physics $sin(x)$ should move from one end to other and diffuse over the time.
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:51
$begingroup$
I have learnt solving advection and diffusion equation in engineering maths book, I havn't learn solving advection-diffusion equation. Could you refer some text book for that?
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:54
$begingroup$
Could you tell $u(x,t) =exp(-0.01 t) sin(x-t)$ is the solution for this. I got this solution by combining heat equation and wave equation solution and it is matching with numerical solution. Please note that the problem is little bit ill-posed because of gradient boundary conditon.
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:59
add a comment |
$begingroup$
Let's rewrite
$$ 2ku_t = u_{xx} - 2ku_x $$
where $k=50$. Separation of variables $u(x,t) = X(x)T(t)$ gives
$$ 2kfrac{T'}{T} = frac{X''-2kX}{X} = -lambda $$
The $X$ part has the characteristic polynomial
$$ r^2 - 2kr + lambda = 0 $$
or $r = k pm sqrt{k^2-lambda} $
Due to the periodicity, we need the roots to have an imaginary part, therefore $k^2-lambda = -mu^2$, or $lambda = mu^2+k^2$. This makes
$$ X(x) = e^{kx}big[Asin(mu x) + Bcos(mu x) big] $$
The periodicity condition $X(x) = X(x+2pi)$ requires $mu = n$ where $n=1,2,3$.
The other boundary condition is $X'(0)=0$. Solving this gives $mu A + kB = 0$, we can write
$$ X_n(x) = e^{kx}big[ksin(nx) - ncos(nx)big] $$
up to a multiplicative constant.
Solving for $T(t)$ and using superposition, we have the general solution
$$ u(x,t) = sum_{n=1}^infty c_n e^{-frac{n^2+k^2}{2k}t}e^{kx}big[ksin(nx) - ncos(nx)big] $$
The initial condition is equivalent to solving for $c_n$ such that
$$ sum_{n=1}^infty c_nbig[ksin(nx) - ncos(nx)big] = e^{-kx}sin x $$
Using orthogonality, we have
$$ c_n = frac{int_0^{2pi} e^{-kx}sin x big[ksin(nx) - ncos(nx)big] dx}{int_0^{2pi} big[ksin(nx) - ncos(nx)big]^2 dx} $$
answered Jan 3 at 11:00
DylanDylan
12.6k31026
$endgroup$
Let's rewrite
$$ 2ku_t = u_{xx} - 2ku_x $$
where $k=50$. Separation of variables $u(x,t) = X(x)T(t)$ gives
$$ 2kfrac{T'}{T} = frac{X''-2kX}{X} = -lambda $$
The $X$ part has the characteristic polynomial
$$ r^2 - 2kr + lambda = 0 $$
or $r = k pm sqrt{k^2-lambda} $
Due to the periodicity, we need the roots to have an imaginary part, therefore $k^2-lambda = -mu^2$, or $lambda = mu^2+k^2$. This makes
$$ X(x) = e^{kx}big[Asin(mu x) + Bcos(mu x) big] $$
The periodicity condition $X(x) = X(x+2pi)$ requires $mu = n$ where $n=1,2,3$.
The other boundary condition is $X'(0)=0$. Solving this gives $mu A + kB = 0$, we can write
$$ X_n(x) = e^{kx}big[ksin(nx) - ncos(nx)big] $$
up to a multiplicative constant.
Solving for $T(t)$ and using superposition, we have the general solution
$$ u(x,t) = sum_{n=1}^infty c_n e^{-frac{n^2+k^2}{2k}t}e^{kx}big[ksin(nx) - ncos(nx)big] $$
The initial condition is equivalent to solving for $c_n$ such that
$$ sum_{n=1}^infty c_nbig[ksin(nx) - ncos(nx)big] = e^{-kx}sin x $$
Using orthogonality, we have
$$ c_n = frac{int_0^{2pi} e^{-kx}sin x big[ksin(nx) - ncos(nx)big] dx}{int_0^{2pi} big[ksin(nx) - ncos(nx)big]^2 dx} $$
answered Jan 3 at 11:00
DylanDylan
12.6k31026
edited Jan 3 at 18:38
answered Jan 3 at 11:00
DylanDylan
12.6k31026
answered Jan 3 at 11:00
DylanDylan
12.6k31026
answered Jan 3 at 11:00
DylanDylan
12.6k31026
12.6k31026
$begingroup$
Sorry for the late response. I think the solution may not a good one. When I was plotting the result, I can't see any advection (travelling wave) in the solution. Based on the physics $sin(x)$ should move from one end to other and diffuse over the time.
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:51
$begingroup$
I have learnt solving advection and diffusion equation in engineering maths book, I havn't learn solving advection-diffusion equation. Could you refer some text book for that?
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:54
$begingroup$
Could you tell $u(x,t) =exp(-0.01 t) sin(x-t)$ is the solution for this. I got this solution by combining heat equation and wave equation solution and it is matching with numerical solution. Please note that the problem is little bit ill-posed because of gradient boundary conditon.
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:59
add a comment |
$begingroup$
Sorry for the late response. I think the solution may not a good one. When I was plotting the result, I can't see any advection (travelling wave) in the solution. Based on the physics $sin(x)$ should move from one end to other and diffuse over the time.
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:51
$begingroup$
I have learnt solving advection and diffusion equation in engineering maths book, I havn't learn solving advection-diffusion equation. Could you refer some text book for that?
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:54
$begingroup$
Could you tell $u(x,t) =exp(-0.01 t) sin(x-t)$ is the solution for this. I got this solution by combining heat equation and wave equation solution and it is matching with numerical solution. Please note that the problem is little bit ill-posed because of gradient boundary conditon.
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:59
$begingroup$
Sorry for the late response. I think the solution may not a good one. When I was plotting the result, I can't see any advection (travelling wave) in the solution. Based on the physics $sin(x)$ should move from one end to other and diffuse over the time.
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:51
$begingroup$
Sorry for the late response. I think the solution may not a good one. When I was plotting the result, I can't see any advection (travelling wave) in the solution. Based on the physics $sin(x)$ should move from one end to other and diffuse over the time.
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:51
$begingroup$
I have learnt solving advection and diffusion equation in engineering maths book, I havn't learn solving advection-diffusion equation. Could you refer some text book for that?
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:54
$begingroup$
I have learnt solving advection and diffusion equation in engineering maths book, I havn't learn solving advection-diffusion equation. Could you refer some text book for that?
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:54
$begingroup$
Could you tell $u(x,t) =exp(-0.01 t) sin(x-t)$ is the solution for this. I got this solution by combining heat equation and wave equation solution and it is matching with numerical solution. Please note that the problem is little bit ill-posed because of gradient boundary conditon.
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:59
$begingroup$
Could you tell $u(x,t) =exp(-0.01 t) sin(x-t)$ is the solution for this. I got this solution by combining heat equation and wave equation solution and it is matching with numerical solution. Please note that the problem is little bit ill-posed because of gradient boundary conditon.
$endgroup$
– Arun Govind Neelan A
Jan 4 at 11:59
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%2f3060304%2fanalytical-solution-for-convection-diffusion-equation%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$
Note, if $u$ solves the above convection-diffusion then $h(x, t)=u(x+t, t)$ solves begin{align} partial_t h =partial_t u + partial_x u= .01 partial_{xx} u = .01partial_{xx}h. end{align} Hence $h$ solves the heat equation.
$endgroup$
– Jacky Chong
Jan 3 at 7:14
$begingroup$
$u(x,t) = exp(-nu t)sin(x-ct)$. Is this the solution to this equation? I found "some" match with the numerical solution but not satisfied. I think I made some mistake. Solution following the physics but I'm not sure whether it is correct or not
$endgroup$
– Arun Govind Neelan A
Jan 3 at 7:29
$begingroup$
In addition to that grid refinement study shows a good match with analytical and numerical result. Could you tell me any way to cross check it? The solution basis is good but I don't know the impact of the boundary condition on the solution and how to cross check it.
$endgroup$
– Arun Govind Neelan A
Jan 3 at 7:43