Numerical solution for the heat equation does not match the exact solution

Question

Erm 2024-9-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2150049-numerical-solution-for-the-heat-equation-does-not-match-the-exact-solution

评论： Torsten 2024-9-5

I write a matlab code for the problem, and i want to compare the numerucal and the exact solution at t=0.5 but the graph of the soultions are very different.

The problem is in the following:

% Parameters
L = 5; % Length of the rod
T = 1;  % Total time
Nx = 100; % Number of spatial points
Nt = 500; % Number of time steps
alpha = 0.01; % Thermal diffusivity
% Discretization
dx = 0.1; % Spatial step size
dt = 0.04; % Time step size
r = alpha * dt / dx^2; % Stability parameter for the scheme
% Initial condition
x = linspace(0, L, Nx);% generates Nx points. The spacing between the points is (L-0)/(Nx-1).
u = cos(x-0.5); %  initial condition
% Preallocate solution matrix
U = zeros(Nx, Nt+1); %  returns an Nx-by-Nt matrix of zeros
U(:,1) = u; %extracts all the elements from the first column of matrix
% Time-stepping loop
for n = 1:Nt
    % Update the solution using explicit scheme for the Heat Equation
    for i = 2:Nx-1
        U(i,n+1) = U(i,n) + r * (U(i+1,n) - 2*U(i,n) + U(i-1,n));
    end
end
Numerical_sol = U
% Plot results
t = linspace(0, T, Nt+1);
[X, T] = meshgrid(x, t);
% Exact solution
Exact= @(x, t) exp(-t) .* cos(pi * x - 0.5 * pi);
figure;
mesh(X, T, U');
xlabel('x');
ylabel('t');
zlabel('u(x,t)');
title('Heat Equation Solution');
figure
plot(x, U(:,1), 'o');

2 个评论
显示无隐藏无

John D'Errico 2024-9-5

You accepted the answer to this question. Then you asked it just again, 42 minutes ago. I closed the duplicate question, since it asks nothing you should not have learned here.

Torsten 2024-9-5

@John D'Errico

I didn't know that you are the Website Master taking the right to delete questions and answers that are actually in flow.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Torsten 2024-9-4

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2150049-numerical-solution-for-the-heat-equation-does-not-match-the-exact-solution#answer_1510724

移动：Torsten 2024-9-4

I cannot believe that the value for alpha has no influence on the solution.

So I guess that the exact solution maybe holds for alpha = 1.

10 个评论
显示 8更早的评论隐藏 8更早的评论

Erm 2024-9-4

在 MATLAB Online 中打开

% Parameters
L = 1; % Length of the rod
T = 1;  % Total time
Nx = 100; % Number of spatial points
Nt = 500; % Number of time steps
alpha = 0.01; % Thermal diffusivity
% Discretization
dx = 0.1; % Spatial step size
dt = 0.04; % Time step size
r = alpha * dt / dx^2; % Stability parameter for the scheme
% Initial condition
x = linspace(0, L, Nx);% generates Nx points. The spacing between the points is (L-0)/(Nx-1).
u = cos(x-0.5); %  initial condition
% Preallocate solution matrix
U = zeros(Nx, Nt+1); %  returns an Nx-by-Nt matrix of zeros
U(:,1) = u; %extracts all the elements from the first column of matrix
% Time-stepping loop
for n = 1:Nt
    % Update the solution using explicit scheme for the Heat Equation
    for i = 2:Nx-1
        U(i,n+1) = U(i,n) + r * (U(i+1,n) - 2*U(i,n) + U(i-1,n));
    end
end
Numerical_sol = U
% Plot results
t = linspace(0, T, Nt+1);
[X, T] = meshgrid(x, t);
% Exact solution
Exact= @(x, t) exp(-t) .* cos(pi * x - 0.5 * pi);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Define the range for x
X = linspace(0, 1, 100); % 100 points between 0 and 1
% Define a specific value for t
t = 1; % You can change this value to see different plots
% Compute the function values
f = exp(-t) .* cos(pi * X - 0.5 * pi);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure;
mesh(X, T, U');
xlabel('x');
ylabel('t');
zlabel('u(x,t)');
title('Heat Equation Solution');
figure
plot(x, U(:,1), 'o');
hold on 
plot(X, f, 'LineWidth', 2);
xlabel('x');
ylabel('f(x)');
title(['Plot of f(x) = exp(-t) .* cos(\pi x - 0.5 \pi) for t = ' num2str(t)]);
Hold off

Torsten 2024-9-4

编辑：Torsten 2024-9-4

在 MATLAB Online 中打开

The exact solution for the differential equation

du/dt = alpha * d^2u/dx^2

with initial condition

u(0,x) = cos(pi*(x-0.5))

and boundary conditions

u(t,0) = u(t,1) = 0

is given by

uexact(t,x) = exp(-alpha*pi^2*t) * cos(pi*(x-0.5))

That's the function you have to compare the numerical solution with.

syms x t alpha pi

u = exp(-alpha*pi^2*t) * cos(pi*(x-0.5));

diff(u,t)

ans =

alpha*diff(u,x,2)

ans =

Erm 2024-9-4

Thanks a lot

请先登录，再进行评论。

Numerical solution for the heat equation does not match the exact solution

2 个评论
显示无隐藏无

采纳的回答

10 个评论
显示 8更早的评论隐藏 8更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

Numerical solution for the heat equation does not match the exact solution

2 个评论 显示 无隐藏 无

采纳的回答

10 个评论 显示 8更早的评论隐藏 8更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

2 个评论
显示无隐藏无

10 个评论
显示 8更早的评论隐藏 8更早的评论