Solving a PDE Using Finite Difference Method

Click For Summary
SUMMARY

The discussion focuses on solving the partial differential equation (PDE) represented by the equation \(\frac{dP}{dt}-A*\frac{{d}^2P}{dx^2}-B*\frac{dC}{dt}=0\) using the finite difference method in MATLAB. The user attempts to combine the finite difference solution for \(P(x,t)\) with the known function \(C(x,t)\) by manipulating matrices, specifically through operations like \(P + B*C\) and \(2*P + B*C\). The conversation clarifies that this is indeed a PDE, despite initial confusion, and emphasizes the importance of boundary and initial conditions for accurate solutions.

PREREQUISITES
  • Understanding of finite difference methods for numerical solutions
  • Familiarity with MATLAB for matrix operations
  • Knowledge of partial differential equations and their characteristics
  • Comprehension of boundary and initial conditions in differential equations
NEXT STEPS
  • Study the implementation of finite difference methods in MATLAB
  • Research boundary and initial conditions for PDEs
  • Explore numerical methods for solving PDEs, such as the Crank-Nicolson method
  • Learn about stability and convergence criteria in finite difference schemes
USEFUL FOR

Mathematicians, engineers, and researchers involved in numerical analysis, particularly those working with partial differential equations and seeking to implement solutions using MATLAB.

kitz2
Messages
4
Reaction score
0
Hi

The equation is:

\frac{dP}{dt}-A*\frac{{d}^2P}{dx^2}-B*\frac{dC}{dt}=0

dP/dt=A*d2P/dx^2 was solved using a finite difference method. If the function C(x,t) is known, is it possible to solve the whole equation by using the finite difference solution as a supplement to the complete solution. I mean solving dP/dt = B*dC/dt, and add them together or something. In case, how do I go about it?

Im using MATLAB and from the first finite difference solution I have an matrix for P(x,t) with all values for different x and t, and I also have a matrix for C(x,t) with all values for x and t. What I did was adding the P matrix with the C matrix multiplied with B , which have equal size. When the result was not as I wanted I tried adding the P matrix with P(matrix)+B*C(matrix), so I had 2*P(matrix) + B*C(matrix). Then I got a result which looked alright, but obviosuly it is just a shot in the dark, and not "math".

Im not very good at math as you might notice :blushing:
 
Physics news on Phys.org
This is not a partial differential equation (pde) it is a second order ordinary differential equation (ode). Are there any boundary or initail conditions that must be satisfied ?
 
Boundary and inital conditions:

C=Co , t=0, 0<=x<=inf.
C=Cdf x-> inf , t >0
P=Po , t=0 0<=x<=inf.
P=Po t>0 x-> inf.
P=Pw x=0 t>0
 
CFDFEAGURU said:
This is not a partial differential equation (pde) it is a second order ordinary differential equation (ode). Are there any boundary or initail conditions that must be satisfied ?

It is a partial differential equation, both t and x are the independent variables and P the dependent function to find. C is a given function, but I don't see how to solve it. Certainly numerical ways should be possible but I can't help you on this.

coomast
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 2 ·
Replies
2
Views
4K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 7 ·
Replies
7
Views
5K
  • · Replies 23 ·
Replies
23
Views
4K