How to solve a system of PDAEs with eigenvalue

  • Context: Graduate 
  • Thread starter Thread starter ktsharp
  • Start date Start date
  • Tags Tags
    Eigenvalue System
Click For Summary
SUMMARY

This discussion focuses on solving a system of Partial Differential Algebraic Equations (PDAEs) numerically. The equations presented include terms for H, W, and c, with specific initial and boundary conditions. The numerical solution requires methods suitable for PDAEs, such as the method of lines or finite difference methods. Participants emphasize the importance of correctly implementing the boundary conditions and selecting appropriate numerical solvers for accurate results.

PREREQUISITES
  • Understanding of Partial Differential Algebraic Equations (PDAEs)
  • Familiarity with numerical methods for solving differential equations
  • Knowledge of boundary value problems and initial conditions
  • Experience with numerical software tools such as MATLAB or Python libraries (e.g., SciPy)
NEXT STEPS
  • Research the method of lines for solving PDAEs
  • Learn about finite difference methods for numerical solutions
  • Explore MATLAB's PDE toolbox for implementing boundary conditions
  • Investigate Python's SciPy library for numerical solvers applicable to PDAEs
USEFUL FOR

Researchers, engineers, and mathematicians involved in numerical analysis, particularly those working with complex systems modeled by Partial Differential Algebraic Equations.

ktsharp
Messages
8
Reaction score
0
I have the following system of partial differential algebraic equations:

[ tex ] \frac{1}{H_p}\frac{\partial H_p}{\partial t} = - \frac{\partial W_p}{\partial x} - \frac{f1(H_p,c_p,W_p)}{H_p}, [ \tex ]
[tex]\frac{1}{H_p}\frac{\partial}{\partial t}(H_p c_p} = - \frac{\partial}{\partial x}(W_p c_p) - \frac{f2(H_p,W_p,c_p)}{H_p}, [\tex]<br /> [tex]0 = f2(H_p,c_p,W_p) + f3(H_p,c_p,W_p). [\tex]<br /> <br /> with the following conditions:<br /> Hp(x,0) = 5<br /> Wp(x,0) = s1(x)<br /> cp(x,0) = s2(x)<br /> <br /> Hp(0,t) = s3(t)<br /> Wp(0,t) = W0<br /> Wp(L,t) = 0<br /> d(cp)/dx (L,t) = 0<br /> <br /> How can I solve this numerically?[/tex][/tex]
 
Physics news on Phys.org
type "/tex" not "\tex" (or "$$") :wink:
ktsharp said:
I have the following system of partial differential algebraic equations:

[tex]\frac{1}{H_p}\frac{\partial H_p}{\partial t} = - \frac{\partial W_p}{\partial x} - \frac{f1(H_p,c_p,W_p)}{H_p},[/tex]
[tex]\frac{1}{H_p}\frac{\partial}{\partial t}(H_p c_p) = - \frac{\partial}{\partial x}(W_p c_p) - \frac{f2(H_p,W_p,c_p)}{H_p},[/tex]
[tex]0 = f2(H_p,c_p,W_p) + f3(H_p,c_p,W_p).[/tex]
 
  • Like
Likes   Reactions: 1 person
[UPDATED]

I have the following system of partial differential algebraic equations:

[tex]\frac{1}{H}\frac{\partial H}{\partial t} = - \frac{\partial W}{\partial x} - \frac{f_1(H,c,W)}{H},[/tex]
[tex]\frac{1}{H}\frac{\partial}{\partial t}(H c) = - \frac{\partial}{\partial x}(W c) - \frac{f_2(H,W,c)}{H},[/tex]
[tex]0 = f_2(H,c,W) + f_3(H,c,W).[/tex]

with the following conditions:
H(x,0) = 5
W(x,0) = s1(x)
c(x,0) = s2(x)

H(0,t) = s3(t)
W(0,t) = W0
W(L,t) = 0
d(c)/dx (L,t) = 0

How can I solve this numerically?
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
4K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 5 ·
Replies
5
Views
4K
  • · Replies 3 ·
Replies
3
Views
4K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 5 ·
Replies
5
Views
3K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 8 ·
Replies
8
Views
5K
  • · Replies 20 ·
Replies
20
Views
3K