Need help in solving 2D Parabolic PDE through PDE toolbox.

3 次查看（过去 30 天）
Aditya
Aditya 2023-12-5
回答： Udit06 2023-12-11
Variables Description:
h is film thickness
p is pressure (Variable)
grid is written in x and z
got the coefficient of pde through ‘pdeCoefficients’
Formatted as a simplified equation.
Code Snippet:
syms x
h = 1 + cos(x);
syms p(x,z)
% pdeeq = diff(p,x,x) + diff(p,z,z) + (diff(h,x)/h)*diff(p,x) - diff(h,x)/h;
pdeeq = laplacian(p,[x, z]) + (diff(h,x)/h)*diff(p,x) - diff(h,x)/h;
coeffs = pdeCoefficients(pdeeq,p)
coeffs.f('show')
coeffs.c('show')
Error:
Warning: After extracting m, d, and c, some gradients remain. Writing all remaining terms to f.
Note:
Seems like this equation can’t be translated in divergence form.
Can this PDE be solved, if yes through PDE toolbox can anyone please guide?

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

回答（1 个）

Udit06
Udit06 2023-12-11
Hi Aditya,
I understand that you are trying to solve a 2-D parabolic PDE that cannot be translated in a form which can be solved using PDE toolbox as specified in the following documentation.
https://www.mathworks.com/help/pde/ug/equations-you-can-solve.html
You can solve such equations using the numerical methods for solving PDE such as finite difference method or finite element method. You can refer to the following resources to understand more about the finite difference method and finite element method respectively.
I hope this helps.

类别

在 Help Center 和 File Exchange 中查找有关 PDE Solvers 的更多信息

标签

产品


版本

R2023a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by